p505.m

这是在AR模型的功率谱估计

%program p505
%blackman-pbt
k=0:1000;
w=(pi/500)*k;
M=k/500;
X=r*(exp(-j*pi/500)).^(m'*k);
magX=abs(X);
plot(M,10*log10(magX));
xlabel('Angular frequency,ω (π)');
ylabel('Power Spectrum Magnitude (dB)');
title('Blackman-Turkey Spectral Estimate')