基于matlab的數(shù)字濾波器設(shè)計濾波仿真

1、基于matlab的數(shù)字濾波器設(shè)計及濾波仿真基于matlab的數(shù)字濾波器設(shè)計及濾波仿真11/11基于matlab的數(shù)字濾波器設(shè)計及濾波仿真山東科技大學電工電子實驗授課中心創(chuàng)新性實驗結(jié)題報告實驗項目名稱基于matlab的數(shù)字濾波器設(shè)計及濾波仿真三個正弦信號用三種方法進行濾波分別編程法、法與FDATool法一、實驗大綱隨著通信行業(yè)和電子計算機技術(shù)的發(fā)展，數(shù)字信號辦理技術(shù)碰到了越來越廣泛關(guān)注，其理論及算法隨著計算機技術(shù)和微電子技術(shù)的發(fā)展獲得了飛速地發(fā)展，數(shù)字濾波器是數(shù)字信號辦理中最重要的組成部分之一，本文詳盡介紹了利用MATLAB信號辦理工具箱SignalProcessingToolbox）快速有效的

2、設(shè)計由軟件組成的老例數(shù)字濾波器的設(shè)計方法。利用Maltab模擬一個基本輸入信號（如三個正弦頻率信號的疊加），自行設(shè)定濾波器的技術(shù)指標，用程序設(shè)計法設(shè)計各種濾波器（低通，高通，帶通），使用Simulink進行濾波仿真最后獲得分其他三個正弦頻率信號，并與理論上得濾波收效進行比較和談?wù)?。二、實驗?zāi)康纳钊肜斫釯IR和FIR數(shù)字濾波器的看法、基根源理和方法，學會用FADtool設(shè)計數(shù)字濾波器，掌握各種數(shù)字濾波器的設(shè)計方法和特點，熟悉MATLAB語言和MATLAB及其數(shù)字信號辦理工具箱的使用方法，經(jīng)過設(shè)計和仿真來加深學生對濾波的理解并提高設(shè)計能力，學習創(chuàng)新性實驗報告的撰寫方法。三、實驗場所及儀器、設(shè)備和資

3、料：實驗室：J11-318實驗設(shè)備：PC機一臺，Matlab軟件四、實驗內(nèi)容（1）、實驗原理:正弦信號：x(n)Asin(wn)巴特沃斯模擬濾波器的設(shè)計：N,Wc=buttord(Wp,Ws,ap,as,s)N,wc=buttord(wp，ws，p，s)用于計算巴特沃斯數(shù)字濾波器的階數(shù)N和3dB截止頻率wc。調(diào)用參數(shù)wp，ws分別為數(shù)字濾波器的通帶、阻帶截止頻率的歸一化值，要求：0wp1，0ws1。1表示數(shù)字頻率pi。p，s分別為通帶最大衰減和組帶最小衰減(dB)。當wswp時，為高通濾波器;當wp和ws為二元矢量時，為帶通或帶阻濾波器，這時wc也是二元向量。N，wc作為butter函數(shù)的調(diào)用

4、參數(shù)。壇N,c=buttord(p，s，p，s，s)B,A=BUTTER(N,Wn,ftype)用于計算巴特沃斯模擬濾波器的階數(shù)N和3dB截止頻率c。p，s，c均為實質(zhì)模擬角頻率BS,AS=butter(N,Wc,s)B,A=BUTTER(N,Wn)，設(shè)計一個階數(shù)為n，頻率為Wn的低通濾波器;可以設(shè)計高通，帶阻濾波器，其中ftype參數(shù)的形式可以指定何種濾波器，ftype為high時，設(shè)計一個階數(shù)為n，頻率為Wn的高通濾波器;ftype為stop時，獲得濾波器階數(shù)為2*n，頻率范圍為Wn=W1,W2的帶阻濾波器;雙線性變換設(shè)計IIR濾波器：BZ,AZ=bilinear(BS,AS,1/T)信號

5、濾波Y=filter(B,A,X)輸入X為濾波前序列，Y為濾波結(jié)果序列，B/A供應(yīng)濾波器系數(shù)，B為分子，A為分母整個濾波過程是經(jīng)過下面差分方程實現(xiàn)的：a(1)*y(n)=b(1)*x(n)+b(2)*x(n-1)+.+b(nb+1)*x(n-nb)-a(2)*y(n-1)-.-a(na+1)*y(n-na)Y,Zf=FILTER(B,A,X,Zi)，輸入X為濾波前序列，Y為濾波結(jié)果序列，B/A供應(yīng)濾波器系數(shù)，B為分子，A為分母，并輸入Zi指定X的初始狀態(tài)，Zf為最后狀態(tài)矢量（2）Simulink的使用方法：模塊庫中的模塊可以直接用鼠標進行拖曳（選中模塊，按住鼠標左鍵不放）而放到模型窗口中進行辦

