p507.m

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

%program p507
lag=255;
[r,z]=xcorr(u,lag,'biased');
stem(z,r);
hndl=stem(z,r);   
set(hndl,'Marker','o')   
set(hndl,'MarkerSize',0.5)
xlabel('n');ylabel('r(n)');
title('Autocorrelation');