ar.m

ar模型 bt算法 ls rls lms music等数字信号处理的源代码

clear;
clc;
a1=-1.595;
a2=0.95;
r(1)=1;
sigma_u=r(1);
sigma_v=sigma_u*(1-a2)*((1+a2)^2-a1^2)/(1+a2);
%%%%%%%%%%%  calculate r(m)    %%%%%%%%%%%%%%
r(2)=-a1*r(1)/(1+a2);
for m=3:256
    r(m)=-a1*r(m-1)-a2*r(m-2);
end
%%%%%%%%%%%%%%%   BT method    %%%%%%%%%%%%%
r1=[fliplr(r(:,2:256)) r];
k=-100:100;
n=-255:255;
P_BT=r1*exp(-j*pi/100).^(n'*k);
plot(pi/100*k,10*log10(P_BT/max(P_BT)))
title('BT method PSD curve')
xlabel('w')
ylabel('normalized PSD(db)')
%%%%%%%%%%%%%%%   AR method    %%%%%%%%%%%%%
w=k*pi/100;
P_AR=(abs(1./(1+a1*exp(-j*w)+a2*exp(-2*j*w)))).^2*sigma_v;
figure
plot(w,10*log10(P_AR/max(P_AR)))
title('AR method PSD curve')
xlabel('w')
ylabel('normalized PSD(db)')