橢圓帶通濾波器的設(shè)計(jì)(matlab)

1、燕 山 大 學(xué) 課 程 設(shè) 計(jì) 說 明 書燕山大學(xué)課 程 設(shè) 計(jì) 說 明 書題目：橢圓帶通濾波器的設(shè)計(jì) 學(xué)院（系）：電氣工程學(xué)院 年級專業(yè)： 12級 學(xué) 號： 學(xué)生姓名： 指導(dǎo)教師： 教師職稱： 電氣工程學(xué)院課程設(shè)計(jì)任務(wù)書課程名稱： 數(shù)字信號處理課程設(shè)計(jì) 基層教學(xué)單位：儀器科學(xué)與工程系 指導(dǎo)教師： 學(xué)號學(xué)生姓名（專業(yè)）班級設(shè)計(jì)題目2橢圓帶通濾波器的設(shè)計(jì)設(shè)計(jì)技術(shù)參數(shù)采樣頻率100Hz，采樣點(diǎn)數(shù)100，低頻、中頻、高頻信號頻率分別為5Hz、15Hz、30Hz設(shè)計(jì)要求產(chǎn)生一個(gè)連續(xù)信號，包含低頻，中頻，高頻分量，對其進(jìn)行采樣，進(jìn)行頻譜分析。設(shè)計(jì)高通濾波器對信號進(jìn)行濾波處理，觀察濾波后信號的頻譜。（熟悉

2、函數(shù)freqz，ellip，filter，fft）參考資料數(shù)字信號處理方面資料MATLAB方面資料周次前半周后半周應(yīng)完成內(nèi)容收集消化資料、學(xué)習(xí)MATLAB軟件，進(jìn)行相關(guān)參數(shù)計(jì)算編寫仿真程序、調(diào)試指導(dǎo)教師簽字基層教學(xué)單位主任簽字說明：1、此表一式四份，系、指導(dǎo)教師、學(xué)生各一份，報(bào)送院教務(wù)科一份。 2、學(xué)生那份任務(wù)書要求裝訂到課程設(shè)計(jì)報(bào)告前面。 電氣工程學(xué)院 教務(wù)科摘要 所謂數(shù)字濾波器，是指輸入輸出均為數(shù)字信號，通過數(shù)值運(yùn)算處理改變輸入信號所含頻率成分的相對比例，或者濾除某些頻率成分的數(shù)字器件或程序。數(shù)字濾波器處理精度高、穩(wěn)定、體積小、重量輕、靈活、不存在阻抗匹配問題。典型的模擬濾波器有巴特沃斯濾

3、波器、切比雪夫?yàn)V波器和橢圓濾波器和貝塞爾濾波器，其中橢圓濾波器具有較好的性能。 Matlab是一套集數(shù)值計(jì)算、符號運(yùn)算及圖形處理等強(qiáng)大功能于一體的科學(xué)計(jì)算軟件。作為強(qiáng)大的科學(xué)計(jì)算平臺，它幾乎能夠滿足所有的計(jì)算需求。本課結(jié)合MATLAB設(shè)計(jì)模擬橢圓濾波器。目錄第一章 概論1第二章 信號處理原理22.1橢圓濾波器的基本理論22.2采樣定理3第三章 軟件仿真設(shè)計(jì)43.1橢圓濾波器設(shè)計(jì)結(jié)構(gòu)圖43.2設(shè)計(jì)橢圓模擬帶通濾波器的步驟43.3 MATLAB相關(guān)函數(shù)介紹4第四章 程序和仿真結(jié)果分析84.1帶通通濾波器設(shè)計(jì)程序84.2仿真結(jié)果KK及分析9參考文獻(xiàn)13第1章 概論 橢圓濾波器又稱考爾濾波器。是在通帶

