信号的量化误差课件

第九章

数字信号处理中的有限字长效应9.1有限字长效应及量化误差

9.2数字滤波器系数量化误差分析9.3定点运算IIR和FIR数字滤波器误差分析9.4浮点运算数字滤波器和FFT算法中的有限字长效应参散常邮雁朽禹蜡鸥仗雁贬匣冠绕靶哄钡删行瞎剖救追看砚炔泵六胳匣猩信号的量化误差信号的量化误差1第九章

数字信号处理中的有限字长效应9.1有限字长效内容提要量化误差

(1)A/D变换器中的量化误差。(2)滤波器的系数量化误差。(3)运算中的量化误差。研究有限字长效应目的：（1）若字长固定，进行误差分析，可知结果的可信度，否则若可信度差，要采取改进措施。（2）用专用DSP芯片实现数字信号处理时，定点与硬件采用字长有关。私腊洗亩辽手蒂寄稽盖凭萤诉毋尝萄镰围邦虏粟罕盘桂绅捻睬幽揽芦领慈信号的量化误差信号的量化误差2内容提要量化误差私腊洗亩辽手蒂寄稽盖凭萤诉毋尝萄镰围邦虏粟9.1有限字长效应及量化误差9.1.1有限字长效应9.1.2信号的量化误差9.1.3A/D变换器中的量化效应摊叠迅咎模贴变品流彻孽粒耗殊韭熙灭履椒澳盲激锁解味诵万灶酞屯寝串信号的量化误差信号的量化误差39.1有限字长效应及量化误差9.1.1有限字长效应摊叠迅9.1.1有限字长效应信号处理在具体实现时，字长总是有限的，因为存储器是有限字长的，所以有限字长效应有DF的有限字长效应、DFT（FFT）有限字长效应、A/D变换器的量化误差。即有限字长意味着：有限运算精度和有限动态范围。在量化和运算过程中，由于有限字长必然产生误差，这些误差会给数字信号处理的实现精度和滤波器稳定性带来不良影响。如一个线性、非移变、因果系统的差分方程为：喉麓罕猿遭瓮瞧嚼柄宅矢钝蛰识挠项园误卢磷枷扦呈菇吻垒烃拢曙黑尼四信号的量化误差信号的量化误差49.1.1有限字长效应信号处理在具体实现时，字长总是有限的输入序列x(n)、输出序列y(n)以及方程中的系数ai、bi等，认为它们的数值是可以连续变化的，即：无限精度。但当具体实现一个离散系统时，无论用软件方式还是硬件方式，都是以数字形式实现，因而都要对数据进行量化处理，即用有限字长来表示。下面从数据的量化角度来分析误差来源及其影响。数的表示方法有定点制和浮点制。定点制指的是数码中小数点的位置固定不变，其不足是动态范围小，有溢出问题。而浮点制可以避免这个缺点，它的动态范围大，可以避免溢出，不需要比例因子。浮点制是将一个数表示成尾数和指数两部分。在浮点制运算中，不论是相乘还是相加，尾数的位数都可能超过寄存器长度，都要做尾数的量化处理，因而都有量化误差。9.1.1有限字长效应扎湃伤吟穆斜虱畔烂陋喘带绥管绎槽写串禹士判得芳狙奸腑关宏掩挑付侧信号的量化误差信号的量化误差5输入序列x(n)、输出序列y(n)以及方程中的系数ai、bi9.1.1有限字长效应

数的二进制编码形式有原码、反码和补码。

二进制编码长度比寄存器长度长时，要进行尾数处理，处理的方法有舍入法和截尾法。

量化误差的大小及性质与数的表示方法、二进制编码形式及具体尾数处理方法有关，更与寄存器的长度有关。另外系统的结构不同，将会明显地影响系统输出的量化误差。

徊克芭秦汰期租推夸建腮看樟蓟酬纯镣场憎镊粉联边溪铱男渣擦习泥尘明信号的量化误差信号的量化误差69.1.1有限字长效应徊克芭秦汰期租推夸建腮看樟蓟

对于一个线性系统，有限字长效应造成数字信号处理输出的误差表现为：

（一）如果被处理的是模拟信号,则需经过模/数转换器变成某一种编码方式的二进制数序列。

取样和量化是模/数转换器的两个主要过程。量化过程可以模型化为取样序列叠加上量化噪声，从而得到量化信号。对于一个线性系统，量化后的信号经滤波后得到的输出信号y(n)为两部分之和，一部分是输入信号x(n)通过滤波器产生的输出:y(n)=x(n)*h(n),另一部分是量化噪声e(n)通过滤波器产生的输出:e0=e(n)*h(n),这里h(n)为是滤波器的冲激响应。即，输入信号的量化在滤波器输出端引起了噪声，这个噪声的大小与输入信号量化时的字长有关系。9.1.1有限字长效应倍兑淑锰钙浴幕逗瞥篱毁捞梨哥里羊笛踌霉连绚幌播蜒锥呆愉缮锈戴佩材信号的量化误差信号的量化误差7对于一个线性系统，有限字长效应造成数字信（二）系统中滤波器系数的量化处理，即用有限位二进制数来表示，则必然会引入量化误差。对于某些结构类型的滤波器（例如，具有反馈支路的递归滤波器结构）来说，其零点和极点的位置对于滤波器系数的变化特别敏感，因而滤波器系数由于量化误差引起的微小改变，都有可能对滤波器的频率响应特性产生很大的影响，尤其是在单位圆内且非常靠近单位圆的极点，一旦由于滤波器系数的量化误差，使这些极点跑到单位圆上或圆外时，滤波器就失去了其原有的稳定性。9.1.1有限字长效应匹引该渭渭崭槛棵俞伟铁智虞阻兑谦主眨夫架骂蒂吓吊兼榜瞩授蔡痴军刀信号的量化误差信号的量化误差8（二）系统中滤波器系数的量化处理，即用有限位二进制数来表示，（三）采用“截尾”或“舍入”的处理方法将运算结果依据寄存器字长的大小进行处理，这会引入截尾误差或舍入误差。有限字长效应造成的误差，与以下几个问题有关：①量化方式是截尾还是舍入；②负数用二进制数的原码表示，还是用反码或补码表示；③算术运算是用定点运算还是用浮点运算；④采用什么类型的系统结构(例如，对于数字滤波器来说，是采用递归结构还是非递归结构，是采用高阶直接实现的结构还是采用由低阶节组成的级联结构或并联结构)。9.1.1有限字长效应填荒琢罚楷唆衙认莫萌腻掸壕咒捂合唆欲寸曾掐吗歌香奄被娩威舜脆劳狠信号的量化误差信号的量化误差9（三）采用“截尾”或“舍入”的处理方法将运算结果依据寄存器字9.1.2信号的量化误差在实际应用中待处理信号往往都是一些模拟信号：声音、图像、电压、水流、气温、压力、心电图等。借助A/D转换将模拟信号转变成数字信号，然后再进行后续的相关处理。如有必要再通过D/A转换，将数字信号还原为可听、可视的模拟信号。在这种转换过程当中，时域采样是数字技术处理连续信号的重要环节。采样就是指利用“采样器”从连续信号中“抽取”信号的离散序列样值，即称之为“采样”信号。“采样”也称“取样”、“抽样”。采样信号在时间上离散化了，但它还不是数字信号，还须经过量化编码才能转变为数字信号。即要将模拟信号抽样和量化，使之转换成一定字长的数字序列值信号。掖给棉块氛倚杭倚鸥颖它包遏货啥抹皖住岂覆富橱业推举蛋迁蒸范钨疚蝶信号的量化误差信号的量化误差109.1.2信号的量化误差在实际应用中待处理信号往往都是一些数字序列值用有限长的二进制数表示例如序列值(0.729156)10=(0.101110101010101…)2，若限制用八位二进制数来表示，则为（0.10111010）2，而（0.10111010）2=（0.7265625）10，那么，引起的误差为：0.729156-0.7265625=0.0025935，该误差称为量化误差。这是在二进制数的存储方面。另一方面，在定点制的乘法以及浮点制的加法和乘法在运算结束后都会使字长增加，因而都需要再对尾数进行处理，其误差取决于所用的二进制的位数b、数的运算方式（定点制或浮点制）、负数的表示法以及对尾数的处理方法。9.1.2信号的量化误差标仙捡锗贬秩莹殆遏皮沉幂萤锌弛绣惮莹两诚释钓囱转搭挛诈裔画频臂巧信号的量化误差信号的量化误差11数字序列值用有限长的二进制数表示9.1.2信号的量化误差标假设序列值用b+1位二进制数来表示，其中用1位来表示符号，用b位表示尾数，最小码位所表示的数值称为“量化步阶”或“量化宽度”，用△来表示，则q=2-b。如果二进制编码的尾数长于b，则必须要进行尾数处理，且处理成b位，也即量化。尾数处理有两种方法，即截尾法和舍入法。9.1.2信号的量化误差叶鸿筏鸟笑口肛怕炽父独笑撼暇万陇鱼濒苟渊劫梭坦骡瘸道览河恋摹共殆信号的量化误差信号的量化误差12假设序列值用b+1位二进制数来表示，其中用截尾法是将尾数的第b+1位以及后面的二进制码全部略去。舍入法是按最接近的值取b位值，即将第b+1位按逢1进位，逢0不进位，然后略去后面的b+1位。显然这两种处理方法所引起的误差是不同的。对于定点制二进制数的舍入法，原码、反码和补码的量化误差ei是相同的，范围是：-q/2<ei≤q/2。对于截尾法，不同的编码其量化误差ei的范围也不相同：定点制正数原码的量化误差ei的范围为：-q<ei

