常微分方程讲义 (2)
高等数学常微分方程讲义,试题,答案
高等数学常微分方程讲义,试题,答案常微分方程§4.1 基本概念和一阶微分方程(甲)内容要点一、基本概念1、常微分方程和阶2、解、通解和特解3、初始条件4、齐次线性方程和非齐次线性方程二、变量可分离方程及其推广1、dyp(x)Q(y)dx(Q(y) 0) 2、齐次方程:dy dxy f x三、一阶线性方程及其推广1、dydyP(x)y Q(x) 2、P(x)y Q(x)y dxdx( 0,1)四、全微分方程及其推广(数学一)1、P(x,y)dx Q(x,y)dy 0,满足Q P2、P(x,y)dx Q(x,y)dy 0,五、差分方程(数学三)(乙)典型例题例1、求y x22Q p (RQ) (RP)但存在R(x,y),使x y x ydydyxy的通解。
dxdx解:y (x xy)22dy0dxydyy2 x d__y x2 y1 x2yduu2令u,则u x udx x(1 u)du 0xdxu 11 udxdu u x C1 ln|xu| u C1例2C1 uce, y cedyy的通解d__ y4uyx求微分方程d__ y4dx1解:此题不是一阶线性方程,但把x看作未知函数,y看作自变量,所得微分方程即x y3是一阶dyydyy11dy 14 dy 133yydy C y Cy 线性方程P(y) ,Q(y) y x e yey 3例3设y e是xy p(x)y x的一个解,求此微分方程满足yx ln2 0的特解xx解:将y e代入微分方程求出P(x) xe先求出对应齐次方程x,方程化为dy(e x 1)y 1 dxx xdy(e x 1)y 0的通解y cex e根据解的结构立刻可得非齐次方程通解y ex cex e dx再由yx ln2 0得2 2ec 0,c e例4设1212故所求解y e exx e x12满足以下件F(x) f(x)g(x),其中f(x),g(x)在( , )内f (x) g(x),g (x) f(x),且f(0) 0,f(x) g(x) 2ex(1)求F(x)所满足的一阶微分方程(2)求出F(x)的表达式解:(1)由F (x) f (x)g(x) f(x)g (x) g2(x) f2(x) [f(x) g(x)]2 2f(x)g(x) (2ex)2 2F(x) 可知F(x)所满足的一阶微分方程为F (x) 2F(x) 4e2x (2)F(x) e2dx4e2xe 2dxdx c e 2x 4e4xdx c e2x ce 2x将F(0) f(0)g(0) 0代入,可知c 1 于是例52F(x) e2x e 2xdy2(1 y)的通解求微分方程(y x) xdxsec2udusec3u 解:令y tanu,x tanv, 原方程化为(tanu tanv)secv2secvdv化简为sin(u v)dudzdudz 1 再令z u v,则1,方程化为sinz 1 sinz dvdvdvdv sinz(sinz 1) 1dz dv c, 1 sinz 1 sinzdz v c,1 sinzv c21 sinz1 sinz z v c 2coszz tanz secz v c z最后Z再返回x,y,v也返回x,即可。
常微分方程讲义精简
例2 求解方程 .解令,有原方程的参数形式为由基本关系式有积分得到从而原方程的参数形式通解为也可以消去参数t ,得到原方程的通积分为通解为例4 求解方程解令原方程的参数形式为(1.72)由基本关系式有或上式又可化为由,代入(1.72)的第三式,得原程的一个特解 .再由,解得,代入(1.72)的第三式,得原方程的通解例5求解方程(1.73)这里,假定是二次可微函数.解 (1.73)的参数形式为(1.74)由基本关系式有整理得由,得,代入(1.74)的第三式,得原方程通解(1.75)由于,由解得隐函数 ,代入(1.74)第三式,得到原方程的一个特解(1.76)(第7讲几种可降阶的高阶方程例1求解方程解令则有通解为从而积分四次,得到原方程的通解第二种可降阶的高阶方程例2求解方程.解令,则代入原方程得或积分后得"其中a"为任意常数. 解出p"得或积分后得其中 b为任意常数. 于是有或其中为任意常数.1.7.3恰当导数方程假如方程( 1.80)的左端恰为某一函数对 x的导数,即(1.80)可化为则(1.80)称为恰当导数方程.这类方程的解法与全微分方程的解法相类似,显然可降低一阶,成为之后再设法求解这个方程.例3求解方程.解易知可将方程写成故有即.积分后即得通解例4 求解方程.解先将两端同乘不为0的因子,则有故,从而通解为参数法第10讲解的延展2.3.1 延展解、不可延展解的定义定义2.1 设是初值问题(2,2)在区间上的一个解,如果(2,2)还有一个在区间上的解,且满足(1)(2)当时,则称解是可延展的,并称是在I2上的一个延展解.否则,如果不存在满足上述条件的解,则称是初值问题(2.2)的一个不可延展解,(亦称饱和解).这里区间I1和I2可以是开的也可以是闭的..3.2 不可延展解的存在性定义2.2设定义在开区域上,如果对于D上任一点,都存在以为中心的,完全属于D的闭矩形域R,使得在R上的关于y满足李普希兹条件,对于不同的点,闭矩形域R的大小以及常数N可以不同,则称在D上关于y满足局部李普希兹条件“柯西收敛准则收敛对,N,使当1.数列,就有,存在对,N,使当2.,时,总有.存在对,A> 0,使当3.,总有.”例1试讨论方程通过点(1,1)的解和通过点(3,-1)的解的存在区间.解此时区域D是整个平面.方程右端函数满足延展定理的条件.容易算出,方程的通解是故通过(1,1)的积分曲线为它向左可无限延展,而当x →2-0时,y →+∞, 所以,其存在区间为(-∞,2),参看图2-10.图 2-10通过(3,-1)的积分曲线为它向左不能无限延展,因为当x →2+0时,y →-∞,所以其存在区间为(2,+∞).顺便指出:这个方程只有解y = 0可以向左右两上方向无限延展.这个例子说明,尽管在整个平面满足延展定理条件,解上的点能任意接近区域D的边界,但方程的解的定义区间却不能延展到整个数轴上去.例2讨论方程解的存在区间.解方程右端函数在无界区域内连续,且对y满足李普希兹条件,其通解为过D1内任一点的初值解.图 2-11在(0,+∞)上有定义,且当x →+0时,该积分曲线上的点无限接近D1的边界线x = 0,但不趋向其上任一点(图2-11).在区域内的讨论是类似的.延展定理是常微分方程中一个重要定理.它能帮助我们确定解的最大存在区间.从推论和上面的例子可以看出,方程的解的最大存在区间是因解而异的.例3考虑方程及在平面上连续,试证明:对于任意及假设,方程满足的解都在(-∞,+∞)上存在.图 2-12证明根据题设,可以证明方程右端函数在整个平面上满足延展定理及存在与唯一性定理的条件.易于看到,为方程在(-∞,+∞)上的解.由延展定理可知,满足任意,的解上的点应当无限远离原点,但是,由解的唯一性,又不能穿过直线,故只能向两侧延展,而无限远离原点,从而这解应在(-∞,+∞)上存在(图2-12).2.4.1 奇解在本章 2.2节的例2中,我们已经看到方程的通解是,还有一解,除解外,其余解都满足唯一性,只有解所对应的积分曲线上每一点,唯一性都被破坏. 这样的解在许多方程中存在.例1求方程的所有解.解该方程的通解是此外还有两个特解和.由于该方程右端函数的根号前只取+号,故积分曲线如图2-13所示,图 2-13显然解和所对应的积分曲线上每一点,解的唯一性均被破坏。
常微分方程的几何解释
(2.2)
a x b, y ,
假设函数 f x, y在给定区域上连续且有界.于是
它在这个区域上确定了一个线素场.下面利用线素场
求出经过 x0, y0 的近似积分曲线.把
x0 ,b n 等分,其分点为:
xk x0 kh, k 0,1, , n
h b x0 , n
xn b
常微分方程
绵阳师范学院
先求出 f x0, y0
用经过 x0, y0 斜率为
y
x1
,
y1
x2
,
y2
f x0, y0 的直线段来近
y0
似积分曲线,其方程为
y y0 f x0, y0 x x0
x0 x1 x2
bx
求出直线上横坐标 x1 处的点的纵坐标
y1 y0 f x0, y0 x1 x0 y0 f x0, y0 h
如果 h 很小 x1, y1 就很接近积分曲线上的点 x1, y x1
因 f x, y 连续.于是由点 x1, y1 出发的斜率为
f x1, y1 的直线段又近似于原积分曲线.它的方程为
了线素场.
y k x
易见在点 x, y 的线素与
过原点与该点的射线重合.
常微分方程
绵阳师范学院
定理2.1 L为(2.1)的积分曲线的充要条件是: 在L 上任一点,L 的切线方向与(2.1)所确定的线 素场在该点的线素方向重合;即L在每间点均与 线素场的线素相切.
证明 必要性 设L为(2.1)的积分曲线,其方程为
20
若初值问题
dy dx
f ( x, y),的解是存在,是否唯一?
清华大学微积分高等数学课件第讲常微分方程二 共32页
29.07.2019
10
给 y(x0)y0 得 C y0
特解
x
x
yex0p(x)d(xy0xx0q(非x)齐e次x0特p(x解)dd x )x
非齐次通解的结构:
设y是y'p(x)y0 (2)的通 解 ,
y(x)是y'p(x)yq(x) (1)的 一 个 ,
则(1)的 通 解 为 y(x)yy(x)
代入方程并计算化简
yC (y) C (y) C (y) yye
C(y)ey
积分得 C(y)eyd yeyC
通解 xCyyey
29.07.2019
14
[例 3]设 a0,f(x)在 [0, )连 续,证 有明 界 方程
dxaxf(t) (t0) dt
每个[0 解 , 在 )有.界
x2 ydx x2ydyd(x2y2) 2
29.07.2019
23
[ 例 1 ]解(x 方 2 y ) d 程 ( x x y ) d 0 y
[解] 凑微分
x 2 d x (xd yyd )x yd 0 y
d(x3)d(x)yd(y2)0
3
2
d(x3 xyy2)0
3
2
通解
x3 xy y2 C
3
2
29.07.2019
24
[例 2] 解方 yd 程 x(y3ln x)d y0
x
[解] 改写为
(ydx lnxd) yy3dy 0 x
(yld n x ln x) d y y 3 d y 0
d(ylnx)d(y4)0 4
通解为
yl nx1 y4 C 4
例如 xd ydx d(x)y
常微分方程 ppt课件
量,x是未知函数,是未知函数对t导数. 现
在,我们还不会求解方程(1.1),但是,如果
考虑k=0的情形,即自由落体运动,此时方程
(1.1)可化为
d2x dt 2
g
(1.2)
将上式对t积分两次得
x(mt)xk12xgt2mgc1t c2
(1.3) (1.1)
ppt课件
11
一般说来,微分方程就是联系自变量、 未知函数以及未知函数的某些导数之间的关 系式. 如果其中的未知函数只是一个自变量 的函数,则称为常微分方程;如果未知函数 是两个或两个以上自变量的函数,并且在方 程中出现偏导数,则称为偏微分方程. 本书 所介绍的都是常微分方程,有时就简称微分 方程或方程.
这样,从定义1.1可以直接验证:
F(x, y, y) 0
(1.8)
如果在(1.8)中能将 y 解出,则得到方程
y f (x, y)
(1.9)
或
M (x, y)dx N(x, y)dy 0
(1.10)
(1.8)称为一阶隐式方程,(1.9)称为一阶显式方程,(1.10)称为微 分形式的一阶方程.
ppt课件
14
n 阶隐式方程的一般形式为
常微分方程
ppt课件
1
常微分方程课程简介
常微分方程是研究自然科学和社会科学中的事物、 物体和现象运动、演化和变化规律的最为基本的数 学理论和方法。物理、化学、生物、工程、航空航 天、医学、经济和金融领域中的许多原理和规律都 可以描述成适当的常微分方程,如牛顿运动定律、 万有引力定律、机械能守恒定律,能量守恒定律、 人口发展规律、生态种群竞争、疾病传染、遗
ppt课件
2
传基因变异、股票的涨伏趋势、利率的 浮动、市场均衡价格的变化等,对这些 规律的描述、认识和分析就归结为对相 应的常微分方程描述的数学模型的研究.
经济数学基础微积分课件 常微分方程
例2 验证函数 y e x e x 是不是方程
y 2 y y 0的解.
解 求 y e x e x 的导数,得 y e x e x , y e x e x
将y、y及y 代入原方程的左边,有
e x e x 2e x 2e x e x e x 0 即函数 y e x e x 不满足原方程,
前页 后页 结束
M1(x) N1(x)
d
x
N2(y) M 2( y)
d
y
0
将(9.2.3)式两边积分后,
(9.2.3)
M1(x) N1(x)
d
x
N2(y) M 2( y)
d
y
C
(C为任意常数)
可验证,此结果即用隐式给出的方程(9.2.3)的通解.
约定:
在微分方程这一章中不定积分式表示被积函数的一
y e p(x)d x q(x)e p(x)d x d x C
即为所求(9.3.1)的通解.
前页 后页 结束
例1 求微分方程 dy 2xy 2xe x2 的通解. dx
解 p(x) 2x, q(x) 2xex2
代入公式
y e2xd x 2xex2 e2xd x d x C
常微分方程
9.1 常微分方程的基本概念 9.2 可分离变量的微分方程 9.3 一阶微分方程与可降阶
的高阶微分方程 9.4 二阶常系数微分方程 9.5 常微分方程的应用举例
结束
9.1 常微分方程的基本概念
定义一 含有未知函数的导数(或微分)的方程称为 微分方程。
常微分方程:未知函数是一元函数的微分方程 偏微分方程:未知函数是多元函数的微分方程 定义二 在微分方程中,所出现的未知函数的最高阶
清华大学微积分高等数学课件第讲常微分方程二教案资料
2020/6/18
15
设 f(x) M
x(t)x0
tea(ts)f(s)ds
0
x0 M
tea(ts)ds
0
x0
M a
2020/6/18
16
二、伯努利(Bernoulli)方程
Bernoulli 方程
dyp(x)yq(x)yn dx
方程两端同除yn
yndyp(x)y1nq(x) dx
令z y1n
例如 xd ydx d(x)y
xdxyd yd(x2y2) 2
xdx y2ydxd(x y)
ydxy2xdyd(xy)
2020/6/18
22
xxd2 yyy2dxd(arctxya)n
yxd2x x y2dyd(arctxya)n
2020/6/18
8
C (x )y 1 (x ) C (x )y 1 '(x ) p (x )C (x )y 1 (x ) q (x )
y1(x)是 ( 2) 的 解 ,
C (x )y 1 '(x ) p (x ) C (x )y 1 (x ) 0 化简得到 C (x)y1(x)q(x)
即 C(x)q(x)ep(x)dx
通解 xCyyey
2020/6/18
14
[例 3] 设 a0,f(x)在 [0, )连 续,证 有明 界 方程
dxaxf(t) (t0) dt
每个[0 解 , 在 )有.界
[证] 设xx(t)是 满 足 初x始 (0)条 x0 件 的 解 .
则 x (t) e a(tx 0 0 teafs (s)d)s(t0 )
2020/6/18
9
积分 C(x) q(x)ep(x)dxC
常微分方程讲义和作业
第四章 常微分方程与数学模型微积分最主要的应用可能就是微分方程了,在物理学、力学、工程技术、经济学和管理科学等实际问题中具有广泛的应用。
一、什么是微分方程例1:含有未知函数的导数或微分的方程称为微分方程,例如()dyu x dx=,其中()y f x =为未知函数,()u x 为已知函数。
满足上述方程的函数()y f x =称为微分方程的解。
求下列微分方程满足所给条件的解: (1)2(2)dyx dx=-,20x y ==; (2)2232d x dt t =,11t dx dt ==,11t x ==。
二、分离变量法※例2:求微分方程y xy '=的通解。
解: 变形为:dy xy dx =, 分离变量:1dy xdx y=(此时漏掉解0y =), 两边同时积分:1dy xdx y =⎰⎰, 得:211ln 2y x C =+, 22111122x C x C y ee e+==,从而22111222x x C y e eC e =±=,其中12CC e =±,为任意非零常数,但0y =亦是方程的解,统一起来,方程的通解为:212x y Ce=,C 为任意常数。
上述求解过程比较繁琐,由于经常出现,为方便计,从分离变量后开始将求解过程简写为:两边同时积分:1dy xdx y =⎰⎰, 得:21ln ln 2y x C =+, 从而 2211ln 22xx C y e e Ce==这个过程严格说是有问题的,但比较简洁,又能得到正确的结果,所以常被采用。
例3:(1)牛顿冷却定律指出:如果物体和周围环境之间的温度相差不是很大的话,物体冷却速度与温差成正比(同样可用于加热的情况)。
命()T t 表示在时刻t 物体的温度,c T 表示周围环境的温度(假定是常数),建立微分方程并求解,得出()T t 的变化规律。
(2)清晨,警察局接到报案,街头发现一具死尸,6:30时测量体温为18℃,7:30时再测一次为16℃,室外温度为10℃(假定不变),人正常体温为37℃,请估计被害人何时死亡?(死亡时刻记为0t ,则0()37T t =,时刻6:30计算时看成6.5)例4:人口预测记时刻t 的人口为()P t ,当考察一个国家或一个较大地区的人口时,()P t 是一个很大的整数,为了利用微积分这一数学工具,将()P t 视为连续、可微函数.记初始时刻(0)t =的人口为0P ,假设人口增长的速度(即增长率)与t 时刻的人口数量()P t 成正比,利用下表中数据为20世纪世界人口建模,增长率是多少,建立的模型与数据相符合吗?解:设比例系数为μ(即增长率),则()P t 满足的微分方程为:0,(0)dPP P P dtμ==. 解出 0()tP t Pe μ= , 表明人口将按指数规律随时间无限增长(0μ>).上式称为人口指数增长模型,也称为马尔萨斯人口模型.以1900年为初始时刻,0(0)=1650P P =,得()1650tP t e μ=, 以1910年数据估计μ,即10(10)16501750P e μ==,解11750l n .0584101650μ=≈,即增长率约为0.6%,增长模型为0.005884()1650t P t e =若以1950年为初始时刻,为20世纪后50年建模,则0=2560P ,得()2560tP t e μ=,以1960年数据估计μ,即10(10)25603040P e μ==,解13040l n 0.017185102560μ=≈,即增长率约为1.7%,增长模型为0.017185()2560t P t e =但是长期来看,任何地区的人口都不可能无限增长,即指数模型不能描述、也不能预测较长时期的人口演变过程,这是因为人口增长率事实上是不断地变化着.排除灾难、战争等特殊时期,一般来说,当人口较少时,其增长较快,即增长率较大;人口增加到一定数量后,增长就会慢下来,即增长率变小.看来,为了使人口预测特别是长期预测能更好地符合实际情况,必须修改人口指数增长模型中关于人口增长率是常数这个基本假设.2.人口阻滞增长模型(Logistic 模型)分析人口增长到一定数量后增长率下降的主要原因,人们注意到,自然资源、环境条件等因素对人口的增长起着阻滞作用,并且随着人口的增加,阻滞作用越来越大.所谓人口阻滞增长模型就是考虑到这个因素,对人口指数增长模型的基本假设进行修改后得到的.阻滞作用体现在对人口增长率μ的影响上,使得μ随着人口数量P 的增加而下降。
常微分方程讲义全文
6、恰当方程
M (x, y)dx + N (x, y)dy = 0
判定:全微分 ⇔ ∂M ≡ ∂N ∂y ∂x
x
y
∫ ∫ 解法: M (x, y)dx + x0
y0 N (x0 , y)dy = C
初值问题: C = 0
例 2xydx + (x2 − y2 )dy = 0
解: ∂M ∂y
≡ ∂N ∂x
uz′ = −(z −1)(z − 2) /(z + 1)
z = 1, z = 2 ⇔ v = u, v = 2u ⇔ y = x + 1, y = 2x
⎛ ⎝⎜
z
3 −
2
−
z
2 −
1
⎞ ⎠⎟
dz
= − du u
⇒
(z − 2)3 (z −1)2
= C /u
( y − 2x)3 = C( y − x −1)2
一阶 线性 二阶 线性 一阶 非线性
齐方程、非齐次方程
在方程中,不含未知函数及其导数的项,称为自由项。 自由项为零的方程,称为齐方程。 自由项不为零的方程,称为非齐方程。
d x = x2 dt
一阶齐线性方程
d2 y d x2
+
b
d d
y x
+
cy
=
sin
x
二阶非齐线性方程
⎜⎛ d x ⎞⎟2 − x2 = t3 ⎝ dt ⎠
一阶非齐非线性方程
微分方程的一般表示形式
n 阶微分方程的一般形式 为 F (x, y′, y′′,L, y(n) ) = 0 。
F
(x,
y′,
y′′)
常微分方程讲义
常微分方程讲义微分方程是数学的一个重要分支,它的本质是求解某个函数的微分(偏导数)方程等式,并得出相应的函数解。
因此,它也被称为“函数微分学”。
微分方程常常用于研究物理和其他科学的解析理论上的问题,比如力学、流体力学、声学、电磁学等方面的研究。
一般来说,微分方程包括微分解析方程、积分方程和偏微分方程,其中,最常用的是偏微分方程。
它是由一个或多个复变量函数的某个变量(或多个变量)的偏导数组成的方程,而它的解就是被偏微分方程包含的函数。
偏微分方程可以分为常微分方程和时滞微分方程,本文讲义主要介绍的是前者,即常微分方程。
常微分方程是由一个复变量函数的某个变量的导数组成的方程式,它的解是一个关于变量的函数。
它的基本思想是:将某些可变量的函数表示为可以用一个或几个未知函数的函数,求解该未知函数,从而求解特定问题所对应的函数,用以描述和分析物理系统的特性。
常微分方程可以通过三种方式求解:第一种是数值方法,即将微分方程的求解转换成一系列的算数计算,它是最常用的解法;第二种是图像方法,它是通过拟合图形来确定方程的解的;第三种是函数解法,即求解方程的解析表达式,它也是研究微分方程的重要方法。
如何求解一般常微分方程?一般来说,要先将原始方程化为不带高次导数(称为常数阶微分方程)或不带高次导数和常数(称为普通微分方程)的形式,然后再进行解算。
这些方程又可以分为线性微分方程、二阶微分方程、高阶微分方程及一类特殊微分方程。
线性微分方程是指形如y′+ay=f(x)的微分方程,它的解可以通过谱解的方法求出,就是将此方程转换为一个定义域上的线性算子的本征方程,再根据本征方程的本质解其求解。
二阶微分方程是指形如y′′+ay′+by=f(x)的微分方程,它的解可以利用解析方法或特殊求解的方法求得,常见的有求解公式或积分方法。
高阶微分方程是指形如y′(n)+ay′(n-1)+…+by=f(x)的微分方程,它的解是求解公式,这种公式只有当所求解的方程满足某些条件时才可以得出,如果不满足,就只能利用特殊的解法来解高阶方程。
考研数学(二)考试大纲解析(常微分方程)【圣才出品】
第 5 章 常微分方程
一、微分方程
表 5-1 微分方程
微分方程
概念 凡表示未知函数、未知函数的导数与自变量之问的关系的方程,称为微分方程.
阶
微分方程中所出现的未知函数的最高阶导数的阶数,称为微分方程的阶.
设函数 y (x) 在区间 I 上有 n 阶连续导数,如果在区间 I 上,有
1 / 13
圣才电子书 十万种考研考证电子书、题库视频学习平台
如果一个一阶微分方程能写成
g y dy f x dx
的形式,则原方程称为可分离变量的微分方程.
2.解法
设 y x 是方程 g y dy f x dx 的解,代入得
三、一阶线性微分方程 1.一阶线性微分方程
dy P(x) y Q(x) dx
2.齐次线性方程 (1)表达式
(2)通解
dy P(x) y 0 dx
2 / 13
圣才电子书 十万种考研考证电子书、题库视频学习平台
y CeP(x)dx (C eC1 ) (C1为常数)
F[x, (x), (x), , (n)(x)] 0
解
则函数 y (x) 就称为微分方程
F(x, y, y, , y(n)) 0
在区间 I 上的解.
通解
如果微分方程的解中含有任意常数,且任意常数的个数与微分方程的阶数相
同,则这样的解称为微分方程的通解.
设微分方程中的未知函数为 y (x) ,则:
方程
y f ( y, y)
令 y p
得
y dp dp dy p dp dx dy dx dy
常微分方程讲义(2)
常微分方程讲义(二)
例:一曲线,其任意一点的切线界于两座标轴之间的部分被切点分成相等的两段,求该曲线满足的微分方程的表达式并解之。
注重几何意义
在什么时候添加积分常量“C ”?
线素场
线素决定曲线的走势;
利用线素场描画曲线走势的好处在于不解出曲线方程也能大致画出曲线的性状(根据Liouville ,1841:不是每个积分都能用初等积分法求出最后的解的)
但线素场的方法不是万能的。
例:用线素场法画x dx
dy
2=的解的图像 用线素场法画y x dx dy +=的解的图像
可分离变量的微分方程的解法
)()(y h x g dx
dy = dx x g dy y h )()
(1=⇒
⎰⎰=⇒dx x g dy y h )()
(1
求微分方程的解,关键在积分——所以娴熟的积分技巧非常必须。
“拆”的思路: 例:)1(2
12y dx
dy
-= y y dx
dy -=2 2by ay dx
dy -= 3242y y dx dy x =+ y
x xy y dx dy 321++=
例:利息、复利与利率期限结构
复计的实际利率≥名义利率
例:关于人口普查
高昂的成本
用于预测的Malthus 方程
我们为什么还需要不定期的人口普查?
例:禽流感与银行倒闭:政府能否不进行干预?。
常微分方程考研讲义第二章 一阶微分方程的初等解法
第二章、一阶微分方程的初等解法[教学目标]1. 理解变量分离方程以及可化为变量分离方程的类型(齐次方程),熟练掌握变量分离方程的解法。
2. 理解一阶线性微分方程的类型,熟练掌握常数变易法及伯努力方程的求解。
3. 理解恰当方程的类型,掌握恰当方程的解法及简单积分因子的求法。
4. 理解一阶隐式方程的可积类型,掌握隐式方程的参数解法。
[教学重难点] 重点是一阶微分方程的各类初等解法 ,难点是积分因子的求法以及隐式方程的解法。
[教学方法] 讲授,实践。
[教学时间] 14学时[教学内容] 变量分离方程,齐次方程以及可化为变量分离方程类型,一阶线性微分方程及其常数变易法,伯努利方程,恰当方程及其积分因子法,隐式方程。
[考核目标]1.一阶微分方程的初等解法:变量分离法、一阶线性微分方程的常数变易法、恰当方程与积分因子法、一阶隐方程的参数解法。
2.会建立一阶微分方程并能求解。
§1 变量分离方程与变量变换1、 变量分离方程1) 变量分离方程形如 ()()dy f x g y dx= (或1122()()()()0M x N y dx M x N y dy +=) (2.1) 的方程,称为变量分离方程,其中函数()f x 和()g y 分别是,x y 的连续函数.2) 求解方法如果()0g y ≠,方程(2.1)可化为,()()dy f x dx g y = 这样变量就分离开了,两边积分,得到 ()()dy f x dx c g y =+⎰⎰ (2.2) 把,()()dy f x dx g y ⎰⎰分别理解为1,()()f x y ϕ的某一个原函数. 容易验证由(2.2)所确定的隐函数(,)y x c ϕ=满足方程(2.1).因而(2.2)是(2.1)的通解.如果存在0y 使0()0g y =,可知0y y =也是(2.1)的解.可能它不包含在方程的通解(2.2)中,必须予以补上.3) 例题例1 求解方程dy x dx y=- 解 将变量分离,得到ydy xdx =-两边积分,即得 22222y x c =-+ 因而,通解为22x y c += 这里的c 是任意的正常数.或解出显式形式y =例2 解方程 2cos dy y x dx= 并求满足初始条件:当0x =时.1y =的特解.解 将变量分离,得到2cos dy xdx y = 两边积分,即得 1sin x c y -=+ 因而,通解为 1sin y x c=-+ 这里的c 是任意的常数.此外,方程还有解0y =.为确定所求的特解,以0x =.1y =代入通解中确定常数c ,得到 1c =-因而,所求的特解为 11sin y x =- 例3 求方程 ()dy P x y dx= (2.3) 的通解,其中()P x 是x 的连续函数.解 将变量分离,得到()dy P x dx y= 两边积分,即得 ln ()y P x dx c=+⎰ 这里的c是任意常数.由对数的定义,即有 ()P x dx c y e +⎰=即()P x dx c y e e ⎰=±令c e c ±= ,得到()P x dx y ce ⎰= (2.4)此外,0y =也是(2.3)的解.如果在(2.4)中允许0c =,则0y =也就包括在(2.4)中,因而,(2.3)的通解为(2.4),其中c 是任意常数.注: 1.常数c 的选取保证(2.2)式有意义.2.方程的通解不一定是方程的全部解,有些通解包含了方程的所有解,有些通解不能包含方程的所有解.此时,还应求出不含在通解中的其它解, 即将遗漏的解要弥补上.3.微分方程的通解表示的是一族曲线,而特解表示的是满足特定条件00()y x y =的一个解,表示的是一条过点00(,)x y 的曲线.2、可化为变量分离方程的类型1).形如 dy y g dx x ⎛⎫= ⎪⎝⎭(2.5) 的方程,称为齐次方程,这里的()g u 是u 的连续函数.另外,ⅰ)对于方程 (,)(,)dy M x y dx N x y = 其中函数(,)M x y 和(,)N x y 都是x 和y 的m 次齐次函数,即对0t >有(,)(,)m M tx ty t M x y ≡ (,)(,)m N tx ty t N x y ≡ 事实上,取1t x=,则方程可改写成形如(2.5)的方程. (1,)(1,)(1,)(1,)m m y y x M M dy x x y y dx x N N x x== ⅱ)对方程 (,)dy f x y dx= 其中右端函数(,)f x y 是x 和y 的零次齐次函数,即对0t >有(,)(,)f tx ty f x y =则方程也可改写成形如(2.5)的方程(1,)dy y f dx x= 对齐次方程(2.5)利用变量替换可化为变量分离方程再求解.令y u x =(2.6) 即y ux =,于是 dy du x u dx dx=+ (2.7) 将(2.6)、(2.7)代入(2.5),则原方程变为 ()du xu g u dx += 整理后,得到 ()du g u u dx x-= (2.8) 方程(2.8)是一个可分离变量方程,按照变量分离法求解,然后将所求的解代回原变量,所得的解便是原方程(2.5)的解.例4 求解方程dy y y tg dx x x=+ 解 这是齐次方程,以,y dy du u x u x dx dx ==+代入,则原方程变为 du xu u tgu dx +=+ 即du tgu dx x= (2.9) 分离变量,即有 dx ctgudu x=两边积分,得到 ln sin ln u x c=+这里的c是任意的常数,整理后,得到 sin u cx = (2.10)此外,方程(2.9)还有解0tgu =,即sin 0u =. 如果(2.10)中允许0c =,则sin 0u =就包含在(2.10)中,这就是说,方程(2.9)的通解为(2.10).代回原来的变量,得到原方程的通解为 sin ycx x =例5 求解方程(0).dyx y x dx +=<解 将方程改写为(0)dyyx dx x =< 这是齐次方程,以,ydyduu x u x dx dx ==+代入,则原方程变为dux dx =(2.11) 分离变量,得到dxx =两边积分,得到(2.11)的通解ln()x c =-+即2[ln()](ln()0)u x c x c =-+-+> (2.12) 这里的c 是任意常数.此外,(2.11)还有解0u =注意,此解不包括在通解(2.12)中.代回原来的变量,即得原方程的通解2[ln()](ln()0)y x x c x c =-+-+>及解0y =.原方程的通解还可表为2[ln()],ln()0,0,x x c x c y ⎧-+-+>=⎨⎩它定义于整个负半轴上. 注:1.对于齐次方程dy y g dx x ⎛⎫= ⎪⎝⎭的求解方法关键的一步是令y u x =后,解出y ux =,再对两边求关于x 的导数得dy du u x dx dx =+,再将其代入齐次方程使方程变为关于,u x 的可分离方程.2.齐次方程也可以通过变换x v y=而化为变量分离方程.这时x vy =,再对两边求关于y 的导数得dx dv v y dy dy =+,将其代入齐次方程dx x f dy y ⎛⎫= ⎪⎝⎭使方程变为,v y 的可分离方程 小结:这一讲我们主要讲解了一阶微分方程的可分离变量法和齐次方程的dy y g dx x ⎛⎫= ⎪⎝⎭形状的解法.而这一齐次方程通过变量替换任然可化为可分离方程,因而,一定要熟练掌握可分离方程的解法.2)形如 111222a xb yc dy dx a x b y c ++=++ (2.13) 的方程经变量变换化为变量分离方程,这里的121212,,,,,a a b b c c 均为常数.分三种情况来讨论(1)120c c ==情形.这时方程(2.13)属齐次方程,有 1122a x b y dy y g dx a x b y x +⎛⎫== ⎪+⎝⎭此时,令y u x=,即可化为变量可分离方程. (2)11220a b a b =,即1122a b a b =的情形.设1122a b k a b ==,则方程可写成22122222()()()k a x b y c dy f a x b y dx a x b y c ++==+++ 令22a x b y u +=,则方程化为 22()du a b f u dx=+ 这是一变量分离方程.(3)1112220,a b c c a b ≠及不全为零的情形.这时方程(2.13)右端的分子、分母都是,x y 的一次式,因此11122200a xb yc a x b y c ++=⎧⎨++=⎩ (2.14) 代表xy 平面上两条相交的直线,设交点为(,)αβ.显然,0α≠或0β≠,否则必有120c c ==,这正是情形(1)(只需进行坐标平移,将坐标原点(0,0)移至(,)αβ就行了,若令X x Y y αβ=-⎧⎨=-⎩(2.15) 则(2.14)化为112200a X b Y a X b y +=⎧⎨+=⎩ 从而(2.13)变为 1122a X bY dY Y g dX a X b Y X +⎛⎫== ⎪+⎝⎭(2.16) 因此,得到这种情形求解的一般步骤如下:(1)解联立代数方程(2.14),设其解为,x y αβ==;(2)作变换(2.15)将方程化为齐次方程(2.16);(3)再经变换Y u X=将(2.16)化为变量分离方程; (4)求解上述变量分离方程,最后代回原变量可得原方程(2.13)的解.上述解题的方法和步骤也适用于比方程(2.13)更一般的方程类型 111222a x b y c dy f dx a x b y c ⎛⎫+== ⎪++⎝⎭此外,诸如 ()dy f ax by c dx++ ()()0y xy dx xg xy dy += 2()dy x f xy dx =2dy y xf dx x ⎛⎫= ⎪⎝⎭ 以及(,)()(,)()0M x y xdx ydy N x y xdy ydx ++-=(其中,M N 为,x y 的齐次函数,次数可以不相同)等一些方程类型,均可通过适当的变量变换化为变量分离方程.例6 求解方程 13dy x y dx x y -+=+- (2.17) 解 解方程组 1030x y x y -+=⎧⎨+-=⎩得1, 2.x y == 令12x X y Y =+⎧⎨=+⎩ 代入方程(2.17),则有dY X Y dX X Y-=+ (2.18) 再令Y u X =即 Y uX = 则(2.18)化为2112dX u du X u u +=-- 两边积分,得 22ln ln 21X u u c=-+-+ 因此22(21)c X u u e +-=±记1,c e c ±= 并代回原变量,就得2212Y XY X c +-=221(2)2(1)(2)(1)y x y x c -+----=此外,易验证2210u u +-=即2220Y XY X +-=也就是(2.18)的解.因此方程(2.17)的通解为22262y xy x y x c +---=其中c 为任意的常数. 3、 应用举例例7 电容器的充电和放电如图(2.1)所示的R C -电路,开始时电容C 上没有电荷,电容两端的电压为零.把开关K 合上“1”后,电池E 就对电容C 充电,电容C 两端的电压C u 逐渐升高,经过相当时间后,电容充电完毕,再把开关K 合上“2”,这时电容就开始放电过程,现在要求找出充、放电过程中,电容C 两端的电压C u 随时间t 的变化规律.解 对于充电过程,由闭合回路的基尔霍夫第二定理,c u RI E += (2.19)对于电容C 充电时,电容上的电量Q 逐渐增多,根据C Q Cu =,得到 ()C C du dQ d I Cu C dt dt dt=== (2.20) 将(2.20)代入(2.19),得到c u 满足的微分方程 c c du RC u E dt+= (2.21) 这里R 、C 、E 都是常数.方程(2.21)属于变量分离方程.将(2.21)分离变量,得到C C du dt u E RC =-- 两边积分,得到 11ln C u E t c RC -=-+ 即 1112t t c RC RC C u E e ec e ---=±= 这里12c c e =±为任意常数.将初始条件:0t =时,0C u =代入,得到2c E =-.所以 1(1)t RC C u E e -=- (2.22)这就是R C -电路充电过程中电容C 两端的电压的变化规律.由(2.22)知道,电压C u 从零开始逐渐增大,且当t →+∞时,C u E →,在电工学中,通常称RC τ=为时间常数,当3t τ=时,0.95C u E =,就是说,经过3τ的时间后,电容C 上的电压已达到外加电压的95%.实用上,通常认为这时电容C 的充电过程已基本结束.易见充电结果C u E =.对于放电过程的讨论,可以类似地进行.例8 探照灯反射镜面的形状在制造探照灯的反射镜面时,总是要求将点光源射出的光线平行地射出去,以保证照灯有良好的方向性,试求反射镜面的几何形状.解 取光源所在处为坐标原点,而x 轴平行于光的反射方向,设所求曲面由曲线()0y f x z =⎧⎨=⎩ (2.23) 绕x 轴旋转而成,则求反射镜面的问题归结为求xy 平面上的曲线()y f x =的问题,仅考虑0y >的部分,过曲线()y f x =上任一点(,)M x y 作切线NT ,则由光的反射定律:入射角等于反射角,容易推知12αα=从而 OM ON =注意到 2dy MP tg dx NPα==及,,OP x MP y OM ===就得到函数()y f x =所应满足的微分方程式dy dx = (2.24) 这是齐次方程.由2.12知引入新变量x u y =可将它化为变量分离方程.再经直接积分即可求得方程的解.对于方齐次方程(2.24)也可以通过变换x v y =而化为变量分离方程也可由x yv =得dx dv v y dy dy=+代入(2.24)得到sgn dv v y v y dy+=+于是sgn dy y y =(2.25)积分(2.25)并代回原来变量,经化简整理,最后得2(2)y c c x =+ (2.26)其中c 为任意常数.(2.26)就是所求的平面曲线,它是抛物线,因此,反射镜面的形状为旋转抛物面22(2)y z c c x +=+ (2.27)小结: 本节我们主要讨论了一阶可分离微分方程和齐次微分方程的求解问题.将各种类型的求解步骤记清楚的同时要注意对解的讨论.§2 线性方程与常数变易法1、一阶线性微分方程 ()()()0dy a x b x y c x dx++= 在()0a x ≠的区间上可以写成 ()()dy P x y Q x dx=+ (2.28) 对于()a x 有零点的情形分别在()0a x ≠的相应区间上讨论.这里假设(),()P x Q x 在考虑的区间上是x 的连续函数.若()0Q x ≡,(2.28)变为 ()dy P x y dx= (2.3) 称为一阶齐线性方程.若()0Q x ≠,(2.28)称为一阶非齐线性方程.2、常数变易法(2.3)是变量分离方程,已在例3中求得它的通解为()P x dx y ce ⎰= (2.4)这里c 是任意的常数.下面讨论一阶非齐线性方程(2.28)的求解方法.方程(2.3)与方程(2.28)两者既有联系又有区别,设想它们的解也有一定的联系,在(2.4)中c 恒为常数时,它不可能是(2.28)的解,要使(2.28)具有形如(2.4)的解, c 不再是常数,将是x 的待定函数()c x ,为此令()()P x dx y c x e ⎰= (2.29)两边微分,得到 ()()()()()P x dx P x dx dy dc x e c x P x e dx dx ⎰⎰=+ (2.30) 将(2.29)、(2.30)代入(2.28),得到()()()()()()()()()P x dx P x dx P x dx dc x e c x P x e P x c x e Q x dx ⎰⎰⎰+=+ 即()()()P x dx dc x Q x e dx -⎰= 积分后得到()()()P x dx c x Q x e dx c -⎰=+⎰(2.31) 这里c是任意的常数..将(2.31)代入(2.29),得到()()()()()() =()P x dx P x dx P x dx P x dx P x dx y e Q x e dx c ce e Q x e dx --⎛⎫⎰⎰=+ ⎪⎝⎭⎰⎰⎰+⎰⎰(2.32) 这就是方程(2.28)的通解.这种将常数变易为待定函数的方法,通常称为常数变易法.实际上常数变易法也是一种变量变换的方法.通过变换(2.29)可将方程(2.28)化为变量分离方程.注: 非齐线性方程的通解是它对应的齐线性方程的通解与它的某个特解之和.例1 求方程1(1)(1)x n dy x ny e x dx++-=+的通解,这里的n 为常数. 解 将方程改写为(1)1x n dy n y e x dx x -=++ (2.33) 先求对应的齐次方程01dy n y dx x -=+ 的通解,得(1)n y c x =+令 ()(1)n y c x x =+ (2.34)微分之,得到 ()(1)(1)()n dy dc x x n x c x dx dx=+++ (2.35) 以(2.34)、(2.35)代入(2.33),再积分,得()x c x e c=+ 将其代入公式(2.34),即得原方程的通解(1)()n x y x e c=++ 这里c是任意的常数. 例2 求方程22dy y dx x y =-的通解. 解 原方程改写为 2dx x y dy y=- (2.36) 把x 看作未知函数,y 看作自变量,这样,对于x 及dx dy 来说,方程(2.36)就是一个线性方程了.先求齐线性方程2dx x dy y= 的通解为2x cy = (2.37)令2()x c y y =,于是2()2()dx dc y y c y y dy dy=+ 代入(2.36),得到 ()ln c y y c=-+ 从而,原方程的通解为 2(ln )x y cy =- 这里c是任意的常数,另外0y =也是方程的解. 特别的,初值问题00()()()dy P x y Q x dx y x y ⎧=+⎪⎨⎪=⎩ 的解为0000()()()=()x x s x x x P d P d P d x x y cee Q s e ds ττττττ-⎰⎰⎰+⎰例3 试证 (1)一阶非齐线性方程(2.28)的任两解之差必为相应的齐线性方程(2.3)之解;(2)若()y y x =是(2.3)的非零解,而 ()y y x =是(2.28)的解,则(2.28)的通解可表为 ()()y cy x y x =+,其中c 为任意常数.(3)方程(2.3)任一解的常数倍或两解之和(或差)仍是方程(2.3)的解.证 (1)设12,y y 是非齐线性方程的两个不同的解,则应满足方程使 1122()(1)()(2)dy py Q x dx dy py Q x dx=+=+ (1)—(2)有1212()()d y y p y y dx-=- 说明非齐线性方程任意两个解的差12y y -是对应的齐次线性方程的解.(2)因为 (()())()()(()()()()d cy x y x dy x d y x c p cy p y Q x p cy y Q x dx dx dx+=+=++=++ 故结论成立.(3)因为12121212()()()(),(),()d y y d y y d cy p cy p y y p y y dx dx dx+-==+=- 故结论成立. 3、Bernoulli 方程形如 ()()n dy P x y Q x y dx=+ ( 0,1n ≠) (2.38) 的方程,称为伯努利(Bernoulli )方程,这里(),()P x Q x 为x 连续函数.利用变量变换可将伯努利方程化为线性方程来求解.事实上,对于0y ≠,用n y -乘(2.38)两边,得到 1()()nn dy y y P x Q x dx --=+ (2.39) 引入变量变换1n z y -= (2.40)从而 (1)n dz dy n y dx dx-=- (2.41) 将(2.40)、2.41)代入(2.39),得到 (1)()(1)()dz n P x z n Q x dx=-+- (2.42) 这是线性方程,用上面介绍的方法求得它的通解,然后再代回原来的变量,便得到(2.38)的通解.此外,当0n >时,方程还有解0y =.例4 求方程26dy y xy dx x=-的通解 解 这是2n =时的伯努利方程,令 1z y -=,得2dz dy y dx dx -=- 代入原方程得到 6dz z x dx x=-+ 这是线性方程,求得它的通解为 268c x z x =+ 代回原来的变量y ,得到 2618c x y x =+ 或者 688x x c y -= 这是原方程的通解. 此外,方程还有解0y =.例5 求方程331dy dx xy x y =+的解 解 将方程改写为 33dx yx y x dy=+ 这是一个自变量为y ,因变量为x 的伯努利方程.解法同上.例6 求方程23y dy e x dx x+=的通解 这个方程只要做一个变换,令,y y du dy u e e dx dx==,原方程改写为 22231du x u u dx x x=+便是伯努利方程.小结;这次主要讨论了一阶线性微分方程的解法.其核心思想是常数变易法.即将非齐线性方程对应的齐线性方程解的常数变易为待定函数,使其变易后的解函数代入非齐次线性方程,求出待定函数()c x ,求出非齐次方程的解.我们还讨论了伯努利方程,求解过程为,先变换,将原方程化为非齐线性方程,再求解.§3 恰当方程与积分因子1、恰当方程的定义将一阶微分方程(,)dy f x y dx = 写成微分的形式(,)0f x y dx dy -=把,x y 平等看待,对称形式的一阶微分方程的一般式为(,)(,)0M x y dx N x y dy += (2.43)假设(,),(,)M x y N x y 在某区域G 内是,x y 的连续函数,而且具有连续的一阶偏导数. 如果存在可微函数(,)u x y ,使得(,)(,)du M x y dx N x y dy =+ (2.44)即 (,), (,)u u M x y N x y x y∂∂==∂∂ (2.45) 则称方程(2.43)为恰当方程,或称全微分方程.在上述情形,方程(2.43)可写成(,)0du x y ≡,于是(,)u x y C ≡就是方程(2.43)的隐式通解,这里C 是任意常数(应使函数有意义).2、 恰当方程的判定准则定理1设(,),(,)M x y N x y 在某区域G 内连续可微,则方程(2.43)是恰当方程的充要条件是, (,)M N x y G y x∂∂=∈∂∂ (2.46) 而且当(2.46)成立时,相应的原函数可取为000(,)(,)(,)x y x y u x y M s y ds N x t dt =+⎰⎰ (2.47) 或者也可取为000(,)(,)(,)yxy x u x y N x t dt M s y ds =+⎰⎰ (2.48) 其中00(,)x y G ∈是任意取定的一点.证明 先证必要性.因为(2.43)是恰当方程,则有可微函数(,)u x y 满足(2.45),又知(,),(,)M x y N x y 是连续可微的,从而有 22M u u N y y x x y x∂∂∂∂===∂∂∂∂∂∂ 下面证明定理的充分性,即由条件(2.46),寻找函数(,)u x y ,使其适合方程(2.45).从(2.47)可知 (,)u N x y y∂=∂ 0000(,)(,) =(,)(,) =(,)(,)(,)y y y x y y y y u M x y N x t dt x x M x y N x t dt M x y M x t dt M x y ∂∂=+∂∂++=⎰⎰⎰ 即(2.45)成立,同理也可从(2.48)推出(2.45).例1. 解方程 21()02x xydx dy y++= (2.49)解 这里21, =()2x M xy N y=+,则y x M x N ==,所以(2.49)是恰当方程.因为N 于0y =处无意义,所以应分别在0y >和0y <区域上应用定理2.3,可按任意一条途径去求相应的原函数(,)u x y .先选取00(,)(0,1)x y =,代入公式(2.47)有 22011()ln 22xyx x u xdx dy y y y =++=+⎰⎰再选取00(,)(0,1)x y =-,代入公式(2.47)有22011()()ln()22xyx x u x dx dy y y y -=-++=+-⎰⎰可见不论0y >和0y <,都有2ln ||2x u y y =+故方程的通解为2ln ||2x y y C +=.3、恰当方程的解法上述定理已给出恰当方程的解法,下面给出恰当方程的另两种常用解法. 解法1. 已经验证方程为恰当方程,从(,)x u M x y =出发,有2(,)(,)()()2x u x y M x y dx y y y φφ≡+=+⎰ (2.50)其中()y φ为待定函数,再利用(,)y u N x y =,有221()22x x y yφ'+=+ 从而1()y y φ'=于是有 ()ln ||y y φ=只需要求出一个(,)u x y ,因而省略了积分常数.把它代入(2.50)便得方程的通解为2ln ||2x u y y C =+= 解法2. 分项组合的方法 对(2.49)式重新组合变为21()02x xydx dy dy y++=于是 2()ln ||02x d y d y +=从而得到方程的通解为 2ln ||2x y y C +=4、积分因子的定义及判别对于微分形式的微分方程(,)(,)0M x y dx N x y dy += (2.43)如果方程(2.43)不是恰当方程,而存在连续可微的函数(,)0x y μμ=≠,使得 (,)(,)0M x y dx N x y dy μμ+= (2.51) 为一恰当方程,即存在函数(,)v x y ,使(,)(,)M x y dx N x y dy dv μμ+≡则称(,)x y μ是方程(2.43)的积分因子.此时(,)v x y C =是(2.51)的通解,因而也就是(2.43)的通解.如果函数(,),(,)M x y N x y 和(,)x y μ都是连续可微的,则由恰当方程的判别准则知道, (,)x y μ为(2.43)积分因子的充要条件是M Ny xμμ∂∂=∂∂ 即 ()M N N M x y y xμμμ∂∂∂∂-=-∂∂∂∂ (2.52) 5、积分因子的求法方程(2.52)的非零解总是存在的,但这是一个以μ为未知函数的一阶线性偏微分方程,求解很困难,我们只求某些特殊情形的积分因子. 定理2 设(,),(,)M M x y N N x y ==和(,)x y ϕϕ=在某区域内都是连续可微的,则方程(2.43) 有形如((,))x y μμϕ=的积分因子的充要条件是:函数(,)(,)(,)(,)(,)(,)y x x y M x y N x y N x y x y M x y x y ϕϕ-- (2.53)仅是(,)x y φ的函数,此外,如果(2.53)仅是(,)x y φ的函数((,))f f x y ϕ=,而()()G u f u du =⎰,则函数((,))G x y e ϕμ= (2.54) 就是方程(2.43)的积分因子.证明 因为如果方程(2.43)有积分因子()μμϕ=,则由(2.52)进一步知()()d M N N M d x y y xμϕϕμϕ∂∂∂∂-=-∂∂∂∂ 即y x x yM N d d N M μϕμϕϕ-=-由()μμϕ=可知左端是ϕ的函数,可见右端y x x yM N N M ϕϕ--也是ϕ的函数,即()y x x yM N f N M ϕϕϕ-=-,于是,有()d f d μϕϕμ=, 从而 ()()f d G e e ϕϕϕμ⎰==反之,如果(2.53)仅是ϕ的函数,即()y x x yM N f N M ϕϕϕ-=-,则函数(2.54)是方程(2.52)的解.事实上,因为()()()()G x y y x NM N M f e M N x yϕμμϕϕϕμ∂∂-=-=-∂∂ 因此函数(2.54)的确是方程(2.43)的积分因子.为了方便应用这个定理,我们就若干特殊情形列简表如下:例2. 解22(31)()0y xy dx xy x dy -++-= 解 这里2231,M y xy N xy x =-+=-,注意y x M N y x -=-所以方程不是恰当的,但是1y xM N Nx-=它仅是依赖与x ,因此有积分因子1dxx ex μ⎰≡=给方程两边乘以因子x μ=得到2223(3)()0xy x y x dx x y x dy -++-=从而可得到隐式通解22321122u x y x y x C ≡-+= 例3. 解方程2()(1)0xy y dx xy y dy ++++=解 这里2,1M xy y N xy y =+=++方程不是恰当的.但是1y xM N My-=-- 它有仅依赖于y 的积分因子 11dyy eyμ-⎰≡=方程两边乘以积分因子1y μ=得到 1()(1)0x y dx x dy y++++=从而可得到隐式通解21ln ||2u x xy y y C ≡+++= 另外,还有特解0y =.它是用积分因子乘方程时丢失的解. 例4. 解方程 223(2)()0y x y dx xy x dy +++=解 这里2232,M y x y N xy x =+=+,不是恰当方程.设想方程有积分因子()x y αβμμ=,其中α,β是待定实数.于是2112111()(2)y xM N y x x N y M x y y x x y x yαβαβαβαβαβαβ----⋅=⋅=--+- 只须取3,2αβ==.由上述简表知原方程有积分因子32x y μ=从而容易求得其通解为:446313u x y x y C ≡+=六、积分因子的其他求法以例4为例,方程的积分因子也可以这样来求:把原方程改写为如下两组和的形式:223()(2)0y dx xydx x ydx x dy +++=前一组有积分因子11yμ=,并且 21()()y dx xydy d xy y+= 后一组有积分因子21xμ=,并且 2321(2)()x ydx x dy d x y x+= 设想原方程有积分因子211()()xy x y y xαβμ==其中α,β是待定实数.容易看出只须3,2αβ==,上述函数确实是积分因子,其实就是上面找到一个.例5. 解方程 1212()()()()0M x M y dx N x N y dy += 其中1M ,2M ,1N ,2N 均为连续函数.解 这里12()()M M x M y =,12()()N N x N y =.写成微商形式就形式上方程是变量可分离方程,若有0y 使得20()0M y =,则0y y =是此方程的解;若有0x 使得10()0N x =,则0x x =是此方程的解;若21()()0M y N x ≠,则有积分因子211()()M y N x μ=并且通解为1212()()()()M x N y u dx dy N x M y ≡+⎰⎰ 例6、试用积分因子法解线性方程(2.28).解 将(2.28)改写为微分方程[()()]0P x y Q x dx dy +-= (2.55)这里()(),1M P x y Q x N =+=-,而()M Ny xP x N∂∂-∂∂=- 则线性方程只有与x 有关的积分因子()P x dxe μ-⎰=方程(2.55)两边乘以()P x dxe μ-⎰=,得()()()()()0P x dx P x dx P x dxxP x e ydx e dy Q x e dx ---⎰⎰⎰-+= (2.56) (2.56)为恰当方程,又分项分组法()()()()0P x dx P x dxd ye Q x e dx --⎰⎰-=因此方程的通解为()()()P x dx P x dxye Q x e dx c --⎰⎰-=⎰即()()[()]P x dxP x dxy e Q x e dx c -⎰⎰=+⎰与前面所求得的结果一样.注:积分因子一般不容易求得可以先从求特殊形状的积分因子开始,或者通过观察法进行“分项分组”法求得积分因子.§4 一阶隐方程与参数表示1、一阶隐方程一阶隐式微分方程的一般形式可表示为:(,,)0F x y y '=如果能解出(,)y f x y '=,则可化为显式形式,根据前面的知识求解. 例如方程2()()0y x y y xy ''-++=,可化为y x '=或y y '=但难以从方程中解出y ',或即使解出y ',而其形式比较复杂,则宜采用引进参数的方法求解.一般隐式方程分为以下四种类型:1) (,)y f x y '= 2) (,)x f y y '= 3) (,)0F x y '= 4)(,)0F y y '=2、求解方法Ⅰ)可以解出y (或)x 的方程1) 讨论形如(,)y f x y '= (2.57)的方程的解法,假设函数(,)f x y '有连续的偏导数,引进参数y p '=,则方程(2.57)变为 (,)y f x p = (2.58)将(2.58) 的两边对x 求导数,得到 f f dpp x y dx∂∂=+∂∂ (2.59) 方程(2.59)是关于,x p 的一阶微分方程,而且属于显式形式.若求得(2.59)的通解形式为(,)p x c ϕ=,将其代入(2.58),于是得到(2.57)通解为 (,(,))y f x x c ϕ=若求得(2.59)的通解形式为(,)x p c ψ=,于是得到(2.57)的参数形式的通解为(,)((,),)x p c y f p c p ψψ=⎧⎨=⎩其中p 为参数, c 是任意常数.若求得(2.59)的通解形式为(,,)0x p c Φ=,于是得到(2.57)的参数形式的通解为(,,)0(,)x p c y f x p Φ=⎧⎨=⎩其中p 为参数, c 是任意常数. 例1 求方程3()20dy dyx y dx dx+-= 的解 解 令dyp dx=,于是有32y p xp =+ (2.60) 两边对x 求导数,得到 2322dp dp p p x p dx dx=++ 即 2320p dp xdp pdx ++= 当0p ≠时,上式有积分因子p μ=,从而32320p dp xpdp p dx ++= 由此可知4234p xp c += 得到42223344c pc x p p p -==-将其代入(2.60),即得43342()c p y p p-=+故参数形式的通解为22334(0) 212c x p p p c y p p ⎧=-⎪⎪≠⎨⎪=-⎪⎩当0p =时,由(2.60)可知0y =也是方程的解.例2 求方程22()2dy dy x y xdx dx =-+的解. 解 令dy p dx =,得到222x y p xp =-+ (2.61) 两边对x 求导数,得到2dp dp p px p x dx dx =--+ 或 (2)(1)0dpp x dx--= 由10dpdx -=,解得p x c =+,于是得到方程的通解为222x y cx c =++ (2.62)由20p x -=,解得2xp =,于是得到方程的一个解为24x y = (2.63)特解(2.63)与通解(2.62)中的每一条积分曲线均相切,因此称为方程的奇解.2) 讨论形如(,)dyx f y dx= (2.64) 的方程的求解方法,方程(2.64)与方程(2.57)的求解方法完全类似,假定函数(,)f y y ' 有连续偏导数.引进参数dyp dx=,则(2.64) 变为 (,)x f y p = (2.65) 将(2.65) 的两边对y 求导数,得到1f f dpp y x dy∂∂=+∂∂ (2.66) 方程(2.66))是关于,y p 的一阶微分方程,而且属于显式形式.设其通解为 (,,)0y p c Φ= 则(2.64)的通解为(,,)0(,)y p c x f y p Φ=⎧⎨=⎩Ⅱ)不显含y (或)x 的方程 3) 讨论形如(,)0F x y '= (2.67) 的方程的解法.记dyp y dx'==,此时(,)0F x p =表示的是xp 平面上的一条曲线,设曲线用参数形式表示为()x t ϕ=,()p t ψ= (2.68) 由于dy pdx =,进而()()dy t t dt ψϕ'= 两边积分,得到()()y t t dt c ψϕ'=+⎰ 于是得到方程(2.67)参数形式的解为()()()x t y t t dt c ϕψϕ=⎧⎪⎨'=+⎪⎩⎰c 是任意常数.例3 求解方程3330x y xy ''+-=解 令y p tx '==,则由方程得 331t x t =+, 2331t p t =+ 于是 23339(12)(1)t t dy dt t -=+ 积分得到 23333329(12)314(1)2(1)t t t y dt c c t t -+=+=+++⎰ 故原方程参数形式的通解为: 3332313142(1)t x t t y c t ⎧=⎪+⎪⎨+⎪=+⎪+⎩4) 讨论形如(,)0F y y '= (2.69) 的方程,其解法与方程(2.67)的求解方法类似.记dy p y dx '==,此时(,)0F y p =表示的是yp 平面上的一条曲线,设曲线用参数形式表示为()y t ϕ=,()p t ψ= 由关系式dy pdx =可知 ()()t dt t dx ϕψ'=,于是0p ≠时,有 ()()t dx dt t ϕψ'=, ()()t x dt c t ϕψ'=+⎰ 故方程(2.69)的参数形式的通解 ()()()t x dt c t y t ϕψϕ'⎧=+⎪⎨⎪=⎩⎰c 是任意常数.此外,不难验证,若(,0)0F y =有实根y k =,则y k =也是方程的解. 例4 求解方程 22(1)(2)y y y ''-=-. 解 令2y yt '-=,则有222(1')y y y t -=由此可以得2'1y t =-,1y t t=+ 代入1dx dy p=,得到 222111(1)1dx dt dt t t t=-+=-- 积分,得到1x c t=+ 故原方程参数形式的通解为 11x c t y t t ⎧=+⎪⎪⎨⎪=+⎪⎩其中c 是任意常数.此外, 当0y '=时原方程变为24y =,于是2y =±也是方程的解. 例5求解方程y '=解 令y p '=,则有p =,取,(,)22p tgt t ππ=∈-,则sin sec tgt x t t === 由dy pdx =得到cos sin dy tgt tdt tdt == 所以cos y t c =-+故原方程参数形式的通解为 sin cos x t y t c =⎧⎨=-+⎩其中c是任意常数.。
常微分方程第三版课本 (2)
常微分方程第三版课本概述“常微分方程第三版课本”是一本由X编写的教材,主要介绍了常微分方程的基本概念、理论和解析方法。
本教材内容丰富、结构清晰,适用于高等院校的常微分方程课程教学,也可以作为自学的参考资料。
目录1.基本概念– 1.1 常微分方程的定义– 1.2 解的定义及存在唯一性定理– 1.3 初值问题和边值问题2.一阶常微分方程– 2.1 可分离变量方程– 2.2 齐次线性方程– 2.3 一阶线性常微分方程– 2.4 Bernoulli 方程和 Riccati 方程– 2.5 可降阶的高阶微分方程3.高阶线性常微分方程– 3.1 高阶常微分方程的一般理论– 3.2 同解、通解和特解– 3.3 常系数齐次线性方程– 3.4 常系数非齐次线性方程及其特解– 3.5 变系数线性方程4.线性常微分方程组– 4.1 二阶齐次线性方程组– 4.2 二阶非齐次线性方程组和线性常系数方程组– 4.3 三阶及三阶以上线性方程组内容简介基本概念本章介绍了常微分方程的基本概念,包括常微分方程的定义、解的定义及存在唯一性定理、初值问题和边值问题。
通过对这些概念的学习,读者可以对常微分方程有一个基本的认识。
一阶常微分方程本章主要介绍了一阶常微分方程的解析方法,包括可分离变量方程、齐次线性方程、一阶线性常微分方程、Bernoulli方程和 Riccati 方程、可降阶的高阶微分方程等。
通过对这些解析方法的学习,读者可以熟练地解决一阶常微分方程的问题。
高阶线性常微分方程本章主要介绍了高阶线性常微分方程的理论和方法。
包括高阶常微分方程的一般理论、同解、通解和特解、常系数齐次线性方程、常系数非齐次线性方程及其特解、变系数线性方程等。
通过对这些理论和方法的学习,读者可以掌握高阶线性常微分方程的解法。
线性常微分方程组本章主要介绍了线性常微分方程组的解法。
包括二阶齐次线性方程组、二阶非齐次线性方程组和线性常系数方程组、三阶及三阶以上线性方程组等。
常微分方程教案(王高雄)第二章
第二章目录内容提要及其它 (1)第二章一阶微分方程的初等解法(初等积分) (2)第一节变量分离方程与变量变换 (2)一、变量分离方程 (2)二、可化为变量分离方程的类型 (6)1、齐次方程 (6)2、可化为变量分离方程 (7)三、应用例题选讲 (10)第二节线性方程与常数变易法 (11)第三节恰当方程与积分因子 (15)一、恰当方程 (15)二、积分因子 (20)第四节一阶隐含方程与参数表示 (23)一、可以解出y(或x)的方程 (24)二、不显含y(或x)的方程 (25)本章小结及其它 (27)内容提要及其它授课题目(章、节)第二章:一阶微分方程的初等解法教材及主要参考书(注明页数)教材:常微分方程(第三版),王高雄等,高等教育出版社,2006年,p30-74主要参考书:[1]常微分方程,东北师范大学微分方程教研室编,高等教育出版社,2005,p1-70[2]常微分方程教程,丁同仁等编,高等教育出版社,1991,p1-20[3]偏微分方程数值解法(第2版),陆金甫关治,清华大学出版社,2004,p1-12[4]常微分方程习题解,庄万主编,山东科学技术出版社,2003,p28-169[5]微分方程模型与混沌,王树禾编著,中国科学技术大学出版社,1999,p15-158[6]差分方程和常微分方程,阮炯编著,复旦大学出版社,2002,p38-124目的与要求:掌握变量分离方程、齐次方程、线性方程、伯努利方程和恰当方程的解法.理解变量变换思想方法和积分因子方法,并能应用于求解一些特殊的常微分方程.掌握四类典型的一阶隐方程的解法.能熟练求解变量分离方程、齐次方程、线性方程、伯努利方程、恰当方程和四类典型的一阶隐方程.领会变量变换思想方法和积分因子方法,并能应用于求解一些特殊的常微分方程.教学内容与时间安排、教学方法、教学手段:教学内容:第1节变量分离方程与变量变换;第2节线性方程与常数变易法;第3节恰当方程与积分因子;第4节一阶隐方程与参数表示:可以解出(或y x)的方程、不显含(或y x)的方程.时间安排:8学时教学方法:讲解方法教学手段:传统教学方法与多媒体教学相结合。
常微分方程(王素云)-第2章
a 1 tan(gat arctan av0 ) 0,
解得
t g a 1 1 arctan av0
第2章 初等积分法
例2.2 解方程
dy
1 y2 .
d x 1 x2
解 对式(2.7)分离变量,得
dy dx 1 y2 1 x2
对上式两端积分后便得隐式通解:
(2.7)
arcsin y arcsinx C
阻力与质点运动速度的平方成正比(比例系数为k>0),求该质 点从抛出至达到最高点的时间.
解 设质点在t时刻的速度为v,v(0)=v0,且有
m d v mg kv2 dt
解得
v a1 tan(gat c), a k / mg
因为v(0)=v0,所以解出c=arctanav0.
第2章 初等积分法 质点达到最高点,即v=0,亦即
第2章 初等积分法 为得出显式通解可从上式解出y :
y sin(arcsinx C) x cosC sin C cos(arcsinx) x cosC sin C 1 x2 x 1 C12 C1 1 x2
(2.8)
其中,常数C1=sinC.容易看出,y=±1也是方程(2.7)的解, 但它们不包含在通解式(2.8)中.注意,这是由于y=±1时式 (2.7)的分子等于零的缘故.
R : a x b, c y d
内的任意一点.为了求方程(2.1)所满足的初值条件y(x0)=y0的 解,可按
y0 G 1 (H (x0 ) C)
以确定常数C,即
C G( y0 ) H (x0 )
代入式(2.5)便得
y G 1 (H (x) G( y0 ) H (x0 ))
(2.6)
现在来分析方程(2.1)的积分曲线,即解在(x,y)平面上 的图像的分布情况.如果g(y)≠0(c<y<d),则由于初值问题的 解存在而且唯一,经过域R内每一点都恰有方程(2.1)的一条 积分曲线,积分曲线在R内彼此不相交.这时,积分曲线的分布 情况很简单.但是,如果g(y)为c<y<d上的某些点,比如 g(y)=g(y1)=0,在这种情况下,方程(2.1)经过直线y=y1上的 点的积分曲线就很可能不止一条
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
常微分方程讲义(一)
课程目标:
掌握常用的常微分方程解题技巧;利用常微分方程的思想建模。
上课方式:
课堂讲授、练习与考试。
课程特点:
承接高数、微积分、数学分析等课程而来,与导数、积分的关系非常紧密,在经济数学中有广泛的应用;常与其他数学工具与方法混合使用。
参考书目:
《常微分方程》,蔡燧林编著,武汉大学出版社,2003;及所有标注有“常微分方程”、“应用”、“经济数学”、“金融数学”的教材与专著。
为什么在模拟经济变化时要引入常微分方程?
注重刻画在无穷小时间段内的变量的动态变化,实现了从“静态”向“动态”的飞跃。
微分方程比初等函数更近于现实,更真于模拟。
什么是方程?)(x
y 。
f
什么是微分方程?
dy的方程;
常微分方程:含有dy、dx、
dx
偏微分方程:含有y ∂、x ∂、x
y ∂∂的方程。
x
y ∆∆的几何含义:割线、割线的斜率 dx
dy 的几何含义:切线、切线的斜率 dx
dy x y x =∆∆→∆0lim :数学上——切线的斜率,导数 经济上——变化率,边际
例:求2x y =与x e y =的导数
应当记下来的等式:
1)'(-=n n nx x ,c x dx nx n n +=⎰-1
x x e e =)'(,c e dx e x x +=⎰
x x 1
)'(ln =,C x dx x +=⎰ln 1
x x cos )'(sin =,⎰+=C x xdx sin cos
x x sin )'(cos -=,⎰+=-C x dx x cos )sin (
x tgx 2sec )'(=,⎰+=C tgx xdx 2sec
x ctgx 2csc )'(-=,⎰+=-C
ctgx dx x )csc (2
0)'(=C
k kx =)'(
'')'(b a b a +=±
'')'(ab b a ab +=
2'
')'(b ab
b a b a -=
'')])'([(g f x g f =
)()')((x p dx x p =⎰
x x x 2121)'(21==-,⎰+=C x dx x 21
a a a x x ln )'(=,C a adx a x x +=⎰ln
211
)'(arcsin x x -=,⎰+=-C x dx x arcsin 112
211
)'(x arctgx +=,⎰+=+C arctgx dx x
211
例:匀速运动与变速运动
例:不良资产的处置
常微分方程的“阶” 考察方程中导数的最高“阶”n n dx
y d dx y d dx dy ......,22, 而不是考察方程中的最高“次方”n dx
dy dx dy dx dy )......()(,2
常微分方程的“解”
通解:曲线族
特解:初值条件
例:检验1121++=
C x C y 是方程0)'(12''2=-+y y
y 的解
例:检验C y y x =+22sin 是方程0')22sin (sin 2=++y y y x y 的解
例:检验由参数方程⎪⎩
⎪⎨⎧+-=+-=C t t y t t x 24321432所决定的函数)(x f y =,是微分方程2)(3+-=dx
dy dx dy x 的解
例:设)(x p 是区间(a ,b )上的连续函数,证明⎰=-dx x p Ce y )(是微分方程0)('=+y x p y 在区间+∞<<<y b x a ,内的解。
例:一曲线经过点(2,0),且其上任意一点的切线界于切点和纵坐标轴之间的部分的长度恒等于2,求此曲线所满足的微分方程的表达式。
文献清单:
● 《常微分方程及其应用》,周义仓等编,科学出版社,2004
● 《微分方程模型》,(美)William F.Lucas 主编,朱煜民等译,国防科技大学,1998
● 《微分方程模型与混沌》,王树禾编著,中国科学技术大学出版社,1999
●《微分方程及其应用》,(美)M.Braun著,张鸿林译,人民教育出
版社,1980
●《常微分方程习题集》,周尚仁等编,人民教育出版社,1980
●《经济应用数学》,万世栋等主编,科学出版社,2002
●《经管财金建模方法及应用》,饶友玲等编著,清华大学出版社,
2004。