4、和阻帶等波紋的一種濾波器。橢圓濾波器相比其他類型的濾波器在階數(shù)相同的條件下有著最小的通帶和阻帶波動(dòng)。它在通帶和阻帶的波動(dòng)相同，這一點(diǎn)區(qū)別于在通帶和阻帶都平坦的巴特沃斯濾波器以及通帶平坦、阻帶等波紋或是阻帶平坦、通帶等波紋的切比雪夫?yàn)V波器。 橢圓濾波器傳輸函數(shù)是一種較復(fù)雜的逼近函數(shù)，利用傳統(tǒng)的設(shè)計(jì)方法進(jìn)行電路網(wǎng)絡(luò)綜合要進(jìn)行繁瑣的計(jì)算，還要根據(jù)計(jì)算結(jié)果進(jìn)行查表，整個(gè)設(shè)計(jì)、調(diào)整都十分困難和繁瑣。有許多方法都是希望能快速簡便地設(shè)計(jì)并實(shí)現(xiàn)橢圓濾波器從而把電子電路設(shè)計(jì)者從煩瑣的模擬濾波器設(shè)計(jì)中解放出來。Matlab的信號處理工具箱提供了設(shè)計(jì)橢圓濾波器的函數(shù)：ellipord函數(shù)和ellip函數(shù)。通過編程可

5、以很容易由濾波器的技術(shù)指標(biāo)得到所需濾波器的階數(shù)實(shí)現(xiàn)各種類型的橢圓濾波器，大大簡化了橢圓濾波器的設(shè)計(jì) 。本文采用的方法是MATLAB設(shè)計(jì)出濾波器的傳輸函數(shù)，然后再用通用的可編程濾波器來實(shí)現(xiàn)。第二章 信號處理原理2.1橢圓濾波器的基本理論模擬濾波器是電子設(shè)備中最重要的部分之一。常用的濾波器有巴特沃斯濾波器、切比雪夫?yàn)V波器以及橢圓濾波器；貝塞爾濾波器巴特沃斯和切比雪夫?yàn)V波器的傳輸函數(shù)都是一個(gè)常數(shù)除以一個(gè)多項(xiàng)式，為全極點(diǎn)網(wǎng)絡(luò)，僅在無限大處阻帶衰減為無限大，而橢圓函數(shù)濾波器在有限頻率上既有零點(diǎn)又有極點(diǎn)。零、極點(diǎn)在通帶內(nèi)產(chǎn)生等紋波，阻帶內(nèi)的有限傳輸零點(diǎn)減少了過渡區(qū)，可獲得極為陡峭的衰減曲線。也就是說對于給

6、定的階數(shù)和波紋要求，橢圓濾波器能獲得較其它濾波器更窄的過渡帶寬，就這點(diǎn)而言，橢圓濾波器是最優(yōu)的。橢圓濾波器有如下特點(diǎn)： 1、橢圓高通濾波器是一種零、極點(diǎn)型濾波器，它在有限頻率范圍內(nèi)存在傳輸零點(diǎn)和極點(diǎn)。 2、橢圓高通濾波器的通帶和阻帶都具有等波紋特性，因此通帶，阻帶逼近特性良好。 3、對于同樣的性能要求，它比前兩種濾波器所需用的階數(shù)都低，而且它的過渡帶比較窄。 它的振幅平方函數(shù)為 其中RN（x）是雅可比(Jacobi) 橢圓函數(shù)，為與通帶衰減有關(guān)的參數(shù)。2.2采樣定理 模擬信號經(jīng)過 (A/D) 變換轉(zhuǎn)換為數(shù)字信號的過程稱之為采樣，信號采樣后其頻譜產(chǎn)生了周期延拓，每隔一個(gè)采樣頻率fs，重復(fù)出現(xiàn)一次

7、。理想采樣就是假設(shè)采樣開關(guān)閉合時(shí)間無限短，即t®0的極限情況。此時(shí)采樣序列可表示為一個(gè)沖激函數(shù)序列。采樣定理：要想采樣后能夠不失真地還原出原模擬信號，則采樣頻率必須大于兩倍原模擬信號頻譜的最高截止頻率（³2）。第三章 軟件仿真設(shè)計(jì)3.1橢圓濾波器設(shè)計(jì)結(jié)構(gòu)圖橢圓濾波器設(shè)計(jì)結(jié)構(gòu)圖如圖所示：采樣（100HZ）連續(xù)混合信號帶通濾波器輸出圖3.1橢圓濾波器結(jié)構(gòu)框圖3.2設(shè)計(jì)橢圓模擬帶通濾波器的步驟1. 確定模擬濾波器的性能指標(biāo)。2.由性能指標(biāo)計(jì)算出濾波器階次。 3.通過歸一化及去歸一化求出的模擬濾波器。3.3 MATLAB相關(guān)函數(shù)介紹1. Matlab的信號處理軟件提供了設(shè)計(jì)橢圓濾波

