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

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,
第九章 常微分方程
30第九章 连续时间:微分方程

• 索罗—斯旺新古典增长模型 新古典生产函数 Y Y (K, L) 边际产品为正但递减
Y K
2Y 0, K 2
0
Y L
0,
2Y L2
0
一次齐次(规模报酬不变)性
Y (K,L) Y (K, L)
人均项目表示为
y (k)
净投资:
K I K S K sY K
同除 L可得
K / L sy k s(k) k
yk a
该非齐次方程的通解为 y(x) y y(0)eax
定义
• y(x) y,y 收敛于y ,y 的时间路径是稳定的
在上例中,当且仅当 a 0时,y(x) y
• 伯努利方程
dy ay cym dx
m 其中a 和 c为常数或者 x 的函数, 为任意除0和1之外的
实数,两边同除 ym 可得
形式 P(D)y 0的通解则非齐次方程 P(D) y f (x) 的通解
为 y yc yp 。
第3节 一阶常系数线性微分方程
最简单形式
dy ay f (x) dx
定理 其非齐次方程的特解为
y(x) eax x eas f (s)ds 0
特殊情形 dy ay k dx
其一个特解(潜在均衡点)为
dt
为常系数的一阶线性微分方程,一特解(潜在均衡点)为
通解为
P
ab
P(t) P cegt
其中c为任意常数而g (b a)
当且仅当 g 0时P(t) P ,因 0条件即为b a
在正常商品时,供给曲线不后仰,条件满足
• 马歇尔供求函数:
PD
a
Q a
PS
b
Q b
动态调整过程:
dQ dt
第九章 连续时间:微分方程

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 • 辅助方程:
第九章-微分方程与差分方程简介市公开课一等奖省赛课获奖课件

x
C2
例3.求解微分方程
y
y2 ,y(0) 1,y(0) 1. y
解: 设
y
p( y) ,则
y
p
dp dy
代入方程得
p dp p2 , dy y
p(
dp dy
p y
)
0
p0
27
第27页
(三)不显含自变量 x 二阶微分方程
2
第2页
第一节 微分方程普通概念
例2.设 s=s(t) 为作自由落体运动物体在 t 时刻
下落距离, 则有
d 2s dt 2 g
s(t) g
s g
ds dt
g
ds dt
gt
C1
s(0) 0
s(0)
0
ds gdt
ds gdt
s gt C1
ds ( gt C1 )dt
ds (gt C1 )dt
于价格P线性函数: QS a bP , QD c dP ,
且 a, b, c, d 都是已知正常数. 当 QS = QD 时, 得
均衡价格 P
ac .
当 QS
> QD 时, 价格将下降,
bd
当 QS < QD 时, 价格将上涨,故价格是时间t 函数.
假设在时刻t价格P(t)改变率与这时过剩需求量
x
因
P(
x)dx
1 x
dx
ln
x
ln
1 x
,
Q(
x)e
P
(
x )dx
dx
1
x 2eln x dx
xdx x2 ,
2
故 y ( x2 C )e(ln x) ( x2 C ) x Cx x3 .
微分方程罗兆富等编第九章非线性偏微分方程Adomian分解法全篇

