armodel_psd.m

用AR模型法进行功率谱估计的各种算法比较

clear all;
close all;

load test;
N=2048;
% 使用自相关法得到功率谱估计；
[xpsd,F]=pyulear(x,8,N,1);
pmax=max(xpsd);
xpsd=xpsd/pmax;
xpsd=10*log10(xpsd+0.000001);
figure('color','w');subplot(321);
plot(F,xpsd);grid on;
%title('pyulear 8');

[xpsd,F]=pyulear(x,11,N,1);
pmax=max(xpsd);
xpsd=xpsd/pmax;
xpsd=10*log10(xpsd+0.000001);
subplot(323);
plot(F,xpsd);grid on;
%title('pyulear 11');

[xpsd,F]=pyulear(x,30,N,1);
pmax=max(xpsd);
xpsd=xpsd/pmax;
xpsd=10*log10(xpsd+0.000001);
subplot(325);
plot(F,xpsd);grid on;
%title('pyulear 30');
%
[xpsd,F]=pburg(x,5,N,1);
pmax=max(xpsd);
xpsd=xpsd/pmax;
xpsd=10*log10(xpsd+0.000001);
subplot(322);
plot(F,xpsd);grid on;
%title('pburg 5');

[xpsd,F]=pburg(x,7,N,1);
pmax=max(xpsd);
xpsd=xpsd/pmax;
xpsd=10*log10(xpsd+0.000001);
subplot(324);
plot(F,xpsd);grid on;
%title('pburg 7');


[xpsd,F]=pburg(x,9,N,1);
pmax=max(xpsd);
xpsd=xpsd/pmax;
xpsd=10*log10(xpsd+0.000001);
subplot(326);
plot(F,xpsd);grid on;
%title('pburg 9');