《數字信號處理課程設計巴特沃斯帶阻濾波器的設計》由會員分享，可在線閱讀，更多相關《數字信號處理課程設計巴特沃斯帶阻濾波器的設計（11頁珍藏版）》請在裝配圖網上搜索。

1、鄭州輕工業(yè)學院 課 程 設 計 任 務 書 題目 模擬巴特沃斯帶阻濾波器的設計 專業(yè)、班級 信息08-1 學號 200807070129 姓名 主要內容、基本要求、主要參考資料等： 1、主要內容 1）設計巴特沃斯模擬低通原型濾波器； 2）按頻率變換設計巴特沃斯模擬帶阻濾波器； 3）對比模擬低通原型濾波器和模擬帶阻濾波器并加以分析 2、基本要求 （1）編制MATLAB下的m文件實現主要內容。 （2）書寫課程設計報告。 3、主要參考資料 楊永雙等編.數字信號處理實驗指導書.鄭州：鄭州輕工業(yè)學院，2007 丁玉美等編著.數字信號處理第三版.西安

2、電子科技大學出版社，2006 完 成 期 限： 指導教師簽名： 課程負責人簽名： 年 月 日 9 第一章、理論簡單介紹 1 1、1 MATLAB概述 1 1、2 濾波器設計 2 1、3 基本定理 4 第二章、設計目的、要求、指標 5 2、1 設計目的 5 2、2 設計要求 5 2、3 設計指標 6 第三章、程序代碼和結果分析 6 3、1 程序代碼 6 3、2結果與分析 7 心得體會 8 參考文獻: 10 第一章、理論簡單介

3、紹 1、1 MATLAB概述 MATLAB 是一個可視化的計算程序，被廣泛地應用在科學運算領域里。它具有功能強大、使用簡單等特點，內容包括：數值計算、符號計算、數據擬合、圖形圖像處理、系統模擬和仿真分析等功能。此外，用Matlab還可以進行動畫設計、有限元分析等。 MATLAB系統包括五個主要部分： 1）開發(fā)環(huán)境：這是一組幫助你使用MATLAB的函數和文件的工具和設備。這些工具大部分是圖形用戶界面。它包括MATLAB桌面和命令窗口，命令歷史，和用于查看幫助的瀏覽器，工作空間，文件和查找路徑。 2）MATLAB數學函數庫：這里匯集了大量計算的算法，范圍從初等函數如：求和，正弦，余弦和

4、復數的算術運算，到復雜的高等函數如：矩陣求逆，矩陣特征值，貝塞爾(Bessel)函數和快速傅立葉變換等。 3）MATLAB語言：這是一種高水平的矩陣/數組語言，含有控制流語句，函數，數據結構，輸入/輸出，和面向對象編程特征。它允許“小型編程”以迅速創(chuàng)立快速拋棄型程序，以及“大型編程”以創(chuàng)立完整的大型復雜應用程序。 4）句柄制圖：這是MATLAB制圖系統。它包括高級別的二維、三維數據可視化，圖像處理，動畫，以及表現圖形的命令。它還包括低級別的命令，這使你不但能在MATLAB的應用中建立完整的圖形用戶界面，而且還能完全定制圖形的外觀。 5）MATLAB應用程序界面(API)：這是使你編寫

5、與MATLAB相合的C或Fortran程序的程序庫。它包括從MATLAB中調用程序（動態(tài)鏈接）， 調用MATLAB為計算引擎，和讀寫MAT-文件的設備。 1、2 濾波器設計 在數字信號處理中，濾波占有極其重要的作用，濾波器是譜分析、雷達信號處理、通信信號處理應用中的基本處理算法。目前常用的濾波器設計方法普遍采用Matlab仿真。 1）巴特沃斯濾波器 巴特沃斯濾波器模擬低通濾波器的平方幅頻響應函數為： (1.1) 式中，為低通濾波器的截止頻率，N為濾波器的階數。 巴特沃斯濾波器的特點：通帶內具有最大平坦的頻率特性，且隨著頻率增大平滑單調下降；階數

