1、課程設計課程設計名稱：用雙線性變換法設計原型低通為橢圓型的數字IIR帶通濾波器專業 班級：xxxxxx學生姓名：xxx學 號： xxxxxxxxxxxx指導教師 ：xxx課程設計時間：2014.6.16 至2014.6.20電子信息工程專業課程設計任務書學生姓名XXX專業班級XXXX學號XXXXXX題目用雙線性變換法設計原型低通為橢圓型的數字IIR低通濾波器課題性質其他課題來源自擬指導教師XX同組姓名主要內容用雙線性變換法設計原型低通為橢圓型的數字IIR低通濾波器，要求通帶邊界頻率為400Hz,阻帶邊界頻率分別為 500Hz,通帶最大哀減1dB,阻帶最 小衰減40dB,抽樣頻率為2000Hz,

2、用MATLAB畫出幅頻特性，畫出并分析濾 波器傳輸函數的零極點；信號 x(t) =x1(t)+x2(t) =sin(2研1t)+sin(2nf2t)經過該濾波器,其中f1 =300Hz, f2 =600Hz,濾波器的輸出y(t)是什么？用 Matlab驗證你的結論 并給出Xi8,x，y的圖形。任務要求1、掌握用雙線性變換法設計原型低通為橢圓型的數字IIR低通濾波器的原理和設計方法。2、求出所設計濾波器的 Z變換。3、用MATLAB畫出幅頻特性圖。4、驗證所設計的濾波器。奔f文獻1、程佩青著，數字信號處理教程，清華大學出版社，20012、Sanjit K. Mitra著，孫洪，余翔宇譯，數字信號

3、處理實驗指導書 (MATLAB 版)，電子工業出版社，2005年1月3、郭仕劍等，MATLAB 7.x數字信號處理，人民郵電出版社，2006年4、胡廣書，數字信號處理理論算法與實現，清華大學出版社，2003年審查意見指導教師簽字：教研室主任簽字：2014年06月12日說明：本表由指導教師填寫，由教研室主任審核后下達給選題學生，裝訂在設計(論文)首頁一需求分析和技術要求1、需求分析濾波器從廣義上來說對特定的頻點或頻點以外的頻率進行有效濾波的電路， 這種電路保留輸入信號中的有用信息，濾除不需要的信息，從而達到信號的檢測、 提取、識別等不同的目的。如果處理的信號是時域離散信號，那么相應的處理系統就稱

4、為數字濾波器，由于在實際工作中被處理的信號都是幅度量化的數字信 號，因此，數字濾波器實際上是用有限精度的算法實現一個線性時不變的時域離 散系統。數字濾波器的種類很多，分類方法也不同，可以從功能上分類，也可以從實 現方法上及設計方法上來分類等等。濾波器在功能上總的可分為四類，即低通 (LP)、高通(HP)、帶通(BP)、帶阻(BS)濾波器等，從實現方法上，由有限長沖激 響應的數字濾波器被稱為FIR濾波器，具有無限長沖激響應的數字濾波器稱作 IIR濾波器。橢圓濾波器(日liptic filter )又稱考爾濾波器(Causer filter ),是在 通帶和阻帶等波紋的一種濾波器。 橢圓濾波器相比

5、其他類型的濾波器，在階數相 同的條件下有著最小的通帶和阻帶波動。 它在通帶和阻帶的波動相同，這一點區 別于在通帶和阻帶都平坦的巴特沃斯濾波器，以及阻帶平坦、通帶等波紋的切比 雪夫濾波器。利用雙線性變換法將模擬傳輸信號 Ha(s)變換為數字傳輸函數G(z),從而是 z域的數字傳輸函數保留s域的模擬傳輸函數的基本性質。 設計成的IIR數字低 通濾波器能夠去掉信號中不必要的高頻成分， 降低采樣頻率，避免頻率混淆，去 掉高頻干擾。2、技術要求1、掌握用雙線性變換法設計原型低通為橢圓型的數字IIR低通濾波器的原理和設計方法。2、求出所設計濾波器的Z變換。3、用MATLAB畫出幅頻特性圖。4、驗證所設計的

