第七章 常微分方程

合集下载

(整理)第七章常微分方程数值解1

(整理)第七章常微分方程数值解1

第七章常微分方程数值解本章介绍求解微分方程数值解的基本思想和方法.●常微分方程含有自变量、未知函数和它的一阶导数和高阶导数的方程.它是描述运动、变化规律的重要数学方法之一,分为两类:1.初值问题,即给出未知函数及导数在初始点的值2.边值问题,即给出未知函数及(或)它的某些导数在区间两个端点的值●常微分方程初值问题00(,),, ||(),,y f x y a x b y y a y y R '=<≤<∞⎧⎨=∈⎩ 其中f 为,x y 的已知函数,0y 为给定的初值. 这里仅讨论一阶标量微分方程初值问题的数值解法.而高阶微分方程通常可化为一阶微分方程组来研究.数值解法寻求微分方程初值问题之解()y x 在一系列离散点012N a x x x x b =<<<<=上的近似值:,,,210y y yN y , 的方法.{}n y : 问题的数值解数值解所满足的离散方程统称为差分格式. 步长:1i i i h x x -=- ,一般取定步长i h h =● 初值问题的适定性(其解是否唯一存在)记(带形)区域:[,]a b R ⨯为G ,即G =[,]a b R ⨯. 设:f G R →为连续映射,若存在常数0L >使得不等式1212|(,)(,)|||f x y f x y L y y -≤- 对一切12(,),(,)x y x y G ∈都成立,则称(,)f x y 在G 上关于y 满足Lipschitz 条件,而式中的常数L 称为Lipschitz 常数.定理 初值问题00(,),, ||(),,y f x y a x b y y a y y R '=<≤<∞⎧⎨=∈⎩,当(,)f x y 在G 上连续,且关于y 满足Lipschitz 条件,则其解存在且唯一.§ 7.1 Euler 方法● Euler 公式将初值问题⎩⎨⎧∈=∞<≤<=',,)(||,),,(00R y y a y y b x a y x f y 的求解区间],[b a N 等分,分点:),,2,1,0(,N n nh a x n =+= ,其中N ab h -=将),(y x f y ='写成等价的积分方程形式:⎰++=+hx x d y f x y h x y τττ))(,()()(,在上式中令n x x =,并用左矩形公式计算右端积分,得到n n n n n R x y x hf x y h x y ++=+))(,()()(, (*)⎰+-=1))(,())(,(n nx x n n n x y x hf dx x y x f R ——余项将(*)中的余项n R 截去,可得))(,()()(1n n n n x y x hf x y x y +≈+则有)(n x y 的近似值n y 的递推公式1,,1,0),,(1-=+=+N n y x hf y y n n n n (*1)----------Euler 公式n R 表示当)(n n x y y =为精确值时,利用(*1)即Euler 公式计算)(1+n x y 时的误差.n n n y x y -=)(ε 为Euler 方法的局部截断误差.例1 用Euler 公式解初值问题⎪⎩⎪⎨⎧=≤<-='1)0(10,2y x y x y y解:取1.0=h ,Euler 公式的具体形式为)2(),(1nnn n n n n n y x y h y y x hf y y -+=+=+其中)10,,1,0,1.0 ===n n nh x n ( 已知10=y ,则有1.11.01)2(0001=+=-+=y x y h y y191818.1)1.12.01.1(1.01.1)2(11112=-+=-+=y x y h y y… … … 依次计算可得109876543,,,,,,,y y y y y y y y其部分结果见下表可见Euler 方法的计算结果精度不太高。

第七章 常微分方程习题课

第七章 常微分方程习题课

7
x3
2
Cx3 .
7
20
例3 求 dy dx
x
y y2 cos
的通解 y
解:将原方程写成
dx 1 x y cos y dy y
x
e
1 y
dy
(
y cos
y
)e
1 y
dy
dy
C
y
(
y
cos
y)
1 y
dy
C
y(C sin y)
21
例4
求通解
2x
y2 3x2
dx
dy 0.
y3
y4
3
2、一阶微分方程的解法
(1) 可分离变量的微分方程
形如 g( y)dy f ( x)dx
解法 g( y)dy f ( x)dx
分离变量法
(2) 齐次方程 形如 dy f ( y) dx x
解法 作变量代换 u y x
4
(3) 可化为齐次的方程
形如 dy f ( ax by c )
通解:y Y y*
求特解的方法 待定系数法
(1) f ( x) ex Pm ( x) 型
0 不是根
设特解 y* xkexQm (x) k 1 是单根 ,
2 是重根
15
(2)
f
(
x)
ex
[
Pl
(1)
(
x)
cos
x
P(2) n
(
x)
sin
x]

设特解
y*
xkex[Rm(1) (x) cosx

2 解得
4b 0,
a 1,
8 b 0,

第七章常微分方程数值解法

第七章常微分方程数值解法

h2 h3 y ( xi 1 ) y ( xi h) y ( xi ) hy '( xi ) y ''( xi ) y '''( xi ) 2! 3!
丢掉高阶项,有
y( xi 1 ) y( xi h) y( xi ) hy '( xi ) yi hf ( xi , yi )
| f ( x, y1 ) f ( x, y2 ) | L | y1 y2 | ,
那么模型问题在 [ a, b] 存在唯一解。
Lipschitz 连续: | f ( x, y1 ) f ( x, y2 ) | L | y1 y2 | .
(1) 比连续性强: y1 y2 可推出 f ( x, y1 ) f ( x, y2 ) ; (2) 比连续的 1 阶导弱:具有连续的 1 阶导,则
f | f ( x, y1 ) f ( x, y2 ) || ( ) || y1 y2 | L | y1 y2 | . y
常微分方程数值解法
目标:计算出解析解 y ( x) 在一系列节点 a x0 x1 xn1 xn b 处的近似值 yi y( xi ) ,即所谓的数值解。节点间距 hi xi 1 xi ,一般 取为等距节点。
常微分方程初值问题的数值解法一般分为两大类: (1)单步法:在计算 yn 1 时,只用到前一步的值,即用到 xn1 , xn , yn ,则给定初
值之后,就可逐步计算。例如 Euler 法、向后欧拉法、梯形公式、龙格-库塔法;
(2) 多步法: 这 类 方 法 在 计算 yn 1 时 , 除 了 用 到 xn1 , xn , yn 外 , 还 要 用到

微分方程

微分方程

dy P ( x ) y Q( x ) dx
dy 2 dx 2 例如 y x , x sin t t , 线性的; dx dt
yy 2 xy 3, y cos y 1,
非线性的.
高等数学(上)
一阶线性非齐次微分方程的通解为:
ye
Ce
P ( x ) dx
过定点的积分曲线; 微分方程的图形
y f ( x , y , y ) 二阶: y x x0 y0 , y x x0 y0
过定点且在定点的切线的斜率为定值的积分曲线.
高等数学(上)
第二节 一阶微分方程
一、可分离变量的微分方程
二、齐次方程
三、一阶线性微分方程
cos x C.
所以原方程通解为
y
1 cos x C . x
高等数学(上)
1 sin x 求方程 y y 的通解. x x
1 解 P( x) , x
sin x Q( x ) , x
sin x y x ln x sin x ln x e e dx C x 1 1 sin xdx C cos x C . x x
高等数学(上)
( x, C1 )
例3 求方程 xy

(5)
y
(4)
0 的通解.
(5)
设y
(4)
P ( x ), y
P ( x )
(4)
代入原方程 分离变量,得
xP P 0, (P 0)
1 2 两端积分,得 y C1 x C 2 , 2
原方程通解为
高等数学(上)

第七章常微分方程练习题(含答案)

第七章常微分方程练习题(含答案)

第7章 常微分方程一、单项选择题1.微分方程3245(''')3('')(')0y y y x -++=阶数是( b )A.4阶 B .3阶 C .2阶 D .1阶2.微分方程222y x dxdy x +=是( b ) A.一阶可分离变量方程 B.一阶齐次方程 C.一阶非齐次线性方程 D.一阶齐次线性方程3.下列方程中,是一阶线性微分方程的是( c )A.0'2)'(2=+-x yy y xB.0'2=-+x yy xyC.0'2=+y x xyD.0)()67(=++-dy y x dx y x4.方程x y xy =-'满足初始条件11==x y 的特解是( a )A.x x x y +=lnB.Cx x x y +=lnC.x x x y +=ln 2D.Cx x x y +=ln 25.微分方程y y x 2='的通解为( c )A .2x y =B . c x y +=2C . 2cx y =D .0=y6.微分方程y y x ='满足1)1(=y 的特解为 ( a )A.x y =B. c x y +=C.cx y =D.0=y8.微分方程05))(sin(2''=+-+x y y xy y 是( a )A 一阶微分方程B 二阶微分方程C 可分离变量的微分方程D 一阶线性微分方程9.微分方程2y xy '=的通解为( c )A .2x y e C =+B . x y Ce =C . 2x y Ce =D .22x y Ce =二、填空题1.微分方程34()"30y y y y '++=的阶数为__2____;2.微分方程0=+y dxdy 的通解是x y ce -=; 3.微分方程02=+'xy y 的通解是2x y ce -=;4.微分方程x y y e +'=的通解是()10,0x ye C e C ++=<; 5. 一阶线性微分方程()()y P x y Q x '+=的通解为()()()()P x dx P x dx P x dx y Ce e Q x e dx --⎰⎰⎰=+⎰; 6. n 阶微分方程的通解含有__n __个独立的任意常数。

常微分方程简介

常微分方程简介

第七章 常微分方程简介我们已经学完一元函数微积分的基本内容.回顾微积分的产生和发展,就会发现它与人们求解微分方程的需要有密切关系.20世纪以前,微分方程问题主要来源于几何学、力学和物理学,而现在它几乎渗透到自然科学和一些社会科学的各个领域,已成为人们研究科学技术,解决实际问题的不可缺少的有力工具.本章我们主要介绍常微分方程的基本概念,一阶微分方程的初等解法,可降阶的高阶方程及常系数线性方程的求解方法,它是本课程的一个重要组成部分.§7.1 基本概念1. 微分方程及其解的定义利用数学手段研究自然现象和社会现象,或解决工程技术问题,一般先要建立数学模型,再对数学模型进行简化和求解,最后结合实际问题对结果进行分析和讨论.数学模型最常见的表达方式是包含自变量和未知函数的方程,在很多情况下未知函数的导数(或微分)也会在方程中出现,于是便自然地称这类方程为微分方程.定义7.1.1 联系着自变量、未知函数及其某些导数的方程称为微分方程. 只含一个自变量的方程称为常微分方程,自变量多于一个的称为偏微分方程.微分方程中实际出现的导数的最高阶数称为微分方程的阶.于是n 阶常微分方程的一般形式是0),,,,()(='n y y y x F , (1.1) 其中F 是2+n 个变元的已知函数,且)(n y 一定出现.(注意,这里我们仅引用了多元函数的记号,它是一元函数记号在形式上的推广.)本章只介绍常微分方程,并简称为微分方程或方程.定义7.1.2 如果方程(1.1)的左边函数F 对未知函数y 和它的各阶导数)(,n y y '的全体而言是一次的,则称它为线性微分方程,否则称它为非线性微分方程.n 阶线性微分方程的一般形式是:)()()()(1)1(1)(x f y x a y x a y x a y n n n n =+'++-- , (1.2)其中)(x a i ),,2 ,1(n i =和)(x f 都是x 的已知函数.例如,下面的方程都是常微分方程:yx y -=', (1.3) 21y y +=', (1.4) 02=+''y y ω (0>ω是常数), (1.5) 它们的阶数分别为1,1,2.方程(1.5)是线性的,而方程(1.3)和(1.4)是非线性的.定义7.1.3 设函数)(x y ϕ=在区间I 上连续,且有直到n 阶的导数,若把)(x y ϕ=及其相应的各阶导数代入方程(1.1),得到关于x 的恒等式,即在I 上0))(),(),(,()(≡'x x x x F n ϕϕϕ ,则称)(x y ϕ=为方程(1.1)在区间I 上的解,若由关系式0),(=y x φ所确定的隐函数是方程(1.1)的解,则称0),(=y x φ为方程(1.1)的隐式解.例如,从定义7.1.3可以直接验证:1) 函数21x y -=和21x y --=都是方程(1.3)在区间)1 ,1(-上的解,而122=+y x 是它的隐式解.2) 函数x y tan =是方程(1.4)在区间)2,2(ππ-上的一个解,而)tan(c x y -=是方程(1.4)在区间)2,2(ππ+-c c 上的解,其中c 为任意常数.3) 函数x y cos 3ω=,x y sin 4ω=都是方程(1.5)在区间),(+∞-∞上的解,而且对任意常数1c 和2c ,x c x c y sin cos 21ωω+=也是方程(1.5)在区间),(+∞-∞上的解.今后对解与隐式解不加区别,统称它们为解.一般情况下也不再指明解的定义区间.从上面的讨论可知,微分方程的解可以包含一个或几个任意常数(与方程的阶数有关),而有的解不含任意常数.为了加以区别,我们给出如下定义:定义7.1.4 方程(1.1)的含有n 个独立的任意常数n c c c ,, ,21 的解) ,, ,,(21n c c c x y ϕ=称为它的通解.不含任意常数的解称为它的特解.这里说n 个任意常数是独立的,其含义是指它们不能合并而使得任意常数的个数减少.例如对于两个任意常数的情形,设函数)(),(x x ψϕ在区间I 上连续,若在I 上≠)()(x x ψϕ常数或≠)()(x x ϕψ常数,则称函数)(),(x x ψϕ在I 上线性无关,这时易知表达式 )()(21x c x c y ψϕ+=中的两个任意常数21,c c 是独立的.例1 验证函数x c x c y sin cos 21ωω+=是方程(1.5)的通解,其中21,c c 为任意常数. 解 x c x c y c o s s i n 21ωωωω+-=',x c x c y s i n c o s 2221ωωωω--='',将y y '',的表达式代入方程(1.5)有) sin cos ( sin cos 21222212x c x c x c x c y y ωωωωωωωω++--=+''0≡, ),(+∞-∞∈x ,所以对任意常数21,c c ,x c x c y sin cos 21ωω+=都是方程(1.5)的解,又由于 ≠xx s i n c o s ωω常数 (πk x ≠,∈k Z ) 即21,c c 是两个独立的任意常数,因此x c x c y sin cos 21ωω+=是方程(1.5)的通解.类似验证 ) sin(B x A y +=ω (A,B 为任意常数)也是方程(1.5)的通解.而x y cos 3ω=和x y sin 4ω=则是方程(1.5)的两个特解.定义7.1.5 为了确定方程(1.1)的特解而给出的附加条件称为定解条件,求方程(1.1)的满足定解条件的特解的问题称为定解问题.方程(1.1)的一种常用的定解条件是初始条件,它的一般提法是)(0x y ,)1(00)(y x y =',…, )1(00)1()(--=n n y x y , (1.6) 其中0x , 0y ,)1(0y ,…,)1(0-n y 是任给的1+n 个常数.求方程(1.1)满足初始条件(1.6)的解的问题称为初值问题或柯西(Cauchy)问题. 例如x y cos 3ω=是初值问题⎩⎨⎧='==+''0)0(,3)0(02y y y y ω 的解,而x y sin 4ω=是初值问题⎩⎨⎧='==+''ωω4)0(,3)0(02y y y y的解.它们都是在求得方程的通解以后,再利用初始条件定出通解中的任意常数而得出.这种做法是具有一般性的.可以证明:对于在一定范围内给出的1+n 个常数:0x , 0y ,)1(0y ,…,)1(0-n y ,利用通解表达式及初始条件(1.6)便可确定通解中的n 个任意常数n c c c ,,,21 ,从而得到相应的初值问题的解.换句话说,在一定范围内,通解包含了方程的所有解,这也是通解这一名词的一种名副其实的解释.2. 微分方程及其解的几何解释考虑一阶微方程),(y x f y =' (1.7) 其中),(y x f 是平面区域D 内给定的连续函数.方程(1.7)的解)(x y ϕ=)(I x ∈在平面上的图形是一条光滑曲线,称它为方程(1.7)的一条积分曲线,记作Γ.任取一点Γ∈),(000y x P ,即I x ∈0,)(00x y ϕ=.由于)(x y ϕ=满足方程(1.7),故按导数的几何意义可知,曲线Γ在点0P 的切线斜率为),())(,()(00000y x f x x f x =='ϕϕ.这说明曲线Γ上任一点处的切线斜率恰好等于方程右边函数),(y x f 在该点的函数值.这样,在区域D 内每一点),(y x P ,都可以作一个以函数),(y x f 在该点的值为斜率的小线段来表明积分曲线(如果存在的话)在该点的切线方向.区域D 连同所有这些小线段称为方程(1.7)的方向场.现在我们可以对微分方程(1.7)及其解作出几何解释:给定方程(1.7),就相当于给定平面区域D 内的一个方向场,反之给定区域D 内的一个方向场,就相当于给定一个形如(1.7)的方程.方程(1.7)的解所对应的积分曲线就是区域D 内这样的一条曲线,在它所经过的每一点都与方向场吻合,即曲线上每一点的切线方向都与方向场在该点的方向一致.求解初值问题⎩⎨⎧=='00)(),(y x y y x f y , 就是求一条经过点),(00y x 并与方向场吻合的光滑曲线.以上这种几何解释,无论在理论上还是在实用上都有很大的价值.从理论上说,它把作为解析对象的微分方程及其解与作为几何对象的方向场及积分曲线沟通起来,从而在微分方程这门学科建立了数与形的联系,这就为我们从几何的角度去分析和思考微分方程的理论问题找到了入口.从实用上说,我们可以通过作出方向场来画出积分曲线的大概图形.这在无法(或无必要)求出解的精确表达式时,使我们能从微分方程本身的特有性质去推断出它的解的某些属性,从而使所讨论的问题在一定程度上获得解决.例2 证明:与微分方程3224xy y y x =-' (1.8) 的积分曲线关于坐标原点(0, 0)成中心对称的曲线,也是方程(1.8)的积分曲线.证 设)(x y ϕ=)(b x a <<是方程(1.8)的一条积分曲线,以x -代x ,y -代y ,得)(x y ϕ=关于原点成中心对称的曲线)(x y -=-ϕ,即)(x y --=ϕ.由于)(x y ϕ=满足方程(1.8),故有)()]([)]([43222x x x x x ϕϕϕ≡-',)(b x a <<.上式中以x -代x ,得3222)]()[()]([)]([)(4x x x x x --≡---'-ϕϕϕ,)(a x b -<<-.或将它改写为3222)]([)]([])([4x x x x x --≡---'--ϕϕϕ,)(a x b -<<-.可见)(x y --=ϕ亦满足方程(1.8).所以它也是方程(1.8)的一条积分曲线.§ 7.2 一阶微分方程的初等解法 本节讲述一阶微分方程的初等解法,即把微分方程的求解问题化为积分问题,因此也称初等积分法.虽然能用初等积分法求解的方程属特殊类型,但它们却经常出现在实际应用中,同时掌握这些方法与技巧,也为今后研究新问题时提供参考和借鉴.1. 变量分离方程形如)()(y g x f dxdy = (2.1) 的方程称为变量分离方程,其中)(x f 和)(y g 都是连续函数.当0)(≠y g 时,把(2.1)改写为dx x f y g dy )()(=.(称为分离变量), 两边积分,得通解(隐式通解)⎰⎰+=c dx x f y g dy )()(. (2.2)这里我们把积分常数C 明确写出来,而把⎰)(y g dy ,⎰dx x f )(分别理解为)(1y g 和)(x f 的一个确定的原函数.在微分方程课程中,我们总是作这样的理解.若存在0y ,使0)(0=y g ,则直接验证可知0y y =也是方程(2.1)的解(称为常数解).一般而论,这种解会在分离变量时丢失,且可能不含于通解(2.2)中,应注意补上这些可能丢失的解.例1 求方程0)(=+y x p dxdy (2.3) 的通解,其中)(x p 为连续函数.解 分离变量dx x p ydy )(-=, 两边积分得 ⎰+-=cdx x p y ~)( ln . 或 ⎰±=-dx x p e cy )(~ . 令cc ~±=,则 ⎰=-dx x p ce y )( )0(≠c .此外0=y 是方程的常数解.若允许0=c ,则此解也含于上式中.所以方程(2.3)的通解为⎰=-dx x p ce y )(. (2.4)其中c 为任意常数.例2 解方程y y x y '=-2231.解 分离变量2231xdx y ydy =-, 两边积分得方程的通解c x y +-=--3112, 或03112=+--c x y . 此外由012=-y 找到原方程的两个特解 1±=y ,但它们不能并入通解.2. 可化为变量分离方程的特殊类型1) 形如)(xy dx dy ϕ= (2.7)的方程称为齐次方程,其中ϕ是连续函数.通过变量代换,可将(2.7)化为变量分离方程,然后按变量分离方程求解. 令u x y=,或ux y =,则 u dx dux dx dy+=,代入(2.7)得)(u u dx dux ϕ=+,或x uu dx du -=)(ϕ.这是一个变量分离方程.例3 解方程x yx ydx dy +=2.解 令ux y =,代入方程得u u u dx dux +=+2,或u dx dux 2=.(2.8) 分离变量并积分,得(2.8)的通解c x u += ln .此外0=u 也是(2.8)的解.代回原变量得原方程的通解c x x y+= ln及特解0=y )0(≠x .例4 解方程2)(1y xy xy ++-='.解 令v y x=,或vy x =,则v dy dvy dy dx +=.代入方程得211v v v dy dvy ++-=+,即21v dydv y +=. 分离变量并积分,有ln ln )1ln(12y c v v =+++ )0(1>c . 从而推出)1(2v v c y ++= )0(1≠±=c c ,或222c cvy y =-, 代回原变量得)2(22c x c y +=, 其中0≠c 为任意常数.2) 形如)(222111c y b x a c y b x a f dx dy ++++= (2.9) 的方程可化为齐次方程或变量分离方程,其中f 是连续函数,i a ,i b ,i c )2 ,1(=i 都是常数,且02221≠+a a ,02221≠+b b ,02221≠+c c .分两种情形讨论:1) 0 2211=b a b a .若02≠a ,则02≠b ,因为如果02=b ,由于0 122211=-=b a b a b a ,推出01=b .与假设21 ,b b 不同时为零相矛盾,从而有k b b a a ==2121 (常数). 令u y b x a =+22,得 )(2122c u c ku f b a dx du +++=,这是变量分离方程.若02=a ,则01≠a ,由0 212211==b a b a b a ,推出02=b .从而01≠b .令u y b x a =+11,得)(2111c c u f b a dx du ++=. 亦化为变量分离方程.2) 0 2211≠b a b a . 这时方程组⎩⎨⎧=++=++00222111c y b x a c y b x a 有唯一解 α=x ,β=y .作平移代换 α-=x X ,β-=y Y ,代入方程(2.9),得)(2211Yb X a Y b X a f dX dY ++=, 这是齐次方程.例5 解方程25--+-=y x y x dx dy . 解 令u y x =-,则251-+-=u u dx du , 即27--=u dx du . 分离变量并积分,得27222c x u u +-=-, 或c x u u =+-1442.代回原变量得通解c y x xy y x =++-+410222.例6 解方程123132-+++=y x y x dx dy . 解 联立⎩⎨⎧=-+=++01230132y x y x .解得 1=x ,1-=y .令1-=x X ,1+=y Y ,得齐次方程YX YX dX dY 2332++=. 又令u XY=,得 uudX du X 2332++=. 或uu dX du X 23)1(22+-=. (2.10) 分离变量013222=-++du u u dX X 由此积分得2ln 11ln 23 1 ln ln 122c u u u X =+-+-+ )0(1>c 从而推出c u u u X =+--3224)11()1( )0(1≠±=c c 或)1()1(54+=-u c u X此外由012=-u 找到(2.10)的两个特解1±=u ,其中1=u 可并入上式(取0=c ).代回原变量.得原方程的通解)()2(5x y c x y +=+-,其中c 为任意常数.而由1-=u 代回原变量找到原方程的一个特解 x y -=.3. 一阶线性方程一阶线性方程的一般形式为)()(x Q y x P dxdy=+, (2.11) 其中P 、Q 为连续函数.当0)(≡x Q 时,(2.11)成为0)(=+y x P dxdy. (2.3) 称它为齐次线性方程,当0)(≡x Q 时,(2.11)称为非齐次线性方程.例1已求出方程(2.3)的通解⎰=-dxx p ce y )(. (2.4)现在对方程(2.11)作变换⎰=-dx x p ue y )(, (2.12)代入(2.11)化简得⎰=dx x p e x Q dxdu)()(, 由此积分,有⎰+⎰=c dx e x Q u dxx p )()(,将它代回到(2.12)即得方程(2.11)的通解))(()()(⎰+⎰⎰=-c dx e x Q e y dxx p dx x p . (2.13)上述求解方法通常称为常数变易法[把(2.4)中c 变易为x 的函数)(x u u =],公式(2.13)也称为方程(2.11)的常数变易公式.具体求解可按上述常数变易法的过程进行,也可直接代公式(2.13).例7 解方程 22yx ydx dy -=. 解 将方程改写为y x ydx dy -=2. 这是以x 为未知量的一阶线性方程.通解为)(22⎰+⎰-⎰=-c dy ye e x dydy .) ln (2y c y -=.4. 可化为一阶线性方程的特殊类型 1)贝努里(Bernoulli)方程 形如n y x Q y x P dxdy)()(=+ (2.14) 的方程称为贝努里方程,其中Q P ,为连续函数,1 ,0≠n 为常数.方程两边同乘以n y -,得 )()(1x Q y x P dxdyy n n=+--, 或)()1()()1(11x Q n y x P n dxdy n n-=-+--. 这里以n y -1为未知量的一阶线性方程.此外,当0>n 时,0=y 也是方程(2.14)的解. 例8 解方程33y x xy dxdy=+. 解 方程两边同乘以3-y ,得 323x xy dxdyy =+--, 或32222x xy dxdy -=---. 所以)2(2232c dx ex e yx x +-=⎰--)(2222c e e x e x x x ++=--212x ce x -++=.方程的通解为1)1(222=++x ce x y .2)形如)()()()()(y f x Q y f x P dxdyy f n =+' 的方程可化为线性方程或贝努里方程,其中Q P ,是连续函数,f 是可微函数,1 ≠n 为常数.例9 解方程y x y dxdyy 2sin cos sin cos =-. 解 原方程即y x y dxyd 2sin cos sin sin =-, 上式两边同乘以y 2sin -,得x y dxdyycos sin sin 12=---, 或x y dxyd cos sin sin 11-=+--. 所以方程的通解为) )(cos (sin 1c dx e x e y x x +-=⎰--x ce x x -++-=)sin (cos 21.此外,0sin =y ,即πk y =∈k (Z )也都是方程的解.上述几种特殊的方程类型,产生于微分方程发展的早期.从中我们可以体会到求解微分方程的一种方法:对于所给微分方程,总是设法通过变形或适当的变量代换将它转化为变量分离方程或一阶线性方程(当作两种基本类型)来求解,以此扩充可求解方程的范围.这种方法通常也称为分离变量法或变量代换法.例10 求解下列微分方程:(1) xey y dx dy -=2; (2) )(2222x y y x y dx dy x -=-; (3))cos(x y dxdy-=; (4) x x x e ye y e y 2212-=-+'-. 解 (1) 将方程改写为xe yy dy dx 211-=. 上式两边同乘以x e -,得211ye y dy dx e x x-=--, 或211ye y dy de x x =+--. 所以)1(112c dy e ye edy dy xy y+⎰⎰=⎰--) (l n 1c y y+=. 方程的通解为x e c y y ) (ln +=.此外,0=y 是方程的一个特解. (2) 原方程即332)21(xy y x xdx dy +-= . 上式两边同乘以3-y ,得 x y x xdx dy y 2)21(233+-=--, 或x y xx dx dy 4)24(232--=--. 所以)4(44222⎰+⋅-=--c dx e x x xe y x x )(442c e xe x x+=-方程的通解为4222x e cy y x =-. 此外,0=y 是方程的一个特解. (3) 令u x y =-,将方程化为1cos -=u dxdu. 分离变量并积分,得c x u+=2c o t .代回原变量得方程的通解c x xy +=-2cot. 此外方程有常数解πk x y 2+= ∈k (Z ).(4) 原方程可改写为2)()(x x x e y e e y --=-. 易知x e y =是它的一个特解.令x e y z -=,得 2z e z x -='. 分离变量并积分,得 c e zx +=1, 或c e z x+=1. 所以原方程的通解为c e e y x x ++=1. 此外方程有特解x e y =. 5. 恰当方程一阶微分方程也常以微分形式出现,即写成0),(),(=+dy y x Q dx y x P . (2.15)如果存在二元可微函数),(y x u ,使得dy y x Q dx y x P y x du ),(),(),(+=, 则称方程(2.15)为恰当方程.这时(2.15)成为0),(=y x du .从而c y x u =),(. 就是它的通解.这里需要用到二元函数微分学的知识,但我们可以利用一元可微函数关于微分形式不变性来理解它.例如在微分方程032223=+dy y x dx xy中,视y 为x 函数,它的左边恰好写成)(32y x d ,所以求得通解为 c y x =32.例11 求下列微分方程的通解: (1) 0)2()(2=-++dy y x dx y x ; (2) 0)1(22=++dy x xydx ;(3) 0)1()1(cos 2=-++dy yxy dx y x .解 (1) 分项组合,得0)(22=++-xdy ydx ydy dx x , 或0)3(23=+-xy y x d 所以方程的通解为c xy y x =+-233. (2) 分项组合,得0)2(2=++dy dy x xydx ,或0)(2=+y y x d .所以通解为c y y x =+2. (3) 分项组合,得 0)1c o s 2=-++y x d y y d x dy y xdx , 或0) ln (sin =++yxy x d .所以通解为c yxy x =++ ln sin . 如果方程(2.15)不是恰当的,我们还可以通过对方程本身变形,结合运用微分运算的法则和技巧,有时还辅以适当的变量代换,将它转化为恰当方程来求解.为此需要掌握一些常见的微分表达式,如)(xy d xdy ydx =+, )2(22y x d y d y x d x+=+,)(2x y d x y d x x d y =-, )(2y xd y x d y y d x =-,) (ln x yd xy ydx xdy =-, )(a r c t a n 22x y d yx y d x x d y =+- 等. 例12 求解下列微分方程: (1) 0)(=-+dy x y ydx ;(2) 0)()(=++-dy y x dx y x ; (3) 0)()(2=-++ydx xdy y ydy xdx x ;(4) 0)()(33=+++dy y x dx y x . 解 (1) 原方程即0=+-y d y x d y y d x 上式两边同乘以21y,得 02=+-y dyyxdy ydx , 或0) ln (=+y yxd . 所以通解为c y yx=+ ln . (2) 分项组合0=-++y d x x d y y d y x d x, 或0)(2122=-++y d x x d y y x d . 以221y x +乘上式,得0)(2)(222222=+-+++y x y d xx d y y x y x d , 由此积分,得c xy y x =++a r c t a n )l n (2122. (3) 以21x 乘方程,得 02=-++xy d xx d y y y d y x d x , 或0)()(2122=++xyyd y x d . 应用极坐标 θcos r x =,θsin r y =,上式变为0tan sin 212=+θθd r dr ,化简得0cos sin =+θθθd dr . 由此积分,得1c o s 1c r =+θ, 代回原变量,得方程的通解12222c xy x y x =+++, 或2222)1)((cx x y x =++,其中21c c =为任意常数.(4) 分项组合0)()(3333=+++y dyx dx y x ydy xdx , 运用微分运算法则,得 0)11(21)(21223322=+-+yx d y x y x d , 或222322)()()(xy y x d xy y x d +=+. 令 22y x u +=,xy v =,上式变为23v u dv du =, 或udv vdu du 2-=.移项合并du v udv )1(2-=,分离变量并积分,得cu v =-2)1(. 通解为)()1(222y x c xy +=-.历史上,有人曾专门按导数形式去求解一阶微分方程,也有人曾试图按微分形式统一处理.但经验表明:单纯采用一种形式总有其不便与困难.求解中我们应特别注意这两种形式的互相转化,灵活运用各种求解方法.例13 求解下列微分方程(1) 252233363224y y y x x xy x dx dy +-+-=; (2) 0)1()(2=++-dy x y dx y x . 解 (1) 改写成微分形式0)363()224(252233=+--+-dy y y y x dx x xy x . 分项组合0)32()36()24(223253=+--++dy y x dx xy dy y y dx x x . 从而有0)(323624=--++y x y y x x d . 所以通解为c y x y y x x =--++323624.(2) 将方程改写为)1(2x y xy dx dy +-=. 或xxy x dx dy +-+=121222. 所以))1(12()1(222c x dxx x x y ++⋅+-+=⎰))1(1()1(22⎰+++=c x xdx )11)1(()1(22c x x x x +++++= 2)1(12x c x +++=.§ 7. 3 一阶微分方程的应用应用微分方程解决实际问题的步骤是:(1) 分析问题,建立相应的微分方程,并提出定解条件; (2) 求定解问题;(3) 利用所得结果解释实际问题.对于步骤 (1) 所涉及到的基本方法有按规律列方程,微元分析法及模拟近似法,下面我们通过举例分别阐述它们的具体运用.1. 按规律列方程在数学、力学、物理、化学等学科中已有许多经过实践或实验检验的规律或定律,如牛顿冷却定理、牛顿运动定律、物质放射性的规律、电路问题中的基尔霍夫(Kirchhoff) 第二定律、曲线的切线性质等,它们都涉及到某些函数的变化率,由此所列出的关系式自然就是包含自变量、未知函数及其某些导数的微分方程式.例1 冷却问题把一个加热到50C 的物体,放到20C 的恒温环境中冷却,求物体的变化规律. 解 根据牛顿冷却定律:温度为u 的物体,在温度为0u 的周围环境中冷却的速率与温差0u u -成正比.在冷却过程中,设物体在时刻t 的温度为)(t u u =,物体冷却的速率就是其温度对时间的变化率dtdu.于是由冷却定律可得 )20(--=u k dtdu, (3.1) 这里0>k 为比例常数,上式右边出现负号,是因为随时间t 的增加,温度u 在减少,即当20>u 时,0<dtdu. 此外,)(t u u =应满足初始条件 20)0(=u . (3.2) 解初值问题(3.1)+(3.2),得所求温度的变化规律kteu -+=3020.可见,物体的冷却是按指数规律变化的(图7—1).当t 增加时,温度开始时下降较快,以后逐渐变慢而趋于环境温度.例2 物体在空气中的下落假设质量为m 的物体在空气中下落,初速度为零.空气阻力与下落速度的平方成正比,阻尼系数为0>k ,求下落速度)(t v v =的变化规律.解 不妨设重力mg (g 为重力加速度)的方向为正,则空气阻力为2kv -,由牛顿第二定律,可得2kv mg dtdvm-=. (3.3) 此外)(t v v =应满足初始条件0)0(=v .将方程(3.3)分离变量dt vmk g dv=-2, 或dt mkv kmg dv =-2.两边积分,得1ln 21ln 21c t m k vv kmg kmg k mg kmg +=-+ )0(1>c .化简得at kmg k mg ce vv 2=-+,其中01≠±=c c 为常数,mkga =. 由条件0)0(=v 可定出1=c .把1=c 代回到上式,并从中解出v ,得所求变化规律1122+-=at at e e k mg v .例3 R —L 电路设有电路如图7—2所示,其中R 、L 、E 都是常数.原来电路中没有电流,当把开关K 拨到1处后,电路中的电流逐渐增加,设0t t =时又将开关K 倒向2,则电路中电流又逐渐减少,试求电路中的电流I 随时间t 的变化规律.解 (1) K 拨到1处)0(0t t ≤≤时,电阻R , 电感L 与电源E 串联成一闭合回路,各元件上电 压降分别为RI ,dtdIL ,E -,由基尔霍夫第二 定律,得图7—20=-+E dtdIL RI , 或LE I L R dt dI =+. (3.3) 此外)(t I I =应满足初始条件0)0(=I . (3.4) 解初值问题(3.3)+(3.4),得)1(tL R e RE I --=.(2) K 倒向2处)(0t t ≥后,回路只由电阻R 与电感L 串联而成,故有0=+I LRdt dI . (3.5) 且满足初始条件00)1()(0I e RE t I t R ∆-=-=. (3.6)解初值问题(3.5)+(3.6),得)(00 t t LR eI I --=.所以)(t I I =的表达式为⎪⎩⎪⎨⎧≥<≤-=---. ,,0 ),1( 0)(000t t e I t t e R E I t t t L RLR例4 几何问题求一曲线,使其上任一点的切线与x 轴的交点 到切点的距离等于该交点到坐标原点的距离.解 设所求曲线的方程为)(x y y =,其 上任一点),(y x P 处的切线方程为图7—3)(x X y y Y -'=-,其中),(Y X 是切线上的动点(图7—3).上式 中令0=Y ,得切线的横截距为y yx '-.依题意,得 222)()(y yx y y y '-=+', 或222y x xyy -='. (3.7)此即所求曲线应满足的微分方程.将(3.7)改写为对称式0)(222=--dy y x xydx ,分项组合0222=+-dy y dy x xydx . 以21y乘上式将它写成 0)(2=+y yx d . 所以c y yx =+2. 或4)2(222c c y x =-+.这是中心在y 轴上且与x 轴相切于原点的圆族.2. 微元分析法用微积分的微元法来建立微分方程式实际上是寻求一些微元之间的关系式.在建立这些关系式时也要用到某已知规律或定律,与前述方法不同之处在于这里是对这些微元来应用规律或定律的.例5 溶液的混合问题一容器内盛有100升盐水,其中含盐10千克,今用每分钟2升的速率将净水注入容器,并不断进行搅拌,使混合液迅速达到均匀,同时混合液以同样速率流出容器,问在任一时刻t 容器内含盐量是多少?解 设在时刻t ,容器内含盐量为)(t Q Q =.经过时间dt 后,流出容器的溶液量为dt 2,从而流出的盐量近似为dt Q 2100⋅,其中100Q为混合液在时刻t 的浓度,而流入容器的盐量为0.于是得含盐量Q 的微元dt QdQ 2100⋅-=, 即Q dt dQ 501-=. (3.8) 此外)(t Q Q =满足初始条件10)0(=Q . (3.9) 解初值问题(3.8)+(3.9),得teQ 110-=.例6 流体流动问题有一盛满了水的圆锥形漏斗,高为10cm ,顶角为60 ,漏斗下面有面积为0.52cm 的孔,求水从小孔流出过程中漏斗中水面高度随时间变化的规律.解 建立坐标系如图7—4所示.设在时刻t ,水面高度为)(t h h =,经过时间dt ,水面高度为dh h +)0(<dh .则在时间间隔],[dt t t +内,相应于漏斗中的水层,其体积近似于以h r 33=为半径,高为)(dh -的圆柱体体积,故这水层的体积微元为 dh h dh r dV 223ππ-=-=.而在这段时间内从小孔流出的水近似为dt h V dV )(5.00=,这里)(0h V 表示水面高度为h 时从小孔中流出的水沿h 图7—4 轴反向的速度,根据水力学中有关定律gh h V 262.0)(0=.其中0.62为流量系数,g 为重力加速度.于是有 dh h dt gh 23231.0π-=,或dh h gdt 23293.0π-=. (3.10)此外)(t h h =应满足初始条件10)0(=h . (3.11) 解初值问题(3.10)+ (3.11),得所求的变化规律)10100(265.45h gt -=π. 3. 模拟近似法在生物、医学、经济等学科的实际问题中,所反映的现象往往是很复杂的.为了研究它们,需要在不同的假设下去模拟实际现象.这个过程当然是近似的,所建立的数学模型,例如微分方程模型,也是作了各种近似和简化.因此,模型的结果是否具有实际意义或满足实际要求,只有通过实践去检验.例7 单种群模型与人口问题动植物种群本身是离散变量,但由于突然增加或减少的只是单一个体或少数几个个体,与全体数量相比,这种增量是很微小的,所以我们可以近似地假设大规模种群随时间是连续地甚至可微地在变化,进而可以应用微分方程这一数学工具来研究.英国人马尔萨斯(Malthus)认为人口增长率(出生率与死亡率之差)为常数,即单位时间内人口增量与人口总量成正比.设在时间t ,人口总数为)(t P P =,则有马尔萨斯人口模型⎪⎩⎪⎨⎧=>=00)()0( P t P a aP dt dP是常数这个初值问题的解为)(00t t a e P P -=,它表明人口按指数曲线增长,这一理论已被实践证明是错误的.1837年,荷兰生物数学专家弗胡斯特(Verhulst)考虑了有影响增长率的竞争项的模拟,得出容易理解的下述单种群数学模型⎪⎩⎪⎨⎧=-=00)()(P t P P bP a dtdP其中正常数a 和b 称为生命系数.解此初值问题,得)(00)(000t t a t t a ebP bP a e aP P --+-=. (3.12) 一些生态学家测得a 的自然值为0.029,而b 的值取决于各国的社会经济条件.据文献记载,美国和法国曾用公式(3.12)预报过人口的变化,结果相当符合实际.根据安徽省统计局统计的数字,芜湖市人口总数在1992年底为206.12万人,假设当时的人口增长率为0.832℅,则00832.00=-bP a ,在公式(3.12)中取19920=t ,60100612.2⨯=P ,对芜湖市人口作出估算,将其中一部分结果与省统计局最近几年公布的芜湖市人口总数列表对照如下从表中看出,上述估算有一定的可信度.按照这个估计,芜湖市人口总数到2020年底将达到245万,到2050年将接近269万.而最终趋势为05.289)(lim ≈=+∞→bat P t (万).§7.4 可降阶的高阶微分方程二阶及二阶以上的微分方程统称为高阶微分方程.求解高阶方程的一种常用的方法就是设法降低方程的阶数.如果能把它降低为一阶方程,我们就有可能运用§7.2所讲的方法.本节介绍几种可降阶的方程类型.1. 不显含未知函数y 的方程0),,,,()()1()(=+n k k y y y x F )1(≥k . (4.1)令z y k =)(,得关于z 的k n -阶方程0),,,,()(='-k n z z z x F , 从而使方程(4.1)降低了k 阶.例1 解方程01)4()5(=-y xy . 解 令z y =)4(,得01=-'z xz . 所以cx ce z dx=⎰=1,即cx y =)4(.对x 积分4次,即得方程的通解54233251c x c x c x c x c y ++++=.2. 不显含自变量x 的二阶方程0),,(='''y y y F . (4.2)令 z y =',并以y 为新方程的自变量,z 为新未知函数,则 dydz z dx dy dy dz dx dz y ===''. 将z y =',dydzz y =''代入方程(4.2),即化为未知函数z 的一阶方程. 例2 解方程02='+''y y y . (4.3) 解 令z y =',则dydzzy ='',代入方程得 02=+z dydzyz , 或0=+z dydzy.(注意0=z 仍是方程的解) 解得 yc z =, 即y c y ='. 分离变量并积分得212c x c y += )2(1c c =.注意,方程(4.3)可改写为0)(=''y y , 这时我们也说(4.3)是恰当方程.进而有0)2(2=''y ,对x 积分两次得222212c x c y +=, 即212c x c y +=.如果一个高阶方程不是恰当的,我们也可以通过变形,看能否将它化为恰当方程来求解. 例3 解方程02='-''y y y . 解 以21y 乘方程,得 022='-''y y y y ,或0)(=''yy , 进而有0) (ln =''y .对x 积分两次得21 ln c x c y +=.此外0=y 也是方程的解.如果把通解表示为x c ce y 1=,其中c ,1c 为任意常数,则特解0=y 也并入其中.例4 一条长度为a 的均匀链条放置在一水平而无摩擦的桌面上,使链条在桌边悬挂下来的长为b )0(a b <<,问链条全部滑离桌面需要多长时间?解 如图7—5所示,设在时刻t ,链条在桌边悬挂下来的长)(t x x =,以ρ表示链条密度(即单位长的质量),按牛顿第二定律,可得g x dt x d a ) ( 22ρρ= (g 为重力加速度). 图7—5或x a g dtx d =22. 令dt dx v =,dx dv v dtx d =22代入上式得 x ag dx dv v=. (4.4) 由假设知 0==b x v . (4.5)解初值问题(4.4)+(4.5),得22b x ag v -=, 即 22b x ag dtdx -=. (4.6) 并且 b x =)0(. (4.7)从(4.6) ,(4.7) 解得)l n (22b b x x g a t -+=. 所求时间为)l n (22bb a a g a T -+=. §7.5 高阶常系数线性方程虽然我们已会求解一阶或高阶的几类特殊类型的方程.但是我们应该知道,即使是一阶方程,例如形式上很简单的黎卡提 (Riccati)方程22y x dxdy +=,我们也不可能用初等积分法求解,这是法国数学家刘维尔(Liouville)于1838年所证明的事实.对于高阶方程,求解的难度会更大,就是高阶线性方程,尽管已从理论上完全掌握其解的结构和性质.但具体求解却无一定规律可循.也就是说,除了一些特殊情形外,是很难求解的.本节我们以二阶为例讨论高阶常系数线性方程,它的求解问题将获得彻底解决.1. 二阶常系数齐次线性方程一般形式为021=+'+''y a y a y .(5.1) 其中1a ,2a 为实数.明显看出0=y )(+∞<<-∞x 是方程(5.1)的解,(称为零解或平凡解).根据方程(5.1)的特点及指数函数t e λ的特性,我们试求(5.1)如下形式的特解x e y λ=.其中λ是待定的(实或复)常数.将x e y λ= 代入(5.1),可得0)(212 =++a a e x λλλ.因为0 ≠x e λ,所以0212=++a a λλ. (5.2) 这样,对于二次代数方程(5.2)的每一个根λ,x e λ就是方程(5.1)的一个解.(5.2)称为方。

计算方法 第七章常微分数值解

计算方法 第七章常微分数值解

1. 整体截断误差和局部截断误差 整体截断误差:数值解 yn和精确解 y(xn) 之差
dy f (x, y) dx
(1)
y(x0 ) y0
en y( xn ) yn
整体截断误差除与 xn步计算有关外,还与 xn1,, x1的计算
有关
分析计算中的某一步,显式单步法的一般形式可写为:
yn1 yn hQ(xn , yn , h)
x0
解 : f ( x, y) y x 1,由Euler公式
yn1 yn h( yn xn 1)
代入h 0.1,有yn1 0.9 yn 0.1( xn 1), 依次算得果如下:
n0 1 2 3
4
5
xn 0 0.1 0.2 0.3 0.4
0.5
yn 1.0 1.0 1.01 1.029 1.0561 1.09049
一、Runge-Kutta法的基本思想(1)
若用p阶Taylor多项式近似函数y( xn1 )有:
yn1
y( xn1 )
y( xn )
hy'( xn )
h2 2!
y"( xn )
hp P!
y( p)( xn )
其中y'( x)
f ( x, y),
y'' ( x)
f
' x
(
x,
y
)
f
' y
(例:考察初值问题来自y( x) 30y( x) 在区间[0, 0.5]上的解。
y(0) 1
分别用欧拉显、隐式格式和改进的欧拉格式计算数值解。
节点 xi 0.0 0.1
0.2 0.3 0.4
0.5
欧拉显式 欧拉隐式

高数 第七章 微分方程常微分方程

高数 第七章 微分方程常微分方程
铁岭师范高等专科学校 理工学院 数学教研室
5.微分方程的初始条件、特解
定义5 用来确定 解n 微分方程 F (x, y, y, ,y(n) ) 的0 通解中任意
常数的条件: y xx0 y0 , y xx0 y0 ,
,y(n1)
y , xx0
( n 1) 0
其中 x0, y0, y0 , , y0(n1都) 是给定的值。上述这种条件叫做初始 条件。
(
y) x
的方程,称为齐次方程。
例如(1)2xy d y (x2 y2 ) d x 0
(2)(x2 y xy2 ) d x (x3 y3) 0
2.齐次方程的求解:
在齐次方程
理有
du
(u)
u
dx x
d d
y x
求出积分后,再以
程的通解。
(
y x
y x
) 中,令 u y ,化简并整
铁岭师范高等专科学校 理工学院 数学教研室
第一讲 微分方程的基本概念;可分离变量的微分方程方 程
授课题目(章节)§7.1 微分方程的基本概念 §7.2 可分离变量的微分方程方程 教学目的与要求: 了解微分方程的阶及微分方程的解、通解、初始条件、特
解等概念; 2.会识别变量可分离的一阶微分方程,熟练掌握可分离
用常数变易法,令:y ueP(x)dx ,其中u u(x) 为待定函数,
带入原非齐次微分方程d y P(x)y Q(x) ,可解
得:u Q(x)eP(x)d x d x C
dx
因此,非齐次微分方程 d y P(x)y
解为: d x
y
e
P(
x)
d
x
(
Q(x)e P(x)d x

第七章 常微分方程数值解法.

第七章  常微分方程数值解法.
12
7.1.1 欧拉法及其截断误差
4、欧拉公式的截断误差是O(h2),公式是1 阶的。
因为 yi+1 ? yi h f (xi , yi ) = y(xi ) + h y¢(xi )
y(x)

yy((xxi )i
1
)y

(yxi()xi
)1y
2
y
(()x(i x) h
xi
) 12 2
2、欧拉公式几何意义: ——用折线代替曲线计算解函数的近似值。

yi1 yi y0 y( x0 )
h
f
(
xi
,
yi )
,
i

0,1,2,

yi1 yi f y0 y(x0 )
(xi
,
yi
)
(
xi1

xi
)
,
i 0,1, 2,
11
7.1.1 欧拉法及其截断误差
(1) f(x,y)在R上连续, (2)在R上关于y满足Lipschitz(李普希兹)条件 即存在常数L,对R上任意点均有以下不等式成立:
|f(x,y1 )–f(x,y2 )|≤L|y1–y2|, x∈[a,b], y1,y2∈R 则上述初值问题存在唯一的连续可微的解函数 y = y(x)。拉法及其截断误差
初值问题:
y f (x, y)

y(
x0
)

y0
1、欧拉公式的构造思想:用差商代替导数
设 x0, x1, x2,鬃?, xn 等距,步长为xi+1 - xi = h, i = 0,1,鬃?,n 1
yⅱ(x) ? y(x + h)- y(x) , y(x h)? y(x) h y (x)= y(x)+ h f (x, y)

常微分方程初值问题的数值解法

常微分方程初值问题的数值解法

第七章 常微分方程初值问题的数值解法--------学习小结一、本章学习体会通过本章的学习,我了解了常微分方程初值问题的计算方法,对于解决那些很难求解出解析表达式的,甚至有解析表达式但是解不出具体的值的常微分方程非常有用。

在这一章里求解常微分方程的基本思想是将初值问题进行离散化,然后进行迭代求解。

在这里将初值问题离散化的方法有三种,分别是差商代替导数的方法、Taylor 级数法和数值积分法。

常微分方程初值问题的数值解法的分类有显示方法和隐式方法,或者可以分为单步法和多步法。

在这里单步法是指计算第n+1个y 的值时,只用到前一步的值,而多步法则是指计算第n+1个y 的值时,用到了前几步的值。

通过对本章的学习,已经能熟练掌握如何用Taylor 级数法去求解单步法中各方法的公式和截断误差,但是对线性多步法的求解理解不怎么透切,特别是计算过程较复杂的推理。

在本章的学习过程中还遇到不少问题,比如本章知识点多,公式多,在做题时容易混淆,其次对几种R-K 公式的理解不够透彻,处理一个实际问题时,不知道选取哪一种公式,通过课本里面几种方法的计算比较得知其误差并不一样,,这个还需要自己在往后的实际应用中多多实践留意并总结。

二、本章知识梳理常微分方程初值问题的数值解法一般概念步长h ,取节点0,(0,1,...,)n t t nh n M =+=,且M t T ≤,则初值问题000'(,),()y f t y t t Ty t y =≤≤⎧⎨=⎩的数值解法的一般形式是1(,,,...,,)0,(0,1,...,)n n n n k F t y y y h n M k ++==-@显示单步法7.2.1 显示单步法的一般形式1(,,),(0,1,...,1)n n n n y y h t y h n M ϕ+=+=-定理7.2.1 设增量函数(,,)n n t y h ϕ在区域00{(,,)|,||,0}D t y h t t T y h h =≤≤<∞≤≤内对变量y 满足Lipschitz 条件,即存在常数K ,使对D 内任何两点1(,,)t u h 和2(,,)t u h ,不等式1212|(,,)(,,)|||t u h t u h K u u ϕϕ-≤-成立,那么,若单步法的局部截断误差1n R +与1(1)p h p +≥同阶,即11()p n R O h ++=,则单步法的整体截断误差1n ε+与p h 同阶,即1()p n O h ε+=。

高等数学第七章常微分方程

高等数学第七章常微分方程

行业PPT模板:/hangye/ PPT素材下载:/sucai/ PPT图表下载:/tubiao/ PPT教程: /powerpoint/ Excel教程:/excel/
高等数学
第七章 常微分方程
因此y=eλ1x是原方程的解。 函数y=C1eλ1x+C2eλ2x的一阶导数和二阶导数分别为 y′=C1λ1eλ1x+C2λ2eλ2x y″=C1λ12eλ1x+C2λ22eλ2x 代入原方程,则 (C1λ12eλ1x+C2λ22eλ2x)-(λ1+λ2)(C1λ1eλ1x+C2λ2eλ2x)+λ1λ2( C1eλ1x+eλ2x)≡0 说明y=C1eλ1x+C2eλ2x也是原方程的解。
微分方程的概念 一阶微分方程 可降阶的高阶微分方程 二阶常系数齐次线性微分方程 二阶常系数非齐次线性微分方程
第一节 微分方程的概念
一、 微分方程的基本概念
例1 已知一条曲线经过点(2,1),且该曲线上任一点
P(x,y)处切线斜率为x,求该曲线的方程.
解 设所求曲线方程为y=y(x).由导数的概念及几何意义
F(x,f(x),f′(x),…,f(n)(x))≡0 则称y=f(x)为微分方程 (7-1-1) 在区间I上的解。
第一节 微分方程的概念
例2 验证函数y=eλ1x和y=C1eλ1x+C2eλ2x均为方程 y″-(λ1+λ2)y′+λ1λ2y=0的解。
解 y=eλ1x的一阶导数和二阶导数分别为 y′=λ1eλ1x, y″=λ12eλ1x 将y,y′,y″代入原方程中,则 λ12eλ1x-(λ1+λ2)λ1eλ1x+λ1λ2eλ1x≡0
dx

7-第七章-常微分方程1

7-第七章-常微分方程1
5/28
-------兔子数量 -------狐狸数量
兔、狐数量 变化相位图
6/28
例7.3 蝴蝶效应 数学家洛伦兹在一次讲演中将大气环流数据对初值的敏感性 形象地解释为“一只蝴蝶在巴西扇动翅膀,会引起德克萨斯 州一场龙卷风”。洛伦兹微分方程模型如下
dx x yz dt dy (y z ) , t [0, 80] dt dz xy y z dt
u x GMx / (x 2 y 2 )3/2 u GMy / (x 2 y 2 )3/2 v v y
初始条件
x (0) (R h ) y(0) 0
u(0) v 0 cos( / 2) v(0) v 0 sin( / 2)
14/28
function u=myfun(t,y) u(1,:)=-8/3*y(1)+y(2)*y(3); u(2,:)=-10*y(2)+10*y(3); u(3,:)=-y(1)*y(2)+28*y(2)-y(3);
8/28
9/28
例7.4 抛射曲线实验 假设阻力与速度成正比,在微分方程中 增加阻力项 x (t ) kx (t ) x (0) 0, x (0) v 0 cos y (t ) g ky (t ) y(0) 0, y (0) v 0 sin 符号计算方法
function z=fun1(t,N) z=0.015*N; ode23('fun1',[1994,2020],12) [T,N]=ode23('fun1',[1994,2020],12)
3/28
例7.2 捕食者与被捕食者问题 海岛上有狐狸和野兔,当野兔数量增多时,狐狸捕食野兔导致 狐群数量增长;大量兔子被捕食使狐群进入饥饿状态其数量下 降;狐群数量下降导致兔子被捕食机会减少,兔群数量回升。 微分方程模型如下

第七章 常微分方程

第七章 常微分方程


1 t RC
,
UC
把初始条件代入得,
E
0 E e0 , 即 E.
所求的充电电路的电压 U C o 的变化规律是
t
U C E (1 e

1 t RC
)
第七章
微 分 方 程
课题二十八 微分方程应用举例
[例 2]生物活体含有少量固定比的放射性 14 C, 其 死亡时存在的 14 C 量按与瞬时存量成比例的速率减 少,其半衰期约为 5730 年,在 1972 年初长沙马王 堆一号墓发掘时,若测得墓中木炭 14 C 含量为原来的 77.2%,试断定马王堆一号墓主人辛追的死亡时间. 解 设墓中木炭 14 C 在 t 时含量为x (t ), 由题设有 dx dx kx(k 0), 分离变量得 k dt, dt x kt 两边积分得 ln x ln C kt, 即x Ce .
第七章
微 分 方 程
课题二十八 微分方程应用举例
【授课时数】 总时数:4 学时. 【学习目标】 1、知道一阶和简单的二阶微分方程的解法; 2、会用一阶微分方程解决实际应用问题; 3、会用简单的二阶微分方程解决实际应用问题. 【重、难点】
重点:一、二阶微分方程的定义和解法,由通过复
习引出.
难点:一些实际应用问题用微分方程构建数学模型
5730 t ln 0.772 2139 (年) ln 2
这表明距1972年初马王堆一号墓发掘时,墓中人已 死亡约2139年,由1971 - 2139 =-168(年).
即墓中人辛追夫人死亡时间约为公元前168年.
第七章
微 分 方 程
课题二十八 微分方程应用举例
[例 3] 有高为1米的半球形容器, 水从它的底部小孔 流出, 小孔横截面积为1平方厘米(如图). 开始时容器内盛 满了水, 求水从小孔流出过程中容器里水面的高度h(水面 与孔口中心间的距离)随时间 t 的变化规律.

高等数学_第7章___常微分方程

高等数学_第7章___常微分方程

第7章 微分方程一、本章提要1. 基本概念微分方程,常微分方程(未知函数为一元函数),偏微分方程(未知函数为多元函数),微分方程的阶数(填空题).齐次方程 :()dy y dxx ϕ=或者()dxxdy yϕ=(计算) 一阶线性微分方程:()()y P x y Q x '+=或者()()x P y x Q y '+=通解公式()d ()d ()e d e P x x P x x y Q x x C -⎡⎤⎰⎰=+⎢⎥⎣⎦⎰ 或者用常数变异法求解.(计算或者填空) 线性相关,线性无关(选择) 可降解(不显含x 或y )的(计算)齐次常系数线性微分方程:特征根法(填空)非齐次常系数线性微分方程:特接用待定系数法. (计算) 微分方程解的结构定理(选择或填空). 换元法也是求解微分方程的重要方法之一. 二、要点解析问题1 常微分方程有通用的解法吗?对本章的学习应特别注意些什么?解析 常微分方程没有通用的求解方法.每一种方法一般只适用于某类方程.在本章 我们只学习了常微分方程的几种常用方法.因此,学习本章时应特别注意每一种求解方法所适用的微分方程的类型.当然,有时一个方程可能有几种求解方法,在求解时,要选取最简单的那种方法以提高求解效率.要特别注意:并不是每一个微分方程都能求出其解析解,大多数方程只能求其数值解.例1 求微分方程 '+=y y 0 的通解.解一 因为 0y y '+= 所对应的特征方程为10r +=,特征根1r =-,所以e xy C -=(C 为任意常数)为所求通解.解二 因为0=+'y y ,所以)0(d d ≠-=y y xy ,分离变量x y y d d -=,两边积分⎰⎰-=x yy d d ,1ln ln y x C =-+, 所以exy C -= (C 为任意常数)三、例题精解例3 求''=y y 4满足初始条件01,2x x yy =='== 的特解.解一 令'=y p ,则d d d d d d d d p p y py pxy x y''==⋅=.将其代入原方程''=y y 4得 y yp p4d d =,分离变量 y y p p d 4d =, 两边积分⎰⎰=y y p p d 4d ,22111422p y C =⋅+, 2224p y C =+,因为001,2x x yp y =='===,所以222241C =⨯+,可得C 2=0.故224p y =,即 p y =±2.这里'=-y y 2 应舍去,因为此时'y 与y 异号,不能够满足初始条件.将2y y '=分离变量便得其解y =23exC +.再由y x ==01,得30C =,于是所求解为2e xy =.上面解法中,由于及时地利用初始条件确定出了任意常数C 1的值,使得后续步骤变得简单,这种技巧经常用到.解二 因为''=y y 4,所以40y y ''-=,特征方程 240r -=, 特征根 122,2r r =-=, 于是其通解为2212e e x x y C C -=+, 由初始条件可得C 1=0 ,C 2=1 ,所求特解为 2e x y =.例4 求方程''+=y y x sin 的通解.解一 该方程为二阶常系数非齐次线性方程,其对应的齐次方程为 ''+=y y 0, 特征方程为 210r +=, 特征根12i,=i r r =-,齐次方程的通解为12cos sin Y C x C x =+,由于方程0sin e sin y y x x ''+==,i i αβ+=(其中0,1αβ==) 恰是特征单根,故设特解为(c o s s i n y x a xb x *=+,代入原方程,可得1,02a b =-= 所以1cos 2y x x *=-,于是所求通解为y C x C x x x =+-1212c o ss i n c o s .上述解法一般表述为:若二阶线性常系数非齐次微分方程 ''+'+=y py qy f x ()中的非齐次项[]()e()c o s ()s i nxnh f x P x x P xx αββ=+,那么该微分方程的特解可设为[]e()c o s ()s i n kxp mm y x P x x Q xx αββ=+,其中(), ()m m P x Q x 均为 m 次待定多项式 {}m h n =m ax ,.如果非齐次项中的αβ,使i αβ±不是特征方程的根,则设0k =;如果i αβ±是特征方程的单根,则取1k =. 例5 求解微分方程x xe y y y 42=+'-''。

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

1.
解 使用常数变易法求解.
将所给的方程改写成下列形式:
y 1 y 1 cos x, xx
则与其对应的线性齐次方程
的通解为
y 1 y 0 x
yC. x
23
设所给线性非齐次方程的通解为
y C(x) 1 . x
将 y 及 y代入该方程,得
于是,有
C( x) 1 1 cos x, xx
C( x) cos xdx sin x C.
种简单易行的方法,即看它们的比是否为常数事, 实 上 ,
y P( x) y Q( x).
则有
C( x) y1 C( x) y1 P( x)C( x) y1 Q( x), 即
C( x) y1 C( x)( y1 P( x) y1 ) Q( x),
17
因 y1 是对应的线性齐次方程的解,故 y1 P( x) y1 0, 因此有
C( x) y1 Q( x),
7
例2
验证方程 y 2 y 的通解为 y = Cx2 (C 为
x
任意常数),并求满足初始条件 y|x = 1 = 2 的特解.
解 由 y = Cx2 得 y = 2Cx,
将 y 及 y 代入原方程的左、右两边,左边有 y= 2Cx,
而右边 2 y 2Cx, x
所以函数 y = Cx2 满足原方程.
5
用来确定通解中的任意常数的附加条件一般称 为初始条件. 通常一阶微分方程的初始条件是
y |xx0 y0 , 即 y( x0 ) y0 . 二阶微分方程的初始条件是 y |xx0 y0 及 y |xx0 y0 , 即 y(x0) = y0 与 y(x0) = y0, 一个微分方程与其初始条件构成的问题,称为 初值问题. 求解某初值问题,就是求方程的特解.
6
例 1 验证函数 y = 3e – x – xe – x 是方程
的解.
y + 2y + y =0
解 求 y = 3e – x – xe – x 的导数, 得
y = - 4e – x + xe - x, y = 5e – x - xe - x,
将 y,y 及 y 代入原方程的左边,有 (5e – x - xe - x) + 2(- 4e – x + xe - x) + 3e – x – xe – x = 0, 即函数 y = 3e – x – xe – x 满足原方程,所以该函数 是所给二阶微分方程的解.
例如,下列方程都是微分方程 (其中 y, v, q 均
为未知函数). (1) y= kx, k 为常数; (2) ( y - 2xy) dx + x2 dy = 0;
3
(4)
y 1 a
1 y2 ;
(5)
d2q
dt 2
g sinq
l
0
(g, l 为常数).
微分方程中出现的未知函数最高阶导数的阶数,
解 将原方程改写为
dx 1 2 y x 1, dy y2
这是一个关于未知函数 x = x(y) 的一阶线性非齐次
方程,
其中
P(
y)
1
2 y2
y
,
它的自由项 Q(y) = 1.
26
代入一阶线性非齐次方程的通解公式,有
x
e
12 y2
y dy
C
e
12 y2
y
dy
dy
1
1
1
y2e y (C e y ) y2(1 Ce y ),
的通解
y P( x) y Q( x)
在运算过程中,我们取线性齐次方程的一个解为
y1 e P( x)dx ,
于是,一阶线性非齐次方程的通解公式,就可写成:
y e P( x)dx C Q( x)e P( x)dxdx.
上述讨论中所用的方法,是将常数 C 变为待 定函数 C(x),再通过确定 C(x) 而求得方程解的方法, 称为常数变易法.
21

P( x)dx
1 dx 2
x 2
,
x
e P ( x )dx e 2 ,
Q( x)e P( x)dxdx
1e xe
x
2 dx
e
x 2
,
2
代入通解公式,得原方程的通解为
xx
x
y (C e 2 )e 2 Ce 2 e x .
22
例 9 求解初值问题.
xy y cos x,
y()
第七章 常微分方程
第一节:常微分方程的基本概念 第二节:一阶微分方程 第三节:一阶微分方程的应用 第四节:二阶梯微分方程的应用
第一节 微分方程的基本概念
一、微分方程 二、微分方程的解
2
一、微分方程
定义 1 凡含有未知函数导数 (或微分) 的方程, 称为微分方程,有时简称为方程,未知函数是一元 函数的微分方程称做常微分方程,未 知 函 数 是 多 元 函数的微分方程称做偏微分方程. 本 教 材 仅 讨 论 常 微分方程,并简称为微分方程.
其中 y1 与 Q(x) 均为已知函数,所以可以通过积分 求得
C
(
x)
Q( x)dx y1
C
,
代入 y = C (x)y1 中,得 y Cy1 y1
Q( x) dx.
y1
容易验证,上式给出的函数满足线性非齐次方程
y P( x) y Q( x),
18
且含有一个任意常数,所以它是一阶线性非齐次方程
即所求通解为
1
x y2(1 Ce y ).
27
第四节 二阶常系数线性微分方 程
一、二阶线性微分方程解的结 构 二、二阶常系数线性微分方程的解法
28
一、二阶线性微分方程解的结构
二阶微分方程的如下形式
y + p(x)y + q(x)y = f (x)
称为二阶线性微分方程,简称二阶线性方程. f (x) 称为自由项,当 f (x) 0 时,称为二阶线性非齐次 微分方程,简称二阶线性非齐次方程. 当 f (x) 恒为 0 时,称为二阶线性齐次微分方程,简称二阶线性 齐次方程. 方程中 p(x)、 q(x) 和 f (x) 都是自变量 的已知连续函数. 这类方程的特点是:右边是已知 函数或零,左边每一项含 y 或 y 或 y,且 每 项 均 为 y 或 y 或 y 的一次项, 例如 y + xy + y = x2 就是二阶线性非齐次方程. 而 y + x(y)2 + y = x2 就不是二阶线性方程.
是可分离变量方程. 分离变量,得
两边积分,得
dy P( x)dx, y
ln y P( x)dx lnC,
所以,方程的通解公式为
y Ce P(x)dx
14
例 6 求方程 y + (sin x)y = 0 的通解. 解 所给方程是一阶线性齐次方程,且 P(x) = sin x, 则
P( x)dx sin xdx cos x,
解 设所求曲线的方程为 y = y(x),根 据 导 数 的几何意义及本题所给出的条件
y = y2,

dx dy
1 y2
,
积分得
x 1 C.
又由于已知曲线过点
(1,
y
2),代入上式,得
C
3 2
.
所以,求此曲线的方程为 x 3 1 .
2y
9
第二节 一阶微分方 程
一、可分离变量方程 二、一阶线性微分方程
称为微分方程的阶. 例如,方程 (1) - (3) 为一阶微 分方程,方程 (4) - (5) 为二阶微分方程. 通常,n 阶微分方程的一般形式为
F(x, y, y, , y(n)) = 0,
其中 x 是自变量, y 是未知函数,F(x, y, y, , y(n)) 是已知函数,而且一定含有 y(n).
4
二、微分方程的解
定义 2 任何代入微分方程后使其成为恒等式 的函数,都叫做该方程的解.若微分方程的解中含有 任意常数的个数与方程的阶数相同, 且任意常数之 间不能合并,则称此解为该方程的通解(或一般解). 当通解中的各任意常数都取特定值时所得到的解,称 为方程的特解.
例如方程 y = 2x 的解 y = x2 + C 中含有一个任意 常数且与该方程的阶数相同,因此,这个解是方程的 通解;如果求满足条件 y(0) = 0 的解,代入通解 y = x2 + C 中,得 C = 0,那么 y = x2 就是方程 y = 2x 的特解.
24
因此,原方程的通解为 y (sinx C) 1 C 1 sin x. x xx
将初始条件 y() = 1 代入,得 C = , 所以, 所求的特解,即初值问题的解为
y 1 ( sin x). x
25
例 10 求方程 y2dx + (x - 2xy - y2)dy = 0 的通解.
12
若 Q (x) 0,则方程成为
y P( x) y 0,

称为一阶线性齐次微分方程,简称线性齐次方程, 若 Q (x) 0,则称方程 ① 为一阶线性非齐次微分 方程,简称线性非齐次方程. 通常方程 ② 称为方 程 ① 所对应的线性齐次方程.
13
1.一阶线性齐次方程的解法
一阶线性齐次方程
y P(x) y 0
由通解公式即可得到方程的通解为 y Cecosx .
15
例 7 求方程 (y - 2xy) dx + x2dy = 0 满足初始
条件 y|x=1 = e 的特解.
解 将所给方程化为如下形式:
dy 1 2x
dx
x2
y 0,
这是一个线性齐次方程,

P(
x)
相关文档
最新文档