数字信号处理第六章 数字信号处理系统的实现

第六章数字信号处理系统的实现

数字信号处理系统的实现方法：a.利用专用计算机或FPGA等专用硬件；b.直接利用计算机和通用软件编程实现。一个数字滤波器的系统函数一般可表示为有理函数形式：

为IIR滤波器形式，{}都为0时就是一个FIR滤波器。对于这样一个系统，也可用差分方程来表示：.IIR、FIR的系统函数网络结构形式软、硬件实现DF.

即一个输出序列是其过去点的线性组合加上当前输入序列与过去点输入序列的线性组合。除了与当前的输入有关，同时还与过去的输入和过去的输出有关，系统是带有记忆的。

对于上面的算式，可以化成不同的计算形式，如直接计算、分解为多个有理函数相加、分解为多个有理函数相乘等等，不同的计算形式也就表现出不同的计算结构，而不同的计算结构可能会带来不同的效果，或者是实现简单，编程方便，或者是计算精度较高等等。另外，数字信号是通过采样和转换得到的，而转换的位数是有限的（一般6、8、10、12、16位），所以存在量化误差，另外，计算机中的数的表示也总是有限的，经此表示的滤波器的系数同样存在量化误差，在计算过程中因有限字长也会造成误差。

.量化误差主要有三种误差：①A/D变换量化效应；②系数的量化效应；③数字运算的有限字长效应。.6.1数字滤波器的结构

6.1.1数字网络的信号流图

差分方程中数字滤波器的基本操作：①加法，②乘法,③延迟。为了表示简单，通常用信号流图来表示其运算结构。对于加法、乘法及延迟这三种基本运算。..只有输出支路的节点称为输入节点或源点；只有输入支路的节点称为输出节点或阱点；既有输入支路又有输出支路的节点叫做混合节点。通路是指从源点到阱点之间沿着箭头方向的连续的一串支路。通路增益是该通路上各支路增益的乘积。如7-1-3-4-8；7-1-2-6-3-4-8；回路是指从一个节点出发沿着支路箭头方向到达同一个节点的闭合通路，它象征着系统中的反馈回路。组成回路的所有支路增益的乘积通常叫做回路增益。如：3-4-5-6-3.梅逊(Mason)公式式中Tk为从输入节点（源点）到输出节点（阱点）的第k条前向通路增益；Δ为流图的特征式Δk是不接触第k条前向通路的特征式余因子为所有不同回路增益之和.为每两个互不接触回路增益之和.例:利用梅逊公式计算图中的系统函数有两条前向通路：

一个回路，其回路增益为则系统函数

.信号流图的转置定理：对于单个输入、单个输出的系统，通过反转网络中的全部支路的方向，并且将其输入和输出互换，得出的流图具有与原始流图相同的系统函数。.信号流图转置的作用：①转变运算结构；②验证计算流图的系统函数的正确与否。运算结构对滤波器的实现很重要，尤其对于一些定点运算的处理机，结构的不同将会影响系统的精度、误差、稳定性、经济性以及运算速度等许多重要的性能。对于无限长单位冲激响应（IIR）数字滤波器与FIR数字滤波器，它们在结构上各有自己不同的特点，因此我们在下面将对它们分别加以讨论。.和每个节点连接的有输入支路和输出支路，节点变量等于所有输入支路的输出之和。在图5.2.2中，(5.2.1).图5.2.2信号流图(a)基本信号流图；(b)非基本信号流图.例求图5.2.2(a)信号流图决定的系统函数H(z)。解将5.2.1式进行z变换，得到经过联立求解得到：.

6.1.2IIR数字滤波器的结构IIR数字滤波器的结构特点：存在反馈环路，递归型结构。同一系统函数，有各种不同的结构形式。其主要结构有：

(1)直接型直接由IIRDF的差分方程所得的网络结构。..图二IIR数字滤波器的网络结构.上述结构缺点：①需要2N个延迟器（z-1），太多。②系数ai、bi对滤波器性能的控制不直接，对极、零点的控制难，一个ai、bi的改变会影响系统的零点或极点分布。③对字长变化敏感（对ai、bi的准确度要求严格）。④易不稳定，阶数高时，上述影响更大。.直接Ⅱ型上面直接型结构中的两部分可分别看作是两个独立的网络(H1(z)和H2(z)),两部分串接构成总的系统函数：由系统函数的不变性（系统是线性的），得

.直接I型的变型.两条延时链中对应的延时单元内容完全相同,可合并，得：直接II型结构.

.

.例IIR数字滤波器的系统函数H(z)为画出该滤波器的直接型结构。解由H(z)写出差分方程如下：.例图.

直接II型(正准型)优缺点：优点：延迟线减少一半，为N个，可节省寄存器或存储单元。缺点：同直接型。通常在实际中很少采用上述两种结构实现高阶系统，而是把高阶变成一系列不同组合的低阶系统（一、二阶）来实现。.（2）级联型（串联）

一个N阶系统函数可用它的零、极点表示，即把它的分子、分母都表达为因子形式

由于系数、都是实数，极、零点为实根或共轭复根，所以有

.

、——实根、——复根且将共轭因子合并为实系数二阶因子，单实根因子看作二阶因子的一个特例，则

