橢圓帶通濾波器課程設計

1、燕 山 大 學 課 程 設 計 說 明 書燕山大學課 程 設 計 說 明 書題目： 橢圓帶通濾波器的設計 學院（系）：電氣工程學院 年級專業(yè)： 09級檢測一班 學 號： 090103020104 學生姓名： 劉 清 指導教師： 王 娜 教師職稱： 講 師 燕山大學課程設計（論文）任務書課程名稱： 數(shù)字信號處理課程設計 基層教學單位： 指導教師：學 號090103020104學生姓名劉清專業(yè)（班級）09級檢測一班設計題目2橢圓帶通濾波器設計設計技術參數(shù) 采樣頻率為100hz,采樣點數(shù)100，低頻、中頻、高頻信號頻率分別為5hz、15hz 、30hz 設計要求產生一個連續(xù)信號，包含低頻率，中頻，高

2、頻分量，對其進行采樣，進行頻譜分析。設計帶通濾波器對信號進行濾波處理，觀察濾波后信號的頻譜。（熟悉函數(shù)freqz,butter,filter,fft）參考資料數(shù)字信號處理方面資料matlab方面資料周次第一周第二周工作計劃收集消化資料、學習matlab軟件，進行相關參數(shù)計算。編寫仿真程序、調試。指導教師簽字基層教學單位主任簽字說明：此表一式四份，學生、指導教師、基層教學單位、系部各一份。年 月 日 目 錄第1章 摘要5第2章 引言5第3章 設計基本原理6 3.1 matlab簡介.6 3.2 橢圓濾波器的特點7 3.3 采樣定理及相關原理8第4章 濾波器設計過程8 4.1 橢圓濾波器設計結構圖

3、8 4.2 橢圓帶通模擬濾波器設計步驟9 4.3 模擬濾波器的matlab實現(xiàn)9第5章 仿真程序和仿真圖105.1帶通濾波器設計程序及其仿真圖11第6章 仿真分析及總結15心得體會16參考文獻16第一章 摘要本文通過對采樣信號進行頻譜分析和利用設計的橢圓濾波器對采樣信號進行濾波處理，并對仿真結果進行分析和處理。詳細介紹了在基于matlab設計橢圓濾波器過程中常用到的工具和命令。利用matlab濾波濾波器設計函數(shù)直接實現(xiàn)橢圓濾波器的設計，介紹了橢圓型濾波器的基本理論和設計思想，給出了基于matlab設計橢圓型帶通濾波器的具體步驟和利用matlab產生一個包含低頻、中頻、高頻分量的連續(xù)信號，并實現(xiàn)

4、對信號進行采樣和分析。第二章 引言現(xiàn)代生活中，數(shù)字信號經過dac轉換獲得的模擬信號的例子太多了，如聲卡中的語音合成輸出，又如試驗室中的合成信號發(fā)生器等，為了濾除諧波干擾，獲得高精度的模擬信號，大多數(shù)就采用本文介紹的衰減特性陡峭的橢圓低通濾波器。因此橢圓濾波器具有廣泛的應用。橢圓濾波器（elliptic filter）是在通帶和阻帶等波紋的一種濾波器。它比切比雪夫方式更進一步地是同時用通帶和阻帶的起伏為代價來換取過渡帶更為陡峭的特性。橢圓濾波器在通帶和阻帶內都具有等波紋幅頻響應特性。由于其極點位置與經典場論中的橢圓函數(shù)有關，所以由此取名為橢圓濾波器。對于給定的階數(shù)和波紋要求，橢圓濾波器能得到較其

5、它濾波器更窄的過渡帶寬，可以獲得對理想濾波器幅頻響應的最好逼近，是一種性價比很高的濾波器，橢圓濾波器相比其他類型的濾波器，在階數(shù)相同的條件下有著最小的通帶和阻帶波動。它在通帶和阻帶的波動相同，這一點區(qū)別于在通帶和阻帶都平坦的巴特沃斯濾波器，以及通帶平坦、阻帶等波紋或是阻帶平坦、通帶等波紋的切比雪夫濾波器。本文將詳細介紹橢圓帶通濾波器的設計過程，并應用其對輸入信號進行濾波，對濾波前后的頻譜作對比。第三章 設計基本原理3.1 matlab簡介matlab是一套集數(shù)值計算、符號運算及圖形處理等強大功能于一體的科學計算語言。作為強大的科學計算平臺，它幾乎能夠滿足所有的計算需求。matlab的信號處理工

