ex18模擬原型濾波器的設(shè)計分解

1、第十次實驗實驗名稱:學生班級:學生姓名:學生學號:指導教師:模擬原型濾波器的設(shè)計電信zgx、實驗目的(1)加深對模擬濾波器基本類型、特點和主要設(shè)計指標的了解(2)掌握模擬低通濾波器原型的設(shè)計方法(3)學習MATLA前言有關(guān)模擬原型濾波器設(shè)計的子函 數(shù)的使用 方法二、實驗原理輸入信號和輸出信號均為連續(xù)時間信號，沖擊響應也是連續(xù)的濾 波器，成為模擬濾波器。模擬濾波器從功能上可以分為低通、高通、帶通、帶阻以及全通濾波器。實際使用中理想濾波器是不可能實現(xiàn)的，必須設(shè)計一個因果可實現(xiàn)的濾波器去逼近。通常通帶和阻帶都允許存在一定誤差容限， 即通 帶不一定完全水平，阻帶也不一定絕對衰減到零。 在通帶和阻帶之間

2、 允許存在一定寬度的過渡帶。三、實驗任務(wù)設(shè)計一個模擬原型低通濾波器，通帶截止頻率fp=6kHz,通帶最大衰減RpW 1dB,阻帶截止頻率fs=15kHz,阻帶最小衰減 As三30dR 要求：分別實現(xiàn)符合以上指標的巴特沃斯濾波器、切比雪夫一型濾波器、切比雪夫二型濾波器、橢圓濾波器，繪制幅頻特性和相頻特 性曲線、零極點分布圖，并列寫傳遞函數(shù)表達式。4 6I391011121314151617131920212223四、實驗過程和結(jié)果1,巴特沃斯濾波器:程序清單:P Editor - C:UsersMQominXDesktopXS,PU V 幅 7t1 ,m 4 fp=600 0：0mgp=2*p

