数字信号处理第三西安科大出高西全丁玉美课后答案.pptx

4.1 教材第教材第5章学习要点章学习要点数字信号处理系统设计完毕后，得到的是该系统的系统函数或者差分方程，要实现还需要按照系统函数设计一种具体的算法。不同的算法会影响系统的成本、运算的复杂程度、运算时间以及运算误差等。教材第5章的学习要点如下：(1)由系统流图写出系统的系统函数或者差分方程。第1页/共114页(2)按照FIR系统的系统函数或者差分方程画出其直接型、级联型和频率采样结构、FIR线性相位结构以及用快速卷积法实现FIR系统。(3)按照IIR系统的系统函数或者差分方程画出其直接型、级联型、并联型。(4)一般了解格型网络结构，包括全零点格型网络结构系统函数、由FIR直接型转换成全零点格型网络结构；全极点格型网络结构及其系统函数。第2页/共114页 4.2 按照系统流图求系统的按照系统流图求系统的 系统函数或者差分方程系统函数或者差分方程具体的网络结构一般用流图表示。掌握教材第5章内容就是必须能根据流图正确地求出系统函数。求系统函数的方法在先修课“信号与系统”中已讲过，这里仅帮助大家复习。求系统函数的方法有两种。一种是先根据流图写出各节点的节点方程，联立节点方程，求出输入和输出之间的关系，得到系统函数；另一种是根据梅荪（Masson）公式直接写出系统函数。显然，后一种简单。下面仅介绍根据Masson公式直接写出系统函数的方法。第3页/共114页 按照梅荪公式写出系统函数为式中，称为流图特征式，其计算公式如下：第4页/共114页式中,表示所有的环路增益之和；表示所有的每两个互不接触的环路增益乘积之和；表示所有的每三个互不接触的环路增益乘积之和；Tk表示从输入节点到输出节点的第k条前向支路的增益；k表示不与第k条前向通路接触的值。下面用例题说明利用梅荪公式直接写系统函数的方法。第5页/共114页例 写出图中流图的系统函数。图第6页/共114页解解：该流图有两个环路，一个是w2w2w2，另一个是w2w2w1w2，环路增益分别为a1z1和a2z2。没有互不接触的环路，这样流图特征式为=1(a1z1a2z2)=1+a1z1+a2z2流图中有三条前向通路：第一条T1是x(n)w2y(n)，它的增益是T1=b0；第二条T2是x(n)w2w2y(n)，它的增益是T2=b1z1；第三条T3是x(n)w2w2w1y(n)，它的增益是T3=b2z2。第7页/共114页流图中的两个环路均与所有的前向通路相接触，因此对应于三条前向通路的1=1，2=1，3=1。这样可以直接写出该流图的系统函数为第8页/共114页 4.3 按照系统函数或者差分方程画系统流图按照系统函数或者差分方程画系统流图 按照系统函数设计系统的实现方法主要依据的是系统函数的特点和要求，画出系统流图，然后根据流图设计用硬件或软件进行实现。系统的网络结构有很多，但最基本的是FIR和IIR网络结构。这两类结构各有特点。FIR结构一般没有反馈回路，单位脉冲响应是有限长的，系统稳定，但相对IIR结构，FIR结构的频率选择性不高，换句话说，要求频率选择性高时，要求FIR有很高的阶数。第9页/共114页FIR中主要有直接型结构、线性相位结构和频率采样结构。IIR网络结构主要有直接型结构、级联型结构和并联型结构。IIR网络结构有反馈回路，单位脉冲响应是无限长的，存在稳定性问题，但频率选择性高。画这些结构的流图时，最好能熟悉前一节介绍的梅荪公式，这样画起来得心应手。第10页/共114页中的线性相位结构中的线性相位结构 FIR线性相位系统具有以下特点:(1)FIR线性相位系统单位脉冲响应满足下式:h(n)=h(Nn1)式中,h(n)是实序列;N表示序列的长度。该式说明h(n)对序列的(N1)/2位置偶对称（公式中取“”号）或奇对称（公式中取“”号）。(2)FIR线性相位系统系统函数满足下面公式:N为偶数第11页/共114页N为奇数(3)FIR线性相位系统零点分布具有四个一组的特点，即如果z1是零点，那么z1*、z11、(z11)*也是零点。以上三点的分析和公式推导请参考教材第5章内容。只要满足上面任意一个特点，就可以判断该系统具有线性相位的特点。按照该系统函数的特点，就可以构成它的线性相位结构，因此并不是所有FIR系统都能形成线性相位结构。线性相位结构的优点是能节约近一半的乘法器。第12页/共114页中的频率采样结构中的频率采样结构 由频率采样定理得到公式：式中,H(k)是在02区间对传数函数等间隔采样N点的采样值，可以对单位脉冲响应h(n)进行DFT得到。这里要注意采样点数必须大于等于h(n)的长度，否则会发生时域混叠现象。因为IIR系统的单位脉冲响应是无限长的，因此不能用频率采样结构实现。第13页/共114页该公式是频率采样结构的基本公式，但它是一个不考虑稳定性，又可以应用复数乘法器的公式。为了稳定，且使用实数乘法器，应使用如下公式：当N为偶数时，第14页/共114页当N为奇数时，式中a0k=2ReH(k)第15页/共114页中的级联结构和并联结构IIR基本结构有直接型、级联型和并联型。一般低阶的用直接型，高阶的用级联型或并联型。在设计级联型结构时，需要将分子式和分母式进行因式分解，阶数高时可借助于MATLAB语言用计算机解决。设计并联结构时要进行部分分式展开。部分分式展开要求分子多项式的阶数低于分母多项式的阶数，否则是一个假分式（分子多项式的阶数不低于分母多项式的阶数），要将其化为整数和真分式之和，然后再对真分式进行部分分式展开。部分分式的各系数通过待定系数法解决。部分分式的一般表达式为第16页/共114页式中,pk是极点l，C是常整数，Ak是展开式中的系数。一般pk、Ak都是复数。为了用实数乘法，将共轭成对的极点放在一起，形成一个二阶网络，公式为第17页/共114页上式中的系数均是实数。总的系统函数为式中，L是(N+1)/2的整数部分。当N为奇数时，Hk(z)中有一个是实数极点。按照上式形成IIR的并联型结构，其中每一个分系统均是一阶网络或者是二阶网络。每个分系统均用直接型结构。第18页/共114页4.4 例例 题题 例例 设FIR滤波器的系统函数为求出其单位脉冲响应，判断是否具有线性相位，画出直接型结构和线性相位结构（如果存在）。第19页/共114页图第20页/共114页解：单位脉冲响应为序列的长度为N=5，序列对n=2对称，因此系统具有线性相位特性。画出其直接型结构和线性相位结构如图和(b)所示。第21页/共114页例例 假设系统函数如下式，画出它的并联型结构。解：上式的分子分母是因式分解形式，再写成下式:上式的第二项已是真分式，可以进行因式分解。第22页/共114页第23页/共114页再根据等式两边同次项系数必须相等的法则确定系数B和C，得到B16，C20最后得到按照上式画出系统并联结构的流图如图所示。第24页/共114页图第25页/共114页例例 为了保证滤波器的因果稳定性，其系统函数的极点必须保证全部集中在单位圆内。如果有极点在单位圆上，则可以形成一个正弦波发生器。利用这一原理试设计正弦波发生器。解解：假设有两个系统函数第26页/共114页令x(n)=A(n)，X(z)=A,得到见配套教材第54页表，上面两式对应的时域信号分别为y1(n)=A sin(0n)u(n)y2(n)=A cos(0n)u(n)第27页/共114页上面两式说明系统H1(z)和H2(z)分别在x(n)=A(n)的激励下可以分别产生正弦波和余弦波。H1(z)和H2(z)的极点为p1,2=ej0,这正是在单位圆上的两个极点，极点的相角为0。这样，H1(z)和H2(z)可以分别称为正弦波和余弦波发生器，画出H1(z)实现结构图如图所示，共需要两个乘法器、两个加法器和两个移位器。运行时要用x(n)=A(n)作激励。也可以令图中的v(n)起始条件为v(0)=A，v(1)=0，v(2)=0，代替输入信号x(n)=A(n)。在实际应用中有时需要两个正交相位正弦波，可以将H1(z)和H2(z)进行组合，同时产生正弦波和余弦波，实现结构如图所示。第28页/共114页图第29页/共114页图第30页/共114页例 研究一个FIR滤波器，其频率响应函数为式中,n0不一定为整数。设该系统的单位脉冲响应h(n)的长度N=15，n0=15/2，且第31页/共114页(1)画出该系统的频率采样结构;(2)求出系统的单位脉冲响应h(n)，并画出直接型结构，要求用最少的乘法器。解解：(1)已知代入N=15，n0=15/2及Hk的值,得到第32页/共114页由频域采样的z域内插公式有第33页/共114页系统频率采样结构如图所示。图第34页/共114页(2)显然，式中,h(0)=0。因为h(n)=h(Nn)，所以其直接型结构的高效形式（乘法运算最少）如图所示。第35页/共114页图第36页/共114页 4.5 教材第教材第9章学习要点章学习要点数字信号处理的实现方法一般有软件实现和硬件实现两种。教材第9章主要学习一般实现中的有关重要问题。教材第9章学习要点如下：(1)数字信号处理的实现中的重要问题是运算误差问题，运算误差主要来自于有限字长效应，表现在数字量化及其量化误差上。量化误差引起量化效应，量化效应主要有A/D变换器中的量化效应、系数量化效应、运算量化误差等。这些量化效应主要和计算中用的寄存器长度有关，寄存器长度愈长，量化效应愈小。第37页/共114页(2)A/D变换器中的量化效应使A/D变换器输出端的信噪比降低，如果不考虑输入信号中的噪声，仅考虑A/D变换器中的量化效应，A/D变换器输出端的信噪比为(3)系数量化效应会影响系统的频率特性，表现在使系统的零、极点位置改变。极点变化严重时，会使系统不稳定。为减少极点位置对于量化效应的敏感程度，应尽量加长寄存器长度，尽量采用阶数低的结构，以及极点不很密集的结构。第38页/共114页(4)运算量化效应主要表现在定点运算中的乘法运算中以及浮点运算中的加法、乘法运算中。运算量化效应会使网络输出端的信噪比降低。运算量化效应的大小主要和寄存器的长度有关，它的长度愈长，运算量化效应愈小。另外，也和网络结构有关，比较起来，一般直接型结构的运算量化效应较大，级联型结构的次之，并联型的最小。(5)注意在加法运算中可能会产生溢出问题，要考虑适当加防溢出的措施。(6)数字信号处理有软、硬件两种实现方法，配套教材中主要介绍了软件实现方法，其中包括如何考虑网络结构的软件实现方法。第39页/共114页 4.6 教材第教材第5章习题与上机题解答章习题与上机题解答1.已知系统用下面差分方程描述:试分别画出系统的直接型、级联型和并联型结构。式中x(n)和y(n)分别表示系统的输入和输出信号。第40页/共114页解：将原式移项得将上式进行Z变换，得到第41页/共114页(1)按照系统函数H(z)，根据Masson公式，画出直接型结构如题1解图（一）所示。题1解图（一）第42页/共114页(2)将H(z)的分母进行因式分解：按照上式可以有两种级联型结构:第43页/共114页画出级联型结构如题1解图（二）（a)所示。画出级联型结构如题1解图（二）（b)所示。第44页/共114页题1解图（二）第45页/共114页(3)将H(z)进行部分分式展开：第46页/共114页第47页/共114页根据上式画出并联型结构如题1解图（三）所示。题1解图（三）第48页/共114页2 设数字滤波器的差分方程为试画出系统的直接型结构。解解：由差分方程得到滤波器的系统函数为画出其直接型结构如题2解图所示。第49页/共114页题2解图第50页/共114页3.设系统的差分方程为y(n)=(a+b)y(n1)aby(n2)+x(n2)+(a+b)x(n1)+ab式中,|a|1，|b|1/2，对上式进行逆Z变换，得到第84页/共114页16.画出题15图中系统的转置结构，并验证两者具有相同的系统函数。解：按照题15图，将支路方向翻转，维持支路增益不变，并交换输入输出的位置，则形成对应的转置结构，画出题15图系统的转置结构如题16解图所示。将题16解图和题15图对照，它们的直通通路和反馈回路情况完全一样，写出它们的系统函数完全一样，这里用Masson公式最能说明问题。第85页/共114页题16解图第86页/共114页题17图17.用b1和b2确定a1、a2、c1和c0，使题17图中的两个系统等效。第87页/共114页解：题17图（a）的系统函数为题16图（b）的系统函数为对比式和式，当两个系统等效时，系数关系为a1=b1,a2=b2c0=2,c1=(b1+b2)第88页/共114页18.对于题18图中的系统，要求：（1）确定它的系统函数；（2）如果系统参数为 b0=b2=1,b1=2,a1=1.5,a2=0.9 b0=b2=1,b1=2,a1=1,a2=2画出系统的零极点分布图，并检验系统的稳定性。解解：（1）第89页/共114页（2）b0=b2=1,b1=2,a1=1.5,a2=0.9零点为z=1(二阶)，极点为 p1,2=0.750.58j,|p1,2|=0.773极零点分布如题18 解图(a)所示。由于极点的模小于1，可知系统稳定。第90页/共114页题18图第91页/共114页题18解图第92页/共114页 b0=b2=1,b1=2,a1=1,a2=2零点为z=1(二阶)，极点为 p1,2=0.51.323j,|p1,2|=1.414极零点分布如题18解图(b)所示。这里极点的模大于1，或者说极点在单位圆外，如果系统因果可实现，收敛域为|z|1.414，收敛域并不包含单位圆，因此系统不稳定。第93页/共114页19*.假设滤波器的系统函数为在单位圆上采样六点，选择r0.95，试画出它的频率采样结构，并在计算机上用DFT求出频率采样结构中的有关系数。解解：第94页/共114页式中，分母分子多项式各有一个零点z=1，相互抵消，因此该系统仍然稳定，属于FIR系统。由系统函数得到单位脉冲响应为h(n)=5(n)+5(n1)+5(n2)+3(n3)+3(n4)+3(n5)H(k)=DFTh(n)k=0,1,2,5第95页/共114页按照上式画出频率采样修正结构如题19*解图所示。图中系数a0k=2ReH(k),a1k=2RerH(k)W6k求系数程序ex519.m如下：%程序ex519.mhn=5，5，5，3，3，3；r=0.95；Hk=fft(hn，6)；for k=1：3，hk(k)=Hk(k)；Wk(k)=exp(j*2*pi*(k1)/6)；endH0=Hk(1)；H3=Hk(4)；r0k=2*real(hk)；r1k=2*real(r*hk.*Wk)第96页/共114页题19*解图第97页/共114页程序运行结果：H(0)=24H(3)=2r0k=48 4 0r1k=45.6000 3.8000 0得到01=48,02=4,11=45.2,12=38 进一步的说明：此题h(n)的长度为6，由单位圆上采样6点得到频率采样结构，满足频率采样定理。但如果采样点数少于6点，则不满足频率采样定理，产生时域混叠现象。第98页/共114页20.已知FIR滤波器的系统函数为：（1）H(z)=1+0.8z1+0.65z2 （2）H(z)=10.6z1+0.825z20.9z3试分别画出它们的直接型结构和格型结构，并求出格型结构的有关参数。解：已知FIR滤波器的系统函数，设计相应的格型结构需要用到的公式如下：ak=h(k)l=1,2,N第99页/共114页式中,N是FIR滤波器的阶数，h(k)是其单位脉冲响应，kl是格型结构的系数。（1）画出直接型结构如题20解图（a）所示。h(n)=(n)+0.8(n1)+0.65(n2)k1=0.485画出格型结构如题20解图（b）所示。第100页/共114页(2)画出直接型结构如题20解图（c）所示。H(z)=10.6z1+0.825z20.9z3h(n)=(n)0.6(n1)+0.825(n2)0.9(n3)k3=0.9第101页/共114页 k1=0.3画出直接型结构如题20解图（d）所示。第102页/共114页题20解图第103页/共114页21.假设FIR格型网络结构的参数k1=0.08,k2=0.217,k3=1.0,k4=0.5，求系统的系统函数并画出FIR直接型结构。解解：用到的公式重写如下：1kl1;l=1,2,N（该题N=3）第104页/共114页最后得到画出它的直接型结构如题21解图所示。系统函数为第105页/共114页题21解图第106页/共114页22.假设系统的系统函数为 H(z)=1+2.88z1+3.4048z2+1.74z3+0.4z4 要求:(1)画出系统的直接型结构以及描述系统的差分方程;(2)画出相应的格型结构，并求出它的系数;(3)判断系统是否是最小相位。解解：(1)系统的差分方程为 y(n)=x(n)+2.88x(n1)+3.4048x(n2)+1.74x(n3)+0.4x(n4)它的直接型结构如题22解图（一）所示。第107页/共114页题22解图（一）第108页/共114页(2)N=4，第109页/共114页由以上得到k1=0.863,k2=1.123,k3=0.684,k4=0.4第110页/共114页题22解图（二）画出其格型结构如题22解图（二）所示。第111页/共114页(3)由系统函数求出系统的零点为1.0429+0.6279i 1.04290.6279i0.3971+0.3350i 0.39710.3350i画出系统的零极点图如题22解图（三）所示。因为系统有两个零点在单位圆外，因此系统不是最小相位系统。第112页/共114页题22解图（三）第113页/共114页感谢您的观看！第114页/共114页