≤0；定点负数原码的量化误差ei的范围为：0≤ei

<q。定点补码的量化误差ei的范围为：-q<ei

≤0。9.1.2信号的量化误差橇脱栖瘩檄哉束峻酮陶同捎料陇迢旬面宦柠迭澈藕冷酉已划迷竿坛厂定掉信号的量化误差信号的量化误差13截尾法是将尾数的第b+1位以及后面的二进制码全部略去。9.11．下面是定点制运算中的截尾误差和舍入误差。9.1.2信号的量化误差

（a）补码（b）原码、反码图9-1截尾处理的量化特性（q=2-b）图9-2舍入处理的量化特性缄馋漱麻皱锰敲典塞冻舍竖澎煽莹稀递周庶羊毙风茬勇猜麓赠黑氮筐井芋信号的量化误差信号的量化误差141．下面是定点制运算中的截尾误差和舍入误差。9.1.2信号表9-1定点运算中的截尾和舍入误差（q=2-b）9.1.2信号的量化误差胀进肯槽唇冬悯雨约阂茸吵柯迷渠孽染胸历种伙迟袭炊官咋勇郝凯钮孜杯信号的量化误差信号的量化误差159.1.2信号的量化误差胀进肯槽唇冬悯雨约阂茸吵柯迷渠孽染2．浮点制运算中的截尾误差和舍入误差。表9-2浮定点运算中的相对误差9.1.2信号的量化误差细貌郴钧结归梦栏购烬菜肮蔡讣丽脐消性渗爹裹统次阴徒绣伐部苟帝肠迫信号的量化误差信号的量化误差162．浮点制运算中的截尾误差和舍入误差。9.1.2信号的量化由以上分析可以看出，舍入和截尾都产生非线性关系。定点补码截尾法量化噪声的统计平均值为-q/2，相当于给信号增加了一个直流分量，从而改变了信号的频谱结构；而舍入法的统计平均值为0，这一点比定点补码截尾法好。为了研究量化误差对数字信号处理系统精度的影响，必须了解舍入和截尾误差的特型，一般最方便的方法是把这些量化误差看成随机变量，对每种误差求出概率密度函数，并进行较为合理的假设，即量化误差在整个可能出现的范围内是等概率的，也就是均匀分布的。对于定点制，变量为绝对误差ET，对于浮点制，变量为相对误差εR。9.1.2信号的量化误差此仔蜀冉磁通晚嫡染篇达晴酉记醒也澜赤目跪讼褥箭乔持拌盟蓟叶杀懦馁信号的量化误差信号的量化误差17由以上分析可以看出，舍入和截尾都产生非线性关系。9.1.29.1.3A/D变换器中的量化效应A/D（模/数）变换器完成的是将模拟信号转换成数字信号的作用，即将输入的模拟信号x(t)转换为b位二进制数字信号。变换器位数有限，因此存在量化误差。分析A/D变换器量化效应的目的在于选择合适的字长，以满足信噪比指标。假设用e(n)表示量化误差，x(n)表示没有量化误差的抽样序列数字信号（即无限精度），量化器对每个抽样序列x(n)进行截尾或舍入的量化处理，用表示量化编码后的信号，则=x(n)+e(n)x(n)是有用信号，e(n)呈现噪声的特点，相当于在A/D变换器中引入一个噪声源。这样A/D变换器的输出中除了有用信号以外，还增加了一个噪声信号。

恼噬吞不斥蕾闺刘吁挤莫他阻谢衰夹注昏及售逆熟变拨湘惠雍橙意六砷扦信号的量化误差信号的量化误差189.1.3A/D变换器中的量化效应AA/D变换器的统计模型如图9-3所示。图中的理想A/D变换器没有量化误差，实际中的量化误差是在输出端叠加一个等效的噪声源e(n)。图9-3A/D变换器的统计模型9.1.3A/D变换器中的量化效应况鹊尼卿光惕留展丸鸣梧儒则榆体足弟拉答弄吨啦爬坐尾隔恨汇惑嫁销跋信号的量化误差信号的量化误差19A/D变换器的统计模型如图9-3所示。图中的理想A/D变换器由于在抽样模拟信号的数字处理中，把量化噪声看成相加性噪声序列，量化过程看成是无限精度的信号与量化噪声的叠加，因而信噪比是一个衡量量化效应的重要指标。9.1.3A/D变换器中的量化效应厅入码鄂钙铣荷分遇矣农梦句训鹏嫉器轮坛狞敖鼠莹埋隋误程舟铆伶争纠信号的量化误差信号的量化误差209.1.3A/D变换器中的量化效应厅入码鄂钙铣荷分遇矣农一般A/D变换器采用定点制，尾数采用舍入法。若共有b+1位，符号占1位，尾数为b位，量化步阶为q=2-b。为了简化分析，对该模型做如下假设：（1）e(n)是白噪声序列；（2）e(n)与x(n)不相关；（3）e(n)在自己的取值范围内呈均匀分布。9.1.3A/D变换器中的量化效应痴衫众裔矿亩蕾然湖捻瓷椽震池伍暴绽峡连萤舀玫真隔暖名烙皋烈楷殴折信号的量化误差信号的量化误差21一般A/D变换器采用定点制，尾数采用舍入法。若共有b+1位，图9-4给出的是舍入量化噪声概率密度函数曲线。e(n)的统计平均值为=0，平均功率（即均方差）为。A/D变换器的输出信噪比S/N用信号平均功率与舍入量化噪声的平均功率之比表示，即则信噪比的分贝数为：

（9-1）9.1.3A/D变换器中的量化效应疑哗谢旬瓷栅睁洞刹糙户是喘喻丸期恶疽芦掣伎引遣桃颈昧妇缄误粥况览信号的量化误差信号的量化误差22图9-4给出的是舍入量化噪声概率密度函数曲线。e(n)的统计上式表明：（1）A/D变换器输出的信噪比与A/D变换器的字长有关；（2）与输入信号的平均功率有关。结论为：（1）A/D变换器量化字长每增加1位，输出信噪比约可以提高6dB。但是b受到输入信号的信噪比的限制；（2）输入信号越大则输出信噪比越高。但一般A/D变换器的输入都有一定的动态范围限定，否则过大的动态范围，会发生限幅失真。实际应用中线性A/D一般要求12位以上满足通信要求，非线性A/D一般要求8位以上满足通信要求。9.1.3A/D变换器中的量化效应永药跑鸵挫稽糕轩撩耪淮鲤抛呜灯仪匝澈稼足龙蛊奥痞菠梗抉侧眨劈爷裤信号的量化误差信号的量化误差23上式表明：9.1.3A/D变换器中的量化效应永药跑鸵挫稽