3、i *fp; fs-15000 ;0i&g5-2*pi*fs;Kp=l;As=3Q;n, OftgcZ=buttotcKOngp, Oings, Rp, As, s ); ZzO* pO, kO=buttap (n)： bO=kO*real(poly(zO): aO-real(poly(pO);Pbpoly2str(bO,? s!)Pa=poly2str(aO,n sn) R, Ong =freqs (bO, aO); subplot (3,1,1). plot (Oing*0rogc/<2*pi), abs(H), grid is (0T 20000,-0,5, L5); Wa

4、beLC幅際特性,)； subplot (3, 2)* plot (0mg*0mgc/ (2*pi), angle(H)3 grid axis (0,20000,-5t 5); yl獨el ('相赫特性)；0mgx0=10mgp, Omgsl/Omgc;Hx=freqs (b0» aQ, OngO); dbHK=-20*logl0(abs(Hx)/max (abs (H)一 subplot(3, 1,3), plat(pO*Omgc, x ); a?c is square, axis equal, grid on ylabeL C零極點分布圖)；|得到輸出:Pb =Pas&q

5、uot;3 + 3-2361 s"4 + 5.2361 s,3 + 5.2361 s,2 + 3,2361 s + 1所以 H(s)=s0 +3.2361s： +5.2361s； +5.2361s；十3.2361so + 1且所求曲線:華實憔s尊矍麋理x 102.切比雪夫一型濾波器:程序清單：'Editor - CAUsersXMaominXDesktopXa.t1 ,m tNm* ，+得到數(shù)據(jù):Pb =0.24565Pa =4 + Q,前281 s"3 + 1, 4539 /2 + 0, 74262 s + 0.27563所以傳遞函數(shù)表達式為：1H(S) 32s0

6、 十0.95281so 十1.4539s。十0.74262so 十0.27563所求曲線圖形：O到舜盛要3 .切比雪夫二型濾波器：程序清單：,Editor - C:UsersMoominDe&ktop37t. 一pub. v 姆通唯 4 O 1 ®得到數(shù)據(jù):Pb =0. 031623 s 4 + 0. 25298 /2 + 0. 25298Pa =s"4 + 1. 7643 s"3 + 1. 5653 s"2 + 0.314S6 s + 0-2529E所以傳遞函數(shù)為:H0.031623§ +0.25298s； +0.25298 一 s4

7、 1.7648s0 1.5653s： 0.81486so 0.25298且所求曲線為:史非糜要1x104Pole-Zero Map4-2024Real Axis (seconds'1)Q熊通/白加沏絲： 。 產(chǎn). <11 J爐加7a威狗.?24 .橢圓濾波器:程序清單：Editor - CAU&ersXMoominXDesktopXS%., 一ED. PUB.V& o| 0 t1 ,m t2,m t3,m t4,m* 貓 +1-fp=6000 ：0nigp=2*pi*fp;2-fs-15000 :0i&gs=2*pi*fs :3 -Rp=l: As-304

8、 -In, OagnI=ellipord(Oagp, OagsKp, As, s7);5 pO,kOZ-ellipap(n Rpt As)6 -b0=k0*real(poly(zO):-aO=real(poly(pO)5-Pb-poly2str (b0T s")9-Pa=poly2str(a0T sn)1011 bl-Zzeros(1, length(aO) -length (bO) s bO'：12 - -sas, g'=tf2sos (bl aO)13 -H, 0mg>freqs(b0, a0);14 - dbH-20*logl0 (tabs (H)+ep5

9、)/raax(ab5(H)'1516 - subplot (2» 2t 1), plot (0mg*0mEn/(2*pi), abs(H), grid17 axis (0; 3000Q, -0,1,1,1);18 -yl曲式(幅瀕特性)；19 -wubpl0t2,2, 2), plot (0mg*0m匕, angle 田)t/rid2Q -a?: is (10, 40000, -4, 41);21-ylabeM, 相頻特性,)；22 -宕ubplot(2, 2, Mh plot (Omg+Omgn/(2竄pi), dbH) , grid23-WabmlC 幅度(dB)'

10、;);24 - subplot (2« 2t 4) * pzaop (bO* aO):25-a?! is square, axis equal, grid on26-ylabeL (零極點分布圖)；27OmgjjO=lOmgp, OagsZ/OiDgn；28-Hx=freqs(bO, a0±OngxO)：29 - dbHx=-20*log (abs (Hx) /max tabs (E)得到數(shù)據(jù):Pb -0.14899 屋2 + 0,56863Pa =s'3 + 0, &7O12 屋2 + 1.2459 s + 0.563632所以傳遞函數(shù)為：H(s)0.14

11、89呢0.56863S3 +0.9701250 +1.2459so + 0.56863且所求曲線為：2邙名醫(yī)矍086 4 o.O 翅袋轅森Pole-Zero Map讓氏。耳。我萬4萬；2 T1a 相 r:二cx勺 Q.ga /2.5：牙1每-1尊掌-；!，一年公七 W11. v, Jr"lrx不2*4。、F; U 詼1-2/0Real Axis (secondsi1)r若匚gs囤將會磔重防o O -24 mp函瞿6020O 80五、實驗思考題1 .什么是模擬濾波器？模擬濾波器是如何分類的？什么是模擬原型濾波器？答：1）輸入信號和輸出信號均為連續(xù)時間信號，沖擊響應也是 連續(xù)的濾波器，成

12、為模擬濾波器。2）模擬濾波器分為有源濾波器和無源濾波器。3）在已知的低通模擬濾波器的基礎(chǔ)上設(shè)計的數(shù)字濾波器稱 為模擬濾波器。2 .試歸納各類模擬濾波器幅頻特性曲線的基本特點。答：1）巴特沃斯濾波器具有單點下降的幅頻特性曲線，通帶和阻帶幅頻都比較平坦，且 N階數(shù)越高，曲線越陡峭。2）切比雪夫一型濾波器的幅頻特性曲線可以看出， N階數(shù) 越低，過渡帶曲線越平緩；系統(tǒng)階數(shù) N越高，過渡帶曲線越陡 峭，在通帶區(qū)間呈現(xiàn)波動狀態(tài)。3）切比雪夫二型濾波器的幅頻特性曲線可以看出， N階數(shù) 越低，過渡帶曲線越平緩；系統(tǒng)階數(shù) N越高，過渡帶曲線越陡 峭，在阻帶區(qū)間呈現(xiàn)波動狀態(tài)。4）橢圓濾波器的幅頻特性曲線可以看出，

13、 N階數(shù)越低，過 渡帶曲線越平緩；系統(tǒng)階數(shù)N越高，過渡帶曲線越陡峭，在阻 帶和通帶區(qū)間均呈現(xiàn)波動狀態(tài)。六、實驗感想通過這次的實驗，使我學到了不少實用的知識，更重要的是，做實 驗的過程，思考問題的方法，這與做其他的實驗是通用的，真正使我們 受益匪淺。1f p=6000 ： 0agp=2*p i *f p ：2-fs-15000 ;0i&gs-2*pi*fs;3-Rp=l; As=30 ;4 -n, Onsn=cheblord(Omgp Omgs, Kp, Ast s* * * * S * 7 * * * * * 13);5ZjO* pO, kOZ=cheblap (n»Rp)：

14、6-bO=kO*real (po(zO);:-aO=real(poly (pO);S-Pb-poly2str (bOT sh)9-Pa=poly2str(aO, s?)1011 blal2eros (1, length(aO)-length(bO) t bO；12 - "sas3 g"=tf2sos (bl, atl)13 - Omg-freqs (bO, aO);14 - dbH-20*logl0 (abs (H)+eps)/max (abs(H);1516 - subplot (2t 2t 1), plot (0mg*0mgn/ (2*pi), abs (H), gri

15、d17axis (0, 30000, -0. I, 1. 11)：13-¥1獨elC幅猱特性)；19 - subplot (2, 2f 2)t plot (0mg*0mgn/(2*pi), angle(HJ)T grid20-a；;is (0,30000,-4, 4)：21-yl金式相頻特性)；22 - subplot (2* 2, 3), plot (0mg+0mgn/ (2*pi), abs(H), grid23 - axis (10, 30000T 0. 1,1.1);24 - ylabel('幅度(dE1');25 - subplot (2, 2r 4), p

16、lot IpOiBgn, ;26 - is square* axis equaL grid on27 -yl&bml零極點分布圖)：2S29 -Omg?;fl= "Oingp, OmgsJ/Oingn：3Q - Hx=f reqs (b0 aOOmgxO);31 - dblb:=-20*log (abs (Ib:)/Dai: labs (H)tLm tZm t3,m*+mJIJ.1 -fp=6Q00 ;Oaigp=2*pi*fp ;2 -fs-15000:0mg5-2*pi*fs;3 -Rp=l;As=30;4 -nf Oagn = cheb2ord(Ongp0ngs, Rp

17、, As,? s"); ZzO* pO* 1 2 3 4 S *, kOZ-cheb2ap (n, As);6 -bO=kO*real(poly(zO)7 -aO-real(poly(pO)S -Pb-p01y2str(bO, sh)9 -Pa=poly2str(aQ,7 * * * * * 13 s)1011 -bl= zercs(lt length(a0)-length (bO), bO：12 - Isos, g；-tf2sos(bl, a0)13 -Owg -f rqs (bQ, aO);14 - dbH=20*Logl0(abs(H)+eps)/nax(abs(H1);1516 - subplot (2* 2, 1) s plat (0mg*0mgn/ (2*pi)* abs (H), grid17 - axis (0, 30000,-0. L, 1.1):13 -WabalC幅猱特性)：19 - subplot (2,2, 2)T plot (0mg*0mgn/ (2*pi), angle(H) P grid20 -axi£(0, 40000,7, 41);21 -yl