你这个是AR(1)模型,可以用SAS很容易模拟出来,不过不会和里面图一模一样,因为这个取决于随机数的产生, 不过应该差不多。你可以改变 rannor 里面种子的数值,得到不同的序列
data a;
xlag=0;
do i=1 to 100;
x1=0.8*xlag+rannor(1);
xlag=x1;
output;
end;
drop xlag;
symbol l=1 w=1 i=join;
proc gplot data=a;
plot x1*i;
run;