山科大信号和系统实验六.docVIP

实验 离散信号与系统的Z变换分析一、?实验目的 1熟悉离散信号Z变换的原理及性质 2熟悉常见信号的Z变换 3了解正/反Z变换的MATLAB实现方法 4了解离散信号的Z变换与其对应的理想抽样信号的傅氏变换和拉氏变换之间的关系 5了解利用MATLAB实现离散系统的频率特性分析的方法 二、?实验原理 1正/反Z变换 Z变换分析法是分析离散时间信号与系统的重要手段。如果以时间间隔对连续时间信号f(t)进行理想抽样，那么，所得的理想抽样信号为： 理想抽样信号的双边拉普拉斯变换F? (s)为： 若令 ，那么的双边拉普拉斯变换F (s)为： 则离散信号f(k)的Z变换定义为：? 从上面关于Z变换的推导过程中可知，离散信号f(k)的Z变换F(z)与其对应的理想抽样信号的拉氏变换F? (s)之间存在以下关系： ? 同理，可以推出离散信号f(k)的Z变换F(z)和它对应的理想抽样信号的傅里叶变换之间的关系为 如果已知信号的Z变换F(z)，要求出所对应的原离散序列f(k)，就需要进行反Z变换，反Z变换的定义为： ????? 其中，C为包围的所有极点的闭合积分路线。 在MATLAB语言中有专门对信号进行正反Z变换的函数ztrans( ) 和itrans( )。其调用格式分别如下： ????????? F=ztrans( f ) ?????对f(n)进行Z变换，其结果为F(z) ????????? F=ztrans(f,v)? ?对f(n)进行Z变换，其结果为F(v) ????????? F=ztrans(f,u,v) ???对f(u)进行Z变换，其结果为F(v) ????????? f=itrans ( F )? ????对F(z)进行Z反变换，其结果为f(n) ????????? f=itrans(F,u)? ?对F(z)进行Z反变换，其结果为f(u) ????????? f=itrans(F,v,u ) ???对F(v)进行Z反变换，其结果为f(u) 注意： 在调用函数ztran( )及iztran( )之前，要用syms命令对所有需要用到的变量（如t,u,v,w）等进行说明，即要将这些变量说明成符号变量。 例① 用MATLAB求出离散序列? 的Z变换 syms k z f=0.5^k Fz=ztrans(f) 运行结果 ：Fz =2*z/(2*z-1) 例已知一离散信号的Z变换式为 ，求出它所对应的离散信号f(k)syms k z FZ=2*z/(2*z-1) fk=iztrans(FZ,k) 运行结果 ：fk =(1/2)^k 2.离散系统的频率特性 同连续系统的系统函数H(s)类似，离散系统的系统函数H(z)也反映了系统本身固有的特性。对于离散系统来说，如果把其系统函数H(z)中的复变量z换成，那么所得的函数就是此离散系统的频率响应特性，即离散时间系统的频率响应为： 其中， 称为离散系统的幅频特性，称为系统的相频特性。同连续系统一样，离散时间系统的幅频特性也是频率的偶函数，相频特性也是频率的齐函数。 由于是频率的周期函数，所以离散系统的频率响应特性也是频率的周期函数，其周期为，或者频率周期为。实际上，这就是抽样系统的抽样频率，而其中的T则是系统的抽样周期。频率响应呈现周期性是离散系统特性区别于连续系统特性的重要特点。因此，只要分析在范围内的情况，便可分析出系统的整个频率特性。 鉴于离散系统频率响应特性的特点，为了表示方便起见，我们通常将其中的用一个变量来代替，即令代入系统函数H(z)中，用函数来表示离散系统的频率响应特性。相应地，用表示幅频特性，而相频特性仍用来表示。应该特别注意的是，虽然这里的变量?仍然称为频率变量，但是它已经不是原来意义上的角频率概念，而实际上是表示角度的概念。我们称之为数字角频率。它与原来角频率的关系为：。也就是说，根据离散系统的系统函数H(z)，令其中的，并且代入0～范围内不同的频率值（实际上是角度值），就可以逐个计算出不同频率时的响应，求出离散系统的频率响应特性。再利用离散系统频率特性的周期性特点（周期为），求出系统的整个频率特性。 离散系统的幅频特性曲线和相频特性曲线能够直观地反映出系统对不同频率的输入序列的处理情况。在函数随的变换关系中，在=0附近，反映了系统对输入信号低频部分的处理情况，而在附近，则反映了系统对输入信号高频部分的处理情况。 一般来说，分析离散系统频率响应特性就要绘制频率响应曲线，而这是相当麻烦的。虽然可以通过几何矢量法来定性画出频率响应特性曲线，但一般来说这也是很麻烦的。值得庆幸的是，MATLAB为我们提供了专门用于求解离散系统频率响应的函数frez() ，其调用格式如下： ??????[H,w]=frez(B,A,N)??? 其中，B和A分别是表示待分析的离散