常微分方程第三版课件2.1.ppt

合集下载

第五讲常微分方程PPT课件

5. 求lim x0
1 cos x
.
1
6.
求
lim
xe
x e
xe
.
7.
设
y
x2
sin
1 x
,
x 0,
存在. 0,
x 0,
求y 0
8. 计算积分
x3 dx.
1 x2
并讨l论im y x x0
是否
第37页/共47页
综合练习
9. 计算下列积分.
1
arctan x
x dx;
2
ln x 1 x2 dx.
任给有理数a,
函数
f(x)满足 f
x
x
0
f
a t dt 1,
求
f(x).
练 (2008年高数二)
求微分方程
d2y dx 2
dy dx
0
的通解.
第26页/共47页
3.掌握二阶常系数非齐次线性微分方程的解法二阶常系数非齐次线性微分方程：
ay by cy f x
的通解为
y Y x y* x
y 4 y 0 的通解.
例: 求齐次方程
4
d2x dt 2
20
dx dt
25 x
0
的通解.
例: 求初值问题
y 4 y 29 y 0
y
x0
0
,
y
x0
15
的解.
第25页/共47页
练 (2006年高数二)
微分方程
y 4 y 5 y 0 的通解为___________
练 (2007年高数一)
第16页/共47页
二阶齐次线性方程解的结构

常微分方程第二章第一讲

2.1.2 可化为变量分离方程的类型
引言有的微分方程从表面上看，不是可分离变量的微分方程，但是，通过适当的变量替换，就可以很容易地化为“变量分离方程”，在这里，介绍两类这样的方程。 1. 齐次方程
1)方程的类型
定义
dy y g ( ) (2.5) 的方程，称为齐次 dx x 微分方程，这里 g (u ) 是 u 的连续函数。 14
dy ( y) f ( x)dx C (2.2)
可以证明这就是方程(2.1)的通解.
2)如果存在 y0， ( y0 ) 0, 则方程（ .1 使 2 ）还有特解
y y0
(**)
微分方程(2.1)的所有解为：式(2.2)和(**).
注意：积分常数C 的相对任意性。
7
3.变量分离方程的解题步骤
即 1 , 2 1 ,
则 ON OM ,
PM 而 tan 2 , OP ON
_____ _____
则有 y'
y x x y
2 2
.
上述方程为齐次微分方程，可用变量变换法求解。
27
小结 1.变量分离方程的形状 dy f ( x) ( y )或M 1 ( x) N1 ( y ) dx M 2 ( x) N 2 ( y ) dy 0 dx 2.变量分离方程的求解：分离变量法步骤：分离变量，两边积分，检查是否有遗漏的特解
2
(*)
23
分离变量，得 dX 1 u du 2 X 1 2u u 两边积分，得 ~ 2 2 ln X ln | u 2u 1 | C
即X (u 2u 1) C1 (C1 e ), 此外容易验证 u 2 2u 1 0 亦为方程（*）的解，因此方程(*)的通解为 X 2 (u 2 2u 1) C1，其中C1为任意常数。

常微分方程(王高雄)第三版 2.1教学教材

（I）齐次方程
ddyxg(yx)
(II) 形如 ddyxfaa21xxbb12yycc12的方,程其中 a1,b1,c1,a2,b2,c2为任意.常数
(I) 形如
dyg(y) dx x
(2.5)
方程称为齐次方程, 这里g(u)是u的连续函. 数
求解方法: 10 作变量代换(引入新变量)u y ,方程化为
x
du g(u)u, (这里d由 yx于 duu)
dx x
dx dx
20 解以上的变量分离方程
30 变量还原.
例4 求解方程 xdy 2xyy dx
(x0)
解: 方程变形为 dy2 yy dx x x
(x0)
这是齐次方程, 令u y 代入得 x
x du u 2 uu 即 x du 2 u
dx
为 (1)的情形,可化为变量分离方程求解.
解的步骤:
10解方程 aa21xx 组 bb1 2yy cc1200,
得解yx
,
20 作变换 YXyx,方程化为
dY a1Xb1Y dX a2Xb2Y
g
(
Y X
)
30再经变 u换 Y,将以上方程化离为方变程量分
X
40 求解
50 变量还原
dx
10 分离变量, 当 (y)0时 ,将 (2.1)写成
dy f (x)dx,
(y)
这样变量就“分离”开了.
20 两边积分得
dy
(y)f(x)d xc (2.2)
1 的某一原函数 f (x)的某一原函数 ( y)
由 (2.2)所确定 y的 (x,c)就函 (2 为 .数 1)的.解
例：
分离变量:

常微分方程----第一章-绪论

目录上页下页返回结束
莱布尼兹(1646 – 1716)
德国数学家, 哲学家. 他和牛顿同为微积分的创始人 , 他在《学艺》杂志上发表的几篇有关微积分学的论文中, 有的早于牛顿, 所用微积分符号也远远优于牛顿 . 他还设计了作乘法的计算机 , 系统地阐述二进制计数法 , 并把它与中国的八卦联系起来 .
c2
cn
则称 y (x,c1,,cn ) 含有n个相互独立的常数。
目录上页下页返回结束
例：y c1 cos x c2 sin x 是 y y 0 的通解。因为 y c1 sin x c2 cos x 而
cos x sin x 1 0
sin x cos x
内容小结
1. 微分方程的基本概念常微分方程，偏微分方程，微分方程的阶
微分方程的解，通解，特解
线性微分方程，非线性微分方程初始条件
作业
P27 2， 3，4， 6，8 （1）（3）（5）
目录上页下页返回结束
牛顿(1642 – 1727)
伟大的英国数学家 , 物理学家, 天文学家和自然科学家. 他在数学上的卓越贡献是创立了微积分. 1665年他提出正流数 (微分) 术 , 次年又提出反流数(积分)术,并于1671 年完成《流数术与无穷级数》一书 (1736年出版). 他还著有《自然哲学的数学原理》和《广义算术》等 .
牛顿在建立微积分的同时，对简单的微分方程用级数来求解。后来瑞士数学家雅各布·贝努利、欧拉、法国数学家克雷洛、达朗贝尔、拉格朗日等人又不断地研究和丰富了微分方程的理论。
常微分方程的形成与发展是和力学、天文学、物理学，以及其他科学技术的发展密切相关的。同时，数学的其他分支的新发展，如复变函数、李群、组合拓扑学等，都对常微分方程的发展产生了深刻的影响，当前计算机的发展更是为常微分方程的应用及理论研究提供了非常有力的工具。

《常微分方程》全套课件(完整版)

捕捉到这种联系，而这种联系，用数学语言表达出来，其结果往往形成一个微分方程.一旦求出这个方程的解，其运动规律将一目了然.下面的例子，将会使你看到微分方程是表达自然规律的一种最为自然的数学语言.
例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 =±１，这两个解不包含在上述解中.
3. 函数
是方程(1.6)在区间(-∞，
+∞)上的解，其中和是独立的任意常数.
4. 函数
是方程(１.７)在区间(-
∞，+∞)上的解，其中和是独立的任意常数.
这里，我们仅验证3，其余留给读者完成.事实上，
（1.13）
显然，方程(1.4)是一阶线性方程；方程(1.5)是一阶非线性方程；方程 (1.6)是二阶线性方程；方程(1.7)是二阶非线性方程.
通解与特解
微分方程的解就是满足方程的函数，可定义如下.
定义1.１设函数在区间I上连续，且有直
到n阶的导数.如果把
代入方程(1.11)，得到在
区间I上关于x的恒等式，

高等数学常微分方程PPT课件

第12页/共35页
【解法】需经过变量代换化为一阶线性微分方程.
除方程两边 , 得
yn d y P( x) y1n Q( x) dx
令 z y1n , 则 dz (1 n) yn d y
dx
dx
dz (1 n) P( x) z (1 n)Q( x) (关于z , x的一阶线性方程) dx
特征方程法
待定
特征方程的根及其对应项
系
数
法 f(x)的形式及其
特解形式
高阶方程可降阶方程
线性方程解的结构
定理1;定理2 定理3;定理4
欧拉方程
第4页/共35页
微分方程解题思路
一阶方程
作变换
降阶
高阶方程
分离变量法全微分方程常数变易法
作变换积分因子
非非变全量微可分
分方离程
特征方程法
［提示］(1)
原方程化为
令u=xy,得 (2) 将方程改写为
d u u ln u (分离变量方程) dx x
d y 1 y y3 (贝努里方程) d x 2x ln x 2x
令 z y2
第17页/共35页
【例3】识别下列一阶微分方程的类型，并求解
1)
【解】
y y x
①可分离变量的微分方程
u e P( x)d x P( x) ue P( x)d x P( x) u e P( x)d x Q( x)
即两端积分得
非齐பைடு நூலகம்方程
dy P(x) y Q(x)
dx
u Q(
对应齐次方程通解
x
)
e
P( x)d
y
x
dx

