主动降噪matlab代码

主动降噪是音频信号处理中一种常用的方法，可以有效地去除噪声，提取出清晰的信号。以下是一个简单的主动降噪的Matlab代码示例。

【首先读取音频文件】
[y, fs] = audioread('input.wav'); % 读取音频文件，得到音频信号y和采样率fs

【预处理】
frameSize = 512; % 设置帧大小
overlap = 256; % 设置帧重叠长度
window = hamming(frameSize); % 使用汉明窗进行加窗

【计算帧数】
totalFrame = floor((length(y) - overlap) / (frameSize - overlap)); % 计算总帧数

【建立Noise Estimate】
noiseEstimate = zeros(frameSize, 1); % 初始化噪声估计
alpha = 0.98; % 设置估计系数
for i = 1:totalFrame
frame = y((i-1)*(frameSize-overlap)+1 : (i-1)*(frameSize-overlap)+frameSize); % 提取当前帧
noiseEstimate = alpha * noiseEstimate + (1 - alpha) * frame.^2; % 更新噪声估计
end

【构建Spectral Subtraction】
SNR = 15; % 设置信噪比
gain = 1 - (1 ./ (1 + (10^(-SNR/10)) .* (noiseEstimate ./ (abs(fft(frame, frameSize)) .^ 2)))); % 计算增益

【主动降噪】
output = zeros(length(y), 1); % 初始化输出信号

for i = 1:totalFrame
frame = y((i-1)*(frameSize-overlap)+1 : (i-1)*(frameSize-overlap)+frameSize); % 提取当前帧
frameWindowed = frame .* window; % 加窗

frameFFT = fft(frameWindowed, frameSize); % 傅里叶变换
frameGain = abs(frameFFT) .* gain; % 根据增益调整频谱
frameProcessed = real(ifft(frameGain .* exp(1j * angle(frameFFT)), frameSize)); % 逆傅里叶变换

output((i-1)*(frameSize-overlap)+1 : (i-1)*(frameSize-overlap)+frameSize) = ...
output((i-1)*(frameSize-overlap)+1 : (i-1)*(frameSize-overlap)+frameSize) + frameProcessed; % 合成输出信号
end

【播放输出信号】
sound(output, fs); % 播放输出信号

该代码实现了主动降噪的基本步骤，包括帧的预处理、噪声估计、信号增益计算和主动降噪处理。根据输入音频信号的特征，可以调整参数和增加其他处理步骤来进一步优化降噪效果。