§3.3 奇解
曲线族的包络?(x,y,c)=0, ?
c
’(x,y,c)=0
'(,,)0, (,,)0.
c x y c x y c
Φ=?
?
Φ=?
c -判別曲线法
?(x,y,c )=0,?c ’(x,y,c )=0
例12
2
2
x
y p xp =-+
2
2
d d .d d 2y y x y x
x x ??
=-+ ???
d d y p x
=
d d 2,
d d p p p p x p x x x =---d 10,d p x
-=2
22
x y cx c
=++d 1(2)0d p p x x ??
--=
???
.p x c =+2
2
()()2x
y x c x x c =+-++
(续)例1
奇解
2
2
2dy dy x y x dx dx ??
=-+
???
2
22
x y cx c
=++1(2)0dp p x dx ??
--= ???
20p x -=2
4
x y =
10dp
dx
-=
(续)例1
奇解2
2
2dy dy x y x dx dx ??
=-+
???
2
2
2
x y cx c =++22,2
02.x y cx c x c ?=
++???=+?
2
4
x y =
求直线轴x cos a+y sin a-p=0 的包
络
奇解
求曲线轴(y-c)2-2/3(x-c)3=0 的包络
p -判別曲线法F (x,y,p )=0, F p ’(x,y,p )=0
0F
p
?=?'(,,)0,(,,)0p F x y p F x y p =???=??
求方程的奇解
2
2
d 10d y y x ??+-= ???
2210,
20.
p y p ?+-=?
=?
求方程
的奇解
2
d d 2d d y y y x x x ??=- ?
??
22,220.
y xp p x p ?=-?
-=?
克莱罗微分方程y=xp+f(p)
求解方程y=xp+1/p
210,
1.x c y cx c ?-=???
?=+??
2
4y x
=
求一曲线,使其上每一点的切线截割坐标轴而成的直角三角形中的面积等于2。
1x y a b
+=d d y b x a
=-2
d d 4
,d d y y y x x x ?
?-=- ??
d d 2.
d d y y y x x x
=±-2
1122y c x c c c x
=±-=-
习题2.5 2.ydy x xdy ydx 2=- 。 解: 2x ,得: ydy x xdy ydx =-2 c y x y d +-=221 即c y x y =+2 2 1 4. xy x y dx dy -= 解:两边同除以x ,得 x y x y dx dy - =1 令u x y = 则dx du x u dx dy += 即 dx du x u dx dy +=u u -=1 得到 ()2ln 2 1 1y c u -=, 即2 ln 21?? ? ??-=y c y x 另外0=y 也是方程的解。 6.()01=-+xdy ydx xy 解:0=+-xydx xdy ydx x d x y x d y y d x -=-2 得到c x y x d +-=??? ? ??2 21
即 c x y x =+2 2 1 另外0=y 也是方程的解。 8. 32 x y x y dx dy += 解:令 u x y = 则: 21u x u dx du x u dx dy +=+= 即2 1u x dx du x = 得到22x dx u du = 故c x u +-=-11 即 21 1x x c y += 另外0=y 也是方程的解。 10. 2 1?? ? ??+=dx dy dx dy x 解:令 p dx dy = 即p p x 2 1+= 而 p dx dy =故两边积分得到 c p p y +-=ln 2 12 因此原方程的解为p p x 21+=,c p p y +-=ln 212 。 12.x y xe dx dy e =?? ? ??+-1 解: y x xe dx dy +=+1
页脚内容1 第 一 章 一阶微分方程的解法的小结 ⑴、可分离变量的方程: ①、形如 )()(y g x f dx dy = 当0)(≠y g 时,得到 dx x f y g dy )()(=,两边积分即可得到结果; 当0)(0=ηg 时,则0)(η=x y 也是方程的解。 例1.1、xy dx dy = 解:当0≠y 时,有xdx y dy =,两边积分得到)(2ln 2为常数C C x y += 所以)(11212 C x e C C e C y ±==为非零常数且 0=y 显然是原方程的解; 综上所述,原方程的解为)(1212 为常数C e C y x = ②、形如0)()()()(=+dy y Q x P dx y N x M 当0)()(≠y N x P 时,可有dy y N y Q dx x P x M ) ()()()(=,两边积分可得结果; 当0)(0=y N 时,0y y =为原方程的解,当0(0=) x P 时,0x x =为原方程的解。 例1.2、0)1()1(22=-+-dy x y dx y x
页脚内容2 解:当0)1)(1(22≠--y x 时,有dx x x dy y y 1 122-=-两边积分得到 )0(ln 1ln 1ln 22≠=-+-C C y x ,所以有)0()1)(1(22≠=--C C y x ; 当0)1)(1(22=--y x 时,也是原方程的解; 综上所述,原方程的解为)()1)(1(22为常数C C y x =--。 ⑵可化为变量可分离方程的方程: ①、形如)(x y g dx dy = 解法:令x y u = ,则udx xdu dy +=,代入得到)(u g u dx du x =+为变量可分离方程,得到)(0),,(为常数C C x u f =再把u 代入得到)(0),,(为常数C C x x y f =。 ②、形如)0(),(≠+=ab by ax G dx dy 解法:令by ax u +=,则b du adx dy +=,代入得到)(1u G b a dx du b =+为变量可分离方程,得到)(0),,(为常数C C x u f =再把u 代入得到)(0),,(为常数C C x by ax f =+。 ③、形如 )(222111c y b x a c y b x a f dx dy ++++= 解法:01、02211 =b a b a ,转化为)(by ax G dx dy +=,下同①; 02、0221 1 ≠b a b a ,???=++=++00222111c y b x a c y b x a 的解为),(00y x ,令???-=-=00y y v x x u
摘要.................................................... 错误!未定义书签。 1.何谓奇解.............................................. 错误!未定义书签。 2.奇解的产生............................................ 错误!未定义书签。 3.包络跟奇解的关系...................................... 错误!未定义书签。 4.理论上证明C-判别曲线与P-判别曲线方法................. 错误!未定义书签。 克莱罗微分方程 ..................................... 错误!未定义书签。 5.奇解的基本性质........................................ 错误!未定义书签。 定理1 ............................................. 错误!未定义书签。 定理2 ............................................. 错误!未定义书签。 定理3 ............................................. 错误!未定义书签。 6.小结.................................................. 错误!未定义书签。参考文献:.............................................. 错误!未定义书签。
习题2.2 求下列方程的解。 1.dx dy =x y sin + 解: y=e ?dx (?x sin e ?-dx c dx +) =e x [- 2 1e x -(x x cos sin +)+c] =c e x -21 (x x cos sin +)是原方程的解。 2.dt dx +3x=e t 2 解:原方程可化为: dt dx =-3x+e t 2 所以:x=e ?-dt 3 (?e t 2 e -? -dt 3c dt +) =e t 3- (5 1e t 5+c) =c e t 3-+5 1e t 2 是原方程的解。 3.dt ds =-s t cos +21t 2sin 解:s=e ?-tdt cos (t 2sin 2 1?e dt dt ?3c + ) =e t sin -(?+c dt te t t sin cos sin ) = e t sin -(c e te t t +-sin sin sin ) =1sin sin -+-t ce t 是原方程的解。 4. dx dy n x x e y n x =- , n 为常数. 解:原方程可化为:dx dy n x x e y n x += )(c dx e x e e y dx x n n x dx x n +??=?- )(c e x x n += 是原方程的解.
5. dx dy +1212--y x x =0 解:原方程可化为:dx dy =-1212+-y x x ?=-dx x x e y 1 2(c dx e dx x x +?-221) )21(ln 2+=x e )(1 ln 2?+--c dx e x x =)1(1 2 x ce x + 是原方程的解. 6. dx dy 234xy x x += 解:dx dy 234xy x x += =23y x +x y 令 x y u = 则 ux y = dx dy =u dx du x + 因此:dx du x u +=2u x 21u dx du = dx du u =2 c x u +=33 1 c x x u +=-33 (*) 将x y u =带入 (*)中 得:3433cx x y =-是原方程的解.
总结一阶微分方程奇解的求法 摘要:利用有关奇解的存在定理,总结出求一阶微分方程奇解的几种方法,并通过一些具体的例题说明这几种方法的应用 Using relevant theorems to develop several methods of finding singular solution of ordinary differential equation. In addition, illustrate the application of these methods through the concrete examples. 关键词:常微分方程 奇解 c-判别式 p-判别式 方法一:利用c-判别式求奇解 设一阶微分方程0, ,=?? ? ?? dx dy y x F ① 可求出方程①的通解为()0,,=c y x φ ② 如果()()???==0 ,,0,,' c y x c y x c φφ ③ 是微分方程①的解,且对③式满足:()()02 '2 '≠+y x φφ ④ 则③是微分方程①的奇解,且是通解②的包络。 例1:方程() 2 2 2 x x y dy dx dy dx + -= 的奇解 解:首先,本具题意求出该微分方程的通解为2 2 2 c cx y x ++= 与4 2 x y = 其中c 为任意常数 当时2 2 2 c cx y x ++= , ()y c cx x c y x -++= 2 2 2 ,,φ 其相应的c -判别式为 ? ??=+=-++02022x 2 c x y c cx 易得到: ? ??=-=2 2c y c x
代入原微分方程,可知? ??=-=2 2c y c x 不是原微分方程的解; 当4 2 x y = 时,易求出2 ,1''x y x ==φφ,则有()()02 '2 '≠+y x φφ 故4 2 x y = 为原微分方程的奇解 例2:试求微分方程() () y y dy dx 9 42 2 1= -的奇解 解:首先,根据题意求出微分方程的通解为:()()0322=---y y c x 其中c 为任意常数 再由相应的c-判别式: ()()()? ??=--=---020 322c x y y c x 易求出:? ??==0y c x 或 ???==3y c x 当???==0y c x 时,代入原微分方程成立; 所以? ??==0y c x 为原微分方程的解 且有()02'=--=c x x φ;()()93232 '-=---=y y y y φ 满足(Φ‘ x )2 +(Φ‘ y )2≠0 易验证???==3y c x 不是原微分方程的解 故x=c, y=0 是元微分方程的奇解。 方法二:利用p-判别法求奇解 在微分方程①中,设y ′=p,则此方程的p-判别式为: ()()?????==0,,0 ,,' p y x F p y x F p ⑤ 消去p 之后得到的函数y=?(x)是微分方程①身为解,
习题2.2 求下列方程的解 1. dx dy =x y sin + 解: y=e ?dx (?x sin e ?-dx c dx +) =e x [- 21 e x -(x x cos sin +)+c] =c e x -2 1 (x x cos sin +)是原方程的解。 2. dt dx +3x=e t 2 解:原方程可化为: dt dx =-3x+e t 2 所以:x=e ? -dt 3 (?e t 2 e -?-dt 3c dt +) =e t 3- (5 1 e t 5+c) =c e t 3-+5 1 e t 2 是原方程的解。 3. dt ds =-s t cos + 21t 2sin 解:s=e ? -tdt cos (t 2sin 2 1 ?e dt dt ? 3c + ) =e t sin -(?+c dt te t t sin cos sin ) = e t sin -(c e te t t +-sin sin sin ) =1sin sin -+-t ce t 是原方程的解。 4. dx dy n x x e y n x =- , n 为常数. 解:原方程可化为: dx dy n x x e y n x += )(c dx e x e e y dx x n n x dx x n +??=?- )(c e x x n += 是原方程的解.
5. dx dy + 1212 --y x x =0 解:原方程可化为: dx dy =-1212 +-y x x ? =-dx x x e y 2 1 2(c dx e dx x x +? -2 21) ) 2 1(ln 2 + =x e )(1ln 2 ?+- -c dx e x x =)1(1 2 x ce x + 是原方程的解. 6. dx dy 2 3 4xy x x += 解: dx dy 2 3 4 xy x x += =2 3y x + x y 令 x y u = 则 ux y = dx dy =u dx du x + 因此:dx du x u += 2 u x 2 1u dx du = dx du u =2 c x u +=3 31 c x x u +=-33 (*) 将 x y u =带入 (*)中 得:3 4 3 3cx x y =-是原方程的解.
摘要 (2) 1.何谓奇解 (2) 2.奇解的产生 (3) 3.包络跟奇解的关系 (4) 4.理论上证明C-判别曲线与P-判别曲线方法 (5) 4.1 克莱罗微分方程 (9) 5.奇解的基本性质 (12) 5.1 定理1 (12) 5.2 定理2 (14) 5.3 定理3 (14) 6.小结 (14) 参考文献: (15)
一阶常微分方程的奇解 摘要 在常微分方程中,我们知道方程的解可以有多种,现在我们来讨论求奇解的方法。我们看到某些微分方程,会存在一些特殊的积分曲线,他并不属于这方程的积分曲线族,但是,在这些特殊的积分曲线上的每一点处,都有积分曲线族中的一条曲线和他在此处相切。从而我们引出了积分曲线族的包络,而为了求微分方程的奇解,,我们应先求出他的通解,然后求通解的包络。 关键词:奇解,包络,C-判别式,P-判别式 1.何谓奇解 设一阶隐式方程),,(,y y x F =0有一特解
)(:x y ψ=Γ,j x ∈ 如果对每一点Γ∈P ,在P 点的任何一个领域,方程),,(,y y x F =0都有一个不同于Γ的解在P 点与Γ相切,则称Γ是微分方程的),,(,y y x F =0的奇解 定义:如果一个一阶微分方程的一个特解的积分曲线上的每一点都至少和这个微分方程的不同的积分曲线相切,并且这相切的积分曲线在切点的任何邻域都不重合,则称这个特解为这个微分方程的奇解 2.奇解的产生 先看一个例子,求方程 033=-?? ? ??y dx dy (1) 或与它等价的方程 3y dx dy = 的解。 经分离变量后,可得(1)的通解 3)(27 1c x y += 容易看出,y=0也是原方程的一个解。现在来研究这个解y=0有什么特殊的地方。由图我们看到,在解y=0上的每一 点)0,(0x 处相切,这种特殊的积分曲线y=0 称为奇积分曲线,他所对应的解就是奇 解,这就是奇解的产生。 我们现在给出曲线族包络的定义 某些微分方程,存在一些特殊的积分 曲线,会存在一些特殊的积分曲线,他并 不属于这方程的积分曲线族,但是,在这些特殊的积分曲线上的每一点处,都有积分曲线族中的一条曲线和他在此处相切。在几何学里,这些特殊的积分曲线称为上述积分曲线族的包络,在微分方程里,这些特殊的积分曲线所对应的解酒称为方程的奇解。
二阶常微分方程解
————————————————————————————————作者: ————————————————————————————————日期:
第七节 二阶常系数线性微分方程 的解法 在上节我们已经讨论了二阶线性微分方程解的结构,二阶线性微分方程的求解问题,关键在于如何求二阶齐次方程的通解和非齐次方程的一个特解。本节讨论二阶线性方程的一个特殊类型,即二阶常系数线性微分方程及其求解方法。先讨论二阶常系数线性齐次方程的求解方法。 §7.1 二阶常系数线性齐次方程及其求解方法 设给定一常系数二阶线性齐次方程为 ?? 22 dx y d +p dx dy +qy=0 (7.1) 其中p 、q 是常数,由上节定理二知,要求方程(7.1)的通解,只要求出其任意两个线性无关的特解y 1,y2就可以了,下面讨论这样两个特解的求法。 我们先分析方程(7.1)可能具有什么形式的特解, 从方程的形式上来看,它的特点是22dx y d ,dx dy ,y 各乘以 常数因子后相加等于零,如果能找到一个函数y,其
22dx y d ,dx dy ,y之间只相差一个常数因子,这样的函数有可能是方程(7.1)的特解,在初等函数中,指数函数e rx ,符合上述要求,于是我们令 y=e r x (其中r 为待定常数)来试解 将y =e rx ,dx dy =re r x,22dx y d =r 2e r x 代入方程(7.1) 得 r 2e rx +pre rx +qerx =0 或 e r x(r 2+pr+q )=0 因为e rx ≠0,故得 ? r 2 +pr +q=0 由此可见,若r 是二次方程 ?? r 2+pr +q=0 (7.2) 的根,那么e r x就是方程(7.1)的特解,于是方程(7.1)的求解问题,就转化为求代数方程(7.2)的根问题。称(7.2)式为微分方程(7.1)的特征方程。 特征方程(7.2)是一个以r 为未知函数的一元二次代数方程。特征方程的两个根r 1,r 2,称为特征根,由代数知识,特征根r 1,r 2有三种可能的情况,下面我们分别进行讨论。 (1)若特证方程(7.2)有两个不相等的实根r 1, r 2,此时e r 1x ,e r2x 是方程(7.1)的两个特解。
常微分方程 2.1 1. xy dx dy 2=,并求满足初始条件:x=0,y=1的特解. 解:对原式进行变量分离得 。 故它的特解为代入得 把即两边同时积分得:e e x x y c y x x c y c y xdx dy y 2 2 ,11,0,ln ,21 2 =====+== ,0)1(.22 =++dy x dx y 并求满足初始条件:x=0,y=1的特解. 解:对原式进行变量分离得: 。 故特解是 时,代入式子得。当时显然也是原方程的解当即时,两边同时积分得;当x y c y x y x c y c y x y dy dx x y ++=====++=+=+≠=+- 1ln 11 ,11,001ln 1 ,11ln 0,1112 3 y xy dx dy x y 32 1++ = 解:原式可化为: x x y x x y x y x y y x y c c c c x dx x dy y y x y dx dy 2 2 2 2 22 2 2 3 22 3 2 )1(1)1)(1(),0(ln 1ln 21ln 1ln 2 1 1 1,0111=++ =++ ≠++-=+ +=+≠+ ? + =+) 故原方程的解为(即两边积分得故分离变量得显然 .0;0;ln ,ln ,ln ln 0 110000 )1()1(4===-==-+=-++=-=+≠===-++x y c y x xy c y x xy c y y x x dy y y dx x x xy x y xdy y ydx x 故原方程的解为即两边积分时,变量分离是方程的解,当或解:由:
10ln 1ln ln 1ln 1,0 ln 0 )ln (ln :931:8. cos ln sin ln 0 7ln sgn arcsin ln sgn arcsin 1 sgn 11,)1(,,,6ln )1ln(2 11 11,11,,,0 )()(:5332 2 22 2 22 2 22 2 c dx dy dx dy x y cy u d u u dx x x y u dx x y dy x y ydx dy y x x c dy y y y y dx dy c x y tgxdx ctgydy ctgxdy tgydx c x x x y c x x u dx x x du x dx du dx du x u dx dy ux y u x y y dx dy x c x arctgu dx x du u u u dx du x u dx du x u dx dy ux y u x y x y x y dx dy dx x y dy x y e e e e e e e e x y u u x y x u u x y x y y x x x +===+=+-===-?-=--+-=-=+-===-=+?=+?=?=--=+===-+=+-=++ =++-++=++===+-==-++-+-- 两边积分解:变量分离:。 代回原变量得:则有:令解:方程可变为:解:变量分离,得 两边积分得:解:变量分离,得::也是方程的解。 另外,代回原来变量,得两边积分得:分离变量得:则原方程化为: 解:令:。两边积分得:变量分离,得:则令解:
常微分方程(第三版) 答案
常微分方程习题答案 2.1 1.?Skip Record If...?,并求满足初始条件:x=0,y=1的特解. 解:对原式进行变量分离得 ?Skip Record If...??Skip Record If...?并求满足初始条件:x=0,y=1的特解. 解:对原式进行变量分离得: ?Skip Record If...?3 ?Skip Record If...? 解:原式可化为: ?Skip Record If...??Skip Record If...??Skip Record If...? ?Skip Record If...? 12.?Skip Record If...? 解?Skip Record If...??Skip Record If...? ?Skip Record If...? 15.?Skip Record If...? ?Skip Record If...?16.?Skip Record If...? 解:?Skip Record If...? ?Skip Record If...?,这是齐次方程,令?Skip Record If...? 17. ?Skip Record If...? 解:原方程化为?Skip Record If...? 令?Skip Record If...? 方程组?Skip Record If...??Skip Record If...? 则有?Skip Record If...? 令?Skip Record If...? 当?Skip Record If...?当?Skip Record If...? 另外 ?Skip Record If...? ?Skip Record If...?
习题3.1 1 求方程dx dy =x+y 2通过点(0,0)的第三次近似解; 解: 取0)(0=x ? 20020012 1)()(x xdx dx y x y x x x ==++=??? 522200210220 121])21([])([)(x x dx x x dx x x y x x x +=+=++=???? dx x x x y x x ])20 121([)(252003+++=?? = 118524400 1160120121x x x x +++ 2 求方程dx dy =x-y 2通过点(1,0)的第三次近似解; 解: 令0)(0=x ? 则 20020012 1)()(x xdx dx y x y x x x ==-+=??? 522200210220 121])21([])([)(x x dx x x dx x x y x x x -=-=-+=???? dx x x x y x x ])20 121([)(252003--+=?? =118524400 1160120121x x x x -+- 3 题 求初值问题: ?????=--=0 )1(22y y x dx dy R :1+x ≤1,y ≤1 的解的存在区间,并求解第二次近似解,给出在解的存在空间的误差估计; 解: 因为 M=max{22y x -}=4 则h=min(a,M b )=4 1 则解的存在区间为0x x -=)1(--x =1+x ≤4 1 令 )(0X ψ=0 ; )(1x ψ=y 0+?-x x x 0)0(2dx=31x 3+31;
)(2x ψ =y 0+])3131([2132?-+-x x x dx=31x 3-9x -184x -637x +4211 又 y y x f ??),(2≤=L 则:误差估计为:)()(2x x ψ-ψ≤32 2 )12(*h L M +=2411 4 题 讨论方程:31 23y dx dy =在怎样的区域中满足解的存在唯一性定理的条件, 并求通过点(0,0)的一切解; 解:因为y y x f ??),(=3221-y 在y 0≠上存在且连续; 而312 3y 在y 0 σ≥上连续 由 3123y dx dy =有:y =(x+c )23 又 因为y(0)=0 所以:y =x 2 3 另外 y=0也是方程的解; 故 方程的解为:y =?????≥00023 x x x 或 y=0; 6题 证明格朗瓦耳不等式: 设K 为非负整数,f(t)和g(t)为区间βα≤≤t 上的连续非负函数,
习题1.2 1. dx dy =2xy,并满足初始条件:x=0,y=1的特解。 解: y dy =2xdx 两边积分有:ln|y|=x 2+c y=e 2 x +e c =cex 2 另外y=0也是原方程的解,c=0时,y=0 原方程的通解为y= cex 2,x=0 y=1时 c=1 特解为y= e 2 x . 2. y 2dx+(x+1)dy=0 并求满足初始条件:x=0,y=1的特解。 解:y 2dx=-(x+1)dy 2y dy dy=-1 1+x dx 两边积分: - y 1 =-ln|x+1|+ln|c| y=|)1(|ln 1+x c 另外y=0,x=-1也是原方程的解 x=0,y=1时 c=e 特解:y= | )1(|ln 1 +x c 3.dx dy =y x xy y 321++ 解:原方程为:dx dy =y y 21+3 1 x x + y y 21+dy=3 1 x x +dx 两边积分:x(1+x 2 )(1+y 2 )=cx 2 4. (1+x)ydx+(1-y)xdy=0 解:原方程为: y y -1dy=-x x 1 +dx 两边积分:ln|xy|+x-y=c 另外 x=0,y=0也是原方程的解。 5.(y+x )dy+(x-y)dx=0
解:原方程为: dx dy =-y x y x +- 令 x y =u 则dx dy =u+x dx du 代入有: -1 12++u u du=x 1dx ln(u 2+1)x 2=c-2arctgu 即 ln(y 2+x 2)=c-2arctg 2x y . 6. x dx dy -y+22y x -=0 解:原方程为: dx dy =x y +x x | |-2)(1x y - 则令 x y =u dx dy =u+ x dx du 2 11u - du=sgnx x 1 dx arcsin x y =sgnx ln|x|+c 7. tgydx-ctgxdy=0 解:原方程为: tgy dy =ctgx dx 两边积分:ln|siny|=-ln|cosx|-ln|c| siny= x c cos 1=x c cos 另外y=0也是原方程的解,而c=0时,y=0. 所以原方程的通解为sinycosx=c. 8 dx dy +y e x y 32 +=0 解:原方程为:dx dy =y e y 2 e x 3 2 e x 3-3e 2 y -=c. 9.x(lnx-lny)dy-ydx=0 解:原方程为: dx dy =x y ln x y
第四章常微分方程数值解 [课时安排]6学时 [教学课型]理论课 [教学目的和要求] 了解常微分方程初值问题数值解法的一些基本概念,如单步法和多步法,显式和隐式,方法的阶数,整体截断误差和局部截断误差的区别和关系等;掌握一阶常微分方程初值问题的一些常用的数值计算方法,例如欧拉(Euler)方法、改进的欧拉方法、龙贝-库塔(Runge-Kutta)方法、阿达姆斯(Adams)方法等,要注意各方法的特点及有关的理论分析;掌握构造常微分方程数值解的数值积分的构造方法和泰勒展开的构造方法的基本思想,并能具体应用它们导出一些常用的数值计算公式及评估截断误差;熟练掌握龙格-库塔(R-K)方法的基本思想,公式的推导,R-K公式中系数的确定,特别是能应用“标准四阶R-K公式”解题;掌握数值方法的收敛性和稳定性的概念,并能确定给定方法的绝对稳定性区域。 [教学重点与难点] 重点:欧拉方法,改进的欧拉方法,龙贝-库塔方法。 难点:R—K方法,预估-校正公式。 [教学内容与过程] 4.1 引言 本章讨论常微分方程初值问题 (4.1.1) 的数值解法,这也是科学与工程计算经常遇到的问题,由于只有很特殊的方程能用解析方法求解,而用计算机求解常微分方程的初值问题都要采用数值方法.通常我们假定(4.1.1)中 f(x,y)对y满足Lipschitz条件,即存在常数L>0,使对,有 (4.1.2) 则初值问题(4.1.1)的解存在唯一. 假定(4.1.1)的精确解为,求它的数值解就是要在区间上的一组离散点 上求的近似.通常取 ,h称为步长,求(4.1.1)的数值解是按节点的顺序逐步
推进求得.首先,要对方程做离散逼近,求出数值解的公式,再研究公式的局部截断误差,计算稳定性以及数值解的收敛性与整体误差等问题. 4.2 简单的单步法及基本概念 4.2.1 Euler法、后退Euler法与梯形法 求初值问题(4.1.1)的一种最简单方法是将节点的导数用差商 代替,于是(4.1.1)的方程可近似写成 (4.2.1) 从出发,由(4.2.1)求得再将 代入(4.2.1)右端,得到的近似,一般写成 (4.2.2) 称为解初值问题的Euler法. Euler法的几何意义如图4-1所示.初值问题(4.1.1)的解曲线y=y(x)过点,从出发,以为斜率作一段直线,与直线交点于,显然有 ,再从出发,以为斜率作直线推进到上一点,其余类推,这样得到解曲线的一条近似曲线,它就是折线.
习题2.2 求下列方程的解 1.dx dy =x y sin + 解: y=e ?dx (?x sin e ?-dx c dx +) =e x [- 2 1e x -(x x cos sin +)+c] =c e x -21 (x x cos sin +)是原方程的解。 2.dt dx +3x=e t 2 解:原方程可化为: dt dx =-3x+e t 2 所以:x=e ?-dt 3 (?e t 2 e -? -dt 3c dt +) =e t 3- (5 1e t 5+c) =c e t 3-+5 1e t 2 是原方程的解。 3.dt ds =-s t cos +21t 2sin 解:s=e ?-tdt cos (t 2sin 2 1?e dt dt ?3c + ) =e t sin -(?+c dt te t t sin cos sin ) = e t sin -(c e te t t +-sin sin sin ) =1sin sin -+-t ce t 是原方程的解。 4. dx dy n x x e y n x =- , n 为常数. 解:原方程可化为:dx dy n x x e y n x += )(c dx e x e e y dx x n n x dx x n +??=?- )(c e x x n += 是原方程的解.
5. dx dy +1212--y x x =0 解:原方程可化为:dx dy =-1212+-y x x ?=-dx x x e y 21 2(c dx e dx x x +?-221) )21(ln 2+=x e )(1 ln 2?+--c dx e x x =)1(1 2 x ce x + 是原方程的解. 6. dx dy 234xy x x += 解:dx dy 234xy x x += =23y x +x y 令 x y u = 则 ux y = d x d y =u dx du x + 因此:dx du x u +=2u x 21u dx du = dx du u =2 c x u +=33 1 c x x u +=-33 (*) 将x y u =带入 (*)中 得:3433cx x y =-是原方程的解.
常微分方程第三版答 案
习题1.2 1. dx dy =2xy,并满足初始条件:x=0,y=1的特解。 解: y dy =2xdx 两边积分有:ln|y|=x 2+c y=e 2 x +e c =cex 2另外y=0也是原方程的解,c=0时,y=0 原方程的通解为y= cex 2,x=0 y=1时 c=1 特解为y= e 2 x . 2. y 2dx+(x+1)dy=0 并求满足初始条件:x=0,y=1的特解。 解:y 2dx=-(x+1)dy 2y dy dy=-1 1+x dx 两边积分: - y 1 =-ln|x+1|+ln|c| y=|)1(|ln 1+x c 另外y=0,x=-1也是原方程的解 x=0,y=1时 c=e 特解:y= | )1(|ln 1 +x c 3.dx dy =y x xy y 321++ 解:原方程为:dx dy =y y 21+3 1 x x + y y 21+dy=3 1 x x +dx 两边积分:x(1+x 2)(1+y 2)=cx 2 4. (1+x)ydx+(1-y)xdy=0 解:原方程为: y y -1dy=-x x 1 +dx 两边积分:ln|xy|+x-y=c 另外 x=0,y=0也是原方程的解。
5.(y+x )dy+(x-y)dx=0 解:原方程为: dx dy =-y x y x +- 令 x y =u 则dx dy =u+x dx du 代入有: -1 12++u u du=x 1dx ln(u 2+1)x 2=c-2arctgu 即 ln(y 2+x 2)=c-2arctg 2 x y . 6. x dx dy -y+22y x -=0 解:原方程为: dx dy =x y +x x | |-2)(1x y - 则令 x y =u dx dy =u+ x dx du 2 11u - du=sgnx x 1 dx arcsin x y =sgnx ln|x|+c 7. tgydx-ctgxdy=0 解:原方程为: tgy dy =ctgx dx 两边积分:ln|siny|=-ln|cosx|-ln|c| siny= x c cos 1=x c cos 另外y=0也是原方程的解,而c=0时,y=0. 所以原方程的通解为sinycosx=c. 8 dx dy +y e x y 32 +=0 解:原方程为:dx dy =y e y 2 e x 3 2 e x 3-3e 2 y -=c.
常微分方程习题及解答
常微分方程习题及解答 一、问答题: 1. 常微分方程和偏微分方程有什么区别?微分方程的通解是什么含义? 答:微分方程就是联系着自变量,未知函数及其导数的关系式。常微分方程,自变量的个数只有一个。偏微分方程,自变量的个数为两个或两个以上。常微分方程解的表达式中,可能包含一个或几个任意常数,若其所包含的独立的任意常数的个数恰好与该方程的阶数相同,这样的解为该微分方程的通解。 2. 举例阐述常数变易法的基本思想。 答:常数变易法用来求线性非齐次方程的通解,是将线性齐次方程通解中的任意常数变易为待定函数来求线性非齐次方程的通解。 例:求()()dy P x y Q x dx =+的通解。 首先利用变量分离法可求得其对应的线性齐次方程的通解为()P x dx y c ?=l ,然后将常 数c 变易为x 的待定函数()c x ,令()()P x dx y c x ? =l , 微分之,得到 ()()()()()P x dx P x dx dy dc x c x P x dx dx ??=+l l ,将上述两式代入
方程中,得到 ()()()()()()()()() P x dx P x dx P x dx dc x c x P x dx c x P x Q x ??+?=+l l l 即 ()() ()P x dx dc x Q x dx -? =l 积分后得到()()()P x dx c x Q x dx c -? =+?%l 进而得到方程 的通解 ()()(()) P x dx P x dx y Q x dx c -? ?=+?%l l 3.高阶线性微分方程和线性方程组之间的联系如何? 答:n 阶线性微分方程的初值问题 ()(1) 11(1)01020()...()()()(),(),....()n n n n n n x a t x a t x a t x f t x t x t x t ηηη---'?++++=??'===?? 其中1 2 ()(),...(),()n a t a t a t f t ,是区间a t b ≤≤上的已知 连续函数,[]0 ,t a b ∈,1 2 ,,...,n ηηη是已知常数。 它可以化为线性微分方程组的初值问题 12100100 00010000010()()()()()()n n n x x a t a t a t a t f t x t η--????????????????'????=+?????? ???? ? ?????----????? =?? L L M M M M M M L L 但是需要指出的是每一个n 阶线性微分方程可化为n 个一阶线性微分方程构成的方程组,反之却不成立。 4.若常系数线性方程组 Ax x ='和Bx x ='有相同的基本解矩阵, 则A
习题1.2 1. dx dy =2xy,并满足初始条件:x=0,y=1的特解。 解: y dy =2xdx 两边积分有:ln|y|=x 2+c y=e 2 x +e c =cex 2 另外y=0也是原方程的解,c=0时,y=0 原方程的通解为y= cex 2,x=0 y=1时 c=1 特解为y= e 2 x . 2. y 2 dx+(x+1)dy=0 并求满足初始条件:x=0,y=1的特解。 解:y 2dx=-(x+1)dy 2 y dy dy=-1 1+x dx 两边积分: - y 1 =-ln|x+1|+ln|c| y=|)1(|ln 1+x c 另外y=0,x=-1也是原方程的解 x=0,y=1时 c=e 特解:y= | )1(|ln 1 +x c 3.dx dy =y x xy y 321++ 解:原方程为:dx dy =y y 21+31 x x + y y 21+dy=31 x x +dx 两边积分:x(1+x 2 )(1+y 2 )=cx 2 4. (1+x)ydx+(1-y)xdy=0 解:原方程为: y y -1dy=-x x 1 +dx 两边积分:ln|xy|+x-y=c 另外 x=0,y=0也是原方程的解。 5.(y+x )dy+(x-y)dx=0 解:原方程为:
dx dy =-y x y x +- 令 x y =u 则dx dy =u+x dx du 代入有: -1 12++u u du=x 1dx ln(u 2+1)x 2=c-2arctgu 即 ln(y 2+x 2)=c-2arctg 2 x y . 6. x dx dy -y+22y x -=0 解:原方程为: dx dy =x y +x x | |-2)(1x y - 则令 x y =u dx dy =u+ x dx du 2 11u - du=sgnx x 1 dx arcsin x y =sgnx ln|x|+c 7. tgydx-ctgxdy=0 解:原方程为: tgy dy =ctgx dx 两边积分:ln|siny|=-ln|cosx|-ln|c| siny= x c cos 1=x c cos 另外y=0也是原方程的解,而c=0时,y=0. 所以原方程的通解为sinycosx=c. 8 dx dy +y e x y 32+=0 解:原方程为:dx dy =y e y 2 e x 3 2 e x 3-3e 2 y -=c. 9.x(lnx-lny)dy-ydx=0 解:原方程为: dx dy =x y ln x y 令x y =u ,则dx dy =u+ x dx du
\ 习题 1. dx dy =2xy,并满足初始条件:x=0,y=1的特解。 解: y dy =2xdx 两边积分有:ln|y|=x 2+c y=e 2 x +e c =cex 2 另外y=0也是原方程的解,c=0时,y=0 原方程的通解为y= cex 2,x=0 y=1时 c=1 特解为y= e 2 x . 2. y 2 dx+(x+1)dy=0 并求满足初始条件:x=0,y=1的特解。 & 解:y 2 dx=-(x+1)dy 2y dy dy=-1 1+x dx 两边积分: - y 1 =-ln|x+1|+ln|c| y=|)1(|ln 1+x c 另外y=0,x=-1也是原方程的解 x=0,y=1时 c=e 特解:y= | )1(|ln 1 +x c 3.dx dy =y x xy y 321++ 解:原方程为:dx dy =y y 21+3 1 x x + y y 21+dy=3 1 x x +dx 两边积分:x(1+x 2)(1+y 2)=cx 2 、 4. (1+x)ydx+(1-y)xdy=0 解:原方程为: y y -1dy=-x x 1 +dx 两边积分:ln|xy|+x-y=c 另外 x=0,y=0也是原方程的解。 5.(y+x )dy+(x-y)dx=0
解:原方程为: dx dy =-y x y x +- { 令 x y =u 则dx dy =u+x dx du 代入有: -1 12++u u du=x 1dx ln(u 2+1)x 2=c-2arctgu 即 ln(y 2+x 2)=c-2arctg 2x y . 6. x dx dy -y+2 2y x -=0 解:原方程为: dx dy =x y +x x | |-2)(1x y - 则令 x y =u dx dy =u+ x dx du 2 11u - du=sgnx x 1 dx } arcsin x y =sgnx ln|x|+c 7. tgydx-ctgxdy=0 解:原方程为: tgy dy =ctgx dx 两边积分:ln|siny|=-ln|cosx|-ln|c| siny= x c cos 1=x c cos 另外y=0也是原方程的解,而c=0时,y=0. 所以原方程的通解为sinycosx=c. 8 dx dy +y e x y 32 +=0 解:原方程为:dx dy =y e y 2 e x 3 ` 2 e x 3-3e 2 y -=c. (lnx-lny)dy-ydx=0