利用MATLAB仿真軟件系統(tǒng)結(jié)合雙線性變換法設計一個數(shù)字巴特沃斯高通IIR濾波器要點

1、 摘要Matlab是一個矩陣設計平臺，傳統(tǒng)數(shù)字濾波器設計需要大量的計算，但是利用Matlab可以快速實現(xiàn)濾波器的設計與仿真，而且頻譜分析功能強大，在數(shù)字信號處理中發(fā)揮了巨大的作用。本次實驗中，用雙線性不變法設計高通巴特沃斯IIR數(shù)字濾波器，介紹了設計步驟，然后在Matlab環(huán)境下進行了仿真與調(diào)試，實現(xiàn)了設計目標。關鍵詞：Matlab 數(shù)字濾波器 雙線性變換法 IIR摘要Matlab是一個矩陣設計平臺，傳統(tǒng)數(shù)字濾波器設計需要大量的計算，但是利用Matlab可以快速實現(xiàn)濾波器的設計與仿真，而且頻譜分析功能強大，在數(shù)字信號處理中發(fā)揮了巨大的作用。本次實驗中，用雙線性不變法設計高通巴特沃斯IIR數(shù)字濾

2、波器，介紹了設計步驟，然后在Matlab環(huán)境下進行了仿真與調(diào)試，實現(xiàn)了設計目標。關鍵詞：Matlab 數(shù)字濾波器 雙線性變換法 IIR摘要Matlab是一個矩陣設計平臺，傳統(tǒng)數(shù)字濾波器設計需要大量的計算，但是利用Matlab可以快速實現(xiàn)濾波器的設計與仿真，而且頻譜分析功能強大，在數(shù)字信號處理中發(fā)揮了巨大的作用。本次實驗中，用雙線性不變法設計高通巴特沃斯IIR數(shù)字濾波器，介紹了設計步驟，然后在Matlab環(huán)境下進行了仿真。關鍵詞：Matlab 數(shù)字濾波器 雙線性變換法 1設計要求和說明利用MATLAB仿真軟件系統(tǒng)結(jié)合雙線性變換法設計一個數(shù)字巴特沃斯高通IIR濾波器。MATLAB工具箱為濾波器的設

3、計應用提供了豐富而簡便的方法，使原來的非常繁瑣復雜的程序設計變成簡單的程序調(diào)用。1.1 設計原理濾波器，顧名思義，就是對系統(tǒng)輸入信號進行濾波。那個數(shù)字濾波器的數(shù)學運算通常用兩種方法來表示。一種是頻域法，即利用FFT快速運算辦法對輸入信號進行離散傅里葉變換，分析其頻譜，然后根據(jù)所希望的頻率特性進行濾波，再利用傅里葉反變換來輸出出時域信號。N階低通巴特沃斯濾波器的特性為： 其中，c為通帶寬度，即截止頻率。當階次N增大時，濾波器的特性曲線變的更加陡峭，其特性就越接近理想的矩形幅頻特性。巴特沃斯濾波器屬于全極點設置，他的極點由以下公式來確定。式中，。所以，在s平面上有2N個極點。這些極點是等間隔的分布

4、在半徑為c的圓周上，并且這些極點都是成復共軛對出現(xiàn)的，極點位置與虛軸是對稱的，且不在虛軸上。雙線性變換法是將平面壓縮變換到某一中介平面的一條橫帶里，再通過標準變換關系將此帶變換到整個z平面上去，這樣就使平面與平面之間建立一一對應的單值關系，消除了多值變換性。為了將平面的軸壓縮到平面的軸上的到一段上，可以通過以下的正切變換來實現(xiàn)：這樣當由經(jīng)變化到時，由經(jīng)過0變化到，也映射到了整個 軸。將這個關系延拓到整個平面和平面，則可以得到再將平面通過標準變換關系映射到平面，即令得到同樣對z求解，得到雙線性變換平面的映射關系如圖一：圖一 雙線性變換平面映射圖這樣的變換叫做雙線性變換。1.2 設計思路數(shù)字濾波器

5、一般設計過程為： 按照實際需要，確定濾波器性能要求； 用一個因果、穩(wěn)定的離散線性時不變系統(tǒng)去逼近這一個性能指標； 用有限精度的計算去實現(xiàn)所設計的系統(tǒng)； 通過模擬驗證所設計的系統(tǒng)是否符合給定的性能要求。 IIR濾波器設計的主要方法是先設計低通模擬濾波器，然后轉(zhuǎn)換為高通、帶通或帶阻數(shù)字濾波器。對于其他如高通，帶通，則通過頻率變換轉(zhuǎn)換為設計相應的高通，帶通等。在設計的全過程的各個步驟，matlab都提供相應的工具箱函數(shù)，使得IIR數(shù)字濾波器設計變得非常簡單?？偟膩碚f，我的設計思路主要有以下兩種：思路一：從歸一化模擬低通原型出發(fā)，先在模擬域內(nèi)經(jīng)頻率變換成為所需類型的模擬濾波器；然后進行雙線性變換，由S