、——为实系数。用若干二阶网络级联构成滤波器，二阶子网络称为二阶节，可用正准型结构实现。

.级联型结构.级联型结构的优缺点：

优点：①简化实现，用一个二阶节，通过变换系数就可实现整个系统；②极、零点可单独控制、调整，调整、可单独调整第对零点，调整、可单独调整第对极点；③各二阶节零、极点的搭配可互换位置，优化组合以减小运算误差；④可流水线操作。缺点：二阶节电平难控制，电平大易导致溢出，电平小则使信噪比减小。.例设系统函数H(z)如下式：试画出其级联型网络结构。解将H(z)分子分母进行因式分解，得到

.例图.（3）并联型

将系统函数展开成部分分式之和，可用并联方式构成滤波器：将上式中的共轭复根成对地合并为二阶实系数的部分分式，上式表明，可用L个一阶网络、M个二阶网络以及一个常数A0并联组成滤波器H（z），结构如下图：.并联型结构.例画出上例题中的H(z)的并联型结构。解将上例中H(z)展成部分分式形式：将每一部分用直接型结构实现，其并联型网络结构如图所示。.例图.特点：①系统实现简单，只需一个二阶节，系统通过改变输入系数即可完成；②极点位置可单独调整；③运算速度快（可并行进行）；④各二阶网络的误差互不影响，总的误差小，对字长要求低。缺点：不能直接调整零点，因多个二阶节的零点并不是整个系统函数的零点，当需要准确的传输零点时，级联型最合适。.[例]已知某三阶数字滤波器的系统函数为试画出其直接型、级联型和并联型结构。.直接型将系统函数H(z)表达为.级联型将系统函数H(z)表达为一阶、二阶实系数分式之积.并联型将系统函数H(z)表达为部分分式之和的形式.6.1.3FIR滤波器的结构FIRDF特点：主要是非递归结构，无反馈，但在频率采样结构等某些结构中也包含有反馈的递归部分。它的系统函数和差分方程一般有如下形式：.基本的结构形式有下几种：（1）直接型（卷积型、横截型）卷积型：差分方程是信号的卷积形式；横截型：差分方程是一条输入x(n)延时链的横向结构。.基本的结构形式有下几种：（1）直接型（卷积型、横截型）直接由差分方程可画出对应的网络结构：.直接型的转置：

.（2）级联型（串联型）当需要控制滤波器的传输零点时，可将系统函数分解为二阶实系数因子的形式：

于是可用二阶节级联构成，每一个二阶节控制一对零点。缺点：①所需要的系数a比直接型的h(n)多；②乘法运算多于直接型。.级联型结构.（3）线性相位型FIR的重要特点是可设计成具有严格线性相位的滤波器，此时满足偶对称或奇对称条件。偶对称时，N为偶数，

N为奇数，.

图N为偶数的线性相位FIR滤波器结构

.

图N为奇数的线性相位FIR滤波器结构

.

...46.由上两式，可得到线性相位FIR滤波器的结构,如图。优点：线相相位型结构的乘法次数减为（N偶数）（N奇数）

(横截型结构乘法次数:N次).（4）频率采样型第二章讨论了有限长序列可以进行频域采样。现是长为的序列，因此也可对系统函数H(z)在单位圆上作等分采样，这个采样值也就是的离散傅里叶变换值H(k)。根据上一章的讨论，用频率采样表达z函数的内插公式为：.H(z)由两部分级联而成，第一部分（FIR部分）

这是一个由N节延时器组成的梳状滤波器，它在单位圆上有N个等分的零点：其频响为.梳状滤波器频响.

第二部分（IIR部分）是一组并联的一阶网络：

此一阶网络在单位圆上有一个极点：该网络在处的频响为，是一个谐振频率为的谐振器。这些并联谐振器的极点正好各自抵消一个梳状滤波器的零点，从而使这个频率点的响应等于。两部分级联后，就得到频率采样型的总结构。

.频率采样型结构.

这一结构的最大特点是它的系数H(k)直接就是滤波器在处的响应，因此，控制滤波器的响应很直接。两个主要的缺点：①所有的系数和都是复数，计算复杂。②所有谐振器的极点都在单位圆上，考虑到系数量化的影响，有些极点实际上不能与梳状滤波器的零点相抵消，使系统的稳定性变差。

.

为了克服这两个缺点，作两点修正：1）将所有零点和极点移到半径为的圆上，略小于1，同时频率采样点也移到该圆上，以解决系统的稳定性。这时.2）共轭根合并，将一对复数一阶子网络合并成一个实系数的二阶子网络。这些共轭根在圆周上是对称点即同样，h(n)因是实数，其DFT也是圆周共轭对称的，.因此可将第k及第N-k个谐振器合并为一个二阶网络其中.

这个二端网络是一个有限Q值的谐振器，谐振频率为。除了以上共轭极点外，还有实数极点，分两种情况：当N为偶数时，有二个实数极点，对应H(0)和H(N/2)，有二个一阶网络：所以有.当为奇数时，只有一个实数极点，对应H(0)，有一个一阶网络：

所以有

..改进后的频率采样型结构如下图.频率采样型特点：1.选频性好，适于窄带滤波，