第7章 数字信号处理中的有效字长效应.

Slides:



Advertisements
Similar presentations
九族文化村兩天一夜遊 組員 : 傅淳鈺 9A0E0019 黃湘蓉 4A 陳誌龍 9A0K0026 潘韋舜 9A0B0951 何奇龍 4A
Advertisements

财务管理 利 润 分 配 利 润 分 配 嘉善中专 杨晓燕. 二、利润分配的项目及顺序 第三节 利润分配 一、利润分配的原则 财务管理 >> 第六章 >> 第三节 三、利润分配政策及影响因素.
不知者無罪嗎 ? 【本報台北訊】國內知名大學胡姓研究 生進口豬籠草在網路上販售,涉嫌違反 植物防疫檢疫法,胡姓研究生表示不知 道豬籠草是違禁品並當場認錯道歉 台北地檢署檢察官念他初犯,昨 天處分緩起訴,但命他繳交六萬 元緩起訴處分金作公益。 豬籠草有潛移性線蟲寄生,一旦植物感 染後,輕則枯萎凋零,重則危害農業經.
一、音调  听过女高音和男低音的歌唱吗?他们的声音 给你的印象是怎样的? 女高音:音调高, 男低音:音调低,比较低沉。
認識食品標示 東吳大學衛生保健組製作.
                                                 伊朗 的今生 与前世 (2)
