matlab中语音信号时频分析代码

以下是一个简单的 MATLAB 语音信号时频分析的代码示例：

% 读入语音信号
[x, fs] = audioread('speech.wav');

% 设置参数
winlen = 512; % 窗口长度
noverlap = round(winlen/2); % 重叠长度
nfft = winlen; % FFT 长度

% 计算短时傅里叶变换
[S, f, t] = stft(x, winlen, noverlap, nfft, fs);

% 绘制时频图
figure;
imagesc(t, f, abs(S));
axis xy;
xlabel('Time (s)');
ylabel('Frequency (Hz)');
title('STFT Magnitude');

% 定义短时傅里叶变换函数
function [S, f, t] = stft(x, winlen, noverlap, nfft, fs)
    % 初始化输出
    S = [];
    t = [];
    f = [];

    % 计算窗口函数
    win = hamming(winlen);

    % 计算帧长度和帧移
    hopsize = winlen - noverlap;
    nframes = fix((length(x)-noverlap)/hopsize);

    % 循环计算短时傅里叶变换
    for i = 1:nframes
        % 计算帧的开始和结束位置
        start = (i-1)*hopsize+1;
        stop = start+winlen-1;

        % 取出帧
        frame = x(start:stop);

        % 应用窗口函数
        frame = frame .* win;

        % 计算短时傅里叶变换
        X = fft(frame, nfft);

        % 取出幅度谱
        S(:,i) = abs(X(1:nfft/2+1));

        % 记录时间和频率轴
        if isempty(t)
            t = start/fs;
        else
            t = [t start/fs];
        end
        if isempty(f)
            f = linspace(0, fs/2, nfft/2+1);
        end
    end
end

该代码使用 `audioread` 函数读取语音信号，然后计算短时傅里叶变换（STFT）。STFT 的参数包括窗口长度、重叠长度和 FFT 长度。函数 `stft` 实现了短时傅里叶变换的计算。最后，该代码绘制了 STFT 幅度谱的时频图。