pc_zfj.m

来自「使用matlab编写的8个现代信号处理程序」· M 代码 · 共 17 行

%************************************%
%     自相关法功率谱估计(间接法)     %
%************************************%
%产生信号序列
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('自相关法功率谱估计(间接法)');