用MATLAB写一段通过识别声音频率的方式提取音频中的人声部分的代码

以下是MATLAB代码，使用基于频率域的人声分离算法：

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

% 短时傅里叶变换
win = hamming(round(fs*0.03)); % 30ms 窗口
hop = round(fs*0.015); % 15ms 帧移
nfft = 2^nextpow2(length(win)); % FFT 点数
[S,F,T] = spectrogram(y,win,hop,nfft,fs);

% 计算能量谱
P = abs(S).^2;

% 估计噪声能量谱
noiseP = mean(P(:,1:10),2); % 前 10 帧作为噪声参考

% 估计信号能量谱
signalP = max(P - noiseP, 0);

% 提取人声部分
theta = 5; % 阈值
mask = signalP ./ (signalP + noiseP); % 二值掩码
mask(mask < theta) = 0; % 低于阈值的部分设置为 0
mask(mask >= theta) = 1; % 高于等于阈值的部分设置为 1
voiceP = mask .* signalP; % 人声能量谱

% 合成音频信号
voiceS = voiceP.^(1/2) .* exp(1i*angle(S)); % 人声频谱
voiceY = real(ifft(voiceS,'symmetric')); % 人声时域信号

% 输出音频文件
audiowrite('voice_file.wav',voiceY,fs);

注：此代码仅供参考，具体实现可能需要根据实际情况进行调整。