《數(shù)字信號處理》課件第7章

第7章數(shù)字信號處理中的有限字長效應(yīng)7.1

數(shù)值表示的有限字長效應(yīng)7.2

A/D變換的有限字長效應(yīng)7.3數(shù)字濾波器系數(shù)量化的有限字長效應(yīng)7.4定點運算對數(shù)字濾波器的影響

7.1數(shù)值表示的有限字長效應(yīng)7.1.1定點數(shù)與浮點數(shù)為了用有限的數(shù)字符號代表所有的數(shù)值，人們通常采用進位制（也稱進制）的方法，即按給定的規(guī)則進位。例如，對于R進制，就表示某一位置上的數(shù)運算時每逢R進一位。R進制下，任何一個數(shù)P可以表示為(7-1)當式（7-1）中的下限m≥0時，P是純整數(shù)，上限n≤0時，P是純小數(shù)?？梢钥闯?，i=0是區(qū)分數(shù)P的整數(shù)部分和小數(shù)部分的關(guān)鍵位置，稱為小數(shù)點位置。根據(jù)小數(shù)點位置是否變化，將數(shù)的表示區(qū)分為定點數(shù)和浮點數(shù)兩種。小數(shù)點位置固定不變的數(shù)稱為定點數(shù)，如11.556，3.189等。定點數(shù)表示數(shù)的格式規(guī)范，他所表示的數(shù)的精度固定。小數(shù)點位置可以變化的數(shù)稱為浮點數(shù)。浮點數(shù)所表示數(shù)的精度是不同的，如111.76，1.2357，1000.1等。同一個數(shù)也可以用不同的形式來表示，如111.76，可以用科學計數(shù)法表示成11.176×101，1.1176×102，0.11176×103等。數(shù)字系統(tǒng)中，所采用的二進制表示法有定點制和浮點制兩種。無論是采用軟件編程還是采用專用硬件實現(xiàn)，存儲和處理過程中所采用的存儲單元的長度都只能是有限的，也就是說數(shù)值表達的精度是有限的。接下來對二進制的定點、浮點表示因寄存器長度限制而導致的誤差進行分析，并就折中的成組浮點制（BlockFloatingPoint，BFP，也稱塊浮點制）進行簡單討論。7.1.2定點制誤差分析

1.數(shù)的定點表示定點制下，一旦確定了小數(shù)點在整個數(shù)碼中的位置，在整個運算過程中即保持不變。因此，根據(jù)系統(tǒng)設(shè)計要求、數(shù)值范圍來確定小數(shù)點處于什么位置很重要，這就是數(shù)的定標。數(shù)的定標有Q表示法和S表示法兩種。Q表示法形如Qn，字母Q后的數(shù)值n表示包含n位小數(shù)。如Q0表示小數(shù)點在第0位的后面，數(shù)為整數(shù)；Q15表示小數(shù)點在第15位的后面，0～14位都是小數(shù)位。S表示法則形如Sm.n，m表示整數(shù)位，n表示小數(shù)位。以16位DSP為例，通過設(shè)定小數(shù)點在16位數(shù)中的不同位置，可以表示不同大小和不同精度的小數(shù)。表7.1列出了一個16位數(shù)的16種Q表示、S表示及它們所能表示的十進制數(shù)值范圍。

2.定點運算定點表示的兩個數(shù)在進行加減法運算前，必須保證這兩個數(shù)的定標值Q嚴格相等。假設(shè)進行加減運算的兩個數(shù)分別為x和y，它們的Q值相等且都為Qc，則進行加減運算的結(jié)果為x±y，結(jié)果的Q值仍為Qc。當x和y的Q值不相等時，假設(shè)它們的Q值分別為Qx和Qy，且有Qx>Qy，運算結(jié)果z的定標值為Qz，進行加減運算的步驟如下：

(1)將Q值較小的數(shù)y的Q值調(diào)整為Qx，即y′=y×2(Qx－Qy)；

(2)計算z′=x±y′的值，其Q值為Qx；