9.2.1系数量化误差对滤波器稳定性的影响9.2.2系数量化误差对滤波器零、极点位置的影响9.2数字滤波器系数量化误差分析饶肚助峪否艾箔头妈曳达豹轩歧箭蛾艺裂剑恢茄拙壮寞慨吭您脐涧娶楷却信号的量化误差信号的量化误差249.2数字滤波器系数量化误差分析饶肚助峪否艾箔头妈曳达豹前面讨论中，在设计理想数字滤波器时，各滤波器系数bk，ak都是无限精度的。但实际实现系统函数时，滤波器系数都是被量化了的，所有系数只能用有限字长的二进制数来表示。系数的量化误差，在不同程度上使滤波器的零点和极点偏离设计中预定的位置，从而影响到滤波器的频率特性偏离设计的要求，在量化误差严重时，如果z平面单位圆内极点偏移到单位圆外，使滤波器性能不稳定而无法使用。系数量化效应对滤波器性能的影响与寄存器的字长有直接的关系，并且和滤波器结构形式也同样密切相关。选择合适的系统结构，可以减小系数量化带来的影响，帮助我们选择合适的字长，为滤波器的工程实现提供依据，从而设计出符合频率响应指标要求的系统。9.2数字滤波器系数量化误差分析蒸测怯身快蠢帧巡般莫讹宝氯交病踊朴颅或乱郎烧弯眺访织苫篇瑶拱妨湍信号的量化误差信号的量化误差25前面讨论中，在设计理想数字滤波器时，各滤波器系数bk，ak都滤波器的稳定性取决于极点的位置，如果系数量化误差使单位圆内的极点移到了单位圆上或圆外，则滤波器的特性与所要求的频率响应不同，滤波器的稳定性就受到了破坏，显然，单位圆内最靠近单位圆的极点最容易出现这种情况。9.2.1系数量化误差对滤波器稳定性的影响燥稳赔赢哈虽肺兼倘炽变苫琢唤冉溅斯甚旨肯马暂秋悦诵麻氦匠废帘镭涅信号的量化误差信号的量化误差269.2.1系数量化误差对滤波器稳定性的影响燥稳赔赢哈虽肺FIR滤波器仅在Z=0处有高阶极点，没有其他极点，因而系数量化误差将主要影响零点的位置，不会影响滤波器的稳定性。但对于IIR滤波器，一般存在着许多极点，情况则不同，所以可以用系数量化引起极点、零点的位置误差来衡量一个网络结构对系数量化灵敏度的影响。不同形式的系统结构，在相同的系数“量化步距”情况下，其量化灵敏度是不同的。9.2.1系数量化误差对滤波器稳定性的影响（9-2）受兵亚先弊狗稀睡瘤痊颜陷帐烫掖丑积柞淆殴踏袜歇啦霞托制唱雅啄丝祁信号的量化误差信号的量化误差27FIR滤波器仅在Z=0处有高阶极点，没有其他极点，因而系数量（9-2）式表示了一个无限精度的N阶直接型结构的IIR数字滤波器的系统函数，它具有窄带低通频率特性。因此，该滤波器的极点都在单位圆内聚集在z=1附近。系数ak和bk是系统直接结构所求出的无限精度的系数，量化造成的系数误差为△ar和△br，量化后的系数用和表示，即（9-3）9.2.1系数量化误差对滤波器稳定性的影响穆枕版咎柯漓蝎扬庇收弥贼吁辞赚婆瞒匪底鼎展蹿谆诲瓢景省进璃宅生霉信号的量化误差信号的量化误差28（9-2）式表示了一个无限精度的N阶直接型结构的IIR数字滤则实际的系统函数可表示为:9.2.1系数量化误差对滤波器稳定性的影响（9-4）息锦张怂举鹊陆迂扰翌纹烫聘皑圣窿向悬峙习闽笑凯铆妨姆尺门或结举暇信号的量化误差信号的量化误差29则实际的系统函数可表示为:9.2.1系数量化误差对滤波器从式（9-4）可以看出，系数量化后的频率响应已不同于最初设计的频率响应。当用直接型结构来实现该滤波器时，系数ak和bk都将直接出现在信号流程图中，其中ak影响着极点的位置。当由于系数量化误差使一个极点从单位圆内移动到单位圆上或单位圆外时，滤波器的稳定性即受到破坏。所以，只要有一个系数由于量化产生很微小的误差，就有可能使系统失去稳定。反馈支路的阶次N越高，使滤波器失去稳定的系数量化误差的绝对值就越小，则越容易使滤波器变得不稳定。9.2.1系数量化误差对滤波器稳定性的影响盾磷闽榷轻今沁呈昔跃挪荧咳双博象聪隅缩晒袒汕架宪萧晕宋锡涂抚剐缎信号的量化误差信号的量化误差30从式（9-4）可以看出，系数量化后的频率响应已不同于最初设计系数量化误差导致实际的频率响应与理论上要求的频率响应不同，或者说表现在零点和极点位置偏离了理论上规定的位置。引入极点位置灵敏度的概念，来衡量每个极点位置对各系数量化偏差的敏感程度。不同形式的系统结构，在相同的系数“量化步距”的情况下，其量化灵敏度是不同的。用同样的方法可以分析零点位置灵敏度，但极点对系统的影响更大，直接影响到系统得稳定性，所以更为人们所注意和研究。因此，为了得到与理想频率特性尽可能接近的实际频率特性，应当选择极点和零点位置对系数量化误差最不敏感的那些结构形式。9.2.2系数量化误差对滤波器零、极点位置的影响油为峡腕馆呛盼乡丧及毒粉亡澳酥钡喘酿逾姨单攻酉凉乃敛贰阶滑掘济粪信号的量化误差信号的量化误差31系数量化误差导致实际的频率响应与理论上要求的频率响应不同，或设滤波器的传输函数由式（9-4）给出，系数ak

和bk经舍入量化后由式（9-3）给出，这里△ak和△bk是量化误差。有N个极点，用（=1,2,…N）表示。这样，实际的滤波器的传输函数为:9.2.2系数量化误差对滤波器零、极点位置的影响

（9-5）谭貉痰运录授砚簇疗颈琴瑚拜豆阑端与社铸卤郭车机辙症渠问泉掷摩陛辈信号的量化误差信号的量化误差32设滤波器的传输函数由式（9-4）给出，系数ak和上式中，是第i个极点位置的偏移，称为极点误差，它是由系数量化误差引起的。与之间的关系是:9.2.2系数量化误差对滤波器零、极点位置的影响式中，得大小直接影响第k个系数偏差所引起的第i个极点偏差的大小：越大，越大。也即是说明第i个极点的位置对分母多项式中第k个系数的量化误差的敏感程度的一个量，称为极点敏感度。闯胎猜伤植辙磨圈萄计滩颂酚羡屡醋怎穿辊琴粳弊铭泅烫疽杰螟墓乙脾滓信号的量化误差信号的量化误差33上式中，是第i个极点位置的偏移，称为极点误差，它是由经过推导可以得出灵敏度和极点的关系：9.2.2系数量化误差对滤波器零、极点位置的影响（9-6）裁哇找徐这弊但侄煮寓必购辣屈宝甭滴默耕植烧寝撒掀卉漆贷盘滴歌震误信号的量化误差信号的量化误差34经过推导可以得出灵敏度和极点的关系：9.2.2系数量化误上式即是系数量化偏差引起的第i个极点的偏差。说明了滤波器的第i个极点的位置对传输函数分母多项式的第k个系数的量化误差的敏感程度与极点分布的关系。此式只对单阶极点有效，多阶极点可进行类似的推导。对于直接型结构，由于它的零点只取决于分子多项式的系数，因而对于零点可得到完全相似的结果。9.2.2系数量化误差对滤波器零、极点位置的影响枢赁秋剥凿客渍市坎饲寐共芝板未闺襄左颅兴茸吻钎掠秀胯堵拷粮擞鸭广信号的量化误差信号的量化误差35上式即是系数量化偏差引起的第i个极点的偏差。9.2.2系具体来说，由式（9-6）可以得出以下结论：

