基于低通滤波器语音信号加噪与去噪含matlab源码

下面是基于低通滤波器的语音信号加噪与去噪的Matlab源码。

加噪部分：

% 读取语音信号
[y, fs] = audioread('speech_signal.wav');

% 生成噪声信号
noise = randn(size(y));

% 设置噪声信噪比（SNR）为10dB
snr = 10;

% 计算噪声功率
noise_power = var(noise);

% 计算信号功率
signal_power = var(y);

% 根据信噪比计算噪声信号的标准差
noise_std = sqrt(signal_power/(10^(snr/10)));

% 缩放噪声信号的标准差
noise = noise_std * noise;

% 添加噪声到语音信号
y_noisy = y + noise;

% 播放加噪后的信号
sound(y_noisy, fs);

% 保存加噪后的信号
audiowrite('speech_signal_noisy.wav', y_noisy, fs);

去噪部分：

% 读取加噪信号
[y_noisy, fs] = audioread('speech_signal_noisy.wav');

% 应用低通滤波器去噪
fc = 1000; % 设置截止频率为1kHz
[b, a] = butter(6, fc/(fs/2)); % 6阶低通滤波器设计
y_denoised = filter(b, a, y_noisy);

% 播放去噪后的信号
sound(y_denoised, fs);

% 保存去噪后的信号
audiowrite('speech_signal_denoised.wav', y_denoised, fs);

以上源码中，首先通过读取原始的语音信号，并生成与其同样长度的随机噪声信号。然后根据设定的信噪比计算噪声信号的标准差，并将其缩放为具有相同标准差的加性噪声。将噪声信号与原始语音信号相加，生成加噪后的语音信号，并保存为`speech_signal_noisy.wav`文件。

在去噪部分，首先读取加噪后的语音信号，然后使用`butter`函数设计一个6阶低通滤波器，该滤波器的截止频率为1kHz。最后，将加噪信号通过低通滤波器进行去噪处理，得到去噪后的语音信号，并将其保存为`speech_signal_denoised.wav`文件。

请注意，以上仅提供了一种简单的基于低通滤波器的语音信号加噪与去噪方法，实际情况可能需要根据不同的应用场景和要求进行适当调整。