我有一个功能,我需要用R每天使用1步绘制3年。
S(t)= S(0)exp(0.06t+0.20w(t)) #(1)
with S(0) =20
w(t) = standard Brownian movement
我有点受阻。我知道这段代码应该可行,但我不知道如何在我的" t"之间保持相同的长度。和我的w(t) 长度(t)= 901,长度(w(t))= 902。如何使它们相似????
t <- seq(0,900,length=901)
length(t)
v = matrix(rnorm(log(20)+0.06*t,sd=sqrt((0.20)^(2)*t)))
z = matrix(NA, ncol=900, nrow=5)
w = function(t)
{c(0,cumsum(v))}
length(w(t))
Bn <-log(20)+ 0.06*t +0.20*w(t)
Br <- log(Bn)
plot(t,Br,type="l",xlab="Temps")
for (i in 1:5) # I need to draw 5 path of the function(1)
{z[i,] = c(0,cumsum(Br[i,]))}
dim(z)
u= apply(z,2,mean) #mean of the 5 path
plot(t,z[1,],xlab="temps",type="l",ylab="Movement Brownian")
for (i in 1:5){lines(t,z[i,])}
lines(t,u,lwd=2,col="red")
我需要帮助。 谢谢你的时间。