clc,clear

randn('state',sum(clock));

elps=rand(1,100);

x(1)=elps(1);

for j=2:100

x(j)=0.39498*x(j-1)+elps(j);

end

但是产生的数据是符合差分后的AR（1）模型的，如何才能产生出类似于差分之前的数据啊？也是就说把产生的数据进行一阶差分后还能建立AR（1）模型。