精选第8章 常微分方程初值问题的数值解法.pptVIP

1 第 8 章 常微分方程初值问 题的数值解法 ? 基础知识 ? 欧拉方法 ? 龙格 — 库塔方法 2 8.1 基础知识 一、问题的提出 很多实际问题都需要求解常微分方程。例如单摆问题。 常微分方程分为线性常微分方程和非线性常微分方程，又可以分为一阶 常微分方程和高阶常微分方程。通过变量的替换，可以把高阶常微分方 程转化为一阶常微分方程再求解。对于一阶常微分方程组，可以写成向 量形式的单个方程，求解方法与一阶常微分方程相似。因此本章只讨论 一阶常微分方程的初值问题： ? ? ? ? ? ? ? ? ? 0 0 ) ( , ) , ( y x y b x a y x f dx dy 目前在常微分方程理论中，只能求出某些特殊类型常微分方程的解析解， 对大部分常微分方程，用解析方法求出常微分方程的精确解非常困难， 甚至不存在解的解析表达式。为满足工程实践的需要，常常用数值解法 求常微分方程的近似解。 在本章中假设讨论的一阶常微分方程的初值问题的解 y(x) 存在、唯一且 足够光滑，方程本身是稳定的，即精确解 y(x) 连续且依赖于初始值及右 端函数。 3 8.1 基础知识 二、数值解法 一阶常微分方程初值问题的数值解法的主要思想，是对区间 [a,b] 上的节点 a=x 0 ＜ x 1 ＜ …… ＜ x n ＜ x n+1 ＜ ……≤b 建立 y(x n ) 的近似值 y n 的某一递推格式，利用初值 y 0 和已计算出的 y 1 ,y 2 ,……, y k-1 递推出 y k ，并且用这一方法反复递推，依次得到 y k+1 ,y k+2 ,……,y n 。这一 求解方法称为步进式求解，相邻 2 个节点的距离称为步长，记为 h i =x i+1 － x i 。 为便于计算，常取成等距节点，称为定步长，这时把步长记为 h 。 一阶常微分方程初值问题的数值解法有多种分类方法。 多步法不能自行启动，必须先用单步法计算出 y 1 ,y 2 ,……,y k-1 ，才能启动一 个 r 步的多步法。 一种分类方法为： ⑴ 单步法：每一轮递推只用到前面一轮的递推结果，递推格式为： y k =y k-1 ＋ hT(x k-1 ,y k-1 ) ⑵ 多步法：每一轮递推要用到前面多轮递推的结果，递推格式为： y k =y k-1 ＋ hT(x k-r ,y k-r ,x k-r+1 ,y k-r+1 ,……, x k-1 ,y k-1 ) ，其中 r ＞ 1 。 4 8.1 基础知识 二、数值解法（续） 另一种分类方法为： ⑴ 显式方法：递推公式的右端都是已知量，可以直接计算出递推的结 果，递推格式为： y k =y k-1 ＋ hT(x k-r ,y k-r ,x k-r+1 ,y k-r+1 ,……, x k-1 ,y k-1 ) ⑵ 隐式方法：递推公式左端的未知量也出现在公式的右端，递推格式 为： y k =y k-1 ＋ hT(x k-r ,y k-r ,x k-r+1 ,y k-r+1 ,……, x k ,y k ) 隐式方法的递推公式其实是一个方程。解方程的运算量可能较大，为避免 解方程，常采用预测 — 校正系统。 ⑶ 预测 — 校正系统：每一轮递推包括预测和校正这 2 个步骤。先用显式方 法计算出 y k ，作为迭代的初值，这一过程称为预测；再把隐式方法的递推 公式作为迭代公式，把预测值 y k 代入迭代公式右端进行迭代，这一过程称 为校正。在校正时往往迭代 1 次或几次，校正值的精度就会有大幅提高。 一阶常微分方程初值问题的数值解法一般是对连续的初值问题进行离散化 处理，把微分方程转化为代数方程来求解。常用的离散化方法有： ⑴ 基 于数值微分的离散化方法，⑵ 基于数值积分的离散化方法，⑶ 基于泰勒 展开的离散化方法。 5 8.2 欧拉方法 一、显式欧拉法 在一阶常微分方程初值问题的数值解法中，显式欧拉（ Euler ）法是最简 单的一种。显式欧拉法有明显的几何含义，缺点是精度不高。对于一阶 常微分方程的初值问题： ? ? ? ? ? ? ? ? ? 0 0 ) ( , ) , ( y x y b x a y x f dx dy 显式欧拉法的递推公式为： y k =y k-1 ＋ hf(x k-1 ,y k-1 ) ， k=1,2,3,…… 。 显式欧拉法每一轮递推只用到前面一轮递推的结果，因此它是单步法。 由基于数值微分的离散化方法、基于数值积分的离散化方法、基于泰 勒展