第九章 常微分方程3-4
合集下载
常微分方程

dy y
P(
x)dx,
ln | y | P( x)dx lnC1 ,(C1为任意常数)
齐次方程的通解为 y Ce P( x)dx (C eC1 )
17
2. 线性非齐次方程 dy P( x) y Q( x) dx
线性齐次方程是线性非齐次方程的特殊情况.
线性齐次方程的通解是 Ce P( x)dx ,
(3)检验改进模型, 观察所得的解能够在多大程度或范围上反映实际问题,
用实际问题检验该模型, 如果存在问题,则需研究, 改进模型.
27
例 冷却问题 将一个温度为50º的物体,放在20º的恒温 环境中冷却,求物体温度变化的规律.
解 冷却定律:“温度为T的物体,在温度为 T0 的环境中 冷却的速率与温差T T0成正比.” 设物体的温度T与时间 t的函数关系为 T T (t),
(t2 x)dt xdx 0 一阶 z x y 一阶
x
未知函数是一元函数的方程为 常微分方程;
未知函数是多元函数的方程为 偏微分方程.
方程中所出现的导数的最高阶数称为 微分方程的阶.
一般的n阶微分方程为 F ( x, y, y,, y(n) ) 0,
或已解出最高阶导数 y(n) f ( x, y, y,, y(n1) ).
9.4 微分方程的应用问题
例 把“大气压随高度变化而降低的速率与所在高度 处的气压成正比”所含关系表示出来.
解:第一步,设未知函数:
设大气压P和高度x之间的函数关系为 P P(x),
大气压随高度变化的速率为 dP
dx
第二步,根据条件写出方程 dP P, 为比例系数,
dx
第三步,取比例系数为正:因 dP 0, 故 0,
第九章 常微分方程
数值分析第九章常微分方程数值解法

高斯-赛德尔迭代法
松弛法
通过迭代更新函数值并逐步放松约束 条件来逼近解,适用于刚性和非刚性 问题。
利用线性组合迭代函数值来逼近解, 具有更高的收敛速度和稳定性。
03
数值解法的稳定性分析
数值解法的稳定性定义
数值解法的稳定性是指当微分方程的初值有微小的扰动时, 其数值解的近似值的变化情况。如果数值解在微小扰动下变 化较小,则称该数值方法是稳定的。
更高的精度和稳定性。
数值逼近法
泰勒级数法
将微分方程的解展开为泰勒级数,通过截断级数来逼 近解。
多项式逼近法
利用多项式来逼近微分方程的解,通过选取合适的基 函数和系数来提高逼近精度。
样条插值法
利用样条函数来逼近微分方程的解,具有更好的光滑 性和连续性。
迭代法
雅可比迭代法
通过迭代更新函数值来逼近微分方程 的解,具有简单易行的优点。
初值和边界条件的处理
根据实际问题,合理设定初值和边界 条件,以获得更准确的数值解。
收敛性和误差分析
对数值解进行收敛性和误差分析,评 估解的精度和稳定性。
数值解法的应用案例分析
人口增长模型
通过数值解法求解人口增长模型,预测未来人口数量,为政策制 定提供依据。
化学反应动力学
利用数值解法研究化学反应的动力学过程,模拟反应过程和结果。
数值分析第九章常微分方 程数值解法
• 引言 • 常微分方程数值解法的基本思想 • 数值解法的稳定性分析 • 数值解法的收敛性和误差分析 • 数值解法的实现和应用案例
01
引言
常微分方程的应用背景
自然科学
描述物理、化学、生物等自然 现象的变化规律。
工程领域
控制系统设计、航天器轨道计 算等。
松弛法
通过迭代更新函数值并逐步放松约束 条件来逼近解,适用于刚性和非刚性 问题。
利用线性组合迭代函数值来逼近解, 具有更高的收敛速度和稳定性。
03
数值解法的稳定性分析
数值解法的稳定性定义
数值解法的稳定性是指当微分方程的初值有微小的扰动时, 其数值解的近似值的变化情况。如果数值解在微小扰动下变 化较小,则称该数值方法是稳定的。
更高的精度和稳定性。
数值逼近法
泰勒级数法
将微分方程的解展开为泰勒级数,通过截断级数来逼 近解。
多项式逼近法
利用多项式来逼近微分方程的解,通过选取合适的基 函数和系数来提高逼近精度。
样条插值法
利用样条函数来逼近微分方程的解,具有更好的光滑 性和连续性。
迭代法
雅可比迭代法
通过迭代更新函数值来逼近微分方程 的解,具有简单易行的优点。
初值和边界条件的处理
根据实际问题,合理设定初值和边界 条件,以获得更准确的数值解。
收敛性和误差分析
对数值解进行收敛性和误差分析,评 估解的精度和稳定性。
数值解法的应用案例分析
人口增长模型
通过数值解法求解人口增长模型,预测未来人口数量,为政策制 定提供依据。
化学反应动力学
利用数值解法研究化学反应的动力学过程,模拟反应过程和结果。
数值分析第九章常微分方 程数值解法
• 引言 • 常微分方程数值解法的基本思想 • 数值解法的稳定性分析 • 数值解法的收敛性和误差分析 • 数值解法的实现和应用案例
01
引言
常微分方程的应用背景
自然科学
描述物理、化学、生物等自然 现象的变化规律。
工程领域
控制系统设计、航天器轨道计 算等。
第九章 连续时间:微分方程

