數(shù)字信號處理（第2版）課件 錢玲第6、7章-FIR數(shù)字濾波器設(shè)計(jì)、數(shù)字信號處理系統(tǒng)的實(shí)現(xiàn)

第6章FIR數(shù)字濾波器設(shè)計(jì)16.1線性相位FIR數(shù)字濾波器的特性6.2窗函數(shù)設(shè)計(jì)法6.3

頻率采樣設(shè)計(jì)法6.4FIR數(shù)字濾波器的等波紋設(shè)計(jì)6.5FIR濾波器與IIR濾波器的比較6.5用Matlab函數(shù)設(shè)計(jì)FIR數(shù)字濾波器第6章FIR數(shù)字濾波器設(shè)計(jì)2優(yōu)點(diǎn)：①很容易獲得嚴(yán)格的線性相位，避免被處理的信號產(chǎn)生相位失真，這一特點(diǎn)在寬頻帶信號處理、陣列信號理、數(shù)據(jù)傳輸?shù)认到y(tǒng)中非常重要；②極點(diǎn)全部在原點(diǎn)（永遠(yuǎn)穩(wěn)定），無穩(wěn)定性問題；③任何一個(gè)非因果的有限長序列，總可以通過一定的延時(shí)，轉(zhuǎn)變?yōu)橐蚬蛄校砸蚬钥偸菨M足；④無反饋運(yùn)算，運(yùn)算誤差小。缺點(diǎn)：①因?yàn)闊o極點(diǎn)，要獲得好的過渡帶特性，需以較高的階數(shù)為代價(jià)；②無法利用模擬濾波器的設(shè)計(jì)結(jié)果，一般無解析設(shè)計(jì)公式，要借助計(jì)算機(jī)輔助設(shè)計(jì)程序完成。6.1線性相位FIR數(shù)字濾波器的特性36.1.1線性相位特性6.1.2線性相位FIR濾波器的幅度函數(shù)6.1.3線性相位FIR濾波器的零點(diǎn)特性6.1線性相位FIR數(shù)字濾波器的特性4

因?yàn)樗且环N線性時(shí)不變系統(tǒng)，可用卷積和形式表示FIR數(shù)字濾波器的差分方程描述對應(yīng)的系統(tǒng)函數(shù)6.1.1線性相位特性51.線性相位的條件線性相位意味著一個(gè)系統(tǒng)的相頻特性是頻率的線性函數(shù)，即式中

為常數(shù)通過這一系統(tǒng)的各頻率分量的時(shí)延為一相同的常數(shù)，系統(tǒng)的群時(shí)延為FIR濾波器的頻響特性為6.1.1線性相位特性6式中，H(

)是正或負(fù)的實(shí)函數(shù)，稱為幅度函數(shù)（amplitudefunction）；

(

)是濾波器的相位函數(shù)(6.1-3)將

(

)=

代入（6.1-3）式得(6.1-4)(6.1-5)式（6.1-5）中，等式兩邊實(shí)部與虛部應(yīng)當(dāng)各自相等，同樣實(shí)部與虛部的比值應(yīng)當(dāng)相等，從而得到6.1.1線性相位特性7將上式兩邊交叉相乘，應(yīng)用三角函數(shù)的恒等關(guān)系(6.1-6)滿足上式的條件是(6.1-7)此時(shí)濾波器的單位脈沖響應(yīng)將偶對稱于中點(diǎn)

。

6.1.1線性相位特性82.線性相位+附加的相位(6.1-8)式中，β為一常數(shù)，上式的群時(shí)延依然是常數(shù)

。(6.1-9)可以得到類似的解為(6.1-10)此時(shí)濾波器的單位脈沖響應(yīng)將奇對稱于中點(diǎn)

。

6.1.1線性相位特性9

偶對稱線性相位特性

奇對稱6.1.2線性相位FIR濾波器的幅度函數(shù)10FIR濾波器為了得到線性相位特性，需要濾波器的單位脈沖響應(yīng)滿足偶對稱或奇對稱，同時(shí)濾波器的長度N可以為偶數(shù)或奇數(shù)，不同的組合可以得到4種類型的FIR線性相位濾波器。1.類型1——h[n]偶對稱，N為奇數(shù)其單位脈沖響應(yīng)滿足：6.1.2線性相位FIR濾波器的幅度函數(shù)相位函數(shù)為幅度函數(shù)為因?yàn)閷τ诔逝紝ΨQ，則對于是偶對稱。故：116.1.2線性相位FIR濾波器的幅度函數(shù)122.類型2——h[n]偶對稱，N為偶數(shù)其單位脈沖響應(yīng)滿足：相位函數(shù)為6.1.2線性相位FIR濾波器的幅度函數(shù)幅度函數(shù)為因?yàn)閷τ诔势鎸ΨQ，則對于是奇對稱。而且當(dāng)時(shí)，，則H(z)在z=

1處必然有一個(gè)零點(diǎn)，136.1.2線性相位FIR濾波器的幅度函數(shù)143.類型3——h[n]奇對稱，N為奇數(shù)其單位脈沖響應(yīng)滿足：相位函數(shù)為，中間項(xiàng)必然為零6.1.2線性相位FIR濾波器的幅度函數(shù)幅度函數(shù)為對于是奇對稱。而且當(dāng)時(shí)，，則H(z)在z=

1

處必然有零點(diǎn)，因?yàn)閷τ诔势鎸ΨQ，上式可以表示為156.1.2線性相位FIR濾波器的幅度函數(shù)164.類型4——h[n]奇對稱，N為偶數(shù)其單位脈沖響應(yīng)滿足：相位函數(shù)為6.1.2線性相位FIR濾波器的幅度函數(shù)幅度函數(shù)為對于是奇對稱。而且當(dāng)時(shí)，，則H(z)在z=1

處必然有零點(diǎn)，176.1.2線性相位FIR濾波器的幅度函數(shù)類型1:h[n]偶對稱、N為奇數(shù)，四種濾波器都可設(shè)計(jì)；類型2:h[n]偶對稱、N為偶數(shù)，可設(shè)計(jì)低、帶通濾波器，不

能設(shè)計(jì)高通和帶阻；類型3:h[n]奇對稱、N為奇數(shù)，只能設(shè)計(jì)帶通濾波器，

其它

濾波器都不能設(shè)計(jì)；類型4:h[n]奇對稱、N為偶數(shù)，可設(shè)計(jì)高、帶通濾波器，不

能設(shè)計(jì)低通和帶阻。186.1.2線性相位FIR濾波器的幅度函數(shù)類型1類型2h[n]h[n]h[n]a[n]b[n]196.1.2線性相位FIR濾波器的幅度函數(shù)類型3類型4h[n]h[n]d[n]c[n]206.1.2線性相位FIR濾波器的幅度函數(shù)例6.1-1N=5,h

[0]=h[1]=h

[3]=h

[4]=

1/2,h

[2]=2，求幅度函數(shù)H(ω)。解:

Ｎ為奇數(shù),并且h[n]滿足偶對稱關(guān)系a[0]=h

[2]=2a[1]=2h

[3]=

1a

[2]=2h

[4]=

1H(ω)=2

cosω

cos2ω

=2

(cosω+cos2ω)216.1.3線性相位FIR濾波器的零點(diǎn)特性將代入式得到

由該式可看出，若z=zi是H(z)的零點(diǎn)，則z=z-1i也一定是H(z)的零點(diǎn)。由于h[n]是實(shí)數(shù)，H(z)的零點(diǎn)還必須共軛成對，所以z=z*i

及z=1/z*也必是零點(diǎn)。226.1.3線性相位FIR濾波器的零點(diǎn)特性所以線性相位濾波器的零點(diǎn)必須是互為倒數(shù)的共軛對，即成四出現(xiàn)，這種共軛對共有四種①既不在單位園上，也不在實(shí)軸上，有四個(gè)互為倒數(shù)的兩組共軛對，zi,z*I,1/zi,1/z*i②在單位圓上，但不在實(shí)軸上，因倒數(shù)就是自己的共軛，所以有一對共軛零點(diǎn)，zi,z*i圖（b）236.1.3線性相位FIR濾波器的零點(diǎn)特性③不在單位圓上，但在實(shí)軸上，是實(shí)數(shù),共軛就是自己，所以有一對互為倒數(shù)的零點(diǎn),zi,1/zi圖（c）④又在單位圓上，又在實(shí)軸上，共軛和倒數(shù)都合為一點(diǎn)，所以成單出現(xiàn)，只有兩種可能，zi=1或zi=-1圖(d)246.1.3線性相位FIR濾波器的零點(diǎn)特性我們從幅度響應(yīng)的討論中已經(jīng)知道，對于類型2濾波器（h[n]偶對稱，N為偶數(shù)），。

