用MATLAB写一段绘制语谱图的代码

以下是一个简单的绘制语谱图的MATLAB代码：

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

% 设置帧长和帧移
frameLength = 256;
frameShift = 128;

% 计算帧数和频谱点数
numFrames = fix((length(x) - frameLength + frameShift) / frameShift);
numFreqs = frameLength / 2 + 1;

% 初始化语谱图矩阵
spectrogram = zeros(numFreqs, numFrames);

% 计算汉明窗
window = hamming(frameLength);

% 计算每一帧的语谱图
for i = 1:numFrames
    % 计算当前帧的起始和终止位置
    start = (i - 1) * frameShift + 1;
    stop = start + frameLength - 1;
    
    % 应用汉明窗
    xw = x(start:stop) .* window;
    
    % 计算快速傅里叶变换
    X = fft(xw);
    
    % 计算功率谱密度
    P = abs(X(1:numFreqs)) .^ 2 / frameLength;
    
    % 将功率谱密度存储到语谱图矩阵中
    spectrogram(:, i) = P;
end

% 绘制语谱图
imagesc(log10(spectrogram));
axis xy;
xlabel('帧数');
ylabel('频率(Hz)');
title('语谱图');
colorbar;

这段代码可以绘制一个简单的语谱图，其中x是音频信号，fs是采样率，frameLength和frameShift分别是帧长和帧移。语谱图矩阵spectrogram存储了每一帧的功率谱密度，使用log10将功率谱密度转换为dB单位后绘制出来。