第四章:长期股权投资 长期股权投资效果 1、控制:50%以上 有权决定对方财务和经营.
第 3 章 聽覺 每章扉頁 3.1 聽覺 3.2 人類耳朵的反應 3.3 聽覺缺陷.
第八章 互换的运用.
新約研讀 彼得前書複習 讀經組
颞下颌关节常见病.
「健康飲食在校園」運動 2008小學校長高峰會 講題:健康飲食政策個案分享 講者:啟基學校-莫鳳儀校長 日期:二零零八年五月六日(星期二)
客家文化的內涵與傳播 潘朝陽 臺灣師大國際與僑教學院院長 臺灣師大東亞系、地理系教授 臺灣師大全球客家文化研究中心主任
第一节 两者之间的差异分析 第二节 总体内部的差异分析 第三节 计算器的使用
南京市国税局国际税务管理处 二00九年二月二十四日
─視覺藝術的元素.
致理科技大學保險金融管理系 實習月開幕暨頒獎典禮
第四章 從分裂到統一 第一節 漢唐之際的大變動
第四章 從分裂到統一 第一節 漢唐之際的大變動
脊柱损伤固定搬运术 无锡市急救中心 林长春.
第一部分 中考基础复习 第一章 声现象.
2013年二手车市场环境分析.
結腸直腸腫瘤的認知.
經歷復活的愛 約翰福音廿一1-23.
往下扎根 向上生长 ——提高学考选考教学实效的几点感想 桐庐中学 李文娟
郭詩韻老師 (浸信會呂明才小學音樂科科主任)
Part I 上海土地市场.
湖北省,简称“鄂”,为中华人民共和国省级行政区。湖北在中国中部、长江中游、洞庭湖以北,介于北纬29°05′至33°20′,东经108°21′至116°07′;北接河南省,东连安徽省,东南和南邻江西、湖南两省,西靠重庆市,西北与陕西省为邻。东西长约740公里,南北宽约470公里,面积18.59万平方公里,占全国总面积的1.95%,居全国第13位。省会是中部地区唯一的副省级城市--武汉市。
现代社会生活中的压力症,是人们身心疾患 发生的根源。在学习企业管理培训课程的时候, 明白了当人们遇上"压力"时,最初的反应便是"
「但圣灵降临在你们身上,你们就必得着能力,
1890年, 一艘名叫“马尔波罗号”的帆船在从新西兰驶往英国的途中,突然神秘地失踪了。 20年后,人们在火地岛海岸边发现了它。奇怪的是:船体原封未动,完好如初;船长航海日记的字迹仍然依稀可辨;就连那些死去多年的船员,也都“各在其位”,保持着当年在岗时的“姿势”; 1948年,一艘名为“乌兰格梅奇号”的荷兰货船,在通过马六甲海峡时,突然遇到海上风暴,当救助人员赶到时,船上所有人员都莫明其妙地死了。
拟动力试验 伪动力试验,计算机加载器联机试验 地震发生和传播的随机性 周期性加载的加载历程是假定的,与实际地震的非周期反应有很大差别
第一章 绪论.
A.李时珍 B.司马迁 C.达尔文 D.袁隆平 D
新时代的劳动者 杜蒙绮.
平湖市当湖高级中学 平湖市教育局教研室 (电话)
2. 戰後的經濟重建與復興 A. 經濟重建的步驟與措施 1.
第四章 时间序列的分析 本章教学目的:①了解从数量方面研究社会经济现象发展变化过程和发展趋势是统计分析的一种重要方法;②掌握时间数列编制的基本要求;③理解和掌握水平速度两方面指标的计算及运用④理解和掌握长期趋势分析和预测的方法。 本章教学重点:现象发展的水平指标和速度指标。 本章教学难点:现象变动的趋势分析。
信号处理与系统课程教学案例 FFT的应用—— 声音信号合成与处理 国防科技大学电子科学与工程学院.
三级综合医院评审解读-生物安全 安徽医科大学第一附属医院检验科 徐元宏.
好好學習 標點符號 (一) 保良局朱正賢小學上午校.
“08高考化学学业水平(必修科目)测试的命题和教学对策研究”
提升國小自然與生活科技領域教師教學智能研習
第1节 光的干涉 (第2课时).
第三章 无限长单位脉冲响应(IIR)滤波器设计
第4章 种群和群落 第3节 群落的结构 自主学习案   合作探究案 课后练习案. 第4章 种群和群落 第3节 群落的结构 自主学习案   合作探究案 课后练习案.
温故知新 1、凸透镜成像的规律有哪些? 2、照相机成像的原理是什么?.
新約概論 台中生命之道靈糧堂 2007年3月4日.
早期的阿拉伯半島 地理環境: 生活情形 (一)三面環海,大多為荒涼貧瘠的沙漠,不利農耕
4. 聯合國在解決國際衝突中扮演的角色 C. 聯合國解決國際衝突的個案研究.
新陸書局股份有限公司 發行 第十九章 稅捐稽徵法 稅務法規-理論與應用 楊葉承、宋秀玲編著 稅捐稽徵程序.
狂賀!妝品系同學美容乙級通過 妝品系三甲 學號 姓名 AB 陳柔諺 AB 陳思妤 AB 張蔡婷安
民法第四章:權利主體 法人 楊智傑.
数字信号处理 Digital Signal Processing(DSP)
如何寫工程計畫書 臺北市童軍會考驗委員會 高級考驗營 版.
数字信号处理 by Zaiyue Yang CSE, ZJU, 2012.
四年級 中 文 科.
青眼究極龍 之 賓果連線 簡豪天、宋華敏製作.
聖誕禮物 歌羅西書 2:6-7.
数字信号处理基础 第7章 FIR数字滤波器的理论和设计
電子白板百萬小學堂 本活動建議搭配電子白板 學生最多可分成2~6組(請按組別按鈕) 老師可以視時間多少,來進行活動 每一組要回答十個問題。
藝術大師-達利.
「但圣灵降临在你们身上,你们就必得着能力,
介入及追蹤紀錄表 編號: 姓/稱謂: 初次103年 月 日 追蹤 月 日 問題型態 (可複選) □ 1. 覺得西藥都很傷胃
我們這班與梵谷有約 授課老師:李美娟 2000/9/18.
依撒意亞先知書 第一依撒意亞 公元前 740 – 700 (1 – 39 章) 天主是宇宙主宰,揀選以民立約,可惜他們犯罪遭
第七章 FIR数字滤波器的设计方法 IIR数字滤波器: FIR数字滤波器: 可以利用模拟滤波器设计 但相位非线性
「但聖靈降臨在你們身上,你們就必得著能力,
第三章时 域 分 析 引言 语音信号的短时处理方法 短时能量和短时平均幅度 短时平均过零率 短时自相关函数 短时时域处理技术应用举例
經文 : 創世紀一章1~2,26~28 創世紀二章7,三章6~9 主講 : 周淑慧牧師
Presentation transcript:

第7章 数字信号处理中的有效字长效应

前言 数字信号处理的实质:一组数值运算。 从设计的角度来讨论:认为数字是无限精度的。 从实现的角度考虑:数字的精度就是有限的 数字系统中的每一个数总是用有限字长的二进制数码表示,运算过程中需要的数字信号的值、系统的系数和运算过程中的结果都是存储在有限字长的存储单元中的,此时数字的精度就是有限的。 从设计时的无限精度到实现时的有限精度,会产生相对于原设计系统的误差,严重时会导致系统崩溃。

在数字系统中三种因有限字长的影响而引起误差的因素 A/D变换的量化效应 系数的量化效应 数字运算过程中的有限字长效应

7.1 二进制数的表示与量化误差 7.1.1 二进制数的表示 1、定点数的表示 2、浮点数的表示 M是数的尾数部分,R是数的基数,P是阶 所有数据小数点的位置是固定不变的 2、浮点数的表示 浮点数的小数点位置是不固定的,它随每个数的大小而变化 M是数的尾数部分,R是数的基数,P是阶 基数确定后,浮点数就完全由尾数和阶决定

一个浮点数可以表示如下: :阶P的符号位; :阶码,阶P的绝对值部分 :尾数M或浮点数N的符号位,尾符或数符; :尾数的绝对值部分,尾码。

规格化浮点数的小数点在数符Sm的后面,且小数点前有一个隐含“1”。其阶码连同阶符统一编码,浮点数的基数为2。 3、IEEE浮点数 规格化浮点数的小数点在数符Sm的后面,且小数点前有一个隐含“1”。其阶码连同阶符统一编码,浮点数的基数为2。 IEEE浮点数格式分单精度和双精度两种,单精度数为32位,双精度数是64位。  

单精度浮点数格式: 双精度浮点数格式:

表达式(a0·a1 a2 a3 … ab)表示的十进制数x为 带符号数的表示 设任意数x的(b十1)位码的形式为 a0·a1 a2 a3 … ab a0表示符号位 1、原码: a0是符号位, a1… ab是小数的绝对值 表达式(a0·a1 a2 a3 … ab)表示的十进制数x为

(2)反码 正数的反码表示与其原码相同,即 [x]反=x (x>0) 负数的反码由其绝对值按位求反后得到  [x]反=2-[x]原 表达式(a0·a1 a2 a3 … ab)表示的十进制数x为

(3)补码 正数的补码表示与其原码相同,即 [x]补=x (x>0) 负数的补码由它的绝对值求反加1后得到: 表达式(a0·a1 a2 a3 … ab)表示的十进制数x为

7.1.2 量化误差 一、定点运算中的截尾误差和舍入误差 1、截尾误差 7.1.2 量化误差 一、定点运算中的截尾误差和舍入误差 1、截尾误差 ①对于正小数x≥0:原码、反码、补码的表示法相同,因而量化影响也相同。 截尾前x有b1位,有 截尾后x有b位,记做 ,有 以ET表示截尾误差,则有

正小数截尾后数值变小,故截尾误差总是负的。当被截位ai (i=b+1到i=b1)均为l时,为最大截尾误差 令 ——量化间距或量化步阶,因此下式成立

②对于负小数x<0,由于a0≠0,不同码制x的表示法各不相同,因而产生的量化误差也不相同。 i)原码:截尾后负小数的绝对值变小,截尾误差为正:

ii)对于补码,

iii)对于反码:

结论: 原码与反码的截尾误差与数的正负有关: 补码的截尾误差皆为负数 正数时误差为负,当x≥0时,-q<ET≤0 负数时误差为正,当x<0时,0≤ET <q 补码的截尾误差皆为负数 对所有的x,-q<ET≤0。

定点制截尾处理的量化特性(q=2-b) (a)补码 (b)原码、反码

