clear all; M=200; nn=randn(); v1(1)=nn; v2(1)=nn; d(1)=sin(0.05*pi); for ii=2:M nn=randn(); v1(ii)=0.8*v1(ii-1)+nn; v2(ii)=-0.6*v2(ii-1)+nn; d(ii)=sin(0.05*pi*ii); end x=d+v1; figure; plot(x);6
d (n) = sin(0.05π n)v1 (n) = 0.8v1 (n 1) + w(n)
v2 (n) = 0.6v2 (n 1) + w(n)w(n)是零均值、单位方差的白噪声 是零均值、
x(n) = d (n) + v1 (n)