(3)將z′的Q值調(diào)整到Qz，即：z=z′×2(Qz－Qx)，得到最終的運算結(jié)果。乘除運算時，假設(shè)進行運算的兩個數(shù)分別為x和y，它們的Q值分別為Qx和Qy，則兩者進行乘法運算的結(jié)果為xy，Q值為Qx+Qy，除法運算的結(jié)果為x/y，Q值為Qx－Qy。在程序或硬件實現(xiàn)中，上述定標值的調(diào)整可以直接通過寄存器的左移或右移完成。若b>0，實現(xiàn)x×2b需將存儲x的寄存器左移b位；若b<0，實現(xiàn)x×2b則需將存儲x的寄存器右移|b|位即可。

3.定點數(shù)的原碼、補碼和反碼表示從前面的分析過程可知，采用定點處理器進行小數(shù)的運算，其小數(shù)點位是由程序員根據(jù)處理問題的需要人為設(shè)定的，在處理的過程中需要時刻關(guān)注Q值的變化及其所表示的物理含義，因而，從處理器的角度來看，定點制下進行小數(shù)的運算實質(zhì)上還是按處理器字長進行整數(shù)的各種運算。同樣，對數(shù)進行存儲和處理的時候，數(shù)的各種碼制表示規(guī)律也同樣適用于定點表示下的小數(shù)處理。計算機基礎(chǔ)課程中，計算機在表達數(shù)值的時候，二進制數(shù)的最高位一般用作符號位，0表示正數(shù)，1表示負數(shù)，一個數(shù)可以有原碼、補碼和反碼三種形式。對于正數(shù)來說，三種碼都一樣，而對于負數(shù)，這三種碼并不一致，用途也各不相同。如原碼適合做乘除法，常被用作設(shè)計串行乘法器，而補碼適合做加減法，加法器的硬件多采用補碼制。表7.2列出了三種碼的特性對比。在數(shù)字信號處理器（DSP）中，通常采用補碼制。數(shù)的轉(zhuǎn)換在MATLAB中的實現(xiàn)有如下幾種方式。

(1)函數(shù)dec2bin可以實現(xiàn)將一個十進制正整數(shù)轉(zhuǎn)換成一個二進制的字符串。調(diào)用格式：

dec2bin(D)

dec2bin(D,N)輸入變量：D是小于252的非負整數(shù)，N為轉(zhuǎn)換后的二進制字符串的長度。

(2)函數(shù)bitcmp可以用于求一個十進制正整數(shù)的補數(shù)，和函數(shù)dec2bin一起使用可以求一個十進制整數(shù)的反碼和補碼。調(diào)用格式：

bitcmp(A,N)輸入變量：A是無符號型整數(shù)，其中A≤2N－1，計算結(jié)果為2N－1－A。

(3)函數(shù)num2bin用于將一個數(shù)值矩陣轉(zhuǎn)化成二進制的字符串。調(diào)用格式：

B=num2bin(Q,X)輸入變量：X是數(shù)值矩陣，Q用于表明X的屬性。在MATLAB里用函數(shù)?quantizer生成量化目標Q，常用的調(diào)用格式為Q=quantizer(［w,f］)，w是字符串B的長度，f是小數(shù)位數(shù)。需要注意的是，如果不是0.xxxx，必須要給整數(shù)位保留兩個比特（包含一位符號位）。

【例7-1】用MATLAB編程求十進制數(shù)-1325的二進制原碼、反碼和補碼表示。

解：因為210<1325<211，需要用11位的二進制數(shù)來表示該十進制數(shù)，假設(shè)用12位二進制數(shù)表示，其中第12位為符號位。利用前面介紹的MATLAB函數(shù)編寫程序如下：

x=-1325;

a=abs(x);

b=dec2bin(a+211,12)%原碼

c=dec2bin(bitcmp(a,12),12)%反碼

d=dec2bin(bitcmp(a,12)+1,12)%補碼運行結(jié)果：

b=110100101101

c=101011010010

d=101011010011■

【例7-2】

用MATLAB編程求十進制數(shù)123.874的二進制原碼表示，小數(shù)位數(shù)保留8位。

解：利用前面介紹的MATLAB函數(shù)編寫程序如下：

Q=quantizer(［16,8］);

num2bin(Q,123.874)運行結(jié)果：

ans=

0111101111011111■