即

是的零點(diǎn)，既在單位圓，又在實(shí)軸，所以，必有單根；類型3濾波器，h[n]奇對稱，N為奇數(shù)，因所以z=1，z=-1都是H(z)的單根；類型4濾波器，

h[n]奇對稱，N為偶數(shù)，H(0)=0，所以z=1是H(z)的單根。所以，h[n]奇對稱→H(0)=0N為偶數(shù)→H(π)=0線性相位濾波器是FIR濾波器中最重要的一種，應(yīng)用最廣。使用時(shí)應(yīng)根據(jù)需用選擇其合適類型，并在設(shè)計(jì)時(shí)遵循其約束條件。256.2窗函數(shù)設(shè)計(jì)法（時(shí)間窗口法）266.2.1窗函數(shù)設(shè)計(jì)法的基本原理6.2.2幾種常用的窗函數(shù)6.2窗函數(shù)設(shè)計(jì)法（時(shí)間窗口法）27

如果希望得到的濾波器的理想頻率響應(yīng)為，那么FIR濾波器的設(shè)計(jì)就在于尋找一個(gè)實(shí)際頻響特性去逼近，逼近方法有三種：窗口設(shè)計(jì)法（時(shí)域逼近）頻率采樣法（頻域逼近）最優(yōu)化設(shè)計(jì)（等波紋逼近）窗函數(shù)設(shè)計(jì)法是從單位脈沖響應(yīng)序列著手，使實(shí)際濾波器的h[n]逼近理想的單位脈沖響應(yīng)序列hd[n]。6.2.1窗函數(shù)設(shè)計(jì)法的基本原理

28

一般來說，理想頻響是分段恒定，在邊界頻率處有突變點(diǎn)，所以，這樣得到的理想單位脈沖響應(yīng)hd[n]往往都是無限長序列，而且是非因果的。但FIR的h

[n]是有限長的，問題是怎樣用一個(gè)有限長的序列去近似無限長的hd[n]

。最簡單的辦法是直接截取一段hd[n]代替h

[n]。這種截取可以形象地想象為h

[n]是通過一個(gè)“窗口”所看到的一段hd[n]，因此，h

[n]也可表達(dá)為h

[n]和一個(gè)“窗函數(shù)”的乘積，即hd[n]可以從理想頻響通過IDTFT變換獲得h

[n]=

hd[n]

w[n]式中w

[n]稱為窗函數(shù)6.2.1窗函數(shù)設(shè)計(jì)法的基本原理

29窗函數(shù)法設(shè)計(jì)的基本步驟（1）構(gòu)造希望逼近的濾波器頻響特性以低通線性相位FIR濾波器設(shè)計(jì)為例（2）求出hd[n]：對進(jìn)行IDTFT得到（3）加窗得到FIR濾波器的單位脈沖響應(yīng)h

[n]h

[n]=

hd[n]w[n]（6.2-3）6.2.1窗函數(shù)設(shè)計(jì)法的基本原理

302.窗函數(shù)法設(shè)計(jì)的性能分析式（6.2-3）表示，的逼近誤差必然與窗函數(shù)直接相關(guān)。所以逼近誤差實(shí)質(zhì)上就是加窗的影響，窗函數(shù)的類型和長度都會(huì)影響逼近誤差。設(shè)窗函數(shù)為矩形窗函數(shù)，即wR[n]=RN[n]（6.2-4）（6.2-5）矩形窗頻譜的線性相位中幅度函數(shù)6.2.1窗函數(shù)設(shè)計(jì)法的基本原理

31N矩形窗幅度函數(shù)的主瓣寬度為，但旁瓣電平較高，第一旁瓣電平為-13.5dB。（6.2-9）6.2.1窗函數(shù)設(shè)計(jì)法的基本原理

32卷積過程分析如下：① 當(dāng)

=0時(shí)的幅度函數(shù)H(0)，根據(jù)（6.2-9）式，此時(shí)幅度函數(shù)應(yīng)該是Hd(θ)與WR(-θ)兩個(gè)函數(shù)乘積的積分，就是WR(-θ)在-

c到

c一段內(nèi)的積分面積；（6.2-9）②

當(dāng)

=

c時(shí)的幅度函數(shù)H(

c)，此時(shí)幅度函數(shù)應(yīng)該是Hd(θ)與WR(

-θ)的一半重疊，因此H(

c)/H(0)=0.56.2.1窗函數(shù)設(shè)計(jì)法的基本原理

33③

當(dāng)

=

c-2π/N時(shí)，WR(

-θ)的主瓣在Hd(θ)的通帶

c之內(nèi)，而右邊具有負(fù)面積的第一旁瓣已全部移出通帶，因此卷積結(jié)果最大值，即H(

c-2π/N)為最大值，幅度函數(shù)出現(xiàn)正肩峰；（6.2-9）④ 當(dāng)

=

c+2π/N時(shí)，WR(

-θ)的全部主瓣在Hd(θ)的通帶

c之外，而通帶內(nèi)第一旁瓣起著主導(dǎo)作用，使得負(fù)值的面積大于正值的面積，因此H(

c+2π/N)為最小值，幅度函數(shù)出現(xiàn)負(fù)肩峰；6.2.1窗函數(shù)設(shè)計(jì)法的基本原理

34⑤ 當(dāng)

>

c+2π/N時(shí)，頻率

繼續(xù)增加，

WR(

-θ)左邊旁瓣與通帶相乘部分的幅度越來越小，因此,

c+2π/N<

≤π區(qū)間，

H(

)的值就在零值附近波動(dòng)，而且波動(dòng)幅度越來越小。由圖6.2-2（e）可看出，

H(

)是對Hd(

)的逼近，它可能在個(gè)別頻率點(diǎn)上取零值，但不在一個(gè)區(qū)間上恒為零，這也是可實(shí)現(xiàn)性的要求。H(

)相對Hd(

)的逼近誤差是波動(dòng)的，通帶和阻帶的最大波動(dòng)誤差相等，即，相當(dāng)于（6.2-10）這個(gè)最大波動(dòng)誤差是矩形窗固有的，與窗函數(shù)長度N無關(guān)，這是吉布斯現(xiàn)象，N增加時(shí)，主瓣寬度變小，幅度的波動(dòng)起伏振蕩變密。6.2.1窗函數(shù)設(shè)計(jì)法的基本原理

35理想低通濾波器的單位脈沖響應(yīng)截短后，對幅度函數(shù)產(chǎn)生的主要影響：① 理想濾波器的Hd(

)在通帶截止頻率

c處的間斷點(diǎn)變成了連續(xù)曲線，從而使H

(

)出現(xiàn)一個(gè)過渡帶，其寬度取決于窗函數(shù)的主瓣寬度，對于矩形窗WR(

)其主瓣寬度等于4π/N

。② 由于窗函數(shù)旁瓣的作用，使幅度函數(shù)出現(xiàn)波動(dòng)。旁瓣所包圍的面積越大，幅度函數(shù)的波動(dòng)越大。過渡帶：通帶邊界頻率

p和帶阻邊界頻率

s通過數(shù)值計(jì)算，可以得到矩形窗設(shè)計(jì)時(shí)的過渡帶為：（6.2-11）6.2.2窗函數(shù)設(shè)計(jì)法的基本原理

363.窗函數(shù)法設(shè)計(jì)的其他濾波器對于其他類型濾波器，只需修改和的表達(dá)式，選擇窗長度和加窗過程是一致的，表6.2-1標(biāo)準(zhǔn)理想濾波器的設(shè)計(jì)公式（適合矩形窗）6.2.2幾種常用的窗函數(shù)37為了加大阻帶衰減，只能改善窗函數(shù)的形狀。一般希望窗函數(shù)滿足兩項(xiàng)要求：①窗函數(shù)幅度函數(shù)的主瓣盡可能地窄，以獲得較窄的過渡帶；②盡量減少幅度函數(shù)的最大旁瓣的相對幅度，也就是能量盡量集中于主瓣，這樣使肩峰和波紋減小，就可增大阻帶的衰減。但是這兩項(xiàng)要求是不能同時(shí)得到滿足，往往是增加主瓣寬度以換取對旁瓣的抑制。需要設(shè)計(jì)可用的窗函數(shù)，其幅度函數(shù)旁瓣電平要較小，而主瓣就會(huì)加寬。這樣窗函數(shù)在邊沿處（n=0和n=N

