第五章有限字长运算对数字滤波器的影响-第三节课件.ppt
- 【下载声明】
1. 本站全部试题类文档,若标题没写含答案,则无答案;标题注明含答案的文档,主观题也可能无答案。请谨慎下单,一旦售出,不予退换。
2. 本站全部PPT文档均不含视频和音频,PPT中出现的音频或视频标识(或文字)仅表示流程,实际无音频或视频文件。请谨慎下单,一旦售出,不予退换。
3. 本页资料《第五章有限字长运算对数字滤波器的影响-第三节课件.ppt》由用户(晟晟文业)主动上传,其收益全归该用户。163文库仅提供信息存储空间,仅对该用户上传内容的表现方式做保护处理,对上传内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知163文库(点击联系客服),我们立即给予删除!
4. 请根据预览情况,自愿下载本文。本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
5. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007及以上版本和PDF阅读器,压缩文件请下载最新的WinRAR软件解压。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 第五 有限 运算 数字滤波器 影响 三节 课件
- 资源描述:
-
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.1 IIR 的有限字长效应的有限字长效应 以一阶IIR滤波器为例,其输入与输出关系可用差分方程表示为:乘积项将引入一个舍入噪声,如图 上述一阶系统的单位脉冲响应为 系统函数为 由于 是迭加在输入端的,故由 造成的输出误差为:)()1()(nxnayny1,0an)()(nuanhnazzzH)()(ne)(ne)(*)()(*)(nuanenhneenf1 输出噪声方差 或 由上两式均可求得 可见字长
3、 越大,输出噪声越小,同样的方法可分析其它高阶DF的输出噪声。0022222)(mmmeefamhcefzdzzHzHj)()(2122)1(122)1(1212222222aaqabefb例4:一个二阶IIR低通数字滤波器,系统函数为 采用定点制算法,尾数作舍入处理,分别计算其直接型、级联型、并联型三种结构的舍入误差。解:直接型 )8.01)(9.01(04.0)(11zzzH)(04.072.07.1104.0)(21zBzzzH直接型结构流图如图)(0ne)(nx)()(nenyf)(2ne)(1ne1z1z04.07.172.01211111()1 1.70.72(10.9)(10.8
4、)B zzzzz 图中 、分别为系数0.04、1.7、-0.72相乘后引入的舍入噪声。采用线性迭加的方法,从图上可看出输出噪声 这三个舍入噪声通过网络 形成的,如图5.28b,因此 是H0(z)的单位脉冲响应,输出噪声的方差为:)(0ne)(1ne)(2ne)(1)(0zBzH)(*)()()()(0210nhnenenenef)(0nh)(nefcefzdzzBzBj)()(1213122输出噪声的方差为:将 代入,利用留数定理得:1222qe2111132(10.9)(10.8)(10.9)(10.8)ecdzjzzzzz 220.90.82223(0.8)(1 0.9)(1 0.8)(0
5、.9)(1 0.9)(1 0.8)42.27619.8422.4fezzzzzzzzzzqqq 级联型 将H(z)分解 结构流图为 图5.29 IIR级联型的舍入噪声分析)(1)(04.08.0119.0104.0)(2111zBzBzzzH)(0ne)(nx)()(nenyf)(2ne)(1ne1z1z9.08.004.0由图中可见,噪声 、通过 网络,噪声 只通过网络 ,即 )(0ne)(1ne)(1zH101211()()()()()HzHzBz BzB z)(2ne)(2zH)(1)(22zBzH)(*)()(*)()()(22110nhnenhnenenef)(1nh)(2nh和 分
6、别是H1(z)和H2(z)的单位脉冲响应,因此:将 代入,得:(思考:如果将H1(z)和H2(z)次序颠倒,结果会怎样)cecefzdzzBzBjzdzzBzBzBzBj)()(12)()()()(12212221211212212,8.01)(,9.01)(221211qzzBzzBe2222222.40.2315.23fqqq因此:将 代入,得:221112122111212()()()()12()()efcecdzjB z B z B zB zzdzjB z B zz12,8.01)(,9.01)(221211qzzBzzBe22215.3715.2fqq2222222.40.4415.
7、373fqqq220.810.4412(10.9)zqqz1z1z0.80.904.0 并联型 将H(z)分解为部分分式 其结构如图:0.36 0.9 -0.32 0.8图 IIR并联型的舍入噪声分析 )(32.0)(36.08.0132.09.0136.0)(2111zBzBzzzH)(0ne)(nx)()(nenyf)(2ne)(1ne1z1z)(3ne 并联型结构有4个系数,有4个舍入噪声,其中 只通过 网络,通过 网络。输出噪声方差为:代入B1(z)和B2(z)及 的值,得:cecefzdzzBzBjzdzzBzBj)()(122)()(1221222111222e2222(0.440
8、.23)1.34fqq)()(10nene)(11zB)()(32nene)(12ZB比较三种结构的误差大小,可知 直接型 级联型 并联型原因:原因:l直接型结构的所有舍入误差都经过全部网络的反馈环节,反馈过程中误差积累,输出误差很大。l级联型结构,每个舍入误差只通过其后面的反馈环节,而不通过它前面的反馈环节,误差小于直接型。l并联型:每个并联网络的舍入误差只通过本身的反馈环节,与其它并联网络无关,积累作用最小,误差最小。该结论对IIR DF有普遍意义。因此,从有限字长效应看,直接型(、型)结构最差,运算误差最大,高阶时避免采用。级联型结构较好。并联型结构最好,运算误差最小。结论:IIR滤波器
9、的有限字长效应与它的结构有关。滤波器的有限字长效应与它的结构有关。5.3.2 FIR的有限字长效应的有限字长效应 IIR的分析方法同样适用于FIR滤波器,FIR滤波器无反馈环节(频率采样型结构除外),不会造成舍入误差的积累,舍入误差的影响比同阶IIR滤波器小,不会产生非线性振荡。以横截型结构为例分析FIR的有限字长效应。舍入噪声舍入噪声 N 阶FIR的系统函数为:无限精度下,直接型结构的差分方程为:有限精度运算时,10)()(NmmzmhzH10)()()(Nmmnxmhny10)()()()()(NmRfmnxmhnenyny每一次相乘后产生一个舍入噪声 故 输出噪声为:如图。)()()()
10、()(nemnxmhmnxmhmR1010)()()()()(NmNmmfnemnxmhneny10)()(Nmmfnene图中可见,所有舍入噪声都直接加在输出端,因此输出噪声是这些噪声的简单和。于是,输出噪声方差与字长有关,与阶数有关,N越高,运算误差越大,或者,在运算精度相同的情况下,阶数越高的滤波器需要的字长越长。12222NqNef例:FIR滤波器,N=10,b=17时N=1024时,因此,滤波器输出中,小数点后只有4位数字是有效的。)103(1085.41234210121122dbNqf)83(1097.412342102412922dbNqf-410*0.705f 动态范围:动态
11、范围:定点运算时,动态范围的限制,常导致FIR的输出结果发生溢出。利用比例因子,压缩信号的动态范围,可避免溢出。FIR输出:10max10)()()()()(NmNmmhxnymnxmhny定点数不产生溢出的条件:为使结果不溢出,对 采用标度因子A,使由此确定A。1)(ny)(nx10max10max)(11)(NmNmmhxAmhAx对于窄带信号略显保守5.3.3 FFT计算中的有限字长效应计算中的有限字长效应 按时间抽取的按时间抽取的8点点FFTM级、N/2个蝶形结1、输出噪声的方差:)()()()()()(1111nxWmxnxnxWmxmxlpNlllpNllxl-1(m)e(l,n)
12、xl(m)xl-1(n)PNW蝶形运算的量化误差-1xl(n),()()(ImReIm)(ReIm)(ImRe)(Re)(1413121111nleWnxenxWeWnxjeWnxeWnxWnxpNllpNpNlpNlpNlRpNl 相乘所引入的舍入误差源是复数,每个复乘包括4个实乘,每个定点实乘产生一个舍入误差源,因此产生四个误差源e1,e2,e3和e4,即)()(),(4321eejeenle2223EBqnle),(当误差源e(l,n)通过后级蝶形结时,加、减运算是不会影响方差的,而通过乘系数 后,对方差的影响为:pNW2222),(E),(E),(EnlenleWWnlepNpN e(
13、l,n)通过所有蝶形结时,加、减和乘的运算对方差均无影响。这样,计算FFT的最后输出误差,只需知道输出结点共连接多少个蝶形结即可。N=8时连接到X(0)的各蝶形结08W08W08W08W08W08W08W0Xx(0)x(4)x(2)x(6)x(1)x(5)x(3)x(7)-1-1-1-1-1-1-1连接到X(k)末端的误差源总数为 112222112NMM222)1(BkkNeEX(k)上叠加的输出噪声 的均方值(即为方差,因均值为零)为:ke3222NqNBk蝶形结输出的最大模值小于等于输入最大模值的两倍)(,)(max2)(,)(max)(,)(max1111nxmxnxmxnxmxlll
14、lll因此,FFT最后输出最大值有可能等于输入值的 倍,即NM2)(max)(max2)(maxnxNnxkXM因而要求 10,1)(NnNnx2、输出的信噪比这意味着为了防止溢出,可以在输入端一次性乘上比例因子1/N。如果假设x(n)在区间内 是均匀等概分布的,因而x(n)的方差为 NN1,122231)(NnxExNNWnxEkXExnkNNn31)()(221022而 因而 22222213)(E)SNR(qNNkXBk从而得到在这种防止溢出方法下输出的信噪比。第二种方法:引入了1/2的比例因子避免输出信噪比过小问题 加比例因子的蝶形运算WNp/21/2e(l,m)xl-1(m)e(l,
展开阅读全文