4.定點量化誤差分析定點制下，當Q值確定以后，即用來表示小數(shù)的寄存器位數(shù)L確定后，其可表示的最小數(shù)的單位也就確定了，記為2－L，這個值稱為量化間距，記作q。超出L位的部分，可以通過直接截斷的方式進行處理，所產(chǎn)生的誤差稱為截尾誤差，也可以通過舍入的方式進行處理，產(chǎn)生舍入誤差。如果數(shù)P的小數(shù)部分是x，通過M位二進制表示，存入Q值定義為L的寄存器中被量化為Q［x］，則其量化誤差e定義為(7-2)e=Q［x］－xx∈［0,1)

1）截尾誤差對于正數(shù)，原碼、補碼和反碼的形式都相同，有(7-3)(7-4)(7-5)顯然，此時有－(2－L－2－M)≤e≤0，當M→∞時，有－2－L≤e≤0。也就是說，截尾誤差最大不超出量化間距q。對于負數(shù)，由于三種碼的表達方式不同，誤差也不同。(7-6)(7-7)(7-8)易知，0≤e≤2－L－2－M。當M→∞時，有0≤e≤2－L。也就是說，用原碼表示負數(shù)時，截尾誤差始終為正，誤差均值為2－(L+1)，且最大誤差不超出量化間距q。當負數(shù)用補碼表示時，有(7-9)(7-10)(7-11)

易知，－(2－L－2－M)≤e≤0。當M→∞時，有－2－L≤e≤0。也就是說，用補碼表示負數(shù)時，其截尾誤差與正數(shù)情況一致，始終為負，均值為－2－(L+1)，最大誤差不超出量化間距q。當負數(shù)用反碼表示時，有(7-12)(7-14)(7-13)

2）舍入誤差當定點數(shù)用于表示小數(shù)的寄存器長度為L，將長度為M(M>L)的數(shù)據(jù)存入該寄存器進行舍入處理時，就是將第L+1位加1，然后再截斷數(shù)據(jù)到第L位。容易理解，無論采用原碼、補碼還是反碼表示，這個過程都調(diào)整了誤差范圍，誤差均值的絕對值由2－(L+1)調(diào)整為0，而誤差范圍調(diào)整為［－2－(L+1)，2－(L+1)］。7.1.3浮點制誤差分析從前面對定點制表示分析的過程可以看出，定點制的優(yōu)點是運算簡便，對處理器要求低。但數(shù)的表達能力有限，所處理數(shù)的動態(tài)范圍較小。同時，寄存器的利用效率低，如表示較小的小數(shù)時，小數(shù)的有效位數(shù)較短，由截尾舍入產(chǎn)生的百分比誤差隨著數(shù)的絕對值的減小而增加?？茖W計數(shù)法中，允許將任意一個數(shù)字表示為一個純小數(shù)與一個指數(shù)相乘的形式

(P)R=S·RcS一般為絕對值介于0~1之間的規(guī)格化尾數(shù)，機器中可用原碼或補碼表示。R為基數(shù)，在數(shù)字處理器中通常取2。c為浮點數(shù)的階碼，也即是指數(shù)，為整數(shù)。采用浮點制表示，可以在兼顧表達范圍的同時保證運算精度。根據(jù)IEEE754國際標準，常用的浮點數(shù)有單精度浮點數(shù)(single)和雙精度浮點數(shù)(double)兩種格式。單精度浮點數(shù)占用4個字節(jié)(32位)存儲空間，在數(shù)字寄存器中，單精度浮點數(shù)的存儲格式如圖7-1所示，包括1位符號位S，8位階碼(指數(shù))E，23位尾數(shù)F。其表示范圍為

(1)當S=0，E=127，F(xiàn)的23位均為1時，表示的浮點數(shù)為最大的正數(shù)：

(01.111..1)2×2127=(2－2－23)×2127≈3.4×1038

(2)當S=1，E=127，F(xiàn)的23位均為0時，表示的浮點數(shù)為絕對值最大的負數(shù)：

(10.000..0)2×2－127=－2×2－127≈－3.4×1038

(3)當S=0，E=－127，F(xiàn)的23位均為0時，表示的浮點數(shù)為最小的正數(shù)：

(01.000..0)2×2－127=1×2－127≈5.9×10－39