6、理。在模型窗口中，選中模塊，則其4個角會出現(xiàn)黑色標記。此時可以對模塊進行以下的基本操作。搬動：選中模塊，按住鼠標左鍵將其拖曳到所需的地址即可。若要走開線而搬動，可按住shift鍵，再進行拖曳。復(fù)制：選中模塊，爾后按住鼠標右鍵進行拖曳即可復(fù)制同樣的一個功能模塊。刪除：選中模塊，按Delete鍵即可。若要刪除多個模塊，可以同時按住Shift鍵，再用鼠標選中多個模塊，按Delete鍵即可。也可以用鼠標采用某地域，再按Delete鍵就可以把該區(qū)域中的全部模塊和線等全部刪除。轉(zhuǎn)向：為了可以序次連接功能模塊的輸入和輸出端，功能模塊有時需要轉(zhuǎn)向。在菜單Format中選擇FlipBlock旋轉(zhuǎn)180度，選擇R

7、otateBlock順時針旋轉(zhuǎn)90度。也許直接按Ctrl+F鍵執(zhí)行FlipBlock，按Ctrl+R鍵執(zhí)行RotateBlock。改變大小：選中模塊，對模塊出現(xiàn)的4個黑色標記進行拖曳即可。模塊命名：先用鼠標在需要更正的名稱上單擊一下，爾后直接更正即可。名稱在功能模塊上的地址也可以變換180度，可以用Format菜單中的FlipName來實現(xiàn)，也可以直接經(jīng)過鼠標進行拖曳。HideName可以隱蔽模塊名稱。顏色設(shè)定：Format菜單中的ForegroundColor可以改變模塊的遠景顏色，BackgroundColor可以改變模塊的背景顏色；而模型窗口的顏色可以經(jīng)過ScreenColor來改變。參

8、數(shù)設(shè)定：用鼠標雙擊模塊，就可以進入模塊的參數(shù)設(shè)定窗口，從而對模塊進行參數(shù)設(shè)定。參數(shù)設(shè)定窗口包括了該模塊的基本功能幫助，為獲得更詳盡的幫助，可以點擊其上的help按鈕。經(jīng)過對模塊的參數(shù)設(shè)定，就可以獲得需要的功能模塊。屬性設(shè)定：選中模塊，打開Edit菜單的BlockProperties可以對模塊進行屬性設(shè)定。包括Description屬性、Priority優(yōu)先級屬性、Tag屬性、Openfunction屬性、Attributesformatstring屬性。其中Openfunction屬性是一個很適用的屬性，經(jīng)過它指定一個函數(shù)名，則當該模塊被雙擊此后，Simulink就會調(diào)用該函數(shù)執(zhí)行，這種函數(shù)在

9、MATLAB中稱為回調(diào)函數(shù)。模塊的輸入輸出信號：模塊辦理的信號包括標量信號和向量信號；標量信號是一種單一信號，而向量信號為一種復(fù)合信號，是多個信號的會集，它對應(yīng)著系統(tǒng)中幾條連線的合成。缺省情況下，大多數(shù)模塊的輸出都為標量信號，對于輸入信號，模塊都擁有一種“智能”的鑒別功能，能自動進行般配。某些模塊經(jīng)過對參數(shù)的設(shè)定，可以使模塊輸出向量信號。（3）SPTool的使用方法:SPTool圖形化信號辦理工具的使用SPTool工具供應(yīng)四個基本的信號辦理圖形用戶界面（GUI）程序，它們分別是（1）信號閱讀器：用于閱讀可視化的信號圖像（2）濾波器設(shè)計器：可用于設(shè)計和編寫FIR和IIR數(shù)字濾波器，絕大多數(shù)MAT

10、LAB信號辦理工具箱供應(yīng)的命令行函數(shù)都可以在這個可視化的濾波器設(shè)計器中被調(diào)用，用戶調(diào)用Pole/Zero編寫器設(shè)計出吻合自己需要的濾波器。3）濾波器閱讀器：這個工具主要用于解析濾波器的特點。4）頻譜閱讀器：用于頻譜解析，使用工具箱供應(yīng)的頻譜估計函數(shù)去解析某個信號序列的功率譜密度。實驗步驟：1、打開SPTool主界面單擊startToolboxessignalprocessionalSPTool，也許直接在命令行輸入SPTool，啟動SPTool工具，出現(xiàn)主界面：主界面有三欄列表，包括信號源列表、濾波器列表和頻譜列表。2、用戶可以從MATLAB主工作空間中導(dǎo)入信號序列、濾波器或頻譜。如從工作空間

