dsp142.m

经过调试的平稳随机信号的AR模型

% question three point one
p=0.0965;
N=2048;
a=sqrt(p);
v=rand(1,N);
v=v/std(v);
v=v-mean(v);
v=a*v;
%产生在第一组数据下的u（n）序列
a1=-0.195;
a2=0.95;
u1=0;
u2=0;
u=[u1,u2];
for n=3:1:2048
    u(n)=-a1*u(n-1)-a2*u(n-2)+v(n);
end
%产生在第一组数据下估计的r（n）序列
r=xcorr(u,1024,'coeff');
r=r(1025:2048);
SBT=fft(r,255);
w=-pi:1/127*pi:pi;
plot(w,SBT)
title('BT method')
ylabel('Sbt(w)')
xlabel('w')