ch09数字信号处理的实现.ppt

　　最后要指出的是按照（9.1.17）式或（9.1.19）式选择衰减因子是比较保守或者说是比较苛刻的。经常用下式计算：  　　  式中，δ是大于1的数，如果输入信号是方差为1的白噪声，可选δ≥5。 （9.1.20） 9.2 数字信号处理技术的软件实现 　　数字信号的处理可以用软件实现，也可以用硬件实现。软件实现指的是在通用计算机上执行数字信号处理程序。这种方法灵活，但一般不能实现实时处理。硬件实现是利用数字信号处理专用集成电路或单片数字信号处理器DSP(Digital Signal Processor)来实现的。目前，这些器件一般都照顾到数字信号处理的特点，内部带有乘法器、累加器，采用流水线工作方法以及并行结构，多总线，速度快，并配有适合数字信号处理的指令等。 　　一些特殊的专用器件有： FFT专用芯片，FIR滤波器、卷积和相关等专用芯片，它的软件算法已在芯片内部用硬件实现，随着超大规模集成电路的发展，DSP芯片成本在不断下降，从而使这种实现方法成为数字信号处理的主导方法。本节主要介绍软件实现方法。 　　（3） 极点偏差与系统函数的阶数N有关，阶数愈高，极点灵敏度愈高，极点偏差也愈大。这样对于一些窄带滤波器，因为要求选择性高，势必要求阶数高，极点的偏差会很大。严重时使极点移到单位圆上或者单位圆外，引起系统不稳定。 　 考虑以上因素，系统的结构最好不用高阶的直接型结构，而将其分解成一阶或者二阶系统， 再将它们进行并联或者串联，以便减小极点偏移量。 　　【例 9.1.2】 按照例6.5.2数字带通滤波的设计指标： 系统采样频率Fs=8 kHz，要求保留2025～2225 Hz频段的频率成分，幅度失真小于1 dB；滤除0～1500 Hz和2700 Hz以上频段的频率成分，衰减大于40 dB。调用MATLAB滤波器设计分析工具FDATool，设计该滤波器，并对其系数用16位字长量化，其中尾数15位。 　　直接型结构系数量化前后的零极点分布图和频响特性曲线分别如图9.1.4(a)和（b）所示。级联型结构系数量化前后的零极点分布图和频响特性曲线分别如图9.1.4(c)和（d）所示。 图9.1.4 窄带滤波器系数量化前后的零极点分布图和 频响特性曲线比较 　　由图可见，六阶椭圆带通滤波器的直接型结构的极点对系数量化误差的敏感度高，相应的极点偏差大，量化误差使频响曲线偏差很大。量化后的滤波器无法使用。但级联型结构的零极点对系数量化误差的敏感度很低，相应的极点几乎看不出偏差，量化前后频响曲线基本重合。所以，工程上高阶滤波器一般都用级联型结构和并联型结构。 　　观察图9.1.4(a)和（c）发现，直接型结构和级联型结构系数量化前后的零点偏移都很小。这是因为该滤波器的各零点之间距离较大（不密集）。而窄带带通滤波器的极点分布密集，所以直接型结构的极点对系数量化误差的敏感度高，相应的极点偏差大。本例的实验结果与理论结论一致。 9.1.4 数字系统中的运算量化效应 　　在数字网络的运算中，其中间结果和最后结果的位数，如果超出了规定的有限位二进制数长度，则需要进行尾数处理，这样便引起了运算量化误差；运算中还可能出现溢出，造成更大的误差；运算误差的大小除了和规定的二进制数的长度有关以外，还和网络结构有关。下面就以上三个问题进行介绍。 　　1. 运算量化效应　 　　在定点制运算中，二进制乘法的结果尾数可能变长，需要对尾数进行截尾或舍入处理；在浮点制运算中无论乘法还是加法都可能使二进制的位数加长，也需要对尾数进行截尾或舍入处理。这样不管是采用定点制还是浮点制，都会因运算差生量化误差，这种误差称为运算量化误差。下面我们仅介绍定点制的乘法量化效应。 　　由于输入信号是随机信号，差生的运算量化误差同样是随机的，需要进行统计分析。运算量化误差在系统中起噪声作用，会使系统的输出信噪比降低。 　　为了分析计算简单，假定运算量化误差具有以下统计特性： ① 系统中所有的运算量化噪声都是平稳的白噪声；② 所有的运算量化噪声之间以及和信号之间均不相关；③ 这些噪声的概率密度都是均匀分布的。 　　假设定点乘法运算按b位进行量化，量化误差用e(n)表示。对于一个乘法支路，如图9.1.5（a）所示，图中节点变量v2(n)=av1(n)，经过量化后用　　　　　　 表示，那么   　　这样量化以后乘法支路的统计模型如图9.1.5（b）所示。因此系统中所有的乘法支路都和图9.1.5（b）一样引入一个噪声源。 图9.1.5 乘法支路及其量化模型 　　图9.1.6中，有两个乘法支路，采用定点制时共引入两个噪声源，即