内容发布更新时间 : 2024/5/21 22:26:38星期一 下面是文章的全部内容请认真阅读。

实 验 报 告

课程名称： 数字信号处理实验 专业班级： 姓 名： 学 号：

实验名称 数字信号处理 实验地点 实验成绩 一、实验目的及任务 ? ? ? ? 学会运用MATLAB求离散时间信号的z变换和z反变换； 学会运用MATLAB分析离散时间系统的系统函数的零极点； 学会运用MATLAB分析系统函数的零极点分布与其时域特性的关系； 学会运用MATLAB进行离散时间系统的频率特性分析。 实验时间 二、实验内容与步骤 4.2.1 z正反变换 序列x?n?的z变换定义为 X?z??Z?x?n???n????x?n?z??n （4-1） 其中，符号Z表示取z变换，z是复变量。相应地，单边z变换定义为 X?z??Z?x?n????x?n?z?n （4-2） n?0?MATLAB符号数学工具箱提供了计算离散时间信号单边z变换的函数ztrans和z反变换函数iztrans，其语句格式分别为 Z=ztrans(x) x=iztrans(z) 上式中的x和Z分别为时域表达式和z域表达式的符号表示，可通过sym函数来定义。 如果信号的z域表示式X(z)是有理函数，进行z反变换的另一个方法是对X(z)进行部分分式展开，然后求各简单分式的z反变换。设X(z)的有理分式表示为 b0?b1z?1?b2z?2???bmz?mB(z) （4-3） X(z)??A(z)1?a1z?1?a2z?2???anz?nMATLAB信号处理工具箱提供了一个对X(z)进行部分分式展开的函数residuez，其语句格式为 [R,P,K]=residuez(B,A) 其中，B，A分别表示X(z)的分子与分母多项式的系数向量；R为部分分式的系数向量；P为极点向量；K为多项式的系数。若X(z)为有理真分式，则K为零。 4.2.2 系统函数的零极点分析 离散时间系统的系统函数定义为系统零状态响应的z变换与激励的z变换之比，即 H(z)?如果系统函数H(z)的有理函数表示式为 Y(z) （4-4） X(z)b1zm?b2zm?1???bmz?bm?1 （4-5） H(z)?a1zn?a2zn?1???anz?an?1那么，在MATLAB中系统函数的零极点就可通过函数roots得到，也可借助函数tf2zp得到，tf2zp的语句格式为 [Z,P,K]=tf2zp(B,A) 其中，B与A分别表示H(z)的分子与分母多项式的系数向量。它的作用是将H(z)的有理分式表示式转换为零极点增益形式，即 H(z)?k(z?z1)(z?z2)?(z?zm) （4-6） (z?p1)(z?p2)?(z?pn)若要获得系统函数H(z)的零极点分布图，可直接应用zplane函数，其语句格式为 zplane(B,A) 其中，B与A分别表示H(z)的分子和分母多项式的系数向量。它的作用是在Z平面上画出单位圆、零点与极点。 4.2.3 系统函数的零极点分布与其时域特性的关系 与拉氏变换在连续系统中的作用类似，在离散系统中，z变换建立了时域函数h(n)与z域函数H(z)之间的对应关系。因此，z变换的函数H(z)从形式可以反映h(n)的部分内在性质。我们仍旧通过讨论H(z)的一阶极点情况，来说明系统函数的零极点分布与系统时域特性的关系。 4.2.4 离散时间LTI系统的频率特性分析