第7部分数字信号理中的有效字长效应

1、第7章数字信号处理中的有效字长效应 前言w数字信号处理的实质：一组数值运算。w从设计的角度来讨论：认为数字是无限精度的。w从实现的角度考虑：数字的精度就是有限的n数字系统中的每一个数总是用有限字长的二进制数码表示，运算过程中需要的数字信号的值、系统的系数和运算过程中的结果都是存储在有限字长的存储单元中的，此时数字的精度就是有限的。w从设计时的无限精度到实现时的有限精度，会产生相对于原设计系统的误差，严重时会导致系统崩溃。w在数字系统中三种因有限字长的影响而引起误差的因素na/d变换的量化效应n系数的量化效应n数字运算过程中的有限字长效应7.1 二进制数的表示与量化误差 1、定点数的表示 n所有

2、数据小数点的位置是固定不变的 2、浮点数的表示 n浮点数的小数点位置是不固定的，它随每个数的大小而变化 prmnm是数的尾数部分，r是数的基数，p是阶 基数确定后，浮点数就完全由尾数和阶决定7.1.1 二进制数的表示一个浮点数可以表示如下： ：阶p的符号位； ：阶码，阶p的绝对值部分 ：尾数m或浮点数n的符号位，尾符或数符； ：尾数的绝对值部分，尾码。 psnpmsrmpsnpmsrm3、ieee浮点数 w规格化浮点数的小数点在数符sm的后面，且小数点前有一个隐含“1”。其阶码连同阶符统一编码，浮点数的基数为2。wieee浮点数格式分单精度和双精度两种，单精度数为32位，双精度数是64位。 单

