實(shí)驗(yàn)十一IIR數(shù)字濾波器設(shè)計(jì)及應(yīng)用

1、數(shù)字信號(hào)處理及實(shí)驗(yàn)實(shí)驗(yàn)報(bào)告實(shí)驗(yàn)題目IIR數(shù)字濾波器設(shè)計(jì)與應(yīng)用姓名組別班級(jí)光電14學(xué)號(hào)【實(shí)驗(yàn)?zāi)康摹考由罾斫釯IR數(shù)字濾波器的特性，掌握IIR數(shù)字濾波器的設(shè)計(jì)原理與設(shè)計(jì)方法，以及IIR數(shù)字濾波器的應(yīng)用?！緦?shí)驗(yàn)原理】IIR數(shù)字濾波器一般為線性時(shí)不變的因果離散系統(tǒng)，N階IIR數(shù)字濾波器的系統(tǒng)函數(shù)可以表達(dá)為z-1的有理多項(xiàng)式: Hz=j=0Mbjz-j1+i=1Naiz-i=b0+b1z-1+bM-1z-(M-1)+bMz-M1+a1z-1+aN-1z-N-1+aNz-N系數(shù)ai；i=1，2，N中至少有一個(gè)非零。對(duì)于因果IIR數(shù)字濾波器，滿足MN。IIR數(shù)字濾波器的設(shè)計(jì)主要通過成熟的模擬濾波器設(shè)計(jì)方法來

2、實(shí)現(xiàn)。首先在頻域?qū)?shù)字濾波器設(shè)計(jì)指標(biāo)轉(zhuǎn)換為模擬濾波器設(shè)計(jì)指標(biāo)，將任意模擬濾波器轉(zhuǎn)換為原型模擬低通濾波器設(shè)計(jì)指標(biāo)，根據(jù)模擬低通濾波器設(shè)計(jì)指標(biāo)設(shè)計(jì)出相應(yīng)的模擬濾波器H(s)，再經(jīng)過脈沖響應(yīng)不變法或雙線性變換法得到所需的IIR數(shù)字濾波器H(z)。由此可見，IIR數(shù)字濾波器設(shè)計(jì)的重要環(huán)節(jié)是模擬低通濾波器的設(shè)計(jì)，主要包括Butterworth、Chebyshev和橢圓等濾波器。 MATLAB 信號(hào)處理工具箱中提供了IIR濾波器設(shè)計(jì)的函數(shù)。常用設(shè)計(jì)函數(shù)有： IIR濾波器階數(shù)選擇 buttord -巴特沃斯(Butterworth)濾波器階數(shù)選擇。 cheb1ord -切比雪夫(Chebyshev)I 型濾

3、波器階數(shù)選擇。 cheb2ord -切比雪夫(Chebyshev)II 型濾波器階數(shù)選擇。 ellipord -橢圓(Elliptic)濾波器階數(shù)選擇。 IIR 濾波器設(shè)計(jì) butter -巴特沃斯(Butterworth)濾波器設(shè)計(jì) cheby1 -切比雪夫(Chebyshev)I 型濾波器設(shè)計(jì) cheby2 -切比雪夫(Chebyshev)II 型濾波器設(shè)計(jì) ellip -橢圓(Elliptic)濾波器設(shè)計(jì) maxflat -通用的巴特沃斯(Butterworth)低通濾波器設(shè)計(jì)yulewalk -Yule-Walker 濾波器設(shè)計(jì)(直接數(shù)字濾波器設(shè)計(jì)法)1.Butterworth濾波器設(shè)

4、計(jì) Butterworth濾波器是通帶、阻帶都單調(diào)衰減的濾波器。 (1)調(diào)用buttord函數(shù)確定巴特沃斯濾波器的階數(shù)，格式為 N,Wc = buttord(Wp,Ws,Ap,As) 輸入?yún)?shù)：Ap,As為通帶最大衰減和阻帶最小衰減，以dB為單位。 Wp,Ws為歸一化通帶截頻和阻帶截頻，歸一化頻率=實(shí)際數(shù)字頻率。因此0<Wp，Ws<1，1對(duì)應(yīng)數(shù)字頻率（rad）、模擬域抽樣頻率的一半fsam2hz。輸出參數(shù)：N為濾波器的階數(shù)；Wc為截頻，0 < Wc < 1。(2)調(diào)用butter函數(shù)設(shè)計(jì)出巴特沃斯濾波器，格式為 b,a = butter(N,Wc,options) 輸入

5、參數(shù)：N和Wc是buttor函數(shù)返回的參數(shù)，含義見上。 Options=low, high, bandpass, stop, 分別對(duì)應(yīng)低通、高通、帶通、帶阻，默認(rèn)情況下為低通或帶通。 輸出參數(shù)：b和a為設(shè)計(jì)出的IIR數(shù)字濾波器H(s)的分子多項(xiàng)式和分母多項(xiàng)式的系數(shù)矩陣。 2.Chebyshev I型濾波器設(shè)計(jì) Chebyshev I型濾波器為通帶紋波控制器：在通帶呈現(xiàn)紋波特性，在阻帶單調(diào)衰減。 N,Wc = cheb1ord(Wp, Ws, Ap, As)b,a = cheby1(N,Ap,Wc,options) 參數(shù)含義與butter中參數(shù)一致。3.Chebyshev II 型濾波器設(shè)計(jì) C

