matlab声谱图代码

以下是一个简单的 Matlab 声谱图代码：

% 读取音频文件并转换为频谱图
filename = 'sample.wav'; % 输入音频文件名
[y, Fs] = audioread(filename); % 读取音频文件
win = hann(512); % 窗函数
noverlap = 256; % 重叠长度
nfft = 1024; % FFT 点数
[S, F, T, P] = spectrogram(y, win, noverlap, nfft, Fs); % 频谱图

% 绘制频谱图
figure;
imagesc(T, F, 10*log10(abs(S).^2)); % 绘制图像
axis xy; % 翻转 y 轴方向
xlabel('时间 (s)');
ylabel('频率 (Hz)');
title('声谱图');
colorbar; % 添加颜色条

代码解释：

1. 读取音频文件 `sample.wav`，并将其转换为 `y` 向量和采样率 `Fs`。
2. 定义窗函数 `win`、重叠长度 `noverlap` 和 FFT 点数 `nfft`。
3. 使用 `spectrogram` 函数将音频信号转换为频谱图，返回频谱图 `S`、频率向量 `F`、时间向量 `T` 和功率谱矩阵 `P`。
4. 使用 `imagesc` 函数绘制频谱图，`10*log10(abs(S).^2)` 是将功率谱转换为分贝单位的公式，`axis xy` 是翻转 y 轴方向。
5. 添加 x 轴和 y 轴标签、标题和颜色条。

注意，该代码只适用于单声道音频文件。如果您的音频文件是立体声，您需要将其转换为单声道或分别处理左右声道并将它们合并为单声道。