2、舍入误差 舍入误差总是处在+q/2之间, 用 表示对x作舍入处理, ER表示舍入误差,有: 定点制舍入处理的量化特性

定点运算中的截尾和舍入误差(q=2-b) 截尾误差 舍入误差 正数 -q <ET≤0 -q /2≤ER≤q /2 负数 原码     截尾误差 舍入误差 正数 -q <ET≤0 -q /2≤ER≤q /2 负数 原码 0≤ET < q 反码 补码

二、浮点运算中的截尾误差和舍入误差 浮点数x=2C·M 的相对误差定义为: 绝对误差可表示成

1、浮点舍入 x=2C·M 中尾数M的误差范围: 而x的阶码为c,则

2、浮点截尾 ①当x>0时,a0=0,三种码制的截尾误差均为-q<ET≤0

2、浮点截尾 ①当x<0时,a0=1,原码和反码的截尾误差为0≤ET<q 对于补码,截尾误差为-q<ET≤0,相应的,可得

  浮点运算中的相对误差(q=2-b)   截尾 舍入 正数 -q<ET≤0 负数 原码 反码 补码

7.2 A/D转换器中的量化误差 A/D转换器 A/D变换器前一般都加一个前置模拟低通滤波器 模拟输入信号必须乘一个比例因子 将输入的模拟信号xa(t)转换为b位二进制数字信号的器件。 b的数值可以是8,12或高至20。 A/D变换器前一般都加一个前置模拟低通滤波器 滤除高于折叠频率(抽样频率之半)的频率, 模拟输入信号必须乘一个比例因子 A/D变换器总是定点制的,必须使信号不超过A/D变换的动态范围

A/D转换器包括 抽样 量化 产生抽样序列x(n)=xa(t)|t=nT=xa(nT), x(n): 可看成是一个无限精度的数字信号 在满足抽样定理的前提下,模拟信号时间离散化的过程是可逆的。 量化 对抽样序列进行幅度上的离散化之后,用某种格式的数字代码来表示。 量化过程是不可逆的 必定要引入量化误差或量化噪声。 量化噪声的大小决定了A/D转换器的动态范围,是恒量A/D转换器性能的一个最重要指标。

7.2.1 量化误差的统计分析 设量化器的输入信号为随机序列x(n),其量化误差也为随机序列e(n): e(n)=Q[x(n)]-x(n) 7.2.1 量化误差的统计分析 设量化器的输入信号为随机序列x(n),其量化误差也为随机序列e(n): e(n)=Q[x(n)]-x(n) 假设e(n)具有下列特性: e(n)是一个平稳随机序列 e(n)与信号也不相关 e(n)本身的任意两个值之间不相关,具有白噪声性质 e(n)在其误差范围内均匀等概分布。

A/D转换原理图 A/D转换的统计模型

:代表了除去直流分量后量化噪声的平均功率 e(n)的两个最重要的统计参数 均值me :代表噪声的直流分量 :代表了除去直流分量后量化噪声的平均功率 其中 :E[·]表示取数字期望 p(e)是误差值e(n)的概率密度

三种误差范围的概率密度函数

1、舍入处理 舍入误差:-q /2≤e(n)≤q /2 e(n)的概率分布密度为1/q 可得均值me 及方差

2、截尾处理 对于x≥0的三种码制和x<0的补码 : 可得均值me 及方差 e(n)的误差为-q <e(n)≤0

对于x<0的原码和反码: e(n)的误差为0≤ET <q 概率分布密度为1/q 可得均值me 及方差

结论: 各种情况的方差 均为 , 不同的只是均值me,分别为 。 量化噪声的方差 与量化间隔的平方 成正比, :量化间隔, 信号处理时选用的字长b+1越长,量化噪声的方差越小。 另:由于截尾噪声具有直流分量,将影响信号的频谱结构,因此一般采用舍入处理。