6、域變換到Z域，而得到所需類型的數(shù)字濾波器。歸一化模擬低通原型數(shù)字高,帶通或帶阻模擬高,帶通或帶阻 模擬域 沖激響應不變法頻率變換 雙線性變換法圖二 先頻率變換再離散思路二：先進行雙線性變換，將模擬低通原型濾波器變換成數(shù)字低通濾波器；然后在Z域內(nèi)經(jīng)數(shù)字頻率變換為所需類型的數(shù)字濾波器。 數(shù)字高,帶通或帶阻數(shù)字原型低通歸一化模擬低通原型 數(shù)字域 雙線性變換法 頻率變換圖三 先離散再頻率變換以上兩種思路都可以，我最后選擇了第一種思路進行設計，即先在模擬域內(nèi)經(jīng)頻率變換成為所需類型的模擬濾波器；然后進行雙線性變換，由S域變換到Z域，而得到所需類型的數(shù)字濾波器。1.3 設計過程下面結(jié)合一個實例來說明巴特沃斯

7、濾波器的設計過程。 例：設計一個巴特沃斯濾波器函數(shù)。該濾波器的性能指標為：通帶的截止頻率為c=10000rad/s， 通帶的最大衰減Ap=3 dB ，阻帶的截止頻率s=40000rad/s ，阻帶的最小衰減 As=35 dB 。 解： （1）求相關參數(shù)（2） 確定參數(shù)N取N=3 ， 根據(jù) 得到所以 ，極點形式可以表示為 即滿足系統(tǒng)性能指標的函數(shù)2 Matlab軟件的程序設計程序設計中用到主要函數(shù)：z,p,k=buttap(N); %創(chuàng)建Buttord低通濾波器原型Bap,Aap=zp2tf(z,p,k); %由零極點轉(zhuǎn)換為傳遞函數(shù)的形式Bbs,Abs=lp2hp(Bap,Aap,Wn); %模

8、擬低通變高通Bbz,Abz=bilinear(Bbs,Abs,FS); %雙線性變換2.1 數(shù)字域指標變換成模擬域指標 其程序為：fp = 400 ; Fs = 300; Rp = 1; Rs = 20; wp =fp*2*pi; ws =fs*2*pi;FS=1000;T=1/FS;程序執(zhí)行結(jié)果為：2.2 數(shù)字域頻率進行預畸變令 Wp=wp/(FS); Ws=ws/(FS); 其畸變程序為：wp2=2*tan(Wp/2)/T; ws2=2*tan(Ws/2)/T; 經(jīng)過預畸變，可以發(fā)現(xiàn)頻率變?yōu)?3 模擬濾波器的設計其程序為N,Wn = buttord(wp2,ws2,Rp,Rs,s) %設計

9、模擬濾波器z,p,k=buttap(N); %創(chuàng)建Buttord低通濾波器原型Bap,Aap=zp2tf(z,p,k); %由零極點轉(zhuǎn)換為傳遞函數(shù)的形式figure(1)freqs(Bap,Aap);%模擬低通濾波器的頻率響應title(模擬濾波器（低通原型）的頻率響應)Bbs,Abs=lp2hp(Bap,Aap,Wn); %模擬低通變高通figure(2)freqs(Bbs,Abs);title(模擬濾波器(高通原型)的頻率響應)程序執(zhí)行后可以發(fā)現(xiàn)其頻率響應為：模擬濾波器的低通原型波形如圖四所示：圖四 模擬濾波器（低通原型）的頻率響應 模擬濾波器高通原型的波形如圖五所示：圖五 模擬濾波器（高

10、通原型）的頻率響應2.4 模擬濾波器變成數(shù)字濾波器其程序為：Bbz,Abz=bilinear(Bbs,Abs,FS); %用雙線性變換法設計數(shù)字濾波器 freqz(Bbz,Abz,512,FS); 程序運行的結(jié)果如圖六：圖六 數(shù)字濾波器的頻率響應由于使用的是雙線性不變法設計的，其相位為非線性。此處主要是基于要獲得嚴格的頻率響應，以及較準確地控制截止頻率的位置，故畫出了詳細的幅頻響應。得到數(shù)字濾波器的詳細頻率響應圖形圖七 詳細的幅頻響應（通帶衰減為5dB）當通帶變?yōu)?dB（之前為1dB）時，阻帶不變時，其波形如下圖十。對比圖七可知，其在通帶處的衰減變?yōu)榱?dB，曲線平滑了一些。從調(diào)試的過程中發(fā)現(xiàn)

