zhangtao 发表于 2011-12-17 10:10 
epoh老师,您好!
程序其他地方我基本能看懂,就是根据brown motion 进行模拟,
下面这部分双循环主要的 ...
for (i in 1:NReps)
代表模拟次数我先设为10,方便理解
for (j in 1:NSteps)
代表模拟天数,我设为1年252天(T=1,NSteps=252)
SPaths[i,j+1] = SPaths[i,j]*exp(nudt + sidt*rnorm(1))
这是stock price 模拟公式
if(i == 1) plot(x,SPaths[i,],"l", ylim=c(30,80),col =sample(colno,1))
先画出第1条模拟路径
if (i != 1) lines(x,SPaths[i,],ylim=c(30,80),col =sample(colno,1))
添加上第2,3,...NReps条模拟路径