第1讲 微分方程
常微分方程第一章
2．微分方程的解
微分方程的解是一个函数，函数就有定义域，设为区 间I . 定义2 设函数 y ( x)在区间 I 上连续，且有直到 n 阶导 ( n) y, y ' ,, y后， ( x), ( x),, ( n ) ( x) 数,若用 分别代替方程(1.8)中的 使(1.8)在 内为关于x的恒等式，即 I
称(1.16)是初值问题
y '' g ' y (0) y0 , y (0) v0
(1.17)
的解，初值问题又叫柯西问题. 由以上简单实例可以看出： 1． 微分方程的求解，与一定的积分运算相联系，因此也常把求解微分方 程的过程称为积分一个微分方程，而把微分方程的解称为这个微分方 程的一个积分.由于每进行一次不定积分运算，会产生一个任意常数， 因此仅从微分方程本身求解（不考虑定解条件），则 n 阶微分方程的 n 个任意常数. 解应该包含 2．微分方程所描述的是物体运动变化的瞬时规律，求解微分方程，就是 从这种瞬时规律出发，去获得运动的全过程.为此，需要给定这一运动 的一个初始状态（即初始条件），并以此为基点去推断这一运动的未 来，同时也可以追朔它的过去. 3． 一般对 n 阶微分方程的初值问题的提法是：
在(1.12)两侧对 t 积分一次，得
y' (t ) gt C1
(1.13) (1.14)
其中 C1 是一个任意常数，再把(1.13)对 t 积分一次，就得
1 y (t ) gt 2 C1t C2 2
其中 C2 是另一个任意常数.可知(1.14)是微分方程(1.12)的通解. 通解(1.14)就表示自由落体的运动规律，在(1.14)中含有两个任意 常数.这说明微分方程(1.12)有无穷多个解. 为了要得到特定的物体运动规律，还必须考虑当运动开始时落体是在什 么地方，且以什么样的速度运动的，即下面的初值条件： (1.15) y(0) y0 ，y' (0) v
常微分方程-第一章-初等积分法
黄丹
danh_m@
第一章
初等积分法
微分方程初值问题
y H = f (x; y )的含义 如果将 y 视为系统状态变量，则导数 y H 就是状态的变化率；如果 将自变量视为时间，微分方程 y H = f (x; y ) 可解释为：
=
y (x) 或 x = x(t); y = y (t)。
有：
C 的 速 度 矢 量 为 (xH (t); y H (t)), 则 b=
=
q
(xH (t))2 + (y H (t))2
xH (t) dy dx
s
1+
dy 2
dx
(1)
另：
=
at y x
(2)
黄丹
danh_m@
黄丹
danh_m@
第一章
初等积分法
微分方程是微积分的自然延续 微积分是人类科学史上一个划时代的重大发现 微积分在几何上的应用产生了微分几何 在物理上广泛和深入的应用产生了微分方程
黄丹
danh_m@
第一章
初等积分法
微分方程是微积分的自然延续 微积分是人类科学史上一个划时代的重大发现 微积分在几何上的应用产生了微分几何 在物理上广泛和深入的应用产生了微分方程 微分方程是表达自然规律的一种最为自然的数学语言
ቤተ መጻሕፍቲ ባይዱ
黄丹
danh_m@
第一章
初等积分法
物体下落问题 设质量为 m 的物体，在时间 t = 0 时，在距地 面高度为H 处以初始速度 v (0) = v0 垂直地面下落，求此物体下 落时距离与时间的关系。
常微分方程 PPT课件
分曲线和积分曲线族的概念，只不过此时积分曲线所在的空间维数不同，我们将
在第4章详细讨论.
最后，我们要指出，本书中按习惯用
代替
而 分别代表
本节要点： 1．常微分程的定义，方程的阶，隐式方程，显式方程，线性方程，非线性方程. 2．常微分方程解的定义，通解，特解，通积分，特积分. 3．初值问题及初值问题解的求法. 4．解的几何意义，积分曲线.
所以它们都是一阶齐次方程．因此，一阶齐次微分方程可以 写为
(1.27)
1.3.1 齐次方程的解法 方程(1.27)的特点是它的右端是一个以为
变元的函数，经过如下的变量变换，它能化 为变量可分离方程.
令 则有 代入方程(1.27)得
（1.28）
方程(1.28)是一个 变量可分离方程，当 时，分离 变量并积分，得到它的通积分 （1.29）
常微分方程课件
第一章 初等积方法 第二章 基本定理 第三章 线性微分方程 第四章 线性微分方程组 第五章 定性与稳定性概念 第六章 一阶偏微方程初步
第1讲 微分方程与解 微分方程
什么是微分方程？它是怎样产生的？这是首先要回答的问题.
300多年前，由牛顿(Newton,1642-1727)和莱布尼兹 (Leibniz,1646-1716)所创立的微积分学，是人类科学史上划
（1.13）
显然，方程(1.4)是一阶线性方程；方程(1.5)是一阶非线性方程；方程 (1.6)是二阶线性方程；方程(1.7)是二阶非线性方程.
通解与特解
微分方程的解就是满足方程的函数，可定义如下.
定义1.１ 设函数 在区间I上连续，且有直
到n阶的导数.如果把
代入方程(1.11)，得到在
区间I上关于x的恒等式，
1 常微分方程的基本知识
常微分方程的基本知识 线性微分方程组理论 高阶线性微分方程
2007年8月 南京航空航天大学 理学院 数学系 1
一. 什么是微分方程？
方 程： 含有未知量的等式. 未知量是数. 代数方程 超越方程： 函数方程： 微分方程:
含有自变量,未知函数及其导数的等式。
未知量是函数.
2007年8月
南京航空航天大学 理学院 数学系
dy dny 则称y (x) 为方程 F(x, y, , , n ) 0 dx dx 在I上的一个(显式)解.
微分方程的解: 代入微分方程能使方程成为恒等式的函数.
2007年8月 南京航空航天大学 理学院 数学系 13
例： 验证y sinx, y cosx都是微分方程
y y 0在(,)上的一个解.
2
微分方程:
联系着自变量,未知函数及其导数的关系式.
为了定量地研究一些实际问题的变化规律,往往是 要对所研究的问题进行适当的简化和假设,建立数学 模型,当问题涉及变量的变化率时,该模型就是微分方 程。
2007年8月
南京航空航天大学 理学院 数学系
3
例 1 一 曲 线 通 过 点 (1,2), 且 在 该 曲 线 上 任 一 点
2007年8月
南京航空航天大学 理学院 数学系
20
定解条件
为了从通解中得到合乎要求的特解,必须根据实 际问题给微分方程附加一定的条件,称为定解条件. 求满足定解条件的求解问题称为定解问题. 常见的定解条件是初始条件,相应的定解问题称 为初值问题。
过定点且在定点的切线 y f ( x, y, y) 二阶: ( x x y01) 的斜率为定值的积分曲线. y x x0 y0 , y 0
第1章微分方程和差分方程
第一章 线性微分方程在讲这部分之前，我们先来看一个非常熟悉的物理问题。
一个一维粒子，初始时刻处于点0x x =，初始速度为0v ，受到阻尼作用，求该粒子的运动轨迹。
解：用()x t 表示粒子在任意时刻t 的位置，根据牛顿第二定律F ma =，有mx F =对于阻尼作用F kx =-，于是，粒子的运动方程mx kx =-这是关于时间t 的常微分方程，非常简单。
求解得12()ek t mx t c c -=+结合初始条件0(0)x x =，0(0)x v =，则010mv c x k =+，02mvc k=- 代入得粒子的运动轨迹0()(1e )kt m mv x t x k-=+-这就是这门课程的第二部分——数学物理方程所要讨论的内容：将物理问题表述成数学方程，然后用各种方法来求解方程。
1.1 常系数齐次线性微分方程方程的阶：微分方程中未知函数导数的最高阶数。
线性方程：微分方程中对于未知函数及其所有导数都是一次的，就称为线性方程，高于一次以上就称为非线性方程。
齐次方程：微分方程不含有不包含未知函数的项。
例如 u = 4 u xx ; 二阶线性，x 2u = u xx ; 二阶线性，(u x )2 + u 2 = 1; 一阶非线性。
一、二阶常系数齐次线性微分方程求解 二阶线性微分方程()()()y P x y Q x y f x '''++=若()0f x ≡为齐次，()0f x ≠为非齐次。
方程y ''+py '+qy =0称为二阶常系数齐次线性微分方程, 其中p 、q 均为常数。
能否适当选取r , 使y =e rx 满足二阶常系数齐次线性微分方程, 为此将y =e rx 代入方程y ''+py '+qy =0得(r 2+pr +q )e rx =0由此可见, 只要r 满足代数方程r 2+pr +q =0, 函数y =e rx 就是微分方程的解。
微分方程1
原方程的通解为 y tan( x C ) x .
2、齐次方程
dy y 形如 f ( ) 的微分方程称为齐次方程. dx x
例如： 方程
( xy y 2 )dx ( x 2 2 xy)dy 0
第九章
常微分方程及其应用
1、微分方程的基本概念
2、一阶微分方程
3、二阶线性微分方程
4、微分方程应用举例
背景
函数是反映客观世界运动过程中量与量之间的一种关系，寻求函数 关系在实践中具有重要意义。许多实际问题，往往不能直接找出需要的
函数关系，却比较容易列出表示未知函数及其导数(或微分)与自变量之
间关系的等式．这样的等式就是微分方程．1676年詹姆士．贝努利致牛 顿的信中第一次提出微分方程，直到十八世纪中期，微分方程才成为一 门独立的学科.微分方程建立后，立即成为研究、了解和知晓现实世界 的重要工具．1846年,数学家与天文学家合作，通过求解微分方程，发 现了一颗有名的新星—-—-海王星．1991年，科学家在阿尔卑斯山发 - 现一个肌肉丰满的冰人，据躯体所含碳原子消失的程度，通过求解微分 方程，推断这个冰人大约遇难于5000年以前，类似的实例还有很多．在 微分方程的发展史中，数学家牛顿、莱布尼兹、贝努利家族、拉格朗日、 欧拉、拉普拉斯等等都做出了卓越的贡献．
向全欧洲数学家挑战，提出一个很艰难的问题：“设在垂直平面内有
任意两点，一个质点受地心引力的作用，自较高点下滑至较低点，不 计摩擦，问沿着什么曲线下滑，时间最短？”
这就是著名的“最速降线”问题。它的难处在于和普通的极大极
小值求法不同，它是要求出一个未知函数（曲线），来满足所给的条 件。这问题的新颖和别出心裁引起了很大兴趣，罗比塔、伯努利兄弟、
常微分方程第一章
第一章一阶微分方程1.1学习目标：1. 理解微分方程有关的基本概念, 如微分方程、方程阶数、解、通解、初始条件、初值问题等的定义和提法. 掌握处理微分方程的三种主要方法: 解析方法, 定性方法和数值方法.2. 掌握变量分离法，用变量替换将某些方程转化为变量分离方程, 掌握一阶线性方程的猜测检验法, 常数变易法和积分因子法, 灵活运用这些方法求解相应方程, 理解和掌握一阶线性方程的通解结构和性质.3. 能够大致描述给定一阶微分方程的斜率场, 通过给定的斜率场描述方程解的定性性质; 理解和掌握欧拉方法, 能够利用欧拉方法做简单的近似计算.4. 理解和掌握一阶微分方程初值问题解的存在唯一性定理, 能够利用存在唯一性定理判别方程解的存在性与唯一性并解决与之相关的问题, 了解解对初值的连续相依性和解对初值的连续性定理, 理解适定性的概念.5. 理解自治方程平衡点, 平衡解, 相线的概念, 能够画出给定自治方程的相线, 判断平衡点类型进而定性分析满足不同初始条件解的渐近行为.6. 理解和掌握一阶单参数微分方程族的分歧概念, 掌握发生分歧的条件, 理解和掌握各种分歧类型和相应的分歧图解, 能够画出给定单参数微分方程族的分歧图解, 利用分歧图解分析解的渐近行为随参数变化的状况.7. 掌握在给定的假设条件下, 建立与实际问题相应的常微分方程模型, 并能够灵活运用本章知识进行模型的各种分析.1.2基本知识： (一) 基本概念1. 什么是微分方程:联系着自变量、未知函数及它们的导数（或微分）间的关系式（一般是 指等式），称之为微分方程. 2. 常微分方程和偏微分方程:(1) 如果在微分方程中，自变量的个数只有一个，则称这种微分方程为常微分方程,例如 )(22t f cy dt dy b dt y d =++, 0)(2=++y dtdyt dt dy .(2) 如果在微分方程中，自变量的个数为两个或两个以上，则称这种微分方程为偏微分方程. 例如 0222222=∂∂+∂∂+∂∂zTy T x T , t T x T ∂∂=∂∂422. 本书在不特别指明的情况下, 所说的方程或微分方程均指常微分方程. 3. 微分方程的阶数: 微分方程中出现的未知函数最高阶导数的阶数. 例如，)(22t f cy dt dyb dty d =++ 是二阶常微分方程； 0222222=∂∂+∂∂+∂∂zTy T x T 与t T x T ∂∂=∂∂422是二阶偏微分方程. 4. n 阶常微分方程的一般形式：(,,,...,)0n n dy d yF t y dt dt=,这里(,,,...,)n n dy d y F t y dt dt 是,,,...,n n dy d y t y dt dt 的已知函数，而且一定含有n n d ydt的项；y 是未知函数，t 是自变量. 5. 线性与非线性：（1） 如果方程(,,,...,)0n n dy d y F t y dt dt=的左端是y 及,...,n n dy d ydt dt 的一次有理式，则称(,,,...,)0n n dy d yF t y dt dt=为n 阶线性微分方程. （2） 一般n 阶线性微分方程具有形式：1111()...()()()n n n n n n d y d y dy a t a t a t y f t dt dt dt---++++= 这里1()a t ,…, ()n a t ,()f t 是t 的已知函数.（3）不是线性方程的方程称为非线性方程. （4） 举例：方程)(22t f cy dt dyb dt y d =++是二阶线性微分方程； 方程0sin 22=+φφl gdtd 是二阶非线性微分方程；方程0)(2=++y dtdy t dt dy 是一阶非线性微分方程. 6. 解和隐式解：如果将函数()y t ϕ=代入方程(,,,...,)0n n dy d yF t y dt dt=后，能使它变为恒等式，则称函数()y t ϕ=为方程的解. 如果关系式,0t yΦ=（）决定的隐函数()y t ϕ=是方程的解，则称,0t yΦ=（）为方程的隐式解. 7. 通解与特解：把含有n 个独立的任意常数n c c c ,...,,21的解 12(,,,...,)n y t c c c ϕ=称为n 阶方程(,,,...,)0n n dy d yF t y dt dt =的通解. 其中解对常数的独立性是指，对ϕ及其 1n -阶导数11,...,n n d d dt dtϕϕ--关于n 个常数 n c c c ,...,,21的雅可比行列式不为0, 即 1212(1)(1)(1)120n n n n n nc c c c c c c c c ϕϕϕϕϕϕϕϕϕ---∂∂∂∂∂∂'''∂∂∂∂∂∂≠∂∂∂∂∂∂L L M M L M L.为了确定微分方程一个特定的解，通常给出这个解所必须满足的条件，称为定解条件.常见的定解条件是初始条件, n 阶微分方程(,,,...,)0n n dy d yF t y dt dt =的初始条件是指如下的n 个条件： 1(1)(1)00001,,...,n n n dy d y t t y y y y dt dt---====，，这里(1)(1)0000,,,...,n t y y y -是给定的n+1个常数. 求微分方程满足定解条件的解，就是所谓定解问题. 当定解条件为初始条件时，相应的定解问题称为初值问题. 把满足初始条件的解称为微分方程的特解. 初始条件不同，对应的特解也不同.(二) 解析方法1．变量分离方程 形如()()dyf t y dtϕ=的方程为变量分离方程，其中(),()f t y ϕ分别为,t y 的连续函数.方程解法如下：若()0y ϕ≠，则()()()()dyf t dt y dyf t dt cy ϕϕ==+⎰⎰上式确定方程的隐式通解. 如果存在0y ，使得()00y ϕ=，则0y y =也是方程的解. 2. 可化为变量分离方程的方程(1) 齐次方程形如 ()dy yg dt t=的方程为齐次方程，()g u 为u 的连续函数. 解法如下：做变量替换y u t =，即y ut =，有dy dut u dt dt=+，从而原方程变为()du t u g u dt +=，整理有()du g u udt t-=，此为变量分离方程，可求解. (2) 形如111222a tb yc dy dt a t b y c ++=++的方程, 其中121212,,,,a a b b c c , 为常数. ●111222a b c k a b c ===的情形. 此时方程化为,dyk dt=可解得y kt c =+. ●11220,a b a b =即1122a bk a b ==的情形: 令 22,u a t b y =+ 则有 122222ku c du dya b a b dt dt u c +=+=++ 此为变量分离方程. ●11220a b a b ≠的情形 对120c c ==的情况, 直接做变量替换y u t=. 当12,c c 不全为零, 求 11122200a t b y c a t b y c ++=⎧⎨++=⎩的解为t y αβ=⎧⎨=⎩. 令 T t Y y αβ=-⎧⎨=-⎩, 则方程组化为112200a T bY a T b Y +=⎧⎨+=⎩. 原方程化为12()a T bY dY Yg dT a T bY T+==+的齐次方程可求解. 3．一阶线性微分方程(1) 一般形式：()()()0dya tb t yc t dt++=，若()0a t ≠，则可写成()()dyP t y Q t dt=+的形式. (2) 一阶齐次线性微分方程：()dyP t y dt =，通解为(),P t dt ce c ⎰ 为任意常数.(3) 一阶非齐次线性微分方程：()()dyP t y Q t dt=+，()0Q t ≠.(4) 齐次线性微分方程的性质性质1 必有零解 0y =;性质2 通解等于任意常数c 与一个特解的乘积; 性质3 任意两个解的线性组合也是该微分方程的解. (5) 非齐次线性微分方程的性质性质1 没有零解;性质2 非齐次方程的解加上对应齐次方程的解仍为非齐次方程的解; 性质3 任意两个非齐次方程的解的差是相应齐次方程的解. (6) 一阶非齐次线性微分方程的解法：(i) 猜测-检验法对于常系数的情形，即 ()P t 为常数, 此时方程为()dyay Q t dt=+, a 为常数. 对应齐次方程的通解为atce , 只需再求一个特解, 这时根据()Q t 为特定的函数,可猜测不同的形式特解. 事实上, 当()BtQ t Ae =, ,A B 为给定常数, 且B a ≠时可设待定特解为BtCe , 而当B a =时, 可设特解形式为BtCte , 后代入方程可确定待定常数C . 当()Q t 为cos ,sin At At 或它们的线性组合时, 其中A 为给定常数. 这时可设待定特解为cos sin B At C At +代入方程后确定,B C 的值. 当()Q t 具有多项式形式1011n n n n a t a t a t a --++++L , 其中01,,n a a a L 为给定常数且00a ≠, 这时可设待定特解为1011n n n n b t b t b t b --++++L 代入方程可求得,0,1,,i b i n = L 的值. 对于()Q t 有上述几种线性组合的形式, 则可设待定特解是上述形式特解的线性组合. (ii) 常数变易法: 令()()P t dty c t e ⎰=，代入方程，求出()c t 后可求得通解为()()(())P t dtP t dty e Q t e dt c -⎰⎰=+⎰.(iii) 积分因子法: 方程改写为()()dyP t y Q t dt-=, 将()P t dt e μ-⎰=, 乘方程两端得 ()()()()()P t dt P t dtP t dt dy e e P t y Q t e dt---⎰⎰⎰-= 即 ()()()()P t dtP t dt d ye Q t e dt--⎰⎰=, 从而通解为 ()()()P t dt P t dt ye Q t e dt c --⎰⎰ =+⎰，即 ()()(())P t dt P t dt y e Q t e dt c -⎰⎰= +⎰.注意, 非齐次线性微分方程通解的结构是: 非齐次线性微分方程的通解等于其对应的齐次线性微分方程的通解加上非齐次线性微分方程的一个特解.4. 伯努利(Bernoulli)方程. 形如()()n dyP t y Q t y dt=+的方程, 其中 n 是常数且0,1,(),()n P t Q t ≠ 是连续函数, 称为伯努利方程. 伯努利方程可通过变量替换 1nz y-=化为(1)()(1)()dyn P t z n Q t dt=-+-, 这是关于未知函数z 的线性方程, 可求其通解.(三) 定性方法与数值方法：1. 斜率场：一阶微分方程(,)dyf t y dt =的解()y t ϕ=代表ty 平面上的一条曲线，称之为微分方程的积分曲线. 微分方程(,)dyf t y dt=的通解()y t ϕ=，c 对应于ty 平面上的一族曲线，称之为微分方程的积分曲线族. 满足初始条件00()y t y =的特解就是通过点00(,)t y 的一条积分曲线. 方程(,)dy f t y dt=的积分曲线上的每一点(,)t y 处的切线斜率dydt 刚好等于函数(,)f t y 在这点的值. 也就是，积分曲线的每一点(,)t y 以及这点上的切线斜率dydt恒满足方程；反之，如果在一条曲线每点上其切线斜率刚好等于函数(,)f t y 在这点的值，则这一条曲线就是方程的积分曲线. 这样，可以用(,)f t y 在ty 平面的某个区域D 内定义过各点的小线段，其斜率为(,)f t y ，一般称这样的小线段为斜率标记. 而对ty 平面上D 内任一点(,)t y , 有这样一个小线段与之对应, 这样在D 内形成一个方向场, 称为斜率场. 斜率场是几何直观上描述解的常用方法2. 欧拉方法：求微分方程初值问题00(,)()dyf t y dty t y⎧=⎪⎨⎪=⎩ 的解，可以从初始条件00()y t y =出发，按照一定的步长t ∆ 依照某种方法逐步计算微分方程的近似解()n n y y t =, 这里0n t t n t =+∆这样求出的解称为数值解. 利用欧拉公式10(,),n n n n n y y f t y t t t n t +=+∆ =+∆,可求初值问题的近似解，这种方法称为欧拉方法.欧拉方法具有一阶误差精度 .如果我们先用欧拉公式求出近似解,再利用梯形公式进行校正, 得到的近似解将具有2阶误差精度, 具体为 预测: 1(,)n n n n y y f t y t +=+∆,校正: 11,11[(,)()]2n n n n n n y y f t y f t y t ++ +=++∆, 这种方法称为改进的欧拉方法.(四) 解的存在性、唯一性及解对初值的连续相依性1. 利普希茨（lipschitz ）条件: 函数(,)f t y 称为在区域2D ⊆R 内关于y 满足利普希茨条件，是指如果存在常数0L >，使得不等式1212(,)(,)f t y f t y L y y -≤-对于所有的12(,),(,)t y t y D ∈都成立, 其中L 称为利普希茨常数. 2. 基本定理(1) 解的存在性定理: 设(,)f t y 在矩形区域2{(,):,}D t y a t b c y d =∈ << <<R 内连续. 如果00(,)t y D ∈, 那么，存在0ε> 和函数()y t , 定义于区间00(,)t t εε-+内，是初值问题00(,)()dyf t y dt y t y⎧=⎪⎨⎪=⎩ 的解. (2) 解的唯一性定理: 设(,)f t y 在矩形区域2{(,):,}D t y a t b c y d =∈ << <<R 内连续且关于y 满足利普希茨条件. 如果00(,)t y D ∈并且12(),()y t y t 是初值问题00(,)()dyf t y dt y t y⎧=⎪⎨⎪=⎩在区间00(,)t t εε-+内的两个解，那么对任意的00(,)t t t εε∈-+,12()()y t y t =,即解是唯一的.注记1: 存在性定理和唯一性定理结合在一起称为初值问题解的存在唯一性定理，叙述如下：设(,)f t y 在矩形区域2{(,):,}D t y a t b c y d =∈ << <<R 内连续且关于y 满足利普希茨条件. 如果00(,)t y D ∈, 那么，存在0ε> 和函数()y t , 定义于区间00(,)t t εε-+内，是初值问题00(,)()dyf t y dt y t y⎧=⎪⎨⎪=⎩ 的唯一解. 因而当我们判断初值问题解的存在唯一性时，要检查(,)f t y 需要满足的条件.注记2: 由于利普希茨条件较难检验，常用(,)f t y 在2{(,):,}D t y a t b c y d =∈ ≤≤ ≤≤R上对y 有连续偏导数来代替. 事实上，如果在D 上y f ∂∂存在且连续，则yf∂∂在D 上有界. 设在D 上L yf≤∂∂, 这时 2121212(,())(,)(,)f t y y y f t y f t y y y yθ∂+--=-∂21y y L -≤,其中 12(,),(,),01t y t y D θ∈ <<. 但反过来满足利普希茨条件的函数(,)f t y 不一定有偏导数存在. 例如(,)||f t y y = 在任何区域内都满足利普希茨条件，但它在0y =处没有导数.(3) 解对初值的连续相依性定理设(,)f t y 在矩形区域2{(,):,}D t y a t b c y d =∈ << <<R 内连续且关于y 满足利普希茨条件. 如果00(,)t y D ∈，00(,,)y t t y ϕ=是初值问题00(,)()dyf t y dt y t y⎧=⎪⎨⎪=⎩在区间00(,)t h t h -+内的解，其中 0h >,那么，对任意给定的0>ε，必能找到正数(,)0h δδε=>，使得 当2220000t t y y δ-+-<（）（）时，初值问题00(,)()dyf t y dt y t y⎧=⎪⎨⎪=⎩的解00(,,)y t t y ϕ=在区间00(,)t h t h -+内也有定义，并且0000|(,,),,|,t t y x t y ϕϕε-<（） 00(,)t t h t h ∈-+. (4) 解对初值的连续性定理设(,)f t y 在矩形区域2{(,):,}D t y a t b c y d =∈ << <<R 内连续且关于y 满足利普希茨条件. 如果00(,)t y D ∈，00(,,)y t t y ϕ=是初值问题00(,)()dyf t y dt y t y⎧=⎪⎨⎪=⎩的解, 那么00(,,)t t y ϕ作为00,,t t y 的三元函数在它存在的范围内是连续的.3. 初值问题的适定性当一个微分方程初值问题的解存在, 唯一并且解连续的依赖于初始条件时, 我们称该问题是适定的. 那么, 对于常微分方程初值问题00(,)()dyf t y dt y t y⎧=⎪⎨⎪=⎩, 只要在00(,)t y 所在的区域内,(,)f t y 连续并且关于y 满足利普希茨条件, 则该初值问题是适定的.(五) 自治方程的平衡点与相线1. 自治方程 当一阶微分方程(,)dy f t y dt=的右端项只是y 的函数而与自变量t 无关, 即()dyf y dt =时, 称为自治方程.2. 平衡解与平衡点 对自治方程()dyf y dt=而言, 若()0f y =有解0y y =, 则称 0()y t y ≡ 是方程的平衡解, 而点0y 称为方程的一个平衡点. 3. 相线相线是仅仅对自治方程()dyf y dt=而言的一种简化的斜率场. 自治方程的斜率场在水平直线上的斜率标记是一样的, 这样只要知道一条竖直直线上的斜率标记, 我们就可以知道整个斜率场. 因而, 在一个竖直的直线上, 我们用向上的箭头表示正的导数, 用向下的箭头表示负的导数. 对于导数为零的点, 用实心圆点来标记它, 则形成该自治方程的相线. 4. 画相线的基本步骤 (1) 画出y -线(竖直线),(2) 找到并在y -线上标记平衡点,不连续点或定义域外的点 (3) 找到()0f y >的区间, 在这些区间上画上向上的箭头, (4) 找到()0f y < 的区间, 在这些区间上画上向下的箭头.5. 初值问题0(),(0)dyf y y y dt= =解的渐近行为 (1) 趋向于平衡点, 如01()(1),2f y y y y =- =;(2) 在无限时间内趋于无穷, 如0(),1f y y y = =;(3) 在有限时间内趋于无穷(爆破), 如20(),1f y y y = =;(4) 在有限时间内停止(导数趋于无穷), 如 01(),1f y y y=- =. 6. 平衡点的分类对于自治方程()dyf y dt=, 如果()f y 在(,)-∞+∞ 内连续, 那么它的解当t 增加时要么(在有限或无限时间里)趋于+∞或-∞, 要么渐近趋于平衡点. 因而,平衡点在自治方程的研究中起着重要的作用. (1) 汇对于初值接近0y 的解, 当t 增加时, 都渐近趋于0y . 对于这样的平衡点0y , 我们称之为汇, 它是稳定的. (2) 源对于初值接近0y 的解, 当t 增加时, 都远离0y . 对于这样的平衡点0y , 我们称之为源,它是不稳定的. (3) 结点既不是源也不是汇的平衡点, 我们称之为结点,它也是不稳定的. 7. 判断平衡点类型的线性化方法 1. 如果0y 是自治方程()dyf y dt=的一个平衡点, 即0()0f y =, 那么 (1) 0y 是源当且仅当()f y 在0y 附近严格单调增加; (2) 0y 是汇当且仅当()f y 在0y 附近严格单调递减. 2. (线性化定理) 如果0y 是自治方程()dyf y dt=的一个平衡点, 即0()0f y =, 并且()f y 是连续可微的, 那么 (1) 若0()0f y '> 则0y 是源; (2) 若0()0f y '<, 则0y 是汇;(3) 若0()0f y '=, 则需要进一步的信息决定其类型.(六) 分歧一阶微分方程解的渐近行为随参数变化发生了类型的变化, 我们称之为分歧现象(或分支, 分叉).1. 分歧发生的条件 对于单参数微分方程族()(,)dy f y f y dtμμ==, 0μμ=是一个分歧值的必要条件是: 存在平衡点0y , 使得 0000(,)(,)0f f y y y μμ∂==∂. 这样我们要找分歧点可以通过求解方程组 (,)0(,)0f y f y y μμ=⎧⎪∂⎨=⎪∂⎩, 得到解 00(,)y μ,0μ为可能的分歧值, 而0y 是可能发生分歧的平衡点. 2. 分歧图解与分歧类型分歧图解是y μ 平面上方程在分歧值附近的所有相线的图, 用以强调当参数经过分歧值时相线所经历的变化.(1) 鞍结点分歧在分歧图解(图1-1)中, 当μ从左到右经过分歧值0μ时, 方程的平衡点从两个变为一个再变为不存在, 这种分歧一般称之为鞍结点分歧. 这类分歧图解在分歧值附近是抛物线的形状(2) 在分歧图解(图1-2)中,当μ从右到左经过分歧值0μ=时, 方程的平衡点由三个变为一个, 这种分歧一般称之为音叉分歧.图 1-1 鞍结点分歧 图 1-2 音叉分歧图 1-3 跨越分歧 图 1-4 复合分歧(3) 在分歧图解(图1-3)中, 当0μ= 时, 方程有一个平衡点; 当0μ≠ 时, 方程有两个平衡点. 0μ=是一个分歧值. 虽然在分歧值的两侧方程都有两个平衡点，但平衡点的稳定性会改变. 当0μ> 时, 0y =是一个汇,它是稳定的; 当0μ<时, 0y =是一个源,它 是不稳定的. 这类分歧一般称为跨越分歧.(4) 在分歧图解(图1-4)中, 当 μ从左到右变化时,相应的方程平衡点依次由一个变为两个,三个,两个再变回一个, 这种分歧一般称之为复合分歧.(七) 一阶微分方程的应用1. 增长和衰减问题设 ()S t 为正在增长或衰减的某研究对象的总量. 如果假设它随时间的变化率dS dt与当前数目成正比, 其比例系数为 k , 则有 dS kS dt =, 或 0dS kS dt-=. 设()S t 可微, 因而是连续函数. Malthus 人口模型满足上述微分方程, 虽然对人口问题, ()S t 是离散的, 只能取整数值, 但该模型系统在一定情况下提供了很好的近似对某一生物种群进行研究时, 该生物种群的增长往往受资源和环境的限制, 引进参量N , 称为最大承载量, 用以表示自然资源和环境条件所能容纳的最大数量, 并且假定 （1）当基数很小时，增长率与当前数成正比；（2）当基数很大，达到资源和环境不能承受的时候，数量开始减少，即增长率为负的. 此时方程可改写为(1)dS S k S dt N=-, 称为具有增长率k 和最大承载量N 的Logistic 模型，该模型最早由荷兰生物学家 Verhulst在1838年提出.2. 温度问题牛顿冷却定律(亦适应于加热的情况)说明物体的温度随时间的变化率与物体所处的周围环境的温差成正比, 设 T 是物体的温度, T 是所处环境的温度, 那么物体温度随时间的变化率为dT dt, 牛顿冷却定律可表示为 ()dT k T T dt=--, 其中k 是正的比例系数, 而负号表示在冷却过程中, 物体温度 T 大于周围环境温度T , 变化率0dT dt <. 在加热过程中0dT dt>, 此时T T <. 3. 稀释问题一容器最初容纳0V 升盐水溶液, 其中含盐 a 克. 每升含盐 b 克的盐水溶液以e 升/分的速度注入,同时, 搅拌均匀的溶液以f 升/分的速度流出, 问在任何时刻 t , 容器中的含盐量.设Q 为任何时刻容器中的含盐量. Q 的变化率dQ dt等于盐的注入率减去流出率. 盐的注入率是 be 克/分. 要决定流出率, 首先计算在时刻t , 容器中的溶液的体积, 它等于最初的体积0V 加上注入的体积 et 后减去流出的体积ft . 因此, 在任一时刻t , 盐水的体积是 0V et ft +-. 在任何时刻的浓度是 0Q V et ft +-, 由此得流出率为 0Qf V et ft +-/分. 于是得到微分方程 0dQ Qf be dt V et ft =-+-, 即 0dQ f Q be dt V et ft+=+-, 这是一个一阶线性方程.4. 电路一个简单的 RC 回路是包含有电阻R (欧姆), 电容C (法拉)和电源V (伏特),如图1-5.图1-5 RC 电路 图1-6 RL 电路由电路学知识，C 的电压()v t 与电阻R 的电压之和应为电源的电压()V t . 电路中的电流I (安培)为 ()dQ dCv t dv I C dt dt dt ===, 其中 Q 为电量从而R 处的电压为 dv RI RC dt=, 由此我们可以建立RC 电路的模型如下：()dv RC v V t dt +=, 即 ()dv V t v dt RC-=. 对于一个包含有电阻R (欧姆), 电感L (亨利)和电源V (伏特)的RL 回路,如图1-6. 电路中的电流应满足的基本方程为 dI R V I dt L L +=.(八) 种群生态学中的模型设()y t 表示一个生物种群的数量, t 为时间, 最简单的种群模型是 Malthus 模型dy ky dt=. Malthus 模型的解()(0)kty t y e =预测了种群数量的指数增长.由于种群数量大的时候，对资源的竞争加剧，因此单位增长率会随种群数目增大而减小，因此更为合理的假设是()dy yf y dt= (*) 这里()f y 是单位增长率，因为dy dt 为增长率，y 是种群数量, 而()/dy f y y dt =. 当考虑种群数量的变化时.对()f y 而言, 其代数形式并不重要, 而关键是其单调性, 凸凹性, 这样我们可以对其进行大致分类:(1) 若()f y 在[0,)+∞上是递减的，称(*)为 Logistic 型；(2) 若()f y 在[0,)+∞上是先增后减的，称(*)为 Allee 效应型；(3) 若()f y 在[0,)+∞上是递减再递增最后递减的，称(*)为 Hysteresis 型.1.3典型例题：例1 考虑微分方程 3220dy y y y dt=--, 问 (1) y 为何值时, ()y t 将保持不变?(2) y 为何值时, ()y t 将增加?(3) y 为何值时, ()y t 将减少?解: 因为当0dy dt =时, ()y t 将保持不变; 当0dy dt >时, ()y t 将增加; 当0dy dt<时, ()y t 将减少. 由3220dy y y y dt=--知, (1) 当32200y y y --=, 即0,4,5y y y = =-=时, ()y t 将保持不变.(2) 当32200y y y -->, 即40y -<< 或5y > 时, ()y t 将增加.(3) 当32200y y y --<, 即4y <- 或05y << 时, ()y t 将减少.例2 假定在鄱阳湖中一种鱼类的数量()S t 随时间的变化按Logistic 模型增长, 增长率为k , 最大承载量为N , 即有 (1)dS S k S dt N=-. 如果每年要从湖中捕获一定量的鱼, 试按下述不同情形对模型做适当修改,(1) 每年捕获10吨?(2) 每年捕获总量的三分之一?(3) 捕获量与总量的平方根成正比?解: (1)(1)10dS S k S dt N=--. (2) 1(1)3dS S k S S dt N =--. (3) (1)dS S k S l S dt N =--, 其中 l 是捕获量与总量平方根的比例系数. 例3 求解方程dy t dt y=- 解：变量分离得 ydt tdy =-.两边积分 22222y t c =-+. 通解为 22t y c +=, c 为任意正常数. 例4 求解方程231dy y dx xy x y+=+ 解：变量分离得 221(1)ydy dx y x x =++, 两边积分 2221()1(1)1ydy dx x dx y x x x x ==-+++⎰⎰⎰.即 22111ln(1)ln ||ln(1)22y x x c +=-++, 1c 为任意常数, 整理得222(1)(1)y x cx ++=, 12c c e =为任意正的常数.例5 求解方程tan dy y xy dx x-=. 解: 将方程改写为 tan dy y y dx x x=+, 这是齐次方程, 做变量替换y u x =，即y ux =，有dy du x u dx dx=+，从而原方程变为 tan du x u u u dx +=+ 即tan du u dx x= 利用分离变量法求得 sin u cx =, 代回原变量得通解为sin y cx x=, c 为任意常数 例6 求解方程22dy x y x y dx=+-. 解: 方程改写为2sgn 1()dy y y x dx x x =+⋅- 令u=y x ，则y ux =，从而2sgn 1du x u u x u dx+=+⋅- 当210u -≠时，2sgn 1dux dx x u =-, arcsin sgn ln u x x c =⋅+, 即 arcsin sgn ln y x x c x=⋅+, c 为任意常数．此外，还有解210u -=，即22y x =．例7 求解方程 13dy x y dx x y -+=+- 解: 解方程组 1030x y x y -+=⎧⎨+-=⎩的解 为 12x y =⎧⎨=⎩. 令 12X x Y y =-⎧⎨=-⎩ , 则原方程化为 dY X Y dX X Y -=+.令 Y u X = ,则可化为变量分离方程 21,12dX u du X u u +=-- 解得 222Y XY X c --=, 代回原变量 有22262y xy x y x c +---=, c 为任意常数.例8 求解方程2()dy y b t dt -=, 其中 (1) 2()1b t t t =++,(2) 4()t b t e =(3) 2()3t b t e =(4) ()cos3b t t =(5) 422()3cos31t t b t e e t t t =+++++解: 对应齐次方程的通解为 2t y ce =, 下面用猜测-检验法求特解(1) 设 21y At Bt C =++ 代入 221dy y t t dt-=++, 有 2222()1At B At Bt C t t +-++=++解得 1,1,12A B C =-=- =-, 从而21112y t t =---, 原方程的通解为 22112t y ce t t =---, c 为任意常数. (2) 设 42t y Ae = 代入 42t dy y e dt-=, 有 44442t t t Ae Ae e -=解得 12A =, 从而4212t y e =, 原方程的通解为 2412t t y ce e =+, c 为任意常数. (3) 不能设2t Ae 形式的特解, 因为它是相应齐次方程的解,不可能是非齐次方程的解,设 23t y Ate = 代入 22t dy y e dt-=, 有 2222223t t t t Ate Ae Ate e +-=解得 3A =, 从而233t y te =, 原方程的通解为2223(3)t t t y ce te c t e =+=+, c 为任意常数.(4) 设 4cos3sin 3y A t B t =+ 代入 2cos3dy y t dt-=, 有 3sin33cos32(cos3sin3)cos3A t B t A t B t t -+-+=有 2310320A B A B -+-=⎧⎨ --=⎩, 解得 23,1313A B =- =, 从而423cos3sin 31313y t t =-+, 原方程的通解为 223cos3sin 31313t y ce t t =-+, c 为任意常数.(5) 根据叠加原理, 由前面4个小题知方程有特解422512313cos3sin 31213132t t y e te t t t t =+-+--- 原方程的通解为242212313cos3sin 31213132t t t y ce e te t t t t =++-+---,c 为任意常数. 例9 求方程22dy y dx x y =-的通解. 解: 将方程改写为222dx x y x y dy y y-==-. 求齐次线性微分方程 2dx x dy y=, 得通解为2x cy =. (常数变易法) 令 2()x c y y =代入原方程 得()1,()ln ||dc y c y y c dy y=- =-+, 从而可得原方程的通解为2(ln ||)x y y c =-+, c 为任意常数.例10 求方程26dy y ty dt t=-的通解. 解: 此为 2n =的伯努利方程. 令 1z y -=可得 6dz z t dt t=-+，此为线性方程可求通解为 268c t z t =-+, 代回原变量得 2618c t y t =-+, 即 688t t c y -=, c 为任意常数. 此外, 原方程还有解0y =.例11 用积分因子法求解方程 32(1)1dy y t dt t =+++. 解: 方程改写为 32(1)1dy y t dt t -=++, 积分因子为 221()(1)dt t t e t μ- -+⎰==+, 乘方程两端得 23(1)2(1)1dy t t y t dt--+-+=+, 即 2(1)1d t y t dt -+=+, 有 421(1)(1)2y t c t =+++, c 为任意常数.例12 若()f t 连续且0()()10t f t f s ds t = , ≠⎰, 试求函数()f t 的一般表达式. 解: 设0()()t F t f s ds =⎰, 则()F t 可导且()()F t f t '=, 这样有1,dF F FdF dt dt= =, 得 2()2,()2F t t c F t t c =+ =±+, 又(0)0F =, 得0c =. 从而 ()2F t t =±,进而 1()()2f t F t t'==±. 例13 求具有性质 ()()()1()()y t y s y t s y t y s ++=- 的函数 ()y t , 已知(0)y '存在. 解: 首先令 0s =, 由已知可得 ()(0)()1()(0)y t y y t y t y +=-, 化简有 2(0)(1())0y y t +=, 知 (0)0y =. 由函数的导数定义00202002()()()lim()()()1()()lim ()(1())lim (1()())()1()lim lim 1()()(0)(1())s s s s s y t s y t y t sy t y s y t y t y s sy s y t s y t y s y s y t s y t y s y y t →→→→→+-'=+-- =+ =-+ = -' = + 变形为 2(0)1()dy y dt y t '=+, 积分得 arctan ()(0)y t y t c ' = +, 由(0)0y =, 知 0c =, 所以满足条件的函数为 ()tan (0))y t y t '= (.例14 下面给定8个微分方程和4个斜率场, 请选出斜率场相应的微分方程, 并说明理由. (1) 2dy t dt =- (2) 24dy y dt=- (3) 2dy y t dt =- (4) 2dy t dt =- (5) 24dy y dt =- (6) 2dy y dt =- (7) dy yt t dt =+ (8) 2dy y t dt=+图1-7 图1-8图1-9 图1-10解: 图1-7对应于(4),图1-8对应于(3),图1-9对应于(2),图1-10对应于(7). 这是因为图1-7的斜率场竖直方向上的斜率标记一样, 知方程的右端项仅是自变量t 的函数()f t , 且当 2t >, ()0f t <, 当2t <时, ()0f t >, 只有(4)满足要求. 图1-8的斜率场知方程右端项为(,)f t y 是 ,t y 的函数, 且当 0y <时,(,)0f t y <, 只有(3)满足.图1-9的斜率场知方程为自治方程有平衡点 2,2y y ==-, 且在 22y -<<时,()0f y <, 知只有(2)满足要求.图1-10的斜率场知方程右端项为(,)f t y 是 ,t y 的函数, 且有平衡解 1y =-, 只有(7)满足要求.例15 利用欧拉方法和改进的欧拉方法, 对步长 0.1t ∆=, 在区间[0,1]上求初值问题21,(0)0dyy y dt=+ =的近似解. 解: 这里 200(,)1,0,0f t y y t y =+==. 利用欧拉公式10(,),n n n n n y y f t y t t t n t +=+∆ =+∆,和 改进的欧拉方法,预测: 1(,)n n n n y y f t y t +=+∆,校正: 11,11[(,)()]2n n n n n n y y f t y f t y t ++ +=++∆,分别计算如下表:欧拉方法改进的欧拉方法n n tn y(,)n n f t y 预测的n y校正的n y 真 解tan y t =0 010 0 1 0.1 0.1000 1.0100 0.1000 0.1005 0.1003 2 0.2 0.2010 1.0404 0.2015 0.2030 0.2027 3 0.3 0.3050 1.0930 0.3072 0.3098 0.3093 4 0.4 0.4143 1.1716 0.4194 0.4234 0.4228 5 0.5 0.5315 1.2825 0.5413 0.5470 0.5463 6 0.6 0.6598 1.4353 0.6769 0.6849 0.6841 7 0.7 0.8033 1.6453 0.8318 0.8429 0.8423 8 0.8 0.9678 1.9366 1.0140 1.0299 1.0296 9 0.9 1.1615 2.34911.2360 1.2592 1.2602 10 11.39642.94991.51791.55371.5574例16 讨论微分方程 233dyy dt=在怎样的区域内满足存在唯一性定理的条件,并求通过点(0, 0) 的一切解.解: 由 23(,)3f t y y =, 知它在全平面内连续, 又由于13(,)2f t y y y-∂=∂, 在除去0y =的区域内连续, 从而在除去0y =的有界闭区域内有界, 进而满足利普希茨条件, 知方程满足初始条件00()0y t y =≠的解在充分小的邻域内存在并且唯一. 当 0y =时, 函数0y =是方程过 (0,0) 的解.当0y ≠时, 方程可变形为 2313y dy dt - =, 积分得 3()y t c =+, c 为任意常数.当0c =时, 得特解 3y t = 是过 (0,0) 的另一个解, 其实, 除零解外, 过(0,0)的所有解可以表示为3111(),0,t c t c y t c ⎧- <=⎨ ≥⎩,3222(),0,t c t c y t c ⎧- >=⎨ ≤⎩, 31132212(),(),0,t c t c y t c t c c t c ⎧- <⎪=- >⎨⎪≤≤⎩,其中12,c c 是满足10c ≤,20c ≥的任意常数, 这些解的定义区间为(,)-∞ +∞, 但本质上在充分小的邻域 (,)εε-内方程所确定的过(0,0)的解只有四个,即 函数30,y y t = =, 3,00,t t y t εε⎧ -<<=⎨ 0≤<⎩及30,0,t y t t εε -<<⎧=⎨ 0≤<⎩.例17 举例说明一阶微分方程初值问题00(,)()dyf t y dt y t y⎧=⎪⎨⎪=⎩解的存在唯一性定理中, 关于(,)f t y 在矩形区域2{(,):,}D t y a t b c y d =∈ << <<R 内连续,关于y 满足利普希茨条件是保证解的存在唯一的非必要条件.解: (1) 当连续条件不满足时, 解也可能是存在唯一的. 如方程1,(,)0,y t dyf t y y t dt =⎧==⎨≠⎩, 显然, (,)f t y 在以原点为心的任何矩形区域内不连续, 间断点为直线y t =, 但过原点的解存在唯一, 这个解就是y t =.(2) 当利普希茨条件不满足时, 解也可能是唯一的. 如ln ||,0(,)0,y y y dyf t y y t dt ≠⎧==⎨=⎩, 由于 11111|(,)(,0)||ln ||0||ln ||||0|f t y f t y y y y -=-=⋅-,当 110,ln ||y y → →-∞无界, 因而(,)f t y 在以原点为心的任何矩形领域内不满足利普希茨条件. 然而方程的所有解为 xce y e =±,c 为任意常数, 及 0y =.过原点(0,0)有唯一解 ()0y t =. 例18 对微分方程(2)(5)dyy y y dt=--而言, 利用存在唯一性定理, 说明满足下列初始条件的解是否存在, 如果存在你能否知道这个解或有关这个解的一些性质.(1) (0)6y =， (2) (0)5y =， (3) (0)1y =, (4) (0)1y =-.解: 由方程的右端项为 ()(2)(5)f y y y y =--仅为 y 的函数在全平面上连续可微,从而由存在唯一性定理, 给定初始条件的解是存在并且是唯一的. 首先由()(2)(5)f y y y y =--知方程有()0,()2,()5y t y t y t = = =三个平衡解.(1) 初始条件为 (0)6y =, 初值位于()5y t =的上方, 由唯一性, 满足这个初始。
第一章常微分方程初值问题数值解法
（1.2.3）
其中rn,k(t)为插值余项。 代到（1.2.2）式中得
u ( tn +1 ) = u ( tn ) +
舍去余项 并用uj代替u(tj)即得
∫
tn+1 tn
Ln , k ( t ) dt + ∫ t rn , k ( t ) dt
tn+1
n
(1.2.4 (1.2.5)
Rn , k = ∫
⎡ ∑ ⎣α u
j =0
j =0
αk ≠ 0
（1.2.1）
j n+ j
⎤ − hβ j f n + j ⎦ = 0（数值解满足的差分方程）
因此称（1.2.1）为多步法 或 k-步法。 又因为(1.2.1)关于 u n + j , f n + j 是线性的，所以称为线性多步法。 为使多步法的计算能够进行，除给定的初值u0 外，还要 知道附加初值u1,u2,…,uk-1 ，这可用其它方法计算。 若 β k = 0 则称(1.2.1)是显式的； 若 β k ≠ 0 则方法（1.2.1）是隐式的。 例如，一般线性二步法可写成：
f ( t , u ( t ) ) = Ln , k +1 ( t ) + rn , k +1 ( t )
其中rn,k+1(t)为插值余项。 同理即
un +1 = un + h ∑ bk +1i f ( tn −i +1 , un −i +1 )
i =0
k +1
其中
bk +1i
=∫ ∏
−1
j =0 j ≠i
0
k +1
第一节 微分方程的基本概念
含有未知函数导数 (或微分) 的方程。 1、微分方程：
常微分方程
偏微分方程
例如： (1) y= kx, k 为常数； (2) ( y - 2xy) dx + x2 dy = 0； 1 2 1 y ; (3) mv(t) = mg - kv(t)； (4) y a d 2 g (5) 2 sin 0 ( g , l 为常数). dt l 2、微分方程的阶—— 微分方程中出现的未知函数最高 阶导数的阶数。
x 2
则
1 1 x P ( x ) , Q( x ) e , 2 2
则
1 x P ( x )dx dx , 2 2
P ( x )dx x 2
e
x 2
P ( x )dx
x 2
e ,
x 2
Q( x)e
1 x dx e e dx e , 2
自由项
当 f (x) 0 时，称为二阶线性非齐次微分方程， 当 f (x)=0 时，称为二阶线性齐次微分方程， 方程中 p(x)、 q(x) 和 f (x) 都是自变量的已知连续函数. 例如 y + xy + y = x2 y + x(y)2 + y = x2 (不是二阶微分方程)
定理 1 如果函数 y1 与 y2 是线性齐次方程的 两个解， 则函数 y = C1 y1 + C2 y2
解 求 y = 3e – x – xe – x 的导数， 得 y = - 4e – x + xe - x, y = 5e – x - xe - x,
有 将 y，y 及 y 代入原方程的左边， (5e – x - xe - x) + 2(- 4e – x + xe - x) + 3e – x – xe – x = 0， 即函数 y = 3e – x – xe – x 满足原方程， 所以该函数是 所给二阶微分方程的解.
微分方程的基本概念与解法
THANK YOU
汇报人：XX
适用范围：对于某些复杂的微分方程，通过代换可以将方程转化为更易 于解决的形式 步骤：选择适当的代换变量，将原方程中的未知函数和其导数表示为代 换变量的函数，从而简化方程 举例：对于形如dy/dx=f(x/y)的微分方程，可以通过令y=xu来将其转 化为关于u和x的方程
积分因子法
定义：积分因子 是使微分方程左 边成为全导数的 因子
添加标题
添加标题
添加标题
添加标题
判定方法：通过求解微分方程的解， 分析解的性质，如解的收敛性和稳 定性
应用：非线性微分方程的稳定性在 物理学、工程学、经济学等领域有 广泛应用
稳定性判据
定义：稳定性是指微分方程的解在初始条件下的行为 判据：如果微分方程的解在初始条件下的行为是收敛的，则称该解是稳定的 判断方法：通过分析微分方程的解的性质，如导数的符号等，来判断解的稳定性 应用：稳定性理论在物理学、工程学等领域有广泛应用
优缺点：步进法简单易行，但精度不易控制，需要选择合适的步长和迭代公式
微分方程的稳定性
线性微分方程的稳定性
定义：如果一个线性微分方程的解在某 个初始条件下保持恒定或随时间有规律 地变化，则称该微分方程是稳定的。
判别方法：通过计算微分方程的特征根或 利用Routh-Hurwitz定理来判断稳定性。
原理：基于泰勒级数展开， 通过迭代逼近精确解
定义：是一种用于求解常微 分方程初值问题的数值方法
步骤：包括预估、校正和更 新三个步骤
优点：精度高，稳定性好， 适用于多维问题
Hale Waihona Puke 步进法定义：通过逐步逼近的方法求解微分方程的数值解法
原理：将微分方程转化为一系列离散点上的代数方程，逐步求解
第一讲 微分方程的基本概念
第一讲 微分方程的基本概念教学目的：了解微分方程的有关概念难 点：微分方程解的分类与判定重 点：常微分方程、通解与特解、初始条件与初值问题我们先通过具体的例子来说明微分方程的有关概念．例1 设曲线y = f (x )在其上任一点(x ，y )的切线斜率为3x 2，且曲线过点(0，-1)，求曲线的方程．解 由导数的几何意义知在点(x ，y )处，有23x dx dy =． (1) 此外，曲线满足条件 .10-==x y (2) (1)式两边积分，得.332c x dx x y +==⎰ (3)其中c 为任意常数．(3)式表示了无穷多个函 数图(6–1)，为得到满足条件(2)的具体曲线，以条件(2)代入(3)，得c = -1．故所求曲线的方程为.13-=x y (4)例2 质量为m 的物体在离地面高为0s 米处，以初速0v 垂直上抛，设此物体的运动只受重力的影响，试确定该物体运动的路程s 与时间t 的函数关系．解 因为物体运动的加速度是路程s 对时间t 的二阶导数，由于物体运动只受重力的影响，所以由牛顿第二定律知所求函数)(t s s =应满足g dts d -=22． (5) 这里g 为重力加速度，取垂直向上的方向为正方向．此外，)(t s 还应满足条件：00(0),(0).s s s v =⎧⎨'=⎩ (6)(5)式两端对t 积分，得1C gt dtds +-=．(7) 再对t 积分，得21221C t C gt s ++-=． (8) 把条件(6)代入(7)和(8)，得0201,s C v C ==，于是有00221s t v gt s ++-=． (9) 关系式(1)与(5)都含有未知函数的导数，它们都称为微分方程．一般地有 定义1 含有未知函数的导数或微分的方程称为微分方程，未知函数为一元函数的微分方程称为常微分方程．微分方程中出现的未知函数的导数的最高阶数称为这个方程的阶．未知函数为多元函数的微分方程称为偏微分方程．本章只介绍常微分方程，在不致混淆的情况下，也称常微分方程为微分方程或简称为方程．可以看出，方程(1)是一阶微分方程，方程(5)是二阶微分方程．而x y y x y 2sin 4='-''+'''． (10)+'''y x 256)(x y ='． (11)都是三阶微分方程．n 阶常微分方程的一般形式为0),,,;()(=⋅⋅⋅'n y y y x F ． (12)其中x 为自变量，y 为未知函数; ),,,;()(n y y y x F ⋅⋅⋅'是)(,,,n y y y x ⋅⋅⋅'的已知函数，且)(n y 的系数不为0．如果方程(12)的左端函数F 为)(,,,,n y y y y ⋅⋅⋅'''的线性函数，则称方程(12)为n 阶线性微分方程．否则称(12)为非线性的．n 阶线性微分方程的一般形式为 )()()()()1(1)(0x f y x a y x a y x a n n n =+⋅⋅⋅++-． (13) 其中)(),(,),(),(10x f x a x a x a n ⋅⋅⋅均为x 的已知函数，且0)(0≠x a ．例如方程(1)为一阶线性方程，方程(5)是二阶线性方程，而方程(11) 是三阶非线性方程．定义2 如果将已知函数)(x y ϕ=代入方程(12)后，能使其成为恒等式，则称函数)(x y ϕ=是方程(12)的解．如果由关系式0),(=Φy x 确定的隐函数)(x y ϕ=是方程(12)的解，则称0),(=Φy x 为方程(14)的隐式解．为今后叙述简便起见，将对微分方程的解和隐式解都不再加以区别，统称为方程的解．定义3 若微分方程的解中所含(独立的)任意常数的个数与微分方程的阶数相等，则称这个解为方程的通解．在通解中给任意常数以确定的值得到的解，称为微分方程的特解．例如，函数(3)、(8)分别是方程(1)、(5)的通解，函数(4)、(9)分别是方程(1)、(5)的特解，它们都由通解得到．通常，为确定n 阶方程(12)的某个特解，需给出该特解应满足的附加条件，称之为定解条件．一般地，n 阶微分方程应有n 个定解条件，才能从通解中确定某个具体的特解．n 阶微分方程(14)常见的定解条件是如下形式的条件：10)1(1000)(,,)(,)(--=='=n n y x y y x y y x y ．其中1100,,,,-n y y y x 为1+n 个给定的常数，通常称这样的定解条件为初始条件．例如，方程(1)满足初始条件(2)的特解是函数(4)，而方程(5)满足初始条件(6)的特解是函数(9)．求微分方程满足某定解条件的解的问题，称为微分方程的定解问题; 求微分方程满足某初始条件的解的问题，称为初值问题．例3 验证： 函数at C at C x sin cos 21+=是微分方程x a dtx d 222+= 0． (14) 的通解．解 求出函数at C at C x sin cos 21+=的导数：,cos sin 21at a C at a C dtdx +-= .sin cos 222122at a C at a C dtx d --= 将以上两式代入方程(14)的左端，得(14)．因此，函数at C at C x sin cos 21+=是方程(14)的解，又此函数中含有两个任意常数，而方程(14)为二阶微分方程，因此，函数at C at C x sin cos 21+=是方程(14)的通解．例4 验证： 由方程C y xy x =+-22所确定的隐函数是微分方程y x y y x -='-2)2(． (15) 的解，并求出满足初始条件11==x y 的特解．解 在方程C y xy x =+-22两边对x 求导，得022='+'--y y y x y x ．即y x y y x -='-2)2(．所以由方程C y xy x =+-22所确定的隐函数是微分方程(15)的解． 以初始条件11==x y 代入方程C y xy x =+-22，得1=C ．于是，所求特解为122=+-y xy x ．小结：微分方程的概念：阶、解、通解、特解、初始条件与初值问题．第二讲 一阶微分方程教学目的：掌握常见一阶微分方程的求解方法难 点：一阶线性非齐次微分方程的通解重 点：可分离变量的微分方程、齐次方程和一阶线性微分方程．一阶微分方程的一般形式0),,(='y y x F ，或),(y x f y ='．本节将介绍某些特殊类型的一阶微分方程的解法，包括可分离变量的微分方程、齐次方程和一阶线性微分方程．1．可分离变量的微分方程如果一阶微分方程能化为dx x M dy y N )()(=(1)的形式，那么原方程称为可分离变量的微分方程．要解这类方程，先把原方程化为(1)式的形式，称为分离变量，再对(1)式两边积分，得⎰⎰=dx x M dy y N )()(，便可得到所求的通解．如果需要求其特解，可由初始条件00y y x x ==代入通解中定出任意常数C 的值，即可得到相应的特解．例1 求解微分方程xy dxdy 2=． 解 原微分方程可以分离变量，分离变量后得xdx dy y21=． 两边积分 ⎰⎰=xdx dy y21． 12ln C x y +=． 2112x c C x e e e y ⋅==+．21x C e e y ⋅±=．因为1C e ±仍是任意常数，把它记作C ，便得原方程的通解为2x Ce y =． 以后为了运算方便起见，把y ln 写成y ln ，以上解答过程简写为：.ln ln 2C x y += 2x Ce y =．只要记住最后得到的任意常数C 可正可负即可．例2 求微分方程0)1()1(22=+-+dy x xy dx y满足初始条件2)1(=y 的特解．解 分离变量，得dx x x dy y y )1(1122+=+． 即dx x x x dy y y ⎪⎭⎫ ⎝⎛+-=+22111． 两边积分，得 C x x y ln 21)1ln(21ln )1ln(2122++-=+． 即 )ln(1)(1ln(222Cx y x =++）． 因此，通解为222)1)(1(Cx y x =++．这里C 为任意常数．把初始条件2)1(=y 代入通解，可得10=C ．于是，所求特解为22210)1)(1(x y x =++．例3 实验得出，在给定时刻t ，镭的衰变速率(质量减少的即时速度)与镭的现存量M = M (t )成正比．又当t = 0时，M = M 0，求镭的存量与时间t 的函数关系．解 依题意，有.0),()(>-=k t kM dt t dM (2) 并满足初始条件.00M M t ==方程(2)是可分离变量的，分离变量后得kdt MdM -=． 两边积分，得C kt M ln ln +-=．即kt Ce M -=． 将初始条件00M M t ==代入上式，得0M C =，故镭的衰变规律可表示为.0kt e M M -=一般地，利用微分方程解决实际问题的步骤为：① 利用问题的性质建立微分方程，并写出初始条件；② 利用数学方法求出方程的通解；③ 利用初始条件确定任意常数的值，求出特解．2．齐次方程可化为形如 ⎪⎭⎫ ⎝⎛=x y f dx dy ． (3) 的微分方程，称为一阶齐次微分方程，简称为齐次方程．例如方程0)2()(22=---dy xy x dx y xy可化为 ⎪⎭⎫ ⎝⎛-⎪⎭⎫ ⎝⎛-=--=x y x y x y xy x y xy dx dy 212222． 它是一阶齐次微分方程．一般地，形如0),(),(=+dx y x N dy y x M 的方程，若),(y x M 与),(y x N 均为x ，y 的m 次齐次函数，则它是可化为形如(3)的齐次方程．齐次方程(3)中的变量x 与y 一般是不能分离的，如果作变量替换xy u =． (4) 就可以把方程(3)化为可分离变量的方程，这是因为ux y =，dxdu x u dx dy +=．将其代入方程(3)，便得)(u f dxdu x u =+． 这是变量可分离的方程，分离变量，并两边积分，得dx x du u u f ⎰⎰=-1)(1． (5) 求出积分后，将u 还原成xy ，便得所给齐次方程的通解． 例4 解微分方程 .tan 2xy x y y =-' 解 原方程可写成 .tan2x y x y y +=' 这是齐次方程．令xy u =，f (u ) = 2tan u + u ．代入(5)得 .tan 2⎰⎰=x dx u du积分得.ln ln ln 2sin ln 2cx c x u =+=.sin 2cx u = 将xy u =代入上式，便得原方程的通解为 .sin 2cx xy = 在微分方程中，一般习惯上把x 看作自变量，但有时若将y 看作自变量，求解时会很简便，如下例．例5 求微分方程023(22=--xydx dy x y ）． 满足初始条件10==x y 的特解．解 原方程可化为y x y x xy x y dy dx ⋅⎪⎪⎭⎫ ⎝⎛-=-=23123222．令yx u =，即uy x =，则dy du y u dy dx +=，代入上式，得 uu dy du y 2512-=． 分离变量，并两边积分，得dy y du u u ⎰⎰=-15122．注左=⎰---2251)51(51u u d (凑微分)即 C y u ln 51ln )51ln(512-=--． 将yx u =代入，得到原方程的通解为 C y x y =-3255 将初始条件10==x y 代入通解中，得到1=C ．于是，所求特解为15325=-y x y ．与齐次方程类似，某些微分方程通过变量替换可化为可分离变量的方程，然后分离变量，经积分可求得通解．变量替换的方法是解微分方程最常用的方法．在后面，我们还会用到这种方法，这里再举一例．例6 求解微分方程11+-=yx dx dy ． 解 令u y x =-，则u x y -=，dx du dx dy -=1，于是 111+=-u dx du ．udx du 1-=． 分离变量，并两边积分，得 C x u +-=22．以y x u -=代回，得C x y x +-=-2)(2．3．一阶线性微分方程可化为形如)()(x Q y x P dxdy =+． (6) 的微分方程，称为一阶线性微分方程，其中)(),(x Q x P 均为x 的已知函数．当0)(≡x Q 时，称方程(6)是齐次的; 当)(x Q 不恒为零时，称方程(6)是非齐次的．设方程(6)是线性非齐次微分方程，把)(x Q 换成零而写出0)(=+y x P dx dy ． (7) 称为对应于方程(6)的线性齐次微分方程．方程(7)是可分离变量的，分离变量后，得dx x P ydy )(-=． 两边积分，得 C dx x P y ln )(ln +-=⎰．于是，方程(7)的通解为⎰=-dx x P Ce y )(． (8)下面求方程(6)的通解．由于方程(7)是(6)的特殊情况，那么方程(6)的通解中必包含着方程(7)的通解．它们的解之间必有某种内在联系，下面我们分析一下方程(6)的解的形式．把方程(6)改写为dx y x Q x P y dy ⎪⎪⎭⎫ ⎝⎛+-=)()(． 两边积分，得1ln )()(ln C dx yx Q dx x P y ++-=⎰⎰． 即⎰⋅⎰=-dx x P dx y x Q e eC y ))(1（． 因为积分dx yx Q ⎰)(中的被积函数含有未知函数y ，因此还不能说得到了方程(6)的解．但是，由于y 是x 的函数，则积分dx yx Q ⎰)(的结果是x 的函数．故可设 )()(1x C eC dx y x Q =⎰．从而有⎰=-dx x P e x C y )()(． (9) 再求未知函数)(x C ．因为(9)式是方程(6)的解，所以(9)式应满足方程(6)，将y 及它的导数⎰-⎰'='--dx x P dx x P e x P x C e x C y )()()()()(． 代入方程(6)，得)()()()()()()()()(x Q e x C x P e x P x C e x C dx x P dx x P dx x P =⎰+⎰-⎰'---．即)()()(x Q e x C dx x P =⎰'-．⎰='dx x P e x Q x C )()()(． 两边积分，得C dx e x Q x C dx x P +⎰=⎰)()()(.把上式代入(9)式，便得方程(6)的通解为⎪⎭⎫ ⎝⎛+⎰⎰=⎰-C dx e x Q e y dx x P dx x P )()()(． (10) 这种将线性齐次方程(7)的通解(8)中的任意常数换成待定函数)(x C ，然后求得线性非齐次方程(6)的通解的方法，叫做常数变易法．将(10)式写成两项之和⎰⎰⎰+⎰=--dx e x Q e Ce y dx x P dx x P dx x P )()()()(． 上式右端第一项是对应的线性齐次方程(7)的通解，第二项是线性非齐次方程(6)的一个特解(即在通解(10)中令0=C ，便得此特解)．因此，一阶线性非齐次方程的通解等于对应的线性齐次方程的通解与线性非齐次方程的一个特解之和．例7 求解微分方程x x x y y sin 2cot =-'．解法1 常数变易法对应齐次方程为.0cot =-'x y y分离变量，得.cot 1xdx dy y =两边积分，得.sin sin ln cot x C Ce Ce y x xdx ⋅==⎰=用常数变易法，把C 换成新的未知函数)(x C ，即令.sin )(x x C y =则.cos )(sin )(x x C x x C y +'='代入原非齐次方程，得x x C 2)(='．两边积分，得C x x C +=2)(．故所求通解为.sin )(2x C x y +=解法2 公式法.sin 2)(,cot )(x x x Q x x P =-=故).(sin )2(sin )sin 1sin 2(sin )sin 2()sin 2(2sin ln sin ln cot cot C x x C xdx x C dx xx x x C dx e x x e C dx xe x e y x x xdx xdx +⋅=+⋅=+⋅⋅=+⋅=+⎰⎰=⎰⎰⎰⎰-- 例8 求微分方程 02)6(2=+'-y y x y 满足初始条件12==x y 的特解． 解 这个方程不是未知函数y 与y '的线性方程，但是可以将它变形为yy x dy dx 262-=． 即23y x y dy dx -=-． (11) 若将x 视为y 的函数，则对于)(y x 及其导数dydx 而言，方程(11)是一个线性方程，由通解公式(10)得 ⎪⎪⎭⎫ ⎝⎛+⎰⎪⎭⎫ ⎝⎛-⎰=⎰-C dy e y e x dy y dy y 332⎪⎪⎭⎫ ⎝⎛+=C y y 213． 以条件2=x 时，1=y 代入，得23=C . 因此，所求特解为 2232y y x +=． 例9 求解微分方程.)(ln 2y x a xy dx dy =+ 解 原方程不是线性方程，但通过适当的变换，可将它化为线性方程．将原方程改写为.ln 112x a y xdx dy y =+-- 即 .ln 111x a y xdx dy =+--- 令1-=y z ，则上式变为.ln 1x a z xdx dz -=- 这是z 关于x 的一阶线性方程．由通解公式(10)，得通解].)(ln 2[2x a C x z -= 所以，原方程通解为.1])(ln 2[2=-x a C xy 一般地，形如 n y x Q y x P dxdy )()(=+ ( 1,0≠n ) (12) 的方程，称为伯努利方程．这类方程可经过变换化为线性方程，方程(12)两边同除以n y 得)()(1x Q y x P dxdy y n n =+--． 再令n y z -=1，则上式化为)()(11x Q z x P dxdz n =+-． 即)()1()()1(x Q n z x P n dxdz -=-+． 这是函数z 关于x 的一阶线性方程，从而可用常数变易法或公式法得出z ，再用n y -1代换z ，即得伯努利方程(12)的解．小结：1．可分离变量的方程：x x f y y g d )(d )(=，两边积分得通解．2．一阶齐次方程：)(x yy ϕ='，令u xy =，得⎰⎰=-x x u u u d )(d ϕ． 注 形如)(c by ax f y ++='的方程可令u c by ax =++转化为可分离变量的方程．3．一阶线性方程：)()(x Q y x P y =+'的通解为]d e )([d )(d )(C x x Q e y x x P x x P +⎰⎰=⎰-． 4．伯努利方程：n y x Q y x P y )()(=+'，令u yn =-1可转化为一阶线性方程．第三讲 可降阶的高阶微分方程教学目的：掌握三种可以降阶的微分方程求解方法重 点：第二类可以降阶的微分方程难 点：第三类可以降阶的微分方程从这节起我们讨论二阶和高于二阶的微分方程，这类方程称为高阶微分方程．有些高阶微分方程可以通过代换化成较低阶的方程来求解．以二阶微分方程而论，如果我们能设法作代换把它从二阶降至一阶，那么就有可能用第二节所讲的方法来求解．下面介绍三种容易降阶的高阶微分方程的求解方法．1．)()(x f y n = 型的微分方程微分方程)()(x f y n = (6-18)的右端仅含有自变量x ，对于这种方程，两端积分便使它降为一个1-n 阶的微分方程()11d )(C x x f y n +=⎰-．再积分可得()[]212d d )(C x C x x f y n ++=⎰⎰-． 依此继续下去，连续积分n 次，便得方程(6-18)的含有n 个任意常数的通解．例 1 求微分方程x y x cos e 2-='''的通解．解 对所给方程连续积分三次，得12sin e 21C x y x +-=''， 212cos e 41C x C x y x +++='， 3221221sin e 81C x C x C x y x ++++=， 这就是所求的通解．2．()y x f y '='',型的微分方程微分方程()y x f y '='', (6-19)中不显含未知函数y ．如果设p y ='，则p xp y '==''d d ，方程(6-19)变成 ),(p x f p ='．这是关于x 和p 的一阶微分方程，设其通解为()1,C x p ψ=． 由于xy p d d =，因此又得到一个一阶微分方程 ()1,d d C x xy ψ=． 对它积分即得(6-19)的通解 ()21d ,C x C x y +=⎰ψ．例 2 求方程()1212='+''+y x y x 的通解．解 所给方程不显含变量y ，令y p '=，则p y '=''，代入原方程得()1212=+'+xp p x ． 它是一阶线性微分方程，化为标准形式221112xp x x p +=++'， 其通解为⎪⎪⎭⎫ ⎝⎛⎰++⎰=⎰++-x x C p x x x x x x d e 11e d 1221d 1222 ()⎥⎦⎤⎢⎣⎡++++=⎰x x x C x d 111112212 211xC x ++=． 将y p '=代入上式，并再积分一次得所求方程的通解()212arctan 1ln 21C x C x y +++=． 例 3 求方程()y x x y '=+''212满足初始条件1|0==x y ，3|0='=x y 的特解． 解 此方程不显含y ，令y p '=，则p y '=''，代入方程得()xp x p 212=+'．分离变量后两边积分得()211x C p +=， 由3|0='=x y 得31=C ，从而()213d d x x y +=． 两边积分得233C x x y ++=，由1|0==x y 得12=C ．故所求特解为133++=x x y ．3．()y y f y '='',型的微分方程微分方程()y y f y '='',中不显含自变量x ，对于这类方程，令y p '=，两边对x 求导得yp p x y y p x p y d d d d d d d d =⋅==''． 则方程(6-20)变成),(d d p y f yp p =． 这是一个关于变量p 和y 的一阶微分方程，设它的通解为()1,C y p y ϕ=='．分离变量并积分，即可得方程(6-20)的通解()21d ,y x C y C =+ϕ⎰．例4 求微分方程()02='-''y y y 的通解． 解 方程中不显含自变量x ，设p y ='，则y p py d d =''，代入原方程得 0d d 2=-p y p yp．如果0≠p ，那么方程中约去p 并分离变量得yy p p d d =． 两端积分并化简，得y C p 1=，即y C y 1='．再分离变量并积分，得21ln ln C x C y +=，即x C C y 1e 2=．如果0=p ，那么C y =，显然它也满足原方程，但C y =已包含在上述解中(令01=C 即得)，所以原方程的通解为x C C y 1e 2=．小结：1．)()(x f y n =型，连续积分n 次，便得方程的含有n 个任意常数的通解．2．()y x f y '='',型（不显含未知函数y ）．令p y ='，方程变成),(p x f p ='．3．()y y f y '='',型（不显含自变量x ）.，令y p '=，方程变成),(d d p y f yp p =．第四讲 二阶常系数线性微分方程教学目的：掌握二阶常系数线性方程的求解方法重 点：二阶常系数齐次和非齐次线性方程的求解难 点：二阶常系数非齐次线性方程的特解二阶常系数线性微分方程的一般形式为)(x f qy y p y =+'+''．这里p 、q 是常数，)(x f 是x 的已知函数．当()f x 恒等于零时，称为二阶常系数齐次线性微分方程，否则称为二阶常系数非齐次线性微分方程．1．二阶常系数齐次线性微分方程定理1 设)(1x y y =与)(2x y y =为二阶常系数齐次线性微分方程0=+'+''qy y p y(1)的相互独立的两个特解(即)()(12x y x y 不恒等于常数)，则2211y C y C y +=为方程(1)的通解，这里1C 与2C 为任意常数．证 按假设)(1x y 与)(2x y 为方程(1)的解，所以有下式成立0111=+'+''qy y p y ，0222=+'+''qy y p y ． 又 2211y C y C y +=， 2211y C y C y '+'='， 2211y C y C y ''+''=''． 代入(1)式左端，得()()()221122112211y C y C q y C y C p y C y C qy y p y ++'+'+''+''=+'+'' 0)()(22221111=+'+''++'+''=qy y p y C qy y p y C ． 即2211y C y C y +=为方程(1)的解． 在)()(12x y x y 不恒等于常数的条件下，2211y C y C y +=中含有两个相互独立的任意常数1C 和2C ，所以2211y C y C y +=是方程(1)的通解．由此定理可知，求方程(1)的通解问题，归结为求(1)的两个相互独立的特解．为了寻找这两个特解，注意到当r 为常数时，指数函数rx y e =和它的各阶导数只相差一个常数因子，因此不妨用rx y e =来尝试．设rx y e =为方程(1)的解，则rx r y e ='，rx r y e 2=''，代入方程(1)得.0)(2=++rx e q pr r由于0e ≠rx ，所以有.02=++q pr r (2) 只要r 满足(2)式，函数rx y e =就是微分方程(1)的解．我们把代数方程(2)称为微分方程(1)的特征方程，特征方程的根称为特征根．由于特征方程是一元二次方程，故其特征根有三种不同的情况，相应地可得到微分方程(1)的三种不同形式的通解．(ⅰ) 当042>-q p 时，特征方程(8-23)有两个不相等的实根1r 和2r ，此时可得方程(1)的两个特解：x r y 1e 1=， x r y 2e 2=，且≠=-x r r y y )(1212e /常数，故x r x r C C y 21e e 21+=是方程(1)的通解．(ⅱ) 当042=-q p 时，特征方程(8-23)有两个相等的实根21r r =，此时得微分方程(1)的一个特解x r y 1e 1=．为求(1)的通解，还需求出与x r 1e 相互独立的另一解2y ．不妨设)(/12x u y y =，则)(e 12x u y x r =， )(e 121u r u y x r +'='， )2(21121u r u r u e y x r +'+''=''． 将22,y y '及2y ''代入方程(1)，得 0])()2[(e 12111=++'++'+''qu u r u p u r u r u x r ．将上式约去x r 1e 并合并同类项，得0)()2(1211=+++'++''u q pr r u p r u ．由于1r 是特征方程(2)的二重根，因此，0121=++q pr r ，且021=+p r ，于是得0=''u ．不妨取x u =，由此得到微分方程(1)的另一个特解x r x y 1e 2=，且≠=x y y 12/常数，从而得到微分方程(1)的通解为x r x r x C C y 11e e 21+=，即)(e 211x C C y x r +=．(ⅲ) 当042<-q p 时，特征方程(2)有一对共轭复根βαi r +=1，βαi r -=2．于是得到微分方程(1)的两个特解x i y )(1e βα+=，x i y )(2e βα-=．但它们是复数形式，为应用方便，利用欧拉公式θθθsin cos e i i +=将1y 和2y 改写成)sin (cos e 1x i x y x ββα+=，)sin (cos e 2x i x y x ββα-=．于是得到两个新的实函数x y y y x βαcos e )(21211=+=， x y y iy x βαsin e )(21212=-=． 可以验证它们仍是(1)的解，且≠=x y y βtan /12常数，故微分方程(1)的通解为)sin cos (e 21x C x C y x ββα+=．综上所述，求微分方程(1)通解的步骤可归纳如下：第一步 写出微分方程(1)的特征方程02=++q pr r ，求出特征根； 第二步 根据特征根的不同形式，按照下表写出微分方程(1)的通解： 表1 特征方程02=++q pr r 的根21,r r 微分方程0'''=++qy py y 的通解两个不等实根21r r ≠ x r x r C C y 21e e 21+=两个相等实根21r r = x r x C C y 1e )(21+=一对共轭复根βαi r ±=2,1 )sin cos (e 21x C x C y x ββα+=例 1 求微分方程043=-'+''y y y 的通解．解 所给微分方程的特征方程为0432=-+r r ．特征根为121, 4.r r ==- 于是，所求微分方程的通解为x x C C y 421e e -+=．例 2 求微分方程044=+'-''y y y 的满足初始条件1|,1|00='===x x y y 的特解．解 所给微分方程的特征方程为0442=+-r r ．特征根221==r r ．故所求微分方程的通解为)(e 212x C C y x +=．求导得x x C x C C y 22212e )(e 2++='．将初始条件1|0==x y 及1|0='=x y 代入以上两式求得.1,121-==C C 故所求特解为)1(e 2x y x -=．例 3 设函数)(x f 可导，且满足⎰⎰-++=xx t t f x t t tf x x f 00d )(d )(21)(． 试求函数)(x f .解 由上述方程知(0)1f =．方程两边对x 求导得⎰-='xt t f x f 0d )(2)(． 由此可得(0)2f '=．上式两边再对x 求导得)()(x f x f -=''．这是二阶常系数齐次线性方程，其特征方程为,012=+r特征根.,21i r i r =-= 于是，所求微分方程的通解为12()cos sin .f x C x C x =+由此得.cos sin )(21x C x C x f +-='由(0)1f =，(0)2f '=得.2,121==C C 所以.sin 2cos )(x x x f +=本节介绍的求二阶常系数齐次线性微分方程通解的原理和方法，也可以用于求解更高阶的常系数齐次线性方程．例 4 求四阶微分方程08)4(='+y y 的通解．解 所给微分方程的特征方程为084=+r r ，即,0)42)(2(2=+-+r r r r 其特征根为.31,2,04,321i r r r ±=-= 于是得方程的通解).3sin 3cos (e e 43221x C x C C C y x x +++=-2．二阶常系数非齐次线性微分方程从第二节的讨论知，一阶非齐次线性微分方程的通解等于对应的齐次线性方程的通解与非齐次线性方程的一个特解之和．而二阶常系数非齐次线性微分方程具有相类似的性质．定理2 设()y y x **=是二阶常系数非齐次线性微分方程)(x f qy y p y =+'+''(3)的一个特解，而Y 为对应于方程(3)的齐次线性微分方程的通解，则y Y y *=+为方程(3)的通解．由此结论可知，二阶常系数非齐次线性微分方程的通解，可按下面三个步骤来求：①求其对应的齐次线性微分方程的通解Y ；②求非齐次线性微分方程的一个特解y *；③原方程的通解为y Y y *=+．求齐次线性微分方程的通解Y 的方法前面已讨论过，所以只要研究一下如何求非齐次方程(3)的一个特解就行．限于篇幅, 这里只讨论)(x f 为以下两种形式的情形．I.,e )()(x m x P x f λ=其中λ是常数，)(x P m 是x 的m 次多项式：m m m m m a x a x a x a x P ++++=--1110)( ；II ．[]()e ()cos ()sin x t n f x P x x P x x λ=ω+ω，其中λ和ω是常数，)(x P t 、)(x P n 分别是x 的t 次和n 次多项式，其中有一个可为零．对于以上两种情形，下面用待定系数法来求方程(3)的一个特解，其基本思想是：先根据)(x f 的特点，确定特解y *的类型，然后把y *代入到原方程中，确定y *中的待定系数．I.x m x P x f λe )()(=型因为方程(3)右端)(x f 是多项式)(x P m 与指数函数x λe 的乘积，而多项式与指数函数乘积的导数仍然是同一类型的函数，因此，我们推测()e x y Q x *λ=(其中)(x Q 是某个多项式)可能是方程(3)的一个解，把y *、()y *'及()y *''代入方程(3)，求出)(x Q 的系数，使()e x y Q x *λ=满足方程(3)即可．为此将()e x y Q x *λ=，[]()e ()()x y Q x Q x *λ''=λ+，2()e ()2()()x y Q x Q x Q x *λ'''''⎡⎤=λ+λ+⎣⎦代入方程(3)并消去x λe ，得)()()()()2()(2x P x Q q p x Q p x Q m =+++'++''λλλ． (4))1( 如果λ不是方程(3)的特征方程02=++q pr r 的根，由于)(x P m 是一个m 次多项式，要使方程(4)的两端恒等，可令)(x Q 为另一个m 次多项式)(x Q m ，即设)(x Q m 为,)(1110m m m m m b x b x b x b x Q ++++=--其中m b b b ,,,10 为待定系数，将)(x Q m 代入(4)，比较等式两端x 同次幂的系数，可得含有m b b b ,,,10 的1+m 个方程的联立方程组，解出),1,0(m i b i =得到所求特解m )2( 如果λ是特征方程02=++q pr r 的单根，即02=++q p λλ，但,02≠+p λ 要使(4)式的两端恒等，)('x Q 必须是m 次多项式，此时可令),()(x xQ x Q m =并且可用同样的方法确定)(x Q m 的系数),1,0(m i b i =．)3( 如果λ是特征方程02=++q pr r 的重根，即02=++q p λλ且,02=+p λ要使式(4)的两端恒等，)(x Q ''必须是m 次多项式，此时可令),()(2x Q x x Q m =并且利用同样的方法可以确定)(x Q m 的系数),1,0(m i b i =．综上所述，我们有以下结论：如果x m x P x f λe )()(=，则二阶常系数非齐次线性微分方程(3)具有形如()e k x m y x Q x *λ=的特解；其中)(x Q m 是与)(x P m 同次(m 次)的多项式，而k 按λ不是特征方程的根、是特征方程的单根或是特征方程的重根依次取为0、1或2．例 5 求方程2106652+-=+'-''x x y y y 的通解．解 所给方程是二阶常系数非齐次线性微分方程，且右端函数形如x m x P λe )(，其中,0=λ.2106)(2+-=x x x P m先求对应齐次方程065=+'-''y y y的通解，其特征方程是0652=+-r r ．特征根,3,221==r r 对应齐次方程的通解为x x C C Y 3221e e +=．因为0=λ不是特征根，因而所求方程有形如的特解．由于()2,y Ax B *'=+()2,y A *''= 将它们代入原方程中得恒等式.2106652)106(622+-=+-+-+x x C B A x A B Ax比较上式两端x 的同次幂的系数可得⎪⎩⎪⎨⎧=+--=-=.2652,10106,66C B A A B A解方程组得.0,0,1===C B A 故所求方程的一个特解为2.y x *=从而所求方程的通解为.e e 23221x C C y x x ++=例 6 求方程x x y y y 2e 24'4=+-''的通解．解 所求方程是二阶常系数非齐次线性微分方程，且右端函数形如x m x P λe )(，其中,2=λ.2)(x x P m =所求解的方程对应的齐次方程044=+'-''y y y 的通解为).(e 212x C C Y x +=由于2=r 是二重特征根，所以设所求方程有形如22()e x y x Ax B *=+的特解．将它代入所求方程可得622.Ax B x +=比较等式两端x 的同次幂的系数，得0,31==B A ．于是得所求方程的一个特解为321e .3x y x *= 最后得所求方程的通解为).31(e 3212x x C C y x ++=II ．]sin )(cos )([e )(x x P x x P x f n l x ωωλ+=型可以推证，如果]sin )(cos )([e )(x x P x x P x f n l x ωωλ+=，则二阶常系数非齐次线性微分方程(3)的特解可设为e [()cos ()sin ],k x m m y x Q x x R x x *λ=ω+ω其中),(x Q m )(x R m 是m 次多项式，},,max{n l m = 而k 按ωλi ±不是特征方程的根或是特征方程的单根依次取0或1．例 7 求方程)sin 7(cos e 2x x y y y x -=-'+''的通解． 解 所求解的方程对应的齐次方程02=-'+''y y y 的特征方程为022=-+r r ，特征根2,121-==r r ，齐次方程的通解为x x C C Y 221e e -+=．因为i i ±=±1ωλ不是特征根，故所求方程具有形如(cos sin )x y e A x B x *=+的特解，求得()e [()cos ()sin ]x y A B x B A x *'=++-，()e [2cos 2sin ]x y B x A x *''=-．代入所求方程并化简得恒等式.sin 7cos sin )3(cos )3(x x x A B x A B -=+--比较上式两端x cos 和x sin 的系数，可得⎩⎨⎧-=--=+-.73,13B A B A 因此,1,2==B A 故e (2cos sin ).x y x x *=+所求通解为.e e )sin cos 2(221x x x C C x x e y -+++=小结：1．特征方程02=++q pr r 的根21,r r 微分方程0'''=++qy py y 的通解两个不等实根21r r ≠ x r x r C C y 21e e 21+= 两个相等实根21r r = xr x C C y 1e )(21+=一对共轭复根βαi r ±=2,1 )sin cos (e 21x C x C y x ββα+=2．x m x P x f λe )()(=型，特解形式()e k x m y x Q x *λ=，其中)(x Q m 是与)(x P m 同次(m 次)的多项式，而k 按λ不是特征方程的根、是特征方程的单根或是特征方程的重根依次取为0、1或2．3．]sin )(cos )([e )(x x P x x P x f n l x ωωλ+=型，特解形式e [()cos ()sin ],k x m m y x Q x x R x x *λ=ω+ω其中),(x Q m )(x R m 是m 次多项式，},,max{n l m = 而k 按ωλi ±不是特征方程的根或是特征方程的单根依次取0或1．4．二阶常系数非齐次线性微分方程的通解，可按下面三个步骤来求： ①求其对应的齐次线性微分方程的通解Y ；②求非齐次线性微分方程的一个特解y *；③原方程的通解为y Y 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 一阶线性偏微分方程
常微分方程的解的表达式中，可能包含一个或者几个常
精品课件-常微分方程(王素云)-第1章
第1章 基 本 概 念
除了初值条件外，另外一种常见的定解条件是边值条件. 最后，我们对n阶微分方程的通解关于n个任意常数的独立 性作一点说明. 一个n阶微分方程的通解包含n个独立的任意常数. 反 之，设y=g(x,C1,C2,…,Cn)是充分光滑的函数族，其中x是自 变量，而C1,C2,…,Cn是n个独立的参数(任意常数),则存在一 个形如式(1.1)的n阶微分方程，使得它的通解恰好是上面的 函数族y=g(x,C1,C2,…,Cn). 我们把这个一般结论的证明留给读者(习题1.1的第4题)， 它的证明方法与例1.8的讨论是类似的.
b f 2 (x) d x 1 [b2 f (b) f (a)]
a
3
第1章 基 本 概 念
上式对于每点b(b>a)都成立，两边对b求导得
f (b) 3 f 2 (b) 2 f (b)
b2
b
改用惯用的符号:
f
( x)
3
f
2 (x) x2
2
f
(x) x
这就是所要建立的微分方程.
第1章 基 本 概 念
v0t
y0
（1.8）
因此它描述了具有初始高度y0和初始速度v0的自由落体运动.
我们称式(1.8)是初值问题式(1.4)与(1.7)的解，亦即初
值问题:
y" g ,
y(0)
y0
,
y'(0) v0
（1.9）
的解.初值问题又叫柯西问题.
第1章 基 本 概 念
再看一例，一曲线通过点(1,2)
M(x,y) 处的切线斜率为2x，求该曲线的方程.
(3) yy ( y)2 1 0,
(4)
d2
dt2
1.1-4微分方程概念及一阶微分方程
此外, s( t ) 还应满足条件 : s(0) = 0, s′(0) = 0.
4
1.1-2
微分方程的初等积分法
求此曲线方程 例 3 设一曲线的曲率半径等 于常数 a ,求此曲线方程 .
解 : 设所求曲线为 y = y ( x ), 则由曲率半径的公式得
(1 + y′ ) y′′
2
3
2
= a,
即
(1 + y′ ) 2 = ± ay′′
12
1.1-2
微分方程的初等积分法
（2） y = C1 cos 3 x + C2 sin 3 x 。
解： y ′ = −3C1 sin 3 x + 3C 2 cos 3 x ，
y ′′ = − 9C 1 cos 3 x − 9C 2 sin 3 x = − 9(C 1 cos 3 x + C 2 sin 3 x ) ，
2
3
由此式解得 y ( x ), 便是所求曲线 .
5
1.1-2
微分方程的初等积分法
1.1 微分方程的基本概念
1.微分方程的定义 1.微分方程的定义
含有自变量、未知函数及未知函数的导数(或微分) 含有自变量、未知函数及未知函数的导数(或微分)的 等式称为微分方程。 等式称为微分方程。
未知函数是一元函数的微分方程称为常微分方程， 未知函数是一元函数的微分方程称为常微分方程，未 常微分方程 知函数是多元函数的微分方程称为偏微分方程 偏微分方程。 知函数是多元函数的微分方程称为偏微分方程。
2
1.1-2
微分方程的初等积分法
例 2 质量为 的物体从静止开始自由 下落, 求物体下落 质量为m
的运动规律 .
微分方程-1
,
江西理工大学理学院 例9 小船从河边 点 0 处 出发驶向对岸(两岸为平 行直线).设 船速为 a ,船行方向始终与河岸垂直, 设河宽 为 h ,河中任意点处的水流速度与该点到 两岸距离的乘积成正比(比例 系数为 k ).求小船 的航行路线 . y 解 设 0 为原点,河岸朝顺水方向为 x 轴 ,
2
江西理工大学理学院
− k 2 (C1 cos kt + C 2 sin kt ) + k 2 (C1 cos kt + C 2 sin kt ) ≡ 0.
故 x = C1 cos kt + C 2 sin kt 是原方程的解 .
Q x t =0
dx = A, = 0, dt t = 0
∴ C1 = A, C 2 = 0.
江西理工大学理学院
第 十 一 章
微 分 方 程
江西理工大学理学院
第 一 节 微分方程基本概念 可分Biblioteka 变量方程江西理工大学理学院
一、问题的提出
例 1 一 曲 线 通 过 点 (1,2), 且 在 该 曲 线 上 任 一 点
M ( x , y ) 处的切线的斜率为 2 x ,求这曲线的方程.
解
设所求曲线为 y = y( x )
ln y = x + C1
2
∴ y = Ce 为所求通解 .
x2
江西理工大学理学院
例5 求方程 f ( xy ) ydx + g ( xy ) xdy = 0 通解 .
解 令u = xy ,
则 du = xdy + ydx , du − ydx f ( u) ydx + g ( u) x ⋅ = 0, x u [ f ( u) − g ( u)] dx + g ( u)du = 0, x dx g ( u) + du = 0, x u[ f ( u) − g ( u)] g ( u) 通解为 ln | x | + ∫ du = C . u[ f ( u) − g ( u)]
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
微 分 方 程
辽宁工程技术大学理学院 应用数学系
求微分方程的数值解
（一）常微分方程数值解的定义
（二）建立数值解法的一些途径
（三）用Matlab软件求常微分方程的数值解
返 回
数学建模实例
1、目标跟踪问题一：导弹追踪问题 2、目标跟踪问题二：慢跑者与狗 3、地中海鲨鱼问题
返 回
（1）
2. 由于弹头始终对准乙舰，故导弹的速度平行于乙舰与导弹头位置的差向量，
（2）
（3）
3．因乙舰以速度v0沿直线x=1运动，设v0=1，则w=5，X=1，Y=t
因此导弹运动轨迹的参数方程为：
4. 解导弹运动轨迹的参数方程 建立m-文件eq2.m如下： function dy=eq2(t,y) dy=zeros(2,1); dy(1)=5*(1-y(1))/sqrt((1-y(1))^2+(t-y(2))^2); dy(2)=5*(t-y(2))/sqrt((1-y(1))^2+(t-y(2))^2); 取t0=0，tf=2，建立主程序chase2.m如下： [t,y]=ode45('eq2',[0 2],[0 0]); Y=0:0.01:2; plot(1,Y,'-')， hold on To Matlab(chase2) plot(y(:,1),y(:,2),'*')
导弹追踪问题 设位于坐标原点的甲舰向位于x轴上点A(1, 0)处的乙舰 发射导弹，导弹头始终对准乙舰.如果乙舰以最大的速度 v0(是常数)沿平行于y轴的直线行驶，导弹的速度是5v0，求 导弹运行的曲线方程.又乙舰行驶多远时，导弹将它击中？ 解法一（解析法）
假设导弹在 t 时刻的位置为 P(x(t), y(t))，乙舰位于 (1, v0 t ) . Q
To Matlab（ff5）
2、取t0=0，tf=12，输入命令： [T,Y]=ode45('rigid',[0 12],[0 1 1]); plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+')
3、结果如图
返 回
图中，y1的图形为实线，y2的图形为“*”线，y3的图形为“+”线.
则微分方程变为一阶微分方程组：
1、建立m-文件vdp1000.m如下： function dy=vdp1000(t,y) dy=zeros(2,1); dy(1)=y(2); dy(2)=1000*(1-y(1)^2)*y(2)-y(1);
-0.5 -1 -1.5 -2 -2.5 0
y1 ' y2 y2 ' 1000 (1 y12 ) y2 y1 y (0) 2, y (0) 0 1 2
解即为导弹的运行轨迹:
5 5 5 y (1 x) (1 x) 8 12 24
4 5
6 5
5 5 ) 处时被导弹击中. 当x 1时 y ，即当乙舰航行到点(1, 24 24 y 5 被击中时间为: t . 若 v0=1, 则在 t=0.21 处被击中. v 0 24 v 0
在chase3.m，不断修改tf的值,分别取tf=5, 2.5, 3.5,…,至3.15时, 狗刚好追上慢跑者.
(2) w=5时 建立m-文件eq4.m如下: function dy=eq4(t,y) dy=zeros(2,1); dy(1)=5*(10+20*cos(t)-y(1))/sqrt ((10+20*cos(t)-y(1))^2+(20+15*sin(t)-y(2))^2); dy(2)=5*(20+15*sin(t)-y(2))/sqrt ((10+20*cos(t)- y(1))^2+(20+15*sin(t)-y(2))^2);
在chase2.m中，按二分法逐步修改tf，即分别取tf=1， 0.5,0.25,…,直到tf=0.21时，得图2. To Matlab(chase2) 结论：时刻t=0.21时，导弹在（1，0.21）处击中乙舰。 返 回
慢跑者与狗
一个慢跑者在平面上沿椭圆以恒定的速率v=1跑步,设椭 圆方程为: x=10+20cost, y=20+5sint. 突然有一只狗攻击他. 这只狗 从原点出发,以恒定速率w跑向慢跑者,狗的运动方向始终指向慢跑者. 分别求出w=20,w=5时狗的运动轨迹. 1. 模型建立
dx dt X x ， 0 即： Y y dy dt w dx ( X x) dt 2 2 ( X x) (Y y) 消去λ 得： dy w (Y y) 2 2 dt ( X x) (Y y )
解 输入命令 ：
[x,y,z]=dsolve('Dx=2*x-3*y+3*z','Dy=4*x-5*y+3*z','Dz=4*x-4*y+2*z', 't')； x=simple(x) % 将x化简 y=simple(y) To Matlab（ff3） z=simple(z)
结 果 为：x = (c1-c2+c3+c2e -3t-c3e-3t)e2t y = -c1e-4t+c2e-4t+c2e-3t-c3e-3t+c1-c2+c3)e2t z = (-c1e-4t+c2e-4t+c1-c2+c3)e2t 返 回
注意:
1、在解n个未知函数的方程组时，x0和x均为n维向量， m-文件中的待解方程组应以x的分量形式写成. 2、使用Matlab软件求数值解时，高阶微分方程必须 等价地变换成一阶微分方程组.
d 2x dx 2 1000 (1 x 2 ) x 0 例 4 dt dt x (0) 2; x' (0) 0 解: 令 y1=x，y2=y1’
2 1.5 1 0.5 0
500
1000
1500
2000
2500
3000
2、取t0=0，tf=3000，输入命令： [T,Y]=ode15s('vdp1000',[0 3000],[2 0]); plot(T,Y(:,1),'-') 3、结果如图
To Matlab（ff4）
例 5 解微分方程组.
y1 ' y2 y3 y2 ' y1 y3 y3 ' 0.51 y1 y2 y1 (0) 0, y2 (0) 1, y3 (0) 1
例1 求
解
d2y解. dt
0 应表达为：D2y=0.
输入命令：dsolve('Du=1+u^2','t')
结
果：u = tg(t-c)
To Matlab（ff1）
例 2 求微分方程的特解.
d 2 y dy 2 4 29 y 0 dx dx y (0) 0, y ' (0) 15
轨迹图见程序chase1
To Matlab(chase1)
解法二(数值解) 令y1=y,y2=y1’，将方程（3）化为一阶微分方程组。
1 (1 x) y ' ' 1 y2 5
1.建立m-文件eq1.m function dy=eq1(x,y) dy=zeros(2,1); dy(1)=y(2); dy(2)=1/5*sqrt(1+y(1)^2)/(1-x); 2. 取x0=0，xf=0.9999，建立主程序ff6.m如下： x0=0，xf=0.9999 [x,y]=ode15s('eq1',[x0 xf],[0 0]); plot(x,y(:,1),’b.') hold on y=0:0.01:2; To Matlab(ff6) plot(1,y,’b*') 结论: 导弹大致在（1，0.2）处击中乙舰
2. 模型求解
(1) w=20时,建立m-文件eq3.m如下: function dy=eq3(t,y) dy=zeros(2,1); dy(1)=20*(10+20*cos(t)-y(1))/sqrt ((10+20*cos(t)-y(1))^2+(20+15*sin(t)-y(2))^2); dy(2)=20*(20+15*sin(t)-y(2))/sqrt ((10+20*cos(t)-y(1))^2+(20+15*sin(t)-y(2))^2); 取t0=0，tf=10，建立主程序chase3.m如下： t0=0;tf=10; [t,y]=ode45('eq3',[t0 tf],[0 0]); T=0:0.1:2*pi; X=10+20*cos(T); Y=20+15*sin(T); To Matlab(chase3) plot(X,Y,'-') hold on plot(y(:,1),y(:,2),'*')
微分方程的解析解
求微分方程（组）的解析解命令: dsolve(‘方程1’, ‘方程2’,…‘方程n’, ‘初始条件’, ‘自变量’)
记号: 在表达微分方程时，用字母 D 表示求微分，D2、D3 等 表示求高阶微分.任何 D 后所跟的字母为因变量，自变量可以指 定或由系统规则选定为确省. 例如，微分方程
5 dx (1 x) dt 2 2 (1 x) (t y ) 5 dy (t y ) 2 2 (1 x) (t y ) dt x(0) 0, y (0) 0