第6章数字信号处理系统的实现.ppt

【例6.1.1】假设某IIR数字滤波器的系统函数为 试画出该数字滤波器的直接II型结构。 解：（1）将系统函数化成式（6.1.1）形式 假设信号量化后用Q[ x(n) ]表示，定义量化误差： 在一般情况下，x(n) 为随机信号，那么e(n) 也是随机信号，因此也经常被称为量化噪声。 要精确知道量化噪声的大小是很困难的，也是没有必要的。因此，通过分析量化噪声的统计特性来描述量化误差，并对其统计特性作如下假定： 首先设计两个MATLAB函数用来处理系数量化的截尾和舍入处理。 函数a2dT(d,n)实现将十进制数据d（在MATLAB中用浮点表示假设为理想精度）按照截尾处理算法处理成n位二进制数，函数的返回值为该n位二进制数的十进制值。函数的代码为： function beq = a2dT(d,n) m = 1; d1 = abs(d); while fix(d1) 0 d1 = abs(d)/(2^m); m = m+1; end beq = fix(d1*2^n); beq = sign(d).*beq.*2^(m-n-1); 下面以6阶椭圆低通滤波器为例进行研究，假设低通滤波器的截止频率为0.4，通带波动0.5dB，阻带最小衰减为50 dB。当采用5位截尾处理时，分析其量化前后的滤波器特性，比较零极点的变化。MATLAB的程序代码如下： [b,a] = ellip(6,0.5,50,0.4); [h,w] = freqz(b,a,512);g = 20*log10(abs(h)); bq = a2dT(b,5); aq = a2dT(a,5); [hq,w] = freqz(bq,aq,512);gq = 20*log10(abs(hq)); plot(w/pi,g,b,w/pi,gq,r:);grid; axis([0 1 -80 5]);xlabel(\omega/\pi); ylabel(Gain, dB); title(实线——量化前,虚线——量化后); zplane(b,a);hold on;plotzp(bq,aq) 当采用10位截尾处理时，分析结果如图所示，量化前后特性相差不大。 对于舍入处理，只要将处理函数改为a2dR就可. 由于 是白色的，各变量之间互不相关，即 代入上式，得 由Parseval定理 H(z)全部极点在单位圆内， 表示沿单位圆逆时针方向的圆周积分。由留数定理： 如为截尾噪声，则输出噪声中还有一直流分量 例6.2.2 一个8位A/D变换器（ ），其输出 作为IIR滤波器的输入，求滤波器输出端的量化噪声功率，已知IIR滤波器的系统函数为： 解：由于A/D的量化效应，滤波器输入端的噪声功率为： 滤波器的输出噪声功率为： 其积分值等于单位圆内所有极点留数的和。单位圆内有一个极点 z=0.999，所以 6.2.4 有限字长运算对数字滤波器的影响 DF的实现，涉及到两种运算：相乘、求和。 定点制运算中，每一次乘法运算之后都要作一次舍入（截尾）处理，因此引入了非线性。 采用统计分析的方法，将舍入误差作为独立噪声e(n)迭加在信号上，因而仍可用线性流图表示定点相乘。 定点相乘运算统计分析的流图表示 对舍入噪声e(n)作如下的假设： 1. e(n) 为平稳随机噪声序列； 2. e(n) 与输入序列 x(n) 不相关，各噪声之间也互不相关。 3. e(n) 为白色噪声； 4.在量化间隔上均匀分布（即每个噪声都是均匀等概率分布）。 有了这些条件，整个系统就可作为线性系统处理。 每一个噪声可用线性离散系统的理论求出其输出噪声，所有输出噪声经线性迭加得到总的噪声输出。 1 IIR 数字滤波器的有限字长乘法运算量化效应 下面以一个例子来讨论IIR 数字滤波器的有限字长乘法运算量化效应。 【例6.2.3】一个二阶IIR低通数字滤波器，系统函数为 ， 采用定点制算法，尾数作舍入处理，分别计算其直接型、级联型、并联型三种结构的舍入误差产生的输出噪声功率。 解：①直接型结构 根据系统函数画出直接型结构流图，并在每一个乘法支路加入一个舍入噪声，得到如图所示流图。 ②级联型 将H(z)分解 结构流图为 将 ????????????按照不同的方式组合，还可以组成其他的级联形式，请读者自己完成其