6、愈高，特性愈接近矩形，過渡帶愈窄，傳遞函數無零點。 這里的特性接近矩形，是指通帶頻率響應段與過渡帶頻率響應段的夾角接近直角。通常該角為鈍角，如果該角為直角，則為理想濾波器。 2）帶阻濾波器的設計規(guī)范 本設計中分別用h(t)、H(s)、Ha(j)表示模擬濾波器的單位脈沖相應、系統函數、頻率響應函數，三者的關系如下： Ha(s) = LT[ha(t)] =edt Ha(j)=FT[ha(t)] =edt 可以用h(t)、H(s)、Ha(j)中任一個描述模擬濾波器，也可以用線性常系數微分方程描述模擬濾波器。但是設計模擬濾波器時，設計指標一般由幅頻相應函數|Ha(j)|給出，而模擬

7、濾波器設計就是根據設計指標，求系統函數H(s)。 工程實際中通常用所謂的損耗函數即衰減函數A（）來描述濾波器的幅頻響應特性，對歸一化幅頻響應函數，A（）定義如下： A（）=-20lg|Ha(j)|=-10lg|Ha(j)| dB 模擬帶阻濾波器的設計指標參數有、、、、。和分別表示帶阻濾波器的通帶下邊界頻率和通帶上邊界頻率；和分別表示帶阻濾波器的阻帶下邊界頻率和阻帶上邊界頻率。1和2分別稱為通帶和阻帶波紋幅度。 帶阻濾波器 1、3 基本原理 1）無失真?zhèn)鬏? 所謂信號無失真?zhèn)鬏斒侵篙斎胄盘柾ㄟ^系統后，輸出信號的幅值和輸入信號的幅值成正比。允許有一定的延時，但沒有波形上的畸變

8、。即系統的幅頻響應||應為常數，相頻響應應與頻率成比例。或者說，濾波器應具有無限寬的定值幅頻與線性相頻。通常定義群延遲為信號系統的延遲時間為td()，用函數表示為： ||=C (常數)且td=C(常數) 2）理想濾波器 理想濾波器應能無失真地傳輸有用信號，而又能完全抑制無用信號。有用信號和無用信號往往占有不同的頻帶。信號能通過濾波器的頻帶稱為通帶,信號被抑制的頻帶稱為阻帶。 由以前所學知識可以知道可能實現的,在具體實現的方面,我們只能想辦法讓實際濾波器的頻率特性只能“逼近”理想濾波器。濾波器的幅頻響應在通帶內不是完全平直的，而是呈波紋變化；在阻帶內，幅頻特性也不為零，而是衰減至某個值

9、；在通帶和阻帶之間存在一個過渡帶，而不是突然下降。 3）傳遞函數 模擬濾波器的設計的理論基礎通常在Laplace域內進行討論，模擬濾波器的技術指標可由平方幅值響應函數的形式給出，而和傳遞函數H(s)存在下面關系： 當給定模擬濾波器的技術指標后，由求出A(-s2)，再適當地選擇分配零極點可求出H(s)。為了使濾波器穩(wěn)定，H(s)的極點必須落在s平面左半平面。濾波器的零點選擇可任取A(-s2)的一半零點，這是因為濾波器對Laplace域表示的傳遞函數并無特殊要求，但如果要求H(s)具有最小相位，零點也必須選擇在s左半平面。 第二章、設計目的、要求、指標

10、 2、1 設計目的 本次數字電路課程設計的主要目的有以下幾點： 1.熟練掌握MATLAB 軟件系統； 2.使學生的數字信號知識系統化、完整化； 3.掌握在MATLAB環(huán)境下進行信號產生、信號變換、濾波器設計的方法； 4. 初步了解和掌握工程實踐中的設計思想和設計方法，具有開發(fā)信號分析、仿真系統的基本能力； 5. 為后繼課程的教學打下扎實基礎； 2、2 設計要求 基于巴特沃斯濾波器設計中心頻率為200Hz，帶寬為150Hz的模擬帶阻濾波器； （1）給定模擬濾波器的性能指標，選擇合適的濾波器類型等； （2）確定濾波器階數； （3）設計模擬低通原型濾波器； （4）按頻