《常微分方程》(第三版)

常微分方程2.11.xy dxdy2=,并求满足初始条件：x=0,y=1的特解. 解：对原式进行变量分离得。

故它的特解为代入得把即两边同时积分得：e e xx y c y x x c y c y xdx dy y22,11,0,ln ,212=====+==,0)1(.22=++dy x dx y 并求满足初始条件：x=0,y=1的特解.解：对原式进行变量分离得：。

故特解是时，代入式子得。

当时显然也是原方程的解当即时，两边同时积分得；当xy c y x y x c y c y x y dy dx x y++=====++=+=+≠=+-1ln 11,11,001ln 1,11ln 0,11123 yxy dx dyxy 321++=解：原式可化为：x x y xx y x yx yyxyc c c c x dx x dy y yx ydxdy 2222222232232)1(1)1)(1(),0(ln 1ln 21ln 1ln 2111,0111=++=++≠++-=++=+≠+∙+=+）故原方程的解为（即两边积分得故分离变量得显然.0;0;ln ,ln ,ln ln 0110000)1()1(4===-==-+=-++=-=+≠===-++x y c y x xy c y x xy c y y x x dy yydx x x xy x y xdy y ydx x 故原方程的解为即两边积分时，变量分离是方程的解，当或解：由：10ln 1ln ln 1ln 1,0ln 0)ln (ln :931:8.cos ln sin ln 07ln sgn arcsinln sgn arcsin 1sgn 11,)1(,,,6ln )1ln(21111,11,,,0)()(:53322222222222c dx dy dx dy xycy ud uu dx x x y u dx xydy x y ydx dy y x x c dy yy yydx dy c x y tgxdx ctgydy ctgxdy tgydx c x x xyc x x u dx xx du xdxdu dxdux u dx dy ux y u x y y dx dy xc x arctgu dxx du u u u dx du x u dxdu xu dx dy ux y u x y x y x y dx dy dx x y dy x y e e e e e e ee x y uu xy x u u x yxyy x xx+===+=+-===-∙-=--+-=-=+-===-=+∙=+∙=∙=--=+===-+=+-=++=++-++=++===+-==-++-+--两边积分解：变量分离：。

常微分方程的几何解释

(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),的解是存在,是否唯一?

常微分方程课件_第01讲(66)

d 2 d g 1 F (t ). 2 dt m dt l ml
(1.11)
20
When we want to determine some particular motion, we should give its initial states(see C/P4): For example, when t=0,
参考书：中文/丁同仁、李承治（北京大学）编，常微分方程教程，第二版，高等教育出版社，2004.7
英文/M. Braun, Differential Equations and Their Applications, Springer-Verlag为：平时（以到课情况、完成作业情况）与期末考试． 2. 考核形式：期末考试采用闭卷考试形式，试题来自试题库． 3. 期末考试题形及分值比例为：填空题(4个)……………………………约占20%，选择题(4个)……………………………约占20%，计算题(5个)……………………………约占44%，证明与应用题(2个)……………………约占16%．
d 2 g 0. 2 dt l
(1.9)
18
(2) Free damped motion(有阻尼的自由运动E/P136): with damping but no external force(有阻尼但无外力) μ---the damping coefficient(阻尼系数) Math. Model(see C/P4):
Solution Function of the curve: y=f(x) Math. Model: Method:
dy 2 x, dx y ( 2) 5.
12
Integral (积分)

常微分方程2.1 线性方程

仍为非齐次方程的解。 4.非齐次方程的两解之差为对应齐次方程的解。 5.非齐次方程的任一解与对应齐次方程的齐次方程
的通解之和是非齐次方程的通解。
11
y e p( x)dx[ g( x)e p( x)dxdx C]
求方程 y 1 y sin x 的通解. p( x) ? g( x) ?
x
x
解
y
e
1 dx x
[
sin x x
e
1 dx
x dx
C
]
cos x
x
C x
考虑：dy dx , yx
对应齐通解：y C , x
设非齐通解：y u( x) , x
u( x) sin x ,
x
x
u( x) cos x C.
12
解初值问题:
( x2 1) y 2xy cos x 0
1
本章的主要内容
2.1 线性方程 2.2 变量可分离方程 2.3 全微分方程 2.4 变量替换法
2.5 一阶隐式方程 2.6 近似解法 2.7 一阶微分方程1 线性方程
一阶线性微分方程 y ' p(x) y g(x)
一、线性齐次方程
若 y ' p(x) y g(x) 中 g(x) 0 时，
先解：dx dy , x y ln y
ln x ln ln y lnC,
设 : x u( y) , u( y) 1 , u( y) 1 ln2 y C,
ln y ln y y
2
x 1 ln y C .
2
ln y
此外, y = 1 也是原方程的解.
16
解微分方程 dy sin y x cos y x 0 dx

常微分方程总结 PPT

2. 可分离变量方程的求解方法: 分离变量后积分; 根据定解条件定常数 . y 3 .齐次方程的求解方法: 令 u , x
8
机动目录上页下页返回结束
3. 解微分方程应用题的方法和步骤
(1) 找出事物的共性及可贯穿于全过程的规律列方程. 常用的方法: 1) 根据几何关系列方程 ( 如: P263,5(2) ) 2) 根据物理规律列方程 ( 如: 例4 , 例 5 )
线性无关概念.
23
机动目录上页下页返回结束
定义: 设 y1 ( x), y2 ( x),, yn ( x) 是定义在区间 I 上的
n 个函数, 若存在不全为 0 的常数使得
则称这 n个函数在 I 上线性相关, 否则称为线性无关. 例如，在( , )上都有
故它们在任何区间 I 上都线性相关;
微分方程的基本概念
含未知函数及其导数的方程叫做微分方程 . 常微分方程 (本章内容)
分类
偏微分方程
方程中所含未知函数导数的最高阶数叫做微分方程
的阶. 一般地 , n 阶常微分方程的形式是
F ( x, y, y,, y ( n ) ) 0
或
y ( n ) f ( x, y, y,, y ( n 1) ) ( n 阶显式微分方程)
y p( x) y q( x) y f ( x) ,
y
( n) ( n 1)
为二阶线性微分方程.
n 阶线性微分方程的一般形式为
a1 ( x) y an 1 ( x) y an ( x) y f ( x) f ( x) 0 时, 称为非齐次方程 ;
f ( x) 0 时, 称为齐次方程.
若 Q(x) 0, 称为非齐次方程 . dy P( x) y 0 1. 解齐次方程 dx

常微分方程_高教出版社_第三版

第一章绪论微分方程: 联系自变量、未知函数以及它的导数间的关系式。