11、中導(dǎo)入信號源數(shù)據(jù)：單擊Fileimport，出現(xiàn)設(shè)置參數(shù)，產(chǎn)生的信號源命名為signal1，點擊OK，可看到在主界面的信號列表中增加一個signal1vector項。3、濾波器的設(shè)計（1）可以從工作空間中直接導(dǎo)入已設(shè)計好的濾波器系數(shù)（如大作業(yè)中用FDATool設(shè)計好的濾波器的系數(shù)）：（2）也可以直接在SPTool中進行設(shè)計產(chǎn)生。如選擇頻譜列表中的filt1design，單擊Edit，或者在頻譜列表中單擊New，出現(xiàn)：為比較兩個信號的頻譜，第一創(chuàng)辦一個功率譜密度對象spect1，在SPTool主窗口中選擇signal1vector，爾后單擊頻譜列表下面的Create按鈕，這時會彈出頻譜閱讀器窗口

12、。即使用缺省的參數(shù)創(chuàng)辦了信號signal1的頻譜，單擊Apply，這是頻譜圖像就會顯示出來。同樣可創(chuàng)辦另一個OUT2頻譜圖soect2。爾后在頻譜列表中使用shift同時選中兩個頻譜，單擊View即可顯示在頻譜閱讀器中，這樣就可以比較濾波前后的頻譜變化。2、實驗內(nèi)容（1）設(shè)計三個頻率混雜的信號y，頻率可設(shè)為100Hz、1000Hz、10000Hz；（2）設(shè)計3個巴特沃斯模擬濾波器，分別為低通、高通、帶通；（3）用雙線性變換法將設(shè)計好的巴特沃斯濾波器變換為數(shù)字濾波器；（4）對混雜信號y進行濾波；（5）比較濾波前后信號各重量的變化；（6）調(diào)試濾波器參數(shù)；（8）Simulink仿真；（9）SPToo

13、l仿真；（10）總結(jié)報告。3、實驗步驟(1)Matlab程序設(shè)計：clc;fs=4000;t=0:1/fs:1;x=sin(2*pi*20*t)+sin(2*pi*100*t)+sin(2*pi*200*t);%低通濾波器設(shè)計：wp=2*30/fs;ws=2*60/fs;Rp=1;As=30;subplot(331);plot(t,x);title(原始信號);N,wc=buttord(wp,ws,Rp,As);B,A=butter(N,wc);H,W=freqz(B,A);y=filter(B,A,x);subplot(334);plot(W,abs(H);title(低通濾波器);subp

14、lot(337);plot(t,y)title(30Hz信號);%高通濾波器設(shè)計：wp=2*170/fs;ws=2*120/fs;Rp=20;As=50;subplot(332);plot(t,x);title(原始信號);N,wc=buttord(wp,ws,Rp,As);B,A=butter(N,wc,high);H,W=freqz(B,A);y=filter(B,A,x);subplot(335);plot(W,abs(H);title(高通濾波器);subplot(338);plot(t,y);title(200Hz信號);%帶通濾波器設(shè)計：fp=70130;fs=30170;wp=2

15、*fp/Fs;ws=2*fs/Fs;rp=;rs=50;subplot(333);plot(t,x);title(原始信號);N,wc=ellipord(wp,ws,rp,rs);B,A=ellip(N,rp,rs,wc);H,w=freqz(B,A);y=filter(B,A,x);subplot(336);plot(W,abs(H);title(帶通濾波器);subplot(339);plot(t,y);title(100Hz信號);(2)Sinmulink仿真：參數(shù)設(shè)計：自上而下分別是頻率為20Hz、200Hz、100Hz，三個濾波器分別為低通濾波器，高通濾波器和帶通濾波器。左邊對原信號

16、機進行觀察，右邊對濾波后的信號進行觀察SPTool和FADtool應(yīng)用：SPTool主界面導(dǎo)入原始信號數(shù)據(jù)原始信號：設(shè)計低通濾波器：設(shè)計高通濾波器：設(shè)計帶通濾波器：五、實驗結(jié)果與解析1、實驗現(xiàn)象、數(shù)據(jù)記錄(1)Matlab程序設(shè)計：(2)Sinmulink仿真：上圖左邊為經(jīng)過低頻濾波器濾波后的低頻信號，右邊為原低頻信號上圖左邊為經(jīng)過高頻濾波器濾波后的低頻信號，右邊為原高頻信號上圖左邊為經(jīng)過帶通濾波器濾波后的中頻信號，右邊為原中頻信號(3)SPTool和FADtool應(yīng)用：低通濾波后的信號：高通濾波后的信號：帶通濾波后的信號：指導(dǎo)教師簽字：年月日2、對實驗現(xiàn)象、數(shù)據(jù)及觀察結(jié)果的解析與談?wù)摚航?jīng)過程序設(shè)計的方法，Sinmulink仿真，SPTool和FADtool三種方法均能獲得理想的結(jié)果；