常微分方程的解法 (二): 欧拉（Euler）方法

浅夏110

§2 欧拉（Euler）方法
2.1 向前 Euler 公式、向后 Euler 公式
Euler 方法就是用差分方程初值问题（3）的解来近似微分方程初值问题（1）的解， 即由公式（3）依次算出  的近似值   。这组公式求问题（1）的数值 解称为向前 Euler 公式。
' [6 D- x9 v3 ^4 h
' W7 r) n+ ~; s
/ k, O# ?5 ~, [" Y  v- |  i
# A- G3 _$ S: Q( b, U2.2 Euler 方法的误差估计
; U" P5 X' K( n& z  i  e0 h: {7 F对于向前 Euler 公式（3）我们看到，当n = 1,2,....时公式右端的  都是近似的， 所以用它计算的 会有累积误差，分析累积误差比较复杂，这里先讨论比较简单的 所谓局部截断误差。
2 o/ z# c$ i' ?& B8 ?4 {# y$ W: |
( d. `9 b: m  ~  V" Y
$ k$ i# q7 u( ~) |/ a- a' @+ l. ]
2 @5 [- C% n4 y* b2 d" F  X7 D6 B/ \

显然 p 越大，方法的精度越高。式（9）说明，向前 Euler 方法是一阶方法，因此 它的精度不高。

§3 改进的 Euler 方法
/ x3 q8 q& s/ E2 I. B3.1 梯形公式
- j, \1 q0 d+ @, n0 \& t, f利用数值积分方法将微分方程离散化时，若用梯形公式计算式(4)中之右端积分， 即
" s' i% \# U- q# T
) [0 j0 u0 L( e; s1 j
- C  F- ^& i/ c7 S
这就是求解初值问题（1）的梯形公式。

' v' _, b$ z- c, \9 t直观上容易看出，用梯形公式计算数值积分要比矩形公式好。梯形公式为二阶方法。 梯形公式也是隐式格式，一般需用迭代法求解，迭代公式为



4 ^2 _: A# p: l7 F如果实际计算时精度要求不太高，用公式（10）求解时，每步可以只迭代一次，由此导 出一种新的方法—改进 Euler 法。

3.2 改进 Euler 法
按式（5）计算问题（1）的数值解时，如果每步只迭代一次，相当于将 Euler 公式 与梯形公式结合使用：先用 Euler 公式求  的一个初步近似值  ，称为预测值，然 后用梯形公式校正求得近似值  ，即


* k: q( E0 E7 A: k/ S
4 s$ Z" k1 L, m' X5 G  K式（11）称为由 Euler 公式和梯形公式得到的预测—校正系统，也叫改进 Euler 法。

为便于编制程序上机，式（11）常改写成

; }7 h& f0 ?$ M  d3 }2 N

改进 Euler 法是二阶方法。
! L8 @4 y$ I* r4 r3 o2 U3 T

————————————————
版权声明：本文为CSDN博主「wamg潇潇」的原创文章，遵循CC 4.0 BY-SA版权协议，转载请附上原文出处链接及本声明。
2 E' n& G0 Y/ O2 @原文链接：https://blog.csdn.net/qq_29831163/article/details/89703276
% D  p; a# h# a. _

$ j& U& `! |4 ^