waveletpsd.m

《Matlab时频分析及其应用》的详细代码

function WaveletPSD(x)
%输入x-各层的小波
nfft=1024;
fs=1;
y=abs(fft(x,nfft));
plot((0:nfft/2-1)/nfft*fs,y(1:nfft/2));
xlabel('频率 f/Hz');
ylabel('功率谱 P/W');