第九章 常微分方程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
引言
常微分方程的应用背景
自然科学
描述物理、化学、生物等自然 现象的变化规律。
工程领域
控制系统设计、航天器轨道计 算等。

第九章 连续时间:微分方程

第九章 连续时间:微分方程

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的恒等式,

高等数学:第九章 常微分方程3-4

高等数学:第九章 常微分方程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

第9章 常微分方程初值问题数值解法

第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 的切线(存在!),则斜率

常微分方程

常微分方程
若存在 ( 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

常微分方程数值解法

常微分方程数值解法

ρ ρ
n+1 n
≤1
三、梯形公式
由 分 径 y ( xn+1) = y ( xn) + 积 途 : xn+1

f ( x, y)dt

积分 梯形 式 且令:yn+1 = y( xn+1), yn = y( xn) 用 公 , h 则 yn+1 = yn + ( f (xn , yn) + f (xn+1 , yn+1)) 得: 2
第九章 常微分方程数值解法
§1 、引言
一 常 分 程 初 问 : 阶 微 方 的 值 题 dy dx = f (x, y) y( x0) = y0
'
a ≤ x ≤b
2 y 例 : 方 程 xy -2 y = 4 x ⇒ y = + 4 x 2 y 令 :f ( x , y ) = + 4 且 给 出 初 值 y (1 )= -3 x 就 得 到 一 阶 常 微 分 方 程 的 初 值 问 题 : 2 y dy = f (x, y) = + 4 dx x y(1) = − 3
n n n n n 2 // n n+1
~
y
n+1
= yn + hf ( xn, yn ) = y(xn) + hf
n+1
~
y
n+1
( x , y( x ))
n n
则 T = y( x ) − = h y (ξ ) x y 2 ~
// n+1 n+1
2
n
< ξ < xn+1

第9章 常微分方程初值问题数值解法

第9章 常微分方程初值问题数值解法

oa
b
a f ( x)dx (b a) f (b)
中矩形公式
b
ab
a f ( x)dx (b a) f ( 2 )
计算方法
梯形公式
bx
右矩形公式 中矩形公式 左矩形公式
§ 欧拉方法几何意义
y y y(x)
y0 y1 y2 0 x0 x1 x2
计算方法
x
§ 隐式欧拉方法
➢隐式欧拉法 /* implicit Euler method */
初 值 问 题 的 解 必 存 在 且唯 一 。
计算方法
§9.1 引言
三. 数值解法含义
所谓数值解法, 就是设法将常微分方程离散化, 建 立差分方程, 给出解在一些离散点上的近似值。
微分方程的数值解: 设方程问题的解y(x)的存在区 间是[a,b], 令a= x0< x1<…< xn =b, 其中hk=xk+1-xk, 如是等距节点h=(b-a)/n, h称为步长。
yi1 yi1 2h f ( xi , yi ) i 1, ... , n 1
计算方法
预估-校正法
三. 预估 — 校正法
/* predictor-corrector method */
方法 显式欧拉 隐式欧拉 梯形公式
中点公式
简单
稳定性最好
精度提高
精度低
精度低, 计算量大
计算量大
精度提高, 显式
在x0 x X上的数值解法。
四. 误差估计、收敛性
和稳定性
计算方法
§9.2 简单的数值方法与基本概念
一. 欧拉(Euler)格式
设 节 点 为xi a ih (i 0,1,2 , n) 方 法 一 :Taylor展 开 法

第九章 常微分方程数值解

第九章 常微分方程数值解



k 0, 1, 2,...
( ( 应用改进欧拉法,如果序列 yn0)1 , yn1)1 , 收敛,它的极限便
满足方程
y n 1 h yn f ( xn , yn ) f ( xn1 , yn1 ) 2
3.公式的截断误差
二元泰勒公式: 设 z=f(x,y) 在点 ( x 0 , y 0 ) 的某一邻域内连续且直到有n+1阶
首先希望能确定系数 1、2、p,使得到的算法格式有2阶 dy f x ( x, y) yi y( x i ) 的前提假设下,使得 f y ( x , y ) 精度,即在 dx 3 Ri y( xi 1 ) yi 1 Oh ( x , y ) f ( x , y ) f ( x , y ) ( f )
2
Q: 为获得更高的精度,应该如何进一步推广?
yi1 yi h[ 1 K1 2 K2 ... m Km] K1 f ( xi , yi ) K2 f ( xi 2 h, yi 21 hK1 ) K3 f ( xi 3 h, yi 31 hK1 32 hK2 )
最常用为四级4阶经典龙格-库塔法
4 阶龙格――库塔法
h y n 1 y n ( k1 2 k 2 2 k 3 k 4 ) 6 k1 f ( x n , y n ) 1 h k 2 f x n h, y n k1 2 2 1 h k 3 f x n h, y n k 2 2 2 k f x h, y hk n n 3 4
y( x ) y0 f (t , y(t ))dt
x0
x
是等价的,当x = x1时,

常微分方程

常微分方程

u
u
例3 R-L-C电路 电路
如图所示的R-L-C电路. 它包含电感L,电阻R,电容C及电源e(t). 设L,R,C均为常数,e(t)是时间t的已知函数.试求当开关K合上后,电 路中电流强度I与时间t之间的关系.
电路的 第二定律: 第二定律 解: 电路的Kirchhoff第二定律 在闭合回路中,所有支路上的电压的代数和为零 在闭合回路中 所有支路上的电压的代数和为零. 所有支路上的电压的代数和为零
三 线性和非线性
dy d y 1.如果方程 F(x, y, , L , n ) = 0 dx dx n dy d y 的左端为y及 , L , n 的一次有理式, dx dx 则称其为n 则称其为n阶线性方程.
如 (1) dy = 2 x
n
dx
(2) xdy − ydx = 0
是线性微分方程.
d 4x d 2x ( 4) + 5 2 + 3 x = sin t 4 dt dt
例1 镭的衰变规律:
设镭的衰变规律与该时刻现有的量成正比, 且已知t = 0时, 镭元素的量为R0克, 试确定在 任意t时该时镭元素的量.
注:镭的变化率与镭的量成正比。
解: 设t时刻时镭元素的量为R(t ),
dR(t ) 由于镭元素的衰变律就是R(t )对时间的变化律 , dt 依题目中给出镭元素的衰变律可得 :
dR = −kR, dt R(0) = R0
这里k > 0, 是由于R(t )随时间的增加而减少.
解之得 : R(t ) = R0 e − kt
即镭元素的存量是指数规律衰减的.
例2 物理冷却过程的数学模型
将某物体放置于空气中, 在时刻 t = 0 时, 测得它的温度为

数值分析教案_常微分方程初值问题数值解法

数值分析教案_常微分方程初值问题数值解法

第九章常微分方程初值问题数值解法图9-1n 作为()n x y 的近似值,得 ()n n y x hf ,)y x ,两边从n x 到1+n x 积分,得()dx x y x f x y x n nx x n n ⎰+=-+1))(,()1 矩形公式计算上式右侧积分,即()()x x x x x d x y x f dx x y x f n nn n⎰⎰++≈11,))(,()n ,得()n n n n y x hf y y ,1+=+,故欧拉法也称为矩形法。

为了达到较高精度的计算公式,对欧拉法进行改进,用梯形公式计算()()([1,2))(,(1++≈+n n n x f x y x f hdx x y x f n 的近似值,得9.2 龙格—库塔法前面讨论的欧拉法与改进的欧拉法都是一步法,即计算y 1+n 时,只用到前一步值。

龙格—库塔(Runge-Kutta)法(简称为R-K 方法)不是通过求导数的方法构造近似公式,而是通过计算不同点上的函数值,并对这些函数值作线性组合,构造近似公式,再把近似公式与解的泰勒展开式进行比较,使前面的若干项相同,从而使近似公式达到一定的阶数。

我们先分析欧拉法与预估—校正法。

对于欧拉法⎩⎨⎧=+=+),(111n n n n y x hf k k y y 每步计算f 的值一次,其截断误差为O (2h )。

对于预估—校正法()()⎪⎪⎩⎪⎪⎨⎧++==++=+121211,,2121k y h x hf k y x hf k k k y y n n n n n n 每步计算f 的值两次,其截断误差为O (3h ).下面对预估—校正法进行改进,将该公式写成更一般的形式()()bh y ah x hf k y x hf k k R k R y y n n n n n n ++==++=+,,2122111 (2.1)其中b a R R ,,,21为待定常数。

选择这些常数的原则是在)(n n x y y =的前提下,使11)(++-n n y x y )的阶尽量高。

文科经管类微积分第九章常微分方程

文科经管类微积分第九章常微分方程

的切线的斜率为 2x,求此曲线 L 的方程.
解 设曲线的方程为y y(x),则有
d y 2x. dx
(1) 微分方程
此外,函数y y(x) 应满足条件
y(x) 2,
(2)
x1
将(1)式两边关于x 积分,得
y 2xd x x2 C
(3)
初始条件 通解
将(2)代入(3),得 C 1, 故所求的曲线方程为
由(1)式,积分一次, 得
v0.4t20,
—(5)
s0.4tC1; —(3) 再积分一次, 得
s0.2t220t. —(6) 在(5)式中令v0, 得t50(s).
s0.2t2 C1tC2, —(4) 再把t50代入(6), 得
这里C1, C2都是任意常数.
s0.25022050500(m).
下页
•微分方程

ln y x3 C1

令 C eC1
说明: 在求解过程中每一步不一定 是同解变形,
减解.
因此可能增、
( C 为任意常数 )
( 此式含分离变量时丢失的解 y≡0 )
作业P172
1. (1)(2)(3)(4) 2. (1)(2)(5) 3. (1)
高等院校非数学类本科数学课程
大 学 数 学(一)
下页
❖几个基本概念
•微分方程的解 满足微分方程的函数叫做该微分方程的解.
提示: 在例1中, 微分方程y2x的解有yx2C和yx21. 在例2中, 微分方程s0.4的解有 s0.2t2 C1tC2, s0.2t2 20tC2和s0.2t220t.
下页
例2 验证函数 y 2sin x 3cos x是方程
y e P( x)dx elnC ,

高等数学基础教材课后答案

高等数学基础教材课后答案

高等数学基础教材课后答案1. 第一章:函数与极限1.1 函数的概念与性质1.2 极限的定义与性质1.3 常用极限和极限运算法则2. 第二章:导数与微分2.1 导数的定义与基本性质2.2 高阶导数与导数的计算2.3 微分的概念与运算3. 第三章:微分中值定理与导数应用3.1 罗尔定理与拉格朗日中值定理3.2 洛必达法则与泰勒公式3.3 极值与最值的判定3.4 应用题:切线与法线、曲率与弧长4. 第四章:不定积分与定积分4.1 不定积分的概念与性质4.2 基本积分表与积分方法4.4 牛顿-莱布尼茨公式与换元积分法5. 第五章:多元函数微分学5.1 多元函数的概念与性质5.2 偏导数与全微分5.3 隐函数与参数方程的求导5.4 高阶导数与泰勒展开5.5 一元函数与多元函数的导数比较6. 第六章:多元函数的极值与条件极值6.1 多元函数的极值判定与求解6.2 条件极值的求解6.3 隐函数的极值7. 第七章:重积分与曲线积分7.1 二重积分的概念与计算7.2 广义积分的概念与性质7.3 三重积分的概念与计算7.4 曲线积分的概念与计算8. 第八章:无界区域上的积分8.1 狄利克雷条件8.2 无界闭区域上的积分8.3 圆周率的计算9. 第九章:常微分方程9.1 一阶常微分方程的解法与应用9.2 高阶常微分方程的解法9.3 变量分离与恰当方程9.4 拉普拉斯变换与常系数线性微分方程10. 第十章:偏微分方程10.1 偏微分方程的基本概念10.2 分离变量方法与特征线法10.3 热传导方程与波动方程10.4 边界值问题与最值问题以上为《高等数学基础教材》课后习题答案的大致内容。

对于每个章节的习题,下面是一些示例题目及其解答作为参考:【第一章:函数与极限】习题1:已知函数f(x)=3x^2+2x-1,求f(-2)的值。

解答:将x=-2代入f(x),得到f(-2)=3*(-2)^2+2*(-2)-1=13。

习题2:证明函数f(x)=x^3+2x^2-3x+5是奇函数。

第九章微分方程

第九章微分方程
1 利用导数的几何意义建立微分方程
例 以点A(0,a)为起点,在第一象限内求一曲线, 使曲线上任一点P处所作切线与x轴交于T,且 |PT|=|OT|
2 利用物理意义建立微分方程
例 某种气体的气压p对于温度T的变化率与气压成 正比,与温度的平方成反比,求函数p(T)满足的微 分方程。
解:dp k p ,(k 0) dT T 2
1
X Y
X
令u Y , X
方程变为u X du 1 u , 分离变量法得 dX 1 u
X 2(u2 2u 1) C, 即Y 2 2XY X 2 C, 将 X x 1,Y y 2 代回, 得原方程的通解
( y 2)2 2( x 1)( y 2) ( x 1)2 C, 或 x2 2xy y2 2x 6 y C1.
dx P( y)x Q( y) dy
此时x看作y的函数,其通解为
x
e
P(
y )dy
(C
Q(
y
)e
P
(
y
)dy
dy)
e tan ydy (C
2
sin
ye
tan
ydy
dy)
cos y(C 2lncos y)
例10 有一个质量为m的质点,从液面由静止 状态开始垂直下沉,设在沉降过程中质 点所受的阻力与沉降速度v成正比,比 例系数为k (k>0),试求质点下沉速度v 及位置x与沉降时间t的关系。
dy g( x) h( y) dx
❖ 一阶线性微分方程
dy P( x) y Q( x) dx
❖ 齐次型方程 ❖ 伯努利方程
dy f ( y ) dx x
dy P( x) y Q( x) yn dx
可分离变量的微分方程 (分离变量法)
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 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
相关文档
最新文档