信号与系统第8章-离散时间系统能够的时域与变换域

1、1第8章 离散时间系统的时域与 变换域分析8.6 应用应用MATLAB分析离散时间系统分析离散时间系统8.5 数字滤波器的一般概念数字滤波器的一般概念8.4 离散系统的频响特性离散系统的频响特性8.3 离散系统的单位样值响应和系统函数离散系统的单位样值响应和系统函数8.2 常系数线性差分方程的求解常系数线性差分方程的求解8.1 离散时间系统与差分方程离散时间系统与差分方程28.1 8.1 离散时间系统与差分方程离散时间系统与差分方程8.1.1 线性时不变离散时间系统线性时不变离散时间系统离散时间系统可以看成为一个离散信号的变换器，当输入信号xn经过该离散系统后，将变换成另一个序列-输出信号yn

2、，其框图如图8.1-1所示。最基本的一类系统：线性时不变离散时间系统线性时不变离散时间系统 线性离散系统线性离散系统是指满足叠加性与均匀性的离散系统。是指满足叠加性与均匀性的离散系统。 38.1.1 8.1.1 线性时不变离散时间系统线性时不变离散时间系统时不变离散系统时不变离散系统是指在同样起始状态下，系统响应与激励施加于系统的时刻无关。即：若激励信号xn产生的响应为yn，则激励信号xn - m产生的响应为yn - m，即发生同步延迟。 48.1.1 8.1.1 线性时不变离散时间系统线性时不变离散时间系统例例8.1-1 设某离散系统激励xn与响应yn之间的关系为 yn = nxn，判断该系

3、统是否为线性时不变系统。解：解： 设y1n和y2n分别为输入x1n和x2n的响应，即 y1n = nx1n，y2n = nx2n（1）当xn= ax1n时，yn = n(ax1n) = anx1n = ay1n， 该系统满足均匀性。（2）当xn = x1n+x2n时， yn = n(x1n +x2n)= nx1n +nx2n= y1n+ y2n 该系统满足叠加性。所以该系统是所以该系统是线性系统线性系统。（3）当xn= x1n-m时，yn = nx1n-m。而y1n-m=(n-m)x1n-m， 显然yn y1n-m， 该系统不是时不变系统不是时不变系统（或称（或称时变系统时变系统）。综合以上讨

4、论，该系统是一个线性时变系统线性时变系统。 58.1.1 8.1.1 线性时不变离散时间系统线性时不变离散时间系统例例8.1-2 判断滑动平均滤波器的线性特性及时不变特性。广义的滑动平均系统的输出yn与输入xn满足以下关系 21121 T 1MkMy nx nx nkMM解：解：假设 y1n =Tx1n和y2n = Tx2n，即y1n和y2n分别为输入x1n和x2n时的输出信号。（1）当输入信号为xn= ax1n时，输出信号为2111121 T 1MkMy nax nax nkMM因而该系统满足均匀性。211112 1MkMax nkay nMM68.1.1 8.1.1 线性时不变离散时间系统

5、线性时不变离散时间系统（2）输入信号为xn= x1n+x2n时，输出信号为212211121212121212121 T ( )11111 MkMMMkMkMynx nx nx n kx n kMMx n kx n kMMMMy ny n21111121 T 1MkMy nx nmx nmky nmMM该系统满足叠加性，所以该系统是线性系统线性系统。（3）假设输入信号为xn= x1n-m，则输出信号为因而该系统是时不变系统时不变系统。综合以上讨论，该系统是一个线性时不变系统线性时不变系统。78.1.2 8.1.2 差分方程差分方程离散时间系统的基本运算单元符号离散时间系统的基本运算单元符号 8

