0

我有一个函数，我需要使用 R 绘制 3 年的每日步长 1。

S(t)= S(0)exp(0.06t+0.20w(t))  #(1)
with S(0) =20
w(t) = standard Brownian movement

我有点受阻。我知道这段代码应该可以工作，但我不知道如何在我的“t”和我的 w(t) length(t)=901 和 length(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")

请我需要帮助。感谢您的时间。

4

0 回答 0