提取音频文件中语音信号的共振峰并画图的Matlab代码

以下是使用MATLAB提取音频文件中语音信号的共振峰并画图的示例代码：

% 读取音频文件
[data, fs] = audioread('sample.wav');

% 预处理信号
data = data(:, 1); % 取其中一路声道
data = data - mean(data); % 去除直流分量
win = hamming(length(data)); % 加窗
data = data .* win;

% 计算信号的频谱
N = length(data);
f = (0:N-1) * fs / N;
spec = abs(fft(data));

% 寻找谷值点
[pks, locs] = findpeaks(spec);
[~, idx] = sort(pks, 'descend');
peaks = locs(idx(1:5)); % 取前5个峰值

% 画图
figure;
subplot(2,1,1);
plot((0:N-1)/fs, data, 'b');
xlabel('时间/s');
ylabel('幅值');
title('原始语音信号');

subplot(2,1,2);
plot(f, spec, 'b');
hold on;
stem(f(peaks), spec(peaks), 'r');
hold off;
xlabel('频率/Hz');
ylabel('幅值');
title('语音信号频谱及共振峰');

解释一下代码的主要过程：

1. 读取音频文件并预处理信号，包括去直流分量和加窗；

2. 计算信号的频谱，使用MATLAB的fft函数；

3. 使用MATLAB的findpeaks函数寻找频谱中的峰值点；

4. 利用MATLAB的plot和stem函数将原始语音信号和共振峰绘制在同一张图上。

注意：示例代码中只取了频谱中前5个峰值，实际应用中可以根据需要调整。