语音增强matlab代码简单

以下是一个简单的MATLAB语音增强代码示例：

% 读取原始音频文件
[x,fs] = audioread('original_audio.wav');

% 设计带通滤波器
N = 100; % 滤波器阶数
low_freq = 1000; % 低截止频率
high_freq = 4000; % 高截止频率
b = fir1(N, [low_freq/(fs/2) high_freq/(fs/2)], 'bandpass');

% 应用滤波器
y = filter(b, 1, x);

% 增加噪声
snr = 10; % 信噪比
noise = randn(size(x)); % 随机高斯噪声
noise = noise / norm(noise) * norm(x) / 10^(snr/20); % 缩放噪声以获得所需的信噪比
y_noisy = y + noise; % 添加噪声

% 应用Wiener滤波器
winlen = 256; % 窗口长度
overlap = winlen / 2; % 窗口重叠
nfft = winlen; % FFT点数
win = hamming(winlen); % 加窗
[S,F,T] = spectrogram(y_noisy, win, overlap, nfft, fs); % 计算噪声信号的短时功率谱
noise_power = mean(abs(S).^2, 2); % 计算噪声功率谱
[S,F,T] = spectrogram(y, win, overlap, nfft, fs); % 计算干净信号的短时功率谱
signal_power = mean(abs(S).^2, 2); % 计算干净信号功率谱
SNR = signal_power ./ noise_power; % 估计信噪比
G = min(1, SNR./(SNR+1)); % 计算Wiener滤波器增益函数
y_enhanced = istft(bsxfun(@times, G, stft(y_noisy, winlen, overlap, nfft)), winlen, overlap, nfft); % 应用Wiener滤波器

% 保存增强后的音频文件
audiowrite('enhanced_audio.wav', y_enhanced, fs);

该代码从一个原始音频文件读取语音信号，并使用一个带通滤波器滤除不需要的频率成分。然后，它添加高斯噪声以模拟噪声环境，并使用一个Wiener滤波器对信号进行增强。最后，增强后的语音信号被保存到一个新的音频文件中。