1附近）比矩形窗變化要平滑而緩慢，以減小由陡峭的邊緣所引起的旁瓣分量，使阻帶衰減增大，但主瓣寬度會(huì)比矩形窗的要寬，這就造成濾波器幅度函數(shù)過渡帶的加寬。6.2.2幾種常用的窗函數(shù)381.常用的固定窗函數(shù)（1）矩形窗（Rectangular）（2）三角窗（Bartlett）窗主瓣寬度為

6.2.2幾種常用的窗函數(shù)39（3）漢寧（Hann）窗,又稱升余弦窗當(dāng)N>>1時(shí)，幅度函數(shù)近似為這三部分之和，使旁瓣相互抵消，能量更集中在主瓣，但代價(jià)是主瓣寬度比矩形窗的主瓣寬度增加一倍。主瓣寬度為6.2.2幾種常用的窗函數(shù)40（4）哈明（Hamming）窗,又稱改進(jìn)的升余弦窗當(dāng)N>>1時(shí)，幅度函數(shù)近似為它是對漢寧窗的改進(jìn)，在主瓣寬度（對應(yīng)第一零點(diǎn)的寬度）相同的情況下，旁瓣進(jìn)一步減小，可使99.96%的能量集中在主瓣內(nèi)。主瓣寬度為6.2.2幾種常用的窗函數(shù)41（5）布萊克曼（Blackman）窗,又稱二階升余弦窗當(dāng)N>>1時(shí)，幅度函數(shù)近似為主瓣寬度為增加一個(gè)二次諧波余弦分量，可進(jìn)一步降低旁瓣，但主瓣寬度進(jìn)一步增加，是矩形窗主瓣的三倍，增加N可減少過渡帶。幾種窗函數(shù)時(shí)域序列包絡(luò)形狀6.2.2幾種常用的窗函數(shù)42圖6.2-5四種窗函數(shù)的歸一化幅度譜（N=49）6.2.2幾種常用的窗函數(shù)43圖6.2-6四種窗口在同一指標(biāo)下設(shè)計(jì)FIR低通濾波器的幅頻特性（N=49，

c=0.4

）6.2.2幾種常用的窗函數(shù)44表6.2-2幾個(gè)窗函數(shù)的性能參數(shù)窗函數(shù)主瓣寬度過渡帶旁瓣峰值衰減/dB阻帶最小衰減/dB矩形4

/N1.8

/N-13-21三角8

/N4.2

/N-26-25漢寧8

/N6.2

/N-31-44哈明8

/N6.6

/N-41-53布萊克曼12

/N11

/N-57-742.可變窗函數(shù)--凱澤（Kaiser）窗凱澤窗可以在主瓣寬度與旁瓣衰減之間自由選擇。式中I0(x)是零階貝塞爾函數(shù)，形狀參數(shù)β可自由選擇，決定主瓣寬度與旁瓣衰減。β越大，w[n]窗越窄，其頻譜的主瓣變寬，旁瓣變小。一般取4<β<9。6.2.2幾種常用的窗函數(shù)45圖6.2-7凱澤窗（a）凱澤窗函數(shù)（b）凱澤窗頻譜當(dāng)n=(N

1)/2時(shí),凱澤窗處于最大值1。β=5.44接近哈明，β=8.5

接近布萊克曼，β=0為矩形。圖6.2-7(a)給出了β=4.5和β=8.5兩種情況下的凱澤窗函數(shù)時(shí)域序列的包絡(luò)圖，圖6.2-7(b)

給出了β=4.5和β=8.5兩種情況下的凱澤窗函數(shù)的歸一化幅度譜。6.2.2幾種常用的窗函數(shù)46凱澤窗設(shè)計(jì)有經(jīng)驗(yàn)公式可供使用，給出過渡帶△

，阻帶衰減

s(dB),則可求凱澤窗FIR濾波器的階數(shù)N和參數(shù)β，即例6.2-1用凱澤窗設(shè)計(jì)一個(gè)FIR低通濾波器，其性能指標(biāo)為：解：根據(jù)上述濾波器的性能指標(biāo)可知：濾波器通帶與阻帶的波動(dòng)均為:（6.2-17）（6.2-18）6.2.2幾種常用的窗函數(shù)47濾波器的過渡帶為：濾波器的截止頻率為：

設(shè)計(jì)理想低通濾波器的頻響特性為：理想低通濾波器的單位脈沖響應(yīng)為：實(shí)際低通濾波器的單位脈沖響為：6.2.2幾種常用的窗函數(shù)48例（補(bǔ)充）用窗函數(shù)設(shè)計(jì)線性相位高通FIRDF，要求通帶邊界頻率，通帶最大衰減，阻帶截止頻率，阻帶最小衰減。

解：（1）選擇窗函數(shù)因?yàn)樽鑾ё钚∷p，可選擇漢寧窗、哈明窗、布萊克曼窗。這里選擇漢寧窗根據(jù)過渡帶寬高通，N為奇數(shù)，N=256.2.2幾種常用的窗函數(shù)49

（2）期望理想濾波器（3）確定期望濾波器的單位脈沖響應(yīng)6.2.2幾種常用的窗函數(shù)50

全通濾波器低通濾波器（4）加窗6.2.2幾種常用的窗函數(shù)51

6.3頻率采樣設(shè)計(jì)法526.3.1頻率采樣法設(shè)計(jì)原理6.3.2頻率采樣法的設(shè)計(jì)步驟6.3頻率采樣設(shè)計(jì)法53工程上，常給定頻域上的技術(shù)指標(biāo)，所以采用頻域設(shè)計(jì)更直接?；舅枷胧顾O(shè)計(jì)的FIR數(shù)字濾波器的頻率特性在某些離散頻率點(diǎn)上的值準(zhǔn)確地等于所需濾波器在這些頻率點(diǎn)處的值，在其它頻率處的特性則有較好的逼近。內(nèi)插公式頻率采樣法則是從頻域出發(fā)，對理想濾波器的頻響特性

進(jìn)行采樣，然后利用采樣值來實(shí)現(xiàn)FIR數(shù)字濾波器的設(shè)計(jì)。6.3.1頻率采樣法設(shè)計(jì)原理54將理想濾波特性經(jīng)單位圓上的N等分采樣后得到頻域采樣值：單位脈沖響應(yīng)：其系統(tǒng)函數(shù)為：令

（6.3-4）（6.3-1）6.3.1頻率采樣法設(shè)計(jì)原理55

（6.3-4）在z平面單位圓上，即，可以得到頻響特性可以表示為：（6.3-5）其中

(

)是內(nèi)插函數(shù)（6.3-6）由式(6.3-5)可以看出，在各頻率采樣點(diǎn)上，實(shí)際濾波器的頻響特性與理想濾波器頻響特性的數(shù)值嚴(yán)格相等。但是，在采樣點(diǎn)之間的頻響特性則是由各采樣點(diǎn)的加權(quán)內(nèi)插函數(shù)疊加形成的，因而存在一定的逼近誤差。6.3.1頻率采樣法設(shè)計(jì)原理56

誤差大小則取決于理想頻響特性的曲線形狀和采樣點(diǎn)的密度。理想頻響特性變化越平緩，則內(nèi)插值越接近理想值，逼近誤差越小。反之，如果采樣點(diǎn)之間的理想頻響特性變化越迅速，則內(nèi)插值與理想值的誤差就越大。因此，在理想頻響特性的不連續(xù)點(diǎn)附近會(huì)形成振蕩特性。采樣點(diǎn)數(shù)愈多，即采樣頻率越高，誤差越小。圖6.3-1頻率采樣的逼近效果從圖中可以看出，其一，由頻率采樣法設(shè)計(jì)所得的頻響特性在采樣點(diǎn)之間出現(xiàn)了起伏振蕩；其二，在通帶和阻帶之間不連續(xù)處，變化較劇烈，出現(xiàn)肩峰。但在通帶和阻帶之間有過渡帶時(shí)，變化比較緩慢，對的逼近較好。6.3.2頻率采樣法的設(shè)計(jì)步驟57