F(u)是非线性项, g是自由项 .
学者们已证明, 无论是从算子方程Lxu还是从Lyu开始
都可得到解
u
un
并且这样得到的解都是等价的并且都
收敛于精确解. n0
然而, 在Lx 和Ly 选用哪一个来求解定解问题则依赖 于下列两个基点:
具(1体)能而使言计之算, 量我达们最考小虑;算子形式的非线性微分方程 (2)具有L使xu解 L级yu数具Ru有加F (速u)收 敛g 的附加条件. (9.2.01)
y
),
Lx
4 x4
.
(9.2.04)
(9.2.01)
14
机动 目录 上页 下页 返回 结束
un
0
Lx1g
Lx1
Ly
un
Lx1
R
un
Lx1
An
n0
n0
n0
n0
(9.2.04)
Adomian分解法指出, 通项un的递推公式是
也就是
u0 0 Lx1g,
uun
0LxL1Lx1ygun1Lx1LLyx1uR(uLnx11R)uLxL1xA1nF1(,un)
t xt2dt 0
0
u(x,t) un (x,t)
n0
uu32.((..xx.,,.tt.)).......LL.ntt.11.0.AA.u12.n..(.x.,..t00.t)t00tddtxtt0013
xt
3
x
Lt 1
(
n0
An
)
xt ■
18
机动 目录 上页 下页 返回 结束
例2. 求解非齐次偏微分方程
机动 目录 上页 下页 返回 结束
例3. 计算F(u)=uux的Adomian多项式.
学者们已证明, 无论是从算子方程Lxu还是从Lyu开始
都可得到解
u
un
并且这样得到的解都是等价的并且都
收敛于精确解. n0
然而, 在Lx 和Ly 选用哪一个来求解定解问题则依赖 于下列两个基点:
具(1体)能而使言计之算, 量我达们最考小虑;算子形式的非线性微分方程 (2)具有L使xu解 L级yu数具Ru有加F (速u)收 敛g 的附加条件. (9.2.01)
y
),
Lx
4 x4
.
(9.2.04)
(9.2.01)
14
机动 目录 上页 下页 返回 结束
un
0
Lx1g
Lx1
Ly
un
Lx1
R
un
Lx1
An
n0
n0
n0
n0
(9.2.04)
Adomian分解法指出, 通项un的递推公式是
也就是
u0 0 Lx1g,
uun
0LxL1Lx1ygun1Lx1LLyx1uR(uLnx11R)uLxL1xA1nF1(,un)
t xt2dt 0
0
u(x,t) un (x,t)
n0
uu32.((..xx.,,.tt.)).......LL.ntt.11.0.AA.u12.n..(.x.,..t00.t)t00tddtxtt0013
xt
3
x
Lt 1
(
n0
An
)
xt ■
18
机动 目录 上页 下页 返回 结束
例2. 求解非齐次偏微分方程
机动 目录 上页 下页 返回 结束
例3. 计算F(u)=uux的Adomian多项式.
第九章--微分方程与差分方程简介

