试图构建一个函数来模拟R中粒子的100条路径

Question

所以基本上可以说我有一个函数X，它将计算一维粒子的随机运动。该函数具有不同的常数，以及一个正常的随机变量W，每条路径每隔0.1ms发生一次。我想模拟100条路径。

X <- 0;
Dt <- 0.0001;
V <- 0.5;
for (j in 0:100){
    W <- rnorm(100, j*Dt*V,1);
    x[0] = 0;
    x[j] = x[j-1] + Dt*V+ W*sqrt(Dt)
}

但是我得到一个错误，说“替换长度为零”，同样在获得了粒子不同位置的位置后，我想对其进行模拟，但是我不确定如何做到这一点。谢谢

Answer 1

X <- array()
Dt <- 0.0001
V <- 0.5
X[1] = 0
for (j in 2:101){
    W <- rnorm(100, j*Dt*V,1)
    X[j] = X[j-1] + Dt*V+ W*sqrt(Dt)
}