基于MATLAB读取wav音频文件提取窄带语谱图并输出结果

您可以使用MATLAB中的`spectrogram`函数来提取窄带语谱图。以下是一个简单的示例代码，可以读取名为`example.wav`的音频文件，提取语谱图并将结果保存为`example_spectrogram.png`文件。

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

% 设置参数
window = hamming(round(fs*0.03)); % 窗函数
noverlap = round(fs*0.02); % 重叠长度
nfft = 1024; % FFT点数
f = linspace(0, fs/2, nfft/2+1); % 频率轴

% 提取语谱图
[S, F, T] = spectrogram(x, window, noverlap, nfft, fs);

% 绘制语谱图
imagesc(T, f, 20*log10(abs(S(1:nfft/2+1,:))));
axis xy;
xlabel('Time (s)');
ylabel('Frequency (Hz)');
title('Spectrogram');
colorbar;

% 保存语谱图
saveas(gcf, 'example_spectrogram.png');

在此示例中，我们使用了一个长度为30毫秒的汉明窗口，20毫秒的重叠长度和1024点FFT。您可以根据需要调整这些参数。注意，语谱图的y轴是频率轴，而不是频率的对数。如果您想要频率的对数刻度，请将y轴改为`logspace(log10(min(f)), log10(max(f)), 100)`。