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

第七章常微分方程的数值解法 计算机学院 陈克建 6 学时 本章内容 §7.1 引言 §7.2 欧拉方法 §7.3 龙格—库塔方法 §7.4 阿达姆斯方法 小结 作业与实验 本章要求 1. 熟悉Euler显格式,梯形法及Euler预校法; 2. 熟悉局部截断误差及绝对稳定性; 3. 掌握龙格—库塔法。 §7.1 引言 本节内容 一. 问题提出 二. 两类求定解问题 三. 数值解法含义 返回章节目录 §7.1 引言 一. 问题提出 有一个或多个导数及其函数的方程式称为微分方程，在工程中常遇到求解微分方程的问题。 §7.1 引言 §7.1 引言 二. 两类求定解问题 实际中求解常微分方程的所谓定解问题有两类:初值问题和边值问题 定解指已知因变量和/或其导数在某些点上是已知的(约束条件)。 1. 边解问题 约束条件为已知，在自变量的任一非初值上，已知函数值和/或其导数值，如 常常可以将边解问题转化为初值问题求解。 §7.1 引言 2. 初值问题 §7.1 引言 §7.1 引言 §7.1 引言 §7.1 引言 §7.1 引言 三. 数值解法含义 所谓数值解法, 就是设法将常微分方程离散化, 建立差分方程, 给出解在一些离散点上的近似值。 微分方程的数值解: 设方程问题的解y(x)的存在区间是[a,b], 令a= x0 x1… xn =b, 其中hk=xk+1-xk, 如是等距节点h=(b-a)/n, h称为步长。 由于y(x)的解析表达式不容易得到或根本无法得到, 我们用数值方法求得y(x)在每个节点xk上y(xk)的近似值, 用yk表示, 即yk≈y(xk), 这样y0, y1, ..., yn称为微分方程的数值解。 §7.2 欧拉方法 本节内容 一. 欧拉格式 二. 梯形格式 三. Euler预估—校正法 四. 误差估计、收敛性和稳定性 返回章节目录 §7.2 欧拉方法 一. 欧拉（Euler）格式 §7.2 欧拉方法 §7.2 欧拉方法 §7.2 欧拉方法 §7.2 欧拉方法 §7.2 欧拉方法 §7.2 欧拉方法 §7.2 欧拉方法 §7.2 欧拉方法 §7.2 欧拉方法 §7.2 欧拉方法 §7.2 欧拉方法 二. 梯形格式 §7.2 欧拉方法 §7.2 欧拉方法 §7.2 欧拉方法 §7.2 欧拉方法 §7.2 欧拉方法 §7.2 欧拉方法 §7.2 欧拉方法 §7.2 欧拉方法 §7.2 欧拉方法 §7.2 欧拉方法 §7.2 欧拉方法 §7.2 欧拉方法 §7.2 欧拉方法 §7.2 欧拉方法 §7.2 欧拉方法 §7.3 龙格 —库塔方法 本节内容 一. 引言 二. 2阶龙格 —库塔格式 三. 高阶龙格 —库塔格式 返回章节目录 §7.3 龙格 —库塔方法 一. 引言 §7.3 龙格 —库塔方法 §7.3 龙格 —库塔方法 §7.3 龙格 —库塔方法 §7.3 龙格 —库塔方法 二. 2阶龙格 —库塔格式 §7.3 龙格 —库塔方法 §7.3 龙格 —库塔方法 §7.3 龙格 —库塔方法 §7.3 龙格 —库塔方法 §7.3 龙格 —库塔方法 三. 高阶龙格 —库塔格式 §7.3 龙格 —库塔方法 §7.3 龙格 —库塔方法 §7.3 龙格 —库塔方法 §7.3 龙格 —库塔方法 §7.3 龙格 —库塔方法 §7.3 龙格 —库塔方法 §7.3 龙格 —库塔方法 §7.3 龙格 —库塔方法 §7.3 龙格 —库塔方法 §7.3 龙格 —库塔方法 §7.4 阿达姆斯方法 本节内容 一. 一般线性多步法 二. 亚当姆斯显式公式 三. 亚当姆斯隐式公式 四. 亚当姆斯预测—校正系统 返回章节目录 §7.4 阿达姆斯方法 一. 一般线性多步法 由于在计算yn+1时, 已经知道yn, yn-1, … , 及?(xn,yn), ?(xn-1,yn-1), … , 利用这些值构造出精度高、计算量小的差分公式就是线性多步法。 即用若干节点处的 y 及 y’ 值的线性组合来近似y(xi+1)。 §7.4 阿达姆斯方法 §7.4 阿达姆斯方法 二. 亚当姆斯显式公式 /* Adams explicit formulae */ §7.4 阿达姆斯方法 §7.4 阿达姆斯方法 §7.4 阿达姆斯方法 §7.4 阿达姆斯方法 三. 亚当姆斯隐式公式 /* Adams implicit formulae */ §7.4 阿达姆斯方法 §7.4 阿达姆斯方法 小结 §7.1 引言 §7.2 欧拉方法 §7.3 龙格—库塔方法 §7.4 阿达姆斯方法 作业与实验 作业（上作