如何解决R中随机过程的模拟
我正在尝试使用不变性原理在 [0,1] 上模拟 Ornstein-Uhlenbeck 过程。 Mu、sigma 和 theta 是选定的常数。我用 R 试过这段代码：
#(1) Invariance principle
n <- 1000
j<- 10
t <- seq(0,1,length=n)
theta <- 1
delta <- diff(exp(2*theta*t))
mu <- 1.2
sigma <- 0.3
U = replicate(j,{
bm <- c(0,cumsum(rnorm(n-1,sqrt(delta))))
ou <- mu + sigma/(2*sqrt(theta))*exp(-theta*t)*bm[exp(2*theta*t)]
})
matplot(U,type = "l")
但是作为输出，我得到了非常奇怪的图表。有人能找出我的解决方案有什么问题吗？
版权声明：本文内容由互联网用户自发贡献，该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容， 请发送邮件至 dio@foxmail.com 举报，一经查实，本站将立刻删除。
