基于頻率抽樣法二維非遞歸濾波器設(shè)計及其簡單應(yīng)(編輯修改稿)
2025-07-24 21:17 本頁面
【文章內(nèi)容簡介】 1）點，另外（N/21）點（從N/2+1到N1）則需根據(jù)線性相位的奇偶對稱條件求得。對N為奇數(shù)，設(shè)計公式為 (35) (36) (37)對N為偶數(shù)，設(shè)計公式為 (38) (39) (310) (311) 最小二乘優(yōu)化設(shè)計單位脈沖響應(yīng)為{h(n1,n2),n1=0,1,…,N11,n2=0,1,…,N21}的二維FIR 濾波器,其頻率響應(yīng)為 (312)二維FIR 濾波器的約束最小二乘設(shè)計是指選擇N1N2個脈沖響應(yīng)系數(shù)h(n1,n2) ,使H(ejω1,ejω2) 逼近某個期望響應(yīng), 在逼近誤差小于預(yù)定值的條件下, 使平方逼近誤差的積分最小.一般情形下, N1N2個脈沖響應(yīng)是獨立的。當滿足一定對稱性時, 獨立脈沖響應(yīng)系數(shù)有所減少, 且濾波器具有大多數(shù)應(yīng)用場合所要求的線性相位。 為描述簡單, 考慮N1 , N2 均為奇數(shù)、脈沖響應(yīng)矩形對稱的情形,即h(n1,n2)=h(N11n1,n2)=h(n1,N21n2)=h(N11n1,N21n2) ………………………………………………………………(313)此時,共有數(shù)目為n=(N1+1)(N2+1)/4的獨立脈沖響應(yīng),即h(n1,n2),n1=0,1,…,r1,n2=0,1,…，r2,其中r1=(N11)/2,r2=(N21)/2；頻率響應(yīng)變?yōu)? (314)其中,相頻響應(yīng)為(r1ω1+r2ω2),是頻率ω1,ω2的線性函數(shù)(線性相位) ,幅頻響應(yīng)P(ω1,ω2)等于 (315)是獨立脈沖響應(yīng)的線性函數(shù)。用{αi|i=1,2,…，n}來表示獨立的脈沖響應(yīng),則幅頻響應(yīng)可表示為： (316)其中是w1,w2的某個函數(shù)。當αi=h(n1,n2),0≤n1r1,0≤n2r2時 (317)當αi=h(n1,r2),0≤n1r1時 (318)當αi=h(r1,n2),0≤n2r2時當αi=h(r1,r2)時, =1。令=[，…，]T，α=[α1,α2,…,αn]T ,則幅頻響應(yīng)可進一步化為 (319)當所設(shè)計的濾波器就是一個線性相位濾波器時,只需考慮幅頻響應(yīng)的逼近問題,問題簡化為:使P(ω1,ω2)按如下約束最小二乘準則逼近期望幅頻響應(yīng)D(ω1,ω2)： (320) (321)其中Ω1及Ω2是頻率區(qū)域{(ω1,ω2)p≤ω1≤p,p≤ω2≤p}的子集,ε(ω1,ω2) 為預(yù)定的逼近誤差上限函數(shù)。容易看到,P(ω1,ω2)關(guān)于ω1,ω2矩形對稱, 因此要求D(ω1,ω2)也具有矩形對稱性,并且只需考慮頻率子區(qū)域Π≡{(ω1,ω2)0≤ω1≤p,0≤ω2≤p}上的逼近問題。此時,逼近問題仍然表示為(1)、(2) ,但其中的Ω1和Ω2是頻率子區(qū)域Π的子集。一般情況下,式(416) 的積分需采用數(shù)值方法計算,式(417)的約束用稠密頻率點集上的約束集代替。為此,我們將頻率子區(qū)域Π離散化為稠密的頻率點集Ω={(ω1i,ω2j) =(ip/L,jp/L)|i,j=0,1,…,L}L是一個充分大的正整數(shù)。將式(416)積分替換為Ω1∩Ω上的求和,將式(417)約束替換為Ω2∩Ω上的約束集。二維線性相位FIR濾波器的約束最小二乘設(shè)計問題可化為:[6] (322) (323) 基于頻率抽樣法二維FIR濾波器的程序?qū)崿F(xiàn)Matlab信號處理工具箱中的firls函數(shù)用來設(shè)計具有最小二乘線性相位FIR濾波器。該濾波器可使指定頻段內(nèi)的理想分段線性函數(shù)與濾波器頻響應(yīng)之間的誤差平方和最小[7]。b=firls(n,f,m)可得到n階FIR濾波器，其幅頻特性與由f和m指定的性能匹配。濾波器系數(shù)矢量b共包含有n+1個系數(shù)，且滿足對稱關(guān)系b(k)=b(n+2k) k=1，…，n+1 對應(yīng)濾波器即為Ⅰ型和Ⅱ型線性相位濾波器。b=firls（n,f,m,w）可利用權(quán)值矢量w對各頻率段進行加權(quán)擬合，w的長度為f和m長度的一半，用以指定各頻率段的權(quán)值。b=firls(n,f,m,`ftype`)和b=firls（f,m,w,`ftype`）可指定濾波器的類型：當ftype=hilbert時，設(shè)計的濾波器為奇對稱的線性相位濾波器（Ⅲ型和Ⅳ型），b的系數(shù)滿足以b(k)=b(n+2k),k=1,…,n+1。這類濾波器包括Hilbert變換器。當ftype=differentiator時，也可設(shè)計出Ⅲ型和Ⅳ型濾波器，但采用了特殊的加權(quán)技術(shù)，對非零幅值的頻段，將誤差乘以(1+f)2，這樣可使低頻段誤差大大小于高頻段誤差。[8]fsamp2函數(shù)是使用頻率抽樣法設(shè)計二維FIR濾波器。h= fsamp2(hd),設(shè)計一個頻率響應(yīng)為hd的二維FIR濾波器，返回濾波器系數(shù)到矩陣h,濾波器h有一個頻率響應(yīng)經(jīng)過點hd。如果hd是在M、N之間，那么H 也在M、N之間。hd是一個包含理想頻率響應(yīng)采樣點的矩陣， 之間的等間隔點上，, 或pi弧度。為了結(jié)果的準確, 使用 FREQSPACE返回頻率點得到hd。h = fsamp2(f1,f2,hd,[M N])產(chǎn)生一個從M到N 的FIR濾波器 ，使濾波器響應(yīng)hd在f1和f2之間。得到的濾波器盡可能接近理想的頻率響應(yīng)。最好的結(jié)果是必須有至少M*N個理想的頻率點。如果您指定的點少于M*N點,fsamp2將出現(xiàn)警告。[3]下面是利用該函數(shù)和二維FIR濾濾器的頻率采樣設(shè)計函數(shù)fsamp2進行設(shè)計：設(shè)計要求：利用頻率抽樣法設(shè)計一20階的具有最小二乘線性相位的FIR帶阻濾波器。首先必須給出一維的相應(yīng)FIR濾波器。f=[0 1]。m=[1 1 1 0 0 0 0 1 1 1]；產(chǎn)生一個20階的具有最小二乘線性相位的一維FIR帶阻濾波器。B=firls(20,f,m,`Hilbert`)。[h0,w]=f
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1