c c 2 4d (b )
m 两根都为负,时间路径收敛,若两根相等, 1 m2 c / d 都小于零,其时间路径也为收敛的。最后,如果为共轭复 数,其实数部分c / d 小于零, P 的时间路径同样为收敛 的,尽管此时带有不断衰减的循环。
第7节 高阶线性微分方程
• 高阶微分方程
Y 2Y 0, 2 0 K K
Y 2Y 0, 2 0 L L
一次齐次(规模报酬不变)性
Y ( K , L) Y ( K , L)
人均项目表示为
y (k )
净投资:
K I K S K sY K
同除 L 可得
• 瓦尔拉斯均衡
• 马歇尔
• 瓦尔拉斯线性模型:
动态调整过程 dP (QD QS ) dt 代入可得
QD aP QD bP
0
dP (a b) P ( ) dt
为常系数的一阶线性微分方程,一特解(潜在均衡点)为
P
• 情形2:相等实根 辅助方程有相等实根m1 m2 r ,则该齐次方程的通解为
y (c1 c2 x)e
rx
• 情形3:共轭复根 辅助方程根为共轭复根m1 a bi 和 m2 a bi ,齐次方程 的通解为
y eax (c1 cos bx c2 sin bx)t
其中 为任意常数而g (b a) 当且仅当 g 0时 P(t ) P ,因 0 条件即为b a 在正常商品时,供给曲线不后仰,条件满足 • 马歇尔供求函数: Q
c
PD
a Q PS b b
a
其中 y1 不是 y 2的常数倍,则该方程的通解为 y c1 y1 c2 y2 , 其中 c1 , c2 为任意常数 • 定义 m2 pm q 0 • 辅助方程:
《常微分方程》全套课件(完整版)

捕捉到这种联系,而这种联系,用数学语言表达出来,其结 果往往形成一个微分方程.一旦求出这个方程的解,其运动规 律将一目了然.下面的例子,将会使你看到微分方程是表达自 然规律的一种最为自然的数学语言.
例1 物体下落问题 设质量为m的物体,在时间t=0时,在距
地面高度为H处以初始速度v(0) = v0垂直地面 下落,求ss此物体下落时距离与时间的关系.
有恒等式
因此,令
,则有
因此,所谓齐次方程,实际上就是方程(1.9)的右端函数 是一个关于变元x,y的零次齐次式.
如果我们把齐次方程称为第一类可化为变量分离的方程,那么我们 下面要介绍第二类这种方程.
1.3.2 第二类可化为变量可分离的方程 形如 (1.30) 的方程是第二类可化为变量可分离的方程.其中, 显然,方程(1.30)的右端函数,对于x,y并不
是方程(1.5)在区间(-1,+1)
上的解,其中C是任意常数.又方程(1.5)有两个明显
的常数解y =±1,这两个解不包含在上述解中.
3. 函数
是方程(1.6)在区间(-∞,
+∞)上的解,其中和是独立的任意常数.
4. 函数
是方程(1.7)在区间(-
∞,+∞)上的解,其中和是独立的任意常数.
这里,我们仅验证3,其余留给读者完成.事实上,
(1.13)
显然,方程(1.4)是一阶线性方程;方程(1.5)是一阶非线性方程;方程 (1.6)是二阶线性方程;方程(1.7)是二阶非线性方程.
通解与特解
微分方程的解就是满足方程的函数,可定义如下.
定义1.1 设函数 在区间I上连续,且有直
到n阶的导数.如果把
代入方程(1.11),得到在
区间I上关于x的恒等式,
例1 物体下落问题 设质量为m的物体,在时间t=0时,在距
地面高度为H处以初始速度v(0) = v0垂直地面 下落,求ss此物体下落时距离与时间的关系.
有恒等式
因此,令
,则有
因此,所谓齐次方程,实际上就是方程(1.9)的右端函数 是一个关于变元x,y的零次齐次式.
如果我们把齐次方程称为第一类可化为变量分离的方程,那么我们 下面要介绍第二类这种方程.
1.3.2 第二类可化为变量可分离的方程 形如 (1.30) 的方程是第二类可化为变量可分离的方程.其中, 显然,方程(1.30)的右端函数,对于x,y并不
是方程(1.5)在区间(-1,+1)
上的解,其中C是任意常数.又方程(1.5)有两个明显
的常数解y =±1,这两个解不包含在上述解中.
3. 函数
是方程(1.6)在区间(-∞,
+∞)上的解,其中和是独立的任意常数.
4. 函数
是方程(1.7)在区间(-
∞,+∞)上的解,其中和是独立的任意常数.
这里,我们仅验证3,其余留给读者完成.事实上,
(1.13)
显然,方程(1.4)是一阶线性方程;方程(1.5)是一阶非线性方程;方程 (1.6)是二阶线性方程;方程(1.7)是二阶非线性方程.
通解与特解
微分方程的解就是满足方程的函数,可定义如下.
定义1.1 设函数 在区间I上连续,且有直
到n阶的导数.如果把
代入方程(1.11),得到在
区间I上关于x的恒等式,
高等数学:第九章 常微分方程3-4

