如何用窗函数来估计功率谱密度？

jakens

怎么能用窗函数估计离散信号的功率谱密度？

jakens

小弟的水平比较低，有一个信号已经用直接法和间接法对信号做了功率谱的分析，还有一个要给信号加一个窗函数再分析一下功率谱，希望有会的大哥帮一下小弟的忙

vincentsuen

clear;
clc;

fs=50;
ns=512;
t=1:ns;
t=t/fs;
z=sin(4*pi*t);

figure(1);
plot(t,z(1:ns));
title('原始时程信号');

[Pxx,w]=periodogram(z,hamming(512),ns,fs,'onesided');
figure(2);
plot(w,Pxx);
title('功率谱函数');
xlabel('频率（Hz）');
ylabel('S(w)');
grid;