(4)當S=1，E=－127，F(xiàn)的23位均為1時，表示的浮點數(shù)為絕對值最小的負數(shù)：

(10.111..1)2×2－127=(－1－2－23)×2－127≈－5.9×10－39雙精度浮點數(shù)占用8個字節(jié)(64位)存儲空間，包括1位符號位、11位階碼、52位尾數(shù)，數(shù)值范圍為1.7E-308～1.7E+308。圖7-1單精度浮點數(shù)存儲格式浮點數(shù)進行加減運算一般需要有五個步驟：

(1)對階：使兩數(shù)的小數(shù)點位置對齊。

(2)尾數(shù)求和(差)：將對階后的兩尾數(shù)按定點加減運算規(guī)則求和(差)。

(3)規(guī)格化:為增加有效數(shù)字的位數(shù)，提高運算精度，必須將求和(差)后的尾數(shù)規(guī)格化。

(4)舍入：為提高精度，要考慮尾數(shù)右移時丟失的數(shù)值位。

(5)判斷結(jié)果：即判斷結(jié)果是否溢出。進行浮點數(shù)的乘除運算分為階碼運算和尾數(shù)運算兩個步驟。階碼運算較為簡便，若是乘法則對階碼寄存器段做加法運算，除法則做減法運算。對于浮點數(shù)的乘法，尾數(shù)運算一般需要四個步驟：

(1)預(yù)處理：檢測兩個尾數(shù)中是否有一個為0，若有一個為0，乘積必為0，不再作其他操作；如果兩尾數(shù)均不為0，則可進行乘法運算。

(2)相乘：兩個浮點數(shù)的尾數(shù)相乘可以采用定點小數(shù)的任何一種乘法運算來完成。

(3)規(guī)格化：相乘結(jié)果可能要進行左規(guī)，左規(guī)時調(diào)整階碼后如果發(fā)生階下溢，則作機器零處理；如果發(fā)生階上溢，則作溢出處理。

(4)尾數(shù)截斷：尾數(shù)相乘會得到一個雙倍字長的結(jié)果，若限定只取1倍字長，則乘積的若干低位將會丟失?？梢酝ㄟ^前面討論的定點制截斷處理方法進行截尾或舍入處理。對于浮點除法，尾數(shù)運算需要先檢測被除數(shù)是否為0，若為0，則商為0；再檢測除數(shù)是否為0，若為0，則商為無窮大，另作處理。若兩數(shù)均不為0，則可進行除法運算。兩浮點數(shù)尾數(shù)相除同樣可采取定點小數(shù)的任何一種除法運算來完成。對已規(guī)格化的尾數(shù)，為了防止除法運算結(jié)果溢出，可先比較被除數(shù)和除數(shù)的絕對值，如果被除數(shù)的絕對值大于除數(shù)的絕對值，則先將被除數(shù)右移一位，其階碼加1，再作尾數(shù)相除。此時所得結(jié)果必然是規(guī)格化的定點小數(shù)。與定點制相比，浮點制一定程度上可以兼顧動態(tài)范圍和運算精度，但其指數(shù)和尾數(shù)都需要參與運算，運算過程復(fù)雜，實現(xiàn)難度大，硬件成本高，在僅提供定點運算的處理器中很難獲得實時的運算結(jié)果。從處理誤差方面看，浮點制所產(chǎn)生的誤差傳播范圍廣，分析起來難度較大，此處不做詳細論述，可參閱相關(guān)書籍。7.1.4成組浮點制(BFP)成組浮點制(BlockFloatingPoint)也稱為塊浮點制，基本思想是兼顧定點制和浮點制的優(yōu)點，將一組數(shù)值相近的數(shù)據(jù)定義成一個具有統(tǒng)一指數(shù)的數(shù)據(jù)塊，換句話說，就是將該組數(shù)據(jù)同時根據(jù)這個共享指數(shù)進行縮放，在保證動態(tài)范圍和精度的同時又不需要考慮彼此間指數(shù)的影響。

BFP是一種在工程實現(xiàn)中較常采用的方法。如Altra公司提供的IP核FFTMegacore中就集成了該算法。

