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

  • 格式:docx
  • 大小:38.24 KB
  • 文档页数:4

下载文档原格式

  / 4
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

拉普拉斯变换是解常系数线性微分方程中经常采用的一种较简便的方法.其基本思想是,先通过拉普拉斯变换将已知方程化成代数方程,求出代数方程的解,再通过逆拉普拉斯变换,得到所求数值问题的解•

一拉普拉斯变换的概念

定义设函数f(t)的定义域为[0,+ R若广义积分jT^f(t)e-Pt dt对于p在某一范

围内的值收敛,则此积分就确定了一个参数为p的函数记作F(p),即F(p)= o F f(t)e-pt dt函数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)=e at(t > 0是常数)的拉氏变换.

解根据定义,有L[e at]= cTe at e-pt dt二(re-(p-a)t dt

这个积分在p > a时收敛所以有

L[e a1]= 0+TO e-(p-a)t dt=1/(p-a) (p> a) (1)

例2求一次函数f(t)=at(t 滾是常数)的拉氏变换.

解L[at]= 0T ate pt dt二-a/p oFtd(e-pt)

=-[at/p e-pt]o+cx+a/p o广-pt dt

根据罗必达法则,有

lim to+ «(-at/p e-pt)=-lim to+ «at/pe pt=-lim to+ «a/p2 e pt

上述极限当p >0时收敛于0,所以有lim to+ «(-at/pe-pt)=0

因此L[at]=a/p o+T e-pt dt

二-[a/p2e-pt]o+cx=a/p2(p> 0)

例3求正弦函数f(t)=sin 3 t(t的拉氏变换.

解L[sin 3 t]=o+°!Sin w-tedt

=[-1/(p2+w2) e-pt(psin 3 t+ 3 coo+°3 t]

=3 /(p+w2) (p > 0) (3)

用同样的方法可求得

L[cos 3 t]=p/(p+ 32) (p > 0) (4)

二拉普拉斯变换的基本性质

三拉普拉斯变换的逆变换

四拉普拉斯变换的应用

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

拉普拉斯变换方法是解线性微分方程的一种简便方法,利用拉普拉斯变换法可以把微分方程变换成为代数方程,在利用现成的拉普拉斯变换表(参见附录一的附表1),即可方便地查得相应的微分方程解。这样就使方程求解问题大为简化。

拉普拉斯变换法的另一个优点是在求解微分方程时,可同时获得的瞬态分量和稳态分量两部分。

有关拉普拉斯变换(简称拉氏变换)的公式见附录一。

应用拉氏变换法得到的解是线性微分方程的全解。用古典方法求解微分方程全解时需要利用初始条件来确定积分常数的值,这一过程比较麻烦。而应用拉氏变换就可省去这一步。因为初始条件已自动地包含在微分方程的拉氏变换式之中了。而且,如果所有初始条件都为零,那么求

取微分方程的拉氏变换式就更为方便,只要简单地用复变量s 来代替微分方程中的 —,s 2代替

dt

…就可得到。

dt

应用拉氏变换法解微分方程的步骤如下: (1)

对线性微分方程中每一项

进行拉氏变换,使微分方程变为复变量 s 的代数方程(称为 变换

方程)

(2) 求解变换方程,得出系统输出变量的象函数表达式。 (3)

将输出的象函数表达式展开成部分分式(部分分式展开法参见附录二)

(4) 对部分分式进行拉氏反变换(可查拉氏变换表),即得微分方程的全解。 举例说明 【例2-7】

设RC 网络如图2-24所示,在开关 K 闭合之前,电容 C 上有初始电压

U c (O)。试求将开关瞬时闭合后,电容的端电压

U c (网络输出)。

U r Ri U c

1 U c c C

idt

消去中间变量i ,得网络微分方程为

可见等式右边由两部分组成,一部分由输入所决定,另一部分由初始值决定。

将输出的象函数U c (s)展成部分分式:

RC U c (s)

1 — U O

s

RC RCs 1U O

RC RCs 1

U c (0)

1

U c (s) -U O

s

r U o

厂 U c (°)

(2-45 )

解 开关K 瞬时闭合,相当于网络有阶跃电压 为

U c (t) U O • 1(t)输入。故网络微分方程 RC 詈…) 对上式进行拉氏变换,得变换方程 (2-44

) RCsU c (s) RCu c (0) U c (s) U r (s)

将输入阶跃电压的拉氏变换式 U r (S )

U 0

—代入上式,并整理得电容端电压的拉氏变换式

s

U c (s)

U O

s(RCs 1)

RC (RCs 1)

U c (0)

RC

等式两边进行拉氏反变换,得

此式表示了 RC 网络在开关闭合后输出电压 u c (t)的变化过程。

比较方程(2-45)和(2-46)可见,方程右端第一项 取决于外加的输入作用U o • 1(t),表示了网络输出响应 u c (t)的稳态分量,也称强迫解;第二项表示 U c (t)的瞬 态分量,该分量随时间变化的规律取决于系统结构参量

R 、C 所决定的特征方程式(即 RCs 1

0 )的根

1 。显然,由于其特征根为负实数,则瞬态分量将随

RC

着时间的增长而衰减至零。第三项为与初始值有关的瞬态 分量,其随时间变化的规律同样取决于特征根,当初始值

U c (0)

0时,则第三项为零,于是就有

RC 网络的阶跃响应U c (t)及其各组成部分的曲线如图 2-25所示。

U c (t) U o u °e RC U c (O)e RC

(2-46 )

U c (t) U o

U o e RC

(2-47 )