信号处理时域离散系统的实现1

第八章时域离散系统的实现8.4格型网络结构8.3IIR网络结构8.2FIR网络结构8.6数字信号处理中的量化效应8.5用软件实现各种网络结构本章内容:8.1引言8.1引言时域离散系统的实现方法:(a)软件实现：按所设计的软件在通用的计算机运行数字信号处理程序。优点：经济，一机可以多用.缺点：处理速度慢.(b)硬件实现：用加法器、乘法器和延时器等组成的专用数字网络设备,以实现信号的处理运算.优点：处理速度快,容易做到实时处理.缺点：不灵活,开发周期较长,且设备只能专用.在实际应用中,通常采用软硬件结合实现.返回数字滤波器的表示方法(a)常系数线性差分方程:（b）数字滤波器的系统函数：01()()()1MkkkNkkkbzYzHzXzaz10()()()NMkkkkynaynkbxnk返回回到本节数字信号处理器中的基本运算单元加法器方框图aa1z数乘器单位延时基本运算单元流图1z返回回到本节本章重点讨论下述内容IIR滤波器的基本结构FIR滤波器的直接型、级联型、线性相位结构，理解频率抽样型结构数字滤波器的格型结构返回回到本节8.2FIR网络结构它的差分方程和系统函数分别为一般称上面两式表示长度为N,阶数为N-1的FIR滤波器.10NkkHzhkz10Nkynhkxnk返回(a)没有反馈支路，即没有环路，非递归型结构。FIR网络结构特点:(b)N-1阶滤波器，N为滤波器的长度，有N-1个零点分布于z平面，z=0处是N-1阶极点。(c)其单位脉冲响应是有限长序列。返回本节主要讲述：8.2.1FIR直接型结构和级联型结构8.2.2线性相位结构8.2.3FIR频率采样结构8.2.4快速卷积法返回FIR滤波器网络结构的五种实现方法（1）直接型结构（2）级联型结构（3）线性相位型结构（4）频率取样型结构（5）快速卷积法返回8.2.1FIR直接型结构和级联型结构1.FIR直接型结构（卷积型、横截型）按照H(z)或者差分方程直接画出结构图。如图8.2.1所示y(n)h(0)h(1)h(2)h(n2)h(n1)z1z1z1x(n)图8.2.1FIR直接型结构流图特点：单位延时器串联，有抽头，称为延时线；简单直观，乘法运算量少，但不易调整零点.返回回到本节2.FIR级联型结构当需要控制滤波器的传输零点时，可将H(z)进行式分解，并将共轭成对的零点放在一起，形成一个系数为实数的二阶形式：11201201()()()MNkiiikiHzhkzaazaz这样级联型网络结构就是由一阶或二阶因子构成的级联结构，其中每一个因式都用直接型实现。返回回到本节例8.2.1设FIR网络系统函数H(z)如下式：画出H(z)的直接型结构和级联型结构。解:将H(z)进行因式分解，得到：它的直接型结构和级联型结构分别如下图所示:-1-2-3H(z)=0.96+2z+2.8z+1.5z-1-1-2H(z)=(0.6+0.5z)(1.6+2z+3z)返回回到本节y(n)x(n)z-1z1z10.9622.81.5（a）直接型结构z1z1z1x(n)0.60.51.623y(n)（b）级联型结构图8.2.3例8.2.1图返回回到本节级联型结构中,每一个一阶网络控制一个零点,调整零点只需调整该因式的两个系数;二阶网络控制一对零点,调整它也只需调整该因式的三个系数.相对于直接型结构来说:FIR级联型结构特点:1）每个基本节控制一对零点，调整零点方便。2）需要对系统函数进行因式分解，系数比直接型多，所需的乘法运算多。返回回到本节FIR滤波器单位抽样响应h(n)为实数，且满足：–第一类偶对称：–第二类奇对称：–对称中心在(N-1)/2处，这种FIR滤波器具有严格线性相位。8.2.2线性相位结构系统函数具有线性相位,它的单位脉冲响应满足下式01nN()(1)hnhNn()(1)hnhNn)()()(jgjeHeHN12N122－第一类线性相位：（）＝－－第二类线性相位：（）＝－返回回到本节10()()NnnHzhnz11202()()NNnnNnnhnzhnz12(1)0()NnNnnhnzz当N为偶数时10()()NnnHzhnz1112(1)201()2NNnNnnNhnzzhz当N为奇数时返回回到本节x(n)y(n)z－1z－1z－1z－1z－1z－1z－1h(0)h(1)h(2)h(N/2－1)x(n)y(n)z－1z－1z－1z－1z－1z－1h(0)h(1)h(2)h((N－1)/2)N＝偶数N＝奇数图8.2.4第一类线性相位网络结构流图返回回到本节x(n)y(n)z－1z－1z－1z－1z－1z－1z－1h(0)h(1)h(2)h(N/2－1)x(n)y(n)z－1z－1z－1z－1z－1z－1h(0)h(1)h(2)h((N－1)/2)N＝偶数N＝奇数－1－1－1－1－1－1－1－1－1图8.2.5第二类线性相位网络结构流图返回回到本节根据线性相位结构流图,和直接型结构比较,如果N取偶数,直接型需要N个乘法器,而线性相位结构需要N/2个,节约了一半的乘法器.如果N取奇数,则乘法器减少到(N+1)/2个,同样也节约了一半的乘法器.返回回到本节8.2.3FIR频率采样结构根据频率采样定理,在频率的区间,对系统的传输函数进行N点等间隔采样,如果N大于等于系统单位脉冲响应的长度M,不会引起信号失真,系统函数和采样值之间服从下面的内插关系02110()(1)1)1(NNkkNHHkWzzzN10(())1ckkNHzHzN2()()0,1,2,1jkNzeHkHzkN返回回到本节()1NcHzz其中子系统：是N节单位延时单元的梳状滤波器2jkNkze0,1,...,1kN在单位圆上有N个等间隔角度的零点：()1jjNcHee22sin2NjNje频率响应：222NNNjjjeee返回回到本节与第k个零点相抵消，使该频率处的频率响应等于H(k)2kN2jkkNkNzWe一阶网络单位圆上有一个极点：'1()()1kkNHkHzWz子系统：返回回到本节由此可知频率采样结构是由一个梳状滤波器N个一阶网络Hk(z)的并联结构进行级联而成.其结构如下图所示x(n)y(n)z－1z－1－z－NH(0)H(1)H(N－1)0NW1NW1NNWz－1N11210N,,,,k,WzkNk返回回到本节频率抽样型结构的优缺点：（1）调整H(k)就可以有效地调整频响特性。（2）若h(n)长度相同，则网络结构完全相同，除了各支路增益H(k)，便于标准化、模块化。（3）系数多为复数，增加了复数乘法和存储量。（4）系统稳定是靠位于单位圆上的N个零极点对消来保证的，导致系统不稳定。返回回到本节频率抽样结构的修正（1）问题：在有限字长情况下，系数量化后极点不能和零点抵消，使FIR系统不稳定。解决方案：将零极点移至半径为r的圆上11rr且1101()()(1)1NNNrkkNHkHzrzNrWz2jkNkzre极点：0,1,...,1kN返回回到本节频率抽样结构的修正（2）H(k)和都是复数H(k)的分布,等间隔采样，关于N/2共轭对称*Nkkzz()*()NkkkNN由对称性：kNW返回回到本节将第k个和第(N-k)个谐振器合并成一个实系数的二阶网络：1()1ˆ()()()()()11kkNkkNkNNHzHzHzHkHNkrWzrWz*1*1()()11()kkNNHkHkrWzrWz101122212cos()kkzzrkrzN012Re[()]2Re[()]kkkNHkrHkW11,2,...,21,2,...,12NkNNkN为奇数为偶数返回回到本节二阶网络都是实系数的,其结构图为当N为偶数时11/2011112211(0)(/2)()(1)[]21112cos()NNNkkkaazHHNHzrzNrzrzrkzrzNz1-r2)2cos(2kNrz1α0kα1k式中，H(0)和H(N/2)为实数。频率采样修正结构由N/2-1个二阶网络和两个一阶网络并联构成，如图所示。返回回到本节当N为奇数时,只有一个采样值H(O)为实数,此时x(n)y(n)z1H(0)zNr1/NH1(z)H2(z)z1rH(N/2))(12zHNNr图8.2.7频率采样结构N=偶数1(1)/201112211(0)()(1)[]2112cos()NNNkkkaazHHzrzNrzrkzrzN返回回到本节例：设计一M阶实系数FIR，已知H(0)=H(1)=1，画出其频率取样型结构。111(1)11111()111NNNNzHzNzWzWz1112222cos()11()2112cos()NzzNHzNzzzN解：频率抽样点数N=M+1由H(N1)=H(1)=1,和1)1(NNNWW返回回到本节实系数频率取样型结构流图)2cos(2N)2cos(2Nx[k]y[k]1/NNzz1z1z1111优点：1.H(m)零点较多时，实现较为简单。2.可以构成滤波器组，实现信号的频谱分析。返回回到本节8.3IIR网络结构IIR网络结构的特点:信号流图中含有反馈支路,即含有环路,递归型结构；其单位脉冲响应序列是无限长的.其网络基本结构有直接型、级联型和并联型三种.01()()()1MkkkNkkkbzYzHzXzaz系统函数：10()()()NMkkkkynaynkbxnk差分方程：返回本节主要讲述：8.3.1IIR直接型网络结构8.3.2IIR级联型结构8.3.3IIR并联型网络结构8.3.4转置型网络结构返回8.3.1IIR直接型网络结构考虑N阶差分方程,即其系统函数01()()()MNkkkkynbxnkaynk01()1MkkkNkkkbzHzaz返回回到本节令,其中取M=N=2,H(z)的实现结构如下图(a)所示12()()()HzHzHz10()MkkkHzbz211()1NkkkHzazH1(z)H2(z)x(n-1)x(n-2)bb2b0Z-11Z-1x(n)Z-1-a1y(n-1)y(n-2)Z-1-a2y(n)（a）返回回到本节zHzHzHzHzH1221将H1(z)和H2(z)交换位置,结点变量w1=w2，即输入结点变量相等，对应延时支路输出结点变量相等,其结构图如下图(b)观察上图,结点w1=w2,前后延时支路合并,可以得到下图8.3.1(c)所示的IIR直接型网络结构由于H(z)2H(z)1x(n)y(n)b0b1b2Z-1Z-1-a1-a2w2w1Z-1Z-1（b）返回回到本节x(n)y(n)-a1-a2b0b1b2z1z1图8.3.1(c)IIR直接型网络结构由上图可以看出，IIR直接型网络结构需要M+N+1次乘法，M+N次加法，延时单元数为M和N中较大的数.返回回到本节IIR直接型结构特点优点：可直接由传输函数或差分方程画出网结构流图，简单直观。缺点:（1）调整零、极点困难；（2）对参数的量化非常敏感，这是由极点对系数的变化过于敏感造成的；（3）容易产生较大误差。返回回到本节例8.3.1设IIR数字滤波器的系统函数H(z)为写出系统的差分方程,并画出