IIR带通滤波器语音去噪要点.docx

1、IIR带通滤波器语音去噪要点摘要语音信号滤波处理时研究用数字信号处理技术和语音学知识对语音信号进 行处理的新兴学科，是目前发展最为迅速的信息科学研究领域的核心技术之一， 通过语音传递信息是人类交流信息最自然、 最有效、 最方便的手段。 本次主要通 过录制一段语音，对其进行时域、频谱分析，并利用 matlab 的信号处理工具箱 对语音进行加噪然后再用 IIR 数字带通滤波器滤除噪声， 最后对比滤波前后的语 音信号的时域、频域特性。关键字：IIR ;双线性变换；模拟低通滤波器；切比雪夫I; MATLAB目录 2前言 1一设计原理 21.数字滤波器简介 22.IIR 数字滤波器的设计原理 33.II

2、R 滤波器的特点 3二 IIR 数字滤波器的设计方法 41. 模拟滤波器 42.双线性变换法 6三 IIR 数字滤波器设计过程 91. 设计步骤 92. 音频信号部分程序 103.程序流程图 104 .仿真结果 11总结 13致谢 14参考文献 15附录： 16、八 、-前言通过语音传递信息室人类最重要、 最有效、最常用和最方便的交换信息的形 式。语音是人类特有的功能， 声音是人类最常用的工具， 是相传递信息的最重要 的手段。因此，语音信号是人类构成思想疏通和感情交流的最重要的途径之一。 并且，由于语言和语音与人的智力活动密切相关，与社会文化和进步紧密相连， 所以它具有最大的信息容量和最高的智

3、能水平。 现在，人类已开始进入了信息化 时代，用现代手段研究语音信号，使人们能够更加有效地产生、传输、存储、获 取和应用信息， 这对于促进社会的发展具有十分重要的意义。 让计算机能听懂人 类的语言，是人类自计算机诞生以来梦寐以求的想法。随着计算机越来越向便携化方向发展， 随着计算环境的日趋复杂化， 人类越 来越迫切要求拜托键盘的束缚而带至以语音输入这样便于使用的、 自然的、人性 化的输入方法。 作为高科技应用领域的热点， 语音信号的采集和分析从理论的研 究到产品的开发已走过了几十个春秋并且取得了长足的进步。它正在直接与办 公、交通、金融、公安、商业、旅游等行业的语音咨询与管理，工业生产部门的

4、语音声控， 电话、电信系统的自动拨号、 辅助控制和查询以及医疗卫生和福利事 业的生活支援系统等各种实际应用领域相接轨， 并且有望成为下一代操作系统的 和应用程序的用户界面。 可见，语音信号采集与分析的研究将是一项极具市场价 值和挑战性的工作。 我们今天进行这一领域的研究与开发就是要让语音信号处理 技术走向人们的日常生活当中，并不断朝更高目标而努力。语音信号采集与分析之所以能够那样长期地、 深深地吸引广大科学工作者对 其进行研究和探讨， 除了它的实用性之外， 另一个重要原因是， 它始终与当时信 息科学中最活跃的前沿的学科保持密切连续， 并且一起发展。 语音信号的采集和 分析是以语音语言学和数字信

5、号处理为基础而形成的一门设计面很广的综合性 学科，与心理、生理学、计算机科学、通信与信息科学以及模式识别和人工智能 等学科都有着非常密切的关系。 对语音信号的采集和分析的研究一直是数字信号 处理技术发展的重要推动力量。 因为许多处理的新方法的提出， 首先是语音信号 处理中获得成功，然后在推广到其他领域。一.设计原理IIR数字滤波器的设计一般是利用目前已经很成熟的模拟滤波器的设计方法 来进行设计，通常采用的模拟滤波器原型有 butterworth 函数、Chebyshev函数、 Bessel函数、椭圆滤波器函数等。1.数字滤波器简介数字滤波器是一种用来过滤时间离散信号的数字系统，通过对抽样数据进

6、行 数学处理来达到频域滤波的目的。 可以设计系统的频域响应，让它满足一定的要 求,从而对通过该系统的信号的某些特定的频率成分进行过滤， 这就是滤波器的 基本原理。如果系统是一个连续系统，则滤波器称为模拟滤波器。如果系统是一 个离散系统，则滤波器称为数字滤波器。数字滤波器的技术指标如图 1-1所示:s 0p:通带截止频率“ s：阻带波动 “ P :通带波动通带衰减（dB）: Ap = -20 lg（ 1 - 云 p ）阻带衰减（dB）: As = -20 lg s信号通过线性系统后，其输出就是输入信号和系统冲击响应的卷积。 除此之 外，输出信号的波形将不同于输入信号的波形。 从频域分析来看，信号