7.2A/D變換的有限字長效應(yīng)7.2.1A/D變換及其量化誤差的統(tǒng)計分析

A/D變換包括采樣和量化兩部分，如圖7-2所示。模擬信號xa(t)經(jīng)過采樣后轉(zhuǎn)換為時域離散信號x(n)，然后對x(n)做截尾或者舍入的量化處理后得到二進制數(shù)字信號x(n)。由于A/D變換總是采用定點制，因此需要將模擬信號乘以比例因子A，以限定其最大值不能超過A/D變換的動態(tài)范圍，即

x(n)=Axa(t))|t=nT=Axa(nT)

(7-15)顯然信號x(n)具有無限精度，受存儲單元的字長限制，必須對其做截尾或者舍入的量化處理，用e(n)表示量化誤差，則^(7-16)圖7-2A/D變換器的非線性模型

A/D變換的量化方式和數(shù)的表示方式直接決定了其量化特性。設(shè)A/D變換的輸出是字長為b+1位的補碼定點小數(shù)，其中b位是小數(shù)部分，若用eT(n)表示采用截尾的誤差，用eR(n)表示采用舍入處理的誤差，由7.1節(jié)可以得到量化誤差的范圍為(7-17)(7-18)這里，q=2－b表示量化階距。式(7-17)、(7-18)給出了量化誤差的范圍，但是要想精確地刻畫出每個量化誤差的值還是非常困難的。通常用統(tǒng)計分析的方法來分析量化誤差的統(tǒng)計特性，研究A/D變換的有限字長效應(yīng)。A/D變換器的統(tǒng)計模型如圖7-3所示。圖7-3A/D變換器的統(tǒng)計模型為了研究其統(tǒng)計特性，首先對量化誤差e(n)作如下假設(shè)：

(1)e(n)是平穩(wěn)隨機序列；

(2)e(n)與采樣信號x(n)互不相關(guān)；

(3)e(n)序列中任意兩個值之間不相關(guān)，即e(n)為白噪聲序列；

(4)e(n)在其取值范圍內(nèi)均勻等概分布。根據(jù)上述假設(shè)，e(n)是與輸入信號完全不相關(guān)的、均勻分布的白噪聲序列，采用截尾和舍入時誤差的概率密度函數(shù)分別如圖7-4(a)、(b)所示。圖7-4量化誤差的概率分布(a)截尾誤差；(b)舍入誤差采用截尾處理時的均值和方差分別為(7-19)(7-20)采用舍入處理時的均值和方差分別為(7-21)(7-22)分析式(7-19)~(7-22)可知：

(1)采用截尾處理時誤差序列的均值不為零，也就是說誤差序列eT(n)中包含直流分量，直流分量的存在會使信號的頻譜在頻率等于0處存在δ函數(shù)，從而影響信號的頻譜結(jié)構(gòu)。而采用舍入處理時誤差序列的均值為0，不存在直流分量。因此，實際應(yīng)用中多采用舍入處理的方式。

(2)采用截尾處理和舍入處理時誤差序列的方差相等，且都為σ2T=σ2R=σ2e，都取決于A/D變換的字長b，字長越長，量化誤差越小。用符號σ2x表示信號功率，則量化信噪比為(7-23)對數(shù)表示為(7-24)從式(7-24)可以看出，信號功率一定的情況下，字長每增加1bit，量化信噪比約增加6個分貝。

【例7-3】

假設(shè)語音信號量化編碼時，選用12bit的A/D，其動態(tài)范圍為0~5V，求系統(tǒng)量化誤差的均方差。

解：量化階距電壓：

Vq=5×2－b=5×2－12=1.2mV

7.2.2量化噪聲通過線性系統(tǒng)這一節(jié)在不考慮系統(tǒng)實現(xiàn)誤差和運算誤差的情況下，將系統(tǒng)近似看作是完全理想的，即具有無限精度的線性系統(tǒng)，討論量化信號通過線性時不變系統(tǒng)的問題。假設(shè)量化序列x(n)=x(n)+e(n)，線性時不變系統(tǒng)的單位脈沖響應(yīng)為h(n)，由7.2.1節(jié)的假設(shè)，信號x(n)和量化噪聲e(n)相互獨立，假設(shè)系統(tǒng)為因果系統(tǒng)，則根據(jù)線性疊加原理，系統(tǒng)的輸出為

