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

拉普拉斯变换是解常系数线性微分方程中经常采用的一种较简便的方法.其基本思想是,先通过拉普拉斯变换将已知方程化成代数方程,求出代数方程的解,再通过逆拉普拉斯变换,得到所求数值问题的解. 一 拉普拉斯变换的概念 定义 设函数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)=Lf(t). 若F(p)是f(t)的拉氏变换,则称f(t)为F(p)的拉氏逆变换(或F(p)的象原函数),记作L-1F(p). 例1 求指数函数f(t)=eat(t0,a是常数)的拉氏变换. 解根据定义,有Leat=0+eate-ptdt=0+e-(p-a)tdt 这个积分在pa时收敛,所以有 Leat=0+e-(p-a)tdt=1/(p-a) (pa) (1) 例2 求一次函数f(t)=at(t0,a是常数)的拉氏变换. 解 Lat=0+ate-ptdt=-a/p0+td(e-pt) =-at/p e-pt0+a/p0+e-ptdt 根据罗必达法则,有 limt0+(-at/p e-pt)=-limt0+at/pept=-limt0+a/p2 ept 上述极限当p0时收敛于0,所以有limt0+(-at/pe-pt)=0 因此Lat=a/p0+e-ptdt =-a/p2e-pt0+=a/p2(p0) (2) 例3 求正弦函数f(t)=sint(t0)的拉氏变换. 解 Lsint=0+sinte-ptdt=-1/(p2+2) e-pt(psint+cost0+=/(p2+2) (p0) (3) 用同样的方法可求得 Lcost=p/(p2+2) (p0)(4) 二 拉普拉斯变换的基本性质 三 拉普拉斯变换的逆变换 四 拉普拉斯变换的应用25 用拉普拉斯变换方法解微分方程拉普拉斯变换方法是解线性微分方程的一种简便方法，利用拉普拉斯变换法可以把微分方程变换成为代数方程，在利用现成的拉普拉斯变换表（参见附录一的附表1），即可方便地查得相应的微分方程解。这样就使方程求解问题大为简化。拉普拉斯变换法的另一个优点是在求解微分方程时，可同时获得的瞬态分量和稳态分量两部分。有关拉普拉斯变换（简称拉氏变换）的公式见附录一。应用拉氏变换法得到的解是线性微分方程的全解。用古典方法求解微分方程全解时需要利用初始条件来确定积分常数的值，这一过程比较麻烦。而应用拉氏变换就可省去这一步。因为初始条件已自动地包含在微分方程的拉氏变换式之中了。而且，如果所有初始条件都为零，那么求取微分方程的拉氏变换式就更为方便，只要简单地用复变量来代替微分方程中的，代替，就可得到。应用拉氏变换法解微分方程的步骤如下：（1）对线性微分方程中每一项进行拉氏变换，使微分方程变为复变量的代数方程（称为变换方程）（2）求解变换方程，得出系统输出变量的象函数表达式。（3）将输出的象函数表达式展开成部分分式（部分分式展开法参见附录二）。（4）对部分分式进行拉氏反变换（可查拉氏变换表），即得微分方程的全解。举例说明【例2-7】设网络如图2-24所示，在开关闭合之前，电容上有初始电压。试求将开关瞬时闭合后，电容的端电压（网络输出）。解开关瞬时闭合，相当于网络有阶跃电压输入。故网络微分方程为消去中间变量，得网络微分方程为（2-44）对上式进行拉氏变换，得变换方程将输入阶跃电压的拉氏变换式代入上式，并整理得电容端电压的拉氏变换式可见等式右边由两部分组成，一部分由输入所决定，另一部分由初始值决定。将输出的象函数展成部分分式：或 （2-45）等式两边进行拉氏反变换，得（2-46）此式表示了网络在开关闭合后输出电压的变化过程。比较方程（2-45）和（2-46）可见，方程右端第一项取决于外加的输入作用1，表示了网络输出响应的稳态分量，也称强迫解；第二项表示的瞬态分量，该分量随时间变化的规律取决于系统结构参量、所决定的特征方程式（即）的根。显然，由于其特征根为负实数，则瞬态分量将随着时间的增长