基于雙線性變換法的IIR數(shù)字高通濾波器設(shè)計(jì)

1、基于雙線性變換法的IIR數(shù)字高通濾波器設(shè)計(jì)基于雙線性變換法的IIR數(shù)字高通濾波器設(shè)計(jì)摘 要隨著信息時(shí)代和數(shù)字世界的到來(lái)，數(shù)字信號(hào)處理已成為當(dāng)今一門(mén)極其重要的學(xué)科和技術(shù)領(lǐng)域。在數(shù)字信號(hào)處理中起著重要的作用并已獲得廣泛應(yīng)用的是數(shù)字濾波器（DF，Digital Filter）。數(shù)字濾波器是一種用來(lái)過(guò)濾時(shí)間離散信號(hào)的數(shù)字系統(tǒng)，通過(guò)對(duì)抽樣數(shù)據(jù)進(jìn)行數(shù)學(xué)處理來(lái)達(dá)到頻域?yàn)V波的目的。實(shí)現(xiàn)IIR濾波器的階次較低，所用的存儲(chǔ)單元較少，效率高，精度高，而且能夠保留一些模擬濾波器的優(yōu)良特性，因此應(yīng)用很廣。Matlab軟件以矩陣運(yùn)算為基礎(chǔ)，把計(jì)算、可視化及程序設(shè)計(jì)有機(jī)融合到交互式工作環(huán)境中，并且為數(shù)字濾波的研究和應(yīng)用提供

2、了一個(gè)直觀、高效、便捷的利器。尤其是Matlab中的信號(hào)處理工具箱使各個(gè)領(lǐng)域的研究人員可以直觀方便地進(jìn)行科學(xué)研究與工程應(yīng)用。本次主要利用matlab的信號(hào)處理工具箱設(shè)計(jì)一個(gè)數(shù)字高通IIR濾波器，并用此濾波器處理一段音頻信號(hào)。此報(bào)告重點(diǎn)介紹了用雙線性不變法設(shè)計(jì)IIR數(shù)字濾波器的基本流程，比較了各種設(shè)計(jì)方法的優(yōu)缺點(diǎn)，總結(jié)了模擬濾波器的性能特征。最后以雙線性不變法設(shè)計(jì)了一個(gè)高通巴特沃斯IIR數(shù)字濾波器，介紹了設(shè)計(jì)步驟，然后在Matlab環(huán)境下進(jìn)行了仿真與調(diào)試，實(shí)現(xiàn)了設(shè)計(jì)目標(biāo)。關(guān)鍵字 IIR數(shù)字濾波器 巴特沃斯 高通 雙線性變換 MATLABIAbstractWith the information

3、era and the advent of the digital world, digitalsignal processing has become one of the important fields of science and technology. In the digital signal processing plays an important role and has been widely applied in the digital filter. Digital fi

4、lter is a discrete time signal is used to filter the digital system,through the mathematical processing of sampling data to reach the goal of frequency domain filterin. Implementation of IIR filter order is low, less storage unit, 

5、;with high efficiency, high precision, and can keep some simulation characteristics of filter,  so a very wide application. Matlab software based on matrix computation, the calculation,visualization and program design of organic integration to interactive environme

6、nt, and provides an intuitive, efficient, convenient tool for the research and application of digital filter.The main use of signal processing toolbox of MATLAB to design a digital high-pass IIR filter, and an audio signal

7、60;with the filter. This report focuses on the basic process of using the bilinearinvariant IIR digital filter design, compares the advantages and disadvantages of various design methods, summarizes thecharacteristics of the analog filter. Finally,&#

8、160;by bilinear invariant to design a high Butterworth IIR digitalfilter, introduces the design steps, and then in the Matlab environment for the simulation and debugging, achieves the design goal.Keywords  IIR  Butterworth high-pass&#

