|
一阶;红色0 20 40 60 80 1000 1时间增长图2:市场投资组合价值过程的累积对数增长Zu,来自(3.2);黑色市值过程的累积对数增长X【n】,来自(3.1);红色模拟的R代码n<-500N<-550T<-1000000G0<-。051SS0<-。04ss1<-。06ss2<-。06Tfac=100/TG<-Tfac*G0if(ss1==ss2)ss<-rep(ss1,N)else ss<-seq(ss1,ss2,(ss2-ss1)/(N-1))ss<-Tfac*ssSS<-Tfac*SS0S<-sqrt(ss)s<-sqrt(ss)x<-log(1:N)x<-x-log(sum(exp(x)))xx<-exp(x)z<-rep(1,N)g<-rep(0,N)gg<-rep(0,N)LT rep(0,N)GR<-NULLGG<-NULLT1<-2*Ttt<-0for(1:T1中的t){mu<-xx/sum(xx[1:N])g<-(.5-mu)*ss+g-.5*SSXX<-xxng2<-sum(g)/2x<-x-g+s*rnorm(N)X1<-xGS<-g+s*rnorm(1)xx0<-log(sum(exp(x+GS)))ii<-顺序(x)x0<-0for(k in 1:(N-1)){x0<-x0+x[ii[k]]if(k*x[ii[k+1]]-x0>ng2)break}xpiv<-2*(x0+ng2)/kx[ii[k:1]]<-xpiv-x[ii[1:k]]x<-排序(-x)xx<-exp(x+GS)yy<-和(xx)xx<-log(xx)if(t%%(t%/%10)==0)打印(c(t,k,g xx[1]))如果(t>T1-t){tt<-tt+1gg<-gg+gGG[tt]<-log(yy)GR[tt]<-xx0z<-z+xxX2<-cumsum(exp(-sort(-X1)))X1<-cumsum(exp(X1))LT<-LT+2*(X2-X1)/xx}z<-z[1:N]/和(z[1:N])gg<-gg/ttlt<-2*cumsum(gg)LT<-LT(1:n中的i)如果(LT[i]<=0){打印(c(i,LT[i],LT[i])LT[i]<-LT[i]}时间=序列(0,tt,tt/1000)时间[1]<-1GR<-cumsum(GR)[时间]GG<-cumsum(GG)[时间]ZZ<-ss/LTZZ<-exp(cumsum(ZZ[n:1])ZZ<-ZZ/sum(ZZ)图形。off()石英(\'1\',宽度=8,高度=6)绘图(时间*Tfac,GR,ylim=范围(GR,GG),类型=\'l\',xlab=\'Time\',ylab=\'Growth\')线条(时间*Tfac,GG,col=2)石英(\'2\',宽度=8,高度=6)绘图(-GG/Tfac,类型=\'l\',xlab=\'Rank\',ylab=\'g\\U k\')石英(\'3\',宽度=8,高度=6)绘图(z,类型=\'l\',ylim=范围(z,ZZ),对数=\'xy\',AB=\'Rank\',ylab=\'Weight\')线条(ZZ,col(2)
|