頻率采樣法設(shè)計(jì)數(shù)字濾波器的步驟：(1)首先根據(jù)理想濾波器的性能指標(biāo)，計(jì)算在通帶、阻帶中的采樣點(diǎn)數(shù)，確定所設(shè)計(jì)濾波器單位脈沖響應(yīng)h[n]的對稱性(奇、偶)；(2)根據(jù)單位脈沖響應(yīng)h[n]的對稱性，計(jì)算各采樣的幅度值和相位值；(3)利用理想濾波器的頻率采樣值Hd[k]=H

[k]ej(k)

，通過傅里葉逆變換IDFT，求所設(shè)計(jì)濾波器的單位脈沖響應(yīng)h[n]，即h[n]=IDFT[Hd[k]]；(4)利用DTFT變換，求濾波器的頻響特性H

(ej

)=DTFT[h[n]]，檢驗(yàn)是否滿足設(shè)計(jì)要求，若不滿足，可以在通帶和阻帶交界處安排一個(gè)或幾個(gè)不等于零的采樣過渡點(diǎn)，重復(fù)(1)，(2)，(3)步驟計(jì)算處理，直到滿足設(shè)計(jì)要求為止。6.3.2頻率采樣法的設(shè)計(jì)步驟58

1.類型1，h[n]=

h[N

1

n]

，N為奇數(shù)濾波器的幅度函數(shù)是偶對稱的，相位函數(shù)幅度及相位的采樣分別為2.類型2，

h[n]=

h[N

1

n]

，N為偶數(shù)濾波器的幅度函數(shù)是奇對稱的，相位函數(shù)幅度及相位的采樣分別為6.3.2頻率采樣法的設(shè)計(jì)步驟59

3.類型3，

h[n]=

h[N

1

n]

，N為奇數(shù)幅度及相位的采樣分別為4.類型4，

h[n]=

h[N

1

n]

，N為偶數(shù)濾波器的幅度函數(shù)是奇對稱的，相位函數(shù)幅度及相位的采樣分別為濾波器的幅度函數(shù)是偶對稱的，相位函數(shù)6.3.2頻率采樣法的設(shè)計(jì)步驟60例6.3-1用頻率采樣法設(shè)計(jì)一個(gè)線性相位FIR數(shù)字低通濾波器。要求單位脈沖響應(yīng)滿足h[n]=

h[N

1

n]

，理想濾波器幅頻特性為要求截止頻率，采樣點(diǎn)數(shù)。解：（1）確定截止頻率所處的位置因?yàn)椴蓸狱c(diǎn)數(shù)33，采樣間隔，則，

,即截頻率位于第6個(gè)采樣與第7個(gè)采樣之間，?。?）確定各采樣點(diǎn)的幅度和相位大小6.3.2頻率采樣法的設(shè)計(jì)步驟61

（3）利用IDFT求FIR數(shù)字濾波器的單位脈沖響應(yīng)

(4)最后計(jì)算所設(shè)計(jì)FIR數(shù)字濾波頻率響應(yīng)6.3.2頻率采樣法的設(shè)計(jì)步驟62為了改善濾波器的頻響特性，可以在通帶和阻帶交界處安排一個(gè)或幾個(gè)在0到1之間的采樣值。首先，在k=8和k=26處，增設(shè)H[k]=0.5，這等效于加寬過渡帶增至4π/33。重新計(jì)算出的如圖6.3-2(b),其阻帶最小衰減約為-30dB。6.3.2頻率采樣法的設(shè)計(jì)步驟63為了進(jìn)一步增加阻帶衰減，又不增加過渡帶帶寬，可以增加采樣點(diǎn)數(shù)。例如，用65點(diǎn)進(jìn)行采樣，在k=14和k=52處，插進(jìn)優(yōu)化采樣值H[k]=0.6，在k=15和k=51處，插進(jìn)優(yōu)化采樣值H[k]=0.1，過渡帶寬為6π/65，并沒有增加，而阻帶最小衰減可達(dá)-60dB以上。見圖6.3-2(c)。所付出的代價(jià)是提高了濾波器的階數(shù)，因而運(yùn)算量也隨之增加。需要說明，在總的采樣點(diǎn)不變的前提下，過渡帶的采樣值不同直接影響到濾波器的頻響特性。在圖6.3-2(b)中，過渡帶采樣點(diǎn)取0.5并非最優(yōu)值。如果采樣最優(yōu)的過渡帶采樣點(diǎn)取0.3904，最小阻帶衰減可以提高10dB。例6.3-1所設(shè)計(jì)的濾波器為類型1線性相位FIR濾波器。6.3.2頻率采樣法的設(shè)計(jì)步驟64

例6.3-2用頻率采樣法設(shè)計(jì)一個(gè)線性相位FIR帶通濾波器，設(shè)N=32，理想濾波器幅頻特性為解:

N=32為偶數(shù)，按類型2線性相位FIR濾波器設(shè)計(jì)頻率間隔為，求得6.3.2頻率采樣法的設(shè)計(jì)步驟65

根據(jù)上面設(shè)計(jì)可以得到帶通濾波器的頻率特性如圖6.3-3所示。6.3.2頻率采樣法的設(shè)計(jì)步驟66

小結(jié)：頻率采樣設(shè)計(jì)法優(yōu)點(diǎn)：①

直接從頻域進(jìn)行設(shè)計(jì)，物理概念清楚，直觀方便；②

適合于窄帶濾波器設(shè)計(jì)，這時(shí)頻率響應(yīng)只有少數(shù)幾個(gè)非零值。典型應(yīng)用：用一串窄帶濾波器組成多卜勒雷達(dá)接收機(jī)，覆蓋不同的頻段，多卜勒頻偏可反映被測目標(biāo)的運(yùn)動(dòng)速度；缺點(diǎn)：截止頻率難以控制。

因頻率取樣點(diǎn)都局限在2π/N的整數(shù)倍點(diǎn)上，所以在指定通帶和阻帶截止頻率時(shí)，這種方法受到限制，比較死板；充分加大N，可以接近任何給定的頻率，但計(jì)算量和復(fù)雜性增加。6.4FIR數(shù)字濾波器的等波紋設(shè)計(jì)67FIR濾波器設(shè)計(jì)的本質(zhì)是尋找一種長度為N的單位脈沖響應(yīng)h[n]，使h[n]的頻響特性與希望逼近的片段恒定常數(shù)濾波器頻響特性的誤差滿足指標(biāo)要求。窗函數(shù)和頻率采樣設(shè)計(jì)方法，逼近誤差在整個(gè)頻域分布極不均勻，如果在誤差最大的頻段剛好達(dá)到指標(biāo)，但是在誤差最小的頻段則遠(yuǎn)遠(yuǎn)優(yōu)于指標(biāo)，從而造成資源浪費(fèi)。

等波紋最佳一致逼近法使最大誤差最小化，并在整個(gè)逼近頻段上均勻分布。用等波紋最佳一致逼近法設(shè)計(jì)的濾波器的幅頻特性在通帶和阻帶都是等波紋的，這就是等波紋的含義，最佳一致逼近指在濾波器長度給定條件下使波紋幅度最小。由于這種設(shè)計(jì)法使誤差均勻分布，所以設(shè)計(jì)的濾波器性能價(jià)格比最高。指標(biāo)相同時(shí)這種設(shè)計(jì)方法使濾波器階數(shù)最低；階數(shù)相同時(shí)，該設(shè)計(jì)法使濾波器的最大逼近誤差最小，即通帶衰減最

小，阻帶衰減最大。68

6.4.1等波紋最佳一致設(shè)計(jì)的基本思想

最優(yōu)化設(shè)計(jì)的前提是最優(yōu)準(zhǔn)則的確定，在FIR濾波器最優(yōu)化設(shè)計(jì)中，常用的準(zhǔn)則有①最小均方誤差準(zhǔn)則②最大誤差最小化準(zhǔn)則。（1）

均方誤差最小化準(zhǔn)則，若以E(ejω)表示逼近誤差，則

均方誤差為69

6.4.1等波紋最佳一致設(shè)計(jì)的基本思想

均方誤差最小準(zhǔn)則就是選擇一組時(shí)域采樣值，以使均方誤差，這一方法注重的是在整個(gè)-π～π頻率區(qū)間內(nèi)總誤差的全局最小，但不能保證局部頻率點(diǎn)的性能，有些頻率點(diǎn)可能會(huì)有較大的誤差，對于窗口法FIR濾波器設(shè)計(jì)，因采用有限項(xiàng)的h[n]逼近理想的hd[n]，所以其逼近誤差為：如果采用矩形窗706.4.1等波紋最佳一致設(shè)計(jì)的基本思想可以證明，這是一個(gè)最小均方誤差。所以，矩形窗窗口設(shè)計(jì)法是一個(gè)最小均方誤差FIR設(shè)計(jì)，根據(jù)前面的討論，我們知道其優(yōu)點(diǎn)是過渡帶較窄，缺點(diǎn)是局部點(diǎn)誤差大，或者說誤差分布不均勻。（2）