6、具箱提供了設計橢圓濾波器的函數(shù)：ellipord函數(shù)和ellip函數(shù)。通過編程可以很容易由濾波器的技術指標得到所需濾波器的階數(shù)，實現(xiàn)各種類型的橢圓濾波器，大大簡化了橢圓濾波器的設計3.2橢圓濾波器的特點幅值響應在通帶和阻帶內都是等波紋的，對于給定的階數(shù)和給定的波紋要求，橢圓濾波器能獲得較其它濾波器為窄的過渡帶寬，就這點而言，橢圓濾波器是最優(yōu)的，其振幅平方函數(shù)為 （3-1） 其中rn（x）是雅可比(jacobi) 橢圓函數(shù)，為與通帶衰減有關的參數(shù)。橢圓濾波器有如下特點：1.橢圓帶通濾波器是一種零、極點型濾波器，它在有限頻率范圍內存在傳輸零點和極點。2.橢圓帶通濾波器的通帶和阻帶都具有等波紋特性，

7、因此通帶、阻帶逼近特性良好。 3.對于同樣的性能要求，它比前兩種濾波器所需用的階數(shù)都低，而且它的過渡帶比較窄。3.3 采樣定理及相關原理 模擬信號經過 (a/d) 變換轉換為數(shù)字信號的過程稱之為采樣，信號采樣后其頻譜產生了周期延拓，每隔一個采樣頻率fs，重復出現(xiàn)一次。理想采樣就是假設采樣開關閉合時間無限短，即t0的極限情況。此時采樣序列可表示為一個沖激函數(shù)序列。采樣定理：要想采樣后能夠不失真地還原出原模擬信號，則采樣頻率必須大于兩倍原模擬信號頻譜的最高截止頻率（2）。第四章 設計過程4.1橢圓濾波器設計結構圖橢圓濾波器設計結構圖如圖所示：采樣（100hz）連續(xù)混合信號帶通濾波器輸出 圖4.1橢

8、圓濾波器結構框圖4.2設計橢圓模擬帶通濾波器的步驟1.確定模擬濾波器的性能指標：,。2.由性能指標計算出濾波器階次n。3.通過歸一化及去歸一化求出的模擬濾波器ha(s)。4.3模擬濾波器的matlab實現(xiàn)1. matlab的信號處理軟件提供了設計橢圓濾波器的函數(shù)：ellipord函數(shù)和ellip函數(shù)。1).ellipord函數(shù)的功能是求濾波器的最小階數(shù)，其調用格式為： n,wn=ellipord(wp,ws,rp,rs,s) （4-1）其中各參量分別為：n-橢圓濾波器最小階數(shù)；wp-橢圓濾波器通帶截止角頻率；ws-橢圓濾波器阻帶起始角頻率；rp-通帶波紋（db）；rs-阻帶最小衰減(db)；本

9、次設計中，用下面程序可確定濾波器階次：clearrp=0.1;rs=40;wp=10 20;ws=8 22;n,wn=ellipord(wp,ws,rp,rs,s);運行可得：n=6，wn=10 20。因此滿足本次設計要求的橢圓濾波器為6階濾波器。2).ellip函數(shù)的功能是用來設計橢圓濾波器，其調用格式：b,a=ellip(n,rp,rs,wp) （4-2） b,a=ellip(n,rp,rs,wp,ftype) （4-3）其中：ftype = high 為高通濾波器， ftype = low為低通濾波器，ftype = stop為帶阻濾波器。返回長度為n+1的濾波器系數(shù)行向量b和a,進而求

10、得傳遞函數(shù)h(z)： （4-4）2.matlab的信號處理工具箱提供了頻譜分析函數(shù):fft函數(shù)、filter函數(shù)和freqz函數(shù)。1).fft函數(shù)功能：對信號進行離散傅里葉變換。其調用格式：fft(x) （4-5） fft(x,n) （4-6） fft(x,dim)或fft(x,n,dim) （4-7）說明：fft(x)是對輸入信號x的離散傅里葉變換。fft(x,n)是n點傅里葉變換，如果x少于n點則補0湊齊位數(shù)，長于n點則截斷。如果x是個矩陣，列的長度將會以同樣的方式調整，fft會對每列進行傅里葉變換，并返回一個相同維數(shù)的矩陣。fft(x,dim)或fft(x,n,dim)是離散傅里葉變換在

11、dim尺度上的應用。dim可適應于任意維度的fft運算。2).filter函數(shù)功能：利用iir濾波器和fir濾波器對數(shù)據(jù)進行濾波。其調用格式：y=filter(b,a,x) y,zf=filter(b,a,x) y=filter(b,a,x,zi)說明：filter采用數(shù)字濾波器對數(shù)據(jù)進行濾波，其實現(xiàn)采用移位直接型結構，因而適用于iir和fir濾波器。濾波器的系統(tǒng)函數(shù)為 （4-8） 即濾波器系數(shù)a=a0 a1 a2 .an,b=b0 b1 .bm,輸入序列矢量為x。這里，標準形式為a0=1，如果輸入矢量a時，a01，則matlab將自動進行歸一化系數(shù)的操作；如果a0=0，則給出出錯信息。y=f

