能否提供一段用于音频信号谱减法处理的MATLAB代码示例？

当然可以。音频信号谱减法通常涉及傅里叶变换将音频信号从时域转换到频域，然后对两个信号的频谱进行操作，最后通过反傅立叶变换返回原始信号。这里是一个简单的MATLAB代码示例，假设我们有两个单声道音频文件（如声源和背景噪音）：

% 加载音频文件（假设src是声源文件，bgd是背景噪音文件）
[src, Fs] = audioread('src.wav');
[bgd, Fs_bgd] = audioread('bgd.wav'); % 确保两个文件采样率一致

% 检查采样率是否匹配
ifFs ~= Fs_bgd
    error('Audio files have different sample rates.');
end

% 投影噪声到源频率域
noise_spec = abs(fft(bgd)); % 取绝对值表示实数
noise_spec = noise_spec./sum(noise_spec); % 归一化

% 对声源信号进行谱减法
src_spec = abs(fft(src));
reduced_spec = src_spec .* (1 - noise_spec);

% 反傅立叶变换回时域
reduced_audio = ifft(reduced_spec);
reduced_audio = real(reduced_audio); % 只保留实部，因为反FFT得到的是复数

% 保存处理后的音频
audiowrite('reduced_audio.wav', reduced_audio, Fs);

这个例子假设背景噪音的幅度在整个频率范围内均匀分布，实际上可能会更复杂，比如需要考虑相位信息。此外，这仅是一个基础的谱减处理，实际应用中可能还需要进一步优化和调整。