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

Question

所以基本上可以说我有一个函数 X 可以计算粒子在一维中的随机运动。该函数有不同的常数和一个正常的随机变量 W，每条路径每 0.1 毫秒发生一次。我想模拟 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)
}

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

标签： rpathsimulate

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)
}