12、ilter(b,a,x)利用給定系數(shù)矢量a和b對x中的數(shù)據(jù)進行濾波，結果放入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ù)功能：離散時間系統(tǒng)的頻率響應。其調用格式：h,w=freqz(b,a,n) （4-9） h,f=freqz(b,a,n,fs) （4-10） h=freqz(b,a,w) （4-11） h=freqz(b,a,f,fs) （4-12） freqz(b,a,n) （4-13）說明: freqz 用于計算數(shù)字濾波器h(z)的

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

14、由矢量f指定。freqz(b,a,n) 用于在當前圖形窗口中繪制幅頻和相頻特性曲線。 第五章 仿真程序和仿真圖5.1帶通濾波器設計程序及其仿真圖1. 由matlab軟件可實現(xiàn)指定信號的輸入fs=100;t=(1:100)/fs;s=sin(2*pi*t*5)+sin(2*pi*t*15)+sin(2*pi*t*30); figure(1);subplot(111);plot(t,s);xlabel(時間(秒);ylabel(幅值); 其波形如下圖圖5.1連續(xù)信號波形圖2.通過ellip函數(shù)設計出橢圓帶通濾波器b,a=ellip(6,0.1,40,10 20*2/fs);h,w=freqz(b,

15、a,512);figure(2);subplot(111);plot(w*fs/(2*pi),abs(h);xlabel(頻率 (hz);ylabel(頻率響應圖);grid;它的頻率響應圖如下：圖5.2 橢圓帶通濾波器頻率響應圖3.信號通過橢圓帶通濾波器的波形圖sf=filter(b,a,s);figure(3);subplot(111);plot(t,sf);xlabel(時間 (s);ylabel(幅值);axis(0 1 -1 1);s=fft(s,512);sf=fft(sf,512);w=(0:255)/256*(fs/2);如下圖所示圖5.3 信號通過濾波器后的波形圖4.信號通過

16、橢圓帶通濾波器的仿真圖subplot(111);plot(w,abs(s(1:256) sf(1:256);xlabel(頻率(hz);ylabel(傅立葉變換圖);grid;legend(濾波前,濾波后);如下圖所示圖5.4 信號通過濾波器后的頻譜圖注：圖中藍色曲線代表濾波前的幅頻曲線，綠線代表濾波后的幅頻曲線。第六章 仿真分析及總結橢圓濾波器能得到較其它濾波器更窄的過渡帶寬，可以獲得對理想濾波器幅頻響應的最好逼近，是一種性價比很高的濾波器。利用matlab語言，其信號處理工具箱提供了豐富的設計方法，可以使得繁瑣的程序設計簡化成函數(shù)的調用，只要以正確的指標參數(shù)調用函數(shù)，就可以正確快捷地得到設

17、計結果從而較方便地設計出橢圓濾波器。橢圓濾波器可以用較少的階數(shù)獲得很高的選擇特性，在設計過程中可以對比濾波器的特性，隨時更改橢圓濾波器通帶截止角頻率wp，阻帶起始角頻率ws，通帶波紋rp，阻帶最小衰減rs等參數(shù)，觀察濾波器的濾波效果。通過在設計中計算出的階次，考慮到實際應用中計算機計算時間的限制，濾波器階次不可能過高，左右浮動試驗多次后得出6階為最佳階數(shù)。通過圖5-3，5-4可知，設計出的濾波器具有較理想的帶通性能，對高頻和低頻信號能有明顯的截止作用，同時也不會對中頻信號造成損失，很好的滿足了設計要求。心得體會經過一個星期的數(shù)字信號處理課程設計，查閱了多方資料，對數(shù)字電子信號這門課程有更深的了

18、解。同時，對于matlab軟件在圖形處理數(shù)學運算的應用有了很好的了解，學會了matab的基本應用，尤其是對于圖形的掌握。在本次課程設計過程中，我深刻地認識到查閱資料和和同學交流對于自己學習的重要性，如果以后有機會我會繼續(xù)的鍛煉自己的能力好好學習參考文獻1高西全,丁玉美.數(shù)字信號處理(第三版)m.西安:西安電子科技大學出版社,2008.2朱義勝,董輝等 譯.信號處理濾波器設計m.北京:電子工業(yè)出版社,2004.3薛年喜.matlab在數(shù)字信號處理中的應用m.北京:清華大學出版社,2003.4肖有平,胡霞.高階橢圓濾波器的設計與仿真j.電子測量技術,2007,(3).5王靖,李永全.橢圓濾波器matlb設計與實現(xiàn)j