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

§5.3有限字长运算对数字滤波器的影响DF的实现，涉及到两种运算：相乘、求和。定点制运算中，每一次乘法运算之后都要作一次舍入（截尾）处理，因此引入了非线性，采用统计分析的方法，将舍入误差作为独立噪声e(n)迭加在信号上，因而仍可用线性流图表示定点相乘。定点相乘运算统计分析的流图表示对舍入噪声e(n)作如下的假设：1.e(n)为平稳随机噪声序列；2.e(n)与输入序列x(n)不相关，各噪声之间也互不相关。3.e(n)为白色噪声；4.在量化间隔上均匀分布（即每个噪声都是均匀等概率分布）。有了这些条件，整个系统就可作为线性系统处理。每一个噪声可用第一章所讲的线性离散系统的理论求出其输出噪声，所有输出噪声经线性迭加得到总的噪声输出。下面分析三种情况下的有限字长效应：IIR滤波器FIR滤波器FFT计算计算输出的量化噪声方差以一阶IIR滤波器为例，其输入与输出关系可用差分方程表示为：乘积项将引入一个舍入噪声，如图上述一阶系统的单位脉冲响应为系统函数为由于是迭加在输入端的，故由造成的输出误差为：321455.3.1IIR的有限字长效应输出噪声方差01或由上两式均可求得可见字长越大，输出噪声越小，同样的方法可分析其它高阶DF的输出噪声。02例4：一个二阶IIR低通数字滤波器，系统函数为采用定点制算法，尾数作舍入处理，分别计算其直接型、级联型、并联型三种结构的舍入误差。解：①直接型直接型结构流图如图图中、、分别为系数0.04、1.7、-0.72相乘后引入的舍入噪声。采用线性迭加的方法，从图上可看出输出噪声这三个舍入噪声通过网络形成的，如图5.28b，因此是H0(z)的单位脉冲响应,输出噪声的方差为：输出噪声的方差为：将代入，利用留数定理得：将H(z)分解贰级联型壹图5.29IIR级联型的舍入噪声分析肆结构流图为叁02和分别是H1(z)和H2(z)的单位脉冲响应,01由图中可见，噪声、通过网络，噪声只通过网络，即因此：将代入，得：（思考：如果将H1(z)和H2(z)次序颠倒，结果会怎样）010203因此：将代入，得：将H(z)分解为部分分式3690.328图IIR并联型的舍入噪声分析③并联型其结构如图：只通过网络，02输出噪声方差为：04并联型结构有4个系数，有4个舍入噪声，其中01通过网络。03代入B1(z)和B2(z)及的值，得：05比较三种结构的误差大小，可知直接型级联型并联型原因：l直接型结构的所有舍入误差都经过全部网络的反馈环节，反馈过程中误差积累，输出误差很大。l级联型结构，每个舍入误差只通过其后面的反馈环节，而不通过它前面的反馈环节，误差小于直接型。l并联型：每个并联网络的舍入误差只通过本身的反馈环节，与其它并联网络无关，积累作用最小，误差最小。该结论对IIRDF有普遍意义。因此，从有限字长效应看，直接型（Ⅰ、Ⅱ型）结构最差，运算误差最大，高阶时避免采用。级联型结构较好。并联型结构最好，运算误差最小。结论：IIR滤波器的有限字长效应与它的结构有关。5.3.2FIR的有限字长效应IIR的分析方法同样适用于FIR滤波器，FIR滤波器无反馈环节（频率采样型结构除外），不会造成舍入误差的积累，舍入误差的影响比同阶IIR滤波器小，不会产生非线性振荡。以横截型结构为例分析FIR的有限字长效应。2019舍入噪声012020N阶FIR的系统函数为：022021无限精度下，直接型结构