www.pudn.com > estimate-ar.rar > p505.m

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