9、160;bilinear transform  MATLABdigital filterII目 錄摘 要IAbstractII1 數(shù)字濾波器11.1 數(shù)字濾波器介紹11.2 IIR數(shù)字濾波器設(shè)計(jì)原理21.2.2雙線性變換法32數(shù)字濾波器設(shè)計(jì)實(shí)現(xiàn)62.1 典型模擬濾波器比較62.2 設(shè)計(jì)步驟62.3 程序流程圖72.4音頻信號(hào)部分程序82.5仿真結(jié)果92.5.1濾波器性能仿真92.5.2濾波器除噪性能103總結(jié)11致謝11參考文獻(xiàn)12附 錄131 數(shù)字濾波器1.1 數(shù)字濾波器介紹數(shù)字濾波器是具有一定傳輸選擇特性的數(shù)字信號(hào)處理裝置,其輸入、輸出均為數(shù)字信號(hào),實(shí)質(zhì)上是一個(gè)由有限精度算法實(shí)現(xiàn)

10、的線性時(shí)不變離散系統(tǒng)。它的基本工作原理是利用離散系統(tǒng)特性對(duì)系統(tǒng)輸入信號(hào)進(jìn)行加工和變換,改變輸入序列的頻譜或信號(hào)波形,讓有用頻率的信號(hào)分量通過(guò),抑制無(wú)用的信號(hào)分量輸出。數(shù)字濾波器和模擬濾波器有著相同的濾波概念,根據(jù)其頻率響應(yīng)特性可分為低通、高通、帶通、帶阻等類(lèi)型,與模擬濾波器相比,數(shù)字濾波器除了具有數(shù)字信號(hào)處理的固有優(yōu)點(diǎn)外,還有濾波精度高(與系統(tǒng)字長(zhǎng)有關(guān))、穩(wěn)定性好(僅運(yùn)行在0與l兩個(gè)電平狀態(tài))、靈活性強(qiáng)等優(yōu)點(diǎn)。時(shí)域離散系統(tǒng)的頻域特性:，其中，分別是數(shù)字濾波器的輸出序列和輸入序列的頻域特性（或稱(chēng)為頻譜特性），是數(shù)字濾波器的單位取樣響應(yīng)的頻譜，又稱(chēng)為數(shù)字濾波器的頻域響應(yīng)。輸入序列的頻譜經(jīng)過(guò)濾波后。

11、因此，只要按照輸入信號(hào)頻譜的特點(diǎn)和處理信號(hào)的目的， 適當(dāng)選擇，使得濾波后的滿足設(shè)計(jì)的要求，這就是數(shù)字濾波器的濾波原理。數(shù)字濾波器根據(jù)其沖激響應(yīng)函數(shù)的時(shí)域特性，可分為兩種，即無(wú)限長(zhǎng)沖激響應(yīng)(IIR)數(shù)字濾波器和有限長(zhǎng)沖激響應(yīng)(FIR)數(shù)字濾波器。IIR 數(shù)字濾波器的特征是，具有無(wú)限持續(xù)時(shí)間沖激響應(yīng)，需要用遞歸模型來(lái)實(shí)現(xiàn)，其差分方程為： 系統(tǒng)函數(shù)為： 設(shè)計(jì)IIR濾波器的任務(wù)就是尋求一個(gè)物理上可實(shí)現(xiàn)的系統(tǒng)函數(shù)H(z)，使其頻率響應(yīng)H(z)滿足所希望得到的頻域指標(biāo)，即符合給定的通帶截止頻率、阻帶截止頻率、通帶衰減系數(shù)和阻帶衰減系數(shù)。1.2 IIR數(shù)字濾波器設(shè)計(jì)原理    IIR數(shù)字

12、濾波器是一種離散時(shí)間系統(tǒng)，其系統(tǒng)函數(shù)為 假設(shè)MN，當(dāng)MN時(shí),系統(tǒng)函數(shù)可以看作一個(gè)IIR的子系統(tǒng)和一個(gè)(M-N)的FIR子系統(tǒng)的級(jí)聯(lián)。IIR數(shù)字濾波器的設(shè)計(jì)實(shí)際上是求解濾波器的系數(shù)和，它是數(shù)學(xué)上的一種逼近問(wèn)題，即在規(guī)定意義上（通常采用最小均方誤差準(zhǔn)則）去逼近系統(tǒng)的特性。如果在S平面上去逼近，就得到模擬濾波器；如果在z平面上去逼近，就得到數(shù)字濾波器。設(shè)計(jì)通數(shù)字濾波器通?？梢詺w納為如圖所示的兩種常用方法。雙線性變換交換頻率數(shù)字高通模擬高通模擬低通方法1數(shù)字低通數(shù)字高通頻率交換雙線性變換模擬低通方法2 圖1-1 數(shù)字帶通濾波器設(shè)計(jì)的兩種方法 方法1:首先設(shè)計(jì)一個(gè)模擬原型低通濾波器，然后通過(guò)頻率變換成所