最大誤差最小化準(zhǔn)則（也叫最佳一致逼近準(zhǔn)則）表示希望逼近的幅度函數(shù)；表示實(shí)際設(shè)計(jì)的濾波器幅度函數(shù)；稱為誤差加權(quán)函數(shù)，用來控制不同頻段(一般指通帶和阻帶)的逼近精度。取值越大，表示相應(yīng)頻段逼近精度越高。等波紋最佳一致逼近法就是在通帶和阻帶以|E(

)|的最大值最小化為準(zhǔn)則，求解濾波器脈沖響應(yīng)h[n]。把通帶和阻帶稱為逼近(或研究)區(qū)域，把過渡帶稱為無關(guān)區(qū)域。應(yīng)當(dāng)注意，設(shè)計(jì)過程中無關(guān)區(qū)寬度不能為零，即不能是理想幅度函數(shù)。71

6.4.1等波紋最佳一致設(shè)計(jì)的基本思想某濾波器的參數(shù)：等波紋濾波器幅頻特性72

6.4.1等波紋最佳一致設(shè)計(jì)的基本思想設(shè)希望逼近濾波器的通帶和阻帶分別為[0,

/4]和[5

/16，

]，對4種不同的性能指標(biāo)參數(shù)，等波紋最佳一致逼近的增益響應(yīng)曲線分別如圖(a)、(b)、(c)和(d)所示。N=33,W[1,10]；(b)N=33,W[10,1]；(c)N=33,W[1,1]；(d)N=63,W[1,1]誤差加權(quán)函數(shù)W(

)和濾波器階數(shù)N對逼近精度的影響73

6.4.2remez函數(shù)與remezord函數(shù)remez函數(shù)與remezord函數(shù)是Matlab工具箱函數(shù)的庫函數(shù)，編程時(shí)可以直接調(diào)用，下面簡單說明其調(diào)用方法。1．remez函數(shù)功能：采用remez算法實(shí)現(xiàn)線性相位FIR數(shù)字濾波器的等波紋最佳一致逼近設(shè)計(jì)。調(diào)用格式：hn=remez(N,f,m,w,’ftype’)返回單位脈沖響應(yīng)向量hn。調(diào)用參數(shù)含義如下：N：FIRDF階數(shù)，設(shè)計(jì)結(jié)果hn長度為N+1；f、m：給出希望逼近的幅特性。f為邊界頻率向量，0≤f≤1，要求f為單調(diào)增向量。

調(diào)用參數(shù)f和m含義圖解w為誤差加權(quán)向量，其長度為f的一半。ftype用于指定濾波器類型。74

6.4.2remez函數(shù)與remezord函數(shù)2．remezord函數(shù)

功能：用于估算FIR數(shù)字濾波器的等波紋最佳一致逼近設(shè)計(jì)的最低階數(shù)N，誤差加權(quán)向量w，歸一化邊界頻率向量fo，從而使濾波器在滿足指標(biāo)的前提下，濾波器階數(shù)最低。其返回參數(shù)作為remez函數(shù)的調(diào)用參數(shù)。調(diào)用格式：[N，fo，mo，w]=remezord(f，m，rip，F(xiàn)s)參數(shù)說明：

返回參數(shù)作為remez函數(shù)的調(diào)用參數(shù)，F(xiàn)s為采樣頻率，省略時(shí)默認(rèn)為2Hz。解：調(diào)用remezord和remez函數(shù)求解。首先要根據(jù)性能指標(biāo)確定remezord函數(shù)的調(diào)用參數(shù)，再直接寫程序調(diào)用remezord和remez函數(shù)得到濾波器的脈沖響應(yīng)h[n]。75

6.4.2remez函數(shù)與remezord函數(shù)76

6.4.2remez函數(shù)與remezord函數(shù)用上述參數(shù)可以設(shè)計(jì)帶通濾波器，如圖6.4-4所示中，N

1=23，即h[n]長度為24，顯示了h[n]序列和增益響應(yīng)曲線。等波紋最佳一致逼近設(shè)計(jì)方法可以使濾波器階數(shù)大大降低。例題可以看出，h[n]滿足偶對稱條件，所以必然具有嚴(yán)格線性相位特性。77

FIRIIR性能嚴(yán)格線性相位，但成本高，運(yùn)算量大，信號延時(shí)也較大存儲(chǔ)單元少，運(yùn)算量小，經(jīng)濟(jì)高效設(shè)計(jì)結(jié)果可得到幅頻特性（可以多帶）和線性相位（最大優(yōu)點(diǎn)）只能得到幅頻特性，相頻特性未知（一大缺點(diǎn)），如需要線性相位，須用全通網(wǎng)絡(luò)校準(zhǔn)，但增加濾波器階數(shù)和復(fù)雜性穩(wěn)定性極點(diǎn)全部在原點(diǎn)（永遠(yuǎn)穩(wěn)定）無穩(wěn)定性問題有穩(wěn)定性問題階數(shù)高

結(jié)構(gòu)非遞歸遞歸結(jié)構(gòu)設(shè)計(jì)工具只有計(jì)算程序可循可借助于AF的設(shè)計(jì)成果快速算法可用FFT實(shí)現(xiàn)，減少運(yùn)算量無快速運(yùn)算方法6.5FIR濾波器與IIR濾波器的比較78

6.6用Matlab函數(shù)設(shè)計(jì)FIR數(shù)字濾波器例6.6-1已知FIR濾波器的單位脈沖響應(yīng)為：，畫出該濾波器的頻響特性，零極點(diǎn)分布圖。解：根據(jù)濾波器的單位脈沖響應(yīng)可知，其h[n]具有對稱性，并且長度為12，所以該濾波器是類型2線性相位濾波器。b=[2,1,-1,3,5,-4,-4,5,3,-1,1,2];n=0:length(b)-1;a=[1];th=0:0.05:2*pi;[H,w]=freqz(b,a);zerophase(b,a,512,'whole');G=tf(b,a);subplot(1,4,1);stem(n,b);xlabel('n');ylabel('h[n]');title('單位脈沖響應(yīng)');subplot(1,4,2);plot(w/pi,abs(H));xlabel('\omega/\pi');ylabel('(dB)');title('濾波器幅頻特性');subplot(1,4,3);plot(w/pi,angle(H)/pi);xlabel('\omega/\pi');ylabel('(rad/\pi)');title('濾波器相頻特性');subplot(1,4,4);pzmap(G);holdon;;title('濾波器零極點(diǎn)分布圖');x=cos(th);y=sin(th);plot(x,y,':');%畫單位圓6.6用Matlab函數(shù)設(shè)計(jì)FIR數(shù)字濾波器例6.6-3設(shè)計(jì)具有任意頻率響應(yīng)的FIR濾波器。本例設(shè)計(jì)一個(gè)“多段濾波器”，它有三個(gè)頻帶，各自具有不同的理想幅度和不同的誤差波動(dòng)，設(shè)計(jì)參數(shù)為：頻帶1：，理想幅度為0.3；頻帶2：，理想幅度為0.7；頻帶3：，理想幅度為1。畫出濾波器的幅頻特性。解：Matlab提供兩個(gè)窗函數(shù)設(shè)計(jì)FIR濾波器的函數(shù)。fir1可以設(shè)計(jì)標(biāo)準(zhǔn)頻率響應(yīng)的濾波器；fir2可以設(shè)計(jì)任意頻率響應(yīng)的多帶濾波器。本例利用fir2函數(shù)設(shè)計(jì)FIR濾波器，使用的窗函數(shù)為Hamming缺省值，fpts=[00.30.40.650.751];mval=[0.30.30.70.71.01.0];b=fir2(100,fpts,mval);[h,omega]=freqz(b,1,512);plot(omega/pi,abs(h));xlabel('\omega/\pi');ylabel('幅度');title('濾波器幅頻特性');gridon;濾波器的幅頻特性呈現(xiàn)多帶特點(diǎn)，該濾波器不屬于標(biāo)準(zhǔn)的低通、高通等濾波特性。

