matlab z域响应,实验七离散时间信号与系统的z域分析及matlab实现.doc

实验七离散时间信号与系统的z域分析及matlab实现.doc

实验七离散时间信号与系统的Z域分析及MATLAB实现实验目的1．熟悉信号Z变换的实现2．熟悉离散时间系统频率响应的实现3．熟悉离散时间系统零极点分布图和系统幅频响应的MATLAB实现实验原理与步骤一．反Z变换的MATLAB实现在MATLAB中，函数RESIDUEZ可以用来计算一个有理分式的留数部分和直接项。设多项式表示为函数RESIDUEZ的调用格式为R,P,CRESIDUEZB,A其中，向量B，A分别决定多项式BZ和AZ的系数，按Z的降幂排列。返回的列向量R包含留数值，列向量P包含极点的位置，行向量包含直接项。一般情况下，用MATLAB求反Z变换的步骤为(1)确定Z变换的收敛域，这是最关键的一步。(2)用MATLAB求出有理分式的极点分布，以及在各个极点上的留数值和直接项的值，从而将有理分式分解为简单的分式之和，然后根据收敛域，求出原始序列。解本例的反Z变换的程序如下B10,1,1A1,0,025确定Z变换表示式R,P,CRESIDUEZB,A运算结果显示为R42P0500005000C4因此，分式XZ可以表示为下面的MATLAB命令将产生频率响应曲线NUMDNDN1D0DENANAN1A0OMEGAPIPI/150PIHFREQZNUM,DEN,OMEGAMAGABSHPHASE180/PIUNWRAPANGLEH则下面的MATLAB程序可以获得频率响应曲线NUM10DEN105OMEGAPIPI/150PIHFREQZNUM,DEN,OMEGASUBPLOT211,PLOTOMEGA,ABSHSUBPLOT212,PLOTOMEGA,180/PIUNWRAPANGLEH此程序产生的曲线如下图1所示三．离散时间系统零极点分布图和系统幅频响应的MATLAB实现试用MATLAB在Z平面中画出HZ的零点和极点，以及系统的幅度响应。解本例的MATLAB程序实现为B1,1,2A1,15,1FIGURESUBPLOT221ZPLANEB,AXLABEL 虚部 YLABEL 实部 TITLE 零极点图 H,WFREQZB,A,250FIGURESUBPLOT222PLOTW,ABSHXLABEL 频率 YLABEL 幅度 TITLE 幅频响应图 系统的零极点图如图2(A)，系统幅频响应曲线如图2(B)。