6、hebyshev II 型濾波器為阻帶紋波控制器:在阻帶呈現(xiàn)紋波特性。在通帶單調(diào)衰減。 N,Wc = cheb2ord(Wp, Ws, Ap, As) b,a = cheby2(N,As,Wc,options) 4.橢圓濾波器設(shè)計(jì) 橢圓濾波器在通阻帶都呈現(xiàn)紋波特性。 N,Wc = ellipord(Wp,Ws,Ap,As) b,a = ellip(N,Ap,As,Wc,options)【實(shí)驗(yàn)結(jié)果與數(shù)據(jù)處理】1、信號(hào)xk=1+cos(4k)+cos(23k) ,確定設(shè)計(jì)指標(biāo)，實(shí)現(xiàn)各種IIR數(shù)字濾波器以實(shí)現(xiàn)以下信號(hào)處理。 （1）設(shè)計(jì)IIR低通濾波器，濾除cos(23k)的成分。MATLAB程序：c

7、lc,clear,close allN,Wc=buttord(0.250,0.677,3,60);b,a=butter(N,Wc,'low');freqz(b,a);axis(0,1,-120,0);title('巴特沃斯低通數(shù)字濾波器')N=5,Wc=0.2702,a =1.0000 -2.2729 2.4660 -1.4359 0.4450 -0.0575b =0.0045 0.0226 0.0452 0.0452 0.0226 0.00452、在“電話撥號(hào)音合成和識(shí)別”實(shí)驗(yàn)中，DTMF信號(hào)是利用譜分析方法實(shí)現(xiàn)解碼。在此討論另一種利用IIR數(shù)字濾波器的解碼方

8、法。如圖所示，接收到的音頻信號(hào)通過低通和高通兩個(gè)濾波器。由于DTMF信號(hào)行頻低于列頻，所以低通和高通可以分離出行和列的信息。再通過帶通濾波器，對(duì)其輸出的振幅大小進(jìn)行判斷，可以確定具體的行頻和列頻。（1）利用DTMF的產(chǎn)生原理合成數(shù)字4，合成公式為 x4k=sin2×770k8000+sin(2×1209k8000)式中：k=0:204;8000為抽樣頻率（Hz）。畫出其時(shí)域波形和頻譜圖。（2）根據(jù)圖中的低通和高通濾波器，畫出這兩個(gè)濾波器的頻譜響應(yīng)；并利用y=filter(b,a,x)對(duì)信號(hào)x4k進(jìn)行濾波；利用fft函數(shù)求得輸出的頻響，繪出幅頻響應(yīng)曲線。MATLAB程序：cl

9、c,clear,close allT0=204;N=205;k=0:T0;x4=sin(2*pi*770*k/8000)+sin(2*pi*1209*k/8000);subplot(2,1,1)plot(k,x4);title('時(shí)域波形')X=fft(x4,N);subplot(2,1,2)plot(k-N/2,abs(fftshift(X);a .低通濾波：clc,clear,close allfigure(1)N,Wc=buttord(0.1925,0.30225,3,60);b,a=butter(N,Wc,'low');freqz(b,a);axis(0

10、,1,-120,0);figure(2)T0=204;N=205;k=0:T0;x4=sin(2*pi*770*k/8000)+sin(2*pi*1209*k/8000);subplot(3,1,1)plot(k,x4);title('時(shí)域波形')X=fft(x4,N);xlabel('f(Hz)');ylabel('magnitude');title('幅度譜')y=filter(b,a,x4);subplot(3,1,2)plot(k,y);title('對(duì)信號(hào)進(jìn)行低通濾波')Y=fft(y,N);subplo

11、t(3,1,3)plot(k-N/2,abs(fftshift(Y);title('低通濾波后的頻譜圖')b.高通濾波：clc,clear,close allfigure(1)N,Wc=buttord(0.1925,0.30225,3,60);b,a=butter(N,Wc,'high');freqz(b,a);axis(0,1,-120,0);figure(2)T0=204;N=205;k=0:T0;x4=sin(2*pi*770*k/8000)+sin(2*pi*1209*k/8000);subplot(3,1,1)plot(k,x4);title('

12、;時(shí)域波形')X=fft(x4,N);xlabel('f(Hz)');ylabel('magnitude');title('幅度譜')y=filter(b,a,x4);subplot(3,1,2)plot(k,y);title('對(duì)信號(hào)進(jìn)行高通濾波')Y=fft(y,N);subplot(3,1,3)plot(k-N/2,abs(fftshift(Y);title('高通濾波后的頻譜圖')【實(shí)驗(yàn)結(jié)論與分析】1、巴特沃思、切比雪夫I型，切比雪夫II型和橢圓模擬圓型濾波器的幅頻特性各有哪些特點(diǎn)？其優(yōu)缺點(diǎn)是什么？答：巴特沃斯：3dB點(diǎn)不變性，通帶內(nèi)有最大平坦的幅度特性，單調(diào)減小 ，過 渡帶及阻帶內(nèi)快速單調(diào)減小 ；切比雪夫I 型：截止頻率，不一定為3dB帶寬，通帶內(nèi)有等間隔的波紋。切比雪夫II 型：阻帶內(nèi)有波紋