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

第7章 常微分方程的数值解法 §1 引言 §2 Euler法和改进的Euler法 §3 Runge-Kutta法 §4 程序示例 本章目的要求: 熟练掌握Euler公式，Runge-Kutta公式以及会用以上公式求常微分方程初值问题的数值解与程序设计. §1 引言 在工程和科学计算中，所建立的各种常微分方程的初值问题，除很少几类的特殊方程能给出解析解，绝大多数的方程是很难甚至不可能给出解析解的，其主要原因在于积分工具的局限性。因此，人们转向用数值方法去解常微分方程，并获得相当大的成功，讨论和研究常微分方程的数值解法是有重要意义的。 §2 Euler法和改进的Euler法 2.1 Euler法(数值微分(折线)法) 若将函数y(x)在点xi处的导数y′(xi)用差商来表示，即 再用yi近似地代替y(xi)，则初值问题(1) 就化为 2.1 Euler法(数值积分法) 设将方程 的两端从 到 求积分，即得 显然，只要能近似的算出其中的积分项，我们就可以得到计算 的公式。 若我们用梯形法计算积分项： 即可得如下计算公式 例1 用Euler法求初值问题 Euler法求解例1的Matlab程序: Euler法的运行结果: 2.2 改进的Euler法 Euler法虽然形式简单，计算方便，但比较粗糙，精度也低。为了达到较高精度的计算公式，对Euler法进行改进. 改进的Euler公式 这样建立的预报－校正系统称为改进的Euler公式： 预报 校正 为方便编程,写成形式： 实践表明，改进的Euler公式明显改善了精度. 改进的Euler公式进行计算: 改进的Euler法求解的Matlab程序: 改进的Euler法运行结果: Runge-Kutta法的设计思想 考察差商 ，根据微分中值定理，存在点 ，利 用所给方程 得 我们称 为区间 上的平均斜率，这样只要对平 均斜率 提供一种算法，相应地我们便导出一种计算格式。 Runge-Kutta方法设计思想就是设法在 内多预报几个点的斜率值，然后把它们加权平均作为平均斜率，以期望构造出更高精度的计算格式。 二阶Runge-Kutta方法 随意考察区间 内一点 ,用两个 点 的斜率 的加权平均代替平均斜率 ， 于是我们就得 到如下计算格式： 其中有两个待定参数 ， 适当选取它们的值，就可使上述格式有 较高的精度。若 ，该格式是二阶的 ，故统称满足这一条件 的一族格式为二阶Runge-Kutta格式。特别地，当 时，上 述格式即为改进的Euler格式，如果取 ，则上述格式称为 变形的Euler格式，亦称为中点格式。 三阶Runge-Kutta方法 为了进一步提高精度，我们可以考虑用三个点 的斜率值 加权平均得出平均斜率 的近似值，其中， 于是就可以构造所谓的三阶Runge-Kutta格式 ，下列Kutta格式是其中 的一种： 四阶Runge-Kutta方法 继续上述过程，我们可以导出四阶Runge-Kutta格式，下列经典 格式是其中的一种： 值得注意的是，Runge-Kutta法的推导基于泰勒展开法，因而它 要求解具有较好的光滑性。如果解的光滑性差，则该方法得到的解 反而不好。 四阶Runge-Kutta方法的Matlab程序: 四阶Runge-Kutta方法运行结果: 内容大纲 1.常微分方程的数值解(调用内置函数); 2.常微分方程的符号解(调用内置函数). 蓝色曲线 ——y（1）； （原方程解） 红色曲线 ——y（2）； 微分方程数值解 解高阶微分方程三步聚: 1.将高阶方程转化为一阶方程组; 2.