matlab维纳滤波语音去噪代码

以下是使用matlab实现维纳滤波语音去噪的代码：

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

% 设置参数
win = 256; % 窗口大小
hop = win/2; % 帧移
nfft = win; % FFT点数
win_type = hamming(win); % 窗口类型
alpha = 0.99; % 平滑系数
SNR = 10; % 信噪比

% 计算STFT
[S, f, t] = stft(x, win, hop, nfft, fs, 'Window', win_type);

% 计算信噪比
noise = x - S;
noise_power = mean(abs(noise(:)).^2);
signal_power = mean(abs(S(:)).^2);
SNR = 10*log10(signal_power/noise_power);

% 计算噪声功率谱密度
N = size(S, 2);
Nk = round(N/2);
noise_power_spectrum = mean(abs(S(:, Nk+1:end)).^2, 2);

% 计算信号功率谱密度
signal_power_spectrum = mean(abs(S(:, 1:Nk)).^2, 2);

% 计算维纳滤波器
G = signal_power_spectrum./(signal_power_spectrum + noise_power_spectrum/SNR);

% 进行维纳滤波
Y = S.*G;

% 计算逆STFT
y = istft(Y, win, hop, nfft, fs, 'Window', win_type);

% 平滑处理
y = filter(1-alpha, [1, -alpha], y);

% 输出结果
audiowrite('clean_speech.wav', y, fs);