Matlab繪制短時(shí)傅里葉變換的頻譜圖和時(shí)間-頻率-幅值三維圖

?

function [t,frequency,f_spectrum]=fft_s(y,windowlength,Fs) % 輸入 : % y-輸入信號(hào) % windowlength-窗長度 0-1的系數(shù)，比如windowlength-0.5；采樣頻率1000，則窗長度為500 % Fs-采樣頻率 windowlength = windowlength*Fs; y = reshape(y,windowlength,[]); % 數(shù)據(jù)分段s = fft(y); % 快速傅里葉變換d = abs(s(1:windowlength/2,:));% 求絕對(duì)值 f_spectrum = 2/windowlength*d; % 求幅值t = (1:size(y,2))*windowlength/Fs; % 計(jì)算時(shí)間 frequency = (0:windowlength/2-1)/windowlength*Fs; % 頻率范圍

總結(jié)

以上是生活随笔為你收集整理的2021-02-28 Matlab绘制短时傅里叶变换的频谱图和时间-频率-幅值三维图的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。