set.seed(137)
time.series = c()
n = 1024
v = 25
w = rnorm(n, 0, 1) / v
x_ini = runif(1,0,1)
for(i in 1:n){
if(i == 1){
time.series[i] = 0.2 + 0.6*cos(2*pi*(i/n))*x_ini + w[i] #
} else{
time.series[i] = 0.2 + 0.6*cos(2*pi*(i/n))*time.series[i-1] + w[i]
}
}
res1 = fix.fit(time.series, 5, 1, type = "Legen")
fix.plot(res1, "Legen")
Run the code above in your browser using DataLab