6、8.1.2 8.1.2 差分方程差分方程例例8.1-3 考察图8.1-5所示的离散系统，它由移位、相加、乘系数等三个基本单元组合而成，试写出其激励xn和响应yn之间关系的差分方程。乘系数单元：输入为yn-1，输出为ayn-1；加法器单元：输入为xn和ayn-1，输出为yn。因此，针对加法器可以写出： yn = xn + ayn-1移项整理可得： yn-ayn-1 = xn(8.1-2)-一阶常系数线性后向差分方程一阶常系数线性后向差分方程 解：解：单位移位（延时）器：输入为yn，输出为yn-1；98.1.2 8.1.2 差分方程差分方程-一阶常系数线性前向差分方程一阶常系数线性前向差分方程 例

7、例8.1-4 将图8.1-5所示的离散系统中的延时器位置稍做调整，组成如图8.1-6所示的系统，试写出其输入、输出关系式。延时器的输出为yn，则输入必为yn+1，即加法器输出为yn+1，针对加法器可写出 yn+1 = xn + ayn解：解：即 yn+1-ayn = xn (8.1-3)108.1.2 8.1.2 差分方程差分方程例例8.1-5 对图8.1-7所示的离散系统，试写出其输入、输出关系式。yn-a1yn-1-a2yn-2 = b0 xn + b1xn-1 + b2xn-2 (8.1-4)解：解：118.1.2 8.1.2 差分方程差分方程 例：例：yn-ayn-1 = xn设输入

8、xn=n，并假设 y-1=0， y0 = x0 + ay-1 = 1y1 = x1 + ay0 = ay2 = x2 + ay1 = a 2yn = xn + ayn-1 = a n此范围仅限于n 0，故 yn = anun差分方程的求解方法差分方程的求解方法 1递推解法（迭代法）2时域经典法3零输入、零状态响应解法4z变换法5状态空间分析法128.2 8.2 常系数线性差分方程的求解常系数线性差分方程的求解00NMkrkra y nkb x nr(8.2-2)8.2.1 8.2.1 线性常系数差分方程的时域经典法求解线性常系数差分方程的时域经典法求解一般地，常系数线性差分方程的解由齐次解和特

9、解齐次解和特解组成。 式(8.2-2)所对应的齐次方程的形式为：00Nkka y nk(8.2-3)10110NNNNaaaa(8.2-8) 特征方程特征方程138.2.1 8.2.1 线性常系数差分方程的时域经典法求线性常系数差分方程的时域经典法求解解其中其中 1， 2， N为差分方程为差分方程(8.2-2)或或(8.2-3)的特征方程的特征方程（8.2.8）特征根，且是单根。）特征根，且是单根。 系数C1，C2，CN是由差分方程边界条件决定的系数。 例例8.2-1 设描写某系统的齐次差分方程为 yn - 0.7yn-1 + 0.1yn-2 = 0， 并设 y-1= -26，y-2= -20

10、2， 求该差分方程的齐次解。齐次差分方程的齐次解为齐次差分方程的齐次解为h1122 nnnNNy nCCC(8.2-9)14差分方程的特征方程为： 2 0.7 + 0.1=0求得特征根为：1 = 0.2，2 = 0.5，于是齐次解为：yn = C1(0.2)n + C2(0.5)n将已知的边界条件y1= 26，y2= 202代入原方程式，通过迭代求出y0和y1：y0 = 0.7 y1 0.1 y2 = 18.2 + 20.2 = 2y1 = 0.7 y0 0.1 y1 = 1.4 + 2.6 = 4将y0和y1分别代入上述解中，得到一组联立方程12122 4 (0.2) (0.5)CCCC由此

11、求得系数, C1 = 10，C2 = 12，则方程的解为 yn = -10 (0.2)n + 12 (0.5)n解：解：8.2.1 8.2.1 线性常系数差分方程的时域经典法求线性常系数差分方程的时域经典法求解解15特征方程有重根特征方程有重根 假定1是特征方程式的K重根，那么，在齐次解中，相应于 1的部分将有K项 121121111KnKnnnKKC nC nCnC(8.2-10)例例8.2-2 求下述差分方程的齐次解。 yn -2yn-1 +2yn-2 -2yn-3+ yn-4 = 0已知边界条件为： y1=1, y2 =0, y3 =1, y5 =1。解解：特征方程为 4-23 +2 2

