数字信号处理第2版张小虹8数字信号处理3课件教学.pptVIP

为了预防溢出，系统往往需要加合适的幅度压缩因子 z?1 z?1 z?1 z?1 z?1 z?1 y1(n) x(n) y(n) yk(n) y2(n) h1(n) hm(n) h2(n) hk(n) hi(n) ——二阶节的冲激响应 ?1m ?2m ?12 ?22 ?11 ?21 A。 只有保证系统中每一个输出节点都不溢出，即每一个 输出节点变量的绝对值都小于1，才能保证系统没有 溢出效应。 令第k个节点的输出为yk(n)，若要yk(n)不溢出，则应满 hk(n)表示从系统输入到第k个节点的单位脉冲响应，则 (8.4-16) ?yk (n)?1 k=1,2,… (8.4-15) 足 m=?? ? yk(n)=? x (n?m)hk(m) = hk(n)? x (n) 上式表明yk(n)的动态范围取决于输入信号x(n)类型及 hk(n)，对不同类型输入信号，yk(n)的动态范围不同。 x(n)是正弦信号，?yk (n)?≤max?Hk(e j?)? ； x(n)是有界信号?x(n)?≤1 ， x(n)是能量有限信号 m?n ? ? x2(n) ?1 ， ? |yk(n)|? ? |hk(n)| ； n=0 n=0 ? |yk(n)|? ? hk(n) 2 1/2 由第 k 个节点的单位脉冲响应可推出一般情况下，任一 输出节点都不溢出的条件为 (8.4-17) 由上式得到输入x(n)的最大幅度限制为 (8.4-18) 式中xmax表示输入序列最大绝对值。 m=?? ? |yk(n)|=? |x (n?m)hk(m)|? xmax ? |hk(m)|1 m=?? ? ? |hk(m)| m=?? ? xmax 1 实际应用中，通常xmax并不能满足这一要求，系统需要 乘以合适的幅度压缩因子A(A1)衰减输入信号幅度，以 保证任意第k个输出节点都不溢出，即(8.4-17)式变为 同样(8.4-18)式变为 A的选取很重要，因为A大了不起作用，而A小了会使信 号的输出功率减小，导致信噪比降低。 (8.4-19) n=?? ? |yk(n)| ? Axmax ? |hk(n)|1 (8.4-20) ? |hk(n)| n=?? ? A 1 xmax 8.4. 5、定点FIR系统舍入运算误差分析 1、输出噪声方差 N阶的FIR系统差分方程为 系统的直接(卷积)形式运算模型如图所示 y(n) h (N-1) h (N-2) h (2) h (1) h (0) A x(n) z?1 z?1 z?1 z?1 z?1 m=0 N?1 y(n)=? h(m)x(n?m) e0(n) eN-1(n) eN-2(n) e2(n) e1(n) 每做一次舍入乘法产生一次误差。在同一时刻n，要产生 N个舍入误差，将其等效为N个噪声源e0(n)，e1(n)，…, eN-1(n)，当n变化时形成随机序列。 系统的直接(卷积)形式舍入运算误差统计模型 y(n)+ f(n) h (N-1) h (N-2) h (2) h (1) h (0) A x(n) z?1 z?1 z?1 z?1 z?1 同样对ei(n) (i=0,1,2…,N?1)作如下假设 (4)误差ei(n)与输入x(n)不相关，且与系统中任何节点 (2)每个误差在量化误差范围内均匀等概分布, 方差为 (1)所有误差ei(n)是平稳的零均值白噪声序列； (3)两个不同噪声源彼此不相关； 变量不相关，从而与输出y(n)也不相关。 ， i=0,1,2…,N?1 ； 2?2b 12 2 ?ei = 此时输出为 结果表明，FIR系统定点舍入运算误差直接到达输出 端，与系统的系数无关；误差与字长b及系统的阶次 N有关， b越大或N 越小，误差越小。 则输出方差 ?en ?f = ? n=0 2 2 N?1 ei(n)= y(n)+f(n) y(n) = y(n)+ ^ ? N?1 n=0 2?2b 12 ? = N 2、输入x(n)的最大幅度限制 因为定点运算要保证各输出节点不溢出，所以与IIR系 统分析相同，FIR系统输入x(n)也有最大幅度限制要求。 由图8 .4-13可见，各支路直接至输出端，所以不溢出 的条件为 |y(n)|=? | h(m)x (n?m) |? xmax ? |h(m)|1 N?1 m=0 N?1 m=0 由此得到输入x(n)的最大幅度限制 若xmax(n)不满足条件，为保证不溢出，在输入端乘以 A(A1)，A的取值为 ? |h(n)| xmax 1 N?1 n=0 ? |h(n)| A 1 N?1 n=0 xmax 8.4.6 IIR