13、需要的模擬高通濾波器，最后再使用沖激不變法或雙線性變換成相應(yīng)的數(shù)字高通濾波器。方法2:先設(shè)計(jì)一個(gè)模擬原型低通濾波器，然后采用沖激響應(yīng)不變法或雙線性變換法將它轉(zhuǎn)換成數(shù)字原型低通濾波器，最后通過(guò)頻率變換把數(shù)字原型低通濾波器變換成所需要的數(shù)字高通濾波器。本課程設(shè)計(jì)采用第一種設(shè)計(jì)方法，先構(gòu)造一個(gè)巴特沃斯模擬低通濾波器，然后將模擬低通濾波器轉(zhuǎn)換成模擬高通濾波器，最后利用雙線性變換將模擬高通濾波器轉(zhuǎn)換成數(shù)字高通濾波器。1.2.2雙線性變換法為了克服沖激響應(yīng)法可能產(chǎn)生的頻率響應(yīng)的混疊失真，這是因?yàn)閺腟平面到平面是多值的映射關(guān)系所造成的。為了克服這一缺點(diǎn)，可以采用非線性頻率壓縮方法，將整個(gè)頻率軸上的頻率范圍壓

14、縮到-/T/T之間，再用z=esT轉(zhuǎn)換到Z平面上。也就是說(shuō)，第一步先將整個(gè)S平面壓縮映射到S1平面的-/T/T一條橫帶里；第二步再通過(guò)標(biāo)準(zhǔn)變換關(guān)系z(mì)=es1T將此橫帶變換到整個(gè)Z平面上去。這樣就使S平面與Z平面建立了一一對(duì)應(yīng)的單值關(guān)系，消除了多值變換性，也就消除了頻譜混疊現(xiàn)象，映射關(guān)系如圖1-3圖1-3雙線性變換的映射關(guān)系為了將S平面的整個(gè)虛軸j壓縮到S1平面j1軸上的-/T到/T段上，可以通過(guò)以下的正切變換實(shí)現(xiàn)式中,T仍是采樣間隔。當(dāng)1由-/T經(jīng)過(guò)0變化到/T時(shí)，由-經(jīng)過(guò)0變化到+，也即映射了整個(gè)j軸。將式（1-9）寫(xiě)成將此關(guān)系解析延拓到整個(gè)S平面和S1平面，令j=s，j1=s1，則得再將S

15、1平面通過(guò)以下標(biāo)準(zhǔn)變換關(guān)系映射到Z平面z=es1T從而得到S平面和Z平面的單值映射關(guān)系為：這兩個(gè)關(guān)系式是S平面與Z平面之間的單值映射關(guān)系，這種變換都是兩個(gè)線性函數(shù)之比，因此稱(chēng)為雙線性變換首先,把可得即S平面的虛軸映射到Z平面的單位圓。其次，將代入，得因此由此看出，當(dāng)<0時(shí)，|z|<1；當(dāng)>0時(shí)，|z|>1。也就是說(shuō)，S平面的左半平面映射到Z平面的單位圓內(nèi)，S平面的右半平面映射到Z平面的單位圓外，S平面的虛軸映射到Z平面的單位圓上。因此，穩(wěn)定的模擬濾波器經(jīng)雙線性變換后所得的數(shù)字濾波器也一定是穩(wěn)定的。雙線性變換法優(yōu)缺點(diǎn)：雙線性變換法與脈沖響應(yīng)不變法相比，其主要的優(yōu)點(diǎn)是避免了