自变量只有一个的称为常微分方程.§1.1 常微分方程模型例1 RLC 电路包括电感L , 电阻R 和电容C 及电源的电路称为RLC电路. 电流I 流经R,L,C 的电压降分别是RI, d d I Lt , Q C,其中Q 为电量, 它与电流的关系为d d QI t =. 基尔霍夫第二定律: 闭合回路中, 所有支路上的电压的代数和为零.如图所示的 RL 电路, 电感L , 电阻R 和电源电压E 为常数.设0t =时, 电路中没有电流. 开关S 合上后电流应满足的微分方程d 0d I E L RI t--=, 即d d I R E I t L L +=, 求出的()I t 应满足: 0t =时, 0I =.如果在0t t =时, 0I I =, 电源E 突然短路, 则E 变为0并且此后一直保持为0, 则电流I 满足方程d 0d I RI t L +=, 及条件0t t =时, 0I I =.再看如图所示的RLC 电路, 电阻R, 电感L 和电容C 都是常数. 电源()e t 是时间t 的已知函数.开关S 合上后, 电流I 应满足的微分方程()d d I Qe t L RI t C=++, 微分上式可得()22d d d 1d d d e t I R I I t L t LC L t++=, 如果()e t =常数, 则有22d d 10d d I R I t L t LC++=. 如果电阻R =0, 则有22d 10d I t LC+= 例2 数学摆解设摆在铅垂线右边时所成夹角ϕ为正. 质点M 沿圆周切向速度v 可表示为d d v l tϕ=. 重力mg 沿圆周切向的分力为MP, 数值为sin mg ϕ-, 于是摆的运动方程为 d sin d vmmg tϕ=-, 即22d sin d gt l ϕϕ=-. 如果是微小振动, 即ϕ比较小时, 可取sin ϕϕ≈, 于是微小振动方程为22d 0d gt lϕϕ+=. 如果摆在一个粘性介质中运动, 设阻力系数为μ, 则摆的运动方程为22d d 0d d gt m t lϕμϕϕ++=. 如果沿摆的运动方向恒有一个外力()F t 作用于它, 则称受迫微小振动, 方程为()22d d 1d d g F t t m t l mlϕμϕϕ++=.摆的初始条件为0t =时, 0ϕϕ=,0d d tϕω=. 例3 人口模型Malthus 假定: 人口出生率是常数r , 则从t 到t t +∆这段时间人口数量()N t 的增长量为()()()N t t N t rN t t +∆-=∆于是人口数量满足d d NrN t = 改写为d d Nr t N= 两边积分可得ln N rt c =+ 这里c 为任意常数, 上式又可变形为rt N ce =这里c c e =, 注意0N =也是解, 所以c 可以是任意常数. 如果设初值条件为0t t =时, ()0N t N =代入上式可得00rt c N e -=, 即方程满足此初始条件的解为()()00r t t N t N e -=.Logistic 模型: 引入环境最大容纳量m N , 假定净相对增长率为()1m N t r N ⎛⎫- ⎪⎝⎭, 则人口模型变为d 1d m N N r N t N ⎛⎫=- ⎪⎝⎭. 例4 传染病模型设某地区在某种传染病传播期间总人数保持不变, 为常数n . 开始感染人数为0x , 在t 时刻的健康人数为()y t ,染病人数为()x t , 则有()()x t y t n +=设单位时间内一个病人能传染的人数和当时健康人数成正比, 比例常数为k , 称之为传染系数, 于是()()()d d x t ky t x t t= 注意到总人数不变, 可得()()0d ,0d xkx n x x x t=-= 此模型称为SI 模型, 即Susceptible, Infective.对无免疫性的疾病, 病人治愈后会再次感染. 设单位时间治愈率为μ, 则SI 模型应修正为 ()()()()()0d ,0d x t ky t x t x t x x tμ=-=, 即()()0d 1,0d x kx n x x kx n x x x t μσ⎛⎫=--=--= ⎪⎝⎭这个称为SIS 模型. 其中1μ是这个传染病的平均传染期,kσμ=是整个传染期内每个病人有效接触的平均人数(接触数).对于免疫性很强的疾病, 病人治愈后不会再被感染, 即在t 时刻的治愈后免疫人数为()r t , 称为移除者(Removed), 设治愈率l 为常数, 即()()d d r t lx t t= 注意到总人口不变, ()()()x t y t r t n ++=, 我们得到d d d d xkxy lx ty kxy t⎧=-⎪⎪⎨⎪=-⎪⎩这个模型称为SIR 模型.例5 两生物种群生态模型某环境中有两种鱼: 被食鱼与捕食鱼. 设t 时刻被食鱼的总数为()x t , 捕食鱼的总数为()y t , 如果没用捕食鱼, 则被食鱼的增长规律为d d xax t=, 设捕食率为b , 则有d d xax bxy t=- 而捕食鱼有一个自然减少率c , 被食鱼供养捕食鱼的能力为d , 则有d d ycy dxy t=-+ 这个称之为Volterra 捕食-被捕食模型.其更一般的模型为()()d d d d d xx a bx cy ty y ex fy t⎧=++⎪⎪⎨⎪=++⎪⎩ 从数学的角度归类:d d I R E I t L L +=, d d N rN t =可以写为d d yay c t+=. 而 ()22d d d 11d d d e t I R I t L t LC L t++= 和()22d d 1d d g F t t m t l ml ϕμϕϕ++= 可以写为()22d d d d y yb cy f t t t++=.§1.2 基本概念和常微分方程的发展历史1.2.1 常微分方程的基本概念(1) 常微分方程和偏微分方程如果在微分方程中自变量的个数只有一个, 则称为常微分方程; 自变量的个数多于一个的微分方程则称为偏微分方程.第一节中的例子都是常微分方程. 以下是偏微分方程2222220T T Tx y z ∂∂∂++=∂∂∂,224T Tx t∂∂=∂∂. 阶数: 微分方程中出现的最高阶导数的阶数称为微分方程的阶数. 一般的n 阶常微分方程具有如下形式:d d ,,,,0d d n n y y F x y x x ⎛⎫= ⎪⎝⎭ ,这里F 是d d ,,,,d d n n y y x y x x 的表达式, 且必含有d d n n yx, y 是未知函数, x 是自变量.此书中常微分方程也简称为微分方程或方程. (2) 线性和非线性如果方程d d ,,,,0d d n n y y F x y x x ⎛⎫= ⎪⎝⎭ 左端为y 及d d ,,d d n n y yx x的一次有理整式, 则称方程为n 阶线性微分方程. 一般n 阶线性微分方程的形式为()()()()1111d d d d d d n n n n n n y y ya x a x a x y f x x x x---++++= , 这里()()()1,,,n a x a x f x 是x 的已知函数.不是线性方程的方程统称为非线性方程. 例如22d sin d gt lϕϕ=- 是二阶非线性方程.(3) 解和隐式解如果函数()y x ϕ=代入方程d d ,,,,0d d n n y y F x y x x ⎛⎫= ⎪⎝⎭ 后能使它变为恒等式, 则称()y x ϕ=为方程的解. 如果关系式(),0x y Φ=决定的隐函数()y x ϕ=是方程的解, 则称(),0x y Φ=为方程的隐式解.例: 一阶微分方程d d y xx y=-的有解y =y =则关系式221x y +=就是此方程的隐式解.解和隐式解统称为方程的解而不加以区别. (4) 通解和特解含有n 个独立的任意常数12,,,n c c c 的解()12,,,,n y x c c c ϕ=称为方程d d ,,,,0d d n n y y F x y x x ⎛⎫= ⎪⎝⎭ 的通解. 同样可定义隐式通解. 它们统称为方程的通解而不加以区分. 为了确定微分方程一个特解所需的条件称为定解条件. 常见的定解条件是初始条件, 方程d d ,,,,0d d n n y y F x y x x ⎛⎫= ⎪⎝⎭ 的初始条件是指当0x x =时, ()()1110001d d ,,,d d n n n y y y y y y x x---=== ,()()110000,,,,n x y y y - 是给定的n +1个常数.求微分方程满足定解条件的解, 就是定解问题. 当定解条件为初始条件时, 称为初值问题, 这也是本书讨论的主要内容.满足初始条件的解称为微分方程的特解. 初始条件不同, 特解也不同.例: 人口模型的解rt N ce =含有一个任意常数c , 所以是d d NrN t=的通解, 而()()00r t t N t N e -=就是满足初始条件0t t =时, ()0N t N =的特解, 它可在通解中令00rt c N e-=得到.例: 二阶微分方程22d d 540d d y yy x x++= 的通解为412x x y c e c e --=+,这里12,c c 是任意常数, 满足初始条件()()d 002,1d y y x== 的特解为43x x y e e --=-. 5) 积分曲线和方向场一阶微分方程()d ,d yf x y x= 的解()y x ϕ=代表xy 平面上一条曲线, 称为微分方程的积分曲线. 而通解(),y x c ϕ=则对应一族曲线, 称为积分曲线族.满足初始条件()00y x y =的解就是过点()00,x y 的积分曲线. 积分曲线上每一点(),x y 的切线斜率正好就是(),f x y . 反之, 如果有某条曲线, 它在点(),x y 的切线斜率是(),f x y , 则它就是一条积分曲线.方向场: 设(),f x y 的定义域为D , 在每个(),x y D ∈上画一个方向, 此方向的斜率等于(),f x y , 这种带有方向的区域称为方程()d ,d yf x y x=确定的方向场. 等斜线: 在方向场中, 方向相同的点的轨迹称为等斜线. ()d ,d yf x y x =的等斜线方程为(),f x y k =,其中k 是参数.例 d 1d yxy x=+. 利用Maple 模拟出的此方程的方向场:6) 微分方程组用两个及两个以上的关系式表示的涉及多个函数的导数的微分方程称为微分方程组.第二章一阶微分方程的初等解法初等解法: 即将微分方程的求解问题转化为积分问题注1 不一定要求用初等函数表示积分.注2 并不是所有的微分方程都有初等解法. §2.1 变量分离方程与变量变换2.1.1 变量分离方程形如()()d d y f x y x ϕ= 的方程, 称为变量分离方程, 这里()f x 和()y ϕ分别是x,y 的连续函数.解法: 如果()0y ϕ≠, 则方程可改写为()()d d y f x x y ϕ=, 两边积分, 得到()()d d y f x x c y ϕ=+⎰⎰,其中c 是任意常数, 而()d y y ϕ⎰和()d f x x ⎰则分别表示()1y ϕ和()f x 的一个原函数. 容易验证方程()()d d y f x x c y ϕ=+⎰⎰所确定的隐函数(),y y x c =就是原微分方程的通解.如果存在0y 使得()00y ϕ=, 则0y y =也是原方程的解, 它不包含的通解中, 须补上.例1 求解方程d d y x x y=-. 解分离变量, 可得d d y y x x =-,两边积分 2222y x c =-+, 化简可得通解为22x y c +=.例2 求解两种群模型 ()()d d d y c x y x x a by -+=-, 0,0x y ≥≥. 解分离变量 d d d c a x b y x y ⎛⎫⎛⎫-=-+ ⎪ ⎪⎝⎭⎝⎭, 积分得 ln ln c x dx a y by k -=-++化简为d c x a by xe y e k --=±考虑条件0,0x y ≥≥以及0y =是解, 可得方程的通解为 d c x a by x e y e k --=这里0k ≥是任意常数.例3 求解Logistic 模型 ()()00d 1,,0d m N N r N N t N N t t N ⎛⎫=-=≥ ⎪⎝⎭. 解分离变量()d d d d m m m N N N N r t N N N N N N ==+--积分可得()ln ln m rt c N N N +=-- 其中c 为任意常数, 化简 ()1rt c m N e N-+=- 即 1m rtN N ce -=+ 这里c c e -=, 代入初值得 001rt m N ce N -=- 最后得到 ()()0011mN r t t mN N N e --=+-.例4 求方程()d d y P x y x=的通解, 其中()P x 是x 的连续函数. 解分离变量 ()d d y P x x y=, 两边积分()ln d y P x x c =+⎰ , 即()d P x x c y e +⎰= , 于是()d P x x c y e e ⎰=±⋅ ,令c c e =± , 于是()d P x x y ce ⎰=. 此外0y =也是方程的解, 但它已包括在上述解中. 故通解为()d P x x y ce ⎰=, 其中c 为任意常数.2.1.2 可化为变量分离方程的类型介绍两种简单情形:1) 形如 d d y y g x x ⎛⎫= ⎪⎝⎭的方程称为齐次方程, 这里()g u 是u 的连续函数.求解方法: 作变量变换 y u x= 于是 d d d d y u x u x x=+, 将上两式代入原方程可得 ()d d u x u g u x+=, 整理后可得 ()d d g u u u x x-=, 这是一个关于u,x 的变量分离方程, 可求解, 再代回原来的变量即可.例5 求解方程d tan d y y y x x x=+. 解这是齐次方程, 令y u x =, d d d d y u x u x x =+, 代入原方程d tan d u x u u u x +=+, 即d tan d u u x x=, 分离变量 d cot d x u u x=, 两边积分 ln sin ln u x c=+ , c 是任意常数, 整理可得sin c u e x =±⋅ ,令c c e =± , 可得sin u cx =, 此外方程还有解tan 0u =, 此解已包括在上式中, 故通解为sin u cx =, c 是任意常数, 代回原来的变量可得到原方程的通解为 sin y cx x=, c 是任意常数. 例6求解方程()d 0d y x y x x+=<. 解将方程改写为d d y y x x= ()0x <, 这是齐次方程, 令y u x =, d d d d y u x u x x =+, 代入原方程得d d u xx=分离变量d x x =, 两边积分()ln x c =-+, 即()()()2ln ln 0u x c x c =-+-+>⎡⎤⎣⎦,这里c 是任意常数, 此外0u =也是方程的解, 它不包括在通解中.代回原来的变量, 得到原方程的通解为()()()2ln ln 0y x x c x c =-+-+>⎡⎤⎣⎦及0y =.或者也可将方程的解表示为 ()()()2ln ,ln 0,0,ln 0.x x c x c y x c ⎧-+-+>⎡⎤⎪⎣⎦=⎨-+≤⎪⎩2) 形如111222d x a x b y c =++ 的方程也可化为变量分离方程. 分三种情形讨论.(1) 120c c ==的情形此时方程可化为 11112222d d yx y x a b a x b y y y g x a x b y x a b ++⎛⎫=== ⎪++⎝⎭ , 从而变为一个齐次方程求解. (2) 11220a b a b =, 即1122a b a b =的情形. 设上面的比值为k , 即1122a b k a b ==, 则方程可写为 ()()22122222d d k a x b y c y f a x b y x a x b y c ++==+++, 令22a x b y u +=, 则有()22d d u a b f u x=+, 这是一个关于u , x 的变量分离方程, 可求解. (3) 11220a b a b ≠, 及1c , 2c 不全为零的情形. 此时方程组11122200a x b y c a x b y c ++=⎧⎨++=⎩ 有解, 设解为(),αβ. 显然()(),0,0αβ≠, 否则与1c , 2c 不全为零矛盾.可通过坐标平移将原点移至(),αβ, 可令,X x Y y αβ=-⎧⎨=-⎩, 从而原方程化为1122d X a X b Y=+, 又转化为情形(1). 对于方程111222d d a x b y c y f x a x b y c ⎛⎫++= ⎪++⎝⎭, 也可用同样方法求解.此外, 下面各种方程也可通过适当的变量变换化为变量分离方程求解.()d d y f ax by c x =++ 令ax by c u ++=. 则()d d d d u y a b a bf u x x=+=+ ()()d d 0yf xy x xg xy y +=令xy u =, 则 ()()()()d d 11d d f u f u u y u y x y x x g u x g u ⎛⎫⎛⎫=+=-=- ⎪ ⎪ ⎪ ⎪⎝⎭⎝⎭. ()2d d y x f xy x= 令xy u =, 则 ()()()d d 11d d u y y x y f u u f u x x x x=+=+=+. 2d d y y xf x x ⎛⎫= ⎪⎝⎭令2y u x =. ()2d d 2d d y u xu x xf u x x=+=，方程变为()2d d f u u u x x-= ()()()(),d d ,d d 0M x y x x y y N x y x y y x ++-=, 其中M, N 为x, y 的齐次函数. 令y u x =. d d d d y u u x x x=+，方程变为12d d 10d d y y y y y y g g x x x x x x ⎛⎫⎛⎫⎛⎫⎛⎫++-= ⎪⎪ ⎪⎪⎝⎭⎝⎭⎝⎭⎝⎭，代入可得()()12d d 10d d u u g u u u x g u u x u x x ⎛⎫⎛⎫⎛⎫++++-= ⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭，整理可得()()()()()21121d d u g u u x x ug u g u +=-+ 例7 求解方程d 1d 3y x y x x y -+=+-. 解解方程组1030x y x y -+=⎧⎨+-=⎩, 解得1,2x y ==, 令12X x Y y =-⎧⎨=-⎩原方程变为 d d Y X Y X X Y-=+ 再令 Y u X=, 即 Y uX = 方程又变为 d 1d 1u u u X X u-+=+, 分离变量得 2d 1d 12X u u X u u+=-- 两边积分 22ln ln 21X u u c=-+-+ 整理得 ()2221c X u u e +-=±令1c c e =± , 则有()22121X u u c +-=代回原变量2212Y XY X c +-=, ()()()()22122121y x y x c -+----= 又2210u u +-=也是原方程的解, 故整理上式可得原方程的通解为22262y xy x y x c +---=这里c 为任意常数.2.1.3 应用举例例7 电容器的充电和放电如图所示的R-C 电路, 开始电容C 没有电荷, 其两端电压为零, 开关合上1后, 电容开始充电, 电压逐渐升高, 充电完毕后, 合上开关2, 电容开始放电, 求充放电过程中电容C 两端的电压c u 随时间t 的变化规律.解对充电过程, 由基尔霍夫第二定律c u RI E +=,由于c Q Cu =, 微分得到 d d d d c u Q I C t t==, 代入可得d d c c u RC u E t+=, 这是c u 满足的微分方程, 分离变量 d d c c u t u E RC=-- 两边积分 11ln c u E t c RC-=-+ 即 1112RC RC t t c c u E e e c e ---=±=,代入初始条件0,0c t u ==可得2c E =-, 于是 ()11RC t c u E e -=-.函数图象如下放电过程类似可讨论.例8 探照灯反射镜面的形状.探照灯要求将点光源射出的光线平行反射出, 求反射镜面的形状.解将点光源设为坐标原点, 设所求曲面为曲线 ()0y f x z =⎧⎪⎨=⎪⎩绕x 轴旋转而成的. 下面求曲线()f x , 如图.过曲线()y f x =上任一点(),M x y 做切线NT , 由反射定律可得12αα=从而 OM ON =切线斜率为 2d tan d y MP x NPα==,又OP x =, MP y =, OM =可得()y f x =满足的微分方程d d y x = 此为齐次方程, 可令y u x =进行求解. 此外, 齐次方程还可令x v y=, 此时x yv =, 微分可得d d d d x v v y y y=+ 代入方程得到d sgn d v v y v y y+=+整理可得d sgn y y y =ln x c ⎛⎫=+ ⎪⎝⎭ 注意到0y >, 可解得(y c v = 代入x v y =可得2y cx =+，整理得 ()22y c c x =+, c 为任意常数.此曲线为抛物线, 反射镜面即为旋转抛物面()222y z c c x +=+.§2.2 线性方程与常数变易法一阶线性微分方程()()()d 0d y a x b x y c x x ++= 当()0a x ≠时可写为 ()()d d y P x y Q x x=+, 下面主要讨论这种形式, 这里()P x 和()Q x 都是连续函数.当()0Q x =时, 方程 ()d d y P x y x= 称为一阶齐线性方程. 若()0Q x ≠称为一阶非齐线性方程.一阶齐线性方程为变量分离方程, 上节例3已求得其通解为()d P x x y ce ⎰=, 其中c 为任意常数.下面讨论一阶非齐线性方程的求解问题. 设想两种方程( 齐次与非齐次) 的解必有某种联系或者形式上的相似. 但是()d P x x y ce ⎰=必不可能是一阶非齐线性方程的解. 设想c 不是常数, 而是一个关于x 的函数()c x , 这是一个待定的函数. 于是, 将()()d P x x y c x e ⎰=代入到一阶非齐线性方程: ()()()()()()()()()d d d d d P x x P x x P x x c x e c x P x e P x c x e Q x x⎰⎰⎰+=+ 即()()()d d d P x x c x Q x e x-⎰= 积分后可得()()()d d P x x c x Q x e x c -⎰=+⎰ 代入得到()()()d d d P x x P x x y e Q x e x c -⎛⎫⎰⎰=+ ⎪⎝⎭⎰ 这就是非齐线性方程的通解. 此法称为常数变易法.例1 求方程()()1d 11d n x y x ny e x x++-=+的通解, 其中n 为常数.解将方程改写为 ()d 1d 1n x y n y e x x x =+++, 先求对应的齐次方程 d d 1y n y x x =+ 的通解, 由公式可知通解为()1n y c x =+.下面应用常数变易法, 令 ()()1n y c x x =+ 微分之可得()()()()1d d 11d d n n c x y x n x c x x x-=+++ 代入原方程()d d x c x e x= 积分可得()x c x e c=+ 因此非齐次方程的通解为()()1n x y x e c =++ 其中c为任意常数. 例2 求方程2d d 2y y x x y =-的通解. 解将方程改写为 d 2d x x y y y=-, 这是一个关于未知函数x 的非齐线性方程.其对应的齐次方程 d 2d x x y y= 的通解是2x cy =.应用常数变易法, 令()2x c y y =并代入方程可得 ()d 1d c y y y=- 积分可得()ln c y y c=-+ 从而原方程的通解为()2ln x y c y =- 这里c是任意常数. 伯努利方程:()()d d n y P x y Q x y x=+ 这里()(),P x Q x 都是连续函数, 且0,1n ≠.可以利用变量变换将伯努利方程化为线性方程. 当0y ≠时, 用n y -乘方程两边()()1d d n n y y y P x Q x x --=+, 令1n z y -=可得 ()d d 1d d n z y n y x x-=- 于是原方程化为 ()()()()d 11d z n P x z n Q x x=-+- 这是一个关于z,x 的线性方程, 可求通解. 此外方程还有解0y =.例3 求方程2d 6d y y xy x x=-的通解. 解这是2n =时的伯努利方程. 令1z y -=可得 2d d d d z y y x x-=- 代入原方程 d 6d z z x x x=-+, 这是线性方程, 它的通解为 268c x z x =+ 代回原变量y 得到2618c x y x =+ 或688x x c y -= 这里c 是任意常数, 此外方程还有解0y =.§2.3 恰当方程与积分因子2.3.1 恰当方程将一阶方程 ()d ,d y f x y x= 写成微分形式(),d d 0f x y x y -=或写成具有对称形式的一阶微分方程()(),d ,d 0M x y x N x y y +=, 这里假设M, N 是x, y 的连续函数, 且具有连续的一阶偏导数.如果方程左端恰好是某个二元函数(),u x y 的全微分, 即()()(),d ,d d ,d d u u M x y x N x y y u x y x y x y ∂∂+≡≡+∂∂ 则称为恰当方程.恰当方程的通解是(),u x y c =, c 是任意常数.下面将解决两个问题(1) 如何判定方程是恰当方程?(2) 如果方程是恰当方程, 如何求出函数(),u x y ? 分析:如果()(),d ,d 0M x y x N x y y +=是恰当方程, 则有 ,u u M N x y ∂∂==∂∂,上两式对y, x 分别再求偏导 22,u M u N y x y x y x∂∂∂∂==∂∂∂∂∂∂ 由假设上面的混合偏导相等, 于是 M N y x∂∂=∂∂ 这是恰当方程的必要条件. 下证这也是恰当方程的充分条件, 即证明当方程满足此条件时能找到函数(),u x y 满足 ,u u M N x y∂∂==∂∂. 首先积分u M x ∂=∂, 得到 ()(),d u M x y x y ϕ=+⎰这里()y ϕ是y 的任意可微函数, 现在选择()y ϕ使u 能满足u N y∂=∂, 即 ()()d ,d d y M x y x N y yϕ∂+=∂⎰, 所以 ()()d ,d d y N M x y x y yϕ∂=-∂⎰. 上式右端与x 无关, 事实上右端对x 的偏导数()()(),d ,d ,d 0N N M x y x M x y x x y x x y N M x y x x y x N M x y ⎡⎤⎡⎤∂∂∂∂∂-=-⎢⎥⎢⎥∂∂∂∂∂⎣⎦⎣⎦∂∂∂⎡⎤=-⎢⎥∂∂∂⎣⎦∂∂=-≡∂∂⎰⎰⎰ 这样()y ϕ就可以积分得到()(),d d y N M x y x y y ϕ⎡⎤∂=-⎢⎥∂⎣⎦⎰⎰ 即求得()(),d ,d d u M x y x N M x y x y y ⎡⎤∂=+-⎢⎥∂⎣⎦⎰⎰⎰, 于是恰当方程的通解即为()(),d ,d d M x y x N M x y x y c y ⎡⎤∂+-=⎢⎥∂⎣⎦⎰⎰⎰ 这里c 是任意常数.例1 求()()222336d 64d 0x xy x x y y y +++=的通解.解这里2236M x xy =+, 2364N x y y =+ 12,12M N xy xy y x∂∂==∂∂ 所以这是一个恰当方程.现求u 使得它同时满足 2236u x xy x∂=+∂ 和 2364u x y y y ∂=+∂ 积分上面第一式可得()3223u x x y y ϕ=++再对y 求导 ()223d 664d y u x y x y y y yϕ∂=+=+∂ 所以()3d 4d y y y ϕ= 积分可得()4y y ϕ=所以得到32243u x x y y =++ 方程的通解为32243x x y y c ++=,这里c 是任意常数.恰当方程可以采用”分项组合”的方法. 此法须熟记一些已知的二元函数的全微分, 如()d d d y x x y xy += 2d d d y x x y x y y ⎛⎫-= ⎪⎝⎭2d d d y x x y y x x -+⎛⎫= ⎪⎝⎭ d d d ln y x x y x xy y ⎛⎫-= ⎪⎝⎭22d d d arctan y x x y x x y y ⎛⎫-= ⎪+⎝⎭ 22d d 1d ln 2y x x y x y x y x y ⎛⎫--= ⎪-+⎝⎭例2 用”分项组合”的办法求例1. 解分组23223d 4d 6d 6d 0x x y y xy x x y y +++= 即342222d d 3d 3d 0x y y x x y +++=再写成()3422d 30x y x y ++= 于是通解为34223x y x y c ++= 这里c 是任意常数.例3 求解方程211cos d d 0x x x y y y y ⎛⎫⎛⎫++-= ⎪ ⎪⎝⎭⎝⎭. 解因为2211,M N y y x y ∂∂=-=-∂∂, 故方程是恰当方程, 分项组合:211cos d d d d 0x x x y x y y y y ⎛⎫++-= ⎪⎝⎭ 即2d d dsin d ln 0y x x y x y y -++= 或 d sin ln 0x x y y ⎛⎫++= ⎪⎝⎭ 所以方程的通解是 sin ln x x y c y ++= 这里c 是任意常数.2.3.2 积分因子这一部分介绍如何将一个非恰当方程转化为恰当方程.如果存在连续可微的函数(),0x y μ≠使得()()()(),,d ,,d 0x y M x y x x y N x y y μμ+= 成为一个恰当方程, 则称(),x y μ为方程()(),d ,d 0M x y x N x y y +=的积分因子.注: 理论上微分方程有解必存在积分因子且不唯一, 从而通解也可能有不同形式.(),x y μ成为方程()(),d ,d 0M x y x N xy y +=的积分因子的充要条件是()()M N y xμμ∂∂=∂∂ 即 M N N M x y y x μμμ⎛⎫∂∂∂∂-=- ⎪∂∂∂∂⎝⎭这是一个关于μ的一阶线性偏微分方程, 事实上, 解这个方程可能会比解原方程更困难. 但对于特殊形式的μ求解会相对容易许多.如果方程存在只与x 有关的积分因子()x μμ=, 则0yμ∂=∂, 这时上述关于μ的一阶线性偏微分方程变为d d M N N x y x μμ⎛⎫∂∂=- ⎪∂∂⎝⎭ 即 d d M N y x x Nμμ⎛⎫∂∂- ⎪∂∂⎝⎭= 于是方程有只与x 有关的积分因子的充要条件是 ()M N y x x Nψ∂∂-∂∂= 这里()x ψ是仅为x 的函数, 如果此条件成立, 则可积分求得方程的一个只与x 有关的积分因子()d x x e ψμ⎰=.同样, 方程有只与y 有关的积分因子的充要条件是 ()M N y x y Mϕ∂∂-∂∂=- 积分可求得方程的一个只与y 有关的积分因子 ()d y y e ϕμ⎰=.例4 试用积分因子法解一阶线性方程.解将一阶线性方程改写为()()d d 0P x y Q x x y +-=⎡⎤⎣⎦设()()(),M x y P x y Q x =+, (),1N x y =-, 计算可得()M N y x P x N ∂∂-∂∂=- 因此方程有只与x 有关的积分因子()d P x x e μ-⎰=, 用它乘以方程两边得()()()()()d d d d d d 0P x x P x x P x x P x e y x e y Q x e x ---⎰⎰⎰-+= 即()()()()d d d d d d 0P x x P x x P x x y e e y Q x e x ---⎰⎰⎰+-=或()()()d d d 0P x x P x x ye Q x e --⎛⎫⎰⎰-= ⎪⎝⎭⎰ 故通解为()()()d d d P x x P x x ye Q x e x c --⎰⎰-=⎰或改写为()()()d d d P x x P x x y e Q x e x c -⎛⎫⎰⎰=+ ⎪⎝⎭⎰. 例5求解方程d d y x x y =-+()0y >.解改写方程d d x x y y x +=即 ()221d 2x y x +=容易看出此方程有积分因子μ=, 用它乘以方程两边22d d x y x +=或写为d x =故通解为x c =+ 或()22y c c x =+.例6 求解方程()d d 0y x y x y +-=解设M y =, N y x =-,1M y ∂=∂, 1N x ∂=-∂, 此方程不是恰当方程.方法1 因为2M N y x M y∂∂-∂∂=--, 故方程有只与y 有关的积分因子 ()2d 2ln 21y y y e e y μ--⎰=== 用它乘以方程的两边得到211d d d 0x x y y y y y+-= 或 2d d d 0y x x y y y y -+= 故通解为 ln x y c y+=. 方法2 将方程改写为d d d y x x y y y -=-由公式知左端有多种积分因子, 其中只和y 有关的积分因子有21y μ=, 用它乘以方程两边可得同样结果. 方法3 改写方程为 d d y y x x y=- 这是一个齐次方程, 令y u x =可求解. 方法4 改写方程为 d 1d x x y y=-这是一个x 作为未知函数的线性方程, 直接用公式可求解. §2.4 一阶隐方程与参数表示一阶隐微分方程的一般形式(),,0F x y y '=如果能够解出y ', 则方程可以采用前面介绍的方法处理. 如果不能解出y '或者解出后形式太复杂, 则可考虑利用变量变换将其变为导数解出的方程. 本节主要介绍一下四种类型1) (),y f x y '= 2) (),x f y y '= 3) (),0F x y '= 4) (),0F y y '=.2.4.1 可以解出y （或x ）的方程1) 先讨论形如 d ,d y y f x x ⎛⎫= ⎪⎝⎭的方程, 假设f 有连续偏导数. 引进参数d d y p x =, 则方程变为 (),y f x p =两边对x 求导可得 d d f f p p x p x∂∂=+∂∂ 这是关于p 的导数解出的方程, 若已求得其通解为 (),p x c ϕ=则原方程的通解为()(),,y f x x c ϕ=.若求得通解为 (),x p c ψ=,则原方程的通解为如下参数形式()()(),,,x p c y f p c p ψψ=⎧⎪⎨=⎪⎩若求得通解为(),,0x p c Φ=,则原方程的通解为()(),,0,x p c y f x p Φ=⎧⎪⎨=⎪⎩ 其中p 是参数, c 是任意常数.例1 求方程3d d 20d d y y x y x x ⎛⎫+-= ⎪⎝⎭的解. 解解出y , 并令d d y p x =, 得到 32y p xp =+两边对x 求导 2d d 322d d p p p p x p x x=++ 整理得23d 2d d 0p p x p p x ++=当0p ≠时, 用p 乘以方程两边3223d d d 0p p x p p x ++=即 423d 04p xp ⎛⎫+= ⎪⎝⎭ 故通解为 4234p xp c += 解出x 并代入32y p xp =+可得 ()43342c p y p p -=+ 因此原方程的参数式的通解为22334212c x p p c y p p ⎧=-⎪⎪⎨⎪=-⎪⎩ 0p ≠ 当0p =时, 直接计算可知0y =也是原方程的解.例2 求方程22d d d d 2y y x y x x x ⎛⎫=-+ ⎪⎝⎭的解. 解令d d y p x =, 得到 222x y p xp =-+ 两边对x 求导 d d 2d d p p p p x p x x x=--+ 即()d 120d p p x x ⎛⎫--= ⎪⎝⎭ 由d 10d p x-= 解得p x c =+, 并得到原方程的通解 222x y cx c =++ 又20p x -=可解得2x p =, 代入后可解得原方程的另一个解 24x y = 此解和通解中每一条曲线相切, 称之为奇解, 下一章将详细介绍.2) 形如 d ,d y x f y x ⎛⎫= ⎪⎝⎭的方程解法与1)类似. 引入参数d d y p x =, 方程变为 (),x f y p =两边对y 求导再以d 1d x y p=代入得 1d d f f p p y p y∂∂=+∂∂ 此为关于p 的导数解出方程, 可求解, 设通解为 (),,0y p c Φ=则原方程的通解为()(),,,0x f y p y p c =⎧⎪⎨Φ=⎪⎩ 例3 求解例1中的方程3d d 20d d y y x y x x ⎛⎫+-= ⎪⎝⎭. 解解出x , 并以d d y p x =代入 ()3,02y p x p p-=≠ 两边对y 求导()()232d d d d 1312p p y y p p y p p p ---= 即3d d 2d 0p y y p p p ++=积分可得42yp p c +=因而 42c p y p-= 代入得 4234c p x p -=, 于是原方程的通解为424342c p x p c p y p ⎧-=⎪⎪⎨-⎪=⎪⎩ 此外还有0y =.2.4.2 不显含y （或x ）的方程3) 形如(),0F x y '=的方程, 可令d d y p y x'==, 则(),0F x p =代表xp 平面上一条曲线, 设这条曲线有参数式 ()(),x t p t ϕψ==因为d d y p x =, 代入上面的参数式可得()()d d y t t t ψϕ'=积分可得()()d y t t t c ψϕ'=+⎰ 于是原方程的参数式通解为()()(),d .x t y t t t c ϕψϕ=⎧⎪⎨'=+⎪⎩⎰ 例4 求解方程3330x y xy ''+-=.解令y p tx '==, 方程变为333230x t x x t +-=可得 331t x t =+ 从而2331t p t =+ 于是()()3233912d d d 1t t y tx x t t -==+积分可得()32331421t y c t +=++ 因此原方程的通解为 ()332331314.21t x t t y c t ⎧=⎪+⎪⎨+=+⎪⎪+⎩4) 形如(),0F y y '=可采用同样方法求解. 令p y '=, 将(),0F y p =表示为参数式()(),y t p t ϕψ==由d d y p x =可得()()d d t t t x ϕψ'=, 所以()()d d t x t t ϕψ'=积分得()()d t x t c t ϕψ'=+⎰于是原方程的通解为()()()d .t x t c t y t ϕψϕ'⎧=+⎪⎨⎪=⎩⎰ 此外若(),00F y =有实根y k =, 则y k =亦是原方程的解.例5 求解方程()()2212y y y ''-=-. 解令2y yt '-=, 代入原方程得 ()2221y yt y t -=即 1y t t=+故21y t '=-所以2d 1d d y x t y t==-', 积分得1x c t=+于是原方程的通解为11x c ty t t ⎧=+⎪⎪⎨⎪=+⎪⎩此外0y '=时原方程变为24y =, 所以2y =±也是原方程的解.§2.5 习题选讲1. 求下列方程的解(1) d sin cos 1d yy x x x+= 解方程变形为sin d cos d d y x x x y x +=, 左端寻找只和x 有关的积分因子2cos x -, 积分可得2cos d tan cos yx x c x c x-=+=+⎰. (2) 2d d d y x x y x y y -=. 解方程两边同乘以21x , 可得 2d d d y x x yy y x -=,即d d y y y x ⎛⎫-= ⎪⎝⎭, 通解为212y y c x +=. (3)d 4sin 1d y ye x x-=- 解方程变为()4sin d d 0y yx e x e y --=, 因为1M Ny x N∂∂-∂∂= 故方程有积分因子x e , 用x e 乘以方程两边可得4sin d d d 0x y x x y e x x e e x e e y --=,即 4sin d x y xe e e x x c =+⎰,所以通解为()2sin cos y x e x x ce -=-+.(5) 22d d 0xxy y xye y x x e y ⎛⎫+-= ⎪⎝⎭解方程可变形为齐次方程22d d xyxyy xye y x x e+=. 令x u y=, 则x uy =, d d d d x u u y y y =+,代入方程可得2d d 1u uu u e u y y ue +=+, 化简并且变量分离1d d u y e u u y ⎛⎫+=- ⎪⎝⎭, 两边积分可得 ln ln uu e y c +=-+代入原来的变量ln x yex c +=.(6) ()1d d 0xy y x x y +-= 解方程改写为2d d d xy x x y y x =-容易看出有积分因子2y -.(7) ()()221d 2d 0x y x x y y +-++-= 解令u x y =+d d d u x y =+, 方程化为()()()()()()()21d 2d 21d 2d d 1d 2d 0u x u yu x u u x u x u u -+-=-+--=++-= 当1u ≠-时, 变量分离可得2d d 1u u x u -=-+, 积分可得3ln 1u u x c -+=-+, ()31u x c e u +-=+代入u 得到方程通解()321x y x y ce +++=,另外1u =-即1x y +=-也是解, 包含于通解中.(8) 23d d y y y x x x=+(伯努利方程)(9)d 32d yy x x=+-(线性方程) (10) 2d d 1d d y y x x x ⎛⎫=+ ⎪⎝⎭(x 可解出的隐方程) (11)2d 1d 3y x y x x y -+=++ 解方程化为()()21d 3d 0x y x x y y -+-++=,可验证M Ny x∂∂=∂∂, 此为恰当方程. (12) d 1d y x y e xe x -⎛⎫+=⎪⎝⎭解方程变形为d d d y y xe x e x xe x --+=容易看出方程有积分因子x e -.(13) ()22d 2d 0x y x xy y +-=(14)d 1d yx y x=++ (15) d d yx y y e x x=+(16) ()d 112d y yx e x-++= 解方程变形为()()12d 1d 0ye x x y --++=方程有只和y 有关的积分因子.(17) ()()2d 1d 0x y x y x y -++=解方程改写为1d 1d 11y x y y x x x -=-++ 此为1n =-的伯努利方程.(18) ()2234d 21d 0x y x x y y +-=提示: 寻找只和y 有关的积分因子.22M Nx y y x∂∂-=∂∂, 12M Ny x M y∂∂-∂∂=--,方程有积分因子12y -.(19) 2d d 240d d y y x y x x x ⎛⎫⎛⎫-+= ⎪ ⎪⎝⎭⎝⎭. (20) 22d 11d y y x ⎡⎤⎛⎫-=⎢⎥ ⎪⎝⎭⎢⎥⎣⎦解令sin y p t '==, 代入得sec y t =±, 由d sin d yt x=可得 sec tan d sin d t t t t x ±=, 整理得2sec d d t t x ±=, 积分得tan x t c =±+消去参数可得()221y x c =++, 此外还有解1y =±.第三章一阶微分方程的解的存在定理在实际应用当中, 如果能够找出方程的通解表达式, 则可以通过它了解和掌握所研究对象的性质. 但是, 很多一阶方程并不能用初等解法求出通解, 而且实际问题中很多情况下都是要求满足初始条件的解, 因此研究初值问题的解的存在和唯一性具有重要的地位.反例解存在而不唯一的例子, 方程d d yx=()0,0的解不止一个,0y =和2y x =都是解.解的存在唯一性的意义: 在解的近似计算中提供的理论依据. 在初值问题中对初值测量偏差所产生的影响.§3.1 解的存在唯一性定理与逐步逼近法3.1.1 存在唯一性定理1) 先考虑导数解出的一阶微分方程()d ,d yf x y x= 这里(),f x y 定义在矩形域00:,R x x a y y b -≤-≤上的连续函数.利普希兹条件如果存在常数0L >使得不等式()()1212,,f x y f x y L y y -≤-对所有()()12,,,x y x y R ∈都成立, 则称函数(),f x y 在R 上满足利普希兹条件, L 成为利普希兹常数.定理1 如果(),f x y 在R 上连续且关于y 满足利普希兹条件, 则方程()d ,d yf x y x=存在唯一的解()y x ϕ=, 定义于区间0x x h -≤上, 连续且满足初始条件()00x y ϕ=这里()(),min ,,max ,x y Rbh a M f x y M∈⎛⎫== ⎪⎝⎭. 证明思路皮卡逐步逼近法首先将为微分方程转化为积分方程()00,d xx y y f x y x =+⎰再任取一个连续函数()0x ϕ代入上面积分方程右端的y , 得到()()()0100,d xx x y f x x x ϕϕ≡+⎰,则()1x ϕ也是连续函数, 如果()()10x x ϕϕ≡, 则()0x ϕ就是积分方程的解, 否则继续把()1x ϕ代入积分方程右端的y ,()()()0201,d xx x y f x x x ϕϕ≡+⎰如果()()21x x ϕϕ≡, 则()2x ϕ就是积分方程的解, 否则可以继续此步骤从而得到一个连续函数列()()()01,,,,n x x x ϕϕϕ可以证明上面的函数列有极限函数()x ϕ, 而它正是积分方程的解. 函数列中的第n 项称为n 次近似解.命题 1 设()y x ϕ=是方程()d ,d yf x y x=的定义于区间00x x x h ≤≤+上, 且满足初始条件()00x y ϕ=的解, 则()y x ϕ=是积分方程()00,d xx y y f x y x =+⎰定义于区间00x x x h ≤≤+上的连续解, 反之亦然.取()00x y ϕ=, 构造皮卡逐项逼近函数列()()()()00001,d x nn x x y x y f ϕϕξϕξξ-⎧=⎪⎨=+⎪⎩⎰ 命题 2 对所有的n , 上式中的()n x ϕ在[]00,x x h +上有定义、连续且满足不等式()0n x y b ϕ-≤.命题3 函数列(){}n x ϕ在[]00,x x h +上是一致收敛的.设()()lim n n x x ϕϕ→∞=, 则()x ϕ连续且()0x y b ϕ-≤. 命题 4 ()x ϕ是积分方程()00,d xx y y f x y x =+⎰定义于[]00,x x h +上的连续解.命题 5 设()x ψ是积分方程()00,d xx y y f x y x =+⎰定义于[]00,x x h +上的一个连续解, 则()()x x ϕψ≡.命题1——5即为定理1的证明.注1 利普希兹条件常用(),f x y 在R 上有对y 的连续偏导代替. 此时, 在R 上,fL y∂≤∂,()()()()212121212,,,f x y y y f x y f x y y y yL y y θ∂+--=-∂≤-注 2 对于线性方程()()d d yP x y Q x x=+, 当()P x 和()Q x 都连续时, 则定理条件就能满足.2) 现在考虑一阶隐方程(),,0F x y y '=根据隐函数定理, 如果在()000,,x y y '的某一邻域内F连续且()000,,0F x y y '=, 而0Fy∂≠'∂, 则y '必可唯一的看成是x, y 的函数(),y f x y '=且导数f FFy y y∂∂∂=-'∂∂∂ 也是连续有界的, 这样(),f x y 即满足利普希兹条件, 于是可得到下面定理.定理2 如果在()000,,x y y '的某一邻域中: 1. (),,F x y y '对所有变元连续, 且存在连续偏导数;2. ()000,,0F x y y '=;3.()000,,0F x y y y '∂≠'∂ 则方程(),,0F x y y '=存在唯一解()0,y y x x x h =-≤满足初始条件()()0000,y x y y x y ''==. 3.1.2 近似计算和误差估计§3.2 解的延拓上节中解的存在唯一性定理是局部性的, 即解只在初值附近较小领域存在. 本节讨论如何延拓解的区间至最大范围.解的延拓定理如果方程()d ,d yf x y x=右端的函数(),f x y 在有界区域G 中连续, 且在G 内关于y 满足局部利普希兹条件, 那么此方程通过G 内任何一点()00,x y 的解()y x ϕ=可以延拓, 直到点()(),x x ϕ任意接近区域G 的边界.推论如果G 是无界区域, 在上面解的延拓定理的条件下, 方程通过点()00,x y 的解()y x ϕ=可以延拓, 以向x 增大的方向来说, 有两种情况: (1) 解()y x ϕ=可以延拓到区间0[,)x +∞;(2) 解()y x ϕ=只可以延拓到区间0[,)x m , m 为有限数,则当x 趋向于m 时, 或者y 无界, 或者()(),x x ϕ趋向于区域的边界.如果函数(),f x y 在整个xy 平面上定义、连续且有界, 同时存在关于y 的一阶连续偏导数, 则方程()d ,d yf x y x=的任一解可以延拓到区间x -∞<<+∞.§3.3 解对初值的连续性和可微性方程()d ,d yf x y x=的解经过初值()00,x y 是唯一的，当初值()00,x y 变化时，解也随之变化. 因此可以把方程的解看成是三元函数()00,,y x x y ϕ=满足()0000,,y x x y ϕ=.解关于初值的对称性设方程()d ,d yf x y x=的经过初值()00,x y 的解是唯一的, 记为()00,,y x x y ϕ=, 则此表达式中(),x y 和()00,x y 可以对调位置, 即成立()00,,y x x y ϕ=。

