短时傅里叶变换语谱图代码

短时傅里叶变换（Short-Time Fourier Transform，STFT）是一种时频分析方法，它将信号分成若干个时间段，在每个时间段内进行傅里叶变换，得到该时间段内信号的频谱。通过对所有时间段内的频谱进行叠加，可以得到整个信号的时频分布情况，即语谱图。

下面是一个MATLAB实现的STFT语谱图代码示例：

% 读取音频文件
[x, fs] = audioread('example.wav');

% 设置参数
win = hamming(256); % 窗函数
noverlap = 128; % 重叠长度
nfft = 512; % FFT长度

% 计算STFT
[S, F, T] = stft(x, win, noverlap, nfft, fs);

% 绘制语谱图
imagesc(T, F, abs(S));
axis xy;
xlabel('Time (s)');
ylabel('Frequency (Hz)');
colorbar;

其中，`stft`函数是MATLAB自带的STFT函数，它的输入参数依次为信号`x`、窗函数`win`、重叠长度`noverlap`、FFT长度`nfft`和采样率`fs`，输出参数为STFT结果`S`、频率向量`F`和时间向量`T`。

需要注意的是，语谱图的颜色表示的是信号的幅度，而不是功率或能量。如果需要绘制功率谱密度图，可以将`abs(S)`替换为`abs(S).^2`。