系统的复频域分析实验报告.docVIP

实验六 系统的复频域分析 信号的拉普拉斯变换 （6.1） 是连续时间傅立叶变换地推广。连续时间傅立叶变换在研究连续时间信号与系统中是很有用的。然而，许多信号不存在傅立叶变换而存在拉普拉斯变换，这使得拉普拉斯变换成为线性时不变系统分析的一种有用方法。对一大类信号来说，它们的拉普拉斯变换可以表示为s的多项式之比，即 这里和分别称作分子和分母多项式。能表示成多项式之比的变换称为有理变换，这里作为满足线性常系数微分方程的LTI系统的系统函数中常常出现。除了一个标量因子外，有理变换是完全由多项式和的根决定的，这些根分别称为零点和极点。由于这些根在LTI系统的研究中起着重要的作用，所以它们以零极点图的方式展现出来的是很方便的。这一章将用拉普拉斯变换在复频域研究LTI系统的一些性质。 §6.1 MATLAB函数lsim（用于系统函数） 目的 用lsim仿真由系统函数表征的因果LTI系统的输出。 相关知识 第二章所讨论的是如何用lsim命令仿真一个输出满足一个线性常系数微分方程的因果LTI连续时间系统。因为系统函数唯一地表征了关联系统输入和输出的微分方程。所以由系统函数表征的因果LTI系统的输出也能够用lsim仿真。如果系统函数给出如下形式： （6.2） 那么，对输入地系统的系统输出就能用lsim(b,a,x,t)仿真，其中MATLAB向量b和a包含了分子分母s多项式的系数。 例如，考虑系统函数，其系数由向量b=[1 1/2]和a=[1 -2]定义。命令lsim(b,a,x,t)将系统对由向量x给出的输入，在向量t给定的时刻上的系统的时间响应应存入向量y中，向量y和输入向量x有相同数量的元素。 基本题 1．定义系数向量a1和b1用以描述由下面系统函数表征的因果LTI系统： a1=[1 -2]; %分子系数 b1=[1 2]; %分母系数 2．定义系数向量a2和b2用以描述由下面系统函数表征的因果LTI系统： a2=[0 3]; %分子系数 b2=[1 0.3]; %分母系数 3．定义系数向量a3和b3用以描述由下面系统函数表征的因果LTI系统： a3=[2 0]; %分子系数 b3=[1 0.8]; %分母系数 4．利用lsim和前面部分定义的向量求这些因果LTI系统对由t=[0:0.1:0.5]，x=cos(t)给出的输入的输出。 代码： a1=[1 -2]; b1=[1 2]; a2=[0 3]; b2=[1 0.3]; a3=[2 0]; b3=[1 0.8]; t=[0:0.1:0.5] x=cos(t); y1=lsim(a1,b1,x,t); y2=lsim(a2,b2,x,t); y3=lsim(a3,b3,x,t); [t,x,y1,y2,y3] 结果： ans = 0 1.0000 1.0000 0 2.0000 0.1000 0.9950 0.6334 0.2948 1.8366 0.2000 0.9801 0.3261 0.5779 1.6667 0.3000 0.9553 0.0692 0.8468 1.4910 0.4000 0.9211 -0.1444 1.0991 1.3105 0.5000 0.8776 -0.3205 1.3323 1.1262 §6.2作连续时间的零极点图 目的 这一节要学习如何在一个零极点图上展现有理系统函数的零极点。 相关知识 一个有理系统函数的零极点能用函数roots计算出。例如，对于系统函数为 的LTI系统，其零极点可用执行如下命令而计算出： b=[1 -1]; a=[1 3 2]; zs=roots(b) zs = 1 ps=roots(a) ps = -2 -1 一个简单的零极点图可以用在复数s平面内，在每个极点的位置放一个x，在每个零点的位置放一个?来完成，也即 plot(real(zs),imag(zs),o); hold on plot(r