第8章 常微分方程数值解法
本章主要内容：
1．欧拉法、改进欧拉法. 2．龙格-库塔法。
3．单步法的收敛性与稳定性。
重点、难点
一、微分方程的数值解法
在工程技术或自然科学中，我们会遇到的许多微分方程的问题，而我们只能对其中具有较简单形式的微分方程才能够求出它们的精确解。对于大量的微分方程问题我们需要考虑求它们的满足一定精度要求的近似解的方法，称为微分方程的数值解法。本章我们主要
讨论常微分方程初值问题⎪⎩⎪⎨⎧==00
)()
,(y
x y y x f dx dy
的数值解法。
数值解法的基本思想是：在常微分方程初值问题解的存在区间[a,b]内，取n+1个节点a=x 0＜x 1＜…＜x N =b （其中差h n = x n –x n-1称为步长，一般取h 为常数，即等步长），在这些节点上把常微分方程的初值问题离散化为差分方程的相应问题，再求出这些点的上的差分方程值作为相应的微分方程的近似值（满足精度要求）。
二、欧拉法与改进欧拉法
欧拉法与改进欧拉法是用数值积分方法对微分方程进行离散化的一种方法。
将常微分方程),(y x f y ='变为()
*+=⎰++1
1))(,()()(n x
n x n n dt
t y t f x y x y
1．欧拉法（欧拉折线法）
欧拉法是求解常微分方程初值问题的一种最简单的数值解法。
欧拉法的基本思想：用左矩阵公式计算（＊）式右端积分，则得欧拉法的计算公式为：N
a
b h N n y x hf y y n n n n -=
-=+=+)1,...,1,0(),(1 欧拉法局部截断误差
11121
)(2
++++≤≤''=n n n n n x x y h R ξξ或简记为O （h 2）。
我们在计算时应注意欧拉法是一阶方法，计算误差较大。
欧拉法的几何意义：过点A 0（x 0，y 0），A 1（x 1，y 1），…，A n （x n ，y n ），斜率分别为f （x 0，y 0），f （x 1，y 1），…，f （x n ，y n ）所连接的一条折线，所以欧拉法亦称为欧拉折线法。
例1用欧拉法解初值问题
⎪⎩
⎪⎨
⎧=≤≤-=1)0()
10(2y x xy dx dy
在x ＝0 (0.2) 1处的近似解。（计算过程保留4位小数）。 【思路】 用欧拉法求解常微分方程的初值问题时，首先熟练掌握欧拉公式的一般形式， 根据具体题目写出找出欧拉公式的迭代式，并根据初始条件和所给步长进行迭代求解。 解 ∵ f (x ，y )＝－2xy ，h ＝0.2，
欧拉公式为：
)
5,4,3,2,1,0()4.01()2(2.0),(1=-=-+=+=+n y x y x y y x hf y y n n n n n n n n n
列表计算如下：
2．改进欧拉法
改进欧拉法比欧拉法的计算准确，是对欧拉法的改进。改进欧拉法的基本思想：用梯形公式计算（＊）式右端积分，则得改进欧拉法的计算公式为：
[]N
a b h N n y x f y x f h
y y n n n n n n -=
-=++
=+++)1,...,1,0(),(),(2
111
利用改进欧拉法计算常微分方程初值问题时，我们应注意此公式为隐式表达式，需要对它进行迭代求解。计算时可以采用一次迭代和多次迭代，因此，就有改进欧拉法预估-校正法公式和反复迭代的改进欧拉法预估-校正法公式。 改进欧拉法预估-校正法公式：
()()[]
⎪⎩
⎪⎨
⎧-=++=+=++++)
1,...,1,0(),(),(2),(011101N n y x f y x f h y y y x hf y y n
n n n n n n n n n
反复迭代的改进欧拉法预估-校正法公式：
()()()[]
⎪⎩⎪⎨
⎧=-=++=+=+++++),...,1,0,1,,0(),(),(2)
,(1111
01m N n y x f y x f h y y y x hf y y m n n n n n m n n n n n
改进欧拉法的局部截断误差1113
1
)(12
++++≤≤'''=n n n n n x x y h R ξξ或简记为O （h 3）。从局部截断误差的形式看，改进欧拉法是二阶方法，因此，它比欧拉法更精确。 例2用预估-校正法求初值问题
⎪⎩⎪⎨⎧=≤≤--='1
)0()
10(2
y x xy y y
在x=0（0.2）1的解。
【思路】掌握预估-校正法的计算公式，根据已知条件迭代求解。 解 步长h=0.2，将2
),(xy
y y x f --=代入预估-校正公式，整理得
⎪⎩⎪⎨
⎧+--=-=+++++2
)0(11)0(1212)0(1
)
((1.01.09.02.08.0n n n n n n n n n n n y x y y x y y y x y y 列表计算如下：
例3用改进欧拉法求解例1的初值问题,要求3)
1()(10--<-m n
m n y y 。
【思路】掌握改进欧拉法的计算公式，根据已知条件迭代求解，并检验迭代解是否满
足精度要求，若满足则确定此解为常微分方程在某点的近似解。 解 将xy y x f 2),(-=代入改进欧拉法的计算公式得：
()()[]()
[]
⎪⎩⎪⎨
⎧+-=++=-=-+=+=-++-++++)
2.0),(),(2)4.01()2(2.0),(1111111
)0(1m n n n n n m n n n n n m n n n n n n n n n n
y x y x y y x f y x f h y y y x y x y y x hf y y
列表计算如下：
