《常微分方程初值问题的数值解法分析》1800字.docxVIP

PAGE16

常微分方程初值问题的数值解法分析

虽然求常微分方程解析解的方法多种多样，但这些方法只能对一些特殊类型的常微分方程进行求解.在实际问题中遇到的常微分方程，往往不能求出其解析解，所以在实际应用时需要通过数值计算求出满足要求的数值解.

在自然科学与社会科学的许多领域中，常常会遇到一阶常微分方程的初值问题

（1.1）

定理：如果方程（1.1）中的右端函数满足：

（1）是实值函数；

（2）函数在矩形区域内连续；

（3）关于满足Lipschitz条件：即存在正常数，使得对任意,下列不等式成立：

，

则方程（1.1）存在唯一的解.

对于方程（1.1），本文给出了三种常见的数值解法：欧拉法、改进的欧拉法、龙格-库塔法.

1.1欧拉法

欧拉法是最简单的方法，虽然其精度不高，但欧拉法的公式简单，也体现了其他方法的基本特征.

欧拉法的具体求解如下：将区间作等分，步长为，节点为：

.

由可以得到在的导数值为：

.

利用泰勒展式：

，

忽略二阶无穷小量，得到：，其中为的近似值.由可以算出，依次进行可计算出在所有节点上的近似值，可以得出递推公式：

（1.2）

欧拉法的几何意义：欧拉法实际上是用过已知点的折线（图(1)中的虚折线）来近似代替过该点的积分曲线.不难看出，越小，则该折线逼近积分曲线越好，因此人们也将欧拉法称为欧拉折线法.由泰勒展式导出的欧拉法是一种单步法，其局部