19
于是非齐次方程的一个特解为:y* =kxa x-1 x
例5 求解差分方程 2y x+1 − 4y x = 2
解:原方程可化为 y x+1 − 2y x = 2 x % 则相应齐方程的通解为 y x =C ⋅ 2 x 由于p=2=a, 所以原方程的特解应设为 y* = Ax 2 x x 代入原方程得: A(x+1)2 x +1 − 2 Ax 2 x = 2 x , 1 ⇒A= 2 1 x * y x = x 2 =x 2 x -1 于是 2 所以原方程的通解为: y x =x 2 x -1 +C ⋅ 2 x
(2)∆(cyx ) = c∆y x (c为常数)
(3)∆ (ay x + bz x ) = a∆y x + b∆z x , b为常数) (a
(4)∆ ( yx z x ) = yx +1∆z x + z x ∆yx = y∆z x + z x +1∆yx
yx z x ⋅ ∆y x − y x ⋅ ∆z x (5) ∆( ) = zx z x ⋅ z x +1
23
1、二阶齐次差分方程的通解 由9.6节可知,要求齐次差分方程的通解,只需找出 两个线性无关的特解即可。仿照一阶齐次差分方程, 设二阶齐次差分方程存在指数形式的解: y x = λ x , (λ ≠ 0) 代入原方程得:
λ x+2 + pλ x+1 + qλ x = 0
即:
λ x + pλ + q = 0
11
9.6、常系数线性差分方程 、
9.6.1 n阶 系 线 差 方 的 本 质 常 数 性 分 程 基 性 n阶 系 线 差 方 的 般 式 : 常 数 性 分 程 一 形 为 yx+n +p1yx+n-1+L+pn-1yx+1+pny1 = f (x) 其 , 1,, n为 知 数 且 n ≠ 0, (x)为 知 数 中 pL p 已 常 , p f 已 函 。 当 (x)=0时 上 方 则 n阶 系 齐 线 差 方 。 , 述 程 为 常 数 次 性 分 程 f 当 (x) ≠ 0时 上 方 则 n阶 系 非 次 性 分 程 , 述 程 为 常 数 齐 线 差 方 。 f
于是非齐次方程的一个特解为:y* =kxa x-1 x
例5 求解差分方程 2y x+1 − 4y x = 2
解:原方程可化为 y x+1 − 2y x = 2 x % 则相应齐方程的通解为 y x =C ⋅ 2 x 由于p=2=a, 所以原方程的特解应设为 y* = Ax 2 x x 代入原方程得: A(x+1)2 x +1 − 2 Ax 2 x = 2 x , 1 ⇒A= 2 1 x * y x = x 2 =x 2 x -1 于是 2 所以原方程的通解为: y x =x 2 x -1 +C ⋅ 2 x
(2)∆(cyx ) = c∆y x (c为常数)
(3)∆ (ay x + bz x ) = a∆y x + b∆z x , b为常数) (a
(4)∆ ( yx z x ) = yx +1∆z x + z x ∆yx = y∆z x + z x +1∆yx
yx z x ⋅ ∆y x − y x ⋅ ∆z x (5) ∆( ) = zx z x ⋅ z x +1
23
1、二阶齐次差分方程的通解 由9.6节可知,要求齐次差分方程的通解,只需找出 两个线性无关的特解即可。仿照一阶齐次差分方程, 设二阶齐次差分方程存在指数形式的解: y x = λ x , (λ ≠ 0) 代入原方程得:
λ x+2 + pλ x+1 + qλ x = 0
即:
λ x + pλ + q = 0
11
9.6、常系数线性差分方程 、
9.6.1 n阶 系 线 差 方 的 本 质 常 数 性 分 程 基 性 n阶 系 线 差 方 的 般 式 : 常 数 性 分 程 一 形 为 yx+n +p1yx+n-1+L+pn-1yx+1+pny1 = f (x) 其 , 1,, n为 知 数 且 n ≠ 0, (x)为 知 数 中 pL p 已 常 , p f 已 函 。 当 (x)=0时 上 方 则 n阶 系 齐 线 差 方 。 , 述 程 为 常 数 次 性 分 程 f 当 (x) ≠ 0时 上 方 则 n阶 系 非 次 性 分 程 , 述 程 为 常 数 齐 线 差 方 。 f
第九章 微分方程 4

而 y C2e 2 x 2(C1 C2 x )e 2 x .
4
2x 2x 而 y C2e 2(C1 C2 x )e .
y(0) C1 1 C1 1 由初始条件得: 解得 : C 2 2, y(0) 2C1 C 2 0,
d x 由牛顿第二定律得: 100 2 gx , dt 2 d x g 即 x 0, 2 100 dt
2
O
且满足: x(0) 20, x(0) 0,
x
g t 100
求得通解: x C1e
g t 100
C 2e
,
x
15
C1 10 x (0) 20 由初始条件 , 得: , x ( 0 ) 0 C 2 10
所求的特解为:
y (1 2 x )e 2 x .
( 3) y 2 y 2 y 0.
解: 特征方程: 2 2 0,
2
特征根为: 1, 2 1 i ,
方程通解: y e (C1 cos x C2 sin x ).
5
x
B. 二阶线性常系数非齐次微分方程的通解
bm ,b m 1 ,...,b1 , b0 . 比较两端同次幂系数 , 即可求出:
( ii ) 如果是特征方程: p q 0 的单根
2
则 2 p q 0, 2 p 0, Q( x )是m次的多项式 ,
可设 : Q( x ) x(bm x m bm 1 x m 1 ... b1 x b0 ),
2 2
3A 9 A 3 14 2 4 A 3B 0 B 4 故 y p ( x ) 3 x 4 x 3 , 2 A 2 B 3C 0, C 14 3 ,
第9章微分方程初值问题的数值解法-1

(x k x k 1 )
y ( x k 1 ) y ( x k ) h y ( ) y ( x k ) h f ( , y ( ) )
记 K*f(,y()) 称为[xk , xk+1]上的平均斜率. 故
y(xk1)y(xk)hK*
当
y(i) k
y(i)(xk)
时,
有
y(xk1)yk1O (hp1). 此时①为
p 阶Taylor方法. p=1时即为Euler公式.
例2: 取步长 h = 0.1, 用一阶、二阶和四阶Taylor方法求解下列初 值问题
y y2
,
y(0) 1
0x1. 2
解: (1) 一阶Taylor法
yk1yk 0.1yk2
Taylor公式推导:
y(xk1)y(xk)hy(xk)h 2 2y(k), xkkxk1
yk1ykhf(xk,yk) k0,1,L,n1
Euler公式几何意义:
y
P2 P1 P0
Pk
也称折线法
x
2. 梯形法
若采用梯形公式计算(★)中的积分项,则有
y(xk1)y(xk)h 2[f(xk,y(xk))f(xk1,y(xk1))]
y ( x k 1 ) y ( x k ) h y ( x k ) h 2 2 !y ( x k ) L h p p !y (p )( x k ) O ( h p 1 )
令
yk 1ykhyk h 22 !yk Lh p p !yk (p)
①
称之为Taylor级数法. 其中 y k (i)y(i)(x k),i 0 ,1 ,2 ,L,p
y(2y3)6y2y6y4
y(4) 24y3y24y5
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
建立微分方程举例
1 利用导数的几何意义建立微分方程
例 以点A(0,a)为起点,在第一象限内求一曲线, 使曲线上任一点P处所作切线与x轴交于T,且 |PT|=|OT|
2 利用物理意义建立微分方程
例 某种气体的气压p对于温度T的变化率与气压成 正比,与温度的平方成反比,求函数p(T)满足的微 分方程。
假设非齐次方程的解为 y C ( x )e
P ( x ) dx C ( x ) Q( x )e dx C
dy P ( x ) y Q( x ), (Q( x )不恒为零) dx
P ( x ) dx
即得通解 y e
Ce
P ( x ) dx
P ( x ) dx
N N 0e
( n m ) t
1 建立共焦抛物线族 y 4C ( x C ) (其中C为任意常 数)所满足的微分方程。
2
2 ydx ( x 2 4 x )dy 0
一阶线性微分方程
定义:形如
y' P ( x ) y Q( x ) 的方程,
称为一阶线性微分方程。 Q( x )称为自由项
f(x,y)可表示成一个
x 的函数与一个 y 的 函数的乘积,
f ( x, y ) g( x )h( y )
dy g( x )dx h( y ) dy h( y ) g( x )dx C
例1
2 ( 2 y cos y ) dy 6 x dx 解:分离变量
( 2 y cos y )dy 6 x dx
y
x
例3
解:设曲线 y=y(x)与椭圆族中的任一椭圆的 交点为M(x,y),则曲线 y=y(x) 在交点 M 处的切线斜率为k1=y’,椭圆在该点处
x 的切线斜率为 k 2 2y
'
2y y 由k1k2=-1.可得 x ' 2y y 即初值问题 x y(a ) b
u
y x
du 2 u u x dx u 1
(1 u)du dx 整理得 2 2 2u u x
d ( 2 2u u ) 2dx 凑微分得 2 2 2u u x
2
两边同时积分得 ln(2 2u u ) 2 ln x ln C
2
( 2 2u u2 ) x 2 C
2
y sin y 2 x C
2 3
一般得到的 y 是 x 的一个隐函数
例2
解:分离变量
y x dy dx 2 2 1 y 1 x
2 2
1 d (1 y ) 1 d (1 x ) 2 2 1 y 2 1 x2
ln(1 y ) ln(1 x ) ln C
dp p 解: k 2 , ( k 0) dT T
3 利用微元法建立微分方程
例 某个地区人口总数 N 是时间 t 的函数,N=N(t).若 这个地区人口的出生率为 n (此时单位时间出生数为 nN),死亡率为m(此时单位时间死亡数为mN).现考察 任一时刻的人口总数. 微元法: [t,t+dt]时间段内, 人口增量=这段时间内出生的人数-死亡的人数 dN=nNdt-mNdt
即2 x 2 2 xy y 2 C
例13
dy y2 2 dx xy x
y u x
y 2 ( ) dy 解: x dx y 1 x
du u2 u x dx u 1
积分得u ln u ln C ln x
1 1 (1 )du dx u x
通解为y Ce
dx x tan y 2 sin y dy
dx P ( y ) x Q( y ) dy
此时x看作y的函数,其通解为
xe
P ( y ) dy
P ( y ) dy (C Q( y )e dy )
tan ydy (C 2 sin ye dy)
e
tan ydy
2 d s dy 2 2 0.4 3x dx dt ' y ( 1 ) 2 s (0) 20 s( 0) 0 称不含任意常数的解为微分方程的 特解。
。
一般地,n阶微分方程 其通解为
F ( x , y, y' , y' ' , , y
1 ln x 解:P ( x ) , Q( x ) 2 x x 1 dx 1 dx ln x x x 则通解为 y e (C 2 e dx ) x ln x ln 1 x e (C 2 xdx) x
1 1 2 (C ln x ) x 2
例9
cos ydx ( x 2 cos y ) sin ydy 0
dy 2dx y x
dy 2dx y x b a
y x ln 2 ln b a b 2 y 2 x a
y
x
例5
dN (n m ) N dt N ( 0) N 0
dN ( n m )dt N N0 0
N
t
N ln ( n m )t N0
dN (n m ) N dt
例 一个容器中装有体积V0 m3的溶液,溶液中含有 某种溶质x0,现以Q m3/s流量向容器中注入清水 (设容器中装有搅拌器使溶解均匀),并以同样流 量从容器排出溶液,求溶液中溶质含量x随时间变化 的规律x(t)。
一阶微分方程
可分离变量微分方程 一阶线性微分方程 齐次型方程 伯努利方程
齐次方程的通解 y Ce
常数变易法,设
x
Байду номын сангаас
2
y C ( x) x
2
代入非齐次方程中得 C ' ( x ) 1
C ( x) x C
得原方程的通解为 y ( x C )x2
例8
x y' xy ln x 0
2
dy 1 ln x y 2 dx x x
1 y(1) 2
P ( x ) dx [ Q( x )e dx c]
P ( x ) dx Q( x )e dx
e
P ( x ) dx
对应齐次 方程通解
非齐次方程特解
例7
dy 2 y x2 dx x
xy'2 y x
2 P( x) x
2 dx x
3
dy 2 解: 齐次 y0 dx x
dv m mg kv dt
v ( 0) 0
kv
.M
mg
ve
k dt m
e
k t m
k t mg m Ce (C ge dt ) k
(C ge
k t m
k dt m
dt )
dx mg Ce v(t ) k dt
x ( 0) 0
dy g( x ) h( y ) dx
dy P ( x ) y Q( x ) dx
dy y f( ) dx x
dy P ( x ) y Q( x ) y n dx
可分离变量的微分方程 (分离变量法)
dy 一阶微分方程 f ( x , y ), dx
dy g( x ) h( y ) dx
当Q( x ) 0时, 称为一阶线性齐次方程,
否则称为一阶线性非齐次方程。
dy P( x) y 0 dx
dy P ( x )dx y
一阶线性齐次微分 方程的通解
ln y P ( x )dx ln C
dy P ( x ) dx y
y Ce
P ( x ) dx
可分离变量的方程
当 f ( u) u 0时, 得
du dx , f ( u) u x
当 u0 ,使 f ( u0 ) u0 0, 则 u u0是新方程的解,
代回原方程 , 得齐次方程的解 y u0 x .
dy 2 x y 例12 dx x y y 2 dy x 解: y dx 1 x
( n)
)0
y y( x, c1 , c2 , , cn )
( n 1 )
方程的初始条件为:
y( x0 ) y0 , y' ( x0 ) y1 , , y
( x 0 ) y n 1
初值 ( n) F ( x , y , y ' , y ' ' , , y )0 问题 (柯 y( x0 ) y0 , y' ( x0 ) y1 , , y ( n1) ( x0 ) yn1 西问 题)
动后,列车行驶了多少时间才停住?且列车
行驶了多少距离?
解:设列车制动后 t秒内行驶了 s s(t )m, 则
d 2s 2 0.4 dt ' s (0) 20 s( 0) 0
定义1
凡表示自变量、未知函数及其未知函数导数 的方程称为微分方程
。
未知函数是一元函数的微分方程称为常微分方程。 未知函数是多元函数的微分方程称为偏微分方程。
2 2
1 y C (1 x )
2 2
把初始条件代入通解得特解 (1+y2)=2(1+x2)
dy g ( x )h( y ) 对于初值问题 dx y ( x 0 ) y0
法1: 法2: 可以先求通解再求特解 可以直接对方程两边同时求变上限定积分
dy g( x )dx h( y ) x0 y0
解:设曲线方程为y y( x ), 则 dy 2 2 y' 3 x , 或 3 x 且y(1) 2 dx 2 3 y 3 x dx x C ,