你这个是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;

追加200.

2# 耕耘使者 感觉该模型是依据某些数据拟合得到的，作者在写书时候去了数据。

不是。
作者不是在判断序列是否平稳，而是在判断模型是否平稳，这是两个概念。
而且，作者明确强调了是用拟合值作的图，而非原始时间序列。

追加300