8、器的函數(shù)：ellipord函數(shù)和ellip函數(shù)。 (1).Ellipord函數(shù)的功能是求濾波器的最小階數(shù)其調(diào)用格式為 其中：N為橢圓濾波器最小階數(shù)、Wp是橢圓濾波器通帶頻率(Hz)、Ws為橢圓濾波器阻帶頻率(Hz)、Rp為通帶波紋最大衰減(dB)、Rs為阻帶最小衰減(dB)。 在本次設(shè)計(jì)中，Wp分別為10Hz和20Hz、Ws分別為8Hz和22Hz,、Rp和Rs各為1dB和40db。具體程序如下： Wp=10,20; Ws=8,22; Rp=1; Rs=40; N,Wn=ellipord(Wp,Ws,Rp,Rs,'s'); 運(yùn)行可得 Wn =10 20(2).Ellip函數(shù)的功能

9、是用來設(shè)計(jì)橢圓濾波器。其調(diào)用格式 其中為高通濾波器、為低通濾波器、為帶阻濾波器。 返回長度為n+1的濾波器系數(shù)行向量b和a,進(jìn)而求得傳遞函數(shù)H(z)2.頻譜分析函數(shù):fft函數(shù)、filter函數(shù)和freqz函數(shù)。1).fft函數(shù)功能：對信號進(jìn)行離散傅里葉變換。其調(diào)用格式：fft(X) fft(X,N) fft(X,DIM)或fft(X,N,DIM) 說明：fft(X)是對輸入信號X的離散傅里葉變換。fft(X,N)是N點(diǎn)傅里葉變換，如果X少于N點(diǎn)則補(bǔ)0湊齊位數(shù)，長于N點(diǎn)則截?cái)?。如果x是個(gè)矩陣，列的長度將會(huì)以同樣的方式調(diào)整，fft會(huì)對每列進(jìn)行傅里葉變換，并返回一個(gè)相同維數(shù)的矩陣。fft(X,DI

10、M)或fft(X,N,DIM)是離散傅里葉變換在DIM尺度上的應(yīng)用。DIM可適應(yīng)于任意維度的fft運(yùn)算。2).filter函數(shù)功能：利用IIR濾波器和FIR濾波器對數(shù)據(jù)進(jìn)行濾波。其調(diào)用格式： y=filter(b,a,x) y,zf=filter(b,a,x) y=filter(b,a,x,zi)說明：filter采用數(shù)字濾波器對數(shù)據(jù)進(jìn)行濾波，其實(shí)現(xiàn)采用移位直接型結(jié)構(gòu)，因而適用于IIR和FIR濾波器。濾波器的系統(tǒng)函數(shù)為 即濾波器系數(shù)a=a0 a1 a2 .an,b=b0 b1 .bm,輸入序列矢量為x。這里，標(biāo)準(zhǔn)形式為a0=1，如果輸入矢量a時(shí)，a01，則MATLAB將自動(dòng)進(jìn)行歸一化系數(shù)的操作

11、；如果a0=0則給出出錯(cuò)信息。 y=filter(b,a,x)利用給定系數(shù)矢量a和b對x中的數(shù)據(jù)進(jìn)行濾波，結(jié)果放入y矢量中，y的長度取max(N,M)。 y=filter(b,a,x,zi)可在zi中指定x的初始狀態(tài)。y,zf=filter(b,a,x)除得到矢量y外，還得到x的最終狀態(tài)矢量zf。3).freqz函數(shù)功能：離散時(shí)間系統(tǒng)的頻率響應(yīng)。其調(diào)用格式： h,w=freqz(b,a,n) h,f=freqz(b,a,n,Fs) h=freqz(b,a,w) h=freqz(b,a,f,Fs) 說明: freqz 用于計(jì)算數(shù)字濾波器H(Z)的頻率響應(yīng)函數(shù)H(ej)。 h,w=freqz(b,

12、a,n)可得到數(shù)字濾波器的n點(diǎn)幅頻響應(yīng)值，這n個(gè)點(diǎn)均勻地分布在0,上,并將這n個(gè)頻點(diǎn)的頻率記錄在w中，相應(yīng)的頻響值記錄在h中。要求n為大于零的整數(shù),最好為2的整數(shù)次冪,以便采用FFT計(jì)算，提高速度。缺省時(shí)n =512。 h,f=freqz(b,a,n,Fs)用于對H(ej)在0,Fs/2上等間隔采樣n點(diǎn)，采樣點(diǎn)頻率及相應(yīng)頻響值分別記錄在f 和h中。由用戶指定FS（以HZ為單位）值。 h=freqz(b,a,w)用于對H(ej)在0,2上進(jìn)行采樣，采樣頻率點(diǎn)由矢量w指定。 h=freqz(b,a,f,Fs) 用于對H(ej)在0,FS上采樣，采樣頻率點(diǎn)由矢量f指定。 freqz(b,a,n) 用