11、：通帶衰減越小，濾波器的性能越好。 阻帶衰減越大，濾波器的性能越好。3.程序調(diào)試中出現(xiàn)的問題 freqs(Bbs,Abs)與freqz(Bbz,Abz,512,FS)剛開始時分得不是很清楚，后來發(fā)現(xiàn)前者是求模擬域頻率響應的，而后者是求數(shù)字域頻率響應的。后者的用法如下：h,w=freqz(b,a,n,whole);或h,f=freqz(b,a,n,whole,Fs);式中，b,a為數(shù)字濾波器分子和分母多項式的系數(shù)，n為復數(shù)頻率的響應點數(shù)，為整數(shù)，最好為2的冪，缺省時為512；Fs為采樣頻率，單位Hz。如果給定該值，則f位置輸出為頻率Hz，若沒有給定，則按角頻率（Angular frequency

12、）給定f的頻率矢量；whole表示返回的頻率f或w值包含z平面整個單位圓頻率矢量，即02；缺省時，頻率f或w值包含z平面上半單位圓（0）之間等間距n個點頻率矢量。h為復頻率響應；w為n點頻率向量（單位rad）；f為n點頻率向量(Hz)。函數(shù)返回值缺省時，繪制幅頻響應和相頻響應圖4小結(jié) 總的來說，這次課程設計讓我對 MATLAB 有了更深刻的了解，對數(shù)字濾波器的設計流程有了大致的了解，掌握了一些設計濾波器的基本方法，提高了理論用于實踐的能力，掌握了更多專業(yè)相關的使用知識與技能。同時，也暴露了我很多的不足，在以后的學習中， 將進一步發(fā)現(xiàn)并克服缺點。參考文獻1 陳懷琛 MATLAB在電子信息課程中的

13、應用（第二版） 電子工業(yè)出版社 20062 郭仕劍 MATLAB7.X數(shù)字信號處理 人民郵電出版社 20063 陳貴明 用MATLAB語言處理數(shù)字信號與數(shù)字圖像 科學出版社 20004 王家文 MATLAB7.0圖形圖像處理 國防工業(yè)出版社 2006附錄：實驗總程序：fp = 400 ; fs = 300; Rp = 1; %通帶最大衰減Rp=1dBRs = 20; %阻帶最小衰減Rs=20dBwp =fp*2*pi; ws =fs*2*pi;FS=1000;T=1/FSWp=wp/(FS);Ws=ws/(FS);wp2=2*tan(Wp/2)/T; %預畸變求濾波器通帶臨界頻率ws2=2*t

14、an(Ws/2)/T; %預畸變求濾波器阻帶臨界頻率 N,Wn = buttord(wp2,ws2,Rp,Rs,s) z,p,k=buttap(N); %創(chuàng)建Buttord低通濾波器原型Bap,Aap=zp2tf(z,p,k); %由零極點轉(zhuǎn)換為傳遞函數(shù)的形式figure(1)freqs(Bap,Aap); %模擬低通濾波器的頻率響應title(模擬濾波器（低通原型）的頻率響應)Bbs,Abs=lp2hp(Bap,Aap,Wn); %模擬低通變高通figure(2)freqs(Bbs,Abs);title(模擬濾波器的頻率響應) Bbz,Abz=bilinear(Bbs,Abs,FS); %雙

15、線性變換figure(3)freqz(Bbz,Abz,512,FS);title(數(shù)字濾波器的頻率響應)hw,w=freqz(Bbz,Abz,512); figure(4)plot(w/pi,20*log10(abs(hw);gridaxis(0,1,-200,10)title( 數(shù)字濾波器詳細的頻率響應)xlabel(w/pi);ylabel(幅度(dB);wp=0.8*pi;ws=0.6*pi;OmegaP=2*1000*tan(wp/2);OmegaS=2*1000*tan(ws/2);lamdas=OmegaP/OmegaS;N=0.5*log10(10.(20/10)-1)/(10.(1/10)-1)/log10(lamdas);Wn= 4.8890e+003az=0 0 0 0 1;bz=1 2.613 3.414,2.613,1;Bbs,Abs=lp2hp(az,bz,Wn) Bbz,Abz=bilinear(Bbs,Abs,1000); hw,w=freqz(Bbz,Abz,512); figure(5)plot(w/pi,20*log10(abs(hw);gridaxis(0,1,-200,10)set(gca,XTickMode,manual,XTick,0,Ws/(pi),Wp/