《常微分方程》PPT课件

两边积分
dxx((tt)) kdt G(y)F(x)C
G(y)
F(x)
P(x) dx
例3 (细菌繁殖模型)在一个理想的环境中，细胞的
繁殖率与细菌的数目成正比，假设t 0时细菌的数
目为 x(t，) 求系统的细菌繁殖规律。
解：设 t示在时x(t刻) 细菌数目，依题意有
t 时y， r k
两边积分
f(x)dx
《常微分方程》PPT课件
本课件PPT仅供大家学习使用学习完请自行删除，谢谢！本课件PPT仅供大家学习使用学习完请自行删除，谢谢！本课件PPT仅供大家学习使用学习完请自行删除，谢谢！本课件PPT仅供大家学习使用学习完请自行删除，谢谢！
6.1 微分方程的根本概念
几何问题引例
物理问题微分方程的根本概念
zxy
2. zaxbyc型方程
作变换 dy abf(z) dx
d y(xy)2 dx 例8. 求方程 n1,令zy1n y 的通解
2
arctan(xy)xC 解：令则 r y(t) krky0 en y0
得方程通解为 VlnV dV lnVadt
将代回ቤተ መጻሕፍቲ ባይዱ原方程通解 y(t)krrky0 en y0
例1 一曲线通过点(1,2),且在该曲线上任意点
。 dz
dx
2处1y 的ddyx 切线斜率为2x，求这曲线的方程
解: 设所求曲线方程为 y = y(x) , 那么有如下关系
式:
r1,2b2ba24ac
①
yx12
②
由 ① 得 y2xdxyf(y) (C为任意常数)
由 ② 得 C = 1, 因此所求曲线方程为 yx21.
dy f (x)g( y) dx