12、 -2 +1 = 0 即 ( -1)2(2 + 1) = 0特征根为 1 = 2 =1（二重根），3 = j, 4 = -j（共轭复根）8.2.1 8.2.1 线性常系数差分方程的时域经典法求线性常系数差分方程的时域经典法求解解16特征方程有重根特征方程有重根 j/2j/212341234 j( j)eennnny nC nCCCC nCCC12=cossin22nnC nCPQ 1 = C1+C2+Q 0 = 2C1+C2-P 1 = 3C1+C2-Q 1 = 5C1+C2+Q由上述方程组解得 C1 = 0, C2 = 1, P = 1, Q = 0从而有 1cos2ny n 8.2.1 8

13、.2.1 线性常系数差分方程的时域经典法求线性常系数差分方程的时域经典法求解解17 例例8.2-3 图8.2-1是一个链形电阻网络，设信号源的电压为x(t)，试确定输出端的电流y(t)（注：本例的链形网络不是一个离散时间系统，因为n不是时间变量，而是表示节点的序号，如前所述差分方程中变量的选取因具体函数而异，并不仅限于时间）。0( )vRy t102( )3( )vvRy tRy t解：根据该网络的结构和电阻值，设从右往左数的第n个节点的对地电压为vn，则有8.2.1 8.2.1 线性常系数差分方程的时域经典法求线性常系数差分方程的时域经典法求解解18 111222v nv nv nv nv

14、nRRR对2 n 100的节点，则有 4 120v nv nv n整理后得 2410 该二阶差分方程式的特征方程为 1223,23特征根为 12 (23)(23)nnv nCC于是节点电压的一般表达式为 8.2.1 8.2.1 线性常系数差分方程的时域经典法求线性常系数差分方程的时域经典法求解解1912120( ) 1(23)(23)3( )vCCRy tvCCRy t根据初始条件，有 123333( ) ( ) 66CRy tCRy t所以 1 (33)(23)(33)(23) ( )6nnv nRy t从而有 1001006( )( )(33)(23)(33)(23)x ty tR显然，当

15、n = 100时，有v100 = x(t)，于是求得8.2.1 8.2.1 线性常系数差分方程的时域经典法求线性常系数差分方程的时域经典法求解解20特解的求解特解的求解 (C0+C1n+C2n2+ Cr-1nr-1+Crnr )n（是方程的r重特征根）nCe nCn（不是方程的特征根）nC0+C1n+C2n2+ Ck-1nk-1+CknkC1sin n+C2cos nsin n（ cos n）C0+C1nA e j n（A为复数）e j nB（常数）特 解 形 式自 由 项特 解 形 式e n（为实数）n knC（常数）自 由项8.2.1 8.2.1 线性常系数差分方程的时域经典法求线性常系数

16、差分方程的时域经典法求解解21h ( 2)ny nC解：（1）齐次解为p12 y nD nD从而特解为 2 1 1y ny nx nx n2 x nn例例8.2-4 求差分方程的完全解。 11y ，且边界条件为其中激励信号为2 x nn21n 代入差分方程的右端，得自由项为（2）将其中，D1和D2为待定系数，代入原方程得12133221D nDDn8.2.1 8.2.1 线性常系数差分方程的时域经典法求线性常系数差分方程的时域经典法求解解22比较两端系数得到 122 /3,1/9DDp21 39y nn所以特解为hp21 ( 2)39ny ny ny nCn完全解为 11y 11( 2)2/3

