畢業(yè)論文matlab實現(xiàn)iir數(shù)字濾波器的仿真(文件)
【正文】 R高通濾波器的仿真由圖可以看出：， 頻率之前迅速衰減。fstop=4。ws=2*pi*fstop/fs。es=sqrt(10^(As/10)1)。for i=1:K theta(i)=pi*(N1+2*i)/(2*N)。endif K(N/2) G0=omega0/(omega0+1)。plot(w/pi,Hw2)。fpb=4。As=。wsb=2*pi*fsb/fs。omegasb=(ccos(wsb))/sin(wsb)。N=ceil(log(es/ep)/log(omegas/omegap))。endfor i=1:K G(i)=omega0^2/(12*omega0*cos(theta(i))+omega0^2)。endfor i=1:K a4(i)=(1+2*omega0*cos(theta(i))+omega0^2)/(12*omega0*cos(theta(i))+omega0^2)。endw=(0+eps):pi/300:pi。IIR高通濾波器39。)。grid。fpb=4。As=。wsb=2*pi*fsb/fs。omegasb=(ccos(wsb))/sin(wsb)。N=ceil(log(es/ep)/log(omegas/omegap))。endfor i=1:K G(i)=omega0^2/(12*omega0*cos(theta(i))+omega0^2)。endfor i=1:K a4(i)=(1+2*omega0*cos(theta(i))+omega0^2)/(12*omega0*cos(theta(i))+omega0^2)。endw=(0+eps):pi/300:pi。IIR帶通濾波器39。)。grid。在本次本科生畢業(yè)設(shè)計的過程中，程老師給予了我許許多多的關(guān)懷和幫助。我也要感謝我的母校嘉應(yīng)學(xué)院，是他提供了良好的學(xué)習(xí)環(huán)境和生活環(huán)境，讓我的大學(xué)生活豐富多姿，為我的人生留下精彩的一筆。7 [M].重慶：重慶大學(xué)出版社,2002。22 徐大鵬，[J].，29 (4) ：1114～1117。16 [J].,(12)：66～69 17 語音低通數(shù)字濾波器[J].,(4)：12～1618 劉衛(wèi)國，陳昭民，[M].北京：高等教育出版社，2002。 2 、現(xiàn)在與未來[J].，13(4)：48～493 [J].，(02)：1～44 李洪偉，、IIR數(shù)字濾波器的應(yīng)用及比較[J].，(15)：79～80 5 Oppenheim A V,Schafer R signal processing[M].Englewood Cliffs,NJ：PrenticeHall,1975。在此論文脫稿之際，我再一次向他致以最誠摯的謝意。首先，我誠摯地感謝我的導(dǎo)師趙新穎老師。幅值39。xlabel(39。plot(w/pi,Hw2)。a0(1)=2*c/(1+omega0)。endfor i=1:K a2(i)=2*(2*c^2+1omega0^2)/(12*omega0*cos(theta(i))+omega0^2)。K=floor(N/2)。ep=sqrt(10^(Ap/10)1)。omegap=abs((ccos(wpb))/sin(wpb))。wpb=2*pi*fpb/fs。fsb=。fs=20。幅值39。xlabel(39。plot(w/pi,Hw2)。a0(1)=2*c/(1+omega0)。endfor i=1:K a2(i)=2*(2*c^2+1omega0^2)/(12*omega0*cos(theta(i))+omega0^2)。K=floor(N/2)。ep=sqrt(10^(Ap/10)1)。omegap=abs((ccos(wpb))/sin(wpb))。wpb=2*pi*fpb/fs。fsb=。fs=20。endw=(0+eps):pi/300:pi。endfor i=1:K a1(i)=2*(omega0^21)/(12*omega0*cos(theta(i))+omega0^2)。omega0=omegap/ep^(1/N)。omegas=cot(ws/2)。As=10。fs=20。)。頻率39。title(39。a0=(omega01)/(omega0+1)。endfor i=1:K G(i)=omega0^2/(12*omega0*cos(theta(i))+omega0^2)。N=ceil(log(es/ep)/log(omegas/omegap))。omegap=tan(wp/2)。Ap=。程序代碼如下：clear。T/ms39。title(39。plot(x)。 %信號頻率值fs=1024。：仿真圖形如圖411和圖412所示：圖411 MATLAB正弦波的仿真圖412 MATLAB正弦波的仿真 如上圖所示，通過編寫程序進(jìn)行正弦波的仿真，可以清楚的看吹正弦波的參數(shù)情況，下面進(jìn)行IIR各種濾波器的仿真。這兩種仿真模式的效果是一樣的，都用同一個軟件。這是因為理想的頻率響應(yīng)是不可能實現(xiàn)的，由于它的幅度響應(yīng)在頻帶之間是突變的，因而其單位抽樣響應(yīng)是非因果的、不可實現(xiàn)的，只能對其逼近。數(shù)字濾波器技術(shù)指標(biāo)指標(biāo)參數(shù)變換模擬濾波器技術(shù)指標(biāo)相應(yīng)模擬濾波器設(shè)計模擬濾波器離散化數(shù)字濾波器圖341 數(shù)字濾波器設(shè)計步驟在設(shè)計一個濾波器之前，必須首先確定一些技術(shù)指標(biāo)。因此，充分利用這些已有的資源將會給數(shù)字濾波器的設(shè)計帶來很大的方便，IIR 數(shù)字濾波器的設(shè)計步驟是：（1）按一定的規(guī)則將給出的數(shù)字濾波器的技術(shù)指標(biāo)轉(zhuǎn)換為模擬濾波器的技術(shù)指標(biāo)；（2）根據(jù)轉(zhuǎn)換后的技術(shù)指標(biāo)設(shè)計模擬低通濾波器。這是很好的性質(zhì)。與FIR濾波器的設(shè)計不同，IIR濾波器設(shè)計時的階數(shù)不是由設(shè)計者指定，而是根據(jù)設(shè)計者輸入的各個濾波器參數(shù)（截止頻率、通帶濾紋、阻帶衰減等），由軟件設(shè)計出滿足這些參數(shù)的最低濾波器階數(shù)。 IIR數(shù)字濾波器在設(shè)計上可以借助成熟的模擬濾波器的成果，如巴特沃斯、契比雪夫和橢圓濾波器等，有現(xiàn)成的設(shè)計數(shù)據(jù)或圖表可查，其設(shè)計工作量比較小，對計算工具的要求不高。 IIR濾波器的設(shè)計步驟IIR數(shù)字濾波器的系統(tǒng)函數(shù)可以寫成封閉函數(shù)的形式。不同的是模擬濾波器的設(shè)計是在S平面上用數(shù)學(xué)逼近法去尋找近似的所需特性H(S),而數(shù)字濾波器則是在Z平面尋找合適的H(z)。對照模擬濾波器的傳遞函數(shù):有了雙線性變換，模擬濾波器的數(shù)字化只須用進(jìn)行置換。這樣s 平面與z 平面是一一對應(yīng)關(guān)系, 消除了多值變換性, 也就消除了頻譜混疊現(xiàn)象。無限沖激響應(yīng)濾波器的單位抽樣響應(yīng)h(n)是無限長的，其差分方程如(32)式所示，是遞歸式的，即結(jié)構(gòu)上存在著輸出信號到輸入信號的反饋，其系統(tǒng)函數(shù)具有(31)式的形式，因此在z平面的有限區(qū)間(0︱z︱∞)有極點存在。圖322 串聯(lián)型（3）IIR濾波器的并聯(lián)型結(jié)構(gòu)，如圖323所示：優(yōu)點： ① 簡化實現(xiàn)，用一個二階節(jié)，通過變換系數(shù)就可實現(xiàn)整個系統(tǒng)； ② 極、零點可單獨控制、調(diào)整； ③ 各二階零、極點的搭配可
電大資料相關(guān)推薦
鄂ICP備17016276號-1