第7章 常微分方程数值解法.pptVIP

第7章 常微分方程数值解法 §1 引言 §2 欧拉法和改进的欧拉法 §3 龙格库塔法 §4 阿达姆斯方法 §5 二阶线性常微分方程边值问题的数值解 §1 引言 在常微分方程中，我们已经掌握了一些典型方程的解法。但许多形式的方程只能用数值方法求近似解，也就是求在某些点上满足一定精度的近似解。现以求一阶常微分方程初值问题 在区间［a，b］上的解为例，介绍数值方法的基本思想。 设f x，y 在带形区域 R：｛a≤x≤b，-∞＜y＜+∞｝ 上为x，y的连续函数，且对任意的y满足李普希茨 Libusize 条件 ｜f x，y1 -f x，y2 ｜≤L｜y1-y2｜ 6―2 其中 x，y1 、 x，y2 ∈R，L为正常数。在求初值问题 6―1 的数值解时，我们通常采用离散化方法 数值微分、数值积分、泰勒展式等 ，求在自变量x的离散点 a x0＜x1＜x2＜…＜xn b 上的准确解y x 的近似值 y0，y1，y2，…，yn 常取离散点x0，x1，x2，…，xn为等距，即 x i+1-xi h，i 0，1，2，…，n-1 h称为步长。图6.1表示为初值问题 6―1 在n+1个离散点上的准确解y x 的近似值。 §2 欧拉法和改进的欧拉法 2.1 欧拉法 折线法 若将函数y x 在点xi处的导数y′ xi 用差商来表示，即 再用yi近似地代替y xi ，则初值问题