基于Matlab的濾波器設(shè)計_第1頁
基于Matlab的濾波器設(shè)計_第2頁
基于Matlab的濾波器設(shè)計_第3頁
全文預覽已結(jié)束

下載本文檔

版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進行舉報或認領(lǐng)

文檔簡介

1、基于Matlab的濾波器設(shè)計 一、巴特沃斯IIR濾波器的設(shè)計在MATLAB下,設(shè)計巴特沃斯IIR濾波器可使用butter函數(shù)。Butter函數(shù)可設(shè)計低通、高通、帶通和帶阻的數(shù)字和模擬IIR濾波器,其特性為使通帶內(nèi)的幅度響應最大限度地平坦,但同時損失截止頻率處的下降斜度。在期望通帶平滑的情況下,可使用butter函數(shù)。butter函數(shù)的用法為:b,a=butter(n,Wn,/ftype/)其中n代表濾波器階數(shù),Wn代表濾波器的截止頻率,這兩個參數(shù)可使用buttord函數(shù)來確定。buttord函數(shù)可在給定濾波器性能的情況下,求出巴特沃斯濾波器的最小階數(shù)n,同時給出對應的截止頻率Wn。buttor

2、d函數(shù)的用法為:n,Wn= buttord(Wp,Ws,Rp,Rs)其中Wp和Ws分別是通帶和阻帶的拐角頻率(截止頻率),其取值范圍為0至1之間。當其值為1時代表采樣頻率的一半。Rp和Rs分別是通帶和阻帶區(qū)的波紋系數(shù)。不同類型(高通、低通、帶通和帶阻)濾波器對應的Wp和Ws值遵循以下規(guī)則:1高通濾波器:Wp和Ws為一元矢量且WpWs;2低通濾波器:Wp和Ws為一元矢量且WpWs;3帶通濾波器:Wp和Ws為二元矢量且WpWs,如Wp=0.1,0.8,Ws=0.2,0.7。二、契比雪夫I型IIR濾波器的設(shè)計在期望通帶下降斜率大的場合,應使用橢圓濾波器或契比雪夫濾波器。在MATLAB下可使用cheb

3、y1函數(shù)設(shè)計出契比雪夫I型IIR濾波器。cheby1函數(shù)可設(shè)計低通、高通、帶通和帶阻契比雪夫I型濾IIR波器,其通帶內(nèi)為等波紋,阻帶內(nèi)為單調(diào)。契比雪夫I型的下降斜度比II型大,但其代價是通帶內(nèi)波紋較大。cheby1函數(shù)的用法為:b,a=cheby1(n,Rp,Wn,/ftype/)在使用cheby1函數(shù)設(shè)計IIR濾波器之前,可使用cheblord函數(shù)求出濾波器階數(shù)n和截止頻率Wn。cheblord函數(shù)可在給定濾波器性能的情況下,選擇契比雪夫I型濾波器的最小階和截止頻率Wn。cheblord函數(shù)的用法為:n,Wn=cheblord(Wp,Ws,Rp,Rs)其中Wp和Ws分別是通帶和阻帶的拐角頻率

4、(截止頻率),其取值范圍為0至1之間。當其值為1時代表采樣頻率的一半。Rp和Rs分別是通帶和阻帶區(qū)的波紋系數(shù)。%*%mode: 1-巴特沃斯低通;2-巴特沃斯高通;3-巴特沃斯帶通;4-巴特沃斯帶阻% 5-契比雪夫低通;6-契比雪夫高通;7-契比雪夫帶通;8-契比雪夫帶阻%fp1,fp2: 通帶截止頻率,當高通或低通時只有fp1有效%fs1, fs2: 阻帶截止頻率,當高通或低通時只有fs1有效%rp: 通帶波紋系數(shù)%as: 阻帶衰減系數(shù)%sample: 采樣率%h: 返回設(shè)計好的濾波器系數(shù)%*%functionb,a=iirfilt(mode,fp1,fp2,fs1,fs2,rp,as,sa

5、mple)wp1=2*fp1/sample;wp2=2*fp2/sample;ws1=2*fs1/sample;ws2=2*fs2/sample;%得到巴特沃斯濾波器的最小階數(shù)N和3bd頻率wnif mode3N,wn=buttord(wp1,ws1,rp,as);elseif mode5N,wn=buttord(wp1 wp2,ws1 ws2,rp,as);%得到契比雪夫濾波器的最小階數(shù)N和3bd頻率wnelseif mode);phase=angle(freq_response);subplot(3,1,2);plot(f,phase);grid; %相頻特性axis(0 sample/2

6、 1.1*min(phase) 1.1*max(phase);ylabel(Phase);xlabel(Frequency-);h=impz(b,a,32); %32點的單位函數(shù)響應t=1:32;subplot(3,1,3);stem(t,h);grid;axis(0 32 1.2*min(h) 1.1*max(h);ylabel(h(n);xlabel(n-);%*%1、ButterWorth模擬和數(shù)字濾波器 (1)butterd函數(shù):ButterWorth濾波器階數(shù)的選擇。 調(diào)用格式:n,Wn=butterd(Wp,Ws,Rp,Rs),在給定濾波器性能的情況下(通帶臨界頻率Wp、阻帶臨界頻

7、率Ws、通帶內(nèi)最大衰減Rp和阻帶內(nèi)最小衰減Rs),計算ButterWorth濾波器的階數(shù)n和截止頻率Wn。 相同參數(shù)條件下的模擬濾波器則調(diào)用格式為:n,Wn=butterd(Wp,Ws,Rp,Rs,s) (2)butter函數(shù):ButterWorth濾波器設(shè)計。 調(diào)用格式:b,a=butter(n,Wn),根據(jù)階數(shù)n和截止頻率Wn計算ButterWorth濾波器分子分母系數(shù)(b為分子系數(shù)的矢量形式,a為分母系數(shù)的矢量形式)。 相同參數(shù)條件下的模擬濾波器則調(diào)用格式為:b,a=butter(n,Wn,s) 例 采樣頻率為1Hz,通帶臨界頻率fp =0.2Hz,通帶內(nèi)衰減小于1dB(p=1);阻帶臨

溫馨提示

  • 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預覽,若沒有圖紙預覽就沒有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負責。
  • 6. 下載文件中如有侵權(quán)或不適當內(nèi)容,請與我們聯(lián)系,我們立即糾正。
  • 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論