告诉你一个思路,估计会误导你,但从另一方面来讲也可以开拓一下眼界。现在的时序分析多半都是在时域中搞,直到把你搞的一个头两个大为止。
验证一个平稳过程是纯随机过程最简单的办法就是看随机过程的频域分布是否形成等腰三角形,如果是,那么过程就是纯随机过程。
代码如下,自己跑一遍。
<br />
set.seed(123456)<br />
x<-rnorm(100,0,1)<br />
par(mfrow=c(1,2))<br />
plot(x,type='l',xlab='时域形式的白噪声',ylab='方差')<br />
h=hist(x,plot=F)<br />
barplot(h$counts,horiz=T,xlab='频域形式的白噪声')<br />
box()<br />
</p>