ex_12_12_1.m

来自「是属于清华大学出的书 数字信号处理及电子课件光盘内容」· M 代码 · 共 57 行

%  ex_12_12_1.m :  习题 12.12


clear all;
close all;

load testdat;
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');