matlab音频去噪代码

### MATLAB 中的音频去噪代码示例

在处理音频信号时，去除噪声是一个常见的需求。下面提供了一个基于频域滤波器设计的方法来实现简单的音频去噪功能。

#### 使用傅里叶变换进行简单谱减法降噪

通过快速傅立叶变换 (FFT)，可以将时间序列转换到频率空间，在这里更容易识别并抑制不需要的成分：

% 加载含噪语音文件
[x, fs] = audioread('noisy_speech.wav');

N = length(x);
nfft = 2^nextpow2(N); % 计算最接近 N 的2次幂作为 FFT 长度

X = fft(x,nfft)/N;    % 对输入数据做归一化后的 FFT 变换
f = fs/2*linspace(0,1,nfft/2+1);

% 设定阈值以区分有用信号和背景噪音
threshold = mean(abs(X)); 

for k=1:nfft
    if abs(X(k)) < threshold * 0.8
        X(k)=0;
    end
end

y = real(ifft(X)*N);     % 进行 IFFT 并取实部恢复原始长度的时间序列
sound(y,fs)

audiowrite('denoised_audio.wav', y, fs);

此脚本读入一段有噪声污染的声音片段 `noisy_speech.wav` ，对其进行离散傅里叶变换，并设置一个幅度阈值来判断哪些部分属于噪声从而将其置零。最后再执行逆向变换得到净化过的新声音文件保存下来[^1]。

对于更复杂的场景，则可能需要用到自适应滤波技术或者其他高级算法比如小波变换、卡尔曼滤波等来进行更加精细有效的处理[^2]。