(7-25)^其中，y(n)是系統(tǒng)對信號x(n)的響應(yīng)，f(n)是噪聲信號e(n)通過線性系統(tǒng)的輸出，有(7-26)量化噪聲通過線性系統(tǒng)的框圖如圖7-5所示。圖7-5量化噪聲通過線性系統(tǒng)(1)若e(n)是舍入誤差，則輸出噪聲信號f(n)的均值為(7-27)由式(7-21)me=0,可求得mf=0。其方差為由前面的假設(shè)e(n)為白噪聲序列，序列中任意兩個值之間不相關(guān)。因此

E[e(n－m)e(n－k)]=σ2eδ(m－k)則有(7-28)從上式可以看出，量化噪聲通過線性時不變系統(tǒng)后，其輸出信號的方差依然和量化字長成反比。在量化字長一定的情況下，其輸出信號的方差取決于系統(tǒng)單位脈沖響應(yīng)的能量。假設(shè)h(n)是實序列，由帕塞伐定理(7-29)代入式(7-28)可得(7-30)

(2)若e(n)是截尾誤差，輸出噪聲信號f(n)的方差仍為式(7-28)，其均值為(7-31)顯然，輸入信號時截尾量化，在輸出端也引入了一個直流分量?！纠?-4】

已知IIR濾波器的系統(tǒng)函數(shù)為假設(shè)其輸入信號x(n)為8位A/D變換器(b=7)的輸出，求濾波器輸出端的量化噪聲功率。

解：由于A/D變換的量化效應(yīng)，濾波器輸入端的噪聲功率為^(7-32)濾波器的輸出噪聲功率為圍線c內(nèi)只有兩個極點z1=0.7，z2=0.8，根據(jù)留數(shù)定理有

7.3數(shù)字濾波器系數(shù)量化的有限字長效應(yīng)7.3.1系數(shù)量化對IIR濾波器性能的影響無限長單位脈沖響應(yīng)(IIR)濾波器的系統(tǒng)函數(shù)可以表示為(7-33)對分子和分母多項式的系數(shù)進行量化。假設(shè)系數(shù)bi和ai的量化值分別為b

i和ai，對應(yīng)的量化誤差為Δbi和Δai，則^^(7-3４)量化后的系統(tǒng)函數(shù)為通過前面章節(jié)的學習可知，系統(tǒng)性能在很大程度上取決于系統(tǒng)的極點。而量化誤差的存在會造成系統(tǒng)極點位置的改變，從而影響系統(tǒng)的穩(wěn)定性。為了衡量系數(shù)量化對極點位置的影響，定義系統(tǒng)中每個極點位置對各系數(shù)偏差的敏感程度為極點位置靈敏度，用極點位置靈敏度來反映系數(shù)量化對濾波器穩(wěn)定性的影響。假設(shè)量化前系統(tǒng)的理想極點為：zi(i=1，2，…，N)，對分母多項式A(z)進行因式分解可得(7-36)系數(shù)量化后，由于量化誤差Δai的存在，對應(yīng)的極點變?yōu)閦

i=zi+Δzi(i=1，2，…，N)。因為極點zi的值和分母多項式中每個系數(shù)ai(i=1，2，…，N)都有關(guān)系，因此可以把極點zi表示成如式(7-37)所示的函數(shù)形式:

zi=zi(a1，a2，…，aN)i=1，…，N

(7-37)由系數(shù)量化帶來的極點的改變量Δzi為(7-38)7.3.2系數(shù)量化對FIR數(shù)字濾波器性能的影響上一節(jié)我們通過分析系數(shù)量化帶來的系統(tǒng)極點位置的變化，研究了系數(shù)量化對IIR數(shù)字濾波器穩(wěn)定性的影響。這一節(jié)我們將以直接型結(jié)構(gòu)為例分析系數(shù)量化對FIR系統(tǒng)性能的影響。對于直接型FIR濾波器而言，系數(shù)量化時濾波器仍然能夠保持線性相位特性，系數(shù)量化主要會對系統(tǒng)的幅頻特性產(chǎn)生影響。由4.3節(jié)可知，F(xiàn)IR濾波器的系統(tǒng)函數(shù)可以表示為(7-44)對系數(shù)h(n)進行量化，假設(shè)量化誤差為e(n)，則量化后的結(jié)果h(n)為

