matlab 实现语音增强（含维纳滤波、谱减法、卡尔曼滤波三种实现方法）

%%三种语音增强方法的测试脚本 %****************************************************** % 在audioread函数中可以设置读入的语音信号 % 改变SNR的值即可改变加入的噪声 % [Input, Fs] = audioread('sp01.wav'); Time = (0:1/Fs:(length(Input)-1)/Fs)'; %取单声道 Input = Input(:,1); %SNR为加入噪声与纯净信号的信噪比(dB) SNR=10; [NoisyInput,Noise] = add_noise(Input,SNR);%NoisyInput为加噪信号，Noise是噪声 %% 三种语音增强方法的实现 [spectruesub_enspeech] = spectruesub(NoisyInput); [wiener_enspeech] = wienerfilter(NoisyInput); [Klaman_Output] = kalman(NoisyInput,Fs,Noise); %% spectruesub绘制 %将信号长度对齐 sig_len=length(spectruesub_enspeech); NoisyInput=NoisyInput(1:sig_len); Input=Input(1:sig_len); wiener_enspeech=wiener_enspeech(1:sig_len); Klaman_Output=Klaman_Output(1:sig_len); Time = (0:1/Fs:(sig_len-1)/Fs)'; % Time= ((0:1/Fs:(sig_len)-1)/Fs)'; figure(1) MAX_Am(1)=max(Input); MAX_Am(2)=max(NoisyInput); MAX_Am(3)=max(spectruesub_enspeech); subplot(3,1,1); plot(Time, Input) ylim([-max(MAX_Am),max(MAX_Am)]); xlabel('Time') ylabel('Amlitude') title('原始信号') subplot(3,1,2); plot(Time, NoisyInput) ylim([-max(MAX_Am),max(MAX_Am)]); xlabel('Time') ylabel('Amlitude') title('加噪信号') subplot(3,1,3); plot(Time, spectruesub_enspeech) ylim([-max(MAX_Am),max(MAX_Am)]); xlabel('Time') ylabel('Amlitude') title('谱减法') %% spectruesub绘制 % Time_wiener = (0:1/Fs:(length(wiener_enspeech)-1)/Fs)'; figure(2) MAX_Am(1)=max(Input); MAX_Am(2)=max(NoisyInput); MAX_Am(3)=max(wiener_enspeech); subplot(3,1,1); plot(Time, Input) ylim([-max(MAX_Am),max(MAX_Am)]); xlabel('Time') ylabel('Amlitude') title('原始信号') subplot(3,1,2); plot(Time, NoisyInput) ylim([-max(MAX_Am),max(MAX_Am)]); xlabel('Time') ylabel('Amlitude') title('加噪信号') subplot(3,1,3); plot(Time, wiener_enspeech) ylim([-max(MAX_Am),max(MAX_Am)]); xlabel('Time') ylabel('Amlitude') title('谱减法') %% Kalman绘制 figure(3) MAX_Am(1)=max(Input); MAX_Am(2)=max(NoisyInput); MAX_Am(3)=max(Klaman_Output); subplot(3,1,1); plot(Time, Input) ylim([-max(MAX_Am),max(MAX_Am)]); xlabel('Time') ylabel('Amlitude') title('原始信号') subplot(3,1,2); plot(Time, NoisyInput) ylim([-max(MAX_Am),max(MAX_Am)]); xlabel('Time') ylabel('Amlitude') title('加噪信号') subplot(3,1,3); plot(Time, Klaman_Output) ylim([-max(MAX_Am),max(MAX_Am)]); xlabel('Time') ylabel('Amlitude') title('Kalman滤波') %% 求语音降噪后的信噪比 SNR(1)=snr(Input,Input-spectruesub_enspeech); SNR(2)=snr(Input,Input-wiener_enspeech); SNR(3)=snr(Input,Input-Klaman_Output);