17、 1/9C （3）将边界条件 代入上式，可得从而C = 8/9hp821 ( 2)939ny ny ny nn所以完全响应为 8.2.1 8.2.1 线性常系数差分方程的时域经典法求线性常系数差分方程的时域经典法求解解238.2.2 8.2.2 线性常系数差分方程的零输线性常系数差分方程的零输 入响应与零状态响应求解入响应与零状态响应求解1NnkkkC 自由响应齐次解 (8.2-11)强迫响应 = ypn = 特解 (8.2-12)p1 NnkkkCy n完全解 (8.2-13)完全响应完全响应 = 零输入响应 + 零状态响应248.2.2 8.2.2 线性常系数差分方程的零输线性常系数差分方

18、程的零输入响应与零状态响应求解入响应与零状态响应求解若系统的激励序列xn = 0，仅由系统的起始状态y-1, y-2, , y-N 引起的响应，称为离散系统的零输入响应，常用yzi n表示。零输入响应：零输入响应：若激励序列xn在n = 0时接入系统，则系统的零状态指 y-1 = y-2 = = y-N = 0，此时，仅由激励序列xn所引起的响应，称为离散系统的零状态响应，常用yzsn表示。零状态响应：零状态响应： yn = yzin + yzsn(8.2-14)258.2.2 8.2.2 线性常系数差分方程的零输线性常系数差分方程的零输入响应与零状态响应求解入响应与零状态响应求解pzizsp

19、111 NNNnnnkkkkkkkkky nCy nCCy n (8.2-17) 强迫响应自由响应零输入响应零状态响应例例8.2-5 已知某系统的差分方程为11 1 23y ny nu n试分别求下面两种起始状态下的完全响应。 （1）y1 = 0， （2）y1 = 1268.2.2 8.2.2 线性常系数差分方程的零输线性常系数差分方程的零输入响应与零状态响应求解入响应与零状态响应求解解：解： （1）齐次解为 h1 2ny nC12 23ny nC当n 0时，方程右端自由项为常数1/3，故可假设特解为D，将其代入差分方程，解得D = 2/3。从而完全解为zs112 323ny nynu n再将

20、y-1 = 0代入yn，可得到C = -1/3。所以完全响应（零状态响应）写为278.2.2 8.2.2 线性常系数差分方程的零输线性常系数差分方程的零输入响应与零状态响应求解入响应与零状态响应求解零输入响应为 zi1 2nynA当n 0时，方程右端自由项为常数1/3，故可假设ypn为D，将其代入差分方程，解得D = 2/3。再将y-1 = 1代入yzin，可得到A = 1/2。所以yzin 为zi1 1 2 2nynp1 2nzsynBy n零状态响应为 （2）齐次解为 h1 2ny nC288.2.2 8.2.2 线性常系数差分方程的零输线性常系数差分方程的零输入响应与零状态响应求解入响应

21、与零状态响应求解考虑零状态条件，即y-1 = 0代入yn，可得到B = -1/3。所以零状态响应为zs112 323nynu nzizs112 623ny nynynu n完全响应为：298.2.3 8.2.3 线性常系数差分方程的线性常系数差分方程的z变换法求解变换法求解00NMkrkra y nkb x nr(8.2-2)100( ) ( )NMklrkrklkra zY zy l zb zX z(8.2-19) x nx n u n输入信号为，, y-2, y-1，两边做单边z变换，得代数方程 ，起始状态为y-N, y-N + 1,10000 ( )( )NMklrkrklkrNNkkk

22、kkka zy l zb zY zX za za z(8.2-21) 301120.7 10.1 10.1 2( )10.70.1yzyyY zzz11222.62 (1.3)1210( )(0.5)(0.2)0.50.210.70.1zz zzzY zzzzzzz例例8.2-6 用z变换法重解例8.2-1的差分方程。解：对差分方程两边取单边z变换。根据时移性质可得 Y(z) - 0.7z-1Y(z) + y-1 + 0.1z-2Y(z) + z-1y-1 + y-2 = 0 代入起始条件y-1 = -26和y2 = 202，并进行部分分式展开求其逆变换，得到零输入响应为12(0.5)10(0

23、.2) nnu nyn = 8.2.3 8.2.3 线性常系数差分方程的线性常系数差分方程的z变换法求解变换法求解318.2.3 8.2.3 线性常系数差分方程的线性常系数差分方程的z变换法求解变换法求解111 1( )( )11yY zX zzz23133( )11(1)(1)zzzzY zzzzzz 1 y ny nnu n例例8.2-7 用z变换法求解差分方程， 其中已知y-1=1。解：解：对差分方程两边取单边z变换。根据时移性质可得 Y(z) - z-1Y(z) + y-1 = X(z) 2( )(1)zX zz代入和起始条件y1 = 1， 得到(1)(1)1 1 22n nn nu