16、頻率響應(yīng)的混疊現(xiàn)象。這是因?yàn)镾平面與Z平面是單值的一一對(duì)應(yīng)關(guān)系。S平面整個(gè)j軸單值地對(duì)應(yīng)于Z平面單位圓一周，即頻率軸是單值變換關(guān)系。這個(gè)關(guān)系重寫(xiě)如下：上式表明，S平面上與Z平面的成非線性的正切關(guān)系，如圖1-4所示。由圖1-4看出，在零頻率附近，模擬角頻率與數(shù)字頻率之間的變換關(guān)系接近于線性關(guān)系；但當(dāng)進(jìn)一步增加時(shí)，增長(zhǎng)得越來(lái)越慢，最后當(dāng)時(shí)，終止在折疊頻率=處，因而雙線性變換就不會(huì)出現(xiàn)由于高頻部分超過(guò)折疊頻率而混淆到低頻部分去的現(xiàn)象，從而消除了頻率混疊現(xiàn)象。圖1-4雙線性變換法的頻率變換關(guān)系但是雙線性變換的這個(gè)特點(diǎn)是靠頻率的嚴(yán)重非線性關(guān)系而得到的，圖1-4所示。由于這種頻率之間的非線性變換關(guān)系，就產(chǎn)

17、生了新的問(wèn)題。首先，一個(gè)線性相位的模擬濾波器經(jīng)雙線性變換后得到非線性相位的數(shù)字濾波器，不再保持原有的線性相位了；其次，這種非線性關(guān)系要求模擬濾波器的幅頻響應(yīng)必須是分段常數(shù)型的，即某一頻率段的幅頻響應(yīng)近似等于某一常數(shù)（這正是一般典型的低通、高通、帶通、帶阻型濾波器的響應(yīng)特性），不然變換所產(chǎn)生的數(shù)字濾波器幅頻響應(yīng)相對(duì)于原模擬濾波器的幅頻響應(yīng)會(huì)有畸變，如圖1-5所示。圖1-5雙線性變換法幅度和相位特性的非線性映射對(duì)于分段常數(shù)的濾波器，雙線性變換后，仍得到幅頻特性為分段常數(shù)的濾波器，但是各個(gè)分段邊緣的臨界頻率點(diǎn)產(chǎn)生了畸變，這種頻率的畸變，可以通過(guò)頻率的預(yù)畸來(lái)加以校正。也就是將臨界模擬頻率事先加以畸變，

18、然后經(jīng)變換后正好映射到所需要的數(shù)字頻率上。2數(shù)字濾波器設(shè)計(jì)實(shí)現(xiàn)2.1 典型模擬濾波器比較1、Butterworth巴特沃斯濾波器：它具有單調(diào)下降的幅頻特性；即最平幅度。2、Chebyshev切比雪夫?yàn)V波器：在通帶或阻帶等波紋，可提高選擇性。3、Bessel貝塞爾濾波器：在通帶內(nèi)有較好的線性相位特性。4、Ellipse橢圓濾波器：其選擇性相對(duì)前三種是最好的。此處選擇巴特沃斯主要是想獲得最平穩(wěn)的幅頻響應(yīng)。而不計(jì)較相位特性，而用雙線性變換法也會(huì)將線性相位變?yōu)榉蔷€性相位。2.2 設(shè)計(jì)步驟根據(jù)以上IIR數(shù)字濾波器設(shè)計(jì)方法，下面運(yùn)用雙線性變換法基于MATLAB設(shè)計(jì)一個(gè)IIR高通濾波器。 (1)確定性能指標(biāo)

19、在設(shè)計(jì)高通濾波器之前,首先根據(jù)工程實(shí)際的需要確定濾波器的技術(shù)指標(biāo): 通帶截止頻率fp=3200Hz；阻帶截止頻率fs=3000Hz；阻帶最小衰減As=100dB和通帶最大衰減Ap=1dB(2)把頻率轉(zhuǎn)化為數(shù)字角頻率 wp=2*pi*fp*T; ws=2*pi*fs*T; (3)頻率預(yù)畸變以及模擬高通性能指標(biāo)轉(zhuǎn)換成模擬低通性能指標(biāo)用=2/T*tan(w/2)對(duì)高通數(shù)字濾波器H(z)的數(shù)字邊界頻率預(yù)畸變,得到高通模擬濾波器H(s)的邊界頻率主要是通帶截止頻率Wp;阻帶截止頻率Ws的轉(zhuǎn)換。抽樣頻率 Fs=44100Hz。通帶截止頻率 Wp2=(2/T)*tan(wp/2)阻帶截止頻率 Ws2=(2/

20、T)*tan(ws/2)(4) 歸一化處理 ； ； ； ； ； 由式2-3，2-4，2-5計(jì)算出N，查表可得模擬低通濾波器的階數(shù)，從而由下式確定模擬高通濾波器的參數(shù) (5)模擬低通濾波器的構(gòu)造借助巴特沃斯濾波器模型得到模擬低通濾波器的傳輸函數(shù)Ha(s)。N,Wn=buttord(wp2,ws2,Ap,As,'s');z,p,k=buttap(N);Bap,Aap=zp2tf(z,p,k); G1=tf(Bap,Aap);(6)模擬低通濾波器轉(zhuǎn)換成模擬高通濾波器調(diào)用lp2bp函數(shù)將模擬低通濾波器轉(zhuǎn)化為模擬高通濾波器。Bbs,Abs=lp2hp(Bap,Aap,Wn); (7)模擬