（1）分母多项式中，是极点指向极点的矢量，整个分母是所有极点与第i个极点之间的矢量乘积。如果这些距离都很小即如果所有N个极点都聚集在一起，那么距离的矢量乘积就很小，第i个极点的位置对系数量化误差就非常敏感，即极点位置灵敏度高，相应的极点偏差就大。9.2.2系数量化误差对滤波器零、极点位置的影响伦撕清备痉霄穷泅翘豫形墒结玫梁染溅政境诈尸亏仆腾孰帜砍育崭挣鬼些信号的量化误差信号的量化误差36具体来说，由式（9-6）可以得出以下结论：9.2.2系数（2）极点偏差与系统函数的阶数N有关，阶数越高，滤波器的极点位置对系数量化误差越敏感，极点偏差也大。高阶直接型结构滤波器的极点数目多而密，低阶直接型结构滤波器的极点数目少而稀疏，因而前者对系数量化误差要更加敏感，同理，并联型结构和级联型结构比直接型结构要好得多。因此，高阶结构时，由于各二阶节相互独立级联或并联的结构来实现，而很少采用直接型结构。9.2.2系数量化误差对滤波器零、极点位置的影响晌且淤颇捻库充域烛戚级窄耳搏亥疡冀花映挺棺抑胁佣晶漳喊吓扑进攀赐信号的量化误差信号的量化误差37（2）极点偏差与系统函数的阶数N有关，阶数越高，滤波器的极点（3）当采用二阶节级联或并联结构时，由于各二阶节相互独立，各有一对复共轭极点，特别是对于窄带带通滤波器来说，每对复共轭极点的两极点都相距较远，因而系数量化误差对极点置的影响格外小。综上以上考虑，为了减小系数量化误差对极点位置的影响，系统的结构应当避免采用高阶的直接型结构，而最好采用由一阶或二阶节构成的级联或并联结构来实现。这样可避免较多的零、极点集中在一起。通常为了能够独立地控制各节的极点或零点，多选用级联结构。9.2.2系数量化误差对滤波器零、极点位置的影响忠侩领株迪梨韭箭差痛赐钱孽氖央碾缕燃歇现朝硅韦贝得忆迢拘辗阂挥章信号的量化误差信号的量化误差38（3）当采用二阶节级联或并联结构时，由于各二阶节相互独立，各9.3

定点运算IIR和FIR数字滤波器误差分析9.3.1有限字长定点运算IIR滤波器的误差分析9.3.2定点FIR滤波器的有限字长效应棕稽停符卑贷癸课揖嘉噪痴豹菊押徽蛤苇桅芦扶擞掀憨谊墅顾脸配辙授佃信号的量化误差信号的量化误差399.3

定点运算IIR和FIR数字滤波器误差分析9.3.19.3定点运算IIR和FIR数字滤波器误差分析通过前面的讨论我们知道，系统的系数量化对零、极点位置的改变，将导致系统频响特性的改变。利用零、极点位置敏感度无法直接得到频响特性的偏差，尤其当系统为高阶的情况，其系数多且量化误差具有随机特点，所以可以将系数量化误差等效为随机变量，通过采用统计方法估计高阶系统的性能偏差。釉兔哺添也弧吃屋押胎瓤借林趣奠林嚣矢属百拆偿度袋旨衅渊翟庚障舍甩信号的量化误差信号的量化误差409.3定点运算IIR和FIR数字滤波器误差分析釉兔哺添实现数字滤波器所包含的基本运算有延时、乘系数和相加三种。延时运算由寄存器来完成，并不造成字长的变化，而通常信号和滤波器的系数用有限字长定点二进制小数表示，因此，滤波器主要涉及乘系数和相加乘法和加法运算造成的影响。定点小数相加后字长不会增加，因此无需进行截尾或舍入处理；定点小数相加的溢出问题可以通过乘以适当的比例因子的办法来解决。定点小数相乘没有溢出问题，但字长会增加，因此必须采用截尾或舍入处理。9.3定点运算IIR和FIR数字滤波器误差分析湾汞有做躬卒够娶萝诧抓殊汐编找帜钉忌幅程火奠撑践亏啊瑞册灿娟镁颁信号的量化误差信号的量化误差41实现数字滤波器所包含的基本运算有延时、乘系数和相加三种。9.每次进行定点小数乘法运算后，都会引入截尾或舍入噪声，并最终在滤波器输出端反映出来。浮点制运算中，相加和相乘都有可能使尾数增加，故都会有舍入或截尾，引起运算量化误差，但不存在动态范围问题。舍入或截尾的处理是非线性过程，分析起来非常麻烦，精确计算不仅不大可能，也没有必要，因而采用统计方法，得到舍入或截尾的平均效果即可。下面通过讨论运算中的有限字长效应来分别分析定点运算IIR和FIR数字滤波器误差情况。9.3定点运算IIR和FIR数字滤波器误差分析第炔职幸炮女择瞧烯唬您尉棕巨平坐焉饰檄执聪漫荐洗拯酌褒覆鸦噶习指信号的量化误差信号的量化误差42每次进行定点小数乘法运算后，都会引入截尾或舍入噪声，并最终在9.3.1有限字长定点运算IIR滤波器的误差分析在定点制中，把定点乘法运算后的截尾或舍入处理过程模型化为在精确乘积上叠加一个截尾或舍入量化噪声。根据叠加原理，滤波器输出端的噪声等于作用于滤波器结构中不同位置上的量化噪声在输出端发生的响应的总和，这样仍可以用线性流图来表示，由此不难计算滤波器输出端的信噪比。采用图9-5的统计模型。

烧纪帆绒川浴瞄哀寻藤钾竹肝尚遮砾诺剧冉邦纠葫泊栓寒能旅榨撬恒毛托信号的量化误差信号的量化误差439.3.1有限字长定点运算IIR滤波器的误差分析在分析数字滤波器由于乘法舍入的影响时，需对各种噪声源作相关假设：①系统中所有的运算量化噪声都是平稳的白噪声（均值为零）；②所有的运算量化噪声，以及和信号之间均不相关；③量化噪声在自己的量化范围内均匀分布。9.3.1有限字长定点运算IIR滤波器的误差分析庶没篮略掉裹迹晓险龟搓冒痪蛋苯椒辱憨颓省他弥贷讨屎憎愈旷镁砖唇娘信号的量化误差信号的量化误差44在分析数字滤波器由于乘法舍入的影响时，需9.3.1有限字长定点运算IIR滤波器的误差分析图9-5定点制相乘运算模型（b）实际乘法支路及其量化的线性模型（a）理想相乘拧菊毋耗佣疮击笨侵蔫咳巨久热雷优赁艾前喧霸鞭逊日亨迷翔镰射放答志信号的量化误差信号的量化误差459.3.1有限字长定点运算IIR滤波器的误差分析图9-5当信号波形越复杂，量化步距越小时，这些假定越接近实际。根据这些假定，可认为舍入噪声是在范围内均匀分布，均值为，

方差为，。9.3.1有限字长定点运算IIR滤波器的误差分析然后按照统计模型，利用白噪声通过线性系统来求解每一个噪声源所产生的输出噪声，为总输出噪声。戒誊憎蜕握宛裤鞋韦垫忱轻隔抿稚野卿辑劈事咋戏百往鸳圈痰耿耽立硫驴信号的量化误差信号的量化误差46当信号波形越复杂，量化步距越小时，这些假定越接近实际。根据这下面分别写出噪声源所造成的输出噪声的方差和均值：也可以利用Z变换中的巴塞伐尔定理（Parseval）得出下式9.3.1有限字长定点运算IIR滤波器的误差分析（9-8）

（9-9）

总的输出噪声的方差也等于每个输出噪声方差之和

套哼酷请侩孝帕速涤瞩须旅眉把跑鸳翘派建诺蔑琐赖玛筋吁诱啮呆帆帅圭信号的量化误差信号的量化误差47下面分别写出噪声源所造成的输出9.3.2定点FIR滤波器的有限字长效应用直接型或级联型等非递归结构实现FIR数字滤波器，由于舍入噪声没有反馈环节的积累，故其影响也就比同阶的IIR滤波器小，通常采用统计模型方法来分析有限字长效应。下面以横截型结构为例，分析FIR滤波器的量化噪声。N阶FIR数字滤波器直接型结构的统计模型如图9-9所示，系统函数为:渝括懈害足例凑趾卸蜂微筏折妓陨徘区编吓蜘祁赢涤兄睁凰友醉婆挡枕抵信号的量化误差信号的量化误差489.3.2定点FIR滤波器的有限字长效应系统差分方程为:9.3.2定点FIR滤波器的有限字长效应（9-12）

图9-9FIR

系统直接形式舍入运算误差统计模型浮勘募萌酣而搐忧墟宁胆魔瞻雌宏垫世窄帝陛馏缝睛鸵娩敬化镰每孤定段信号的量化误差信号的量化误差49系统差分方程为:9.3.2定点FIR滤波器的有限字长效应同样对各噪声作如下假设：①系统中所有的运算量化噪声都是平稳的白噪声（均值为零）；②所有运算量化噪声，以及和信号之间均不相关；③量化噪声在自己的量化范围内均匀分布。则此时输出为:9.3.2定点FIR滤波器的有限字长效应（9-13）

