常微分方程差分方程解法归纳
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
‘
P(x)dx
C (x) =Q(x)e ,
,再对其两边积分得
fP(x) dx
C(x)二.Q(x)e dx C ,于是将其回代入
常微分方程解法归纳
1. 一阶微分方程部分
①可分离变量方程(分离变量法) 如果一阶微分方程 d
^ = f (x, y)中的二元函数 f (x, y)可表示为f (x, y)二
g(x)h(y) dx 的形式,我们称 3 =g(x)h(y)为可分离变量的方程。
dx 对于这类方程的求解我们首先将其分离变量为 -dy g(x)dx 的形式,再对此式两边积 h(y)
分得到 型 g(x)dx C 从而解出 3二g(x)h(y)的解,其中C 为任意常数。
' h(y) ' dx 具体例子可参考书本 P10 — P11的例题。
②一阶线性齐次、非齐次方程(常数变易法) 如果一阶微分方程史=f (x, y)中的二元函数f (x, y)可表示为 dx f(x, y) =Q(x) - P(x)y 的形式,我们称由此形成的微分方程 dy P(x)y =Q(x)为一阶线 dx
性微分方程,特别地,当 Q(x) =0时我们称其为一阶线性齐次微分方程,否则为一阶线性 非齐次微分方程。
对于这类方程的解法,我们首先考虑一阶线性齐次微分方程
裂P(x)厂0,这是可 —P(x)dx
分离变量的方程,两边积分即可得到 y 二Ce • ,其中 C 为任意常数。
这也是一阶线性 非齐次微分方程的特殊情况,两者的解存在着对应关系,设 C(x)来替换C ,于是一阶线性 非齐次微分方程存在着形如 y=C(x)e - …P(x)dx …P(x)dx
得至U C (x)e —P(x)C(x)e
-P(x)dx dy 的解。
将其代入 P(x)y 二Q(x)我们就可 dx
…P(x)dx
P(x)C(x)e • 二Q(x)这其实也就是 —'P(x)dx y = C(x)e 即得一阶线性
微分方程鱼,P(x)y =Q(x)的通解 dx
-P(x)dx
y =e .
Q(x)e
P(x)dx
dx + C
I 。
具体例子可参照书本P16 —P17的例题。
③一阶齐次型微分方程(变量代换)
如果一阶微分方程矽二f (x, y)中的二元函数f (x, y)满足对于一切非零实数t都有等
dx
式f (tx,ty) = f (x, y)成立,我们称一阶微分方程为一阶齐次型微分方程。
dx
对于此类微分方程的解法,我们一般利用变量代换的方法将其化为一阶可分离变量的方
程然后再相应求解。
事实上,如果我们令t二丄于是f(x,y)二f(1,#) = (y)。
于是一阶齐次型微分方程
x x x
dy = f (x, y)可表示为包二(丫)然后令u二y将其化为一阶可分离变量微分方程。
具体dx dx x
x
过程如下:令,则y=xu,鱼=厂x du,代入方程鱼=(')可得
x dx dx dx x
u x d^ = (u)也就是虫=(u) 7,它的通解是易求得的,求出它的通解之后将u =丄dxdx x x
回代就可得到一阶齐次型微分方程鱼二f (x, y)的通解。
dx
当然,有时候我们令t二丄于是f(x,y) = f (-,1H (y) °于是一阶齐次型微分方程
y x x
dy dy “ x、,亠,dx 1 “ 「人x dx dv
f (x, y)可表示为(一)也就是此时令v ,贝V v y -,
dx dx y dy ‘..气冬) y dy dy
y
dx 1 dv 1
代入方程可得v y 然后再依次求解。
有时候后者的代换方法会更dy屮(仝) dy屮(v)
y
简洁,当然两者的解法本质上是没有区别的,具体求解时可以灵活地运用。
具体例子可参看书本P20 —P22的例题。
④伯努利方程(变量代换)
当厶二ab i - a i b = 0时,〉「有唯一解,可以化上面的方程为齐次方程
,求解此方程,并将上=x -〉J ; = y - 1代回就得到齐次方程
的解。
当厶二ad -aR =0时要分两种情况讨论。
如果一阶微分方程 3 = f (x, y)中的二元函数f (x, y)满足
dx
f (x, y) = Q(x)y n
-P(x)y,(n 严0,1),我们就称由此形成的微分
dy
P(x)y = Q(x) y n
,(n = 0,1)为伯努利方程。
dx
对于此类方程的求解,我们可以通过变量替换将其转化为一阶线性微分方程求解。
我们
可以在方程少 P(x)y =Q(x)y n
,( n = 0,1)两边同除以y n
,可以将方程变形为
dx
y 』巴 P(x) y 1
』=Q(x)即-—】P(x)y 1jn
= Q(x)。
我们令 z = y 1J
,于是方 dx 1 一 n dx
程即dz
• (1 - n)P(x)z =(1 - n)Q(x)利用一阶线性微分方程 dy
- P(x)y =Q(x)的通解
dx
dx
-P(x)dx
y 二 e
Q(x)e
P(x)dx
dz
dx+C [可得 ——+(1 - n)P(x)z = (1 - n)Q(x)的通解,
再将 丿 dx
z = y 1
』回代就得到了伯努利方程 3 P(x)y 二Q(x)y n
,(n = 0,1)的通解。
dx
具体例子可参照书本 P22 — P23的例题。
⑤变量代换方法的应用----其他类型的齐次微分方程
ax by
的方程也是齐次方程。
对于这种类型的方程通过简单的代换就
可以化为一阶齐次型微分方程来进行求解。
我们讨论更一般的情形,对于形如
ax by c
的齐次方程,我们令x= : j 二• 1 ,其中:-/■为待定常数,
a b a -> } b- c
&占 +b d +时 +
t hP +G /
,可以选取适当的使得」
ax by c
情况一:若,则—=—=k。
原方程可以化为巴=f k(ai X + biy) + c。
令
a i
b i dx J + by + G /
^a i x b i y,则y = ^(z -印刈得到变量可分离的方程
b i
按照相应的解法即可求解。
情况二:若d=o,则3“与b中至少有一个为0.当b = 0时,原方程为d^= f ax+ c
dx gx + G / 是可变量分离的方程,按照相应的解法即可求解。
当b = 0时,可以令
z = ax+by,型 =「空—a [ 原方程就变为了" C这是可变量分
dx b Idx 丿 b J dx 丿g x + G ,
离的方程,按照相应的解法即可求解。
具体例子可参看书本P24 —P25的例题。
2.可降阶的高阶微分方程部分(主要讨论二阶微分方程)
①形如y(n)= f(x)的微分方程
对于形如y(n)= f (x)的微分方程,我们可以连续对等式两边积分n次便可以求得其含
有n个任意常数的通解为y =•…f (x)dx 少◎... C n。
(n -1)! (n-2)!
n个积分符号
具体例子可参看书本P28例题。
②形如y = f(x,y)的微分方程
般二阶微分方程可以表示为y = g(x, y, y ),当因变量y不显含时形成了如
f (x, y )的不显含因变量y的二阶微分方程。
我们可以通过变量代换来进行降阶。
我
们令,于是方程可化为-^ = f (x, p),这是一个以p为未知函数,以x为
dx dx
自变量的一阶微分方程,我们可以容易求得。
设其通解为p= (x,G),则史二(x,C i)。
dx
两边积分就得到原方程的通解为y二(x,C1)dx C2。
其中C1,C2为任意常数。
具体可参看书本P28 —P30例题(注意例4! !)
③ 形如y = f (y, y )的微分方程
与不显含因变量 y 的二阶微分方程的定义类似,我们把形如 y"= f ( y, y )的微分方
程
称为不含自变量 X 的二阶微分方程。
我们仍然通过变量代换来求解此类方程。
我们令
p = y, y -上=—P
y = p —P ,于是方程可化为 p —P = f (y, p),这是一个关于 p,
y
dx dy
dx dy
dy
dx ,两边积分就得到原方程的通解为
'■ (y,C i )
常数。
具体例子可参看书本 P32 — P34例题。
注:在可降阶的微分方程求解问题中,在消去所设的变元如 p 时我们 一定要注意是否会丢失p=0的解。
3.线性微分方程
在介绍线性微分方程的解法之前有必要先介绍线性微分方程解的结构与性质。
我们直接
介绍n 阶线性微分方程y
(n)
- a i (x)y (nJ) - ... - a n 4(x)/ a n (x)y 二f (x)的解的结构与性
质。
对于区间[a,b ]上的n 个函数y^x)、y 2(x)、y 3(x)、……、y n (x), 若存在n 个不全为0
n
的常数k 1> k 2、k 3、 .. 、k n 使得在[a,b ]上有工K % (x) = 0 ,我们就称这n 个函数在区间
i =1
[a,b ]上是线性相关的,否则就是线性无关的。
的一阶微分方程,我们可以容易求得。
设其通解为
p= (y,C i ),则由 p = y = —丫 可得
dx
二 X C 2
其中C i ,C 2为任意
此外对于n阶线性微分方程y(n)' a1(x)y(n4) ' ... ' a n4(x)/ a n(x)y = f (x)的系数
a1(x)> a2(x)、...a n4(x)、a n(x)都为常数是我们称该方程为n阶线性常系数方程,否则为n 阶线性变系数方程。
进一步细分,对于自由项 f (x),若f (x)二0就称原方程为
n阶线性齐次方程,否则为n阶线性非齐次方程。
若函数yjx)、y2(x)、y3(x)、……、y n(x)是n
阶线性齐次方程的n个线性无关的特解,
则y = G y i(x) • C2y2(x) •…• C n y n(x)为
n阶线性齐次方程的通解。
若函数y i(x)、y2(x)、y3(x)、……、y n (x)是n
阶线性非齐次方程的n个线性无关的特
解,此外函数y p (x)是
n阶线性非齐次方程的1个线性无关的特解,则^y p(x) C1y1(x)
C2y2(x) ... - C n y n(x)为n阶线性齐次方程的通解。
①二阶常系数齐次线性微分方程
我们把形如y、py 'qy=0的微分方程称为二阶常系数齐次线性微分方程.其中p、q均为常数.我们知道如果 y i、y2是二阶常系数齐次线性微分方程的两个线性无关解.那么y£i y i C2y2就是它的通解.现在先尝试能否适当选取 r.使yR 满足二阶常系数齐次线性微分方程.为此将丫壬伙代入方程y'py'qy=O得(r 2卩r q)e rx =0 .由此可见.只要r满足代数方程r2pr q =0 .函数y=e rx就是微分方程的解,
接下来介绍一般的解法,我们把方程r2pr ^0叫做微分方程y'py「qy=0的特征方程,
特征方程的两个根 r i、r2可用公式r12=―P——一、卩-- 4^ 求出,
2
特征方程的根与通解的关系:
(1)特征方程有两个不相等的实根门、r2时•函数=e rix、y2二e r2x是方程的两个线性
v e r1x
无关的解,这是因为•函数y^e rix、y2二e r2x是方程的解.又Z =吕二* 不是常数
V2 e r2x
因此方程的通解为 y =C i e
riX
C z e
护.
⑵特征方程有两个相等的实根 rim 时.函数
y 1 =e rix
、y 2 =xe rix
是二阶常系数齐次线
性微分方程的两个线性无关的解.这是因为.y^e rix 是方程的解.又
(xe riX
) p(xe rix
) q(xe rix
)二(2口 x 『)e rix
p(i xr i )e rix
qxe
rix
=e rix (2r i p) xe rix (r i 2 pr i q^ 0 ,
所以y 2 =xe rix 也是方程的解.且里二嘍=x 不是常数.因此方程的通解为
r i
x
y =C i e"x
C 2xe"x
= (C i C 2X)e
rix
.
(3)特征方程有一对共轭复根 r i, 2=二匸 时.函数y=e (e)x
、y=e (:+)x
是微分方程的两个
线性无关的复数形式的解 .函数y=e :x cos x 、y 壬:x
sin x 是微分方程的两个线性无关的实数 形式的解.函数y i 壬(J )x
和y 2=e (*:)x
都是方程的解.而由欧拉公式.得
y i =e (: i )x =e :x (cos F x i sin F x).
故e x cos -x> y^e
x
sin :x 也是方程解,可以验证.y i=e-x cos :x 、y^e
x
sin :x 是方程的线性
无关解.因此方程的通解为 y=e x
(C i cos :x C 2Sin x ).
求二阶常系数齐次线性微分方程 y \py 「qy =0的通解的步骤为:
第一步写出微分方程的特征方 r 2
pr ^0
第二步求出特征方程的两个根 r i 、「2,
第三步
根据特征方程的两个根的不同情况
.写出微分方程的通解•
②二阶线性常系数非齐次方程
我们把形如y 'py qy ^f a )的微分方程称为二阶常系数非齐次线性微分方程 .其中p 、q 是
(cos
常数,二阶常系数非齐次线性微分方程的通解是对应的齐次方程的通解y=Y(x)与非齐次方程本身的一个特解y=y*(x)之和:yM(x) - y*(x).
当f(x)为两种特殊形式时.方程的特解的求法:
一、f(x)二P m(x)e 'X 型
当f(x)护m(x)e x时.可以猜想.方程的特解也应具有这种形式.因此.设特解形式为y* q(x)e x.将其代入方程.得等式
Q (x) (2 ■ p)Q (x) (2 p,q)Q(x)=P m(x).
(1)如果■不是特征方程r2卩r・q=O的根.则,2・p「q^0,要使上式成立.Q(x)应设为m
次多项式:
Q m(X)=b o X m b i X m十亠亠b m」x b m .
通过比较等式两边同次项系数.可确定b o .b i . .b m .并得所求特解
y* =Q m(x)e x.
⑵如果■是特征方程r2 pr ^0的单根.则 Zp'pn.但2'pa要使等式
Q (x) (2,p)Q (x) (2p・:;q)Q(x)=P m(x).
成立Q(x)应设为m 1次多项式:
Q(x)=xQ m(x).
Q m(x) =b o x m b i x m"亠亠b m4X b m .
通过比较等式两边同次项系数.可确定b o .b i.b m .并得所求特解
y* =xQ m(x)e x,
(3)如果九是特征方程r2巾r+q=0的二重根.则丸2十p?+q=O. 2^+p=0 .要使等式
Q (x) (2,p)Q (x) (2p・:;q)Q(x)二P m(x).
成立Q(x)应设为m 2次多项式:
Q m(x)龙0X m b i x m° 亠亠b m」x b m .
通过比较等式两边同次项系数.可确定bo .bi . • • • .bm .并得所求特解
y* =x2Qm(x)e x,
综上所述.我们有如下结论:如果f(x)二Pm(X)e'x.则二阶常系数非齐次线性微分方程y py qy』(x)有形如
y* =x k Q m(x)e x
的特解.其中Q m(x)是与P m(x)同次的多项式.而k按,不是特征方程的根、是特征方程的
单根或是特征方程的的重根依次取为0、1或2 .
二、f(x) =e 'x [P i(x)cos .x P n(x)sin,x]型
方程 y py qy =e x[P i(x)cos,x P n(x)sin x]的特解形式应用欧拉公式可得
e x[P i(x)cos x P n(x)sin ・x]
a /恥+小护e^x_e%x
=e x[P(x卢 2 P n(x)匕e ]
P n(x)
=2[P(x)-iP n(x)]e(」)x2【P(x) iP n(x)]e(z)x
=P(x)e('十')x • P(x)e('』')x
1 — 1
其中P(x)=2(P-PJ) . P(x^1(F l P n i) •而m=max{l.n}.
设方程 y py qy =P(x)e r i )x的特解为 y1* =x k Q m(x)e(' ' )x
则%* =x k Q m(x)e(' $ >必是方程目 p/ q^P(x)e r 4 0的特解
其中k按• _r ■不是特征方程的根或是特征方程的根依次取0或1 .
于是方程 y py、qy 弋x[P i(x)co^ x - P n(x)sin「x]的特解为
y* =x k Q m(x)e(' J )x• x k Q m(x)e('十)x
= x k e 'X[Q m(x)(cos x isin x) Q m(x)(cos x-i sin x)
=x k e x[R⑴m(x)cos ,x R(::2)m(x)sin ,x],
综上所述.我们有如下结论:
如果f(x)=e x[P i(x)cos • .x P n(x)sin • .x].则二阶常系数非齐次线性微分方程 y py qy =f(x) 的特解可设为 y*=x k e x[R⑴m(x)cos x R2)m(x)sin x]
其中R⑴m(x)、R(2)m(x)是m次多项式.m=max{ l . n}.而k按• h .(或,-i .)不是特征方程的根或是特征方程的单根依次取0或1 .
②高阶线性常系数微分方程
对于n阶线性微分方程y(n)- a nJ y(nJ) - ... - a i y a0y二f (x)的解。
我们首先讨论n 阶线性常系数微分方程y(n)■ a nd y(nJ)• qy ' a0y =0的解。
它的特征方程为
n _1
■n4」=0,与二阶的情况类似,故可按解得的情况按下表写出微分方程所对应的解。
i =0
根据上表写出函数n个线性无关的特解yjx)、y2(x)、y3(x)、……、y n(x),则
y二Gy i(x) • C2y2(x) •… C n y n(x)为n阶线性常系数齐次方程的通解。
其次对于y(n)' an4y(n4) '... ai/ a°y二f (x),我们主要是讨论其特解的形式若f (x) =P m(x)e,,其中P m(x)是m次待定的多项式,我们可设特解为
y p(x)二x k Q m(x)e x,其中Q m(x)是m次待定的多项式,k为〉作为特征方程根的重数。
即若〉不是特
征方程根,则k=0 ;若:•是L重特征方程根,则k=L
若f(x) =e:X[P m(x)cos -x Q h(x)sin -x],其中P m(x)、Q h(x)是m 次、h 次待定的多项式,我们可设特解为y p(x)二乂怙订& cos x T n sin -x],其中S n(x)、T n(x)是n次待定的多项式,n=max{m,h} , k为复数根作为特征方程根的重数。
单纯的课本内容,并不能满足学生的需要,通过补充,达到内容的完善
教育之通病是教用脑的人不用手,不教用手的人用脑,所以一无所能。
教育革命的对策是手脑联盟,结果是手与脑的力量都可以大到不可思议。