实验三Z变换.pptVIP

一、实验目的 1、学会运用MATLAB求离散时间信号的z变换和z反变换； 2、学会运用MATLAB分析离散时间系统的系统函数的零极点分布与时频特性分析； 二、实验设备 三、实验原理 三、实验原理 三、实验原理 三、实验原理 三、实验原理 三、实验原理 如果信号的z域表示式是有理函数，进行z反变换的另一个方法是对进行部分分式展开，然后求各简单分式的z反变换.如果X(z)的有理分式表示为： 三、实验原理 三、实验原理 【例3-3】试用MATLAB命令进行部分分式展开，并求出其z反变换。 三、实验原理 三、实验原理 三、实验原理 三、实验原理 三、实验原理 三、实验原理 三、实验原理 三、实验原理 三、实验原理 三、实验原理 三、实验原理 三、实验原理 三、实验原理 三、实验原理 三、实验原理 三、实验原理 三、实验原理 三、实验原理 三、实验原理 三、实验原理 四、实验内容 四、实验内容 四、实验内容 四、实验内容 四、实验内容 五、实验报告要求 【1】 试用ztrans函数求下列函数的z变换。 zx=sym(2^(n-1)-(-2)^(n-1)); Z=ztrans(x); simplify(Z) ans= z^2/(z-2)/(z+2) 【2】 试用iztrans函数求下列函数的z反变换。 Z=sym(z*(2*z^2-11*z+12)/(z-1)/(z-2)^3); x=iztrans(Z); simplify(x) ans= -3+3*2^n-1/4*2^n*n-1/4*2^n*n^2 【3】试用MATLAB的residuez函数，求出X(z)的部分分式展开和。 【4】试用MATLAB画出下列因果系统的系统函数零极点分布图,以及对应的时域单位抽样响应h(n)的波形。 【5】试用MATLAB绘制以下系统的频率响应曲线。 记录实验主要程序及结果。 * * 实验三 z变换 2、MATLAB6.5 软件 1、计算机 (1)??? 序列的正反Z变换 ???其中，符号表示取z变换，z是复变量。相应地，单边z变换定义为： MATLAB实现 MATLAB符号数学工具箱提供了计算离散时间信号单边z变换 的函数ztrans和z反变换函数iztrans，其语句格式分别为 Z=ztrans(x) x=iztrans(z) 上式中的x和Z分别为时域表达式和z域表达式的符号表示， 可通过sym函数来定义。 【例3-1】 试用ztrans函数求下列函数的z变换。 x=sym(a^n*cos(pi*n)); Z=ztrans(x); simplify(Z) ans= z/(z+a) 【例3-2】 试用iztrans函数求下列函数的z反变换。 Z=sym((8*z-19)/(z^2-5*z+6)); x=iztrans(Z); simplify(x) ans= -19/6*charfcn[0](n)+5*3^(n-1)+3*2^(n-1) -19/6*charfcn[0](n)+5*3^(n-1)+3*2^(n-1) charfcn[0](n)是函数在MATLAB符号工具箱中的表示， 反变换后的函数形式为： MATLAB信号处理工具箱提供了一个对进行部分分式展开的函数residuez，其语句格式为： [R,P,K]=residuez(B,A) 其中，B，A分别表示X(z)的分子与分母多项式的系数向量，分子与分母多项式按照z^-1升幂排列。 R为部分分式的系数向量； P为极点向量； K为多项式的系数。 若X(z)为有理真分式，则K为零。 解：MATLAB源程序为 B=[18]; A=[18,3,-4,-1]; [R,P,K]=residuez(B,A) R= 0.3600 0.2400 0.4000 P= 0.5000 -0.3333 -0.3333 K= [] 从运行结果可知， 表示系统有一个二重极点。 所以，X(z)的部分分式展开为 2、系统函数的零极点分析 离散时间系统的系统函数定义为系统零状态响应的z变换与 激励的z变换之比: 如果系统函数 的有理函数表示式为 在MATLAB中系统函数的零极点就可通过函数roots得到，也可 借助DSP工具箱中的函数tf2zp得到，tf2zp的语句格式为： [R,P,K]=tf2zp(B,A) 其中，B与A分别表示分子与分母多项式的系数向量。 它的作用是将H(z)的有理分式表示式转换为零极点增益形式： MATLAB实现 【实例3-4】 已知一离散因果LTI系统的系统函数为： 试用MATLAB命令求该系统的零极点。 。 B=[1,0.3