《常微分方程》(王高雄)第三版课后

y= 1 。 1 + ln1 + x
3
dy = 1 + y2 dx xy + x3 y
解：原式可化为：
dy = 1 + y2 •
1
1+ 显然
y2
≠
0, 故分离变量得
y
dy =
1
dx
dx y x + x3
y
1+ y2
x + x3
两边积分得 1 ln1 + 2
y2
=
ln
x
−
1 ln1 + 2
x2
+ ln c (c
c x2 , y
=
0也包含在此通解中。
故原方程的解为原
x2
y2 y2 +
2
=
c
x2,
x
=
0.
解 (2)令xy = u，则原方程化为 du = 1 (u 2 + u 2 + u) = 1 4u
dx x 2 − u 2
x 2−u2
分离变量得 2 − u 2 du = 1 dx，两边积分得 ln y = x 2 y 2 + c，这也就是方程的解。
dx dx
dx t 2
变量分离
t
t2 2 +1
dt
=
dx，两边积分t
−
arctgt
=
x
+
c，代回变量
x + y − arctg(x + y) = x + c
13. dy = 2x − y − 1 dx x − 2 y + 1
解：方程组2x − y −1 = 0, x − 2 y + 1 = 0;的解为x = − 1 , y = 1 33

常微分方程全册ppt课件

z z (5) z ; x y
2u 2u (6) 2 x y uz 0 . 2 x y
都是偏微分方程注: 本课程主要研究常微分方程，同时把常微分方程简称为微分方程或方程
微分方程的阶定义微分方程中出现的未知函数的最高阶导数或微分的阶数称为微分方程的阶数.
z z (5) z ; x y
2 3
(2) xdy ydx 0 ;
d 4x d 2x (4) 5 2 3x sin t ; 4 dt dt
2u 2u (6) 2 x y uz 0 . 2 x y
常微分方程如果在一个微分方程中，自变量的个数只有一个，则这样的微分方程称为常微分方程
两种群竞争模型
Lorenz方程
Lorenz吸引子，蝴蝶效应
对初值的敏感性
分形（fractal）
吸引盆
总结
微分方程反映量与量之间的关系，与时间有关，是一个动态系统从已知的自然规律出发，考虑主要因素，构造出由自变量、未知函数及其导数的关系史，即微分方程，从而建立数学模型数学模型的建立有多种方式研究微分方程的解和解结构的性质，检查是否与实际相吻合，不断改进模型由微分方程发现或预测新的规律和性质
如:
dy (1) 2x dx
是一阶微分方程
(2) xdy ydx 0
d 2x dx (3) tx x 0 2 dt dt
d 4x d 2x (4) 5 2 3x sin t 4 dt dt
3
是二阶微分方程
是四阶微分方程
n阶微分方程的一般形式为
此ppt下载后可自行编辑
教学课件
常微分方程

常微分方程(第三版)课件第一章

2u 2u 2u 8. 2 2 4 xy y
§1.1 Sketch of ODE n阶隐式方程 n阶显式方程方程组
偏微分方程偏微分方程不是微分方程
9. f 2 ( x) sin x
§1.1 Sketch of ODE
微分方程模型举例/Modeling of ODE/
CH.1 Introduction
本章要求/Requirements/
能快速判断微分方程的类型；
掌握高阶微分方程及其初值问题的一般形式；
理解微分方程解的意义。
§1.1 Sketch of ODE
§ 1.1 微分方程概述/ Sketch of ODE/
微分方程理论起始于十七世纪末，是研究自然现象强有力的工具，是数学科学联系实际的主要途径之一。
§ 1.2 基本概念/Basic Conception/
1. 常微分方程和偏微分方程 2. 一阶与高阶微分方程 3. 线性和非线性微分方程 4. 解和隐式解 5. 通解和特解 6. 积分曲线和积分曲线族 7. 微分方程的几何解释-----方向场
§1.2 Basic Conception
常微分方程与偏微分方程/ODE and PDE/
电子课件
常微分方程
Ordinary differential equation
王高雄周之铭朱思铭王寿松编
常微分方程
Ordinary differential equation
• • • • • • • 第一章第二章第三章第四章第五章第六章第七章绪论一阶微分方程的初等解法一阶微分方程的解的存在定理高阶微分方程线性微分方程组定性理论初步1 2 一阶线性偏微分方程
常微分方程的解的表达式中，可能包含一个或者几个常