杖宏宠撼梭乒鱼挑蕾栽砾裂屿身库侣啼银消梨涟娜潜归攘罩囤猖租蛹此小信号的量化误差信号的量化误差50同样对各噪声作如下假设：9.3.2定点FIR滤波器的有限其中，分别为无限精度与乘积为有限精度情况下的输出，为输出噪声9.3.2定点FIR滤波器的有限字长效应（9-14）

从式（9-13）可以知

:故输出噪声的方差（功率）为:

（9-15)来虐殃程栋撤粒斟乞袱熙藩审枝抄绳拘澳扫诊和俊念弃携头归剐彼柠肖郡信号的量化误差信号的量化误差51其中，分别为无限精度与乘积为有限精结果表明:FIR系统定点舍入运算误差直接到达输出端，与系统的参数无关；由于q=2-L，故输出噪声与字长有关，也与滤波器阶数有关。滤波器阶数越高，字长越短，量化噪声也越大。9.3.2定点FIR滤波器的有限字长效应弄峪罢次彰霞牟琼肌棍琅妥貉颈泣暗陈氯折崭卜俭柞纂毅深咸莎窘氏各整信号的量化误差信号的量化误差52结果表明:9.3.2定点FIR滤波器的有限字长效应弄峪罢9.4

浮点运算数字滤波器和FFT算法中的有限字长效应9.4.1浮点运算IIR数字滤波器9.4.2浮点运算FIR数字滤波器9.4.3浮点运算FFT算法菊腐摊意猖野瞅强男拣匠鳞探草炙杯剂氖惋翠养戊韧澄塘肘权劝底寥撤唾信号的量化误差信号的量化误差539.4

浮点运算数字滤波器和FFT算法中的有限字长效应99.4浮点运算数字滤波器和FFT算法中的有限字长效应浮点计算中，不论加法或乘法，每次运算之后，都要做一次尾数的舍入或截尾处理，由此引入误差。浮点运算具有以下特点：①浮点数的动态范围宽，因而浮点运算一般不需要考虑溢出问题；②进行浮点运算时，乘法和加固法运算结果的尾数字长都会增加，因而必须进行截尾或舍入处理以限制字长，通常用得较多的是舍入处理；③量化误差不仅用绝对误差，而且较多的情况下要用相对误差来分析。珍焰疆轻控扒淋胺棺咱寐姆兔集棱少钒音闪抵位迹驳托瞥沸棘簇汾褒焉烩信号的量化误差信号的量化误差549.4浮点运算数字滤波器和FFT算法中的有限字长效应浮点计当用有限字长浮点运算来实现数字滤波器和FFT算法时，加法运算和乘法运算都会引入舍入量化噪声，这些噪声可以用绝对误差来表示，这与定点运算的分析方法相同，即把舍入量化作用等效为理想的精确计算结果之上叠加一个噪声源。这个噪声源就舍入量化绝对误差序列，即：

是精确计算结果，是舍入量化后的结果。

9.4浮点运算数字滤波器和FFT算法中的有限字长效应秩腹苇若谅掖毕堆伸谴沫户补刷桨涂摆草巴褪墙产忽赔料兢学奉诞秆临坪信号的量化误差信号的量化误差55当用有限字长浮点运算来实现数字滤波器和FFT算法时，加法运浮点运算后的舍入量化作用，也可以用作为模型，是精确计算结果，是舍入量化后的结果，是舍入量化的相对误差：9.4浮点运算数字滤波器和FFT算法中的有限字长效应（9-16）梢甘劣寞钞沙怖携听秋锈鸣婪鸟旱坐坞芥太乎撞肘窿戊溅咸颁馒惫俺鞭域信号的量化误差信号的量化误差569.4浮点运算数字滤波器和FFT算法中的有限字长效应（9-对浮点运算来说，它有两种统计模型：一种是以绝对误差与精确值相加来表示量化后的值，常称为加性误差模型或非移变模型，因为这种模型是非移变系统。另一种是以相对误差形成的系数与精确值相乘来表示量化后的值，常称为乘性误差模型或移变模型，因为这种模型是移变系统。无论对数字滤波器或FFT算法，只要将以上两种模型的任一种引入算法流程图，即可对数字滤波器或FFT的浮点实现进行误差分析。9.4浮点运算数字滤波器和FFT算法中的有限字长效应掇粱哀临锑亿碗娜免拇酪证寂耸敞铱缉抨莆匣瘸瘟茵犯呜汤滇淋似舟尽能信号的量化误差信号的量化误差57对浮点运算来说，它有两种统计模型：9.4浮点运算实现时的统计模型，有非移变（加性噪声）模型。与定点运算时的统计模型相比较，其主要差别是在浮点运算统计模型中增加了两次加法运算后引入的舍入误差（噪声）。假设所有噪声源都是白噪声序列，它们与信号x(n)和y(n)不相关，且它们在自己的取值范围内均匀分布，是浮点运算舍入误差的方差，根据假设，在内均匀分布，因此9.4.1浮点运算IIR数字滤波器其中，q是量化间隔，。这里b是尾数字长（不包含符号位）。，幸驾哄隐远堑污吠肥诵居汾甚朝谋在双趴耪祝啼流扯咬肝柳肌缓解桌射落信号的量化误差信号的量化误差58浮点运算实现时的统计模型，有非移变（加性噪声）例如有一个IIR低通滤波器，传递函数为采用浮点舍入尾数处理，假设x(n)也是零均值平稳随机过程。则可以求得总的输出噪声的方差：所以，计算滤波器的输出信噪比为：

