MATLAB中显示X、Y、Z三个方向信号的时频特征，使用时频分析方法

---

方法1：使用短时傅里叶变换（STFT）

% 示例数据（假设已有x, y, z三个方向的信号）
fs = 1000;          % 采样频率
t = 0:1/fs:1;       % 时间向量
x = chirp(t, 0, 1, 100);  % X方向信号（线性调频）
y = sin(2*pi*50*t);       % Y方向信号（50Hz正弦波）
z = randn(size(t));       % Z方向信号（噪声）% 时频分析
figure
subplot(3,1,1)
spectrogram(x, 256, 250, 256, fs, 'yaxis');
title('X方向时频图');subplot(3,1,2)
spectrogram(y, 256, 250, 256, fs, 'yaxis');
title('Y方向时频图');subplot(3,1,3)
spectrogram(z, 256, 250, 256, fs, 'yaxis');
title('Z方向时频图');

---

方法2：使用连续小波变换（CWT）- 更推荐

% 小波变换（需要Wavelet Toolbox）
figure
subplot(3,1,1)
cwt(x, 'amor', fs);
title('X方向小波时频图');subplot(3,1,2)
cwt(y, 'amor', fs);
title('Y方向小波时频图');subplot(3,1,3)
cwt(z, 'amor', fs);
title('Z方向小波时频图');

---

方法3：使用Hilbert-Huang变换（适合非平稳信号）

% 经验模态分解+希尔伯特谱（需要Signal Processing Toolbox）
figure
subplot(3,1,1)
hht(emd(x), fs);  % 对X方向信号进行EMD分解后显示Hilbert谱
title('X方向Hilbert谱');subplot(3,1,2)
hht(emd(y), fs);
title('Y方向Hilbert谱');subplot(3,1,3)
hht(emd(z), fs);
title('Z方向Hilbert谱');

---

方法4：自定义STFT绘图（更灵活控制）

% 自定义STFT参数
window = hamming(256);
noverlap = 200;
nfft = 1024;figure
subplot(3,1,1)
[s,f,t] = spectrogram(x, window, noverlap, nfft, fs);
imagesc(t, f, 10*log10(abs(s)));
axis xy; colorbar; ylabel('频率 (Hz)'); title('X方向时频图');subplot(3,1,2)
[s,f,t] = spectrogram(y, window, noverlap, nfft, fs);
imagesc(t, f, 10*log10(abs(s)));
axis xy; colorbar; ylabel('频率 (Hz)'); title('Y方向时频图');subplot(3,1,3)
[s,f,t] = spectrogram(z, window, noverlap, nfft, fs);
imagesc(t, f, 10*log10(abs(s)));
axis xy; colorbar; ylabel('频率 (Hz)'); title('Z方向时频图');

---

关键参数说明：

• 窗函数：hamming/hann/kaiser等
• 重叠点数：影响时频图平滑度
• FFT点数：影响频率分辨率
• 颜色映射：使用colormap(jet)可改变颜色

推荐选择：

1. 对于一般振动/生物信号 → 使用小波变换（方法2）
2. 需要快速实现 → 使用spectrogram（方法1）
3. 非线性非平稳信号 → 使用HHT（方法3）