796.6用Matlab函數(shù)設(shè)計(jì)FIR數(shù)字濾波器例6.6-4用等波紋切比雪夫法，設(shè)計(jì)一個(gè)線性相位數(shù)字帶通濾波器，設(shè)N=36，三段不同頻帶的加權(quán)值分別為1，2和0.5，邊界頻率如下:解：在窗函數(shù)設(shè)計(jì)中，不能精確指定通帶與阻帶頻率，并且通帶與阻帶區(qū)間上誤差不是均勻分布。靠近邊緣處誤差大；反之，遠(yuǎn)離邊緣處誤差小。等波紋切比雪夫設(shè)計(jì)法是通帶與阻帶上誤差均勻分布，因而可以得到較低階數(shù)的濾波器。f=[00.20.30.60.71];m=[001100];N=36;W=[120.5];b=remez(N-1,f,m,W);[h,omega]=freqz(b,1,512);%UsingHammingWindowsubplot(1,2,1);plot(omega/pi,abs(h));xlabel('\omega/\pi');ylabel('幅度');title('濾波器幅頻特性');axis([0101.1]);gridon;subplot(1,2,2);plot(omega/pi,20*log10(abs(h)));xlabel('\omega/\pi');ylabel('(dB)');itle('濾波器增益響應(yīng)');axis([01-701]);gridon;806.6用Matlab函數(shù)設(shè)計(jì)FIR數(shù)字濾波器f=[00.20.30.60.71];m=[001100];N=36;W=[120.5];b=remez(N-1,f,m,W);[h,omega]=freqz(b,1,512);%UsingHammingWindowsubplot(1,2,1);plot(omega/pi,abs(h));xlabel('\omega/\pi');ylabel('幅度');title('濾波器幅頻特性');axis([0101.1]);gridon;subplot(1,2,2);plot(omega/pi,20*log10(abs(h)));xlabel('\omega/\pi');ylabel('(dB)');itle('濾波器增益響應(yīng)');axis([01-701]);gridon;由于權(quán)系數(shù)不同，濾波器高頻段誤差較大。81小結(jié)線性相位FIR濾波器線性相位濾波器的條件4種類型線性相位濾波器的特性窗函數(shù)設(shè)計(jì)FIR濾波器設(shè)計(jì)步驟窗函數(shù)類型與特性頻率響應(yīng)設(shè)計(jì)FIR濾波器設(shè)計(jì)原理設(shè)計(jì)步驟822024/10/28dsp-chap7-201883第7章

數(shù)字信號處理系統(tǒng)的實(shí)現(xiàn)7.4數(shù)字濾波器的系數(shù)量化效應(yīng)7.3輸入信號的量化誤差7.2數(shù)的表示與量化誤差7.1數(shù)字濾波器結(jié)構(gòu)7.6FFT實(shí)現(xiàn)中的量化效應(yīng)7.5數(shù)字濾波器的運(yùn)算量化效應(yīng)7.7用Matlab分析數(shù)字信號處理系統(tǒng)與量化效應(yīng)2024/10/28dsp-chap7-201884學(xué)習(xí)目標(biāo)：熟練應(yīng)用系統(tǒng)框圖或信號流圖描述數(shù)字信號處理的過程

（a）IIR數(shù)字濾波器的三種結(jié)構(gòu)：直接型，級聯(lián)型和并聯(lián)型；

（b）FIR數(shù)字濾波器的結(jié)構(gòu)：直接型，級聯(lián)型，線性相位型，

頻率取樣型，頻域快速卷積型；第7章

數(shù)字信號處理系統(tǒng)的實(shí)現(xiàn)2.了解數(shù)字信號處理中有限字長造成的影響和效應(yīng)

（a）數(shù)的二進(jìn)制表示法；

（b）A/D轉(zhuǎn)換器的量化誤差及量化效應(yīng)統(tǒng)計(jì)分析；（c）數(shù)字濾波器系數(shù)量化所產(chǎn)生的誤差及統(tǒng)計(jì)分析模型；（d）數(shù)字濾波器實(shí)現(xiàn)時(shí)的運(yùn)算量化誤差及統(tǒng)計(jì)分析模型；

極限環(huán)效應(yīng)（e）FFT實(shí)現(xiàn)中的誤差分析2024/10/28dsp-chap7-201885數(shù)字濾波器可以用差分方程、單位脈沖響應(yīng)和系統(tǒng)函數(shù)來描述，但對于研究系統(tǒng)的實(shí)現(xiàn)，即它的運(yùn)算結(jié)構(gòu)來說，用方框圖（或信號流圖）表示最為直接。7.1數(shù)字濾波器結(jié)構(gòu)E1(z1)x[n]x[n1]單位延時(shí)器x1[n]

x2[n]x1[n]+x2[n]加法器x[n]aax[n]數(shù)乘器X(z)Y(z)=z

1X(z)z

1X(z)Y(z)=aX(z)aX1(z)Y(z)=X1(z)+X2(z)X2(z)2024/10/28dsp-chap7-201886例1某離散系統(tǒng)的信號流圖如下圖所示，分析系統(tǒng)的系統(tǒng)函數(shù)，

并求解其沖激響應(yīng)。7.1數(shù)字濾波器結(jié)構(gòu)z1Y(z)W(z)z1z1423X(z)53Y(z)解：W(z)=X(z)

(3z1+5z2+3z3)W(z)Y(z)=(2z2+4z3)W(z)如果收斂域：2024/10/28dsp-chap7-201887IIR數(shù)字濾波器：其沖激響應(yīng)h[n]是無限長度序列

；

或者系統(tǒng)函數(shù)H(z)至少有一個(gè)非零的極點(diǎn)。7.1.2IIR數(shù)字濾波器結(jié)構(gòu)IIR數(shù)字濾波器的結(jié)構(gòu)主要有三種形式：直接型，級聯(lián)型，并聯(lián)型。2024/10/28dsp-chap7-2018887.1.2IIR數(shù)字濾波器結(jié)構(gòu)1、IIR數(shù)字濾波器的直接I型結(jié)構(gòu)（設(shè)M=N）圖7.1-3直接I型結(jié)構(gòu)方框圖

b0x[n]

z

1

b1

z

1

b2

z

1bM

1++++

bM++

z

1

a1y[n]+

z

1

a2

z

1+

aN

1

aN圖7.1-4直接I型信號流圖aNy[n]a1z

1z

1z

1a2aN

1x[n]b0b1z

1z

1z

1b2bM

1bM特點(diǎn)：需要M+N個(gè)延時(shí)器

（存儲(chǔ)器）2024/10/28dsp-chap7-2018897.1.2IIR數(shù)字濾波器結(jié)構(gòu)2、IIR數(shù)字濾波器的直接II型結(jié)構(gòu)（設(shè)M=N）直接II性結(jié)構(gòu)x[n]y[n]b0b1a1z

1z

1z

1b2bN

1bNa2aN

1aNW(z)特點(diǎn)：需要max（M，N）

個(gè)延時(shí)器（存儲(chǔ)器）2024/10/28dsp-chap7-2018907.1.2IIR數(shù)字濾波器結(jié)構(gòu)3、轉(zhuǎn)置型結(jié)構(gòu)流圖將方框圖或信號流圖中所有支路方向顛倒成反向，且輸入輸出的位置互相調(diào)換一下，得到系統(tǒng)的轉(zhuǎn)置流圖。以前述的直接II型為例：圖7.1-6轉(zhuǎn)置的直接II型信號流圖

x[n]y[n]b0b1a1z

1z

1z

1b2bM

1bMa2aN

1aN2024/10/28dsp-chap7-2018917.1.2IIR數(shù)字濾波器結(jié)構(gòu)3、轉(zhuǎn)置型結(jié)構(gòu)流圖直接I、II型結(jié)構(gòu)都稱為直接型結(jié)構(gòu)。其優(yōu)點(diǎn)是簡單直觀。缺點(diǎn)是：系數(shù)ai和bi對濾波器性能（零點(diǎn)和極點(diǎn)）的控制關(guān)系不直接，也就是說，當(dāng)系數(shù)ai中有一個(gè)發(fā)生變化，則所有極點(diǎn)位置都會(huì)變化，系數(shù)bi對零點(diǎn)的影響也是如此。而且階數(shù)N越大，影響也越大。所以高階的IIR濾波器一般會(huì)采用級聯(lián)或并聯(lián)結(jié)構(gòu)以減小上述影響。2024/10/28dsp-chap7-2018927.1.2IIR數(shù)字濾波器結(jié)構(gòu)4、IIR數(shù)字濾波器的級聯(lián)型結(jié)構(gòu)圖7.1-7級聯(lián)型的二階基本節(jié)Hi(z)1

