用拉普拉斯变换方法解微分方程

拉普拉斯变换是解常系数线性微分方程中经常采用的一种较简便的方法.其基本思想是,先通过拉普拉斯变换将已知方程化成代数方程,求出代数方程的解,再通过逆拉普拉斯变换,得到所求数值问题的解.一拉普拉斯变换的概念定义设函数f(t)的定义域为[0,+∞),若广义积分∫0+∞f(t)e-ptdt对于p在某一范围内的值收敛,则此积分就确定了一个参数为p的函数,记作F(p),即F(p)=∫0+∞f(t)e-ptdt函数F(p)称为f(t)的拉普拉斯变换(或称为f(t)的象函数),表示为F(p)=L[f(t)].若F(p)是f(t)的拉氏变换,则称f(t)为F(p)的拉氏逆变换(或F(p)的象原函数),记作L-1[F(p)].例1求指数函数f(t)=eat(t≥0,a是常数)的拉氏变换.解根据定义,有L[eat]=∫0+∞eate-ptdt=∫0+∞e-(p-a)tdt这个积分在p＞a时收敛,所以有L[eat]=∫0+∞e-(p-a)tdt=1/(p-a)(p＞a)(1)例2求一次函数f(t)=at(t≥0,a是常数)的拉氏变换.解L[at]=∫0+∞ate-ptdt=-a/p∫0+∞td(e-pt)=-[at/pe-pt]0+∞+a/p∫0+∞e-ptdt根据罗必达法则,有limt0+∞(-at/pe-pt)=-limt0+∞at/pept=-limt0+∞a/p2ept上述极限当p＞0时收敛于0,所以有limt0+∞(-at/pe-pt)=0因此L[at]=a/p∫0+∞e-ptdt=-[a/p2e-pt]0+∞=a/p2(p＞0)(2)例3求正弦函数f(t)=sinωt(t≥0)的拉氏变换.解L[sinωt]=∫0+∞sinωte-ptdt=[-1/(p2+ω2)e-pt(psinωt+ωcosωt]0+∞=ω/(p2+ω2)(p＞0)(3)用同样的方法可求得L[cosωt]=p/(p2+ω2)(p＞0)(4)二拉普拉斯变换的基本性质三拉普拉斯变换的逆变换四拉普拉斯变换的应用2–5用拉普拉斯变换方法解微分方程拉普拉斯变换方法是解线性微分方程的一种简便方法，利用拉普拉斯变换法可以把微分方程变换成为代数方程，在利用现成的拉普拉斯变换表（参见附录一的附表1），即可方便地查得相应的微分方程解。这样就使方程求解问题大为简化。拉普拉斯变换法的另一个优点是在求解微分方程时，可同时获得的瞬态分量和稳态分量两部分。有关拉普拉斯变换（简称拉氏变换）的公式见附录一。应用拉氏变换法得到的解是线性微分方程的全解。用古典方法求解微分方程全解时需要利用初始条件来确定积分常数的值，这一过程比较麻烦。而应用拉氏变换就可省去这一步。因为初始条件已自动地包含在微分方程的拉氏变换式之中了。而且，如果所有初始条件都为零，那么求取微分方程的拉氏变换式就更为方便，只要简单地用复变量s来代替微分方程中的dtd，2s代替22dtd，…就可得到。应用拉氏变换法解微分方程的步骤如下：（1）对线性微分方程中每一项进行拉氏变换，使微分方程变为复变量s的代数方程（称为变换方程）（2）求解变换方程，得出系统输出变量的象函数表达式。（3）将输出的象函数表达式展开成部分分式（部分分式展开法参见附录二）。（4）对部分分式进行拉氏反变换（可查拉氏变换表），即得微分方程的全解。举例说明【例2-7】设RC网络如图2-24所示，在开关K闭合之前，电容C上有初始电压)0(cu。试求将开关瞬时闭合后，电容的端电压cu（网络输出）。解开关K瞬时闭合，相当于网络有阶跃电压0)(utuc·)(1t输入。故网络微分方程为idtCuuRiuccr1消去中间变量i，得网络微分方程为)(tuudtduRCrcc（2-44）对上式进行拉氏变换，得变换方程)()()0()(sUsURCusRCsUrccc将输入阶跃电压的拉氏变换式susUr0)(代入上式，并整理得电容端电压的拉氏变换式)0()1()1()(0ccuRCsRCRCssusU可见等式右边由两部分组成，一部分由输入所决定，另一部分由初始值决定。将输出的象函数)(sUc展成部分分式：)0(111)(00ccuRCsRCuRCsRCussU或)0(11111)(00ccuRCsuRCsussU（2-45）等式两边进行拉氏反变换，得tRCctRCceueuutu1100)0()(（2-46）此式表示了RC网络在开关闭合后输出电压)(tuc的变化过程。比较方程（2-45）和（2-46）可见，方程右端第一项取决于外加的输入作用0u·1)(t，表示了网络输出响应)(tuc的稳态分量，也称强迫解；第二项表示)(tuc的瞬态分量，该分量随时间变化的规律取决于系统结构参量R、C所决定的特征方程式（即01RCs）的根RC1。显然，由于其特征根为负实数，则瞬态分量将随着时间的增长而衰减至零。第三项为与初始值有关的瞬态分量，其随时间变化的规律同样取决于特征根，当初始值0)0(cu时，则第三项为零，于是就有tRCceuutu100)(（2-47）RC网络的阶跃响应)(tuc及其各组成部分的曲线如图2-25所示。