数字信号处理及Matlab实现第7章.pptVIP

字长至少要14位。 可见一个三阶系统对字长的要求已经非常严格了，如果阶数再高对量化误差的要求将更加苛刻。 如果采用定点二进制小数表示，则应满足： 结论： 系数量化对零、极点位置的影响与零、极点位置的分布以及滤波器的结构均有密切的关系 高阶滤波器：避免采用直接型的结构，而应分解为最低阶的级联结构或并联结构 对于极点灵敏度很高的场合，可以来用双精度的系数，以便有效的达到精度的要求。 2.系数量化对二阶子系统极点位置的影响 高阶系统 级联型和并联型优于直接型。 级联型和并联型的基本子系统是二阶节 一个具有共轭极点的二阶系统有各种不同的结构 不同结构对于系数量化的敏感度也不同 二阶IIR系统的系统函数为 具有一对共轭对称的复极点 有 得到 共轭对称极点对组成的基本二阶网络的直接型实现 说明: 对于二阶网络，其极点的半径 完全由系数 决定，极点在实轴上的坐标值 取决于系数 。 如果 、 用三位字长表示，b＝3(不算符号位)， 只能表达8种半径 值和 之间的15种实轴坐标值 a1三位二进码β1,β2,β3 所表达的 的值 |a1 | 极点横坐标 a2三位二进码 所表达的 a2的值 极点半径 0.00 0.00 0.000 0.000 0.000 0.000 0.01 0.25 0.125 0.001 0.125 0.354 0.10 0.50 0.250 0.010 0.250 0.500 0.11 0.75 0.375 0.011 0.375 0.612 1.00 1.00 0.500 0.100 0.500 0.707 1.01 1.25 0.625 0.101 0.625 0.789 1.10 1.50 0.750 0.110 0.750 0.865 1.00 1.75 0.875 0.111 0.875 0.935 极点位置如右图中的网眼节点 如果所需要的极点位置不在这些网眼节点上时，就只能以最靠近的一个节点来代替这一极点位置，这样就引入了极点位置误差 系数量化使零极点位置的取值范围由一个连续域变为一个离散的平面点阵，从而造成零极点的漂移，导致系统特性的改变。 在平面上量化位置的分布密度是不均匀的 在实轴附近分布得稀，在虚轴附近分布得密； 半径小的地方分布得稀，半径大的地方分布得密。 这祥就会使实轴附近的极点(例加低通、高通滤波器)量化误差较大、而对虚轴附近的极点(例如带通滤波器)量化误差较小 这种分布只是二阶直接型结构的情况，不同结构的滤波器，系数对零极点位置的影响是不一样的 结论： 7.4.2 频率响应偏差的统计分析 1、IIR数字滤波器系数量化的统计分析 N阶IIR直接型结构为例，其理想精度的系统函数为 其中 是系数ai、bi的量化结果： 对系数ai、bi量化后，其实际传递函数为： 故系数量化后，实际的系统函数为 得到系统函数的偏差为 式中 系数量化造成的系统频率响应的偏差 系数量化的统计分析模型 2、FIR数字滤波器系数量化的统计分析 设FIR滤波器的系统函数为 对各系数 进行量化成为 有： 令 ， ，则 所以 FIR系数量化时，系统函数产生的误差不会超过 所以 当作定点舍入处理时，因为 7.5数字滤波器运算中的有限字长效应 实现数字滤波器所包含的基本运算有延时、乘系数和相加三种。因为延时不造成字长的变化，所以只讨论乘系数和相加运算造成的影响。 为了便于用统计方法分析这些量化误差的平均效应，我们假定 所有这些噪声都是平稳的白噪声序列； 所有噪声都与信号不相关，并且各噪声之间也互不相关 每个误差噪声都在其误差范围内呈均匀等概率分布 7.5.1 IIR滤波器定点运算舍入误差的统计分析 定点相乘运算的流图表示 采用统计分析方法后，实际的输出可以表示为 a理想相乘 b实际相乘的非线性流图 c统计模型的线性流图 方差 而每一个噪声源 所造成的输出噪声方差及均值： 最后将所有的输出噪声线性叠加就得到总的输出噪声 7.5.2 IIR滤波器定点加法运算的溢出问题 所有的相加运算，均可能产生溢出，应避免 使每个相加点的输入端都引入比例因子A加以限制，使相加点的数值绝对值保持小于1。 为使第k个相加节点上的输出信号 yk(n)不发生溢出，需要在输入端引入比例因子Ak。 要使yk(n)不发生溢出，则需使| yk(n) |＜1，所以 若 表示输入信号的最大幅度 因此 选择其中最小的比例