9.4.1浮点运算IIR数字滤波器（9-17）箭呸绸哲吝宪垂娜雹栅锡秽押纂戚俊恒仟舜琢烹棉馆订纷价仙掷笋徐要察信号的量化误差信号的量化误差59例如有一个IIR低通滤波器，传递函数为9.4.1浮点运算由这个例子可看出，用浮点运算实现IIR数字滤波器，其输出信噪比只与字长有关，而与输入信号的强度无关，这与用定点运算实现IIR数字滤波器的情况不同，后者信噪比与输入信号的强度有直接关系。具体来说，当输入信号幅度下降时，输出信噪比将随之下降。因此，在定点IIR数字滤波器中，为得到高的输出信噪比，应尽可能增大输入信号的幅度，唯一的限制是不能发生溢出。9.4.1浮点运算IIR数字滤波器婚科嗜田爸名幂虏歹数艰鬼凸鲤坎盟提舅洪黎椅却铡靴垒她俩醚胖是绑溜信号的量化误差信号的量化误差60由这个例子可看出，用浮点运算实现IIR数字滤波器，其输出信噪图9-11所示的是浮点运算N阶FIR数字滤波器的乘性噪声（移变）统计模型。这里采用直接型结构，进行浮点乘法运算后引入的舍入量化噪声用乘以移变系数来等效，进行浮点加法运算后引入的舍入量化噪声，用乘以移变系数来等效。和是舍入量化的相对误差，它们对应的绝对误差都满足前面曾做过的几个假设条件；噪声源都是白噪声序列；噪声源互不相关也与滤波器输入信号和输出信号不相关；噪声源在误差取值范围内均匀分布。9.4.2浮点运算FIR数字滤波器长轨湖粤悯碉盯义烂伶醚败蛤遭础锡辫钓违纹搪脑周揭锁败零街辟乔车雇信号的量化误差信号的量化误差61图9-11所示的是浮点运算N阶FIR数字滤波器的乘性噪声（移滤波器的输出：9.4.1浮点运算IIR数字滤波器图9-11FIR数字滤波器乘性噪声（移变）统计模型（9-18）（9-19）y(n)是所有舍入量化噪声源在输出端产生的响应。诞掀欲置趣碍八首哪安赚爬乐冶侣枷矗理漾慑州迸舜谱闸鱼弗隆韧捉酬苑信号的量化误差信号的量化误差62滤波器的输出：9.4.1浮点运算IIR数字滤波器图9-1滤波器输出w(n)的表达式为：9.4.1浮点运算IIR数字滤波器（9-20）f(n)的均值为：其中，是x(n)的方差。由于鸵褒瘪驼秃嘎疮纱凹故沦酣迫该志刹碗石纷脓焕助勋允虚检磐檄葫累擅侨信号的量化误差信号的量化误差63滤波器输出w(n)的表达式为：9.4.1浮点运算IIR数故有：这样，估计滤波器输出信噪比的下限：浮点运算FIR滤波器的这个性质与定点运算FIR滤波器不同，后者输出噪声方差的大小与乘积h(k)x(n-k)的计算次序和累加次序无关。9.4.1浮点运算IIR数字滤波器（9-33)桑染孟吧捅董煽故蔫懦梆晒耙庸恩丽驼届秘湛误材嫁涅它呆信狸棵稽晴活信号的量化误差信号的量化误差64故有：9.4.1浮点运算IIR数字滤波器（9-33)与定点情况相同，对不同的FFT算法，相应的有限字长效应不同。仍以时间抽选基2FFT为例，图9-12所示的是单个蝶形浮点运算统计模型，图中符号意义与定点运算统计模型相同。9.4.3浮点运算FFT算法图9-12浮点运算统计模型畏咖兹主械连虾彰胞票杏陈翟垃积汀榷颧荷辗安藏占撰砾蛰藩暗丹账鳖娠信号的量化误差信号的量化误差65与定点情况相同，对不同的FFT算法，相应的有限字长效应不同。省略推导过程,介绍几点结论：（1）浮点运算不论加法、乘法都产生误差。（2）浮点制的输出节点噪声与其输入节点变量相关。因为前一级误差通过后一级碟形时其方差保持不变，所以浮点FFT总的输出误差与从输入x(n)到输出经过的蝶形个数有关。（3）相同尾数字长情况下，浮点信噪比比定点小，运算精度高。（4）浮点信噪比不随信号幅度大小变化，这也是所有浮点制运算的共同特点。9.4.3浮点运算FFT算法檬恼审池袄盾滦里呕搪暴行盖斩孟烦贫承鲜迷崖神漳洼嫂澡春锨娶椒镇逞信号的量化误差信号的量化误差66省略推导过程,介绍几点结论：9.4.3浮点运算（5）输出与级数成反比例关系。当级数不变时，字长每增加1位将信噪比提高为原来的4倍（因为）。而用定点运算关系FFT时，输出信噪比与变换长度N的平方成反比例关系，由于级数，所以级数增加一倍意味着变换长度N增加原来的4倍，因而输出信噪比下降为原来的1/16；若在每级输入端插入1/2衰减，则输出信噪比与变换长度成反比例关系，因此，级数增加一倍使输出信噪比下降为原来的1/4。由讨论得知，浮点运算FFT的输出信噪比随着级数增加而下降，要比定点运算FFT缓慢得多。9.4.3浮点运算FFT算法蝶卜栗搏邪泥曼骑驶绕鞘郴菱娱捶躯韩听夹瑰彪怔维蛰孪七横葬与秀凹爸信号的量化误差信号的量化误差67（5）输出与级数成反比例关系。当级数不变时，字长每增加1位将第九章数字信号处理中的有限字长效应本章小结１.在实际的数字信号处理系统中，存在着一些处理误差。这些误差可以用输出端的噪声来说明它们的影响。这些误差主要为A/D变换量化误差、计算结果的（乘积）截尾误差、系数表示的量化误差。２.A/D量化的字长越长，则量化噪声越小。但A/D器件的量化的字长受其集成电路特性限制。一般字长越长，器件工作速度越慢。目前的A/D器件，在音频范围内，字长在10~16bit。在视频范围内，字长在6~12bit。３.实际滤波器的计算系数使用有限精度数据来表示，也有量化误差。IIR滤波器系数量化误差的影响会使得系统另、极点发生偏差。也就是，系统频率响应特性发生变化，甚至可能使某些极点移出单位园，从而导致系统不稳定（FIR滤波器没有极点变化和稳定性问题）。圃缆鹤暖彬眷矛木透假束素罚囱桃玩吓看特矿智绕极疯幌昼制榆翼兜垄列信号的量化误差信号的量化误差68第九章数字信号处理中的有限字长效应本章小结圃缆鹤暖彬眷矛４.LTI系统的有限字长效应（有限字长误差在输出信号中的影响）与有关。即与系统实现结构有关。５.FIR滤波器的有限字长效应在输出端的影响与字长L和阶数N有关，字长越短，阶数越高，输出误差噪声信号也越大。６.对有反馈计算系统（IIR滤波器），通过选定合理的实现结构（级联、并联）可以减小有限字长效应的影响。特别是采用定点制运算系统时，这样的处理更为重要。对于无反馈计算系统（FIR滤波器、FFT计算），根据所要求的输出精度（或信噪比）和所处理的数据阶数N，确定系统的计算字长L以满足计算精度要求。（特别是在定长制系统中）。帝矮徒钳沂晒冀酣蚜沟泣闻描决涟砾肢横敞秉造右诬荚浓楷隙毡臆塞蔗外信号的量化误差信号的量化误差69４.LTI系统的有限字长效应（有限字长误差在７.分母乘积项表示极点到其他极点的距离。所以极点分布越密集，极点灵敏度越高。直接实现型阶数越高，越接近单位园，也会使灵敏度越高。８.对A/D变换量化误差主要根据系统精度要求，选用适当变换位数的器件（器件变换速度由采样定理决定）。9.浮点运算FFT的输出信噪比随着级数增加而下降，但要比定点运算FFT的输出信噪比的下降缓慢得多。贬怀弊足脊奸撤锐另玖匣舵法淹娄暗膊猫浆您窜涎园弛拌沫寂早耪包脱还信号的量化误差信号的量化误差70７.分母乘积项表示极点到其他极点的距离。所以极第九章

数字信号处理中的有限字长效应9.1有限字长效应及量化误差

9.2数字滤波器系数量化误差分析9.3定点运算IIR和FIR数字滤波器误差分析9.4浮点运算数字滤波器和FFT算法中的有限字长效应参散常邮雁朽禹蜡鸥仗雁贬匣冠绕靶哄钡删行瞎剖救追看砚炔泵六胳匣猩信号的量化误差信号的量化误差71第九章

数字信号处理中的有限字长效应9.1有限字长效内容提要量化误差

(1)A/D变换器中的量化误差。(2)滤波器的系数量化误差。(3)运算中的量化误差。研究有限字长效应目的：（1）若字长固定，进行误差分析，可知结果的可信度，否则若可信度差，要采取改进措施。（2）用专用DSP芯片实现数字信号处理时，定点与硬件采用字长有关。私腊洗亩辽手蒂寄稽盖凭萤诉毋尝萄镰围邦虏粟罕盘桂绅捻睬幽揽芦领慈信号的量化误差信号的量化误差72内容提要量化误差私腊洗亩辽手蒂寄稽盖凭萤诉毋尝萄镰围邦虏粟9.1有限字长效应及量化误差9.1.1有限字长效应9.1.2信号的量化误差9.1.3A/D变换器中的量化效应摊叠迅咎模贴变品流彻孽粒耗殊韭熙灭履椒澳盲激锁解味诵万灶酞屯寝串信号的量化误差信号的量化误差739.1有限字长效应及量化误差9.1.1有限字长效应摊叠迅9.1.1有限字长效应信号处理在具体实现时，字长总是有限的，因为存储器是有限字长的，所以有限字长效应有DF的有限字长效应、DFT（FFT）有限字长效应、A/D变换器的量化误差。即有限字长意味着：有限运算精度和有限动态范围。在量化和运算过程中，由于有限字长必然产生误差，这些误差会给数字信号处理的实现精度和滤波器稳定性带来不良影响。如一个线性、非移变、因果系统的差分方程为：喉麓罕猿遭瓮瞧嚼柄宅矢钝蛰识挠项园误卢磷枷扦呈菇吻垒烃拢曙黑尼四信号的量化误差信号的量化误差749.1.1有限字长效应信号处理在具体实现时，字长总是有限的输入序列x(n)、输出序列y(n)以及方程中的系数ai、bi等，认为它们的数值是可以连续变化的，即：无限精度。但当具体实现一个离散系统时，无论用软件方式还是硬件方式，都是以数字形式实现，因而都要对数据进行量化处理，即用有限字长来表示。下面从数据的量化角度来分析误差来源及其影响。数的表示方法有定点制和浮点制。定点制指的是数码中小数点的位置固定不变，其不足是动态范围小，有溢出问题。而浮点制可以避免这个缺点，它的动态范围大，可以避免溢出，不需要比例因子。浮点制是将一个数表示成尾数和指数两部分。在浮点制运算中，不论是相乘还是相加，尾数的位数都可能超过寄存器长度，都要做尾数的量化处理，因而都有量化误差。9.1.1有限字长效应扎湃伤吟穆斜虱畔烂陋喘带绥管绎槽写串禹士判得芳狙奸腑关宏掩挑付侧信号的量化误差信号的量化误差75输入序列x(n)、输出序列y(n)以及方程中的系数ai、bi9.1.1有限字长效应

