MATLAB中从语音信号获取基频的代码

获取语音信号的基频可以采用自相关法或基于峰值拟合的方法。以下是一个基于自相关法的示例代码：

% 读取语音文件
[s, fs] = audioread('speech.wav');

% 设置自相关函数的参数
frame_length = round(20 * fs / 1000); % 每帧的长度为20ms
frame_shift = round(10 * fs / 1000); % 相邻帧的重叠长度为10ms

% 计算自相关函数
acf = zeros(frame_length, length(s));
for i = 1 : length(s) - frame_length
    frame = s(i : i + frame_length - 1);
    acf(:, i) = xcorr(frame);
end

% 找到自相关函数的峰值
[max_acf, max_lag] = max(acf(1 : round(frame_length / 2), :), [], 1);

% 计算基频
f0 = fs ./ max_lag;
f0(isnan(f0)) = 0; % 处理没有检测到基频的帧

% 显示基频轨迹
t = (0 : length(s) - 1) / fs;
subplot(2, 1, 1);
plot(t, s);
xlabel('Time (s)');
ylabel('Amplitude');
title('Speech signal');
subplot(2, 1, 2);
plot(t(1 : end - frame_length), f0);
xlabel('Time (s)');
ylabel('Fundamental frequency (Hz)');
title('F0 contour');

需要注意的是，这个方法并不是完全准确的，因为语音信号的基频会受到许多因素的影响，例如说话人的性别、年龄、情感等等。因此，如果需要更精确的基频估计，可以尝试其他方法，例如基于倒谱法的方法。