nonlin_sim.r

来自「是例子滤波方面的一个很好的算法模块」· R 代码 · 共 15 行

sn <- 1 ##observation sd
ss <- 10 ## system sd
sp<- 1 ##prior sd
sig<-ss 

n<-50 #number of observations
#generate the data
alpha<-rep(0,n)
alpha[1]<-sp*rnorm(1)
for(t in 2:n){
alpha[t]<- alpha[t-1]*0.5+25*alpha[t-1]/(1+alpha[t-1]^2)+8*cos(1.2*(t-1))+ss*rnorm(1)
}

z<- alpha*alpha/20+sn*rnorm(n)