拉氏变换与反变换.doc

word2.5 拉氏变换与反变换机电控制工程所涉与的数学问题较多，经常要解算一些线性微分方程按照一般方法解算比拟麻烦，如果用拉普拉斯变换求解线性微分方程，可将经典数学中的微积分运算转化为代数运算，又能够单独地明确初始条件的影响，并有变换表可查找，因而是一种较为简便的工程数学方法2.5.1 拉普拉斯变换的定义如果有一个以时间 为自变量的实变函数 ，它的定义域是，那么的拉普拉斯变换定义为(2.10) 式中，是复变数，〔σ、ω均为实数〕，称为拉普拉斯积分；是函数的拉普拉斯变换，它是一个复变函数，通常也称 为的象函数，而称为 的原函数；L是表示进展拉普拉斯变换的符号 式〔2.10〕明确：拉氏变换是这样一种变换，即在一定条件下，它能把一实数域中的实变函数变换为一个在复数域内与之等价的复变函数 2.5.2 几种典型函数的拉氏变换1.单位阶跃函数 的拉氏变换单位阶跃函数是机电控制中最常用的典型输入信号之一，常以它作为评价系统性能的标准输入，这一函数定义为单位阶跃函数如图2.7所示，它表示在 时刻突然作用于系统一个幅值为1的不变量单位阶跃函数的拉氏变换式为当，如此所以〔2.11〕图2.7 单位阶跃函数的拉氏变换指数函数也是控制理论中经常用到的函数，其中 是常数。

