2021-02-28 Matlab绘制短时傅里叶变换的频谱图和时间-频率-幅值三维图

Matlab绘制短时傅里叶变换的频谱图和时间-频率-幅值三维图

 

function [t,frequency,f_spectrum]=fft_s(y,windowlength,Fs)
% 输入 : 
% y-输入信号
% windowlength-窗长度 0-1的系数，比如windowlength-0.5；采样频率1000，则窗长度为500
% Fs-采样频率
windowlength = windowlength*Fs;
y = reshape(y,windowlength,[]); % 数据分段

s = fft(y); % 快速傅里叶变换

d = abs(s(1:windowlength/2,:));% 求绝对值
f_spectrum = 2/windowlength*d; % 求幅值

t = (1:size(y,2))*windowlength/Fs; % 计算时间
frequency = (0:windowlength/2-1)/windowlength*Fs; % 频率范围