21、高通濾波器轉(zhuǎn)換成數(shù)字高通濾波器利用雙線性變換法將模擬高通濾波器Ha(s)轉(zhuǎn)換成數(shù)字高通濾波器H(z)。Bbz,Abz=bilinear(Bbs,Abs,FS); 2.3 程序流程圖首先確定性能指標(biāo)，把頻率轉(zhuǎn)化為數(shù)字角頻率，進(jìn)而在進(jìn)行頻率預(yù)畸變，用=2/T*tan(w/2)對(duì)高通數(shù)字濾波器H(z)的數(shù)字邊界頻率預(yù)畸變,得到高通模擬濾波器H(s)的邊界頻率主要是通帶截止頻率Wp;阻帶截止頻率Ws的轉(zhuǎn)換。抽樣頻率Fs=44100Hz。上述準(zhǔn)備工作做好之后，就先把模擬高通性能指標(biāo)轉(zhuǎn)換成模擬低通性能指標(biāo)，然后設(shè)計(jì)模擬低通濾波器，借助巴特沃斯濾波器模型得到模擬低通濾波器的傳輸函數(shù)Ha(s)。然后將模擬低通

22、濾波器轉(zhuǎn)化為模擬高通濾波器。最后利用雙線性變換法將模擬高通濾波器Ha(s)轉(zhuǎn)換成數(shù)字高通濾波器H(z)。開(kāi)始讀入數(shù)字濾波器技術(shù)指標(biāo) 將指標(biāo)轉(zhuǎn)換成歸一化模擬低通濾波器的指標(biāo) 設(shè)計(jì)歸一化的模擬低通濾波器階數(shù)N和3db截止頻率 模擬低通濾波器變換成模擬高通濾波器 模擬高通濾波器轉(zhuǎn)換成數(shù)字高通濾波器 結(jié)束 圖2-1程序流程圖2.4音頻信號(hào)部分程序（1）寫(xiě)入聲音信號(hào) fnose =10;T=length(x)/fs;t=T/length(x):T/length(x):T; j=1:1:length(x);（2）產(chǎn)生噪聲并加到聲音中 nose =sin(2*pi*fnose*T*j);y=x+5*nose

23、'b,a=butter(N,Wn,'high');Y=fft(y);2.5仿真結(jié)果2.5.1濾波器性能仿真源程序設(shè)計(jì)了模擬低通濾波器、模擬高通濾波器與數(shù)字高通濾波器，對(duì)數(shù)字高通濾波器的性能仿真如下：2.5.2濾波器除噪性能為了實(shí)現(xiàn)濾波器的應(yīng)用，程序中加入了有噪聲的音頻信號(hào)，通過(guò)對(duì)其濾波處理，來(lái)顯示數(shù)字高通濾波器的功能，下面顯示未加入噪聲，加入噪聲和濾波后的頻譜波形圖。 濾波前加噪信號(hào)的波形和頻譜濾波后信號(hào)的波形和頻譜3總結(jié)這次課設(shè)下來(lái)，對(duì)設(shè)計(jì)高通數(shù)字濾波器的整個(gè)過(guò)程有了很好的掌握，懂得了設(shè)計(jì)濾波器的基本方法，對(duì)雙線性變換法，巴特沃斯濾波器有了一定了解，同時(shí)呢也熟悉了MA