令如此与求单位阶跃函数同理，就可求得〔〕设，，如此由欧拉公式，有所以(2.13〕同理 (〕4.单位脉冲函数 δ(t)的拉氏变换单位脉冲函数是在持续时间期间幅值为的矩形波其幅值和作用时间的乘积等于1，即如图2.8所示图2.8 单位脉冲函数单位脉冲函数的数学表达式为其拉氏变换式为此处因为时，，故积分限变为2.15)单位速度函数，又称单位斜坡函数，其数学表达式为见图所示图2.9 单位速度函数单位速度函数的拉氏变换式为利用分部积分法令如此所以当时，,如此〔〕单位加速度函数的数学表达式为图2.10 单位加速度函数其拉氏变换式为〔〕2.5.3 拉氏变换的主要定理根据拉氏变换定义或查表能对一些标准的函数进展拉氏变换和反变换，但利用以下的定理，如此对一般的函数可以使运算简化拉氏变换也服从线性函数的齐次性和叠加性〔1〕齐次性 设，如此〔2.18〕式中——常数〔2〕叠加性设,，如此〔2.19〕两者结合起来，就有这说明拉氏变换是线性变换设如此式中——函数在 时刻的值，即初始值同样，可得的各阶导数的拉氏变换是〔〕式中,，…——原函数各阶导数在时刻的值如果函数与其各阶导数的初始值均为零〔称为零初始条件〕，如此各阶导数的拉氏变换为〔2.21〕假如 可以进展拉氏变换，如此除了在 的极点以外，〔〕式中，。

同样有一般地，有〔〕设 ，如此〔〕式中——积分 在 时刻的值当初始条件为零时，〔〕对多重积分是〔〕当初始条件为零时，如此〔〕设，且时，，如此〔〕函数为原函数沿时间轴延迟了，如图2.11所示图2.11 函数在控制理论中，经常遇到一类的函数，它的象函数只需把用代替即可，这相当于在复数坐标中，有一位移设，如此〔〕例如 的象函数，如此的象函数为它明确原函数在 时的数值〔2.30〕即原函数的初值等于 乘以象函数的终值设，并且存在，如此〔〕即原函数的终值等于乘以象函数的初值这一定理对于求瞬态响应的稳态值是很有用的设，，如此有〔〕即两个原函数的卷积分的拉氏变换等于它们象函数的乘积式〔2.32〕中，为卷积分的数学表示，定义为〔〕式中——比例系数例如,的象函数,如此的象函数为在某些情况下,在处有一个脉冲函数这时必须明确拉普拉斯积分的下限是还是，因为对于这两种下限，的拉氏变换是不同的为此，可采用如下符号予以区分：假如在处包含一个脉冲函数，如此因为在这种情况下显然，如果在处没有脉冲函数，如此有2.5.4 拉普拉斯反变换拉普拉斯反变换的公式为〔〕式中——表示拉普拉斯反变换的符号通常用局部分式展开法将复杂函数展开成有理分式函数之和，然后由拉氏变换表一一查出对应的反变换函数，即得所求的原函数。

在控制理论中，常遇到的象函数是的有理分式为了将写成局部分式，首先将的分母因式分解，如此有式中，，，…，是的根的负值，称为的极点，按照这些根的性质，可分为以下几种情况来研究2. 的极点为各不一样的实数时的拉氏反变换〔〕式中，是待定系数，它是处的留数，其求法如下〔2.38〕再根据拉氏变换的迭加定理，求原函数[例 2.1] 求的原函数解:首先将的分母因式分解，如此有 即得3.含有共轭复数极点时的拉氏反变换如果 有一对共轭复数极点,，其余极点均为各不一样的实数极点将展成式中,和可按下式求解即〔〕因为(或〕〕两边都应是复数，令等号两边的实部、虚局部别相等，得两个方程式，联立求解，即得，两个常数[例 2.2] ,试求其局部分式解:因为〔〕含有一对共轭复数极点,和一个极点，故可将式〔2.40〕因式分解成〔2.41〕以下求系数、和由式〔2.40〕和式〔2.41〕相等，有〔2.42〕 用乘以上式两边，并令，得到上式可进一步写成由上式两边实部和虚局部别相等，可得联立以上两式，可求得为了求出系数,用乘方程〔2.42〕两边，并令,将代入，得 将所求得的,, 值代入〔2.41〕，并整理后得的局部分式查拉氏变换表便得， 结果见式〔3.16〕。

[例 2.3] 求解: 将的分母因式分解，得利用方程两边实部、虚局部别相等得解得，所以这种形式再作适当变换：查拉氏变换表得4.中含有重极点的拉氏反变换 设有r个重根，如此将上式展开成局部分式〔2.43〕 式中，，，…，的求法与单实数极点情况下一样…，的求法如下：……如此〕[例 2.4] 设，试求的局部分式 解: 〔2.45〕含有2个重极点，可将式〔2.45〕的分母因式分解得〔2.46〕以下求系数、和将所求得的、、值代入式〔2.46〕，即得的局部分式查拉氏变换表可得[例 2.5] 求的拉氏反变换 解: 将展开为局部分式上式中各项系数为于是查拉氏变换表，得(1) 概述MATLAB是美国Math Works公司的软件产品，是一个高级的数值分析、处理与计算的软件，其强大的矩阵运算能力和完美的图形可视化功能，使得它成为国际控制界应用最广的首选计算机工具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=[1 11 39 52 26]den= [1 10 35 50 24]命令[r,p,k]=residue(num,den)于是得到如下结果[r,p,k]=residue(num,den) r=2.5000 p=k= 1如此得如果F(s)中含重极点，如此局部分式展开式将包括如下诸项式中，p(j)为一个q重极点[例2.7] 试将如下函数展开成局部分式解：对于该函数有num=[0 1 4 6]den =[1 3 3 1]命令[r,p,k]=residue(num,den)将得到如下结果：[r,p,k]=residue(num,den)r=p=k=[ ]所以可得注意，本例的余项k为零。

2.5.5 应用拉氏变换解线性微分方程应用拉氏变换解线性微分方程时，采用如下步骤：(1) 对线性微分方程中每一项进展拉氏变换，使微分方程变为 的代数方程；(2) 解代数方程，得到有关变量的拉氏变换表达式；(3) 用拉氏反变换得到微分方程的时域解整个求解过程如图2.12所示2.8 设系统微分方程为假如，初始条件分别为 、，试求解: 对微分方程左边进展拉氏变换利用迭加定理将上式逐项相加，即得方程左边的拉氏变换对方程右边进展拉氏变换得写成一般形式应该强调指出是微分方程的特征方程，也是该系统的特征方程利用局部分式将展开为求待定系数、、、、:代入原式得查拉氏变换表得当初始条件为零时，得 / 。