7、通过线性 系统后，输出信号的频谱将是输入信号的频谱与系统传递函数的乘积。 除非系统 函数是常数，否则，输出信号的频谱将不同于输入信号的频谱。2.IIR数字滤波器的设计原理（1） 按照一定规则把给定的滤波器技术指标转换为模拟低通滤波器的技术指标；（2） 根据模拟滤波器的技术指标设计为相应的模拟低通滤波器；（3） 根据双线性不变法把模拟滤波器转换为数字滤波器；（4） 根据要设计的带通滤波器，首先把它们的技术指标转化为模拟低通滤波器 的技术指标，设计为相应的模拟低通滤波器，最后通过频率转化的方法来得到所 要的滤波器如下图1-2所示：图1-2 IIR 数字滤波器的设计原理3.IIR滤波器的特点（1）

8、IIR数字滤波器的系统函数可以写成封闭函数的形式。（2） IIR数字滤波器采用递归型结构，及结构上带有反馈环路。IIR滤波器运算 结构通常由延时、乘以系数和相加的基本运算组成，可以组合成直接型、并联型、 级联型、正准型四种结构形式，都具有反馈回路，由于运算中的舍入处理，使误 差不断累积，有时会产生微弱的几声震荡。（3） IIR数字滤波器在设计上可以借助成熟的模拟滤波器的成果，如：巴特沃 斯、切比雪夫和椭圆滤波器等，有现成的设计数据和图表可查，其设计工作量比 较小，对计算工具的要求不高。在设计一个IIR数字滤波器时，我们根据指标先 写出模拟滤波器的公式，然后通过一定的变化，将模拟滤波器的公式转化

9、为数字 滤波器的公式。（4） IIR数字滤波器的相位特性不好控制，对相位要求较高时，需加相位校准 网络。二. IIR数字滤波器的设计方法IIR数字滤波器是一种离散时间系统，其系统函数为公式（2-1 ）:M M、bkZ 上 IT （z-Zr）H（Z）二占 k罟 （2-1）_k1 akZ ii （z-Zp）k=1 k=1（1） 当 M=N,H（z）:N 阶 IIR 系统 +（ M-N）阶 FIR 系统。（2） 以上两种表示等价，部分分式形式和零极点增益形式。（3） iir系统的逼近，就是找到滤波器的系数 ak, bk，或者是系统的零极点和 增益（乙p,k ）。1.模拟滤波器（一）技术指标如图2-1

10、所示：图2-1模拟滤波器的技术指标p :通带波动通带衰减（dB）: Ap 二-20lg（1 -）阻带衰减（dB）:滤波器的衰减函数:As 20 lg、sA w = -20lg jw dB（二）切比雪夫I型的模拟低通滤波器的幅度响应如下图 2-2所示:A H （ j 图2-2 切比雪夫I型的模拟低通滤波器的幅度响应幅度响应模方为式（2-2）：1 （2-2）-值越大通i ；2cN( / 匕)式中wc为有效通带截止频率，I二表示与通带波纹有关的参量, 带不动愈大，N表示滤波器阶数。Cn（x）是N阶切比雪夫多项式，定义为（2-3）所示：(2-3)cosN arccos(x) coshN arccos

