一阶常微分方程初值问题
微分方程中的初值问题理论
微分方程中的初值问题理论微分方程是数学领域中的重要分支，它描述了一种变量与其变化率之间的关系。
在实际问题中，经常会遇到需要确定微分方程的解的具体形式，并以给定的初值条件作为起点进行求解的情况，这就是初值问题。
初值问题理论是微分方程研究的基础之一，本文将介绍微分方程中初值问题的理论基础和解法。
一、初值问题的定义初值问题是指给定一个微分方程及其解空间上一点的值，通过求解微分方程，确定解空间上满足给定初值条件的特定解。
初值问题的一般形式可以表示为：̇= (, )= ₀= ₀其中，表示未知函数，是自变量，是因变量，表示关于和的函数关系。
是关于和的函数，是任意给定实数。
初值问题的目标是找到满足上述方程和初值条件的特定解。
二、初值问题的解法解决初值问题的方法有很多种，常见的有解析解法和数值解法。
1. 解析解法解析解法是通过一系列数学手段，直接求得微分方程的解的公式，从而得到满足初值条件的特定解。
这种方法适用于某些特定形式的微分方程，例如线性微分方程、可分离变量的微分方程等。
解析解法的优势在于可以得到精确的解析表达式，从而能够准确描述问题的性质和变化规律。
但是，对于一些复杂的非线性微分方程，往往无法找到解析解，这时需要采用数值解法。
2. 数值解法数值解法是通过近似计算，利用离散的数值方法求解微分方程并得到数值近似解。
这种方法的思路是将微分方程转化为差分方程，并利用离散的计算方法逼近微分方程的解。
数值解法的优势在于适用性广，能够处理各种类型的微分方程，并能够得到任意精度的解。
常见的数值解法包括欧拉法、龙格-库塔法、改进欧拉法等。
三、初值问题的存在唯一性定理对于一阶常微分方程，初值问题存在唯一性定理是指在一定条件下，初值问题的解是存在且唯一的。
存在性定理：设 (, ) 是微分方程 , µ区间上的解且在 µ上连续，则初值问题在 [a,b] 上存在解。
唯一性定理：设 (, ) 和 (, ) 是微分方程在一定区域上的两个解，如果对于 µ [a,b] 上的某个点 x₀， ̇ (x₀) = ̇ (x₀)，那么在整个区域上µ， (x) = (x)，这就是说，在初值问题存在的条件下，初值问题的解是唯一的。
常微分方程初值问题的解法及应用
常微分方程初值问题的解法及应用常微分方程是数学中非常重要的一部分，它涉及了许多领域的模型建立和问题求解。
本文将介绍常微分方程初值问题的解法及其应用。
一、常微分方程初值问题的定义常微分方程初值问题是指给定一个常微分方程，以及它在某一点上的初始条件，求解该方程的解曲线。
通常，一个常微分方程初值问题可以表示为:y'(x) = f(x,y), y(x0) = y0,其中，y(x)是未知函数，f(x,y)是已知函数，y(x0) = y0是初始条件。
二、常微分方程初值问题的解法常微分方程初值问题的解法有多种，下面我们将介绍几种常用的方法。
1.欧拉法欧拉法是最简单的一种求解常微分方程初值问题的方法。
该方法基于初始条件，通过不断迭代计算得到近似解曲线。
具体步骤如下：步骤1：设定步长h，确定求解区间[x0, xn]，计算步数n。
步骤2：初始化，即确定初始点(x0, y0)。
步骤3：根据方程dy/dx = f(x,y)和初始点(x0, y0)，计算斜率k = f(x0, y0)。
步骤4：根据已知的斜率和步长h，计算下一个点的坐标(xi+1,yi+1)。
步骤5：重复步骤3和步骤4，直到达到步数n。
步骤6：得到近似解曲线。
2.改进的欧拉法（改进欧拉法）改进的欧拉法是对欧拉法的改进，其求解精度比欧拉法更高。
具体步骤如下：步骤1：设定步长h，确定求解区间[x0, xn]，计算步数n。
步骤2：初始化，即确定初始点(x0, y0)。
步骤3：根据方程dy/dx = f(x,y)和初始点(x0, y0)，计算斜率k1 =f(x0, y0)。
步骤4：根据已知的斜率k1和步长h/2，计算中间点的坐标(x0+h/2, y0+k1*h/2)。
步骤5：根据方程dy/dx = f(x,y)和中间点的坐标(x0+h/2, y0+k1*h/2)，计算斜率k2= f(x0+h/2, y0+k1*h/2)。
步骤6：根据已知的斜率k2和步长h，计算下一个点的坐标(xi+1,yi+1)。
常微分方程初值问题数值解法
数值解法的必要性
实际应用需求
许多实际问题需要求解常微分方程初值问题，如物理、 化学、生物、工程等领域。
解析解的局限性
对于复杂问题，解析解难以求得或不存在，因此需要 采用数值方法近似求解。
数值解法的优势
未来发展的方向与挑战
高精度算法
研究和发展更高精度的算法，以提高数值解的准确性和稳定性。
并行计算
利用并行计算技术，提高计算效率，处理大规模问题。
自适应方法
研究自适应算法，根据问题特性自动调整计算精度和步长。
计算机技术的发展对数值解法的影响
1 2
硬件升级
计算机硬件的升级为数值解法提供了更强大的计 算能力。
它首先使用预估方法（如欧拉方法）得到一个 初步解，然后使用校正方法（如龙格-库塔方法） 对初步解进行修正，以提高精度。
预估校正方法的优点是精度较高，且计算量相 对较小，适用于各种复杂问题。
步长与误差控制
01
在离散化过程中，步长是一个重要的参数，它决定 了离散化的精度和计算量。
02
误差控制是数值逼近的一个重要环节，它通过设定 误差阈值来控制计算的精度和稳定性。
能够给出近似解的近似值，方便快捷，适用范围广。
数值解法的历史与发展
早期发展
早在17世纪，科学家就开始尝 试用数值方法求解常微分方程。
重要进展
随着计算机技术的发展，数值 解法在20世纪取得了重要进展， 如欧拉法、龙格-库塔法等。
当前研究热点
目前，常微分方程初值问题的 数值解法仍有许多研究热点和 挑战，如高精度算法、并行计
软件优化
软件技术的发展为数值解法提供了更多的优化手 段和工具。
普通微分方程的初值问题
普通微分方程的初值问题微分方程是描述自然现象和工程问题中物理量的变化规律的数学工具。
其中，最常见的是普通微分方程，也就是仅仅涉及一个自变量和一个因变量的微分方程。
普通微分方程的初值问题，是指给定微分方程的初始值或初值条件，求解能够满足该条件的解函数的过程。
我们以一般的一阶线性微分方程为例，尝试深入了解普通微分方程的初值问题。
一、普通微分方程的求解一阶线性微分方程通常采用分离变量法求解，具体步骤如下：1. 将方程变为 dy/dx = f(x)g(y) 的形式；2. 将变量 y 和 x 分别表示为一个积分变量和一个常数变量；3. 将 dx 和 dy 带入方程后进行积分，得到解的通解函数 y = F(x,C)（C为常数）；4. 根据初始条件确定常数 C 的值，得到特解 y = f(x)。
二、一阶线性微分方程的初值问题一阶线性微分方程的一般形式为：dy/dx + p(x)y = q(x)其中， p(x) 和 q(x) 都是已知的函数。
对这个微分方程，我们设它的解为 y(x)，它的初值条件为 y(x0) = y0。
那么，这个方程的解就是 y(x)。
其实，我们求解一阶线性微分方程的初值问题的过程，就是要找到满足这个方程的 y(x) 和初始条件 y(x0) = y0 的函数。
步骤如下：1. 根据一阶线性微分方程，将已知的 p(x) 和 q(x) 带入通解公式y(x) = e^-P(x) * ∫（e^P(x)*q(x)dx + C）其中，P(x) = ∫p(x)dx2. 利用初始条件 y(x0) = y0，确定常数 C 的值，得到特解 y(x) = f(x)。
在求解一阶线性微分方程的初值问题时，需要注意以下几个方面：1. 判断通解是否存在如果微分方程解不唯一，可能会出现通解不存在的情况。
例如，当 p(x) = 0 时，方程就变为 dy/dx = q(x)，它的通解是 y(x) = C +∫q(x)dx，前提是该积分存在。
解常微分方程初值问题
解常微分方程初值问题常微分方程初值问题是求解一个确定初始值条件下的常微分方程的解。
解常微分方程的方法有很多种，下面将介绍几种常用的方法和相关参考内容。
1. 变量分离法：将微分方程中的变量分离，然后进行分离变量的积分。
这是解常微分方程最常用的方法之一。
相关参考内容：《普通微分方程教程》（陈英席著）、《普通微分方程》（王永乐著）2. 齐次方程法：对于齐次方程 dy/dx = f(x,y)（其中 f(x,y) 是关于 x 和 y 的函数），通过引入新的变量 u = y/x，将其转化为一个关于 u 的单变量方程。
然后再解这个方程。
相关参考内容：《普通微分方程与应用》（杨万明、杨卓玲著）、《数学物理方程》（尤伯杯著）3. 线性方程法：对于形如 dy/dx + P(x)y = Q(x) 的线性方程，可以使用积分因子法将其转化为一个可解的方程。
相关参考内容：《普通微分方程讲义》（陈方正、李学勤著）、《分析数学基础讲义》（包维楷等著）4. 变换法：通过进行适当的变量变换，将原方程转化为易于求解的形式。
相关参考内容：《常微分方程讲义》（李鼎立著）、《常微分方程教程》（张世忠、赵寿明著）5. 解特殊的微分方程：一些特殊的微分方程有相应的解法，例如 Bernoulli 方程、Riccati 方程等。
相关参考内容：《常微分方程教程》（孙士焜著）、《微分方程教程》（刘川著）此外，常微分方程的初值问题可以利用数值方法进行求解，例如 Euler 方法、Runge-Kutta 方法等。
相关参考内容：《数值分析》（李庆扬、褚国新著）、《常微分方程数值解法》（赵义、余长星著）解常微分方程初值问题需要动用到微积分、线性代数等数学知识，因此具备扎实的数学基础是解题的前提。
上述参考内容对于理解和掌握常微分方程的解法都具有很好的帮助，读者可以根据自己的实际情况选择适合的参考教材进行学习。
此外，还可以通过参考数学相关的学术论文和网络资源来进一步深入了解常微分方程的解法。
一阶常微分方程的数值求解
作业 利用Euler方法和R-K方法求解一个 常微分初值问题，并比较数值结 果，计算数值解和解析解的误差。
在 [ xk , xk 1 ] 内多取几个点，将它们的导数加权平均代 替 f ( x, y( x)) ，设法构造出精度更高的计算公式。
常用的是经典的 四阶R-K方法
y0 y( x0 ), xk 1 xk h yk 1 yk h (L1 2 L2 2 L3 L4 )/6
若 f 在 D {a x b,| y | } 内连续，且满足 Lip 条件：
dy f ( x , y) , y( x0 ) y0 , x [a, b] dx
L 0, s.t.| f ( x, y1 ) f ( x, y2 ) || y1 y2 | , 则上述问题的连续可
Matlab函数数值求解
[T,Y] = solver(odefun,tspan,y0)
其中 y0 为初值条件，tspan为求解区间；Matlab在数值求解 时自动对求解区间进行分割，T (向量) 中返回的是分割点的 值(自变量)，Y (向量) 中返回的是解函数在这些分割点上的函 数值。
solver 为Matlab的ODE求解器（可以是 ode45、ode23、
其中
L1 L2 L3 L4
f ( xk , yk ) f ( xk h / 2, yk hL1 / 2) f ( xk h / 2, yk hL2 / 2) f ( xk h, yk hL3 )
例3：利用四阶R-K方法求解例1与例2，并与Euler方法的 数值解进行比较。
y( xk 1 ) y( xk ) y( xk 1 ) y( xk ) dy O ( h) dx x h h k
第八章_常微分方程初值问题的单步法
为使局部截断误差为O(h ) ,应取
则
k j f ( xi c j h, y( xi c j h)) f ( xi c j h, y( xi ) hc j y( xi )) f ( xi c j h, y( xi ) hc j f ( xi , y( xi )))
f ( xi c j h, yi h a jm km )
f xx f xx ( xi , yi ), f yy f yy ( xi , yi ), f xy f xy ( xi , yi ),
由此得
k1 f , k2 f ( xi c2 h, yi c2 hk1 ) f ( xi , yi ) c2 h( f x k1 f y ) O(h ),
一. Euler方法
a x0 x1 x2 xN 1 xN b, ba x j x0 jh, h , j 1, 2, , N . N
y0 y ( x0 ), yi 1 yi hf ( xi , yi ), i 0,1, , N 1
一阶常微分方程初值问题的数值方法
------单步法
武汉大学数学与统计学院
一阶常微分方程初值问题的一般形式是:
y f ( x, y ), a x b (1) y(a) y0 D {( x, y ) a x b, c y d }
称f(x,y)在区域D上对y满足Lipschitz 条件是指:
定理显然得证.
8.1.2 一阶常微分方程初值问题的 Runge-Kutta方法
考虑一阶常微分方程初值问题
y f ( x, y ), a x b, y (a ) y0 ,
常微分方程初值问题数值解法
0.4 1.3582 1.3416 0.9 1.7178 1.6733
0.5 1.4351 1.4142 1.0 1.7848 1.7321
7
初值问题（2.2）有解 y ，1按2这x 个解析式子
算出的准确值 y(x同n )近似值 一y起n 列在表9-1中，两者 相比较可以看出欧拉方法的精度很差.
17
所以，局部截断误差可理解为用方法（2.10）计算一步的 误差，也即公式（2.10）中用准确解y(x代) 替数值解产生
的公式误差.
根据定义，显然欧拉法的局部截断误差
Tn1 y( xn1) y( xn ) hf ( xn , y( xn ))
y(xn h) y(xn ) hy(xn )
y(2) n1
yn
hf
( xn1,
y (1) n1
).
11
如此反复进行，得
y (k 1) n1
yn
hf
( xn1,
y(k) n1
),
(k 0,1, ).
（2.6）
由于 f (x,对y) 满足y 利普希茨条件（1.3）. 由（2.6）减 （2.5）得
y (k 1) n 1
yn1
h
f
( xn1,
y(k) n 1
积分曲线上一点 (x的, y切)线斜率等于函数 值.
的f (x, y)
如果按函数 f (在x, y) 平x面y上建立一个方向场，那 么，积分曲线上每一点的切线方向均与方向场在该点的方 向相一致.
基于上述几何解释，从初始点 P0 (x出0 ,发y0，) 先依 方向场在该点的方向推进到 x 上x1一点 ，P然1 后再从 P1 依方向场的方向推进到 x 上x2一点 ，循P2此前进做出
《常微分方程指导与实验》第2章：一阶微分方程的解的存在定理
第二章 一阶微分方程的解的存在定理§2.1 一阶微分方程解的基本理论主要内容一 导数已解出方程初值问题解的存在唯一性定理 考虑导数已解出的一阶DE 的初值问题()()00,y f x y y x y '=⎧⎪⎨=⎪⎩(2.1)(2.2)这里()y x f ,是在闭矩形域R ： a x x ≤-0，b y y ≤-0上的连续函数。
定义2.1 如果存在常数0>L ，使得对于所有的点()1,y x ，()2,y x R ∈，都有不等式()()2121,,y y L y x f y x f -≤-成立，则称函数()y x f ,在R 上关于y 满足李普希兹（Lipschitz ）条件。
1定理2.1 （毕卡存在唯一性定理） 如果()y x f ,在R 上满足条件： 1）连续；2）关于y 满足李普希兹条件，则初值（2.1）和（2.2）在区间h x x ≤-0上存在唯一解()x y y =，其中()M b a h ,m in=，()y x f M R y x ,max ),(∈=。
注1 取数h 的意义。
注意到()y x f M R y x ,max ),(∈=，从而积分曲线()x y y =在任一点()()R x y x ∈,处的切线斜率()M x y ≤'。
于是从点()o y x p ,0引两条斜率分别为M 和M -的直线1l 和2l ，便知过点P 的积分曲线必限制在图2.1和图2.2的阴影区域内。
而直线1l 和2l 相交情形有如下两种可能。
（i ）若相交成如图 2.1所示的情况，则a Mb>，积分曲线()x y y =在a x x ≤-0上不越出R ,从而应取a h =。
（ii ）若相交成如图 2.2所示的情况，则a Mb >，积分曲线()x y y =在Mb x x ≤-0上不越出R ,从而应取Mb h =。
总之，取()M ba h ,min =，就是为了使初值问题（2.1）和（2.2）的解在h x x ≤-0上总存在。
数值分析(25) 常微分方程初值问题的
(9-1)
只要 f (x, y) 在[a, b] R1 上连续，且关于 y 满足 Lipschitz 条 件，即存在与 x, y 无关的常数 L 使 | f ( x, y1 ) f ( x, y2 ) | L | y1 y2 | 数值解 对任意定义在 [a, b] 上的 y1(x) 和 y2(x) 都成立，则上述IVP存 在唯一解。
3、 中点欧拉公式 /* midpoint formula */
中心差商近似导数
y( x 2 ) y( x0 ) y( x1 ) 2h
y( x2 ) y( x0 ) 2h f ( x1 , y( x1 ))
y n1 yn1 2hf (xn , y n ) n 1, 2,
h yn1 yn ( f ( xn , yn ) f ( xn 1 , yn hf ( xn , yn )) 2
可得
h yn 1 yn h (1 )( xn yn ) 1 2
yn 1 0.905 yn 0.095 xn 0.1
x0
x1
y1
x0 x1
过点 ( x0 , y0 ) ，以 f ( x0 , y0 ) 为切线斜率的
切线方程为
y y0 f ( x0 , y0 )( x x0 )
用 y1 y0 f ( x0 , y0 )( x1 x0 ) y0 hf ( x0 , y0 ) 近似代替 y ( x1 )
若数值公式的局部截断误差 En 1 O ( h p 1 ) ，则称此数值 公式是 p 阶的， p 是正整数。
欧拉法的局部截断误差,由Taylor展开: En 1 y( xn 1 ) yn 1
第6章常微分方程初值问题的解法
ykh 2 k[ (ykx k 1 ) ( yk 1x k 1 1 )]
yk11 29 1yk1k05110
预估－校正Euler方法：
y k 1 0 .90 y k 5 0 .00 k 9 0 .1 5
20
Euler方法
xk
yk
yk y(xk)
0.0 1.000000
0.0
梯形方法
yk
yk y(xk)
1.000000
0.0
续
预估－校正方法
yk
yk y(xk)
1.000000
0.0
0.1 1.000000 0.2 1.010000
4.8×10－3 8.7×10－3
1.004762 1.018594
y(0) 1
其解析解为： y1xe-t2dt x[0,1] 0 很难得到其解析解
4
例如：
y=x+y , x[0,1]
y(0) 1
其解析解为 yx12ex
只有一些特殊类型的微分方程问题能够得到用解析表达式 表示的函数解，而大量的微分方程问题很难得到其解析解。
因此，只能依赖于数值方法去获得微分方程的数值解。
例如：
y=x+y , x[0,1]
y(0) 1
其解析解为：yx12ex
3
但是, 只有一些特殊类型的微分方程问题能够得到用解析 表达式表示的函数解，而大量的微分方程问题很难得到其解 析解。
因此，只能依赖于数值方法去获得微分方程的数值解。
例如：
y =e-x2 ,
x[0,1]
7.5×10－5 1.4×10－4
数学物理第六章-特征线法
x at c1, x at c2
解二 考虑弦振动方程的Cauchy问题
utt u(
x,
a2uxx 0,
0) (x),
x ut (x, 0)
,t 0
(x),
x
(1) (2)
特征线族
dx dt
2
a2
0
(3)
(3)称为特征方程。
即
dx a 0, dx a 0
dt
dt
:x (
x(s),t dx , dt , ds ds
t(s),u u(x(s),t
du ds
)
|P0
。显然，向量
(s)) n与
在点 P0 的切向量为 在点 P0 相互垂直。
记 (a(x, t, u), b(x, t, u), c(x, t, u)) |P0 , 则方程(9)式 恰好表示法向量 n 与 在点 P0 处相互垂直。因此，在曲面 u u(x,t) 上，以向量 为切向量且与曲线 相交，交点为
(x t 1)et 1
所以
u(x,t) 1 (x t 1) et 1 (x t 1)e2t
2
2
下面考虑一阶拟线性方程，即一阶导数的系数与未知函数 u 有关。
一阶拟线性方程柯西问题的一般形式为
a( u(
x, x,
t, u)ux
0)
b(x,t,u)ut c( (x), x
( , 0,( )) 的曲线方程为
dx
ds
a( x, t, u),
x(0)
dt
ds
b( x, t, u),
t(0) 0
(11)
du ds
c(
x,
t,
u),
一阶常微分方程初值问题的数值解fortron
一阶常微分方程初值问题的数值解是微分方程数值求解中的基础问题，对于工程、物理、生物等领域的科学计算和数值模拟具有重要意义。
本文将从常微分方程初值问题的数值解的基本原理和数值方法入手，详细介绍使用Fortran语言进行一阶常微分方程初值问题的数值解的实现过程和注意事项。
1. 常微分方程初值问题的数值解基本原理常微分方程初值问题的数值解是通过数值方法来逼近微分方程的解。
对于一阶常微分方程初值问题：dy/dx = f(x, y)y(x0) = y0其中，f(x, y)是给定的函数，y(x0) = y0是给定的初值条件。
求解该初值问题即是要找到一个函数y(x)近似地满足该微分方程，并且在点x = x0处与给定的初值条件相符。
2. 常微分方程初值问题的数值解的数值方法常见的数值方法包括欧拉方法、改进的欧拉方法、四阶龙格-库塔方法等。
其中，四阶龙格-库塔方法是最常用和最经典的数值方法之一。
该方法通过取若干个函数值点上的斜率的加权平均值来逼近微分方程的解，具有较高的数值精度和稳定性。
3. 使用Fortran语言实现一阶常微分方程初值问题的数值解Fortran是一种古老但经典的科学计算语言，以其高效的数值计算和科学工程计算而闻名。
下面将结合Fortran语言的特点，介绍如何使用Fortran语言实现一阶常微分方程初值问题的数值解。
（1）定义常微分方程的函数f(x, y)在Fortran程序中，首先需要定义常微分方程的函数f(x, y)，并将其定义为一个子程序或函数。
这里以一个简单的一阶线性常微分方程为例：f(x, y) = x + y则在Fortran程序中可以这样定义：```function f(x, y)real :: x, y, ff = x + yend function f```（2）实现四阶龙格-库塔方法在Fortran程序中，可以实现四阶龙格-库塔方法来数值解常微分方程初值问题。
具体做法是按照龙格-库塔方法的算法，在程序中编写相应的代码实现。
9、常微分方程初值问题数值解法
( k +1) yn +1
xn +1 ∫xn
− yn +1 |≤
hL 2
|
(k ) yn +1
− yn +1 |,
( 只要 hL < 1，则( 2.8)的ynk +1)收敛到(2.7)的yn +1. +1 2
三、单步法的局部截断误差与阶
一阶常微分方程初值问题(1.1)(1.2)的单步法的一般形式 yn +1 = yn + hϕ ( xn , yn , yn +1, h).
clear x=0,yn=1 %初始化 for n=1:10 yp=yn+0.1*(yn-2*x/yn); %预测 x=x+0.1; yc=yn+0.1*(yp-2*x/yp) ; yn=(yp+yc)/2 %校正 end
( 2.2)
作业： 作业：P381, 1, 2(1).
龙格—库塔 库塔(Runge-Kutta)法 §3 龙格 库塔 法
进一步 y ( xn +1 ) = y ( xn ) + ∫
xn +1 xn
f ( x, y ( x))dx,
(Байду номын сангаас.3)
∫
⇒ 其中
xn +1 xn
f ( x, y ( x))dx ≈ h ∑ ci f ( xn + λi h, y ( xn + λi h)).
yn +1 = yn + hϕ ( xn , yn , h),
i =1
r
(3.4) (3.5) 欧拉法r = 1, p = 1.改进
一阶常微分方程的初值问题
printf( " %12.5f %15.6e \n", t_new, y );
do{
for( j = 1; j <= nstep_pr; j++ ){
t_old = t_new;
t_new = t_new + h;
yn = y;
t_mid = t_old + hh;
yn = y;
k1 = h*fun( yn, t_old );
scanf( "%d", &nstep_pr );
printf( "Maximum t?\n" );
scanf( "%f", &t_limit );
y = 1.0; /* Setting the initial value of the solution */
h = t_pr/nstep_pr;
printf( "h=%g \n", h );
用经典4阶 Runge-Kutta 法求解的格式为
h
ui1 ui 6 k1 2k2 2k3 k4
h
vi1 vi 6 c1 2c2 2c3 c4
k1 f ( xi , ui , vi )
c1 g( xi , ui , vi )
h1
1
k2 f ( xi 2 , ui 2 hk1 , vi 2 hc1 )
它每一步计算f(x,y)两次，截断误差为O(h3)
dy y dt
y(0) 1, 0 t 1
精确解： y e t
function [t,y] = Heun(ode,tspan,h,y0) t = (tspan(1):h:tspan(end))'; n = length(t); y = y0*ones(n,1); for i=2:n
常微分方程--第三章 一阶微分方程的解的存在定理(3.1-3.2)_OK
x
L x0 1( ) 0 ( )d
L
x x0
M (
x0 )d
ML 2
(x
x0 )2
其中第二个不等式是由Lipschitz条件得到的,
由Lipschitz条件
17
设对于正整数n, 有不等式
n (x) n1(x)
MLn1 n!
(x
x0
)n
,
则当x0 x x0 h时,由Lipschitz条件有
dy dx
f
(x, y), (3.1)
y(x0 ) y0
证明: 若y (x)为(3.1)的连续解,则
d ( x)
dx
f
( x, ( x)),
(x0 ) y0
对第一式从x0到x取定积分得
x
即
x (x) (x0 ) x0 f (x,(x))dx (x) y0 x0 f (x,(x))dx
x
f ( , ( )) f ( ,( )) d x0
x
x
L ( ) ( ) d L g( )d
x0
x0
令u(x) L
x
g( )d ,
x0
则u(x)是定义于[x0, x0 h]上连续可微函数,
且u(x0 ) 0,0 g(x) u(x), u'(x) Lg(x),于是
u(x) Lu(x), (u(x) Lu(x))eLx 0,
(4) (x)是积分方程(3.5)定义于[x0 h, x0 h]上连续解
且唯一.
9
下面分五个命题来证明定理,为此先给出
积分方程
如果一个数学关系式中含有定积分符号且在定积分符 号下含有未知函数, 则称这样的关系式为积分方程.
一类一阶常微分方程初值问题的无穷多解性
中 图 分 类 号 Ｏ１ ２ ｌ ７ －
一
阶常微 分方 程初值 问题
（）一 ， ＞ ０， （ ）＝ ０ ｔ ０
（１ ｆ∈ ｃ ０ 。 ）， Ｏ 一 ０ ５ ＞ ０ Ｖ５ ０ ｆ） Ｅ ，。 ， （） ， ） ， ＞ ． ２存 （ ） 则 问题 （ ） 在无 穷 多个 解 的充分 必要 条件 是 厂满足 １
充 分性 先求 解 正解 ． ３ ， 用 分离 变量法 ， 由（ ） 应 并 对 Ｓ ０到 ｔ 从 积分 ， 可得
和 正 解 （ ｔ ０时 ， （）＞ Ｏ 当 ＞ “￡ ）
（）一 ÷ ， ＞ ０ ｔ ，
而忽视 了这样 的无 穷 多个 非负 解 ： 对任 意 的 ｃ ０ ＞ ，
Ｆ， ０ ０≤ ｔ ｃ ≤ ，
的所 有无穷 多个 解． 定理 １ 设 ｂ － 和 厂分别 满足
再求 非负 解 ． 对任 意 的 ｃ ０ 由（ ）和 ＞ ， ３
（ ） ＝ ０， ｃ
（ １ ｂ ｃ ｏ ｏ ） （） ｘ＜ ｏ ， ｂ ） ∈ （ ，。 ，ｌ ｄ ６ 。
收稿 日期 ：０ ８— １ — ２ ｝ 改 日期 ：０ ０— ０ ２０ ０ １修 ２１ ５— ０ ． ４ 基 金 项 目 ： 家 自然 科 学 基 金 项 目 （０ ７ ０ ６ ， 东 省 自然 科 学 基 金 国 １０１６） 山
Ｕ （ ）＞ ０，ｔ ｌ￡ ＞ ｔ． １
０≤ ｔ ｃ ≤ ；
ｌ（一 Ｃ０ －Ｂｃ）４， ≥ ｃ Ｑ （１ Ｂ） （ ） ｔ ＞ ）ｒ
项 目（ ０ ９ Ｒ Ｏ ７ ５ ． ２０Ｚ Ｂ １９ ）
对方程 （） 用分离变量法 ， ２应 并对 从 ｆ ｔ 到 积分 ， 可得
Ｊ ｆ（ Ｊ… … ｏ ）一 一
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
y = 1.7178 x = 0.9000
y = 1.7848 x = 1.0000
§2 Euler方法
➢Euler方法的导出
f ( x) f ( x0 ) f '( x0 )( x x0 ) f ''( x0 ) / (2!)( x x0 )2 ...
将 y( xn1) 在点 xn处进行Taylor展开
称上述公式为向前Euler 公式。
若将 y( xn ) 在点 xn1 处进行Taylor展开
y( xn )
y( xn1)
hy( xn1)
h2 y(n )
2!
略去h2项：y( xn1) y( xn ) hf ( xn1, y( xn1))
然后用 yn代替 y( xn ) ，即得
可用拟合方法求该组数据 ( xi ,的yi近) 似曲线
欧拉(Euler)公式： (yn1-yn )/(xn1-xn )=f(xn,yn ), yn1 yn hf (xn ,yn ), (n 0,1, 2,L ).
(2.1)
例1 取h 0.1,利用Euler公式求解
y y 2x（, 0 x 1)
解：向前Euler公式： yn1 yn hf ( xn , yn )
yn1 0.1xn 0.9 yn 0.1
向后Euler公式： yn1 yn hf ( xn1, yn1)
yn1
1 1.1 (0.1xn1
yn
0.1)
y( x0 ) y0
常用的一些解析解法： 分离变量法、变量代换、 常数变易法、Lapalace变换等
Th9.1.1（解的存在唯一性）
对于初值问题 (，) 如果 f (在x,下y)列区域内连续：
G {a x b;| y | }
且关于 y满足Lipschitz条件，即存在常数 L， 使0
y
( 2.2)
y(0) 1,
解：yn1
yn
h( yn
2xn ) yn
1.1yn
0.2xn . yn
计算及结果如下.
clc;clear; y=1, x=0, %初始化
y= 1
x= 0
y = 1.1000 x = 0.1000
for n=1:10
y = 1.1918 x = 0.2000
第九章 常微分方程数值解法
/*Numerical Method for Ordinary Differential Equations*/
本章主要介绍一阶常微分方程初值问题的数值解法。
§1 引言
➢初值问题及其数值解的概念
dy 一阶常微分方程初值问题： dx
f ( x, y);a
xb ()
yn1 yn hf ( xn1, yn1) n 0,1, 2,L , N 1
称上述公式为向后Euler 公式。
向后Euler 公式为隐式格式，需要利用迭代法求解
例2：分别利用向前和向后Euler
dy
x
y 1
方法求解初值问题的数值 dx
解（取步长为 h 0.1 ） y(0) 1
| f ( x, y1) f ( x, y2 ) | L | y1 y2 |;x, y G
则初值问题 (存) 在唯一解，且解是连续可微的。
Def 1所谓数值解是指：在解的存在区间上取一系列点
x0 x1 x2 ... xn ...
逐个求出 y( xi )的近似值 yi (i 1, 2,3, ...)
y=1.1*y-0.2*x/y, x=x+0.1,
y = 1.2774 x = 0.3000
end
y = 1.3582 x = 0.4000
y = 1.4351 x = 0.5000
几何意义
y = 1.5090 x = 0.6000 y = 1.5803 x = 0.7000
y = 1.6498 x = 0.8000
y( xn1)
y( xn )
hy( xn )
h2 y(n )
2!
略去h2项： y( xn1) y( xn ) hf ( xn , y( xn ))
然后用 yn代替 y( xn ) ，即得
Hale Waihona Puke yn1 yn hf ( xn , yn ) n 0,1, 2,L , N 1
等距节点：xi x0 ih; h : 步长
➢初值问题 (的) 解析解及其数值解的几何意义：
y
(
x0
,
y
0
( )
x1, y1)
••
(
x2 ,
••
y
2)
•
( xn , y n )
•• y y( x)
积分曲线
o x0 x1 x2
xn
x
初值问题 (的) 解表示过点 ( x0的, y一0 )条曲线 初值问题 (的) 数值解表示一组离散点列 ( xi , yi )