畢業(yè)論文matlab實(shí)現(xiàn)iir數(shù)字濾波器的仿真(參考版)
2025-06-22 17:47本頁(yè)面
【正文】 22 徐大鵬，[J].，29 (4) ：1114～1117。16 [J].,(12)：66～69 17 語(yǔ)音低通數(shù)字濾波器[J].,(4)：12～1618 劉衛(wèi)國(guó)，陳昭民，[M].北京：高等教育出版社，2002。7 [M].重慶：重慶大學(xué)出版社,2002。 2 、現(xiàn)在與未來(lái)[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。我也要感謝我的母校嘉應(yīng)學(xué)院，是他提供了良好的學(xué)習(xí)環(huán)境和生活環(huán)境，讓我的大學(xué)生活豐富多姿，為我的人生留下精彩的一筆。在此論文脫稿之際，我再一次向他致以最誠(chéng)摯的謝意。在本次本科生畢業(yè)設(shè)計(jì)的過(guò)程中，程老師給予了我許許多多的關(guān)懷和幫助。首先，我誠(chéng)摯地感謝我的導(dǎo)師趙新穎老師。grid。幅值39。)。xlabel(39。IIR帶通濾波器39。plot(w/pi,Hw2)。endw=(0+eps):pi/300:pi。a0(1)=2*c/(1+omega0)。endfor i=1:K a4(i)=(1+2*omega0*cos(theta(i))+omega0^2)/(12*omega0*cos(theta(i))+omega0^2)。endfor i=1:K a2(i)=2*(2*c^2+1omega0^2)/(12*omega0*cos(theta(i))+omega0^2)。endfor i=1:K G(i)=omega0^2/(12*omega0*cos(theta(i))+omega0^2)。K=floor(N/2)。N=ceil(log(es/ep)/log(omegas/omegap))。ep=sqrt(10^(Ap/10)1)。omegasb=(ccos(wsb))/sin(wsb)。omegap=abs((ccos(wpb))/sin(wpb))。wsb=2*pi*fsb/fs。wpb=2*pi*fpb/fs。As=。fsb=。fpb=4。fs=20。grid。幅值39。)。xlabel(39。IIR高通濾波器39。plot(w/pi,Hw2)。endw=(0+eps):pi/300:pi。a0(1)=2*c/(1+omega0)。endfor i=1:K a4(i)=(1+2*omega0*cos(theta(i))+omega0^2)/(12*omega0*cos(theta(i))+omega0^2)。endfor i=1:K a2(i)=2*(2*c^2+1omega0^2)/(12*omega0*cos(theta(i))+omega0^2)。endfor i=1:K G(i)=omega0^2/(12*omega0*cos(theta(i))+omega0^2)。K=floor(N/2)。N=ceil(log(es/ep)/log(omegas/omegap))。ep=sqrt(10^(Ap/10)1)。omegasb=(ccos(wsb))/sin(wsb)。omegap=abs((ccos(wpb))/sin(wpb))。wsb=2*pi*fsb/fs。wpb=2*pi*fpb/fs。As=。fsb=。fpb=4。fs=20。plot(w/pi,Hw2)。endw=(0+eps):pi/300:pi。endif K(N/2) G0=omega0/(omega0+1)。endfor i=1:K a1(i)=2*(omega0^21)/(12*omega0*cos(theta(i))+omega0^2)。for i=1:K theta(i)=pi*(N1+2*i)/(2*N)。omega0=omegap/ep^(1/N)。es=sqrt(10^(As/10)1)。omegas=cot(ws/2)。ws=2*pi*fstop/fs。As=10。fstop=4。fs=20。 IIR高通濾波的仿真目的：用程序編寫(xiě)一個(gè)低通濾波器，仿真圖形如圖431和圖432所示：圖431 IIR高通濾波器的仿真圖432 IIR高通濾波器的仿真由圖可以看出：， 頻率之前迅速衰減。)。ylabel(39。頻率39。)。title(39。Hw2=1./(1+(tan(w/2)/omega0).^(2*N))。a0=(omega01)/(omega0+1)。endfor i=1:K a2(i)=(1+2*omega0*cos(theta(i))+omega0^2)/(12*omega0*cos(theta(i))+omega0^2)。endfor i=1:K G(i)=omega0^2/(12*omega0*cos(theta(i))+omega0^2)。K=floor(N/2)。N=ceil(log(es/ep)/log(omegas/omegap))。ep=sqrt(10^(Ap/10)1)。omegap=tan(wp/2)。wp=2*pi*fpass/fs。Ap=。fpass=4。程序代碼如下：clear。 IIR低通濾波器的仿真目的：用程序編寫(xiě)一個(gè)低通濾波器。T/ms39。)。title(39。振幅39。plot(x)。x=sin(2*pi*Frq*t)。 %信號(hào)頻率值fs=1024。clear all。：仿真圖形如圖411和圖412所示：圖411 MATLAB正弦波的仿真圖412 MATLAB正弦波的仿真 如上圖所示，通過(guò)編寫(xiě)程序進(jìn)行正弦波的仿真，可以清楚的看吹正弦波的參數(shù)情況，下面進(jìn)行IIR各種濾波器的仿真。在此，我們只對(duì)扯關(guān)系仿真進(jìn)行說(shuō)明。這兩種仿真模式的效果是一樣的，都用同一個(gè)軟件。根據(jù)這個(gè)描述就可以分析其頻率特性和相位特性，以驗(yàn)證設(shè)計(jì)結(jié)果是否滿足指標(biāo)要求；或者利用計(jì)算機(jī)仿真實(shí)現(xiàn)設(shè)計(jì)的濾波器，再分析濾波結(jié)果來(lái)判斷。這是因?yàn)槔硐氲念l率響應(yīng)是不可能實(shí)現(xiàn)的，由于它的幅度響應(yīng)在頻帶之間是突變的，因而其單位抽樣響應(yīng)是非因果的、不可實(shí)現(xiàn)的，只能對(duì)其逼近。在很多實(shí)際應(yīng)用中，例如語(yǔ)音或音頻信號(hào)處理中，數(shù)字濾波器常用來(lái)實(shí)現(xiàn)選頻操作。數(shù)字濾波器技術(shù)指標(biāo)指標(biāo)參數(shù)變換模擬濾波器技術(shù)指標(biāo)相應(yīng)模擬濾波器設(shè)計(jì)模擬濾波器離散化數(shù)字濾波器圖341 數(shù)字濾波器設(shè)計(jì)步驟在設(shè)計(jì)一個(gè)濾波器之前，必須首先確定一些技術(shù)指標(biāo)。簡(jiǎn)而言之：IIR的技術(shù)指標(biāo)→模擬濾波器的技術(shù)指標(biāo)→設(shè)計(jì)模擬低通濾波器H(s) →IIR數(shù)字低通濾波器H1(z) →IIR數(shù)字濾波器數(shù)字濾波器的設(shè)計(jì)步驟如圖341所示。因此，充分利用這些已有的資源將會(huì)給數(shù)字濾波器的設(shè)計(jì)帶來(lái)很大的方便，IIR 數(shù)字濾波器的設(shè)計(jì)步驟是：（1）按一定的規(guī)則將給出的數(shù)字濾波器的技術(shù)指標(biāo)轉(zhuǎn)換為模擬濾波器的技術(shù)指標(biāo)；（2）根據(jù)轉(zhuǎn)換后的技術(shù)指標(biāo)設(shè)計(jì)模擬低通濾波器。目
電大資料相關(guān)推薦
鄂ICP備17016276號(hào)-1