【正文】 hd(n)直接截?cái)嘁鸬?，也稱為 截?cái)嘈?yīng) 。 第 7章 有限脈沖響應(yīng)數(shù)字濾波器的設(shè)計(jì) 第二類線性相位 FIR數(shù)字濾波器的相位 :θ(ω)=π/2ωτ， 類似推導(dǎo)可得到 : 1j j j ( / 2 )g0( e ) ( ) e ( ) eNnnH h n H? ? ????? ? ? ????? 2) 第二類線性相位對(duì) h(n)的約束條件 函數(shù) h(n)cos［ ω(n－ τ)］關(guān)于求和區(qū)間的中心 (N－ 1)/2奇對(duì)稱，是滿足式（ ）的一組解，因?yàn)?cos［ ω(n－ τ)］關(guān)于 n=τ偶對(duì)稱，所以要求 τ和 h(n)滿足如下條件： 10( ) c o s [ ( ) ] 0Nnh n n???????（ ） 第 7章 有限脈沖響應(yīng)數(shù)字濾波器的設(shè)計(jì) 第二類線性相位特性： h(n)應(yīng)當(dāng)關(guān)于 n=(N－ 1)/2點(diǎn)奇對(duì)稱。 1． 線性相位 FIR 第 7章 有限脈沖響應(yīng)數(shù)字濾波器的設(shè)計(jì) ()dd?? ?? ??稱兩種情況為線性相位 ????? ，??)(起始相位， 00)( ?????? ??線性相位 FIR濾波器是指 θ(ω)是 ω的線性函數(shù)，即 嚴(yán)格地說，此時(shí) θ(ω)不具有線性相位特性，但以上兩種情況都滿足群時(shí)延是一個(gè)常數(shù)， 第二類線性相位 第一類線性相位 第 7章 有限脈沖響應(yīng)數(shù)字濾波器的設(shè)計(jì) 1) 第一類線性相位對(duì) h(n)的約束條件 第一類線性相位 FIR數(shù)字濾波器的相位函數(shù) θ(ω)=ωτ， 由式（ ）和（ ）得到 : 1j j jg0( e ) ( ) e ( ) e NnnH h n H? ? ??????????（ ） 1g0( ) ( c os j sin ) ( ) ( c os j sin ) Nnh n n n H? ? ? ?? ????? ? ?? 2. 線性相位 FIR的時(shí)域約束條件 第 7章 有限脈沖響應(yīng)數(shù)字濾波器的設(shè)計(jì) 由式（ ）得到 : 1010( ) c o sc o ss in( ) s inNnNnh n nh n n?????????????將（ ）式中兩式相除得到： 1g01g0( ) c os ( ) c os ( ) sin ( ) sinNnNnH h n nH h n n? ?? ?? ?? ????????????????（ ） 1100( ) c os sin ( ) sin c osNNnnh n n h n n? ?? ? ?????????第 7章 有限脈沖響應(yīng)數(shù)字濾波器的設(shè)計(jì) 簡(jiǎn)得到 : 函數(shù) h(n)sinω(n－ τ)關(guān)于求和區(qū)間的中心 (N－ 1)/2奇對(duì)稱，是滿足（ ）式的一組解。 第 7章 有限脈沖響應(yīng)數(shù)字濾波器的設(shè)計(jì) ★★ 呈偶對(duì)稱呈奇對(duì)稱對(duì),)(??????? 20gH★ 不適合低通，帶阻濾波器 第 7章 有限脈沖響應(yīng)數(shù)字濾波器的設(shè)計(jì) )(?gH表 線性相位 FIR數(shù)字濾波器的時(shí)域和頻域特性一覽 第 7章 有限脈沖響應(yīng)數(shù)字濾波器的設(shè)計(jì) 第 7章 有限脈沖響應(yīng)數(shù)字濾波器的設(shè)計(jì) 四種線性相位 FIR濾波器 相位 特性 幅度特性 適用 特點(diǎn) 1?)(??0)n1N(h)n(h ????????? 2 1N)(N奇 呈偶對(duì)稱，對(duì) ???? 202N偶 偶對(duì)稱呈奇對(duì)稱對(duì)????20 ,??)n1N(h)n(h ????22 1N)( ?????????)(??03N奇 呈奇對(duì)稱，對(duì) ???? 204N偶 偶呈對(duì)稱對(duì)呈奇對(duì)稱，對(duì)?????? 20低通 高通 帶通 帶阻 低通 帶通 帶通 高通 帶通 )nc o s ()n()()(ng??? －＋ ???? 102M hhHg 0( ) 2 ( ) c o s [ ( ) ]MnH h n n? ? ?????)]n([in)n()(ng??? －shH M???? 102g0( ) 2 ( ) s in [ ( ) ]MnH h n n? ? ?????3. 線性相位 FIR數(shù)字濾波器的零點(diǎn)分布特點(diǎn) （ ） 10( ) ( )NnnH z h n z???? ?11001( 1 ) ( 1 ) 10( ) ( ) ( 1 ) ( ) ( )NNnnnnNN m NmH z h n z h N n zh m z z H z???????? ? ? ? ? ??? ? ? ? ?? ? ? ???? 將 h(n)=177。這時(shí) 有 c?? ? )( ?? ?RgW)(H d ?)0(H/)(H c ??第 7章 有限脈沖響應(yīng)數(shù)字濾波器的設(shè)計(jì) 3） 時(shí)， 的主瓣全部在 的通帶內(nèi)，這時(shí)應(yīng)出現(xiàn)正的肩峰。 第 7章 有限脈沖響應(yīng)數(shù)字濾波器的設(shè)計(jì) 以上兩點(diǎn)就是對(duì) hd(n)用矩形窗截?cái)嗪螅陬l域的反映，稱為吉布斯效應(yīng)。 圖 矩形窗函數(shù)長(zhǎng)度的影響 例如：設(shè)計(jì)一個(gè)截止頻率為 的數(shù)字低通濾波器 4/??cw第 7章 有限脈沖響應(yīng)數(shù)字濾波器的設(shè)計(jì) ★★ 對(duì)于同一個(gè)的窗函數(shù)當(dāng)隨著 N的增加，阻帶衰減越快，但是第一旁瓣值不變， N越大，過渡帶越窄，旁瓣越多，通帶越寬。 Bg=8π/N。 ? s=－ 53 dB 第 7章 有限脈沖響應(yīng)數(shù)字濾波器的設(shè)計(jì) 5． 布萊克曼（ Blackman）窗 其頻譜函數(shù)為 （ ） )(1π4c o o )(Bl nRN nN nn N?????????????????????????????????????????????????????????????????????）（）（）（1π4jR1π4jR1π2jR1π2jRjRjBl e)e()e()e(NNNNWWWWWW??????第 7章 有限脈沖響應(yīng)數(shù)字濾波器的設(shè)計(jì) 其幅度函數(shù)為 這樣其幅度函數(shù)由五部分組成，它們都是移位不同，且幅度也不同的 WRg(ω)函數(shù)，使旁瓣再進(jìn)一步抵消。一般 ?加大，主瓣加寬，旁瓣幅度減小 ，典型數(shù)據(jù)為 4 ? 9。 第 7章 有限脈沖響應(yīng)數(shù)字濾波器的設(shè)計(jì) 表 凱塞窗參數(shù)對(duì)濾波器的性能影 第 7章 有限脈沖響應(yīng)數(shù)字濾波器的設(shè)計(jì) 表 6種窗函數(shù)的基本參數(shù) 第 7章 有限脈沖響應(yīng)數(shù)字濾波器的設(shè)計(jì) 例如：設(shè)計(jì)一個(gè)截止頻率為 的數(shù)字低通濾波器 4/??cw第 7章 有限脈沖響應(yīng)數(shù)字濾波器的設(shè)計(jì) ★ 矩形窗過渡帶最窄，布拉克曼窗過渡帶最寬，但是對(duì)于阻帶最小衰減性矩形窗最差，只有－ 21dB，布萊克曼窗最好，達(dá)到－ 75dB。 先按照阻帶衰減選擇窗函數(shù)類型。 解 (1） 選擇窗函數(shù) w(n)，計(jì)算窗函數(shù)長(zhǎng)度 N。因而將所設(shè)計(jì)的 窗函數(shù)法的 MATLAB 第 7章 有限脈沖響應(yīng)數(shù)字濾波器的設(shè)計(jì) hn=fir1(M, wc) 返回 6 dB截止頻率為 wc的 M階（單位脈沖響應(yīng) h(n)長(zhǎng)度 N=M+1） FIR低通（ wc為標(biāo)量）濾波器系數(shù)向量 hn，默認(rèn)選用哈明窗。主要用于設(shè)計(jì)幅度特性形狀特殊的濾波器（如數(shù)字微分器和多帶濾波器等）。high39。所以 h(n)長(zhǎng)度為 N=M+1=24。 % alph=*(rs21)^+*(rs21)。wus=*pi。另外，對(duì)于高通和帶阻濾波器，這里 N只能取奇數(shù)。 N=15和 N=75兩種情況下的幅度內(nèi)插波形 Hg(ω)如圖，圖中的空心圓和實(shí)心圓點(diǎn)分別表示 N=15和N=75時(shí)的頻域幅度采樣。當(dāng) N確定時(shí)， m越大，過渡帶越寬。wp=pi/3。 ％加一個(gè)過渡采樣 thetak=－ pi*(N－ 1)*(0:N－ 1)/N。N=ceil(m+1)*2*pi/Bt)+1。過渡帶采樣值可以設(shè)置為經(jīng)驗(yàn)值，或用累試 （ 5） 對(duì) H(k)進(jìn)行 N點(diǎn) IDFT，得到第一類線性相位 FIR數(shù)字濾波器的單位脈沖響應(yīng)： 1 01( ) I D F T [ ( ) ] ( ) , 0 , 1 , 2 , , 1N knNkh n H k H k W n NN??? ? ? ?? …（ ） （ 6） 檢驗(yàn)設(shè)計(jì)結(jié)果。頻率采樣法同樣滿足這一規(guī)律。這樣，由于時(shí)域混疊及截?cái)?，?h(n)與 hd(n)有偏差。 % B=wlswlp。 %根據(jù) ()式計(jì)算 kaiser窗所需階數(shù) N 第 7章 有限脈沖響應(yīng)數(shù)字濾波器的設(shè)計(jì) wc=(wp+ws)/2/pi。 理想低通濾波器的通帶截止頻率 ωc=(ωs+ωp)/2=，所以得到 : 式中， w(n)是長(zhǎng)度為 24（ ? = st8 4 0 8 2 2 . 2 8 8 72 . 2 8 5 2 . 2 8 5 0 . 2M B? ? ?? ? ???ds i n [0 . 4 ( ) ] 1( ) ( ) ( ) ( ) , 1 1 . 5( ) 2nNh n h n w n w nn? ??? ? ?? ? ? ???第 7章 有限脈沖響應(yīng)數(shù)字濾波器的設(shè)計(jì) 實(shí)現(xiàn)本例設(shè)計(jì)的 MATLAB程序?yàn)? ％ : 例 用凱塞窗函數(shù)設(shè)計(jì)線性相位低通FIR fp=1500。 ％調(diào)用 fir1計(jì)算高通 FIR數(shù)字濾波器的 h(n) ％略去繪圖部分 第 7章 有限脈沖響應(yīng)數(shù)字濾波器的設(shè)計(jì) 運(yùn)行程序得到 h(n)的 25個(gè)值 : h(n)=［ － － － － － － － － － － － － ］ 高通 FIR數(shù)字濾波器的 h(n)及損耗函數(shù)如圖 。 ws=pi/4。當(dāng) ftype=high時(shí)，設(shè)計(jì)高通 FIR濾波器；當(dāng) ftype=stop，且 wc=［ wcl, wcu］時(shí)，設(shè)計(jì)帶阻 FIR濾波器。 過渡帶： 本例中過渡帶寬度 Bt≤ωp－ ωs=π/4， 漢寧窗的精確過渡帶寬度 Bt=，所以要求 Bt=≤π/4N≥。然后根據(jù)過渡帶寬度估計(jì)窗口長(zhǎng)度 N。 N=20 ?改變 N只能改變過渡帶寬窄，不能改變阻帶的衰減性； ?阻帶的衰減性由窗函數(shù)類型決定，矩形窗最差，布萊克曼窗最好 。 ?=，窗函數(shù)接近布萊克曼窗。 布萊克曼窗的四種波形如圖 . ?????????????????????????????????????????????????1π41π4 1π21π2)()(RgRgRgRgRgB l gNWNWNWNWWW??????第 7章 有限脈沖響應(yīng)數(shù)字濾波器的設(shè)計(jì) 圖 布萊克曼窗的四種波形 參數(shù)為 : ?n=－ 57 dB。 漢寧窗的四種波形如圖 第 7章 有限脈沖響應(yīng)數(shù)字濾波器的設(shè)計(jì)