(7-45)量化以后的系統(tǒng)函數(shù)為^(7-46)式中E(z)是誤差函數(shù)e(n)的z變換，對應(yīng)的頻率響應(yīng)為(7-47)量化字長一定時，系數(shù)的量化誤差是個有界量。假設(shè)量化字長為b，量化階距為q，舍入時|e(n)|≤q/2=2－b/2。則量化誤差的幅度響應(yīng)為(7-48)從式(7-48)可以看出，濾波器階數(shù)固定的情況下，量化誤差的幅度響應(yīng)和量化字長成反比，可以通過增加量化字長來減少量化誤差的影響。從概率統(tǒng)計的角度，由7.2.1節(jié)可知，采用有限精度的舍入運算時，e(n)的均值為0，方差等于q2/12。假設(shè)N個系數(shù)的量化誤差之間是不相關(guān)的，則量化誤差的頻率響應(yīng)的方差σ2E是這N項誤差的方差之和，即(7-49)

7.4定點運算對數(shù)字濾波器的影響在工程實現(xiàn)中，為了滿足系統(tǒng)信噪比要求，需要分析濾波器的運算誤差，以選擇合適的濾波器運算位數(shù)。數(shù)字濾波器的實現(xiàn)需要乘法、加法和延遲三種基礎(chǔ)運算，其中延遲運算不會造成字長的變化。定點制運算中，加法運算不會造成尾數(shù)字長的變化，誤差主要來自于數(shù)的溢出；而兩個尾數(shù)長度為n的數(shù)相乘運算后，乘積的尾數(shù)長度為2n，受存儲器長度的限制，要對乘積做截尾或者舍入運算，這勢必會引入誤差，影響濾波器的性能。因此，這里僅分析乘法和加法，特別是乘法運算的舍入誤差對數(shù)字濾波器性能的影響。圖7-6定點相乘運算的算法流圖(a)理想相乘；(b)實際相乘的非線性流圖；(c)統(tǒng)計分析的流圖7.4.1IIR濾波器運算中的有限字長效應(yīng)一個N階的IIR濾波器的輸入輸出關(guān)系可以用式(7-51)所示的N階的線性常系數(shù)差分方程來描述:(7-51)(7-52)對式(7-52)兩邊作Z變換可得根據(jù)式(7-33)，有(7-53)從式(7-53)可以看出舍入處理后，濾波器的輸出y(n)等于y(n)和f(n)之和，其中f(n)是系統(tǒng)總誤差e(n)通過傳輸函數(shù)為He(z)=1/A(z)的系統(tǒng)的輸出。^根據(jù)前面的假設(shè)，每一個乘法支路中的舍入誤差是均值為零的平穩(wěn)白噪聲序列，系統(tǒng)總誤差e(n)是M+N+1個乘法支路的舍入誤差之和，因此e(n)也是均值為零的平穩(wěn)白噪聲序列，其方差(功率)等于各個乘法支路中舍入誤差的方差之和。假設(shè)e(n)的方差為σ2e，則(7-54)式中，σ2ε為單個乘法支路中舍入誤差的方差。由7.2.2節(jié)，e(n)通過系統(tǒng)He(z)后輸出信號f(n)的方差σ2f為(7-55)式(7-55)也可表示為(7-56)從式(7-55)和式(7-56)可以看出，量化字長(量化階距)固定的情況下，σ2f取決于IIR濾波器的系統(tǒng)函數(shù)。顯然，濾波器的結(jié)構(gòu)不同時，σ2f的值也不相等。

【例7-5】一個因果的IIR系統(tǒng)的系統(tǒng)函數(shù)試分別在直接型、級聯(lián)型和并聯(lián)型情況下，求解由乘法舍入誤差所產(chǎn)生的輸出噪聲的方差。

解：

1)直接型結(jié)構(gòu)由有A(z)=1－0.9z－1+0.2z－2=(1－0.5z－1)(1－0.4z－1)M=0，N=2直接型結(jié)構(gòu)