clear all;
Fs=1000;
n=0:1/Fs:0.3;
w0=200*pi;
w1=400*pi;
xn=cos(w0*n)+sin(w1*n)+randn(size(n));
subplot(211)
plot(n,xn);
xlabel('時間（s）')
ylabel('幅度')
title('xn=cos(w0*n)+sin(w1*n)+randn(size(n))');
p=floor(length(xn)/3)+1;
nfft=1024;
%range=whole;
[xpsd,f]=pyulear(xn,p,nfft,Fs);
%繪功率譜估計
pmax=max(xpsd);
xpsd=xpsd/pmax;
xpsd=10*log10(xpsd+0.000001);
subplot(212)
plot(f,xpsd);
title('基于yule-walker方法的功率譜估計')
ylabel('功率譜密度dB')
xlabel('頻率Hz')
grid on;
ymin=min(xpsd)-2;
ymax=max(xpsd)+2;
axis([0,Fs/2,ymin,ymax]);