1i

1iz

1z

1

2i

2i

H1(z)

H2(z)

HL(z)

y[n]

x[n]

A圖7.1-8IIR數(shù)字濾波器的級聯(lián)型結(jié)構(gòu)2024/10/28dsp-chap7-2018937.1.2IIR數(shù)字濾波器結(jié)構(gòu)5、IIR數(shù)字濾波器的并聯(lián)型結(jié)構(gòu)r0ir1i

1iz

1z

1

2iAigiz

1

y[n]

x[n]

C0圖7.1-9IIR數(shù)字濾波器的并聯(lián)型結(jié)構(gòu)H1(z)

H2(z)

HP(z)

2024/10/28dsp-chap7-2018947.1.2IIR數(shù)字濾波器結(jié)構(gòu)例7.1-1分別畫出下式所描述系統(tǒng)的直接II型、并聯(lián)型和級聯(lián)型結(jié)構(gòu)。解：(1)直接II型，如右圖所示120.75z

1z

10.125x[n]y[n]1(2)并聯(lián)型：x[n]0.25z

1y[n]

0.5z

1

251882024/10/28dsp-chap7-2018957.1.2IIR數(shù)字濾波器結(jié)構(gòu)(2)并聯(lián)型，也可以將上式的并聯(lián)型變換一下：x[n]0.25z

1y[n]

0.5z

1

2526

4(3)級聯(lián)型：x[n]0.25z

1y[n]

0.5z

1注意：（a）級聯(lián)和并聯(lián)的子系統(tǒng)階數(shù)必須低于整個(gè)系統(tǒng)的階數(shù)；

（b）級聯(lián)和并聯(lián)子系統(tǒng)的各項(xiàng)系數(shù)必須為實(shí)數(shù)。2024/10/28dsp-chap7-201896FIR數(shù)字濾波器：其沖激響應(yīng)h[n]是有限長度序列，

通常取0

n

N

1

；

從而其系統(tǒng)函數(shù)H(z)的極點(diǎn)全是0。7.1.3FIR數(shù)字濾波器結(jié)構(gòu)另外，如果FIR數(shù)字濾波器的單位沖激響應(yīng)滿足對稱或反對稱特性，即h[n]=

h[N

1

n]，則FIR濾波器的相位響應(yīng)具有線性特點(diǎn)。直接型，級聯(lián)型，線性相位型，快速卷積型，頻率取樣型。所以，F(xiàn)IR濾波器一般具有如下結(jié)構(gòu)：2024/10/28dsp-chap7-2018971、IIR數(shù)字濾波器的直接I型結(jié)構(gòu)（設(shè)M=N）

h[0]x[n]

z

1

h[1]

z

1h[2]

z

1h[N2]++++

h[N1]y[n]y[n]x[n]h[0]h[1]z

1z

1z

1h[2]h[N2]h[N1]7.1.3FIR數(shù)字濾波器結(jié)構(gòu)y[n]x[n]h[0]h[1]z

1z

1z

1h[2]h[N2]h[N1]2024/10/28dsp-chap7-2018982、FIR數(shù)字濾波器的級聯(lián)型結(jié)構(gòu)

1iz

1z

1

2i

0i

H1(z)

H2(z)

HL(z)

y[n]

x[n]

A7.1.3FIR數(shù)字濾波器結(jié)構(gòu)圖7.1-16FIR濾波器的級聯(lián)型結(jié)構(gòu)x[n]

01z

1z

1z

1

11

21z

1

02

12

22y[n]

z

1z

1

2L

1L

0L2024/10/28dsp-chap7-2018993、FIR線性相位數(shù)字濾波器的線性相位性結(jié)構(gòu)7.1.3FIR數(shù)字濾波器結(jié)構(gòu)h[n]=

h[N

1

n]，0

n

N

1當(dāng)N為奇數(shù)時(shí)，當(dāng)N為偶數(shù)時(shí)，2024/10/28dsp-chap7-20181003、FIR線性相位數(shù)字濾波器的線性相位性結(jié)構(gòu)7.1.3FIR數(shù)字濾波器結(jié)構(gòu)h[n]=

h[N

1

n]，0

n

N

1當(dāng)N為奇數(shù)時(shí)，當(dāng)N為偶數(shù)時(shí)，x[n]z

1z

1h[0]z

1

1y[n]

z

1z

1(a)z

1z

1z

1h[1]h[2]h[(N

1)/2]h[(N

3)/2]

1

1

1x[n]z

1z

1h[0]z

1

1y[n]

z

1z

1(b)z

1z

1z

1h[1]h[2]

1

1

1z

1

1h[(N/2)2]h[(N/2)1]2024/10/28dsp-chap7-20181014、FIR數(shù)字濾波器的快速卷積型結(jié)構(gòu)7.1.3FIR數(shù)字濾波器結(jié)構(gòu)假設(shè)輸入信號x[n]長度為N1，如果取N

N1+N2

1，并假設(shè)x[n]和h[n]的N點(diǎn)DFT分別為X[k]和H[k]，則根據(jù)DFT的圓周卷積定理，有y[n]

Y[k]=H[k]X[k]也即：

y[n]=IDFT[H[k]X[k]]x[n]0

n

N1

1補(bǔ)零N點(diǎn)DFTh[n]0

n

N2

1he[n]0

n

N

1H[k]0

k

N

1

N點(diǎn)IDFTy[n]0

n

N1+N2

2補(bǔ)零N點(diǎn)DFTxe[n]0

n

N

1X[k]0

k

N

1H[k]X[k]0

k

N

12024/10/28dsp-chap7-20181025、FIR數(shù)字濾波器的頻率采樣型結(jié)構(gòu)7.1.3FIR數(shù)字濾波器結(jié)構(gòu)設(shè)h[n]的N點(diǎn)DFT為H[k]，則圖7.1-18FIR濾波器的頻率采樣型結(jié)構(gòu)x[n]y[n]

z

NH[N

1]z

1WN

(N

1)H[1]z

1WN

1H[0]z

1WN01/N（a）原理：2024/10/28dsp-chap7-20181035、FIR數(shù)字濾波器的頻率采樣型結(jié)構(gòu)7.1.3FIR數(shù)字濾波器結(jié)構(gòu)（b）特點(diǎn)：（1）Hc(z)是一個(gè)由N節(jié)延時(shí)單元組成的梳狀濾波器，它在單位圓上有N個(gè)等分的零點(diǎn)，而一階濾波器Hk(z)在單位圓上有一個(gè)對應(yīng)的極點(diǎn)；（2）系數(shù)H[k]和WN

k都是復(fù)數(shù)，增加了運(yùn)算的復(fù)雜性；（3）所有諧振器的極點(diǎn)都在單位圓上，由于量化效應(yīng)的影響，系統(tǒng)的穩(wěn)定性不能保證。2024/10/28dsp-chap7-20181045、FIR數(shù)字濾波器的頻率采樣型結(jié)構(gòu)7.1.3FIR數(shù)字濾波器結(jié)構(gòu)（c）修正的頻率采樣型結(jié)構(gòu)：當(dāng)N為偶數(shù)時(shí)，當(dāng)N為奇數(shù)時(shí)，其中：r0k=2Re[H[k]]，r1k=

2Re[rH[k]WNk]2024/10/28dsp-chap7-2018105例7.1-2已知FIR濾波器的系統(tǒng)函數(shù)為：

H(z)=(1+0.5z

1)(1+2z

1)(1

0.25z

1)(1

4z

1)

分別畫出它的直接型、級聯(lián)型、線性相位型和頻率采樣型結(jié)構(gòu)。解：(1)直接型，如右圖所示(2)級聯(lián)型（一階節(jié)級聯(lián)）：7.1.3FIR數(shù)字濾波器結(jié)構(gòu)H(z)=1

1.75z

1

8.625z

2

1.75z

3+z

4x[n]z

1y[n]

z

1z

1z

1

1.75

8.625

1.75(a)

x[n]z

1z

10.52y[n]

z

1

4z

1

0.25(b-1)2024/10/28dsp-chap7-2018106(2)級聯(lián)型，也可以將上式的級聯(lián)型變換為二階節(jié)的級聯(lián)：H(z)=(1+2.5z

1+z

2)(1

4.25z

1+z

2)----線性相位7.1.3FIR數(shù)字濾波器結(jié)構(gòu)(b-2)x[n]z

