继续上述代码，绘制基频曲线

在计算出基音频率之后，我们可以继续绘制基频曲线。下面是一个基于Matlab的绘制基频曲线的代码示例：

% 定义帧长和帧移
frame_len = round(0.03*Fs);
frame_shift = round(0.01*Fs);

% 对语音信号进行分帧
num_frames = floor((length(y)-frame_len)/frame_shift)+1;
frames = zeros(frame_len, num_frames);
for i = 1:num_frames
    idx = (i-1)*frame_shift+1;
    frames(:,i) = y(idx:idx+frame_len-1);
end

% 计算每一帧的基音频率
f0 = zeros(num_frames, 1);
for i = 1:num_frames
    r = xcorr(frames(:,i), max_delay, 'coeff');
    [max_val, max_idx] = max(r);
    p = max_idx;
    f0(i) = Fs/p;
end

% 绘制基频曲线
t = (0:num_frames-1)*frame_shift/Fs;
figure;
plot(t, f0);
xlabel('Time (s)');
ylabel('Fundamental Frequency (Hz)');

首先，我们定义了帧长和帧移，分别为30ms和10ms。然后，对语音信号进行分帧，计算每一帧的基音频率。最后，通过绘制基频曲线来展示基音频率随时间的变化。