matlab 基音

Matlab中计算音频信号的基音可以使用自相关函数。具体步骤如下：

1. 读入音频信号，使用`audioread`函数。

2. 将音频信号转化为单声道信号，如果是多声道信号，则使用`mean`函数将多个声道合并为一个。

3. 对单声道信号进行预加重处理，使用`filter`函数。

4. 对预加重后的信号进行帧分析，将信号分成多个重叠的帧，使用`buffer`函数。

5. 对每一帧信号进行自相关分析，使用`xcorr`函数。

6. 找到自相关函数的第一个峰值，即为基音周期。

7. 计算基音频率，即为采样率除以基音周期。

代码示例：

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

% 将音频信号转化为单声道信号
x = mean(x, 2);

% 预加重处理
a = 0.97;
b = [1 -a];
x = filter(b, a, x);

% 帧分析
frame_len = 0.03; % 帧长为30ms
frame_shift = 0.01; % 帧移为10ms
frame_len_sample = round(frame_len * fs);
frame_shift_sample = round(frame_shift * fs);
frames = buffer(x, frame_len_sample, frame_len_sample-frame_shift_sample);

% 自相关分析
N = size(frames, 1);
pitch = zeros(N, 1);
for i = 1:N
    r = xcorr(frames(:, i));
    [~, locs] = findpeaks(r);
    pitch(i) = locs(1);
end

% 计算基音频率
pitch_freq = fs ./ pitch;