MATLAB读取.mat文件短时傅里叶变换和小波变换绘制时频图代码

以下是MATLAB读取.mat文件并进行短时傅里叶变换和小波变换绘制时频图的示例代码：

% 读取.mat文件中的信号
load signal.mat;
x = signal; % 信号
fs = 1000; % 采样率

% 短时傅里叶变换
window_length = 256; % 窗口长度
overlap = window_length/2; % 重叠长度
nfft = 1024; % 傅里叶变换点数
[S, F, T] = spectrogram(x, window_length, overlap, nfft, fs);

% 绘制时频图
figure;
imagesc(T, F, abs(S));
set(gca, 'YDir', 'normal');
xlabel('时间 (s)');
ylabel('频率 (Hz)');
title('短时傅里叶变换时频图');

% 小波变换
wname = 'db4'; % 小波基
level = 5; % 分解层数
[C, L] = wavedec(x, level, wname);
c = detcoef(C, L, level); % 提取细节系数
t = 1:length(c);
scales = 2.^(0:level); % 尺度
coefs = cwt(c, scales, wname); % 连续小波变换

% 绘制时频图
figure;
imagesc(t, scales, abs(coefs));
set(gca, 'YDir', 'normal');
xlabel('时间 (s)');
ylabel('尺度');
title('小波变换时频图');

其中，`signal.mat` 文件中存放的是信号数据，可以自行准备或者使用 MATLAB 自带的数据集。`window_length` 表示窗口长度，`overlap` 表示重叠长度，`nfft` 表示傅里叶变换点数，在短时傅里叶变换中用于控制时频图分辨率和平滑度。`wname` 表示小波基名称，`level` 表示小波变换的分解层数，在小波变换中用于控制时频图分辨率和频带宽度。