




版權(quán)說(shuō)明:本文檔由用戶(hù)提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡(jiǎn)介
1、實(shí)驗(yàn)五 有限長(zhǎng)單位脈沖響應(yīng)濾波器設(shè)計(jì)1、 實(shí)驗(yàn)?zāi)康?、 掌握用窗函數(shù)法、頻率采樣法以及優(yōu)化設(shè)計(jì)法設(shè)計(jì)FIR濾波器的原理及方法,熟悉相應(yīng)的MATLAB編程。2、 熟悉線(xiàn)性相位FIR濾波器的幅頻特性和相頻特性。3、 了解各種不同窗函數(shù)對(duì)濾波器性能的影響。 2、 實(shí)驗(yàn)原理window=ones(1, N): 產(chǎn)生N點(diǎn)矩形窗,行向量。window=hann(N): 產(chǎn)生N點(diǎn)漢寧窗,列向量。window=hanning(N): 產(chǎn)生N點(diǎn)非零漢寧窗,列向量。等價(jià)于去除hann(N+2)的第一個(gè)零元素和最后一個(gè)零元素,得到的N點(diǎn)非零窗函數(shù)。window=hamming(N): 產(chǎn)生N點(diǎn)海明窗,列向量。win
2、dow=blackman(N): 產(chǎn)生N點(diǎn)布萊克曼窗,列向量。window=kaiser(N, beta): 產(chǎn)生參數(shù)為beta的N點(diǎn)凱塞窗,列向量。M, Wd, beta, ftype=kaiserord(f, a, dev, fs): 凱塞窗參數(shù)估計(jì)。f為一組邊界頻率,最高頻率為fs/2。a為f中各個(gè)頻帶的幅度值,通帶取1,阻帶取0。如果f中有2個(gè)元素,則形成3個(gè)頻帶,其中第1個(gè)和第3個(gè)是通帶或阻帶,第2個(gè)是過(guò)渡帶,a中也有2個(gè)元素,指明第1個(gè)和第3個(gè)頻帶是通帶還是阻帶;如果f中有4個(gè)元素,則形成5個(gè)頻帶,其中1,3和5是通帶或阻帶,2和4是過(guò)渡帶,a中有3個(gè)元素,指明1,3和5是通帶還是
3、阻帶。dev的維數(shù)與a相同,指明每個(gè)頻帶上的波動(dòng)值。fs為采樣頻率。M為FIR濾波器的階數(shù),M=N-1。Wd為歸一化邊界頻率,等于數(shù)字邊界角頻率除以,或者邊界頻率除以fs/2。beta就是凱塞窗的參數(shù)。ftype為濾波器的類(lèi)型。b = fir1(M, Wd, 'ftype', window): 用窗函數(shù)法求FIR濾波器的系數(shù)b(單位脈沖響應(yīng))。M為濾波器的階數(shù),M=N-1。Wd為一組歸一化邊界頻率,通帶和阻帶間隔分布,無(wú)過(guò)渡帶;只有一個(gè)元素,表示低通或高通濾波器;有兩個(gè)元素表示帶通和帶阻濾波器;有三個(gè)及以上元素,表示多帶濾波器。'ftype'表示濾波器類(lèi)型,
4、39;high'表示高通濾波器,'stop'表示帶阻濾波器,'DC-0'表示多帶濾波器的第一個(gè)頻帶為阻帶,'DC-1'表示多帶濾波器的第一個(gè)頻帶為通帶。window為窗口類(lèi)型,缺省為海明窗。b = fir2(M, f, m, window): 用頻率采樣法求FIR濾波器的系數(shù)b。M為濾波器的階數(shù),M=N-1。f為一組歸一化頻率,第一個(gè)元素必須為0,最后一個(gè)元素必須為1(對(duì)應(yīng)奈奎斯特頻率,即采樣頻率的一半),中間的元素按升序排列。m的維數(shù)與f相同,指明f中每個(gè)頻率上的理想幅度。window為窗口類(lèi)型,缺省為海明窗。Fir2可以實(shí)現(xiàn)任意幅度特
5、性的濾波器。三、實(shí)驗(yàn)內(nèi)容1、用窗函數(shù)法設(shè)計(jì)一個(gè)線(xiàn)性相位FIR高通濾波器,通帶邊界頻率為0.7,阻帶邊界頻率為0.5,要求阻帶衰減不小于50dB,窗函數(shù)從矩形窗、漢寧窗、海明窗和布萊克曼窗中選取,且要求濾波器的階數(shù)最小。給出該濾波器的單位脈沖響應(yīng)h(n)的解析式,并用MATLAB繪出時(shí)域波形和幅頻特性。2、用窗函數(shù)法設(shè)計(jì)一個(gè)線(xiàn)性相位FIR帶通濾波器,通帶邊界頻率為3kHz和7kHz,阻帶邊界頻率為2kHz和9kHz,采樣頻率為20kHz,要求阻帶衰減不小于40dB,窗函數(shù)從矩形窗、漢寧窗、海明窗和布萊克曼窗中選取,且要求濾波器的階數(shù)最小。給出該濾波器的單位脈沖響應(yīng)h(n)的解析式,并用MATLA
6、B繪出時(shí)域波形和幅頻特性。3、用窗函數(shù)法設(shè)計(jì)一線(xiàn)性相位FIR帶阻濾波器,通帶邊界頻率為600Hz和1400Hz,阻帶邊界頻率為700Hz和1200Hz,采樣頻率為4000Hz,要求阻帶衰減不小于50dB,窗函數(shù)從矩形窗、漢寧窗、海明窗和布萊克曼窗中選取,且要求濾波器的階數(shù)最小。給出該濾波器的單位脈沖響應(yīng)h(n)的解析式,并用MATLAB繪出時(shí)域波形和幅頻特性。4、分別用矩形窗、漢寧窗、海明窗和布萊克曼窗設(shè)計(jì)FIR低通濾波器,理想濾波器的邊界頻率為0.8,N=51,用Fir1函數(shù)設(shè)計(jì)。根據(jù)所設(shè)計(jì)濾波器的阻帶最小衰減值確定實(shí)際阻帶邊界頻率。5、用頻率采樣法設(shè)計(jì)一個(gè)線(xiàn)性相位低通濾波器,N=40,通帶
7、邊界頻率為/4,過(guò)渡帶設(shè)置一個(gè)采樣點(diǎn)|H(k)|=0.5,給出單位脈沖響應(yīng)h(n)的序列值,并作出幅頻響應(yīng)。6、用頻率采樣法設(shè)計(jì)一個(gè)線(xiàn)性相位高通濾波器,通帶邊界頻率為2/3,過(guò)渡帶設(shè)置一個(gè)采樣點(diǎn)|H(k)|=0.5,分別求N=31和N=32時(shí)的單位脈沖響應(yīng)h(n)的序列值,并分別作出幅頻響應(yīng)。7、用頻率采樣法設(shè)計(jì)一個(gè)多帶線(xiàn)性相位濾波器,理想幅頻特性如下圖所示,N=50,過(guò)渡帶設(shè)置一個(gè)采樣點(diǎn),過(guò)渡點(diǎn)的值設(shè)為0.5,窗函數(shù)選矩形窗。給出該濾波器的單位脈沖響應(yīng)h(n)的解析式,并用MATLAB繪出時(shí)域波形和幅頻特性。4、 實(shí)驗(yàn)結(jié)果與分析 1. 解:求通帶和阻帶的數(shù)字邊界頻率:求理想高通濾波器的邊界頻
8、率。理想高通濾波器的邊界頻率是實(shí)際濾波器幅度函數(shù)正負(fù)肩峰頻率的中心,但這里不知道正負(fù)肩峰的頻率,因此用通帶和阻帶邊界頻率的中心近似:選擇窗函數(shù)和窗口長(zhǎng)度。阻帶衰減不小于50dB,因此選擇海明窗(其阻帶衰減為53dB,滿(mǎn)足衰減要求)。濾波器的過(guò)渡帶寬為,因此窗口長(zhǎng)度為:線(xiàn)性相位延遲常數(shù)為:根據(jù)理想邊界頻率和線(xiàn)性相位延遲常數(shù),求理想單位脈沖響應(yīng):窗函數(shù)與理想單位脈沖響應(yīng)相乘,得到線(xiàn)性相位FIR低通濾波器的單位脈沖響應(yīng):時(shí)域波形:幅頻特性:程序:N=33;n=0:1:N-1;wn=0.6*pi;a=(N-1)/2;hd=(sin(pi*(n-a)-sin(n-a)*wn)./(pi*(n-a);hd
9、(a+1)=0.4;win=hamming(N);h=win'.*hd;figure;stem(n,h);xlabel('n');ylabel('h(n)');grid;title('線(xiàn)性相位高通濾波器FIR的單位脈沖響應(yīng)h(n)');H,w=freqz(h,1);H=20*log10(abs(H);figure;plot(w/pi,H);axis(0 1 -100 0);xlabel('頻率/Hz');ylabel('幅度/dB');title('線(xiàn)性相位高通濾波器FIR,海明窗,N=33'
10、;);2.解:求通帶和阻帶的數(shù)字邊界頻率:求理想帶濾波器的邊界頻率。理想帶通濾波器的邊界頻率是實(shí)際濾波器幅度函數(shù)正負(fù)肩峰頻率的中心,但這里不知道正負(fù)肩峰的頻率,因此用通帶和阻帶邊界頻率的中心近似:選擇窗函數(shù)和窗口長(zhǎng)度。阻帶衰減不小于40dB,且階數(shù)最小,因此選擇漢寧窗(其阻帶衰減為44dB,滿(mǎn)足衰減要求)。濾波器的過(guò)渡帶寬為,選擇較小過(guò)渡帶寬,因此窗口長(zhǎng)度為:線(xiàn)性相位延遲常數(shù)為:根據(jù)理想邊界頻率和線(xiàn)性相位延遲常數(shù),求理想單位脈沖響應(yīng):窗函數(shù)與理想單位脈沖響應(yīng)相乘,得到線(xiàn)性相位FIR帶通濾波器的單位脈沖響應(yīng):時(shí)域波形:幅頻特性:程序:N=31;n=0:1:N-1;w1=0.25*pi;w2=0.
11、8*pi;a=(N-1)/2;hd=(0.5-0.5.*cos(pi/15*n).*(sin(w2*(n-a)-sin(w1*(n-a)./(pi*(n-a);hd(a+1)=0.55;win=hanning(N);h=win'.*hd;figure;stem(n,h);xlabel('n');ylabel('h(n)');grid;title('線(xiàn)性相位帶通濾波器FIR的單位脈沖響應(yīng)h(n)');H,w=freqz(h,1);H=20*log10(abs(H);figure;plot(w/pi,H);axis(0 1 -100 0);x
12、label('頻率/Hz');ylabel('幅度/dB');title('線(xiàn)性相位帶通濾波器FIR,漢寧窗,N=31');3.解:求通帶和阻帶的數(shù)字邊界頻率:求理想帶濾波器的邊界頻率。理想帶阻濾波器的邊界頻率是實(shí)際濾波器幅度函數(shù)正負(fù)肩峰頻率的中心,但這里不知道正負(fù)肩峰的頻率,因此用通帶和阻帶邊界頻率的中心近似:選擇窗函數(shù)和窗口長(zhǎng)度。阻帶衰減不小于50dB,且階數(shù)最小,因此選擇海明窗(其阻帶衰減為53dB,滿(mǎn)足衰減要求)。濾波器的過(guò)渡帶寬為,選擇較大過(guò)渡帶寬,由于選擇的是帶阻濾波器,N必須為奇數(shù),因此窗口長(zhǎng)度為:線(xiàn)性相位延遲常數(shù)為:根據(jù)理想邊界頻
13、率和線(xiàn)性相位延遲常數(shù),求理想單位脈沖響應(yīng):窗函數(shù)與理想單位脈沖響應(yīng)相乘,得到線(xiàn)性相位FIR帶通濾波器的單位脈沖響應(yīng):時(shí)域波形:頻域特性:程序:N=67;n=0:1:N-1;w1=0.35*pi;w2=0.65*pi;a=(N-1)/2;hd=(0.54-0.46.*cos(pi/33*n).*(sin(pi*(n-a)+sin(w1*(n-a)-sin(w2*(n-a)./(pi*(n-a);hd(a+1)=0.7;win=hanning(N);h=win'.*hd;figure;stem(n,h);xlabel('n');ylabel('h(n)');
14、grid;title('線(xiàn)性相位帶阻濾波器FIR的單位脈沖響應(yīng)h(n)');H,w=freqz(h,1);H=20*log10(abs(H);figure;plot(w/pi,H);axis(0 1 -90 0);xlabel('頻率/Hz');ylabel('幅度/dB');title('線(xiàn)性相位帶阻濾波器FIR,海明窗,N=33');4.實(shí)驗(yàn)結(jié)果:程序:N=51;M=N-1;n=0:1:M;Wd=0.8;w1=ones(N,1);w2=hanning(N);w3=hamming(N);w4=blackman(N);b1=fir
15、1(M,Wd,'low',w1);b2=fir1(M,Wd,'low',w2);b3=fir1(M,Wd,'low',w3);b4=fir1(M,Wd,'low',w4);figure;stem(n,b1);xlabel('n');ylabel('h(n)');grid;title('FIR低通濾波器的單位脈沖響應(yīng)h(n),矩形窗');H,w=freqz(b1,1);H=20*log10(abs(H);figure;plot(w/pi,H);axis(0 1 -100 0);xlab
16、el('頻率/Hz');ylabel('幅度/dB');title('FIR低通濾波器,矩形窗');figure;stem(n,b2);xlabel('n');ylabel('h(n)');grid;title('FIR低通濾波器的單位脈沖響應(yīng)h(n),漢寧窗');H,w=freqz(b2,1);H=20*log10(abs(H);figure;plot(w/pi,H);axis(0 1 -100 0);xlabel('頻率/Hz');ylabel('幅度/dB');
17、title('FIR低通濾波器,漢寧窗');figure;stem(n,b3);xlabel('n');ylabel('h(n)');grid;title('FIR低通濾波器的單位脈沖響應(yīng)h(n),海明窗');H,w=freqz(b3,1);H=20*log10(abs(H);figure;plot(w/pi,H);axis(0 1 -100 0);xlabel('頻率/Hz');ylabel('幅度/dB');title('FIR低通濾波器,海明窗');figure;stem(n,
18、b4);xlabel('n');ylabel('h(n)');grid;title('FIR低通濾波器的單位脈沖響應(yīng)h(n),布萊克曼窗');H,w=freqz(b4,1);H=20*log10(abs(H);figure;plot(w/pi,H);axis(0 1 -100 0);xlabel('頻率/Hz');ylabel('幅度/dB');title('FIR低通濾波器,布萊克曼窗');5.解:由于(2/40)×5<0.25,因此下邊界在4和5之間,上邊界在6和35之間。又因?yàn)?/p>
19、N為偶數(shù),屬于第二類(lèi)線(xiàn)性相位濾波器,幅度函數(shù)在02上以為中心呈奇對(duì)稱(chēng)分布。因此頻率采樣值的幅度為:頻響采樣的相位為:頻率采樣值H(k)為:?jiǎn)挝幻}沖響應(yīng)h(n)為:實(shí)驗(yàn)結(jié)果:程序:clear;N=40;k=0:N-1;Hk=ones(1,5) 0.5 zeros(1,29) -0.5 -ones(1,4);h1=real(ifft(Hk.*exp(-j*pi*(N-1)*k/N);h2=hanning(N)'.*h1;H1,w=freqz(h1,1);H2,w=freqz(h2,1);figure;plot(w/pi,20*log10(abs(H1),':r',w/pi,
20、20*log10(abs(H2),'-b');axis(0 1 -70 10); grid; xlabel('omega/pi'); ylabel('幅度/dB');legend('矩形窗','漢寧窗'); title('頻率采樣法,N=40,過(guò)渡帶設(shè)一個(gè)采樣點(diǎn)');6.解:(1)由于(2/31)×10<2/3,因此下邊界在9和10之間,上邊界在11和20之間。又因?yàn)镹為奇數(shù),高通,屬于第一類(lèi)線(xiàn)性相位濾波器,幅度函數(shù)在02上以為中心呈偶對(duì)稱(chēng)分布。因此頻率采樣值的幅度為:頻響采樣的相位
21、為:頻率采樣值H(k)為:?jiǎn)挝幻}沖響應(yīng)h(n)為:(2)由于(2/32)×10<2/3,因此下邊界在9和10之間,上邊界在11和20之間。又因?yàn)镹為偶數(shù),高通,屬于第四類(lèi)線(xiàn)性相位濾波器,幅度函數(shù)在02上以為中心呈偶對(duì)稱(chēng)分布。因此頻率采樣值的幅度為:頻響采樣的相位為:頻率采樣值H(k)為:?jiǎn)挝幻}沖響應(yīng)h(n)為:實(shí)驗(yàn)結(jié)果:程序:N=31;n=0:N-1;k=0:N-1;Hk=zeros(1,10) 0.5 ones(1,10) 0.5 zeros(1,9);h1=real(ifft(Hk.*exp(-j*pi*(N-1)*k/N);h2=hamming(N)'.*h1;f
22、igure;stem(n,h1);xlabel('n');ylabel('h(n)');grid;title('線(xiàn)性相位高通濾波器的單位脈沖響應(yīng)h(n),N=31');H1,w=freqz(h1,1);H2,w=freqz(h2,1);figure;plot(w/pi,20*log10(abs(H1),':r',w/pi,20*log10(abs(H2),'-b');axis(0 1 -70 10); grid; xlabel('omega/pi'); ylabel('幅度/dB')
23、;legend('矩形窗','海明窗'); title('頻率采樣法,N=31,過(guò)渡帶設(shè)一個(gè)采樣點(diǎn)');clear;N=32;n=0:N-1;k=0:N-1;Hk=zeros(1,10) 0.5 ones(1,11) 0.5 zeros(1,9);h1=ifft(Hk.*exp(-j*pi*(N-1)*k/N);h2=hamming(N)'.*h1;figure;stem(n,h1);xlabel('n');ylabel('h(n)');grid;title('線(xiàn)性相位高通濾波器的單位脈沖響應(yīng)h(n
24、),N=32');H1,w=freqz(h1,1);H2,w=freqz(h2,1);figure;plot(w/pi,20*log10(abs(H1),':r',w/pi,20*log10(abs(H2),'-b');axis(0 1 -120 10); grid; xlabel('omega/pi'); ylabel('幅度/dB');legend('矩形窗','海明窗'); title('頻率采樣法,N=32,過(guò)渡帶設(shè)一個(gè)采樣點(diǎn)');7.實(shí)驗(yàn)結(jié)果:程序:clear;N=50; n=0:N-1;k=0:N-1;Hk=zeros(1,5) 0.5 ones(1,4) 0.5 zeros(1,4) 0.5 ones(1,4) 0.5 zeros(1,9) 0.5 ones(1,4) 0.5 zeros(1,4) 0.5 ones(1,4) 0.5 zeros(1,4);h1=real(ifft(exp(-j*pi*(N-1)*k/N-j*pi/2).*Hk);figure;stem(n,h1);xlabel('n')
溫馨提示
- 1. 本站所有資源如無(wú)特殊說(shuō)明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶(hù)所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁(yè)內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒(méi)有圖紙預(yù)覽就沒(méi)有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫(kù)網(wǎng)僅提供信息存儲(chǔ)空間,僅對(duì)用戶(hù)上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶(hù)上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶(hù)因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。
最新文檔
- 教育技術(shù)在商業(yè)領(lǐng)域的特殊應(yīng)用探討
- 心理輔導(dǎo)在學(xué)科融合中的重要作用
- 情感智能在教育領(lǐng)導(dǎo)力中的應(yīng)用
- 探索游戲化教育提升學(xué)習(xí)動(dòng)力與效率的途徑
- 學(xué)生心理健康與校園欺凌的預(yù)防策略
- 智能教育工具在醫(yī)療培訓(xùn)中的實(shí)踐與探索
- 教育改革的新方向基于VR的智能教學(xué)模式研究
- 以虛擬帶實(shí)踐深化教育教學(xué)創(chuàng)新改革研究與實(shí)踐案例分析
- 智能化教育工具在提升教學(xué)質(zhì)量中的作用
- 2025撫州幼兒師范高等專(zhuān)科學(xué)校單招《職業(yè)適應(yīng)性測(cè)試》通關(guān)題庫(kù)及參考答案詳解【達(dá)標(biāo)題】
- 爭(zhēng)分奪秒備戰(zhàn)期末主題的學(xué)習(xí)班會(huì)
- 門(mén)店?duì)I銷(xiāo)課件 完整版
- 中國(guó)全部城市名及拼音
- 未成年人紋身治理-主題班會(huì)
- 中國(guó)紅色文化精神學(xué)習(xí)通超星課后章節(jié)答案期末考試題庫(kù)2023年
- 平方差公式公開(kāi)課一等獎(jiǎng)?wù)n件市公開(kāi)課一等獎(jiǎng)?wù)n件省賽課獲獎(jiǎng)?wù)n件
- 2023年公司商業(yè)機(jī)密保密管理
- 2023山西焦煤集團(tuán)有限責(zé)任公司井下操作工招聘2000人筆試模擬試題及答案解析
- 中國(guó)西部國(guó)際汽車(chē)城可行性研究報(bào)告(上)
- 通信工程施工組織方案
- 火災(zāi)報(bào)警系統(tǒng)操作及維護(hù)保養(yǎng)規(guī)程(保得威爾)
評(píng)論
0/150
提交評(píng)論