%************************************%
%     自相關法功率譜估計(間接法)     %
%************************************%
%產生信號序列
n=linspace(0,2*pi,1024);
xn=2*sin(n*40)+cos(n*100);
figure,plot(xn);
title('輸入信號');
%計算序列的自相關函數
cxn=xcorr(xn,'unbiased');
%計算序列的PSD
CXk=fft(cxn);
Pcw=abs(CXk);
%繪圖
figure,plot(Pcw);
grid on
title('自相關法功率譜估計(間接法)');