pp_zfj.m

来自「文件包含现代信号处理的常用算法的matlab源代码」· M 代码 · 共 17 行

%产生信号序列
n=linspace(-pi,pi,1024);
x=2*sin(n*40)+cos(n*100);
figure,plot(x);
title('输入信号');
%计算序列的DFT
y1=fft(x);
figure,plot(abs(y1));
grid on
title('频谱');
%计算序列的PSD
y2=conj(y1);
Ppw=y1.*y2/length(n);
%绘图
figure,plot(Ppw);
grid on
title('周期图法功率谱估计(直接法)');