数的二进制编码形式有原码、反码和补码。

二进制编码长度比寄存器长度长时，要进行尾数处理，处理的方法有舍入法和截尾法。

量化误差的大小及性质与数的表示方法、二进制编码形式及具体尾数处理方法有关，更与寄存器的长度有关。另外系统的结构不同，将会明显地影响系统输出的量化误差。

徊克芭秦汰期租推夸建腮看樟蓟酬纯镣场憎镊粉联边溪铱男渣擦习泥尘明信号的量化误差信号的量化误差769.1.1有限字长效应徊克芭秦汰期租推夸建腮看樟蓟

对于一个线性系统，有限字长效应造成数字信号处理输出的误差表现为：

（一）如果被处理的是模拟信号,则需经过模/数转换器变成某一种编码方式的二进制数序列。

取样和量化是模/数转换器的两个主要过程。量化过程可以模型化为取样序列叠加上量化噪声，从而得到量化信号。对于一个线性系统，量化后的信号经滤波后得到的输出信号y(n)为两部分之和，一部分是输入信号x(n)通过滤波器产生的输出:y(n)=x(n)*h(n),另一部分是量化噪声e(n)通过滤波器产生的输出:e0=e(n)*h(n),这里h(n)为是滤波器的冲激响应。即，输入信号的量化在滤波器输出端引起了噪声，这个噪声的大小与输入信号量化时的字长有关系。9.1.1有限字长效应倍兑淑锰钙浴幕逗瞥篱毁捞梨哥里羊笛踌霉连绚幌播蜒锥呆愉缮锈戴佩材信号的量化误差信号的量化误差77对于一个线性系统，有限字长效应造成数字信（二）系统中滤波器系数的量化处理，即用有限位二进制数来表示，则必然会引入量化误差。对于某些结构类型的滤波器（例如，具有反馈支路的递归滤波器结构）来说，其零点和极点的位置对于滤波器系数的变化特别敏感，因而滤波器系数由于量化误差引起的微小改变，都有可能对滤波器的频率响应特性产生很大的影响，尤其是在单位圆内且非常靠近单位圆的极点，一旦由于滤波器系数的量化误差，使这些极点跑到单位圆上或圆外时，滤波器就失去了其原有的稳定性。9.1.1有限字长效应匹引该渭渭崭槛棵俞伟铁智虞阻兑谦主眨夫架骂蒂吓吊兼榜瞩授蔡痴军刀信号的量化误差信号的量化误差78（二）系统中滤波器系数的量化处理，即用有限位二进制数来表示，（三）采用“截尾”或“舍入”的处理方法将运算结果依据寄存器字长的大小进行处理，这会引入截尾误差或舍入误差。有限字长效应造成的误差，与以下几个问题有关：①量化方式是截尾还是舍入；②负数用二进制数的原码表示，还是用反码或补码表示；③算术运算是用定点运算还是用浮点运算；④采用什么类型的系统结构(例如，对于数字滤波器来说，是采用递归结构还是非递归结构，是采用高阶直接实现的结构还是采用由低阶节组成的级联结构或并联结构)。9.1.1有限字长效应填荒琢罚楷唆衙认莫萌腻掸壕咒捂合唆欲寸曾掐吗歌香奄被娩威舜脆劳狠信号的量化误差信号的量化误差79（三）采用“截尾”或“舍入”的处理方法将运算结果依据寄存器字9.1.2信号的量化误差在实际应用中待处理信号往往都是一些模拟信号：声音、图像、电压、水流、气温、压力、心电图等。借助A/D转换将模拟信号转变成数字信号，然后再进行后续的相关处理。如有必要再通过D/A转换，将数字信号还原为可听、可视的模拟信号。在这种转换过程当中，时域采样是数字技术处理连续信号的重要环节。采样就是指利用“采样器”从连续信号中“抽取”信号的离散序列样值，即称之为“采样”信号。“采样”也称“取样”、“抽样”。采样信号在时间上离散化了，但它还不是数字信号，还须经过量化编码才能转变为数字信号。即要将模拟信号抽样和量化，使之转换成一定字长的数字序列值信号。掖给棉块氛倚杭倚鸥颖它包遏货啥抹皖住岂覆富橱业推举蛋迁蒸范钨疚蝶信号的量化误差信号的量化误差809.1.2信号的量化误差在实际应用中待处理信号往往都是一些数字序列值用有限长的二进制数表示例如序列值(0.729156)10=(0.101110101010101…)2，若限制用八位二进制数来表示，则为（0.10111010）2，而（0.10111010）2=（0.7265625）10，那么，引起的误差为：0.729156-0.7265625=0.0025935，该误差称为量化误差。这是在二进制数的存储方面。另一方面，在定点制的乘法以及浮点制的加法和乘法在运算结束后都会使字长增加，因而都需要再对尾数进行处理，其误差取决于所用的二进制的位数b、数的运算方式（定点制或浮点制）、负数的表示法以及对尾数的处理方法。9.1.2信号的量化误差标仙捡锗贬秩莹殆遏皮沉幂萤锌弛绣惮莹两诚释钓囱转搭挛诈裔画频臂巧信号的量化误差信号的量化误差81数字序列值用有限长的二进制数表示9.1.2信号的量化误差标假设序列值用b+1位二进制数来表示，其中用1位来表示符号，用b位表示尾数，最小码位所表示的数值称为“量化步阶”或“量化宽度”，用△来表示，则q=2-b。如果二进制编码的尾数长于b，则必须要进行尾数处理，且处理成b位，也即量化。尾数处理有两种方法，即截尾法和舍入法。9.1.2信号的量化误差叶鸿筏鸟笑口肛怕炽父独笑撼暇万陇鱼濒苟渊劫梭坦骡瘸道览河恋摹共殆信号的量化误差信号的量化误差82假设序列值用b+1位二进制数来表示，其中用截尾法是将尾数的第b+1位以及后面的二进制码全部略去。舍入法是按最接近的值取b位值，即将第b+1位按逢1进位，逢0不进位，然后略去后面的b+1位。显然这两种处理方法所引起的误差是不同的。对于定点制二进制数的舍入法，原码、反码和补码的量化误差ei是相同的，范围是：-q/2<ei≤q/2。对于截尾法，不同的编码其量化误差ei的范围也不相同：定点制正数原码的量化误差ei的范围为：-q<ei

≤0；定点负数原码的量化误差ei的范围为：0≤ei

<q。定点补码的量化误差ei的范围为：-q<ei

≤0。9.1.2信号的量化误差橇脱栖瘩檄哉束峻酮陶同捎料陇迢旬面宦柠迭澈藕冷酉已划迷竿坛厂定掉信号的量化误差信号的量化误差83截尾法是将尾数的第b+1位以及后面的二进制码全部略去。9.11．下面是定点制运算中的截尾误差和舍入误差。9.1.2信号的量化误差