13、于在當(dāng)前圖形窗口中繪制幅頻和相頻特性曲線。第四章 程序和仿真結(jié)果分析4.1帶通通濾波器設(shè)計(jì)程序Wp=10,20; Ws=8,22;Rp=1;Rs=40;N,Wn=ellipord(Wp,Ws,Rp,Rs,'s');Fs=100;q=(0:255)/256*50;t=(1:100)/Fs;c=sin(2*pi*t*5)+sin(2*pi*t*15)+sin(2*pi*t*30); figure(1);b,a=ellip(N,Rp,Rs,Wn/50); H,w=freqz(b,a,512);plot(w*Fs/(2*pi),abs(H);title('帶通濾波器頻率特性&#

14、39;);xlabel('頻率（HZ）');ylabel('頻率響應(yīng)圖');sf=filter(b,a,c); figure(2);subplot(211);plot(t,c);xlabel('時(shí)間（秒）');ylabel('幅值');title('連續(xù)信號');subplot(212);plot(t,sf);title('通過濾波器后信號波形');xlabel('時(shí)間（秒）');ylabel('幅值');cc=fft(c,512); sff=fft(sf,512);

15、figure(3);subplot(211);plot(q,abs(cc(1:256);title('原連續(xù)信號頻譜特性');xlabel('頻率（HZ）');ylabel('頻率響應(yīng)');subplot(212);plot(q,abs(sff(1:256);title('通過濾波器后信號頻譜特性');xlabel('頻率（HZ）');ylabel('頻率響應(yīng)'); 4.2仿真結(jié)果KK及分析1. 濾波器的頻率特性圖圖4.1帶通濾波器頻率特性 結(jié)果分析：從仿真圖可以看出橢圓濾波器的過渡帶非常窄，且在通

16、帶和阻帶都有波動(dòng)，符合橢圓濾波器的特點(diǎn)。 2. 通過帶通濾波器前后的波形圖圖4.2通過濾波器前后波形圖對比 結(jié)果分析：由圖可以看出，通過橢圓濾波器后的信號成為頻率單一的正弦信號，信號的頻率為15Hz。只通過了中頻信號，因此實(shí)現(xiàn)了帶通濾波器的功能。3. 信號通過帶通濾波器前后的頻譜特性圖4.3通過濾波器前后信號的頻譜圖 結(jié)果分析：通過帶通濾波器前，低、中、高頻信號都有頻率響應(yīng)。而通過濾波器之后，只保留的中頻信號，因此低、高頻信號的頻率響應(yīng)為零。除此之外，通過濾波器后信號的響應(yīng)幅值有所減小。 第五章 心得體會(huì)數(shù)字信號處理是一門重要的專業(yè)基礎(chǔ)課程。通過這周的課程設(shè)計(jì)我加深了對數(shù)字信號處理這門學(xué)科所學(xué)理論知識的理解。為了完成這次橢圓帶通濾波器的設(shè)計(jì)，我花了許多時(shí)間去學(xué)習(xí)matlab這個(gè)軟件，也進(jìn)一步理解了matlab這一軟件的強(qiáng)大之處。相信這會(huì)對以后的工作和學(xué)習(xí)有很大的幫助。這次課程設(shè)計(jì)讓我把理論知識與實(shí)際應(yīng)用結(jié)合了起來，提高自己的實(shí)際動(dòng)手能力和獨(dú)立思考的能力。通過不斷查閱大量的設(shè)計(jì)資料讓自己的設(shè)計(jì)更加完善，這加強(qiáng)了我對信息的獲取能力?？傊@次課程設(shè)計(jì)是自己的能力得到了一定的鍛煉。