MATLAB短时傅里叶变换并绘制时频图

### 绘制短时傅里叶变换的三维时频图 为了在 MATLAB 中绘制短时傅里叶变换（STFT）的三维时频图，可以采用两种主要方法：一种是利用内置函数 `spectrogram` 或者较新版本中的 `stft` 函数来获取 STFT 的数据；另一种则是基于这些数据构建三维图形。 #### 使用 `spectrogram` 获取 STFT 数据并绘图 通过设置合适的参数选项，可以从 `spectrogram` 函数获得用于创建三维表面图的数据矩阵。具体来说： ```matlab % 设置信号和采样率 Fs = 1000; % 采样频率 (Hz) t = 0:1/Fs:2-1/Fs; x = chirp(t,0,1,150); % 线性扫频信号作为例子 % 计算短时傅立叶变换 windowSize = 128; noverlap = windowSize / 2; nfft = 512; [s,f,t,p] = spectrogram(x,hamming(windowSize),noverlap,nfft,Fs); % 将功率谱密度转换成dB单位 P = 10*log10(p); ``` 接着使用 `surf()` 来制作三维图表，并调整视角以便更好地观察结果: ```matlab figure; surf(t,f,P,'EdgeColor','none'); axis tight; view(-45,60); % 修改视点角度以优化显示效果 xlabel('Time (seconds)'); ylabel('Frequency (Hz)'); zlabel('Power/frequency (dB/Hz)'); title('Spectrogram of Chirp Signal Using surf() Function'); colorbar; ``` 此段代码会生成一个具有时间轴、频率轴以及强度表示的颜色条形图，其中颜色深浅代表不同位置上的能量大小[^1]。 #### 利用 `stft` 函数简化流程 对于拥有更新版 MATLAB 的用户而言，可以直接运用更为简便高效的 `stft` 函数来进行相同的操作。这不仅减少了编写额外辅助脚本的需求，而且提高了计算效率。 ```matlab % 对同一组测试数据再次执行操作 [S,F,T] = stft(x,'Window',hamming(windowSize),'OverlapLength',noverlap,... 'FFTLength',nfft,'SampleRate',Fs); % 显示三维图像 figure; waterfall(T,F,abs(S)); shading interp; colormap jet; xlabel('Time (sec)') ylabel('Frequency (kHz)') zlabel('|STFT|') title({'Waterfall Plot Showing Magnitude Spectra Over Time'; ... '(Using |stft| function)'}) c=colorbar; c.Label.String='Magnitude'; set(gca,'YScale','log') % 可选：使y轴变为对数刻度提高可读性 ``` 上述两部分展示了如何分别借助传统方式 (`spectrogram`) 和现代工具 (`stft`) 实现目标——即构造出能够直观反映音频或其他类型波形特征变化过程的三维可视化作品[^3]。