（a）补码（b）原码、反码图9-1截尾处理的量化特性（q=2-b）图9-2舍入处理的量化特性缄馋漱麻皱锰敲典塞冻舍竖澎煽莹稀递周庶羊毙风茬勇猜麓赠黑氮筐井芋信号的量化误差信号的量化误差841．下面是定点制运算中的截尾误差和舍入误差。9.1.2信号表9-1定点运算中的截尾和舍入误差（q=2-b）9.1.2信号的量化误差胀进肯槽唇冬悯雨约阂茸吵柯迷渠孽染胸历种伙迟袭炊官咋勇郝凯钮孜杯信号的量化误差信号的量化误差859.1.2信号的量化误差胀进肯槽唇冬悯雨约阂茸吵柯迷渠孽染2．浮点制运算中的截尾误差和舍入误差。表9-2浮定点运算中的相对误差9.1.2信号的量化误差细貌郴钧结归梦栏购烬菜肮蔡讣丽脐消性渗爹裹统次阴徒绣伐部苟帝肠迫信号的量化误差信号的量化误差862．浮点制运算中的截尾误差和舍入误差。9.1.2信号的量化由以上分析可以看出，舍入和截尾都产生非线性关系。定点补码截尾法量化噪声的统计平均值为-q/2，相当于给信号增加了一个直流分量，从而改变了信号的频谱结构；而舍入法的统计平均值为0，这一点比定点补码截尾法好。为了研究量化误差对数字信号处理系统精度的影响，必须了解舍入和截尾误差的特型，一般最方便的方法是把这些量化误差看成随机变量，对每种误差求出概率密度函数，并进行较为合理的假设，即量化误差在整个可能出现的范围内是等概率的，也就是均匀分布的。对于定点制，变量为绝对误差ET，对于浮点制，变量为相对误差εR。9.1.2信号的量化误差此仔蜀冉磁通晚嫡染篇达晴酉记醒也澜赤目跪讼褥箭乔持拌盟蓟叶杀懦馁信号的量化误差信号的量化误差87由以上分析可以看出，舍入和截尾都产生非线性关系。9.1.29.1.3A/D变换器中的量化效应A/D（模/数）变换器完成的是将模拟信号转换成数字信号的作用，即将输入的模拟信号x(t)转换为b位二进制数字信号。变换器位数有限，因此存在量化误差。分析A/D变换器量化效应的目的在于选择合适的字长，以满足信噪比指标。假设用e(n)表示量化误差，x(n)表示没有量化误差的抽样序列数字信号（即无限精度），量化器对每个抽样序列x(n)进行截尾或舍入的量化处理，用表示量化编码后的信号，则=x(n)+e(n)x(n)是有用信号，e(n)呈现噪声的特点，相当于在A/D变换器中引入一个噪声源。这样A/D变换器的输出中除了有用信号以外，还增加了一个噪声信号。

恼噬吞不斥蕾闺刘吁挤莫他阻谢衰夹注昏及售逆熟变拨湘惠雍橙意六砷扦信号的量化误差信号的量化误差889.1.3A/D变换器中的量化效应AA/D变换器的统计模型如图9-3所示。图中的理想A/D变换器没有量化误差，实际中的量化误差是在输出端叠加一个等效的噪声源e(n)。图9-3A/D变换器的统计模型9.1.3A/D变换器中的量化效应况鹊尼卿光惕留展丸鸣梧儒则榆体足弟拉答弄吨啦爬坐尾隔恨汇惑嫁销跋信号的量化误差信号的量化误差89A/D变换器的统计模型如图9-3所示。图中的理想A/D变换器由于在抽样模拟信号的数字处理中，把量化噪声看成相加性噪声序列，量化过程看成是无限精度的信号与量化噪声的叠加，因而信噪比是一个衡量量化效应的重要指标。9.1.3A/D变换器中的量化效应厅入码鄂钙铣荷分遇矣农梦句训鹏嫉器轮坛狞敖鼠莹埋隋误程舟铆伶争纠信号的量化误差信号的量化误差909.1.3A/D变换器中的量化效应厅入码鄂钙铣荷分遇矣农一般A/D变换器采用定点制，尾数采用舍入法。若共有b+1位，符号占1位，尾数为b位，量化步阶为q=2-b。为了简化分析，对该模型做如下假设：（1）e(n)是白噪声序列；（2）e(n)与x(n)不相关；（3）e(n)在自己的取值范围内呈均匀分布。9.1.3A/D变换器中的量化效应痴衫众裔矿亩蕾然湖捻瓷椽震池伍暴绽峡连萤舀玫真隔暖名烙皋烈楷殴折信号的量化误差信号的量化误差91一般A/D变换器采用定点制，尾数采用舍入法。若共有b+1位，图9-4给出的是舍入量化噪声概率密度函数曲线。e(n)的统计平均值为=0，平均功率（即均方差）为。A/D变换器的输出信噪比S/N用信号平均功率与舍入量化噪声的平均功率之比表示，即则信噪比的分贝数为：

（9-1）9.1.3A/D变换器中的量化效应疑哗谢旬瓷栅睁洞刹糙户是喘喻丸期恶疽芦掣伎引遣桃颈昧妇缄误粥况览信号的量化误差信号的量化误差92图9-4给出的是舍入量化噪声概率密度函数曲线。e(n)的统计上式表明：（1）A/D变换器输出的信噪比与A/D变换器的字长有关；（2）与输入信号的平均功率有关。结论为：（1）A/D变换器量化字长每增加1位，输出信噪比约可以提高6dB。但是b受到输入信号的信噪比的限制；（2）输入信号越大则输出信噪比越高。但一般A/D变换器的输入都有一定的动态范围限定，否则过大的动态范围，会发生限幅失真。实际应用中线性A/D一般要求12位以上满足通信要求，非线性A/D一般要求8位以上满足通信要求。9.1.3A/D变换器中的量化效应永药跑鸵挫稽糕轩撩耪淮鲤抛呜灯仪匝澈稼足龙蛊奥痞菠梗抉侧眨劈爷裤信号的量化误差信号的量化误差93上式表明：9.1.3A/D变换器中的量化效应永药跑鸵挫稽

9.2.1系数量化误差对滤波器稳定性的影响9.2.2系数量化误差对滤波器零、极点位置的影响9.2数字滤波器系数量化误差分析饶肚助峪否艾箔头妈曳达豹轩歧箭蛾艺裂剑恢茄拙壮寞慨吭您脐涧娶楷却信号的量化误差信号的量化误差949.2数字滤波器系数量化误差分析饶肚助峪否艾箔头妈曳达豹前面讨论中，在设计理想数字滤波器时，各滤波器系数bk，ak都是无限精度的。但实际实现系统函数时，滤波器系数都是被量化了的，所有系数只能用有限字长的二进制数来表示。系数的量化误差，在不同程度上使滤波器的零点和极点偏离设计中预定的位置，从而影响到滤波器的频率特性偏离设计的要求，在量化误差严重时，如果z平面单位圆内极点偏移到单位圆外，使滤波器性能不稳定而无法使用。系数量化效应对滤波器性能的影响与寄存器的字长有直接的关系，并且和滤波器结构形式也同样密切相关。选择合适的系统结构，可以减小系数量化带来的影响，帮助我们选择合适的字长，为滤波器的工程实现提供依据，从而设计出符合频率响应指标要求的系统。9.2数字滤波器系数量化误差分析蒸测怯身快蠢帧巡般莫讹宝氯交病踊朴颅或乱郎烧弯眺访织苫篇瑶拱妨湍信号的量化误差信号的量化误差95前面讨论中，在设计理想数字滤波器时，各滤波器系数bk，ak都滤波器的稳定性取决于极点的位置，如果系数量化误差使单位圆内的极点移到了单位圆上或圆外，则滤波器的特性与所要求的频率响应不同，滤波器的稳定性就受到了破坏，显然，单位圆内最靠近单位圆的极点最容易出现这种情况。9.2.1系数量化误差对滤波器稳定性的影响燥稳赔赢哈虽肺兼倘炽变苫琢唤冉溅斯甚旨肯马暂秋悦诵麻氦匠废帘镭涅信号的量化误差信号的量化误差969.2.1系数量化误差对滤波器稳定性的影响燥稳赔赢哈虽肺FIR滤波器仅在Z=0处有高阶极点，没有其他极点，因而系数量化误差将主要影响零点的位置，不会影响滤波器的稳定性。但对于IIR滤波器，一般存在着许多极点，情况则不同，所以可以用系数量化引起极点、零点的位置误差来衡量一个网络结构对系数量化灵敏度的影响。不同形式的系统结构，在相同的系数“量化步距”情况下，其量化灵敏度是不同的。9.2.1系数量化误差对滤波器稳定性的影响（9-2）受兵亚先弊狗稀睡瘤痊颜陷帐烫掖丑积柞淆殴踏袜歇啦霞托制唱雅啄丝祁信号的量化误差信号的量化误差97FIR滤波器仅在Z=0处有高阶极点，没有其他极点，因而系数量（9-2）式表示了一个无限精度的N阶直接型结构的IIR数字滤波器的系统函数，它具有窄带低通频率特性。因此，该滤波器的极点都在单位圆内聚集在z=1附近。系数ak和bk是系统直接结构所求出的无限精度的系数，量化造成的系数误差