1z

1z

12.5z

1

4.25y[n]

(3)線性相位型結(jié)構(gòu)：：H(z)=(1+z

4)

1.75(z

1+z

3)

8.625z

2(c)x[n]z

1z

1

y[n]

z

1

1.75

8.625

z

12024/10/28dsp-chap7-2018107(4)頻率采樣型結(jié)構(gòu)：取H(z)的8等分點(diǎn)采樣值為：H[k]={

10.125，j11.10，10.625，

j6.15，

3.125，j6.15，10.625，

j11.10}，0

k

77.1.3FIR數(shù)字濾波器結(jié)構(gòu)設(shè)修正半徑r=1（即不修正極點(diǎn)位置）y[n]

x[n]

z

81/821.251.414

115.695

8.696z

1z

1

1z

1z

1

1

3.125z

1

10.125z

1

1.414

1z

1z

12024/10/28dsp-chap7-2018108數(shù)的記錄與表示：7.2數(shù)的表示與量化誤差如：123.456=1102+2101+3100+4101+5102+6103日常記數(shù)法：十進(jìn)制，包含十個(gè)不同的數(shù)碼0，1，2，…，9；基：10；高位低位整數(shù)小數(shù)（尾數(shù)）權(quán)位數(shù)字信號處理中：二進(jìn)制，包含二個(gè)不同的數(shù)碼0，1；基：2；如：11001.1101=124+123+022+021+120

+121+122+023+124為縮短二進(jìn)制表示的長度或快速讀寫二進(jìn)制數(shù)：八進(jìn)制，十六進(jìn)制如：11001.1101=（31.64）8=（1A.E）162024/10/28dsp-chap7-20181091.定點(diǎn)二進(jìn)制數(shù)：小數(shù)點(diǎn)位置固定不動(dòng)7.2.1二進(jìn)制數(shù)的表示再如：（11001.1101）2=（25.8125）10，用定點(diǎn)小數(shù)（二進(jìn)制）表示：（25.8125/32）10=（0.806640625）10=（0.110011101）2但為了運(yùn)算方便，通常定點(diǎn)制都把數(shù)

a限制在

1之間，即

1<a<1。這時(shí)，小數(shù)點(diǎn)固定在分?jǐn)?shù)的第一位二進(jìn)碼之前，而整數(shù)位定義為“符號位”，代表數(shù)的正負(fù)(“0”表示正數(shù)，“1”表示負(fù)數(shù))，小數(shù)部分稱為“尾數(shù)”。如：（0.101）2=（0.625）10定點(diǎn)數(shù)有三種編碼方式：（1）原碼：也稱“符號-幅度碼”；（2）補(bǔ)碼：也稱“2的補(bǔ)碼”；（3）反碼：也稱“1的補(bǔ)碼”；三種編碼對于正數(shù)來說，完全相同；但是對于負(fù)數(shù)，則有不同的表示。例：（

0.625）10=（

1.101）2原=（

1.011）2補(bǔ)=（

1.010）2反2024/10/28dsp-chap7-20181101.定點(diǎn)二進(jìn)制數(shù)：小數(shù)點(diǎn)位置固定不動(dòng)7.2.1二進(jìn)制數(shù)的表示特點(diǎn)：定點(diǎn)數(shù)做加減法運(yùn)算其結(jié)果可能或超出

1的范圍，稱為“溢出”。定點(diǎn)數(shù)做乘法運(yùn)算不會(huì)造成溢出，但是字長卻要增加一倍，一般在定點(diǎn)乘法運(yùn)算以后需要對尾數(shù)做截尾或舍入處理，以保證字長不變，但是這樣處理后會(huì)帶來截尾誤差或舍入誤差。2.浮點(diǎn)二進(jìn)制數(shù)：x=M

2c，M是它的尾數(shù)部分，2c是它的指數(shù)部分，c是階數(shù)，稱為“階碼”。尾數(shù)和階碼都用帶符號的定點(diǎn)數(shù)表示，通常M的范圍是0.5

|M|<1特點(diǎn)：浮點(diǎn)數(shù)做加減法運(yùn)算：需要對階，然后再做加減；其動(dòng)態(tài)范圍很大，不會(huì)出現(xiàn)“溢出”。浮點(diǎn)數(shù)做乘法運(yùn)算時(shí)：階碼相加，尾數(shù)相乘但不論做加減或乘法運(yùn)算，都存在尾數(shù)的截尾或舍入處理而帶來誤差的問題。2024/10/28dsp-chap7-20181111.截尾誤差7.2.2定點(diǎn)制的量化誤差當(dāng)定點(diǎn)系統(tǒng)采取截尾處理時(shí)，對于正數(shù)，三種碼的表示法是相同的，因而量化的影響也是相同的。設(shè)

x是(b1+1)位的正數(shù)：即其十進(jìn)制的數(shù)值為：截尾后字長為b位，顯然b<b1，對x截尾后的量化值用QT[x]表示為截尾的量化誤差用ET表示，為當(dāng)

i=0，i=b+1，b+2，

，b1時(shí)，取等號，即ET=0當(dāng)

i=1，i=b+1，b+2，

，b1時(shí)，誤差達(dá)到最大值，即2024/10/28dsp-chap7-20181121.截尾誤差7.2.2定點(diǎn)制的量化誤差對于負(fù)數(shù)，三種碼的表示法各不相同，因而量化產(chǎn)生的誤差也各不相同。設(shè)

x是(b1+1)位的負(fù)數(shù)。（a）負(fù)數(shù)原碼：截尾的量化誤差用ET表示，為令q=2

b，它是(b+1)位(含符號位)字長正數(shù)的最小值，這個(gè)值稱為“量化寬度”或“量化階”。定點(diǎn)正數(shù)的截尾誤差為：

q<ET

0定點(diǎn)負(fù)數(shù)原碼的截尾誤差為：

0

ET<q2024/10/28dsp-chap7-20181131.截尾誤差7.2.2定點(diǎn)制的量化誤差對于負(fù)數(shù)，三種碼的表示法各不相同，因而量化產(chǎn)生的誤差也各不相同。設(shè)

x是(b1+1)位的負(fù)數(shù)。（b）負(fù)數(shù)反碼：截尾的量化誤差用ET表示，為定點(diǎn)負(fù)數(shù)反碼的截尾誤差為：

0

ET<q2024/10/28dsp-chap7-20181141.截尾誤差7.2.2定點(diǎn)制的量化誤差對于負(fù)數(shù)，三種碼的表示法各不相同，因而量化產(chǎn)生的誤差也各不相同。設(shè)

x是(b1+1)位的負(fù)數(shù)。（c）負(fù)數(shù)補(bǔ)碼：截尾的量化誤差用ET表示，為定點(diǎn)負(fù)數(shù)反碼的截尾誤差為：

q<ET

0歸結(jié)起來，定點(diǎn)表示法的截尾誤差為：正數(shù)及負(fù)數(shù)補(bǔ)碼的截尾誤差為負(fù)數(shù)，其范圍：

q<ET

0；負(fù)數(shù)原碼與負(fù)數(shù)反碼的截尾誤差為正數(shù)，其范圍：0

ET<q。2024/10/28dsp-chap7-20181152.舍入誤差7.2.2定點(diǎn)制的量化誤差對字長為(b1+1)位的定點(diǎn)數(shù)作舍入處理時(shí)，是通過對尾數(shù)的第(b+1)位加1，然后截取到b位實(shí)現(xiàn)的。舍入之后的量化間距為q=2

b，即兩個(gè)數(shù)間最小非零差是q。無論是原碼、反碼還是補(bǔ)碼，其誤差總是在

q/2之間。以QR[x]表示對x作舍入處理，ER表示舍入誤差，則有

q/2<ER

q/2例7.2-1

若寄存器字長b=2的數(shù)字系統(tǒng)，當(dāng)經(jīng)過運(yùn)算處理后字長增加到b1=4?，F(xiàn)有運(yùn)算后的三個(gè)數(shù)：

負(fù)數(shù)原碼x1B=1.1001，負(fù)數(shù)反碼x2反=1.1010，和負(fù)數(shù)補(bǔ)碼x3補(bǔ)=1.1010（1）若采用截尾處理，分別求出上述三數(shù)所引起的誤差；（2）若采用舍入處理，分別求出上述三數(shù)所引起的誤