x0
2
6
.......
| (x) (x) |
ALn | x x0
|n
Ln h n A
n
0.
n!
n!
lim
n
|
(
x)
(
x)
|
0
当x
[
x0
h,
x0
h]时,
(
x)
(
x)。
推论: 考虑微分方程
y' f (x, y),(x, y) D.
若函数 f(x,y), fy(x,y) 在区域 D 上连续,则过 D 内任 一点 (x0,y0), 有且只有一条方程的通积分通过。
x
y1(x) y0 x0 f (x, y0 )dx,
估计 | y1(x) y0 | 后发现
x [x0 h, x0 h]
x
| y1(x) y0 | x0 f (x, y0 )dx M | x x0 | Mh b.
可见(x, y1(x)) R.
h
min(
a,
b M
),
M
max{|
上连续,且对 y 满足李氏条件,则初值问题在 区间 [x0-h,x0+h] 上有且只有一个解,其中常数
h
min( a,
b M
), M
max{|
f (x, y) |, (x, y) R}
y' f (x, y)
y(x0 )
y0
在 [x0-h,x0+h] 上有唯一解.
h
min(
a,
b M
),
M
max{|
f (x, y1) f (x, y2 ) f y (x, )(y1 y2 ), ( y1, y2 ).
常微分方程基本概念

其中c1,, cn为相互独立的任常数 .
注1:称函数y (x, c1,, cn )含有n个独立常数,是指
存在(x, c1,, cn )的某一邻域,使得行列式
c1
(, ',, (n1) )
(c1, c2 ,, cn )
'
c1
(n1)
c1
c2
cn
'
c2
'
cn 0
(n1) (n1)
为了从通解中得到合乎要求的特解,必须根据实 际问题给微分方程附加一定的条件,称为定解条件.
求满足定解条件的求解问题称为定解问题.
常见的定解条件是初始条件,n阶微分方程的初
始条件是指如下的n个条件:
当x
x 0时,
y
y0 ,
dy dx
y (1) 0
,,
d (n1) y dxn1
y (n1) 0
这里x0 , y0 , y0(1) ,, y0(n1)是给定的 n 1个常数.
定义6 在通解中给任意常数以确定的值而得到的解 称ቤተ መጻሕፍቲ ባይዱ方程的特解.
例如 y sinx, y cosx都是方程 y" y 0的特解. 可在通解y c1sinx c2cosx中分别取 c1 1, c2 0,得到 : y sinx,
c1 0, c2 1,得到 : y cosx.
3 定解条件
tx
dx dt
3
x
0;
d4x d2x (4) dt4 5 dt2 3x sin t;
都是常微分方程
2.偏微分方程 如果在一个微分方程中,自变量的个数为两个或两 个以上,称为偏微分方程.
如 (5) z z z ;
x y
注1:称函数y (x, c1,, cn )含有n个独立常数,是指
存在(x, c1,, cn )的某一邻域,使得行列式
c1
(, ',, (n1) )
(c1, c2 ,, cn )
'
c1
(n1)
c1
c2
cn
'
c2
'
cn 0
(n1) (n1)
为了从通解中得到合乎要求的特解,必须根据实 际问题给微分方程附加一定的条件,称为定解条件.
求满足定解条件的求解问题称为定解问题.
常见的定解条件是初始条件,n阶微分方程的初
始条件是指如下的n个条件:
当x
x 0时,
y
y0 ,
dy dx
y (1) 0
,,
d (n1) y dxn1
y (n1) 0
这里x0 , y0 , y0(1) ,, y0(n1)是给定的 n 1个常数.
定义6 在通解中给任意常数以确定的值而得到的解 称ቤተ መጻሕፍቲ ባይዱ方程的特解.
例如 y sinx, y cosx都是方程 y" y 0的特解. 可在通解y c1sinx c2cosx中分别取 c1 1, c2 0,得到 : y sinx,
c1 0, c2 1,得到 : y cosx.
3 定解条件
tx
dx dt
3
x
0;
d4x d2x (4) dt4 5 dt2 3x sin t;
都是常微分方程
2.偏微分方程 如果在一个微分方程中,自变量的个数为两个或两 个以上,称为偏微分方程.
如 (5) z z z ;
x y
第9章 常微分方程初值问题数值解法

2
数值分析
第9章 常微分方程初值问题数值解法
《常微分方程》中介绍的微分方程主要有:
(1)变量可分离的方程 (2)一阶线性微分方程(贝努利方程) (3)可降阶的一类高阶方程 (4)二阶常系数齐次微分方程 (5)二阶常系数非齐次微分方程 (6)全微分方程 本章主要介绍一阶常微分方程初值问题的数值解法。
进一步: 令
y n1 y n
xn 1 xn
y n 1 y( x n 1 ) , y n y( x n )
f ( x , y( x ))dx h f ( x n , y n )
宽
9
高
实际上是矩形法
数值分析
第9章 常微分方程初值问题数值解法
(3)
用Taylor多项式近似并可估计误差
解决方法:有的可化为显格式,但有的不行 18
数值分析
第9章 常微分方程初值问题数值解法
与Euler法结合,形成迭代算法 ,对n 0,2, 1,
( yn0 )1 yn hf x n , yn ( k 1) h ( yn1 yn f x n , yn f x n1 , ynk )1 2
7
数值分析
第9章 常微分方程初值问题数值解法
建立数值解法的常用方法
建立微分方程数值解法,首先要将微分方程离散 化. 一般采用以下几种方法: (1) 用差商近似导数
dy yx yx x x dx x y
n 1 n n 1 n
n
,
n
进一步: 令
yn1 y( xn1 ) , yn y( xn )
由 x0 , y0 出发取解曲线 y y x 的切线(存在!),则斜率
数值分析
第9章 常微分方程初值问题数值解法
《常微分方程》中介绍的微分方程主要有:
(1)变量可分离的方程 (2)一阶线性微分方程(贝努利方程) (3)可降阶的一类高阶方程 (4)二阶常系数齐次微分方程 (5)二阶常系数非齐次微分方程 (6)全微分方程 本章主要介绍一阶常微分方程初值问题的数值解法。
进一步: 令
y n1 y n
xn 1 xn
y n 1 y( x n 1 ) , y n y( x n )
f ( x , y( x ))dx h f ( x n , y n )
宽
9
高
实际上是矩形法
数值分析
第9章 常微分方程初值问题数值解法
(3)
用Taylor多项式近似并可估计误差
解决方法:有的可化为显格式,但有的不行 18
数值分析
第9章 常微分方程初值问题数值解法
与Euler法结合,形成迭代算法 ,对n 0,2, 1,
( yn0 )1 yn hf x n , yn ( k 1) h ( yn1 yn f x n , yn f x n1 , ynk )1 2
7
数值分析
第9章 常微分方程初值问题数值解法
建立数值解法的常用方法
建立微分方程数值解法,首先要将微分方程离散 化. 一般采用以下几种方法: (1) 用差商近似导数
dy yx yx x x dx x y
n 1 n n 1 n
n
,
n
进一步: 令
yn1 y( xn1 ) , yn y( xn )
由 x0 , y0 出发取解曲线 y y x 的切线(存在!),则斜率
常微分方程

若存在 ( x, c1 ,
, cn ) 的一个邻域,使得
, c1 ′ , c1 , c2 ′ , c2 , , cn ′ cn ( n 1) cn
≠0
( n 1) ( n 1) , , c1 c2
,
则称 y = ( x, c1 ,
, cn ) 含有n个相互独立的常数。
y 例: = c1 cos x + c2 sin x 是 y′′ + y = 0 的通解。 因为 y′ = c1 sin x + c2 cos x 而
是
在 (∞, +∞)上的解。
2
y = tan(t )
例:xdx +
x = 1+ x
'
在 (
π π
, ) 上的解。 2 2
ydy = 0 有隐式解 x 2 + y 2 = C ( C 为任意常数)。
n 阶方程的通解: 把含有 n 个相互独立的任意常数
c 称为 c1,c 2, , n 的解 y= x1,c1, ,c n) n 阶方程的通解。 (
耦合摆的动态演示
摆长减小的单摆
我们只研究这样一个方程:
θ( t ) 2 2 t 10 θ ( t ) + 2 θ( t ) =0 t 10 t 10 t
用Maple7编写的单摆模型的动态示意图
1.1.2 微分方程的基本概念
凡含有自变量、未知函数以及未知函数的导数(或微分)的方程称为 微分方程。例如:
用maple 7解双摆的运动微分方程
2 2 θ1 ( t ) = 10 θ2 ( t ) 20 θ1 ( t ) t
2
2 2 θ2 ( t ) = 20 θ1 ( t ) 20 θ2 ( t ) t
, cn ) 的一个邻域,使得
, c1 ′ , c1 , c2 ′ , c2 , , cn ′ cn ( n 1) cn
≠0
( n 1) ( n 1) , , c1 c2
,
则称 y = ( x, c1 ,
, cn ) 含有n个相互独立的常数。
y 例: = c1 cos x + c2 sin x 是 y′′ + y = 0 的通解。 因为 y′ = c1 sin x + c2 cos x 而
是
在 (∞, +∞)上的解。
2
y = tan(t )
例:xdx +
x = 1+ x
'
在 (
π π
, ) 上的解。 2 2
ydy = 0 有隐式解 x 2 + y 2 = C ( C 为任意常数)。
n 阶方程的通解: 把含有 n 个相互独立的任意常数
c 称为 c1,c 2, , n 的解 y= x1,c1, ,c n) n 阶方程的通解。 (
耦合摆的动态演示
摆长减小的单摆
我们只研究这样一个方程:
θ( t ) 2 2 t 10 θ ( t ) + 2 θ( t ) =0 t 10 t 10 t
用Maple7编写的单摆模型的动态示意图
1.1.2 微分方程的基本概念
凡含有自变量、未知函数以及未知函数的导数(或微分)的方程称为 微分方程。例如:
用maple 7解双摆的运动微分方程
2 2 θ1 ( t ) = 10 θ2 ( t ) 20 θ1 ( t ) t
2
2 2 θ2 ( t ) = 20 θ1 ( t ) 20 θ2 ( t ) t
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
x0
x
x
以此类推, 得到一串函数序列 y1 ( x) y0 f ( x, y0 ) dx,
x0
x [ x0 h, x0 h] x [ x0 h, x0 h]
y2 ( x) y0 f ( x, y1 ( x))dx,
x0
x
.......... ....... yn ( x) y0 f ( x, yn -1 ( x ))dx,
x x0
两边积分得 y ' dx f ( x, y )dx
x0 x0 x x0
x
x
y ( x) y ( x0 ) f ( x, y )dx y ( x) y0 f ( x, y )dx.
若y ( x)是积分方程y y0 f ( x, y )dx的解, 则两边
| f ( x, y1 ) f ( x, y2 ) | L | y1 y2 |, ( x, yi ) D, i 1,2.
若f y ( x, y )在区域D内连续, 则函数f ( x, y )在D 内的任意一个闭矩形R中关于y满足李氏条件.
证明: 因为f y ( x, y )在区域D内连续,
.......
n n n n | x x | L h n 0 | ( x ) ( x ) | AL A 0. n! n!
lim | ( x) ( x) | 0 当x [ x0 h, x0 h]时, ( x) ( x)。
n
推论:
( x) y0 f ( x, ( x))dx,
x0
x
( x) y0 f ( x, ( x))dx,
x0
x
| ( x) ( x) |
x
x0
[ f ( x, ( x)) f ( x, ( x))]dx
(4)证明初值问题的解是唯一的。
假设问题有解 y ( x)外, 还有另一个解 y ( x), 其中x [ x0 h, x0 h]. | ( x) ( x) |
上连续,且对 y 满足李氏条件,则初值问题在 区间 [x0-h,x0+h] 上有且只有一个解,其中常数
b h min(a, M ), M max{| f ( x, y) |, ( x, y) R}
推论:
考虑微分方程
y' f ( x, y), ( x, y) D.
若函数 f(x,y),fy(x,y) 在区域 D 上连续,则过 D 内任 一点 (x0,y0), 有且只有一条方程的通积分通过。
x
x0
[ f ( x, ( x)) f ( x, ( x))]dx
f ( x, y)满足李氏条件
x
x0
f ( x, ( x)) f ( x, ( x)) dx
x x0
L ( x) ( x) dx LA dx LA | x x0 |,
x0 x
( x), ( x)在[ x0 h, x0 h]上连续
有局部唯一解的充分条件: 定理:设初值问题中的函数 f(x,y) 在闭矩形域
R {( x, y) || x x0 | a, | y y0 | b}
上连续,且对 y 满足李氏条件,则初值问题在 区间 [x0-h,x0+h] 上有且只有一个解,其中常数
b h min(a, M ), M max{| f ( x, y) |, ( x, y) R}
x
x ( x) ( x) | L LA | x x0 | dx AL2 . x0 2 2 3 x | x x | | x x | 0 0 | ( x) ( x) | L AL2 dx AL3 . x0 2 6
0
x
0
f ( x, yn1 ( x))dx,x [ x0 h, x0 h] f ( x, ( x))dx,
可见(x)既是积分方程的解也是初值问题的解。
(4)证明初值问题的解是唯一的。
b h min(a, M ), M max{| f ( x, y) |, ( x, y) R}
A
x0 a x0
b M
D
x0 a
O
x0 a
x1
x2
x
x0 a
定理证明的步骤:
(1)作皮卡序列;
yn ( x) ( x), x [ x0 h, x0 h]; (2)证明 lim n
(3)由 yn ( x) y0 x x 得 ( x) y0 x
( y1 , y2 ).
注意
a. 对已经得到的初值
问题的解 y y ( x), x [ x0 h, x0 h], 通过延拓可以得到 更大范围的积分曲线.
O
y0 b y0 b
y ' f ( x, y ) , y ( x0 ) y 0
y
( x0 , y0 )
x0
x
y1 ( x) y0 f ( x, y0 ) dx,
x
x [ x0 h, x0 h]
代入积分方程,得 y2 ( x) y0 f ( x, y1 ( x))dx,
x0
x [ x0 h, x0 h]
估计 | y2 ( x) y0 | 后发现 | y2 ( x) y0 |
x0 x
x [ x0 h, x0 h]
并且( x, yn ( x )) R, n 1,2,....
(2)证明
lim yn ( x) ( x),
n
x [ x0 h, x0 h];
(3)由 yn ( x) y0 x f ( x, yn1 ( x))dx,
(1)作皮卡序列.
将常数函数
y y0 f ( x, y)dx,
x0
x
y y0, x [ x0 h, x0 h] 代入积分方程,得 y1 ( x) y0 f ( x, y0 )dx,
x0 x
x [ x0 h, x0 h]
估计 | y1 ( x) y0 | 后发现 | y1 ( x) y0 |
则f y ( x, y )在D内的任意一个闭矩形 R中也连续, 从而f y ( x, y )在R上有界. f ( x, y )关于y有偏导数, 则由 拉格朗日中值定理 , 对( x, y1 ), ( x, y2 ) R, f ( x, y1 ) ( x, y2 ) f y ( x, )( y1 y2 ) M | y1 y2 |,
考虑微分方程
y' f ( x, y), ( x, y) D.
若函数 f(x,y), fy(x,y) 在区域 D 上连续,则过 D 内任
一点 (x0,y0), 有且只有一条方程的通积分通过。
定理:设初值问题中的函数 f(x,y) 在闭矩形域
R {( x, y) || x x0 | a, | y y0 | b}
0
x
x [ x0 h, x0 h]
两边取极限,得
( x) y0 f ( x, ( x))dx,
x0
x
可见(x)既是积分方程的解也是初值问题的解。
(4)证明初值问题的解是唯一的。
y y0 f ( x, y)dx,
x0
x
假设问题有解 y ( x)外, 还有另一个解 y ( x), 其中x [ x0 h, x0 h]. 则y ( x)和y ( x)都满足积分方程,即
x1
x0 a
x2
x
x0 a
注意
b.由皮卡序列可找到 y=y(x) 的一串近似解。 c.若 f(x,y) 在区域 D 上连续,但不满足李氏 条件,这时过 D 内任一点,微分方程仍然 有解,但解可能不唯一。
| f ( x, y1 ) f ( x, y2 ) | L | y1 y2 |, ( x, yi ) D, i 1,2.
若函数f ( x, y )在凸区域D内关于y有偏导数,且 | f y | 在D内有界,则f ( x, y )在D内关于y满足李氏条件。
证明: 因为f ( x, y )在凸区域D内关于y有偏导数,则由
x0
x
求导得y ' f ( x, y ( x)), 且y ( x0 ) y0 f ( x, y )dx y0 .
x0
x0
满足初值问题,可见此时 y ( x) 也是初值问题的解。
x 3. 初值问题 y ' f ( x, y ) , 或 y y0 x f ( x, y)dx, 0 y ( x0 ) y 0
2. 初值问题
y ' f ( x, y ) , y ( x0 ) y 0
与积分方程
y y0 f ( x, y)dx,
x0
x
等价。
y ' f ( x, y ) , y ( x0 ) y 0
y y0 f ( x, y)dx,
x0
x
y ' f ( x, y ) 由 , y ( x0 ) y0
(4)证明初值问题的解是唯一的。
假设问题有解 y ( x)外, 还有另一个解 y ( x), 其中x [ x0 h, x0 h].
2 | x x | 0 | ( x ) ( x ) | L LA | x x0 | dx AL2 . x0 2 2 3 x | x x | | x x | 0 0 | ( x ) ( x ) | L AL2 dx AL3 . x0 2 6 x
x
x
以此类推, 得到一串函数序列 y1 ( x) y0 f ( x, y0 ) dx,
x0
x [ x0 h, x0 h] x [ x0 h, x0 h]
y2 ( x) y0 f ( x, y1 ( x))dx,
x0
x
.......... ....... yn ( x) y0 f ( x, yn -1 ( x ))dx,
x x0
两边积分得 y ' dx f ( x, y )dx
x0 x0 x x0
x
x
y ( x) y ( x0 ) f ( x, y )dx y ( x) y0 f ( x, y )dx.
若y ( x)是积分方程y y0 f ( x, y )dx的解, 则两边
| f ( x, y1 ) f ( x, y2 ) | L | y1 y2 |, ( x, yi ) D, i 1,2.
若f y ( x, y )在区域D内连续, 则函数f ( x, y )在D 内的任意一个闭矩形R中关于y满足李氏条件.
证明: 因为f y ( x, y )在区域D内连续,
.......
n n n n | x x | L h n 0 | ( x ) ( x ) | AL A 0. n! n!
lim | ( x) ( x) | 0 当x [ x0 h, x0 h]时, ( x) ( x)。
n
推论:
( x) y0 f ( x, ( x))dx,
x0
x
( x) y0 f ( x, ( x))dx,
x0
x
| ( x) ( x) |
x
x0
[ f ( x, ( x)) f ( x, ( x))]dx
(4)证明初值问题的解是唯一的。
假设问题有解 y ( x)外, 还有另一个解 y ( x), 其中x [ x0 h, x0 h]. | ( x) ( x) |
上连续,且对 y 满足李氏条件,则初值问题在 区间 [x0-h,x0+h] 上有且只有一个解,其中常数
b h min(a, M ), M max{| f ( x, y) |, ( x, y) R}
推论:
考虑微分方程
y' f ( x, y), ( x, y) D.
若函数 f(x,y),fy(x,y) 在区域 D 上连续,则过 D 内任 一点 (x0,y0), 有且只有一条方程的通积分通过。
x
x0
[ f ( x, ( x)) f ( x, ( x))]dx
f ( x, y)满足李氏条件
x
x0
f ( x, ( x)) f ( x, ( x)) dx
x x0
L ( x) ( x) dx LA dx LA | x x0 |,
x0 x
( x), ( x)在[ x0 h, x0 h]上连续
有局部唯一解的充分条件: 定理:设初值问题中的函数 f(x,y) 在闭矩形域
R {( x, y) || x x0 | a, | y y0 | b}
上连续,且对 y 满足李氏条件,则初值问题在 区间 [x0-h,x0+h] 上有且只有一个解,其中常数
b h min(a, M ), M max{| f ( x, y) |, ( x, y) R}
x
x ( x) ( x) | L LA | x x0 | dx AL2 . x0 2 2 3 x | x x | | x x | 0 0 | ( x) ( x) | L AL2 dx AL3 . x0 2 6
0
x
0
f ( x, yn1 ( x))dx,x [ x0 h, x0 h] f ( x, ( x))dx,
可见(x)既是积分方程的解也是初值问题的解。
(4)证明初值问题的解是唯一的。
b h min(a, M ), M max{| f ( x, y) |, ( x, y) R}
A
x0 a x0
b M
D
x0 a
O
x0 a
x1
x2
x
x0 a
定理证明的步骤:
(1)作皮卡序列;
yn ( x) ( x), x [ x0 h, x0 h]; (2)证明 lim n
(3)由 yn ( x) y0 x x 得 ( x) y0 x
( y1 , y2 ).
注意
a. 对已经得到的初值
问题的解 y y ( x), x [ x0 h, x0 h], 通过延拓可以得到 更大范围的积分曲线.
O
y0 b y0 b
y ' f ( x, y ) , y ( x0 ) y 0
y
( x0 , y0 )
x0
x
y1 ( x) y0 f ( x, y0 ) dx,
x
x [ x0 h, x0 h]
代入积分方程,得 y2 ( x) y0 f ( x, y1 ( x))dx,
x0
x [ x0 h, x0 h]
估计 | y2 ( x) y0 | 后发现 | y2 ( x) y0 |
x0 x
x [ x0 h, x0 h]
并且( x, yn ( x )) R, n 1,2,....
(2)证明
lim yn ( x) ( x),
n
x [ x0 h, x0 h];
(3)由 yn ( x) y0 x f ( x, yn1 ( x))dx,
(1)作皮卡序列.
将常数函数
y y0 f ( x, y)dx,
x0
x
y y0, x [ x0 h, x0 h] 代入积分方程,得 y1 ( x) y0 f ( x, y0 )dx,
x0 x
x [ x0 h, x0 h]
估计 | y1 ( x) y0 | 后发现 | y1 ( x) y0 |
则f y ( x, y )在D内的任意一个闭矩形 R中也连续, 从而f y ( x, y )在R上有界. f ( x, y )关于y有偏导数, 则由 拉格朗日中值定理 , 对( x, y1 ), ( x, y2 ) R, f ( x, y1 ) ( x, y2 ) f y ( x, )( y1 y2 ) M | y1 y2 |,
考虑微分方程
y' f ( x, y), ( x, y) D.
若函数 f(x,y), fy(x,y) 在区域 D 上连续,则过 D 内任
一点 (x0,y0), 有且只有一条方程的通积分通过。
定理:设初值问题中的函数 f(x,y) 在闭矩形域
R {( x, y) || x x0 | a, | y y0 | b}
0
x
x [ x0 h, x0 h]
两边取极限,得
( x) y0 f ( x, ( x))dx,
x0
x
可见(x)既是积分方程的解也是初值问题的解。
(4)证明初值问题的解是唯一的。
y y0 f ( x, y)dx,
x0
x
假设问题有解 y ( x)外, 还有另一个解 y ( x), 其中x [ x0 h, x0 h]. 则y ( x)和y ( x)都满足积分方程,即
x1
x0 a
x2
x
x0 a
注意
b.由皮卡序列可找到 y=y(x) 的一串近似解。 c.若 f(x,y) 在区域 D 上连续,但不满足李氏 条件,这时过 D 内任一点,微分方程仍然 有解,但解可能不唯一。
| f ( x, y1 ) f ( x, y2 ) | L | y1 y2 |, ( x, yi ) D, i 1,2.
若函数f ( x, y )在凸区域D内关于y有偏导数,且 | f y | 在D内有界,则f ( x, y )在D内关于y满足李氏条件。
证明: 因为f ( x, y )在凸区域D内关于y有偏导数,则由
x0
x
求导得y ' f ( x, y ( x)), 且y ( x0 ) y0 f ( x, y )dx y0 .
x0
x0
满足初值问题,可见此时 y ( x) 也是初值问题的解。
x 3. 初值问题 y ' f ( x, y ) , 或 y y0 x f ( x, y)dx, 0 y ( x0 ) y 0
2. 初值问题
y ' f ( x, y ) , y ( x0 ) y 0
与积分方程
y y0 f ( x, y)dx,
x0
x
等价。
y ' f ( x, y ) , y ( x0 ) y 0
y y0 f ( x, y)dx,
x0
x
y ' f ( x, y ) 由 , y ( x0 ) y0
(4)证明初值问题的解是唯一的。
假设问题有解 y ( x)外, 还有另一个解 y ( x), 其中x [ x0 h, x0 h].
2 | x x | 0 | ( x ) ( x ) | L LA | x x0 | dx AL2 . x0 2 2 3 x | x x | | x x | 0 0 | ( x ) ( x ) | L AL2 dx AL3 . x0 2 6 x