数字信号处理及Matlab实现李辉第7章课件教学.pptVIP

PPT研究院 POWERPOINT ACADEMY * * 其次，为了保证不溢出，即要满足 如上图中的阴影部分． 由此式求出 、 的关系 若 必定期足 系统不产生溢出极限环振荡的充要条件： 由于溢出极限环的振荡幅度为1，数字滤波器无法对数据进行滤波，可以通过饱和型加法运算的办法加以消除。 饱和加法器的特点是，当相加器的输入之和大于1或小于－1时，就分别以1和－1代表相加结果，其特性如图所示，这样就能克服溢出振荡 补码饱和加法器特性 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级计算完为止。 如果在整个变换过程中，一次溢出也没有发生，输出信号方差最大，当输入为白噪声型信号时，其值为 输出信噪比为 最坏的情况是各级都产生了衰减，这相当于第二种方法。所以成组浮点法快速傅立叶变换的信噪比的变化范围为 馋死 研究当 时，仅仅由一个系数 的量化所引起的极点 的变化 若字长为8位，由量化引起的误差q/2可达大约0.002，求得 ，于是 极点远远超出单位圆，这样的变化显然是太大了。 字长至少要14位。 可见一个三阶系统对字长的要求已经非常严格了，如果阶数再高对量化误差的要求将更加苛刻。 如果采用定点二进制小数表示，则应满足： 结论： 系数量化对零、极点位置的影响与零、极点位置的分布以及滤波器的结构均有密切的关系 高阶滤波器：避免采用直接型的结构，而应分解为最低阶的级联结构或并联结构 对于极点灵敏度很高的场