11、h(x)L（三）切比雪夫低通滤波器的频域特性如式（2-4）:、2 1（1） 0 _w _Wc时，-jw 在俐 2间振荡1 +呂2（2） wwc时,H（jw由单调下降（N增大，下降加速）2 1（3） IIjWc 2控制了通带衰减N为奇数时,|H( jO=12 iN为偶数时,|H( jO|j =21（四）切比雪夫低通滤波器的设计步骤（1） 由通带截频Wp确定c Wc=Wp（2） 由通带衰减确定； ；=.100.1Ap-1（3） 由通带、阻带指标确定Narccosh(丄 j100.1As -1)2.双线性变换法（一）基本思想将非带限的模拟滤波器映射为最高频率的带限模拟滤波器7）：2w arctan(

12、w-wT ) 飞_oo +旳 -n /T 兀 /T】（二） s域到z域的映射关系如式（2-810）:wS(2-4)(2-5),其方法为式(2-6(2-6)(2-7)(2-8)式中,T仍是采样间隔=j-ta n匚)T 2Q2 sin(巧) =i 2j t 111 cos()2(2-9)T 1 e-j,?s. .jw、z-e J2 1 - z-1S = Tl zJ s-s(2-10)（三）稳定性分析2 1 -zT 1 z，2 T + s二 z =2/T - s令s - ；- jw，则有i 2 2(2/T Y)2 Z2(2-11)(1) CT 0, z 1(2) r=0, z =1(3) 0, z

13、1(2/T -二)2S域左半平面映射到z域单位圆内S域虚轴映射到z 域单位圆上S域右半平面映射到z域单位圆外因此，稳定的模拟滤波器经双线性变换后所得的数字滤波器也一定是稳定的。（四）双线性变换法优缺点双线性变换法与脉冲响应不变法相比，其主要的优点是避免了频率响应的混叠 现象。缺点是幅度响应不是常数时会产生幅度失真。（五）Q和3的关系如下图2-3所示：w = 2tan =T 12丿图2-3 模拟频率3与数字频率Q的对应关系由图2-3看出，在零频率附近，模拟角频率Q与数字频率3之间的变换关系 接近于线性关系；但当Q进一步增加时，3增长得越来越慢，最后当时,3终止在折叠频率3 =冗处，因而双线性变换

14、就不会出现由于高频部分超过折叠 频率而混淆到低频部分去的现象，从而消除了频率混叠现象。但是双线性变换的这个特点是靠频率的严重非线性关系而得到的，如式（2-10）及图2-3所示。由于这种频率之间的非线性变换关系， 就产生了新的问 题。首先，一个线性相位的模拟滤波器经双线性变换后得到非线性相位的数字滤 波器，不再保持原有的线性相位了；其次，这种非线性关系要求模拟滤波器的幅 频响应必须是分段常数型的，即某一频率段的幅频响应近似等于某一常数 （这正 是一般典型的低通、高通、带通、带阻型滤波器的响应特性） ，不然变换所产生的数字滤波器幅频响应相对于原模拟滤波器的幅频响应会有畸变，如图 2-4所示。图2-

15、4 双线性变换法幅度和相位特性的非线性映射对于分段常数的滤波器，双线性变换后，仍得到幅频特性为分段常数的滤波器，但是各个分段边缘的临界频率点产生了畸变， 这种频率的畸变，可以通过频率的预畸来加以校正。也就是将临界模拟频率事先加以畸变， 然后经变换后正好 映射到所需要的数字频率上。三 IIR 数字滤波器设计过程1. 设计步骤根据以上 IIR 数字滤波器设计方法，下面运用双线性变换法基于 MATLAB 设计一个 IIR 带通滤波器。(1)确定性能指标在设计带通滤波器之前 , 首先根据工程实际的需要确定滤波器的技术指标：通带截止频率fp仁1200 , fp2=3000;阻带截止频率fs仁1000 ;

16、 fs2=3200;阻带最小衰减As=100d环口通带最大衰减Ap=1dB;(2)把频率转化为数字角频率wp1=2*pi*fp1*T; wp2=2*pi*fp2*T;ws1=2*pi*fs1*T; ws2=2*pi*fs2*T;(3)频率预畸变用Q =2/T*tan(w/2)对带通数字滤波器 H(z)的数字边界频率预畸变，得 到带通模拟滤波器H(s)的边界频率主要是通带截止频率 Wp1,Wp2阻带截止 频率Ws1 Ws2的转换。抽样频率fs=10KHz。通带截止频率 Wp1=(2/T)*tan(wp1/2) ; Wp2=(2/T)*tan(wp2/2) ;阻带截止频率 Ws1=(2/T)*ta

17、n(ws1/2) ; Ws2=(2/T)*tan(ws2/2)(4)模拟带通性能指标转换成模拟低通性能指标WP=Wp1,Wp2;W0=sqrt(Wp1*Wp2);B=Wp2-Wp1; WS=Ws1,Ws2;(5)模拟低通滤波器的构造借助切比雪夫 (Chebyshev) 滤波器得到模拟低通滤波器的传输函数Ha(s)。N,Wc=cheb1ord(WP,WS,Ap,As,s); % 求阶数和边缘频率 z0,p0,k0=cheb1ap(N,Ap); % 求极点,零点和增益num=k0*real(poly(z0); % 模拟低通滤波器系统函数的分子多项式 den=real(poly(p0);(6)模拟低

18、通滤波器转换成模拟带通滤波器调用 lp2bp 函数将模拟低通滤波器转化为模拟带通滤波器。numt,dent=lp2bp(num,den,W0,B);%模拟带通滤波器系统函数的分子多项式和分母多项式(7)Ha(s) 转换成数字带通滤波器模拟带通滤波器转换成数字带通滤波器 利用双线性变换法将模拟带通滤波器H(z)。numd,dend=bilinear(numt,dent,fs);2.音频信号部分程序1)写入声音信号y,fs,nbits=wavread (1.wav);sound(y,fs,nbits);%回放语音信号N= length (y) ;%求出语音信号长度Y=fft(y,N);%傅里叶变换

19、产生噪声并加到声音中noise=0.01*randn(N,1);%随机函数产生噪声S=y+noise;%语音信号加入噪声sound(S);Si=fft(S);%滤波前傅里叶变换y1=filter(numd,dend,S);sound(y1);y2=fft(y1);3.程序流程图首先确定性能指标， 把频率转化为数字角频率， 进而在进行频率预畸变， 用Q =2/T*tan(w/2)对带通数字滤波器 H(z)的数字边界频率预畸变,得到带通模拟 滤波器H(s)的边界频率主要是通带截止频率 Wp1,Wp2阻带截止频率 Ws1 Ws2的转换。抽样频率 fs=10KHz。上述准备工作做好之后，就先把模拟带通

20、性能指标转换成模拟低通性能指标，然后设计模拟低通滤波器，借助切比雪夫 (Chebyshev)滤波器得到模拟低通 滤波器的传输函数 Ha(s) 。然后调用 lp2bp 函数将模拟低通滤波器转化为模拟带 通滤波器。最后利用双线性变换法将模拟带通滤波器 Ha(s)转换成数字带通滤波 器 H(z) 。；流程图如下图 3-1 所示：开始在Window下录制语音将语音格式改为 wav对语音信号进行频谱分析，画出时域和频域波形图加入噪声用切比雪夫设计IIR带通滤波器画出其频率响应用IIR滤波器对语音信号进行滤波画出语音信号滤波前后的波形并进行比较结束图3-1 IIR 带通滤波器去噪流程4 .仿真结果源程序设

21、计了模拟低通滤波器、模拟带通滤波器与数字带通滤波器，对数字带通滤波器的性能仿真如图 3-2 :图3-2 滤波器的辐频相频特性曲线为了实现IIR数字带通滤波器的应用，程序中加入了有噪声的音频信号， 通过对其滤波处理，来显示数字带通滤波器的功能， 下面显示未加入噪声如图 3-3所示频谱及加入噪声之后和加噪后经滤除噪声的频谱和波形如图 3-4所示:0.50-0.5-1滤波前信号频谱滤波前信号波形1滤波后信号频谱0 2 4 64x 10滤波后信号波形图3-3语音信号的波形及频谱图3-4 滤波前后的频谱与噪声总结在现代通信系统中， 由于信号中经常混有各种噪声和干扰， 所以很多信号分 析都是基于滤波器而进

22、行的， 而数字滤波器一定运算关系改变输入输出信号所含 频率成分的相对比例或者滤出某些频率成分的器件， 具有处理精度高、 稳定、灵 活、不从在阻抗匹配问题等优点， 可以实现模拟滤波器无法实现的特殊滤波功能。 本次试验的主要利用 MATLAB件实现基于语音信号去噪处理的IIR滤波器的设 计。此次过程中主要有： 一、根据给定的指标设计得到 IIR 数字带通滤波器； 二、 采集语音信号并对其进行时域、 频域分析； 三、给采集好的信号叠加单频正弦噪 声；四、将叠加噪声之后的语音信号通过 IIR 带通滤波器得到去噪语音信号， 并 将去噪前后的语音信号进行时域、频域的分析对比，故而可以的到实验结果。 最后，

23、虽然基于语音信号去噪处理的 IIR 滤波器的课程设计得以完成， 但是也存 在一些问题， 如选用不同的语音信号时， 经滤波后的不到与原始信号类似的信号 频谱，有事甚至什么语音信号都没有。因此，在以后的学习中要多加注意，并尽 可能找到解决的方法。致谢通过基于语音信号去噪处理的 IIR 滤波器的设计我学习掌握了好多知识， 首 先是对MATLA有了一个全新的认识，其次是对MATLAB勺更多操作和命令的使用 有了更高的掌握，最重要的是对 MATLA处理数字信号处理的相关能力有了更高 的飞跃，就对MATLA中涉及到数字信号处理的相关命令来说，通过这次课程设 计的亲身操作和实践，学习掌握了较多的原来不知道的

24、或不熟悉的命令。在此次基于语音信号去噪处理的 IIR 带通滤波器设计仿真过程中， 滤波器的 设计对于这次课设的我来说完全是实践中的新知识，虽然网上有很多关于 IIR带通滤波器的完整设计， 能有借鉴的模板和方法固然很令人放松， 但自己觉得这 样做有点浪费青春， 很没意思。 因为我借鉴过来后， 自己对基于语音信号去噪处 理的 IIR 带通滤波器设计既没有记忆， 更谈不上理解。所以一个人暗暗下定决心， 要把这块硬骨头啃下来。 因此图书馆便是不可少去的地方， 翻阅多本参考书， 进 行多重比较， 才在心中有了设计的基本概念。 在多次爬山涉水中， 才明白怎样进 行完整的设计，怎样用MATLA实现相关的仿真

25、。最后，非常感谢陈海燕老师和同学在这次基于语音信号去噪处理的 IIR 滤波 器课程设计中对我的帮助， 因为你们的不懈努力和耐心教导， 才使在这次课程设 计中收获颇丰。参考文献1陈后金.数字信号处理 M. 北京：高等教育出版社 ,2005.12程佩青.数字信号处理 M. 北京：清华大学出版社 ,2007.23从玉良.数字信号处理原理及其 MATLAB实现M.北京：电子工业出版社.2009.74胡广书. 数字信号处理理论、 算法与实现 M. 北京：清华大学出版社 .2003.85丁玉美，高西全 . 数字信号处理 M. 西安电子科技大学出版社 .20016刘泉. 数字信号处理原理与实现 M. 电子工

26、业出版社（第 2 版） .20097陈怀琛.数字信号处理教程MATLAB?义与实现M.电子工业出版社.20048张圣勤.MATLAB7.0实用教程M.北京：机械工业出版社.2006.模拟滤波器的通带截止频率模拟滤波器的阻带截止频率 带通滤波器的通带宽度 带通滤波器的中心频率求阶数和边缘频率 求极点，零点和增益 模拟低通滤波器系统函数的 分子多项式 模拟低通滤波器系统函数的分母 多项式模拟带通滤波器系统函数的 分子多项式和分母多项式 双线性变换 , 有模拟滤波器 转为数字滤波器数字带通滤波器的幅频特性附录：%技术指标fp1=1200; fp2=3000;fs1=1000; fs2=3200;Ap

27、=1; As=100;fs=10000; T=1/fs;wp1=2*pi*fp1*T; wp2=2*pi*fp2*T; ws1=2*pi*fs1*T ； ws2=2*pi*fs2*T; %带通到低通的频率转换Wp1=(2/T)*tan(wp1/2);Wp2=(2/T)*tan(wp2/2);WP=Wp1,Wp2; %Ws1=(2/T)*tan(ws1/2);Ws2=(2/T)*tan(ws2/2);WS=Ws1,Ws2; %B=Wp2-Wp1; %W0=sqrt(Wp1*Wp2); % %切比雪夫模拟低通原型滤波器的设计 N,Wc=cheb1ord(WP,WS,Ap,As,s); % z0,p

28、0,k0=cheb1ap(N,Ap); % num=k0*real(poly(z0); %den=real(poly(p0); %模拟低通原型滤波器转换为模拟带通滤波器 numt,dent=lp2bp(num,den,W0,B); %numd,dend=bilinear(numt,dent,fs); %h,w=freqz(numd,dend); % dbH=20*log10(abs(h)+eps)/max(abs(h);%以分贝表示带通滤波器的幅频特性figure(1);subplot(2,1,1);plot(w/pi,db(h); grid on; xlabel(w(rad);ylabel(

29、|H(z)|); axis(0,1,-100,25); title( 幅频特性曲线 (db);subplot(2,1,2);plot(w/pi,angle(h);grid on;xlabel(w(rad); ylabel(H(z); axis(0,1,-5,5);title( 相频特性曲线 ); y,fs,nbits=wavread ( sound(y,fs,nbits);N= length (y) ;Y=fft(y,N);figure(2);subplot(2,1,1);plot(y);title( subplot(2,1,2);plot(abs(Y); title( 原始信号频谱 ); t

30、=0:1/fs :1;noise=0.5*sin(2000*pi*t);noise=noise(:,1);S=y+noise;sound(S);Si=fft(S);y1=filter(numd,dend,S);sound(y1);y2=fft(y1);figure(3); subplot(2,2,1); plot(abs(Si),g);title( 滤波前信号频谱 ); subplot(2,2,2);公鸡打鸣 .wav);%回放语音信号 求出语音信号长度 傅里叶变换原始信号波形 );%随机函数产生噪声语音信号加入噪声滤波前傅里叶变换滤波后傅里叶变换plot(abs(y2),r);title( 滤波后信号频谱 ); subplot(2,2,3);plot(S);title( 滤波前信号波形 );subplo