成都理工大学2012届本科毕业设计（论文）
致谢
在此论文撰写过程中，要特别感谢我的导师老师的指导与督促，同时感谢她的谅解与包容。没有老师的帮助也就没有今天的这篇论文。求学历程是艰苦的，但又是快乐的。感谢我的辅导员段老师，谢谢他在这四年中为我们全班所做的一切，他不求回报，无私奉献的精神很让我感动，再次向她表示由衷的感谢。在这四年的学期中结识的各位生活和学习上的挚友让我得到了人生最大的一笔财富。在此，也对他们表示衷心感谢。
本文参考了大量的文献资料，在此，向各学术界的前辈们致敬！
28
成都理工大学2012届本科毕业设计（论文）
参考文献
[1] Boll S F.Suppression of Acoustic Noise in Speech Using Spectral Subtraction[J].IEEE Trans.on Acoustics,Speech,and Signal Processing,1979,27(2):113-120. [2] Berouti M,Schwartz R,Makhoul J.Enhancement of SpeechCorrupted by Acoustic Noise[J].IEEE Trans.on Acoustics,Speech,and Signal Processing,1979,4:208-211. [3] 胡航，语音信号处理，哈尔滨工业大学出版社，2000 年 5 月
[4]ThomsonDJ.SpectrumEstimationand armonicAnalysis[J].Proc.IEEE,1982,70(9):1 055-1 096.
[5]皇甫堪,陈建文,楼生强.现代数字信号处理[M].北京:电子工业出版社,2003.
[6]Hu Yi,Loizou P C.Speech Enhancement Based on WaveletThresholding the Multitaper Spectrum[J]. IEEE Trans.onSpeech and Audio Processing,2004,12(1):59-67. [7]吴红卫,吴镇扬,赵力.基于多窗谱的心理声学语音增强[J].声学学 报,2007,32(3):275-281. [8]潘欣裕，童兴法，赵鹤鸣，基于谱能比例加权的谱减法语音增强研究[J]. 中 国电子学会第十五届信息论学术年会暨第一届全国网络编码学术年会论文集
[9] 程正,赵鹤鸣. 基于多频带谱减法的语音增强算法的研究[J]. 002-8331(2007)36-0040-03.
[10] 曹瑜镠,方元,吕勇.基于最小统计及谱减法的语音增强[J]. 语音技术. 002-8684(2006)12-0043-04
[11］白文雅，黄健群，陈智伶.基于维纳滤波语音增强算法的改进实现[J].电声 技术，2007,31（1） ：44-46.
[12］蔡斌.一种改进型 MMSE 语音增强方法[J].信号处理，2004,20（1） ：70-74. [13］陈俊，孙洪，董航.基于 MMSE 先验信噪比估计的语音增强[J].武汉大学学 理学版，2005，51（5） ：638-642.
[14］樊昌信，曹丽娜.通信原理[M].北京：国防工业出版社，2008. [15］姚天任.数字语音处理[M].武汉：华中科技大学出版社，2005. 附录 1
Ft=8000; Fp=1000;
29
成都理工大学2012届本科毕业设计（论文）
Fs=1200; wp=2*Fp/Ft; ws=2*Fs/Ft; rp=1; rs=50;
p=1-10.^(-rp/20); %通带阻带波纹 s=10.^(-rs/20); fpts=[wp ws]; mag=[1 0]; dev=[p s];
[n21,wn21,beta,ftype]=kaiserord(fpts,mag,dev);%kaiserord求阶数截止频率 b21=fir1(n21,wn21,Kaiser(n21+1,beta)); %由fir1设计滤波器 [h,w]=freqz(b21,1); %得到频率响应 plot(w/pi,abs(h));
title('FIR低通滤波器'); 附录2
Ft=8001; Fp=4000; Fs=3500; wp=2*Fp/Ft; ws=2*Fs/Ft; rp=1; rs=50;
p=1-10.^(-rp/20); %通带阻带波纹 s=10.^(-rs/20); fpts=[ws wp]; mag=[0 1]; dev=[p s];
[n23,wn23,beta,ftype]=kaiserord(fpts,mag,dev);
b23=fir1(n23,wn23,'high',Kaiser(n23+1,beta)); %由fir1设计滤波器 [h,w]=freqz(b23,1); %得到频率响应 plot(w*12000*0.5/pi,abs(h)); title('FIR高通滤波器'); axis([3000 6000 0 1.2]); 附录3
Ft=8000; Fp=1000; Fs=1200;
wp=2*pi*Fp/Ft; ws=2*pi*Fs/Ft; fp=2*Ft*tan(wp/2); fs=2*Fs*tan(wp/2);
[n11,wn11]=buttord(wp,ws,1,50,'s'); %求低通滤波器的阶数和截止频率 [b11,a11]=butter(n11,wn11,'s'); %求S域的频率响应的参数
30
成都理工大学2012届本科毕业设计（论文）
[num11,den11]=bilinear(b11,a11,0.5); %双线性变换实现S域到Z域的变换 [h,w]=freqz(num11,den11); %根据参数求出频率响应 plot(w*8000*0.5/pi,abs(h)); title('IIR低通滤波器'); legend('用butter设计'); grid; 附录4
Ft=8000; Fp=4000; Fs=3500;
wp1=tan(pi*Fp/Ft); %高通到低通滤波器参数转换 ws1=tan(pi*Fs/Ft); wp=1;
ws=wp1*wp/ws1;
[n13,wn13]=cheb1ord(wp,ws,1,50,'s'); %求模拟的低通滤波器阶数和截止频 [b13,a13]=cheby1(n13,1,wn13,'s'); %求S域的频率响应的参数 [num,den]=lp2hp(b13,a13,wn13); %将S域低通参数转为高通的
[num13,den13]=bilinear(num,den,0.5); %利用双线性变换实现S域到Z域转 [h,w]=freqz(num13,den13); plot(w*21000*0.5/pi,abs(h)); title('IIR高通滤波器'); legend(用cheby1设计'); 附录5
[y,fs,bits]=wavread('wangqingtian.wav'); sound(y,fs) n=length(y) y_p=fft(y,n); f=fs*(0:n/2-1)/n; figure(1)
subplot(2,1,1); plot(y);
title('原始语音信号采样后时域波形'); xlabel('时间轴') ylabel('幅值 A') subplot(2,1,2);
plot(f,abs(y_p(1:n/2)));
title('原始语音信号采样后的频谱图'); xlabel('频率Hz'); ylabel('频率幅值');
noise=1*sin(2*pi*3000*n); y_z=y+noise; sound(y_z,fs) L=length(y_z); y_zp=fft(y_z,L);
31
成都理工大学2012届本科毕业设计（论文）
f=fs*(0:L/2-1)/L; figure(2)
subplot(2,1,1); plot(y_z);
title('加噪语音信号时域波形'); xlabel('时间轴') ylabel('幅值 A') subplot(2,1,2);
plot(f,abs(y_zp(1:L/2))); title('加噪语音信号频谱图'); xlabel('频率Hz'); ylabel('频率幅值');
附录6
[y,fs,bits]=wavread('wangqingtian.wav'); sound(y,fs) n=length(y) y_p=fft(y,n); f=fs*(0:n/2-1)/n; figure(1)
subplot(2,1,1); plot(y);
title('原始语音信号采样后的时域波形'); xlabel('时间轴') ylabel('幅值A') subplot(2,1,2);
plot(f,abs(y_p(1:n/2)));
title('原始语音信号采样后的频谱图'); xlabel('频率Hz'); ylabel('频率幅值'); L=length(y)
noise=0.1*randn(L,2); y_z=y+noise; sound(y_z,fs) n=length(y); y_zp=fft(y_z,n); f=fs*(0:n/2-1)/n; figure(2)
subplot(2,1,1); plot(y_z);
title('加噪语音信号时域波形'); xlabel('时间轴') ylabel('幅值A') subplot(2,1,2);
plot(f,abs(y_zp(1:n/2)));
32