11、率變換設計模擬帶阻濾波器； 2、3 設計指標 模擬低通濾波器的設計指標有Wp, Wp, αp和αs。 Wp：通帶截止頻率； Ws：阻帶截止頻率；αp：通帶中最大衰減系數；αs：阻帶最小衰減系數。αp和αs一般用dB數表示。對于單調下降的幅度特性，可表示成： ( 如果Ω=0處幅度已歸一化到1，即|Ha(j0)|=1,αp和αs表示為： 第三章、程序代碼和結果分析 3、1 程序代碼 %設計中心頻率為200Hz，帶寬為150Hz的模擬帶阻濾波器 for n=1:6; %設計階數為1~6的Butterworth模擬

12、原型濾波器 [z,p,k]=buttap(n); %[z0，p0，k0]=buttap(N) 用于計算N階巴特沃斯歸一化(3dB截止頻率Ωc=1)模擬低通原型濾波器系統函數的零、極點和增益因子 [b,a]=zp2tf(z,p,k); %將零點極點增益形式轉化傳遞函數形式 [H,w]=freqs(b,a,0:0.01:2); %計算指定頻率點的復數頻率響應 subplot(1,1,1),plot(w,abs(H)); %繪制幅度響應 xlabel(w/wc);ylabel(Butterworth ︳H(jw) ︳); title(Butterworth 低通原型濾波器(wc=1))

13、 text(0.4,0.85,n=1) %作必要的標記 text(0.7,0.8,n=2) text(0.89,0.7,n=4) text(1.18,0.1,n=6) hold on; figure(1) grid on; end figure(2) hold on; for n=1:6; %設計階數為1~6的Butterworth模擬帶阻濾波器 [z,p,k]=buttap(n); %[z0，p0，k0]=buttap(N) 用于計算N階巴特沃斯歸一化(3dB截止頻率Ωc=1)模擬低通原型濾波器系統函數的零、極點和增益因子 [b,a]=zp2tf(z,

14、p,k); %將零點極點增益形式轉化傳遞函數形式 f0=200;bw=150; %中心點頻率和阻帶帶寬 [bt,at]=lp2bs(b,a,f0,bw); %頻率轉換（低通變帶阻） [Ht,ft]=freqs(bt,at,0:2:400); %計算帶阻濾波器指定頻率點的復數頻率響應 plot(ft,abs(Ht));%繪制幅度響應 xlabel(f/hz);ylabel(Butterworth ︳H(jw) ︳); title(Butterworth 帶阻濾波器(f0=200,bw=150)); text(68,0.9,n=1) %作必要的標記 text(115,0.8,n

15、=2) text(305,0.8,n=4) text(306,0.9,n=6) grid on; end 3、2 結果與分析 1） 仿真結果 圖 5.1 Butterworth型模擬原型低通濾波器 圖 5.2 Butterworth型模擬帶阻濾波器 2） 結果分析 從圖5.2可以看到，中心頻率為200hz，帶寬為150hz，滿足設計要求。 進一步分析我們可以得到：階數越高，邊沿越陡峭，越接近于理想情況（近似方波邊沿），但是這就意味著設備的復雜程度加大，成本提高。從圖5.1和圖5.2

16、可知，當n=4和n=6就已經很接近了，兩者的收斂性相差無幾，一階系統最為簡單，但效果極差，遠遠達不到設計效果。 心得體會 經過此次課程設計我學到了很多知識，尤其是對matlab軟件的功能有了進一步的了解?！稊底中盘柼幚怼氛n程是一門理論和技術發(fā)展十分迅速、應用非常廣泛的前沿性學科，它的理論性和實踐性都很強，它的特點是：要求的數學知識多，包括高等代數、數值分析、概率統計、隨機過程等。要求掌握的基礎知識強，網絡理論、信號與系統是本課程的理論基礎。與其他學科密切相關，即與通信理論、計算機、微電子技術不可分，又是人工智能、模式識別、神經網絡等新興學科的理論基礎之一。 在課程設計中，matlab

17、中的一些函數不是很清楚，可以通過命令窗口直接點擊“help”，你就可以得到本地機器上matlab的基本的幫助信息。 對于某些不是很明確的命令，只知道大體所屬范圍，譬如說某個工具箱，直接在 命令窗口中點擊 help toolboxname，一般可以得到本工具箱有關的信息：版本號，函數名等。 知道函數名，直接用help funname就可以得到相應的幫助信息。同時我認為我們的工作是一個團隊的工作，團隊需要個人，個人也離不開團隊，必須發(fā)揚團結協作的精神。這一點，是非常寶貴的、重要的。對我而言，知識上的收獲重要，精神上的豐收更加可喜。 參考文獻: 丁玉美，高西全編著．數字信號處理．西安：西安電子科技大學出版社，2006． 鄭君里等編．信號與系統．北京：高等教育出版社，1981． 程佩青編著．數字信號處理教程．北京：清華大學出版社，2007． 鄭阿奇等編．MATLAB實用教程．北京：電子工業(yè)出版社，2004． 楊永雙等編.數字信號處理實驗指導書.鄭州：鄭州輕工業(yè)學院，2007．