拉氏变换与反变换

2.5拉氏变换与反变换机电控制工程所涉及的数学问题较多，经常要解算一些线性微分方程。按照一般方法解算比较麻烦，如果用拉普拉斯变换求解线性微分方程，可将经典数学中的微积分运算转化为代数运算，又能够单独地表明初始条件的影响，并有变换表可查找，因而是一种较为简便的工程数学方法。2.5.1拉普拉斯变换的定义如果有一个以时间为自变量的实变函数，它的定义域是，那么的拉普拉斯变换定义为(2.10)式中，是复变数，（σ、ω均为实数），称为拉普拉斯积分；是函数的拉普拉斯变换，它是一个复变函数，通常也称为的象函数，而称为的原函数；L是表示进行拉普拉斯变换的符号。式（2.10）表明：拉氏变换是这样一种变换，即在一定条件下，它能把一实数域中的实变函数变换为一个在复数域内与之等价的复变函数。2.5.2几种典型函数的拉氏变换1.单位阶跃函数的拉氏变换单位阶跃函数是机电控制中最常用的典型输入信号之一，常以它作为评价系统性能的标准输入，这一函数定义为单位阶跃函数如图2.7所示，它表示在时刻突然作用于系统一个幅值为1的不变量。单位阶跃函数的拉氏变换式为当，则。所以ttf0ttf0edstFsLftfttsjs0est)(sF)(tf)(sF)(tf)(tf)(sF)(sF)(1t)0(1)0(0)(1ttt0t0e1de)(1)](1[)(0stststttLsF0)Re(s0elimstt（2.11）图2.7单位阶跃函数2.指数函数的拉氏变换指数函数也是控制理论中经常用到的函数，其中是常数。令则与求单位阶跃函数同理，就可求得（2.12）3.正弦函数与余弦函数的拉氏变换设，，则由欧拉公式，有所以ssstLst1)1(00e1)(1ttsFsttsttdeedeej21)(0j0j1(2.13）同理(2.14）4.单位脉冲函数δ(t)的拉氏变换单位脉冲函数是在持续时间期间幅值为的矩形波。其幅值和作用时间的乘积等于1，即。如图2.8所示。图2.8单位脉冲函数单位脉冲函数的数学表达式为其拉氏变换式为此处因为时，，故积分限变为。ttsttstsdeedej210)j(0)j(0ej10ej1j21)j()j(tstsss22j1j1j21sss(2.15)5.单位速度函数的拉氏变换单位速度函数，又称单位斜坡函数，其数学表达式为见图2.9所示。图2.9单位速度函数单位速度函数的拉氏变换式为利用分部积分法令则所以当时，,则（2.16）6.单位加速度函数的拉氏变换单位加速度函数的数学表达式为如图2.10所示图2.10单位加速度函数其拉氏变换式为（2.17）2.5.3拉氏变换的主要定理根据拉氏变换定义或查表能对一些标准的函数进行拉氏变换和反变换，但利用以下的定理，则对一般的函数可以使运算简化。1.叠加定理拉氏变换也服从线性函数的齐次性和叠加性。（1）齐次性设，则（2.18）式中——常数。（2）叠加性设,，则（2.19）两者结合起来，就有这说明拉氏变换是线性变换。2.微分定理设则式中——函数在时刻的值，即初始值。同样，可得的各阶导数的拉氏变换是（2.20）式中,，…——原函数各阶导数在时刻的值。如果函数及其各阶导数的初始值均为零（称为零初始条件），则各阶导数的拉氏变换为（2.21）3.复微分定理若可以进行拉氏变换，则除了在的极点以外，（2.22）式中，。同样有一般地，有（2.23）4.积分定理设，则（2.24）式中——积分在时刻的值。当初始条件为零时，（2.25）对多重积分是（2.26）当初始条件为零时，则（2.27）5.延迟定理设，且时，，则（2.28）函数为原函数沿时间轴延迟了，如图2.11所示。sFsttfLddsFstftL222ddd11,2,3,dnnnnLtftFsns图2.11函数6.位移定理在控制理论中，经常遇到一类的函数，它的象函数只需把用代替即可，这相当于在复数坐标中，有一位移。设，则（2.29）例如的象函数，则的象函数为7.初值定理它表明原函数在时的数值。（2.30）即原函数的初值等于乘以象函数的终值。8.终值定理设，并且存在，则（2.31）即原函数的终值等于乘以象函数的初值。这一定理对于求瞬态响应的稳态值是很有用的。9.卷积定理设，，则有（2.32）即两个原函数的卷积分的拉氏变换等于它们象函数的乘积。式（2.32）中，为卷积分的数学表示，定义为10.时间比例尺的改变（2.33）式中——比例系数例如,的象函数,则的象函数为11.拉氏变换的积分下限在某些情况下,在处有一个脉冲函数。这时必须明确拉普拉斯积分的下限是还是，因为对于这两种下限，的拉氏变换是不同的。为此，可采用如下符号予以区分：若在处包含一个脉冲函数，则因为在这种情况下显然，如果在处没有脉冲函数，则有2.5.4拉普拉斯反变换拉普拉斯反变换的公式为（2.36）jj1de)(πj21)(ccstssFsFLtf式中——表示拉普拉斯反变换的符号通常用部分分式展开法将复杂函数展开成有理分式函数之和，然后由拉氏变换表一一查出对应的反变换函数，即得所求的原函数。1.部分分式展开法在控制理论中，常遇到的象函数是的有理分式为了将写成部分分式，首先将的分母因式分解，则有式中，，，…，是的根的负值，称为的极点，按照这些根的性质，可分为以下几种情况来研究。2.的极点为各不相同的实数时的拉氏反变换（2.37）式中，是待定系数，它是处的留数，其求法如下（2.38）再根据拉氏变换的迭加定理，求原函数[例2.1]求的原函数。解:首先将的分母因式分解，则有1L)(tf即得3.含有共轭复数极点时的拉氏反变换如果有一对共轭复数极点,，其余极点均为各不相同的实数极点。将展成式中,和可按下式求解即（2.39）因为(或）是复数，故式（2.39）两边都应是复数，令等号两边的实部、虚部分别相等，得两个方程式，联立求解，即得，两个常数。[例2.2]已知,试求其部分分式。解:因为（2.40）含有一对共轭复数极点,和一个极点，故可将式（2.40）因式分解成（2.41）以下求系数、和。由式（2.40）和式（2.41）相等，有（2.42）用乘以上式两边，并令，得到上式可进一步写成由上式两边实部和虚部分别相等，可得联立以上两式，可求得为了求出系数,用乘方程（2.42）两边，并令,将代入，得![endif]将所求得的,,值代入（2.41），并整理后得的部分分式查拉氏变换表便得，结果见式（3.16）。[例2.3]已知求。解:将的分母因式分解，得利用方程两边实部、虚部分别相等得解得，所以这种形式再作适当变换：查拉氏变换表得4.中含有重极点的拉氏反变换设有r个重根，则将上式展开成部分分式（2.43）式中，，，…，的求法与单实数极点情况下相同。，，…，的求法如下：……则(2.44）[例2.4]设，试求的部分分式。解:已知（2.45）含有2个重极点，可将式（2.45）的分母因式分解得（2.46）以下求系数、和。将所求得的、、值代入式（2.46），即得的部分分式查拉氏变换表可得。[例2.5]求的拉氏反变换。解:将展开为部分分式上式中各项系数为于是查拉氏变换表，得5.用MATLAB展开部分分式(1)概述MATLAB是美国MathWorks公司的软件产品，是一个高级的数值分析、处理与计算的软件，其强大的矩阵运算能力和完美的图形可视化功能，使得它成为国际控制界应用最广的首选计算机工具。SIMULINK是基于模型化图形的动态系统仿真软件，是MATLAB的一个工具箱，它使系统分析进入一个崭新的阶段，它不需要过多地了解数值问题，而是侧重于系统的建模、分析与设计。其良好的人机界面及周到的帮助功能使得它广为科技界和工程界所采用。(2)用MATLAB进行部分分式展开MATLAB有一个命令用于求B(s)/A(s)的部分分式展开式。设s的有理分式为式中（i=）和(j=）的某些值可能为零。在MATLAB的行向量中，num和den分别表示F(s)分子和分母的系数，即num=[]den=[1]命令[r,p,k]=residue(num,den)MATLAB将按下式给出F(s)部分分式展开式中的留数、极点和余项：上式与式（2.37）比较，显然有p(1)=-p1，p(2)=-p2，…，p(n)=-pn；r(1)=A1,r(2)=A2，…，r(n)=An；k（s）是余项。[例2.6]试求下列函数的部分分式展开式解：对此函数有num=[111395226]den=[110355024]命令[r,p,k]=residue(num,den)于是得到下列结果[r,p,k]=residue(num,den)r=1.00002.5000-3.00000.5000p=-4.0000-3.0000-2.0000-1.0000k=1则得如果F(s)中含重极点，则部分分式展开式将包括下列诸项式中，p(j)为一个q重极点。[例2.7]试将下列函数展开成部分分式解：对于该函数有num=[0146]den=[1331]命令[r,p,k]=residue(num,den)将得到如下结果：[r,p,k]=residue(num,den)r=1.00002.00003.0000p=-1.0000-1.0000-1.0000k=[]所以可得注意，本例的余项k为零。2.5.5应用拉氏变换解线性微分方程应用拉氏变换解线性微分方程时，采用下列步骤：(1)对线性微分方程中每一项进行拉氏变换，使微分方程变为的代数方程；(2)解代数方程，得到有关变量的拉氏变换表达式；(3)用拉氏反变换得到微分方程的时域解。整个求解过程如图2.12所示。2.8设系统微分方程为若，初始条件分别为、，试求。解:对微分方程左边进行拉氏变换利用迭加定理将上式逐项相加，即得方程左边的拉氏变换对方程右边进行拉氏变换得写成一般形式应该强调指出是微分方程的特征方程，也是该系统的特征方程。利用部分分式将展开为s求待定系数、、、、:代入原式得查拉氏变换表得当初始条件为零时，得