数字信号处理及Matlab实现课件作者李辉第7章节.ppt

7.6 FFT算法的有限字长效应 从运算的角度分析，FFT与数字滤波器一样都是系统，因此，二者的有限字长效应分析方法基本上是相同的 不同的是FFT中的运算是一个复数运算，因此误差源是一个复数，同时FFT运算有多个输入端及多个输出端，各误差源对不同的输出端的影响是不同的。 以时间抽取DIT为例进行有限字长效应分析，并且针对的是舍入情况，其他FFT算法及截尾运算结果是相似的 7.6.1 蝶形运算的统计模型 定点制运算中只有乘法引入量化误差，加法运算不引入误差 蝶形运算中只在乘系数 时引入一个误差源 是复数相乘，所产生误差 就是一个复数 一个复数相乘要由四个实数相乘来构成、每个实数相乘都将引入一个相应的误差。因而共有四个误差 假设是舍入量化，则 复数 的方差 为 分析乘以 后的方差影响 通过所有后级蝶形时，其方差保持不变 每一个输出端都与N－1个蝶形相连， 即有N－1个量化噪声源对每个输出端有贡献 在终端，在离散傅里叶变换 上叠加的输出噪声的均方值为 和FIR滤波器的直接型实现一样， 输出噪声的总方差正比于N 7.6.2 防止溢出和FFT输出的信噪比 分析蝶形运算 同样 所以 蝶形运算的输出最大值不超过输入端最大值的两倍。 从前一级到后一级，最大模值是逐级非减的，只要最后一级不出现溢出，则前一级计算一定不会溢出。 三种防止溢出的办法 1、输入端—次衰减法 由蝶形公式可以证明 为使 不溢出，即 ， 所以输入 应满足条件 若 在 内等概率分布，则 的方差为 输出信噪比为 由于 ，输出信号的方差为 信噪比与N2成反比，N增加一倍，SNR下降4倍 若保持运算精度不变，每增加一级运算，字长也需相应增加一位。这种防止溢出的办法，使得输入幅度被限制得过小，造成输出信号／噪声比值过小 2、逐级衰减法 对每个蝶形结的两个输入支路都乘上个比例因子 每—输入信号均乘以1/2，因此每个蝶形有二个乘法运算，即有二个误差源 、 。若每个误差源的方差为 则蝶形运算总的误差方差为 由于系数1/2存在，各误差源到输出端的传输系数不再是±1或 的连乘积， 而是 或 的连乘积，连乘次数与蝶形位置有关。 因此第m级蝶形的误差源在输出端的误差方差 为 总的输出方差为 输出信噪比为 与输入端一次衰减法相比，信噪比有了很大的提高 为保证运算精度不变，N增加4倍，字长只需增加一位 3、成组浮点运算 保证 ， 将原始数据用成组浮点制表示为 为某个整数，称为共阶数 对 作定点FFT，当某个蝶形结出现溢出，则将整个这一级的序列(运算过的，未运算的)全部左移一位(即乘1／2因子)，并在阶码 p 上加1，然后运算从发生溢出的蝶形继续下去。 当以后某级蝶形运算又出现溢出时，再对该级的输入乘上1/2衰减因子，并在阶码p上加1继续运算，直到L级计算完为止。 如果在整个变换过程中，一次溢出也没有发生，输出信号方差最大，当输入为白噪声型信号时，其值为 输出信噪比为 最坏的情况是各级都产生了衰减，这相当于第二种方法。所以成组浮点法快速傅立叶变换的信噪比的变化范围为 结论： 系数量化对零、极点位置的影响与零、极点位置的分布以及滤波器的结构均有密切的关系 高阶滤波器：避免采用直接型的结构，而应分解为最低阶的级联结构或并联结构 对于极点灵敏度很高的场合，可以来用双精度的系数，以便有效的达到精度的要求。 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