24、nu nu nyn =328.3 8.3 离散系统的单位样值响应和系统函数离散系统的单位样值响应和系统函数8.3.1 8.3.1 单位样值响应单位样值响应定义定义：离散时间系统受单位样值信号n激励而产生的零状态响应，称为单位样值响应（也称为单位冲激响应），一般以hn表示。hn在离散时间系统中的作用，完全类似于连续系统中的由 (t)引起的冲激响应h(t)。)()(ttx系统连续系统连续系统：( )( )zsyth t x nn系统离散系统：离散系统： zsynh nhn的求法：的求法： 1、迭代法；、迭代法； 2、z变换法变换法338.3.1 8.3.1 单位样值响应单位样值响应例例8.3-1：

25、已知已知yn-1/3yn-1=xn, 试求其单位样值响应试求其单位样值响应 hn。对于因果系统，对于因果系统， h-1=0 , x-1=-1=0210 100113111013311212( )3311 1 ( )33nhhhhhhh nh nn1,01 330,0nnnh nu nn hn - 1/3hn-1 = n- 齐次解的形式yn - 1/3 yn-1 = xn解：解：348.3.1 8.3.1 单位样值响应单位样值响应 ) 3)(2(365365131)()(222212zzzzzzzzzzXzYzH 解：解：对上述差分方程取z变换 zXzzYzz)31 ()651 (2213222

26、121)3)(2(3)(2zzzzzzzzzH例例(补充补充): 已知因果系统的差分方程为已知因果系统的差分方程为 5 16 2 3 2y ny ny nx nx n求求hn。358.3.1 8.3.1 单位样值响应单位样值响应3222121)(zzzzzH11 2 2 3 22nnh nnu nu n 11 (2 32) 2nnnu n 11 (2 32) 12nnnnu n 1 (2 32) 1nnnu n368.3.1 8.3.1 单位样值响应单位样值响应0 mu nnm例例8.3-4 试利用线性时不变系统的特性，讨论离散系统对单位阶跃信号un的零状态响应单位阶跃响应gn与单位样值响应h

27、n之间的关系。解：解：由式(7.1-4)第一个等号，知 1nu nu nu n 另一方面，根据式(7.1-5)，得从而有 1h ng ng ng n (8.3-2)0 mg nh nm利用线性时不变系统的时不变特性和叠加性，有 (8.3-1)378.3.1 8.3.1 单位样值响应单位样值响应单位样值响应单位样值响应h(n)与系统稳定性及因果性的关系与系统稳定性及因果性的关系:离散稳定系统 nnh| |离散因果系统 0, 0nnhnunhnh或8.3.2 8.3.2 线性时不变系统的时域分析线性时不变系统的时域分析卷积和卷积和 mx nx mnm因为根据线性时不变系统的性质， 有 zsmynx

28、 m h nm x nh n(8.3-3)388.3.2 8.3.2 线性时不变系统的时域分析线性时不变系统的时域分析卷积和卷积和式(8.3-3)说明，对于线性时不变离散时间系统，其零状态响应是激励信号和其单位样值响应的卷积和。这一特点类似于连续时间系统零状态响应的意义，而且其计算方法也类似于卷积积分的四个步骤，即包括反四个步骤，即包括反褶、时移、相乘、求和褶、时移、相乘、求和。当然也可以用变换域方法实现。在7.1.3节我们已给出了卷积和的定义及其计算的一个简单示例。 例例1：某系统某系统hn=an un, 0a1 xn=un-un-N, 求求yn=xn*hn n mmmy nx mhn mu

29、 mu m N au n m398.3.2 8.3.2 线性时不变系统的时域分析线性时不变系统的时域分析卷积和卷积和hn或hmxn或xmn或mn或m1)当n0，hn恒为正值，且1 1 1 mmmph nph nph n递减恒定递增 （若pm=1为二阶极点，则hn线性增长）（2）pm 0时的情况相同。(3) 当pm为复数时，一对共轭复数极点对应于hn的一项振幅按nmp规律变化的正弦序列分量。558.3.5 8.3.5 离散时间系统的因果性和稳定性离散时间系统的因果性和稳定性离散系统稳定的充分必要条件是单位样值响应绝对可和。 nh n nnH zh nh n zZ Z nnnnH zh n zh

30、n z 令 则：1z nH zh n 离散系统稳定的充分必要条件是离散系统稳定的充分必要条件是H(z)的收敛域必须包含单位圆。的收敛域必须包含单位圆。568.3.5 8.3.5 离散时间系统的因果性和稳定性离散时间系统的因果性和稳定性因果系统的充分必要条件是单位样值响应是因果序列，即 , 0 0, 0h nnh nn从而系统函数的收敛域为 | z | R+下面讨论几种不同的系统。（1）因果系统（）因果系统（hn为因果信号）为因果信号）R+1)Re(z)Im(zj收敛域为zR+则则 H(z)的全部极点落在单位圆内。的全部极点落在单位圆内。因为收敛域应包含单位圆，则1R+因而系统也是稳定的。578

31、.3.5 8.3.5 离散时间系统的因果性和稳定性离散时间系统的因果性和稳定性（2）若）若hn为左边序列为左边序列收敛域为zR-则则 H(z)的全部极点落在单位圆外。的全部极点落在单位圆外。因而系统也是稳定的。收敛域为RzR+-则则 H(z)有些极点落在单位圆内，有些极点落在单位圆内，有些极点落在单位圆外。有些极点落在单位圆外。因为收敛域应包含单位圆，则1,1RR+-（3）若）若hn为双边序列为双边序列R+1)Re(z)Im(zjR-588.3.5 8.3.5 离散时间系统的因果性和稳定性离散时间系统的因果性和稳定性解：解：(1) 0.5 nh nu n因果、稳定5 . 01)Re(z)Im(

32、zj21)Re(z)Im(zj(3) 21nh nun (2) 2 nh nu n21)Re(z)Im(zj因果、非稳定非因果、稳定例例8.3-12：检验下列系统的因果性和稳定性。检验下列系统的因果性和稳定性。5 . 05 . 0)() 1 (zzzzH22)()2(zzzzH22)()3(zzzzH(4)( )0.52(0.5)(2)zH zzzz598.3.5 8.3.5 离散时间系统的因果性和稳定性离散时间系统的因果性和稳定性(4)( )0.52(0.5)(2)zH zzzz非因果、稳定0.51)Re(z)Im(zj22 0.5 213nnh nu nun 608.3.5 8.3.5 离

33、散时间系统的因果性和稳定性离散时间系统的因果性和稳定性例例 对于下列差分方程所表示的离散因果系统对于下列差分方程所表示的离散因果系统 0.2 1 0.24 2 1y ny ny nx nx n(1)求系统函数求系统函数H(z), 并说明它的收敛域及系统的稳定性；并说明它的收敛域及系统的稳定性；(2)求单位样值响应求单位样值响应hn；(3)当激励当激励 xn为单位阶跃序列时，求零状态响应为单位阶跃序列时，求零状态响应yzsn解：解：(1)将差分方程两边取z变换，得)()()(24. 0)(2 . 0)(121zXzzXzYzzYzzY于是)6 . 0)(4 . 0() 1()()()(zzzzz

34、XzYzH618.3.5 8.3.5 离散时间系统的因果性和稳定性离散时间系统的因果性和稳定性H(z)的极点为0.4和0.6，因为是因果系统，所以收敛域|z|0.6，因此该系统是稳定的。(2) 将H(z)/z展成部分分式，得到6 . 04 . 04 . 04 . 1)(zzzzzH 1.4(0.4)0.4( 0.6) nnh nu n(3)若激励xn=un,则1)(zzzX)6 . 0)(4 . 0)(1() 1()()()(2zzzzzzXzHzY2.080.930.1510.40.6zzzzzz 2.080.93(0.4)0.15( 0.6) nny nu n628.4 8.4 离散系统的

35、频响特性离散系统的频响特性 8.4.1 8.4.1 离散系统频率响应特性的意义离散系统频率响应特性的意义 离散系统的频率响应特性（简称频响特性），反映了离散系统在正弦序列激励下的稳态响应随频率变化的情况。 jn Tx nAe设式中：jAAe ()() jnTjTjTy nh nx nAeH eH ex n系统的稳态响应为：系统的稳态响应为：- 离散系统频响特性离散系统频响特性其中：其中：() ( )j Tj Tj Tnz enH eh n eH z （8.4-3）638.4.1 8.4.1 离散系统频率响应特性的意义离散系统频率响应特性的意义() ( )jjj nz enH eh n eH z

36、（8.4-4） 离散系统的频率响应是离散系统的频率响应是 的连续的周期函数，其的连续的周期函数，其周期是周期是2 。这是离散系统不同于连续系统的一个突出。这是离散系统不同于连续系统的一个突出特点。特点。()jH e- 幅频特性幅频特性- 相频特性相频特性( ) T 所以式（8.4-3）可改写为：( )()()jjjH eH ee （8.4-5）64|H(ej)| 2 c2 2 2 2658.4.2 8.4.2 频响特性的几何作图法频响特性的几何作图法频率响应也可以由H(z)在z平面上零极点分布来判断。12100121()()()()( )()()()()MrMrNNmmzzzzzzzzH zH

37、Hzpzpzpzp( )101()()|()|()MjrjjjrNjmmezH eHH eeep 分子分子：rjjrrezAe分母分母：mjjmmepB e668.4.2 8.4.2 频响特性的几何作图法频响特性的几何作图法幅频特性 j(e)H11MrrNmmAB= | H0 |(8.4-6)1Mrr1Nmm () = (8.4-7)相频特性( )101()()|()|()MjrjjjrNjmmezH eHH eeep 678.4.2 8.4.2 频响特性的几何作图法频响特性的几何作图法图8.4-2 离散系统频率响应H(ej)的几何确定法Show688.4.2 8.4.2 频响特性的几何作图法

38、频响特性的几何作图法例例8.4-1 求下图所示一阶因果离散系统的频率响应。求下图所示一阶因果离散系统的频率响应。xna1zyn解解 : 该一阶系统的差分方程为 1 y nay nx n 它是因果系统，其系统函数为( )zH zza频率响应为1()(1cos )sinjjjeH eeaaja698.4.2 8.4.2 频响特性的几何作图法频响特性的几何作图法幅频特性为j(e)H2112 cosaa= sin1cosaa () = arctan相频特性为（1）当 0 a 1时：Show708.4.2 8.4.2 频响特性的几何作图法频响特性的几何作图法（2）当 1 a 0时：j(e)H2112 c

39、osaa= sin1cosaa () = arctanShow718.4.2 8.4.2 频响特性的几何作图法频响特性的几何作图法例例8.4-2 分析图8.4-5所示二阶离散系统的幅频特性。解：解： 0.9 10.81 21y ny ny nx n81. 09 . 081. 09 . 01)(2211zzzzzzzH323119 . 0,9 . 0, 0jjepepz728.4.2 8.4.2 频响特性的几何作图法频响特性的几何作图法112jjj()j1jj /3jj /3e12e(e )( )e(e0.9e)(e0.9e)zAHH zBB j(e )H112AB B121B B112( )=

40、() Show738.4.2 8.4.2 频响特性的几何作图法频响特性的几何作图法例例8.4-3 证明图8.4-7(a)所示系统为二阶全通离散系统。其中bl, b2为实数，且满足b12 4b2 0。 证明证明：图示系统的差分方程为 yn + b1yn 1 + b2yn 2 = b2xn+ b1xn 1 + xn 2748.4.2 8.4.2 频响特性的几何作图法频响特性的几何作图法yn + b1yn 1 + b2yn 2 = b2xn+ b1xn 1 + xn 2zbzbzbbzzbzbzzbbzH112211221121121)(1221122122()coscos sinsincoscos

41、 sinsinjjjjjebb eH eb ebebbj bbbj b()1jH e-全通系统全通系统758.4.2 8.4.2 频响特性的几何作图法频响特性的几何作图法122212112212121( )1bb zzb zb zH zb zb zzb zb全通系统的特点：21111()()()()b zzzzzpzp212(40)bb2121142bjbbp21211242bjbbzb212|pb2121|zb768.4.2 8.4.2 频响特性的几何作图法频响特性的几何作图法11,jjprepre1111,jjzezerr令： 778.5 8.5 数字滤波器的一般概念数字滤波器的一般概念8

42、.5.1 8.5.1 数字滤波器原理数字滤波器原理数字滤波器数字滤波器 hn或或H(z) x n y n y nh nx n)()()(zXzHzY()()()jjjY eH eX e（8.5-1）输入序列xn的频谱X(e j)经过滤波器后变为输出序列yn的频谱Y(e j)。 788.5.1 8.5.1 数字滤波器原理数字滤波器原理数字滤波器对连续时间信号处理的一般过程798.5.2 8.5.2 数字滤波器结构数字滤波器结构解：解：（1）直接形式）直接形式例例8.5-1：分别用直接形式、级联形式和并联形式模拟系统：32132353142)(zzzzzzH)(zY3)(zX1z51z1z4238

43、08.5.2 8.5.2 数字滤波器结构数字滤波器结构)321)(1 ()2(2353142)(21121132132zzzzzzzzzzzzH（2）级联形式）级联形式11112)(zzzH212123212)(zzzzzH)(zY)(zX31z221z21z1818.5.2 8.5.2 数字滤波器结构数字滤波器结构2311212311224( )1 3531123zzzzzH zzzzzzz（3）并联形式）并联形式,1)(111zzzH21212321)(zzzzzH)(zY)(zX31z1z21z11828.6 8.6 应用应用MATLABMATLAB分析离散时间系统分析离散时间系统fil

44、ter- 计算对于指定时间范围的激励序列的响应conv-求两个有限时间区间非零的离散时间序列卷积和freqz-分析系统频率响应tf2zp，zp2tf-提供系统函数的两种不同形式的转换zplane-可以绘制零极点分布图的函数impz- 计算系统单位样值响应的函数838.6 8.6 应用应用MATLABMATLAB分析离散时间系统分析离散时间系统0.2z 4j10.8ep4j20.8ep例例8.6-1 已知离散时间系统的系统函数的零极点分别为：绘制系统的零极点分布图，并绘出系统的单位样值响应hn的时域波形。解： z=0.2;p=0.8*exp(pi*i/4);0.8*exp(-pi*i/4);k=1;subplot(121);zplane(z,p);a,b=zp2tf(z,p,k);subplot(122);impz(a,b,20);title(hn); 848.6 8.6 应用应用MATLABMATLAB分析离散时间系统分析离散时间系统858.6 8.6 应用应用MATLABMATLAB分析离散时间系统分析离散时间系统2 3 12 2 1y ny ny nx nx n 4sin 5nx nu n例8.6-2 已知