24、TLAB的環(huán)境，鞏固了相關(guān)知識(shí)。初步掌握了MATLAB語(yǔ)言在數(shù)字信號(hào)處理中一些基本庫(kù)函數(shù)的調(diào)用和編寫(xiě)基本程序等應(yīng)用；熟悉了濾波器設(shè)計(jì)的一般原理，對(duì)濾波器有了一個(gè)感性的認(rèn)識(shí)；學(xué)會(huì)了數(shù)字高通濾波器設(shè)計(jì)的一般步驟；加深了對(duì)濾波器設(shè)計(jì)中產(chǎn)生誤差的原因以及雙線性變換法優(yōu)缺點(diǎn)的理解和認(rèn)識(shí)。總之，使理論聯(lián)系了實(shí)際，鞏固并深化了對(duì)課本基本知識(shí)的認(rèn)識(shí)和理解，使理論得以升華。在做本次課程設(shè)計(jì)的過(guò)程中，我深深地感受到了自己所學(xué)到知識(shí)的有限，明白了只學(xué)好課本上的知識(shí)是不夠的，要通過(guò)圖書(shū)館和互聯(lián)網(wǎng)等各種渠道來(lái)擴(kuò)充自己的知識(shí)。在實(shí)驗(yàn)過(guò)程中我們?cè)?jīng)遇到過(guò)問(wèn)題。但是從中我們學(xué)習(xí)到了如何對(duì)待遇到的困難，進(jìn)一步培養(yǎng)了我們一絲不茍

25、的科學(xué)態(tài)度和不厭其煩的耐心。所有的這些心得會(huì)對(duì)我以后的學(xué)習(xí)和工作有幫助作用，忠心感謝學(xué)校給我們提供這次實(shí)驗(yàn)機(jī)會(huì)。致謝本設(shè)計(jì)過(guò)程中得到老師的悉心指導(dǎo)，設(shè)計(jì)過(guò)程中老師多次幫我分析思路，開(kāi)拓視角，在我遇到困難想要放棄的時(shí)候給予我最大的支持和鼓勵(lì)。老師那嚴(yán)謹(jǐn)求實(shí)的治學(xué)態(tài)度、踏實(shí)堅(jiān)韌的工作精神，將使我終身受益；另外同學(xué)和朋友的熱心幫助也使我學(xué)到了不少的知識(shí)。再多華麗的語(yǔ)言也顯得蒼白，在此，謹(jǐn)向老師和同學(xué)們致以誠(chéng)摯的謝意和崇高的敬意。因?yàn)樗麄兾覍?duì)這次課程設(shè)計(jì)有了深刻的印象，雖然只有短暫的一周，但在這期間，卻讓我受益匪淺，我深刻地認(rèn)識(shí)到了“理論聯(lián)系實(shí)際”這句話的重要性與真實(shí)性。參考文獻(xiàn)1 程佩青，數(shù)字信號(hào)處

26、理.北京：清華大學(xué)出版社,2007.22 景振毅，張澤兵，董霖.北京：中國(guó)鐵道出版社，2008.123 從玉良.數(shù)字信號(hào)處理原理及其MATLAB實(shí)現(xiàn)M.北京：電子工業(yè)出版社.2009.74 胡廣書(shū).數(shù)字信號(hào)處理理論、算法與實(shí)現(xiàn)M.北京：清華大學(xué)出版社.2003,85 劉波MATLAB信號(hào)處理北京：電子工業(yè)出版社，2006 6 施陽(yáng)等MATLAB語(yǔ)言工具箱西安：西北工業(yè)大學(xué)出版社，1999 7 丁玉美數(shù)字信號(hào)處理西安電子科技大學(xué)出版社 第二版8 倪養(yǎng)華數(shù)字信號(hào)處理與實(shí)現(xiàn)上海：上海交通大學(xué)出版社，1998 9 張葛祥，李娜MATLAB仿真技術(shù)與應(yīng)用北京：清華大學(xué)出版社，2003 10 飛思科技產(chǎn)品研發(fā)中心MATLAB7輔助信號(hào)處理技術(shù)與應(yīng)用北京：電子工業(yè)出版社， 2005 附 錄源程序FS=8000;T=1/FS;