有限字长运算对数字滤波器的影响第三节

《有限字长运算对数字滤波器的影响第三节》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、§5.3有限字长运算对数字滤波器的影响DF的实现，涉及到两种运算：相乘、求和。定点制运算中，每一次乘法运算之后都要作一次舍入（截尾）处理，因此引入了非线性，采用统计分析的方法，将舍入误差作为独立噪声e(n)迭加在信号上，因而仍可用线性流图表示定点相乘。定点相乘运算统计分析的流图表示对舍入噪声e(n)作如下的假设：1.e(n)为平稳随机噪声序列；2.e(n)与输入序列x(n)不相关，各噪声之间也互不相关。3.e(n)为白色噪声；4.在量化间隔上均匀分布（即每个噪声都是均匀等概率分布）。有了这些条件，整个系统就可作为线性系统处理。每一个噪声可用第一章所讲的线性离散系统的理论求出其输出噪声，所

2、有输出噪声经线性迭加得到总的噪声输出。下面分析三种情况下的有限字长效应：IIR滤波器FIR滤波器FFT计算计算输出的量化噪声方差5.3.1IIR的有限字长效应以一阶IIR滤波器为例，其输入与输出关系可用差分方程表示为：乘积项将引入一个舍入噪声，如图上述一阶系统的单位脉冲响应为系统函数为由于是迭加在输入端的，故由造成的输出误差为：输出噪声方差或由上两式均可求得可见字长越大，输出噪声越小，同样的方法可分析其它高阶DF的输出噪声。例4：一个二阶IIR低通数字滤波器，系统函数为采用定点制算法，尾数作舍入处理，分别计算其直接型、级联型、并联型三种结构的舍入误差。解：①直接型直接型结构流图如图图中、

3、、分别为系数0.04、1.7、-0.72相乘后引入的舍入噪声。采用线性迭加的方法，从图上可看出输出噪声这三个舍入噪声通过网络形成的，如图5.28b，因此是H0(z)的单位脉冲响应,输出噪声的方差为：输出噪声的方差为：将代入，利用留数定理得：②级联型将H(z)分解结构流图为图5.29IIR级联型的舍入噪声分析由图中可见，噪声、通过网络，噪声只通过网络，即和分别是H1(z)和H2(z)的单位脉冲响应,因此：将代入，得：（思考：如果将H1(z)和H2(z)次序颠倒，结果会怎样）因此：将代入，得：③并联型将H(z)分解为部分分式其结构如图：0.360.9-0.320.8图IIR并联型的舍入噪声分

4、析并联型结构有4个系数，有4个舍入噪声，其中只通过网络，通过网络。输出噪声方差为：代入B1(z)和B2(z)及的值，得：比较三种结构的误差大小，可知直接型>级联型>并联型原因：l直接型结构的所有舍入误差都经过全部网络的反馈环节，反馈过程中误差积累，输出误差很大。l级联型结构，每个舍入误差只通过其后面的反馈环节，而不通过它前面的反馈环节，误差小于直接型。l并联型：每个并联网络的舍入误差只通过本身的反馈环节，与其它并联网络无关，积累作用最小，误差最小。该结论对IIRDF有普遍意义。因此，从有限字长效应看，直接型（Ⅰ、Ⅱ型）结构最差，运算误差最大，高阶时避免采用。级联型结构较好。并联型结构最好

5、，运算误差最小。结论：IIR滤波器的有限字长效应与它的结构有关。5.3.2FIR的有限字长效应IIR的分析方法同样适用于FIR滤波器，FIR滤波器无反馈环节（频率采样型结构除外），不会造成舍入误差的积累，舍入误差的影响比同阶IIR滤波器小，不会产生非线性振荡。以横截型结构为例分析FIR的有限字长效应。①舍入噪声N阶FIR的系统函数为：无限精度下，直接型结构的差分方程为：有限精度运算时，每一次相乘后产生一个舍入噪声故输出噪声为：如图。图中可见，所有舍入噪声都直接加在输出端，因此输出噪声是这些噪声的简单和。于是，输出噪声方差与字长有关，与阶数有关，N越高，运算误差越大，或者，在运算精度相同的

6、情况下，阶数越高的滤波器需要的字长越长。例：FIR滤波器，N=10，b=17时N=1024时，因此，滤波器输出中，小数点后只有4位数字是有效的。②动态范围：定点运算时，动态范围的限制，常导致FIR的输出结果发生溢出。利用比例因子，压缩信号的动态范围，可避免溢出。FIR输出：定点数不产生溢出的条件：为使结果不溢出，对采用标度因子A，使由此确定A。对于窄带信号略显保守5.3.3FFT计算中的有限字长效应按时间抽取的8点FFTM级、N/2个蝶形结1、输出噪声的方差：xl-1(m)e(l,n)xl(m)xl-1(n)蝶形运算的量化误差-1xl(n)相乘所引入的舍入误差源是复数，每个复乘包括4个实

7、乘，每个定点实乘产生一个舍入误差源，因此产生四个误差源e1，e2，e3和e4，即当误差源e(l,n)通过后级蝶形结时，加、减运算是不会影响方差的，而通过乘系数后，对方差的影响为：e(l,n)通过所有蝶形结时，加、减和乘的运算对方差均无影响。这样，计算FFT的最后输出误差，只需知道输出结点共连接多少个蝶形结即可。N=8时连接到X(0)的各蝶形结x(0)x(4)x(2)x(6)x(1)x(5)x(3)x(7)-1-1-1-1-1-1-1