3、精度浮点数格式： 双精度浮点数格式： 数符阶码（8位）尾数（23位）b31b22b23b0b30数符阶码（8位）尾数（23位）b63b51b52b0b62带符号数的表示 设任意数x的(b十1)位码的形式为 a0a1 a2 a3 ab a0表示符号位 1、原码： a0是符号位， a1 ab是小数的绝对值 01)1(110 xxxxxx或原biiiaax12) 1(0表达式（a0a1 a2 a3 ab）表示的十进制数x为 (2)反码 正数的反码表示与其原码相同，即x反x（x0）负数的反码由其绝对值按位求反后得到 x反2-x原 表达式（a0a1 a2 a3 ab）表示的十进制数x为 biiibaax

4、10221(3)补码 正数的补码表示与其原码相同，即x补x（x0） 负数的补码由它的绝对值求反加1后得到： 01110 xxxxx反补biiibiiiaaax110212表达式（a0a1 a2 a3 ab）表示的十进制数x为 7.1.2 量化误差 一、定点运算中的截尾误差和舍入误差 1、截尾误差 对于正小数x0：原码、反码、补码的表示法相同，因而量化影响也相同。 截尾前x有b1位，有 112biiiax截尾后x有b位，记做 ，有 xqt biiitaxq12以et表示截尾误差，则有 11111222bbiiibiiibiiittaaaxxqe 正小数截尾后数值变小，故截尾误差总是负的。当被截位

5、ai (ib1到ib1)均为l时，为最大截尾误差 11111222bbiiibiiibiiittaaaxxqe111max222bbbbiite0221tbbe有021tbeqq00 xeqt或令 量化间距或量化步阶，因此下式成立 bq 2对于负小数x0，由于a00，不同码制x的表示法各不相同，因而产生的量化误差也不相同。 i）原码:截尾后负小数的绝对值变小，截尾误差为正： 11110221biiibiiiaaax biiitaxq12 112bbiiittaxxqe所以1220bbte得00 xqet满足ii）对于补码， 11110212biiibiiiaaax biiitaxq121 11

6、2bbiiittaxxqe所以00 xeqt有iii）对于反码： 111110221221biiibbiiibaaax biiibtaxq1221 111222bbiiibbttaxxqe所以00 xqet有w结论：w原码与反码的截尾误差与数的正负有关：n正数时误差为负，当x0时，qet0n负数时误差为正，当x0时，0et qw补码的截尾误差皆为负数n对所有的x，qet0。 定点制截尾处理的量化特性（q2b） qqxqtxqqxqtx(a)补码 (b)原码、反码2、舍入误差 xqr舍入误差总是处在q/2之间，用 表示对x作舍入处理，er表示舍入误差，有： xxqerr22qeqr定点制舍入处理

7、的量化特性 qqxqrx定点运算中的截尾和舍入误差（q2b） 截尾误差舍入误差正数q et0q /2erq /2负数原码0et q反码0et q补码q 0时，a00，三种码制的截尾误差均为qet0 得，取由于,222/21mincccxx002xxqrc002xqr2、浮点截尾 当x0时，a01，原码和反码的截尾误差为0etq 得，取由于,22/221mincccxx020 xqxcr002xqr对于补码，截尾误差为qet0，相应的，可得 020 xqr浮点运算中的相对误差（q2b） qqr02rqqqrqr20 截尾舍入正数qet0负数原码反码补码02rq7.2 a/d转换器中的量化误差 w

8、a/d转换器n将输入的模拟信号xa(t)转换为b位二进制数字信号的器件。nb的数值可以是8，12或高至20。wa/d变换器前一般都加一个前置模拟低通滤波器n滤除高于折叠频率(抽样频率之半)的频率，w模拟输入信号必须乘一个比例因子na/d变换器总是定点制的，必须使信号不超过a/d变换的动态范围a/d转换器包括w抽样n产生抽样序列x(n)xa(t)|t=ntxa(nt)，nx(n)： 可看成是一个无限精度的数字信号n在满足抽样定理的前提下，模拟信号时间离散化的过程是可逆的。w量化n对抽样序列进行幅度上的离散化之后，用某种格式的数字代码来表示。n量化过程是不可逆的n必定要引入量化误差或量化噪声。n量

9、化噪声的大小决定了a/d转换器的动态范围，是恒量a/d转换器性能的一个最重要指标。 7.2.1 量化误差的统计分析 设量化器的输入信号为随机序列x(n)，其量化误差也为随机序列e(n)：e(n)=qx(n)-x(n) 假设e(n)具有下列特性：1. e(n)是一个平稳随机序列2. e(n)与信号也不相关3. e(n)本身的任意两个值之间不相关，具有白噪声性质4. e(n)在其误差范围内均匀等概分布。 a/d转换原理图 采样量化编码)(txa)(ntxa)( nx理想a/dc)(txa nx +)(nx)(nea/d转换的统计模型 e(n)的两个最重要的统计参数 n均值me ：代表噪声的直流分量

10、 e2：代表了除去直流分量后量化噪声的平均功率 me e nep e dee ( )( )eeee e nmemp e de222( ( ) ()( )其中 :e表示取数字期望 p(e)是误差值e(n)的概率密度 三种误差范围的概率密度函数 1、舍入处理 w舍入误差：q /2e(n)q /2we(n)的概率分布密度为1/q 可得均值me 及方差e20211)(22222qqqqreeqedeqdeeepm222322222121311)(qeqdeeqdeepmeqqqqree2、截尾处理 w对于x0的三种码制和x0的补码 ：ne(n)的误差为q e(n)0n概率分布密度为1/q 可得均值me

11、 及方差e2qeqedeqdeeepmqqte21211)(0202022212121)(qdeqeqdeepmeqteew对于x0的原码和反码：ne(n)的误差为0et qn概率分布密度为1/q 可得均值me 及方差e2qeqedeqdeeepmqqte21211)(0202022212121)(qdeqeqdeepmeqteew结论：各种情况的方差 均为 ，e22121q220qq、不同的只是均值me，分别为 。 另：由于截尾噪声具有直流分量，将影响信号的频谱结构，因此一般采用舍入处理。量化噪声的方差 与量化间隔的平方 成正比， ：量化间隔，信号处理时选用的字长信号处理时选用的字长b1越长

12、，量化噪声的方差越小。越长，量化噪声的方差越小。e22qbq 27.2.2 量化信噪比与所需字长的关系 w量化的信噪比 信号的平均功率 与量化噪声的平均功率 之比 2xe222222221212xbxexq信噪比用对数表示时记作snr，单位db )lg(1079.1002. 6lg10)(222xexbdbsnr信号功率 越大，信噪比越高；随着字长b的增加，信噪比也增大，字长b每增加一位，则信噪比增加约6db。字长越长，a/d变换的信噪比越高 2x 为了使信号不超过定点制运算所允许的动态范围，用一个小于1的正数a（0a1去乘x(n) 信号的平均功率 与量化噪声的平均功率 之比 2xae2aba

13、dbsnrxexlg20)lg(1079.1002. 6lg10)(2222选择a值：使 12xa25. 102. 6)lg(2079.1002. 6lg10)(222babadbsnrxexw若需得到信噪比大于70db，至少需要满足b12bitw字长过长也无必要：n输入信号xa(t)本身也有一定的信噪比，字长长到a/d变换器的量化噪声比xa(t)的噪声电平更低则没有意义w为提高信噪比，n可以增大输入信号，但这受到a/d变换器动态范围的限制，n还可以增加字长b，但这又受到输入信号xa(t)信噪比的限制 7.2.3 量化噪声通过线性非时变系统 量化噪声通过线性系统 +)()(zhnh)(nx)(

14、 nx)()()( nenynyf)(ne)()()()()()()( )( nenynhnenxnhnxnyf设e(n)是定点补码舍入误差，e(n)的均值为me 、方差为 e2则系统量化噪声的输出 的均值mf和方差 计算如下： 2f)(nef0)()()()(0meffmhmnhneeneem022222)()()(meffffmhneemnee根据parseval定理， 也可以用下式表示： 2f或者在单位圆上计算 cemefzdzzhzhjmh)()(2)(120222dehdehehjejjef2222| )(|2)()(2如果e(n)是补码截尾白噪声：输出噪声的方差 不变，输出的均值m

15、f如下： 2f)()()()()()()(000jememffehmmhmmnemhenhneeneem7.3 系数量化对数字滤波器的影响 w7.3.1 极点位置灵敏度 n每个极点位置对各系数偏差的敏感程度n决定系统的稳定性n对系统性能的影响亦较大 理想数字滤波器的系统函数为 )()(1)(10zazbzazbzhniiimiiiniiimiiizazbzh101)( 其中 是系数ai、bi的量化结果：iiba、对系数ai、bi量化后，其实际传递函数为：，iiiaaaiiibbb下面讨论系数量化误差对极点的影响：nkkniiizzzaza111)1 (1)(n21，kzzkz n21，kzzz

16、kk系数量化后的极点为：原系统函数h(z)的分母多项式为令a(z)0，得到h(z)的极点：极点位置的偏差量，由各个系数偏差 引起的，因此kzian211，kaazziniikk结论： 就是极点 对系数 变化的灵敏度v 越大， 对 的影响也越大v 越小， 对 的影响就越小iakzkzikaziakzikaziaikaz下面根据a(z)来求这个极点位置灵敏度 的表达式： ikazkkzziikzzkazaazzza)()(kzzkiikzzaazaaz)()(iizaza)(nklllnnklllkzzzzzzzza1111)()1 ()(而：得极点位置灵敏度为：nklllkinkikzzzaz1

17、)(分母中的每一个因子 代表着某一极点 指向当前极点 的矢量，而整个分母正是所有极点指向当前极点 的矢量积 )(lkzz lzkzkz高阶直接型结构滤波器对系数量化误差的敏感性高于低阶直接型结构滤波器并联型结构及级联型结构每对极点只受与之有关的两个系数的影响，且每个子系统的极点密集度要稀疏得多，因而极点位置受系数量化的影响比直接型结构要小得多 rezimzj1-10rezimzj1-10a 极点间距离长b 极点间距离短极点位置灵敏度与极点间距离成反比 例1 设一低通滤波器的传递函数如下，分析计算系数量化对极点位置的影响。332211321119508. 08934. 29425. 211)(z

18、azazazzzzh解 经计算求得 的极点分别是)(zh0052352198. 098. 099. 0jjezzezz，）（）（151510098. 0198. 0199. 011)(zezezzhjj于是：研究当 时，仅仅由一个系数 的量化所引起的极点 的变化 031 aa2a1z1z）（00552312121198. 099. 098. 099. 099. 0)(jjeeazzzzazz若字长为8位，由量化引起的误差q/2可达大约0.002，求得 ，于是2578. 01z248. 12578. 099. 0111zzz极点远远超出单位圆，这样的变化显然是太大了。再研究 ，即将极点移到单位圆

19、上需要的字长 01. 01z00007578. 02计算得： a510578. 7二进制数中为2-14 2-13，字长至少要14位。510578. 72a510578. 7如果 2.945的量化误差等于 ，就会使量化后系统的极点移到单位圆上可见一个三阶系统对字长的要求已经非常严格了，如果阶数再高对量化误差的要求将更加苛刻。w考察用三个一阶的环节级联或并联组成这个系统。n每一个环节中极点从0.99变到1.00，允许变化0.01n而且这个环节的稳定性并不受另一环节影响n所以其所需字长为7位w结论：n系数量化对零、极点位置的影响与零、极点位置的分布以及滤波器的结构均有密切的关系n高阶滤波器：避免采用

20、直接型的结构，而应分解为最低阶的级联结构或并联结构n对于极点灵敏度很高的场合，可以来用双精度的系数，以便有效的达到精度的要求。7.3.2 系数量化对二阶子系统极点位置的影响 w高阶系统n级联型和并联型优于直接型。w级联型和并联型的基本子系统是二阶节n一个具有共轭极点的二阶系统有各种不同的结构n不同结构对于系数量化的敏感度也不同 二阶iir系统的系统函数为 221111)(zazazhjrez， 21221112211cos21)1)(1 (1zrzrzrezrezazajj2cos122arar，具有一对共轭对称的复极点有得到共轭对称极点对组成的基本二阶网络的直接型实现)(zx)(zy1z1z

21、cos2r2r221111)(zazazh2cos122arar，2cos122arar，说明:对于二阶网络，其极点的半径 完全由系数 决定，极点在实轴上的坐标值 取决于系数 。2acosr1ar1a2a如果 、 用三位字长表示，b3(不算符号位)，r87cosr只能表达8种半径 值和 之间的15种实轴坐标值三位二进码1,2,3 所表达的的值 |a1 |极点横坐标 a2三位二进码 所表达的a2的值 极点半径0.000.000.0000.0000.0000.0000.010.250.1250.0010.1250.3540.100.500.2500.0100.2500.5000.110.750.3

22、750.0110.3750.6121.001.000.5000.1000.5000.7071.011.250.6250.1010.6250.7891.101.500.7500.1100.7500.8651.001.750.8750.1110.8750.9352cos1ar321、2arw极点位置如下图中的网眼节点n如果所需要的极点位置不在这些网眼节点上时，就只能以最靠近的一个节点来代替这一极点位置，这样就引入了极点位置误差rexjimx10-1a1=1.751.501.251.000.750.500.25a1=0.250.500.751.001.251.501.75a2=0.1250.2500

23、.3750.5000.6250.7500.875w系数量化使零极点位置的取值范围由一个连续域变为一个离散的平面点阵，从而造成零极点的漂移，导致系统特性的改变。w在平面上量化位置的分布密度是不均匀的n在实轴附近分布得稀，在虚轴附近分布得密；n半径小的地方分布得稀，半径大的地方分布得密。w这祥就会使实轴附近的极点(例加低通、高通滤波器)量化误差较大、而对虚轴附近的极点(例如带通滤波器)量化误差较小w这种分布只是二阶直接型结构的情况，不同结构的滤波器，系数对零极点位置的影响是不一样的 结论：)(nx)(ny1z1zcosrcosrsinrsinr该网络的系统函数为 ：221)(cos211)(zrz

24、rzh当系数量化时，是对 及 进行量化，因而所得到的网格点子在z平面是均匀分布的cosrsinr这里系数量化后对z平面的所有区域，所产生的误差是相同的。基本二阶网络的另一种实现具有共轭极点对的二阶数字网络耦合形式实现情况下，系数量化为三比特时极点的可能位置 7.3.3 频率响应偏差的统计分析 1、iir数字滤波器系数量化的统计分析 n阶iir直接型结构为例，其理想精度的系统函数为 )()(1)(10zazbzazbzhniiimiii其中 是系数ai、bi的量化结果：iiba、对系数ai、bi量化后，其实际传递函数为：，iiiaaaiiibbb故系数量化后，实际的系统函数为niiimiiiza

25、zbzh101)( )()()()()()()()()()()()()()()()()()()()(11100zhzhzazazhzazbzhzazbzazazbzbzazbzazazbzbzazazbzbzheniiiniiimiiimiii得到系统函数的偏差为 )()()()()()()()(zazazhzazbzhzhzheniiizaza1)( miiizbzb0)( 式中)()()(jjjeeheheh系数量化造成的系统频率响应的偏差系数量化的统计分析模型 可以用频响的均方偏差来描述系数量化所引起的频率特性偏差 ceejezdzzhzhjdeh)()(21)(21122均方偏差 也是

26、一个随机变量它的均值即为频响偏差 22cccnhkcmrrzdzzazazhzhjnzdzzazajmqzdzzazazhzhjaezdzzazajbe)()()()(2)()(12112)()()()(21)()(1211112111210222、fir数字滤波器系数量化的统计分析 设fir滤波器的系统函数为nnnznhzh0)()( )(nh)(nh对各系数 进行量化成为)()()(nenhnhn0n0n0)()()()(nnnnnnzneznhznhzh有：令 ， ，则n0)()(nnznezejez )()()(zhzhzennjnieneee0)()(nnnnjnnnjnineene

27、eneee000| )(|)(|)()(所以nnne0|)(|fir系数量化时，系统函数产生的误差不会超过2) 1(| )(|)(0qnneeenni所以当作定点舍入处理时，因为2| )(|qne7.4 数字滤波器的运算量化效应 w为了便于用统计方法分析这些量化误差的平均效应，我们假定n所有这些噪声都是平稳的白噪声序列；n所有噪声都与信号不相关，并且各噪声之间也互不相关n每个误差噪声都在其误差范围内呈均匀等概率分布 7.4.1 iir滤波器定点运算舍入误差的统计分析a理想相乘 b实际相乘的非线性流图 c统计模型的线性流图)(nx)(nya)(nx)(nya )( ny)(nxa)( ny)(n

28、e定点相乘运算的流图表示采用统计分析方法后，实际的输出可以表示为)()()( nenyny7.4.1 iir滤波器定点运算舍入误差的统计分析a理想相乘 b实际相乘的非线性流图 c统计模型的线性流图)(nx)(nya)(nx)(nya )( ny)(nxa)( ny)(ne定点相乘运算的流图表示采用统计分析方法后，实际的输出可以表示为)()()( nenyny方差 meceeefmhzdzzhzhj)()()(2221220)(neefnhmm kffnene1)()( kifif022)(nef而每一个噪声源 所造成的输出噪声方差及均值：最后将所有的输出噪声线性叠加就得到总的输出噪声重要结论：

29、有限字长效应与滤波器的结构型式有密切关系。例 2 已知一个iir滤波器的系统函数如下，用定点制算法，尾数舍入，分别求出直接型、级联型和并联型实现系统时量化误差的方差 。2f)8 . 01)(9 . 01 (2 . 0)(11zzzh、 1、直接型结构)(172. 07 . 112 . 0)(21zazzzh)(nx1z1z0.2)(0ne)(1ne)(2ne1.70.72)()(nenyf1z1z1.70.72)(nef)()()(210nenenea 相乘引入的舍入噪声b 三个舍入噪声通过相同的传输网络直接型的舍入噪声：系数0.2、1.7以及0.72相乘后的舍入噪声)()()(210nene

30、ne、)(1)(1zazh它们均经过相同的传输网络 输出噪声的方差是)()()()()(1210nhnenenene4123222qqe 22222112245.2280.89)8 . 01)(8 . 09 . 01)(9 . 08 . 0(8 . 0)9 . 01)(9 . 08 . 01)(8 . 09.0(9 . 0)8 . 01)(9 . 01)(8 . 01)(9 . 01 (12qzdzzzzzjeecefd)(1)(2 . 08 . 0119 . 012 . 0)(2111zazazzzh)(nx1z1z0.2)(0ne)(1ne)(2ne0.90.8)()(nenyf2、级联型

31、：可以有几种排列形式)9 . 01 ()(11zza先令 ，并把0.2置于第一级，即 )()(1)(211zazazh)(1)(21zazh222112211212178. 280.89)8 . 01)(8 . 01 (12)8 . 01)(9 . 01)(8 . 01)(9 . 01 (12eececefezdzzzjzdzzzzzj)(0ne)(1ne)(1zh)(2ne)(2zh、 通过网络 、 通过网络 输出噪声的方差是 61222221qqe12222qe22102.15qfe第一级有两个误差源，第二级有一个误差源，故有 所以 再令 ，并把0.2置于第二级，有 )8 . 01 ()(

32、11zza)(2 . 0)(19 . 012 . 08 . 011)(2111zazazzzh)(nx1z1z0.2)(1ne)(2ne0.90.8)()(nenyf)(0ne)(0ne)(1ne)(1zh)(2ne)(2zh、 通过网络 、 通过网络 输出噪声的方差)(1)(11zazh)()(1)(212zazazh2221122112212226. 5592. 3)9 . 01)(9 . 01 (12)8 . 01)(9 . 01)(8 . 01)(9 . 01 (2 . 02eececefezdzzzjzdzzzzzj第一级有两个误差源，第二级有一个误差源，故有 所以 12221qe6

33、1222222qqe221176. 1qfe结论：系数bi和a(k)的排列不同，输出误差也有所不同。 距单位圆远的极点和系数bi置于较后级的误差会小些。3、并联型 )(4 . 6)(2 . 78 . 014 . 69 . 012 . 7)(2111zazazzzh1z)(0ne)(1ne0.9)(nx7.2)()(nenyf1z)(2ne)(3ne0.86.4并联型结构需要4个系数，因此共有四个舍入噪声 、 通过网络)(0ne)(1ne)(1zh)(2ne)(3ne)(2zh、 通过网络)(1)(11zazh)(1)(22zazh21221212)26. 578. 2()9 . 01)(9 .

34、 01 (12)8 . 01)(8 . 01 (12ececefpzdzzzjzdzzzj因此输出的方差为则第一级、第二级都有两个误差源，故有 所以 61222222212qqeee2234. 1qfq由此得出结论：对iir滤波器，从有限字长效应来看，不论是哪一种型式的直接型结构都是最差的，运算误差最大，特别在高阶时应避免采用。级联型结构较好，而并联型结构具有最小的运算误差。7.4.2 iir滤波器定点加法运算的溢出问题 w所有的相加运算，均可能产生溢出，应避免w使每个相加点的输入端都引入比例因子a加以限制，使相加点的数值绝对值保持小于1。为使第k个相加节点上的输出信号 yk(n)不发生溢出，

35、需要在输入端引入比例因子ak。 mkkkmnxamhny)()()(要使yk(n)不发生溢出，则需使| yk(n) |1，所以 max| )(|xnxmkkkmhxany)(| )(|maxmkkmhxa)(1maxmaxx若 表示输入信号的最大幅度因此选择其中最小的比例因子作为系统比例因子的最终选择，即 min1kkkaa令ak=1，得到保证第k个节点上不出现溢出时的最大输入值上限 mkmhx)(1|max7.4.3 极限环振荡 wiir滤波器在无限精度的情况下，当它的所有极点均位于单位圆内时，系统肯定是稳定的。n当去掉输入信号后随着n的增加，系统输出逐渐衰减趋向于零。w在有限字长情况下，由

36、于量化过程的非线性作用，系统输出将不随n的增加而趋于零w零输入极限环振荡：系统输出不随n的增加而趋于零，而是衰减到某一非零的范围幅度后呈现振荡特性w以舍入处理的一阶iir滤波器为例来说明这一现象 设一阶iir系统的系统函数为 在定点运算中，每次乘法运算以后均要对尾数进行舍入处理。因此，实际的非线性差分方程可表示为 在无限精度运算下，其差分方程为111)(azzh)() 1()(nxnayny)() 1( )( nxnyaqnyr)(nx1za)( ny kq一阶iir网络的非线性流图 在无限精度情况下，如果输入信号)(87)(nnx只要|a|1，系统稳定，y(n)将逐渐衰减为零。 nany87

37、)(输出有限精度运算时，系数a=0.5=0.100b的情况n00.1110.0000.0000000.0000.87510.0000.1110.0111000.1000.50020.0000.1000.0100000.0100.25030.0000.0100.0010000.0010.12540.0000.0010.0001000.0010.125)(nx) 1( ny) 1( nyarnya) 1( )( ny -2 -1 0 1 2 3 4 5 6 7 n)( ny -2 -1 0 2 4 6 n)( ny 1 3 5时21ab时21aa以下也称为“死带”区域 81)( ny零输入极限环振

38、荡 2q2| ) 1( )1( |qnyanyaqr2| ) 1( |)1( |qnyanyaqr2| ) 1( | ) 1( |qnyany下面分析振荡幅度和字长的关系由于舍入误差的绝对值在 以内，因此可得或者表明：极限环幅度与量化间隔成正比， 增加字长(减小量化间隔)将使极限环振荡减弱 如b=3，故有 ， ，有8123q5 . 0|a125. 05 . 01161| ) 1( | )( |nyny一阶iir网络的死带范围与表7-4-1结果一致7.4.4 定点运算的溢出振荡 溢出可能产生比之更大的误差，并在输出的最大幅度界限内振荡，称为溢出极限环振荡 溢出极限环振荡产生的原因：定点加法运算存

39、在着溢出1-10123)(vfv补码加法运算的非线性 以定点补码运算的二阶iir滤波器为例进行讨论此二阶滤波器的差分方程为)()2() 1()(21nxnyanyany2121222211)(11)(ppzppzzzazazh24221121aaap，系统函数为首先考虑系统的稳定性1|21，p)(zh1p2p稳定的条件是 的两个极点 、 在单位圆内1a2a可得 、 的值域为图示的大三角形1a2a1|21，p即 、 落于大三角形之内才能满足 (0,1)(2,-1)(-2,-1)iii1a2a其次，为了保证不溢出，即要满足1| )2() 1(| )(|21nyanyany1|21 aa1| )(|

40、ny1|21 aa如上图中的阴影部分1a2a由此式求出 、 的关系若必定期足系统不产生溢出极限环振荡的充要条件由于溢出极限环的振荡幅度为1，数字滤波器无法对数据进行滤波，可以通过饱和型加法运算的办法加以消除。饱和加法器的特点是，当相加器的输入之和大于1或小于1时，就分别以1和1代表相加结果，其特性如图所示，这样就能克服溢出振荡1-101输出)(vf输入v-1补码饱和加法器特性 7.4.5 浮点运算中的有限字长效应 w乘法和加法之后都需要进行尾数量化。n流图中乘法支路之后或加法节点处都将叠加个量化噪声 w舍入或截尾处理只影响尾数的字长，不影响阶码，但所产生的误差值却与阶码值有关。n分析中用相对误

41、差比用绝对误差更适合 w系统的输出误差大小与系统的结构有关n本节上一部分系统结构对输出误差影响的结论同样成立 w取相同的尾数字长时，浮点运算的误差要比定点运算的误差小。n浮点制中可以充分利用尾数的字长，使每一个数不论其数值大小，都能保持同样的相对精度 w浮点制系统的字长一定时，其输出的信噪比为一常数；n定点制系统中，输出噪声的方差与信号无关，因此信号越大，输出的信噪比越大 7.5 fft算法的有限字长效应 w从运算的角度分析，fft与数字滤波器一样都是系统，因此，二者的有限字长效应分析方法基本上是相同的w不同的是fft中的运算是一个复数运算，因此误差源是一个复数，同时fft运算有多个输入端及多

42、个输出端，各误差源对不同的输出端的影响是不同的。w以时间抽取dit为例进行有限字长效应分析，并且针对的是舍入情况，其他fft算法及截尾运算结果是相似的7.5.1 蝶形运算的统计模型 )()()(1jxwixixmrnmm)()()(1jxwixixmrnmm定点制运算中只有乘法引入量化误差，加法运算不引入误差蝶形运算中只在乘系数 时引入一个误差源rnwme)( jxm)(ixm)(1ixm)(1ixmrnwme-1 是复数相乘，所产生误差 就是一个复数me)( jxwmrn4321,eeee一个复数相乘要由四个实数相乘来构成、每个实数相乘都将引入一个相应的误差。因而共有四个误差假设是舍入量化，

43、则re)(imim)(reim)(imre)(re)()(4321ewjxewjxjewjxewjxewjxwjxqrnmrnmrnmrnmmrnmrnmr)(4221eejeeem34|222423222122qeeeeeeeeeeembme2b复数 的方差 为分析乘以 后的方差影响 rnw2222|bmrnrnmeewwee通过所有后级蝶形时，其方差保持不变 me每一个输出端都与n1个蝶形相连，即有n1个量化噪声源对每个输出端有贡献 在终端，在离散傅里叶变换)(kx上叠加的输出噪声的均方值为和fir滤波器的直接型实现一样，输出噪声的总方差正比于n3) 1(2222nqnnbbf=7.5.2 防止溢出和fft输出的信噪比 分析蝶形运算)()()(1jxwixixmrnmm+=+| )(| )(| )(| )(| )(|1jxixjxwixixmmmrnmm| )(| )(| )(|1jxixjxmmm