ar1_yulewalker.m

来自「对信号在AR模型下进行的的yuler_walker,和burg算法谱估计」· M 代码 · 共 28 行

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]);