武汉理工大学《信息处理课群综合训练与设计》课程设计说明书
subplot(2,4,1)
plot(source); %画出原音波形 title('原音波形')
source1=fft(source,length(source)); subplot(2,4,5);
plot(abs(source1)); %画出原音频谱 title('原音频谱') %axis([0 8000 0 80]);
[y1,Fs1,bits1]=wavread('RLSprimsp.wav'); subplot(2,4,2)
plot(y1) %画出麦克风主语音波形 title('麦克风主语音波形'); Y1=fft(y1,length(y1));
subplot(2,4,6); plot(abs(Y1)); %画出麦克风主语音频谱 title('麦克风主语音输入频谱'); %axis([0 8000 0 80]);
subplot(2,4,3); plot(fref); %画出噪声语音波形 title('噪声语音波形');
[y2,Fs2,bits2]=wavread('RLSrefns.wav'); Y2=fft(y2,length(y2));
subplot(2,4,7); plot(abs(Y2)); %画出噪声语音频谱 title('噪声语音频谱');
subplot(2,4,4); plot(output); %画降噪后的语音波形 title('降噪后的语音波形');
wavwrite(output,fs,'asd') %生成降噪后的语音 [y,Fs3,bits3]=wavread('asd.wav'); Y=fft(y,length(y)); subplot(2,4,8);
plot(abs(Y)); %画出降噪后语音频谱
11
武汉理工大学《信息处理课群综合训练与设计》课程设计说明书
title('降噪后语音输出频谱'); %axis([0 8000 0 80]); f=1; while f==1
flag=input('请选择输出语音；1：原式语音 2：加噪主语音 3：噪声语音 4：降噪后语音');
switch flag case 1
sound(source,fs,bits); case 2
sound(y1,Fs1,bits1); case 3
sound(y2,Fs2,bits2); case 4
sound(output,Fs3,bits3); end
f=input('是否重新选择？ 是 输入1 否 输入0'); end
12
武汉理工大学《信息处理课群综合训练与设计》课程设计说明书
5仿真结果与分析
5.1原始音频信号
（1）原始语音波形
图5-1原始语音波形图
（2）原始语音频谱
图5-2原始语音频谱
13
武汉理工大学《信息处理课群综合训练与设计》课程设计说明书
5.2麦克风主语音信号
（1）麦克风主语音波形
图5-3增强的语音信号波形图
（2）麦克风主语音频谱
图5-4麦克风主语音频谱
14
武汉理工大学《信息处理课群综合训练与设计》课程设计说明书
5.3噪声语音信号
（1）噪声语音波形
（2）噪声语音频谱
图5-5噪声语音波形
图5-6噪声语音频谱
15