6、濾波器。設計原理及設計思路1、設計原理數字濾波器是一種用來過濾時間離散信號的數字系統，通過對抽樣數據進行 數學處理來達到頻域濾波的目的。數字濾波器也是具有一定傳輸選擇特性的數字 信號處理裝置，其輸入、輸出均為數字信號，實質上是一個由有限精度算法實現 的線性時不變離散系統。IIR數字濾波器采用遞歸型結構，即結構上帶有反饋環路。IIR濾波器運算結構通常由延時、乘以系數和相加等基本運算組成，可以組合成直接型、正準型、 級聯型、并聯型四種結構形式，都具有反饋回路。數字濾波器根據其沖激響應函 數的時域特性，可分為兩種，即無限長沖激響應(IIR)數字濾波器和有限長沖激響 應(FIR)數字濾波器。IIR數字

7、濾波器的特征是，具有無限持續時間沖激響應，需 要用遞歸模型來實現，其差分方程為：HNj(n) = 2 %” (一力十 t；y ( n j )(2.1)M之H =1+之口高-，e(2.2)設計IIR濾波器的任務就是尋求一個物理上可實現的系統函數H(z),使其頻率響應H(z)滿足所希望得到的頻域指標，即符合給定的通帶截止頻率、阻帶截止頻 率、通帶衰減系數和阻帶衰減系數。本次課程設計采用雙線性變換法。雙線性變換法設計IIR數字濾波器雙線性變換法主要是采用非線性頻率壓縮方法，將整個頻率軸上的頻率范圍 壓縮到-兀/T兀/T之間，再用z=esT轉換到Z平面上。也就是說，第一步先將 整個S平面壓縮映射到S平

8、面的-冗/ T兀/T 一條橫帶里；第二步再通過標準變 換關系z=es1T將此橫帶變換到整個Z平面上去。這樣就使S平面與Z平面建立了 一一對應的單值關系，消除了多值變換性，也就消除了頻譜混疊現象，映射關系 如圖1所示圖1雙線性變換的映射關系為了將S平面的整個虛軸j Q壓縮到S1平面j Q1軸上的-九/T到冗/T段上，可以通過以下的正切變換實現2c 1Tl（2 3）ta = tan i（2.3）一，一T 2 2 J式中，T仍是米樣間隔。當Q 1由-冗/T經過0變化到冗/T時，。由-8經過0變化到+oo,也即映射了整個j Q軸。將式（2.3）寫成2 ej 1T/2 _ej1T/2（2.4）T ejQ

9、T/2. ej干/2將此關系解析延拓到整個 S平面和S1平面，令j Q=s, j Q1=s1,則得S1T/2-S1T/22 e- eQ =，SS1T/2-S1T/2T e e2 Sj2= tanhl = 一T 2 Ts1T1 -eS1T1 e再將S1平面通過以下標準變換關系映射到 Z平面z=es1T從而得到S平面和Z平面的單值映射關系為:2 1 -zs 二jT 1 zdT21 S Sz=2 TI_dT21 S S2T(2.(5)(2.(6)(2.(7)(2.(8)式（2.7）與式（2.8）是S平面與Z平面之間的單值映射關系，這種變換都 是兩個線性函數之比，因此稱為雙線性變換式（2.6）與式（2

10、.7）的雙線性變換符合映射變換應滿足的兩點要求。首先，把z=ej可得s=T=j2n 但jT1+e0T 12 J(2.9)即S平面的虛軸映射到Z平面的單位圓。其次，將s=(7 +j Q代入式(2.8 ),得2.一.,- j1 z = T(2.10)因此工|z-T(2.11)由此看出，當(r0時，|z|1。也就是說，S平面的左半平面映射到Z平面的單位圓內，S平面的右半平面映射到Z平面的單位圓外, S平面的虛軸映射到Z平面的單位圓上。因此，穩定的模擬濾波器經雙線性變換后所得的數字濾波器也一定是穩定的。2、設計思路具體實現方法如下:(1)確定數字低通濾波器的技術指標：通帶截止頻率 cop、通帶衰減Rp

11、阻帶截止頻率cos、阻帶衰減Rs。(2)將數字低通濾波器的技術指標轉換成模擬低通濾波器的技術指標。如果采用脈沖響不變法，邊界頻率的轉換關系為(2.12)如果采用雙線性變換法，邊界頻率的轉換關系為(2.13)(3)按照模擬低通濾波器的技術指標設計模擬低通濾波器。(4)利用雙線性變換法將模擬濾波器 Ha(s),從s平面轉換到z平面，得到數字低通濾波器系統函數H(z) 0(5)數字低通技術指標為： p=0.4 兀 rad , Rp=1dB; s=0.5 兀 rad , Rs=40dB(6)模擬低通的技術指標為：通帶截止頻率 p=(2/T)*tan( p/2)=tan(0.392)=0.0068阻帶截

12、止頻率 s=(2/T)*tan( s/2)=tan(0.628)=0.011通帶最小衰減Rp=1dB和阻帶最大衰減 Rs=40dB;歸一化截止角頻率 p=2pi*Fs/Ft ; w s=2pi*Fs/Ft(7)利用模擬橢圓濾波器設計數字濾波器。通帶截止頻率為：p=0.4*pi ；阻帶截止頻率為：s=0.5*pi ;通帶最大衰減為：Rp=1;阻帶最大衰減為：Rs=40；設定周期為1s；模擬低通濾波器的生成：num1,den1=ellip(N,Rp,Rs,Wn,s) ；最后實現輸入輸出、幅頻特性、相頻特性的圖形。三程序流程圖開始確定數字低通濾波器的各項性能指標數字低通濾波器的技術指標轉換成模擬低通

13、濾波器的技術指標按照模擬低通濾波器的技術指標設計模擬低通濾波器利用雙線性變換法將模擬濾波器轉變為數字低通濾波器利用模擬橢圓濾波器設計數字濾波器結束圖2濾波器設計流程方框圖四程序源代碼MATLAB程序如下：clear allclose allCoeff_LF = 1;% 低頻信號系數F1(300 Hz)Coeff_HF = 1;% 高頻信號系數F2(600 Hz)Coeff_NOISE = 0;% 噪聲系數，如0(沒有噪聲),1, sqrt(2), 3等等.F1 = 300; % HzF2 = 600; % HzFS = 2000; % Hz抽樣頻率%時間定義EndTime = 1/FS * 1

14、023;Time = 0:1/FS:EndTime;N = length(Time);% 產生正弦信號 x1(t),x2(t)和x(t)SinSignal_x1= Coeff_LF*sin(2*pi*F1*Time)% x1(t)SinSignal_x2= Coeff_HF*sin(2*pi*F2*Time)% x2(t)SinSignalSum = SinSignal_x1 + SinSignal_x2+ Coeff_NOISE*randn(1, N);SinSignal = SinSignalSum;% x(t)%繪制信號x1(t)figuresubplot(3,1,1)xlabel(時間

15、)xlim(0 0.5)title(時域信號 x1(t)subplot(3,1,2)plot(Time, SinSignal_x1)xlabel(時間)xlim(0 0.05)title(時域信號 x1(t)FFTSinSignal_x1 = 2/N*fft(SinSignal_x1)subplot(3,1,3)FnSin_x1 = abs(FFTSinSignal_x1);FW = 2*pi*(Time)/EndTime;plot(FW(2:length(FW)*FS/pi/2, FnSin_x1(2:length(FW)hold onplot(FW(2:length(FW)-2*pi)*F

16、S/pi/2, FnSin_x1(2:length(FW)xlim(-FS/2 FS/2)xlabel(頻率)%繪制x2(t)figuresubplot(3,1,1)plot(Time, SinSignal_x2)xlabel(時間)xlim(0 0.5)title(時域信號 x(2) subplot(3,1,2)plot(Time, SinSignal_x2)xlabel(時間)xlim(0 0.05)title(時域信號 x2(t)FFTSinSignal_x2 = 2/N*fft(SinSignal_x2)subplot(3,1,3)FnSin_x2 = abs(FFTSinSignal

17、_x2);plot(FW(2:length(FW)*FS/pi/2, FnSin_x2(2:length(FW)hold onplot(FW(2:length(FW)-2*pi)*FS/pi/2, FnSin_x2(2:length(FW) xlim(-FS/2 FS/2)xlabel(頻率)%繪制信號x(t)figuresubplot(3,1,1)plot(Time, SinSignal)xlabel(時間)xlim(0 0.5)title(時域信號)subplot(3,1,2)plot(Time, SinSignal)xlabel(時間)xlim(0 0.05)title(時域信號 x(t

18、)FFTSinSignal = 2/N*fft(SinSignal)subplot(3,1,3)FnSin = abs(FFTSinSignal);plot(FW(2:length(FW)*FS/pi/2, FnSin(2:length(FW)hold onplot(FW(2:length(FW)-2*pi)*FS/pi/2, FnSin(2:length(FW) xlim(-FS/2 FS/2)xlabel(頻率)%設計低通濾波器Rp = 1;%通帶衰減/dBRs = 40; % 阻帶衰減/dBOmegaP1_1=400; %通帶邊界頻率OmegaS1_1=500; % 阻帶邊界頻率FS=2

19、000; %抽樣頻率Wp1=2*pi*OmegaP1_1/FS; %將模擬頻率變為數字角頻率Ws1=2*pi*OmegaS1_1/FS; %將模擬頻率變為數字角頻率OmegaP1=tan(Wp1/2); % 非線性變換OmegaS1=tan(Ws1/2); % 非線性變換Eta_P=OmegaP1/OmegaP1; % 歸一化Eta_S=OmegaS1/OmegaP1; % 歸一化Lemta_P_EquivalentLowPass=1/Eta_P; % 變為等效低通參數Lemta_S_EquivalentLowPass=1/Eta_S; % 變為等效低通參數%估計濾波器階數N, Wn = El

20、lipord(Lemta_P_EquivalentLowPass, Lemta_S_EquivalentLowPass,Rp,Rs,s)；%時間濾波器num1,den1 = ellip(N,Rp,Rs,Wn,s);num2,den2=lp21P(num1,den1,OmegaP1);num,den=bi1inear(num2,den2,0.5);%計算增益響應w = 0:pi/255:pi;h = freqz(num,den,w);g = 20*log10(abs(h);%繪制增益響應figure;plot(w/pi,g);gridaxis(0 1 -60 5);xlabel(omega 八p

21、i); ylabel(增益 /dB);title(橢圓低通濾波器的增益響應);%繪制零點和極點figure;z,p,k=tf2zp(num,den);zplane(z,p);title(零點和極點)figure subplot(2,1,1)YLP = filter(num,den, SinSignal);plot(Time, YLP)xlim(0 0.05)xlabel(時間)title(低通濾波后的時域信號)FFTSinSignalSum = 2 / N * fft(YLP) subplot(2,1,2)FnSinSum = abs(FFTSinSignalSum);plot(FW(2:le

22、ngth(FW)*FS/pi/2, FnSinSum(2:length(FW)hold onplot(FW(2:length(FW)-2*pi)*FS/pi/2, FnSinSum(2:length(FW) xlim(-FS/2 FS/2)ylim(0 1)xlabel(頻率)%returnEnd Program五仿真結果圖時域信號x1(t)10-100.050.10.150.20.250.30.350.40.450.5時間時域信號x1(t) 10-100.0050.010.0150.020.0250.030.0350.04 0.0450.05時間10.50-1000-800-600-400-

23、20002004006008001000頻率時域信號x(2)時間時域信號x2(t)0-1 00.0050.010.0150.020.0250.030.0350.04 0.0450.05時間0.50-1000-800-600-400-20002004006008001000頻率時域信號0-200.050.10.150.20.250.30.350.40.450.5時間時域信號x(t)20-200.0050.010.0150.020.0250.030.0350.04 0.0450.05時間-1000-800-600-400-20002004006008001000頻率 /二10.80.60.40.20-0.2-0.4-0.6-0.8-1零點和極點-1-0.500.51Real Part10.5 0-0.5-1低通濾波器后的時域信號00.0050.010.0150.020.0250.030.0350.04 0.0450.05時間-I0c-1000-800-600-40