first1.m

这是数字信号处理的作业的第一部分包括了lms算法的基本步骤

a1=-0.195;
a2=0.95;
vt=(1-a2)*[(1+a2)^2-a1^2]/(1+a2)
r(1)=-a1/(1+a2)；
r(2)=-a2+a1^2/(1+a2)
p1=1/2*[-a1+(a1^2-4*a2)^(1/2)]
p2=1/2*[-a1-(a1^2-4*a2)^(1/2)]
for m=3:255
    r(m)=-a1*r(m-1)-a2*r(m-2);
end;        %计算r(m)
Pbt=r0;
Par=0;
w=0:2*pi/999:2*pi;
for n=1:255
    Pbt=Pbt+r(n)*[exp(j*w*n)+exp(-j*w*n)];
end;
Par=vt./[abs(1+a1*exp(-j*w)+a2*exp(-j*2*w))].^2;
subplot(2,1,1);
plot(w,Pbt);
xlabel('w');
title('Pbt ');        %在第一子图绘制BT法估计的PSD
subplot(2,1,2);
plot(w,Par);
xlabel('w');
title('Par ');