7.2.2 量化信噪比与所需字长的关系 量化的信噪比 信号的平均功率 与量化噪声的平均功率 之比 信噪比用对数表示时记作SNR,单位dB 7.2.2 量化信噪比与所需字长的关系 量化的信噪比 信号的平均功率 与量化噪声的平均功率 之比 信噪比用对数表示时记作SNR,单位dB 信号功率 越大,信噪比越高; 随着字长b的增加,信噪比也增大,字长b每增加一位,则信噪比增加约6dB。字长越长,A/D变换的信噪比越高

为了使信号不超过定点制运算所允许的动态范围,用一个小于1的正数A(0<A<1=去乘x(n) 信号的平均功率 与量化噪声的平均功率 之比 选择A值:使

若需得到信噪比大于70dB,至少需要满足b=12bit 字长过长也无必要: 输入信号xa(t)本身也有一定的信噪比,字长长到A/D变换器的量化噪声比xa(t)的噪声电平更低则没有意义 为提高信噪比, 可以增大输入信号,但这受到A/D变换器动态范围的限制, 还可以增加字长b,但这又受到输入信号xa(t)信噪比的限制

7.2.3 量化噪声通过线性非时变系统 量化噪声通过线性系统

设e(n)是定点补码舍入误差,e(n)的均值为me 、方差为 则系统量化噪声的输出 的均值mf和方差 计算如下:

根据Parseval定理, 也可以用下式表示: 或者在单位圆上计算 如果e(n)是补码截尾白噪声: 输出噪声的方差 不变,输出的均值mf如下:

7.3 系数量化对数字滤波器的影响 7.3.1 极点位置灵敏度 每个极点位置对各系数偏差的敏感程度 决定系统的稳定性 对系统性能的影响亦较大

理想数字滤波器的系统函数为 对系数ai、bi量化后,其实际传递函数为: 其中 是系数ai、bi的量化结果:

下面讨论系数量化误差对极点的影响: 原系统函数H(z)的分母多项式为 令A(z)=0,得到H(z)的极点: 系数量化后的极点为:

:极点位置的偏差量,由各个系数偏差 引起的,因此 结论: 就是极点 对系数 变化的灵敏度 越大, 对 的影响也越大 越小, 对 的影响就越小

下面根据A(z)来求这个极点位置灵敏度 的表达式: 而:

分母中的每一个因子 代表着某一极点 指向当前极点 的矢量,而整个分母正是所有极点指向当前极点 的矢量积 得极点位置灵敏度为: 分母中的每一个因子 代表着某一极点 指向当前极点 的矢量,而整个分母正是所有极点指向当前极点 的矢量积 高阶直接型结构滤波器对系数量化误差的敏感性高于低阶直接型结构滤波器 并联型结构及级联型结构每对极点只受与之有关的两个系数的影响,且每个子系统的极点密集度要稀疏得多,因而极点位置受系数量化的影响比直接型结构要小得多

极点位置灵敏度与极点间距离成反比

例1 设一低通滤波器的传递函数如下,分析计算系数量化对极点位置的影响。 解 经计算求得 的极点分别是 于是:

研究当 时,仅仅由一个系数 的量化所引起的极点 的变化 若字长为8位,由量化引起的误差q/2可达大约0.002,求得 ,于是 极点远远超出单位圆,这样的变化显然是太大了。

再研究 ,即将极点移到单位圆上需要的字长 如果 =2.945的量化误差等于 ,就会使量化后系统的极点移到单位圆上 二进制数中为2-14< <2-13,字长至少要14位。 可见一个三阶系统对字长的要求已经非常严格了,如果阶数再高对量化误差的要求将更加苛刻。

考察用三个一阶的环节级联或并联组成这个系统。 每一个环节中极点从0.99变到1.00,允许变化0.01 而且这个环节的稳定性并不受另一环节影响 所以其所需字长为7位 结论: 系数量化对零、极点位置的影响与零、极点位置的分布以及滤波器的结构均有密切的关系 高阶滤波器:避免采用直接型的结构,而应分解为最低阶的级联结构或并联结构 对于极点灵敏度很高的场合,可以来用双精度的系数,以便有效的达到精度的要求。

7.3.2 系数量化对二阶子系统 极点位置的影响 高阶系统 级联型和并联型的基本子系统是二阶节 级联型和并联型优于直接型。 7.3.2 系数量化对二阶子系统 极点位置的影响 高阶系统 级联型和并联型优于直接型。 级联型和并联型的基本子系统是二阶节 一个具有共轭极点的二阶系统有各种不同的结构 不同结构对于系数量化的敏感度也不同

二阶IIR系统的系统函数为 具有一对共轭对称的复极点 有 得到

共轭对称极点对组成的基本二阶网络的 直接型实现

说明: 对于二阶网络,其极点的半径 完全由系数 决定,极点在实轴上的坐标值 取决于系数 。 如果 、 用三位字长表示,b=3(不算符号位), 只能表达8种半径 值和 之间的15种实轴坐标值

三位二进码β1,β2,β3 所表达的 的值 |a1 | 极点横坐标 a2三位二进码 a2的值 极点半径 0.00 0.000 0.01 0.25 0.125 0.001 0.354 0.10 0.50 0.250 0.010 0.500 0.11 0.75 0.375 0.011 0.612 1.00 0.100 0.707 1.01 1.25 0.625 0.101 0.789 1.10 1.50 0.750 0.110 0.865 1.75 0.875 0.111 0.935

极点位置如下图中的网眼节点 如果所需要的极点位置不在这些网眼节点上时,就只能以最靠近的一个节点来代替这一极点位置,这样就引入了极点位置误差

结论: 系数量化使零极点位置的取值范围由一个连续域变为一个离散的平面点阵,从而造成零极点的漂移,导致系统特性的改变。 在平面上量化位置的分布密度是不均匀的 在实轴附近分布得稀,在虚轴附近分布得密; 半径小的地方分布得稀,半径大的地方分布得密。 这祥就会使实轴附近的极点(例加低通、高通滤波器)量化误差较大、而对虚轴附近的极点(例如带通滤波器)量化误差较小 这种分布只是二阶直接型结构的情况,不同结构的滤波器,系数对零极点位置的影响是不一样的

基本二阶网络的另一种实现 该网络的系统函数为 : 当系数量化时,是对 及 进行量化,因而所得到的网格点子在z平面是均匀分布的 这里系数量化后对z平面的所有区域,所产生的误差是相同的。

具有共轭极点对的二阶数字网络耦合形式实现情况下,系数量化为三比特时极点的可能位置

7.3.3 频率响应偏差的统计分析 1、IIR数字滤波器系数量化的统计分析 N阶IIR直接型结构为例,其理想精度的系统函数为 7.3.3 频率响应偏差的统计分析 1、IIR数字滤波器系数量化的统计分析 N阶IIR直接型结构为例,其理想精度的系统函数为 对系数ai、bi量化后,其实际传递函数为: 其中 是系数ai、bi的量化结果:

故系数量化后,实际的系统函数为

得到系统函数的偏差为 式中 系数量化造成的系统频率响应的偏差

系数量化的统计分析模型

可以用频响的均方偏差来描述 系数量化所引起的频率特性偏差 均方偏差 也是一个随机变量.它的均值即为频响偏差

2、FIR数字滤波器系数量化的统计分析 设FIR滤波器的系统函数为 对各系数 进行量化成为 有:

令 , ,则 所以 FIR系数量化时,系统函数产生的误差不会超过 当作定点舍入处理时,因为 所以

7.4 数字滤波器的运算量化效应 为了便于用统计方法分析这些量化误差的平均效应,我们假定 所有这些噪声都是平稳的白噪声序列; 7.4 数字滤波器的运算量化效应 为了便于用统计方法分析这些量化误差的平均效应,我们假定 所有这些噪声都是平稳的白噪声序列; 所有噪声都与信号不相关,并且各噪声之间也互不相关 每个误差噪声都在其误差范围内呈均匀等概率分布

7.4.1 IIR滤波器 定点运算舍入误差的统计分析 定点相乘运算的流图表示 a理想相乘 b实际相乘的非线性流图 c统计模型的线性流图 采用统计分析方法后,实际的输出可以表示为 定点相乘运算的流图表示

7.4.1 IIR滤波器 定点运算舍入误差的统计分析 定点相乘运算的流图表示 采用统计分析方法后,实际的输出可以表示为 a理想相乘 b实际相乘的非线性流图 c统计模型的线性流图 定点相乘运算的流图表示 采用统计分析方法后,实际的输出可以表示为

而每一个噪声源 所造成的输出噪声方差及均值: 最后将所有的输出噪声线性叠加就得到总的输出噪声 方差

重要结论:有限字长效应与滤波器的结构型式有密切关系。 例 2 已知一个IIR滤波器的系统函数如下,用定点制算法,尾数舍入,分别求出直接型、级联型和并联型实现系统时量化误差的方差 。 。 、 1、直接型结构

直接型的舍入噪声 b 三个舍入噪声 通过相同的传输网络 a 相乘引入的舍入噪声 :系数0.2、1.7以及—0.72相乘后的舍入噪声 它们均经过相同的传输网络

输出噪声的方差是

2、级联型:可以有几种排列形式 先令 ,并把0.2置于第一级,即

、 通过网络 、 通过网络 输出噪声的方差是

第一级有两个误差源,第二级有一个误差源,故有 所以

再令 ,并把0.2置于第二级,有

、 通过网络 、 通过网络 输出噪声的方差

第一级有两个误差源,第二级有一个误差源,故有 所以 结论:系数bi和A(k)的排列不同,输出误差也有所不同。 距单位圆远的极点和系数bi置于较后级的误差会小些。

3、并联型

并联型结构需要4个系数,因此共有四个舍入噪声 、 通过网络 、 通过网络 因此输出的方差为

则第一级、第二级都有两个误差源,故有 所以 由此得出结论:对IIR滤波器,从有限字长效应来看,不论是哪一种型式的直接型结构都是最差的,运算误差最大,特别在高阶时应避免采用。级联型结构较好,而并联型结构具有最小的运算误差。

7.4.2 IIR滤波器定点加法运算的溢出问题 所有的相加运算,均可能产生溢出,应避免 使每个相加点的输入端都引入比例因子A加以限制,使相加点的数值绝对值保持小于1。 为使第k个相加节点上的输出信号 yk(n)不发生溢出,需要在输入端引入比例因子Ak。

若 表示输入信号的最大幅度 因此 要使yk(n)不发生溢出,则需使| yk(n) |<1,所以 选择其中最小的比例因子作为系统比例因子的最终选择,即

令AK=1,得到保证第k个节点上不出现溢出时的最大输入值上限

7.4.3 极限环振荡 IIR滤波器在无限精度的情况下,当它的所有极点均位于单位圆内时,系统肯定是稳定的。 7.4.3 极限环振荡 IIR滤波器在无限精度的情况下,当它的所有极点均位于单位圆内时,系统肯定是稳定的。 当去掉输入信号后.随着n的增加,系统输出逐渐衰减趋向于零。 在有限字长情况下,由于量化过程的非线性作用,系统输出将不随n的增加而趋于零 零输入极限环振荡:系统输出不随n的增加而趋于零,而是衰减到某一非零的范围幅度后呈现振荡特性 以舍入处理的一阶IIR滤波器为例来说明这一现象

设一阶IIR系统的系统函数为 在无限精度运算下,其差分方程为 在定点运算中,每次乘法运算以后均要对尾数进行舍入处理。因此,实际的非线性差分方程可表示为

一阶IIR网络的非线性流图 在无限精度情况下,如果输入信号 输出 只要|a|<1,系统稳定,y(n)将逐渐衰减为零。

有限精度运算时,系数a=0.5=0.100b的情况 n 0.111 0.000 0.000000 0.875 1 0.011100 0.100 0.500 2 0.010000 0.010 0.250 3 0.001000 0.001 0.125 4 0.000100 ┇

零输入极限环振荡 以下也称为“死带”区域

下面分析振荡幅度和字长的关系 由于舍入误差的绝对值在 以内,因此 或者 可得 一阶IIR网络的死带范围 表明:极限环幅度与量化间隔成正比, 增加字长(减小量化间隔)将使极限环振荡减弱 如b=3,故有 , ,有 与表7-4-1结果一致

7.4.4 定点运算的溢出振荡 溢出可能产生比之更大的误差,并在输出的最大幅度界限内振荡,称为溢出极限环振荡 7.4.4 定点运算的溢出振荡 溢出可能产生比之更大的误差,并在输出的最大幅度界限内振荡,称为溢出极限环振荡 溢出极限环振荡产生的原因:定点加法运算存在着溢出 补码加法运算的非线性

以定点补码运算的二阶IIR滤波器为例进行讨论 此二阶滤波器的差分方程为 系统函数为

首先考虑系统的稳定性 稳定的条件是 的两个极点 、 在单位圆内 可得 、 的值域为图示的大三角形 即 、 落于大三角形之内才能满足

其次,为了保证不溢出,即要满足 由此式求出 、 的关系 若 系统不产生溢出极限环振荡的充要条件 必定期足 如上图中的阴影部分.

由于溢出极限环的振荡幅度为1,数字滤波器无法对数据进行滤波,可以通过饱和型加法运算的办法加以消除。 饱和加法器的特点是,当相加器的输入之和大于1或小于-1时,就分别以1和-1代表相加结果,其特性如图所示,这样就能克服溢出振荡 补码饱和加法器特性

7.4.5 浮点运算中的有限字长效应 乘法和加法之后都需要进行尾数量化。 7.4.5 浮点运算中的有限字长效应 乘法和加法之后都需要进行尾数量化。 流图中乘法支路之后或加法节点处都将叠加—个量化噪声 舍入或截尾处理只影响尾数的字长,不影响阶码,但所产生的误差值却与阶码值有关。 分析中用相对误差比用绝对误差更适合 系统的输出误差大小与系统的结构有关 本节上一部分系统结构对输出误差影响的结论同样成立 取相同的尾数字长时,浮点运算的误差要比定点运算的误差小。 浮点制中可以充分利用尾数的字长,使每一个数不论其数值大小,都能保持同样的相对精度 浮点制系统的字长一定时,其输出的信噪比为一常数; 定点制系统中,输出噪声的方差与信号无关,因此信号越大,输出的信噪比越大

7.5 FFT算法的有限字长效应 从运算的角度分析,FFT与数字滤波器一样都是系统,因此,二者的有限字长效应分析方法基本上是相同的 以时间抽取DIT为例进行有限字长效应分析,并且针对的是舍入情况,其他FFT算法及截尾运算结果是相似的

7.5.1 蝶形运算的统计模型 定点制运算中只有乘法引入量化误差,加法运算不引入误差 蝶形运算中只在乘系数 时引入一个误差源

是复数相乘,所产生误差 就是一个复数 一个复数相乘要由四个实数相乘来构成、每个实数相乘都将引入一个相应的误差。因而共有四个误差 假设是舍入量化,则

复数 的方差 为 分析乘以 后的方差影响 通过所有后级蝶形时,其方差保持不变

每一个输出端都与N-1个蝶形相连, 即有N-1个量化噪声源对每个输出端有贡献 在终端,在离散傅里叶变换 上叠加的输出噪声的均方值为 和FIR滤波器的直接型实现一样, 输出噪声的总方差正比于N

7.5.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级计算完为止

如果在整个变换过程中,一次溢出也没有发生,输出信号方差最大,当输入为白噪声型信号时,其值为 输出信噪比为 最坏的情况是各级都产生了衰减,这相当于第二种方法。所以成组浮点法快速傅立叶变换的信噪比的变化范围为

7.5.3 浮点FFT算法中的量化效应 采用的分析方法基本相同 与定点FFT算法误差分析不同之处有 对输入信号和舍入噪声做同样的假设 对任一输出节点,考虑与其相连的那些蝶形 与定点FFT算法误差分析不同之处有 浮点加法和浮点乘法之后都需做尾数处理,都会引入噪声源,因此模型设计与定点FFT系统有所不同 模型中的绝对误差换需算成相对误差,再进行分析 浮点制FFT系统任一输出节点处的输出信噪比为