%  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');