第八章 常微方程数值解.docVIP

第八章 常微分方程数值解 由常微分方程理论可知，我们只能求一些特殊类型的常微分方程。而实际上许多常微分方程求解非常困难。 本章主要讨论一阶常微分方程的初值问题： （8－1） 从理论上讲只要方程中的连续且关于满足李普希兹（Lipschitz）条件，即存在常数L，使 则常微分方程存在唯一解。 微分方程数值解：就是求微分方程的解在一系列离散节点 处的近似值（i=1，2，…，n） 称为由到的步长，通常取为常数h。 求数值解，首先将微分方程离散化，常用方法有： 用差商代替微商 若用向前差商代替微商，即 （i=1，2，…，n） 则得 即 数值积分法 利用数值积分法左矩形公式 = 可得同样算法 用泰勒（Taylor）公式 得离散化计算公式 §1 欧拉（Euler）方法 1.1欧拉方法 对一阶微分方程（8—1），等分区间为份，， 则 由以上讨论可知，无论用一阶向前差商，还是用数值积分法左矩形公式，或者用泰勒公式取前两项都可得到同样的离散化计算公式 代入初值则得到数值算法： （i=1，2，…，n－1） （8－2） 称其为欧拉方法。 几何上欧拉方法就是用一条折线近似表示曲线。（如图）： P0