第8章微分方程数值解法.pptVIP

例1: 取步长 h = 2/10, 2/20, 2/30, 2/40, 分别用欧拉法、改进的欧拉法和梯形法求解. 解: 记 f (x, y) = y－ x y2, xk= k h (k = 0, 1, 2,···, n ) (1) . Euler法: yk+1 = yk + h( yk－ xk yk2 ) (k = 0, 1, ···,n) y0 = 1 当 h = 2/10时, n=10. 由Euler公式可得: k 0 1 2 3 4 yk+1 1.2 1.3824 1.506 1.53504 1.46503 k 5 6 7 8 9 yk+1 1.32877 1.17077 1.02113 0.89169 0.783788 (2) . 改进的Euler法: k 0 1 2 3 4 yk+1 1.1912 1.34384 1.42348 1.41905 1.3473 k 5 6 7 8 9 yk+1 1.23726 1.11424 0.994151 0.884751 0.788666 (3) . 梯形法(计算过程略) n 10 20 30 40 h 0.2 0.1 0.0667 0.05 误差 0.1059 0.0521 0.0342 0.0256 Euler法误差: 改进的Euler法误差: n 10 20 30 40 h 0.2 0.1 0.0667 0.05 误差 0.0123 0.0026 0.0011 5.9612e-004 预-校方法, h=0.2时 误差最大值: 0.0123 欧拉方法, h=0.2时 误差最大值: 0.1059 解析解: 第8 章 微分方程数值解法 基础教学部数学教研室 彭 晓 华 内容提纲 引言 Euler法及其改进 Runge-Kutta方法 线性多步法 误差分析 数值解法的收敛性、相容性和稳定性 边值问题数值解法简介 包含自变量、未知函数及未知函数的导数或微分的方程称为微分方程。 在微分方程中, 自变量的个数只有一个, 称为常微分方程。 自变量的个数为两个或两个以上的微分方程叫偏微分方程。 微分方程中出现的未知函数最高阶导数的阶数称为微分方程的阶数。 如果未知函数y及其各阶导数 8.1 引 言 8.1.1 微分方程知识回顾 都是一次的,则称它是线性的,否则称为非线性的。 在高等数学中，对于常微分方程的求解，给出了一些典型方程求解析解的基本方法： 一阶：可分离变量法、齐次方程、一阶线性方程以及伯努利方程 高阶：可降阶方程、 常系数齐次线性方程的解法、常系数非齐次线性方程的解法 但能求解的常微分方程仍然是有限的，大多数的常微分方程是求不出解析解的。 这个一阶微分方程就不能用初等函数及其积分来表 达它的解。 例如 从实际问题当中归纳出来的微分方程，通常主要依靠数值解法来解决。本章主要讨论一阶常微分方程初值问题 ( 8.1 ) 在区间 可以证明,如果函数在带形区域 R={a≤x≤b,-∞＜y＜∞}内连续，且关于y满足李普希兹(Lipschitz)条件，即存在常数L(它与x,y无关)使 对R内任意两个 都成立,则方程( 8.1 )的解 在?a, b?上存在且唯一。 上的数值解法。 是假定h为定数，称为定步长，这时节点可表示为 8.1.2 数值方法的基本思想 1、数值解法需要把连续性的问题加以离散化，从而求出离散节点处的数值解。 常微分方程初值问题(8.1)式的数值解法，首先要算出精确 解y(x)在区间?a,b?上的一系列离散节点 处的函数值 的近似值: 相邻两个节点的间距 称为步长，本章总 描述这类算法，要求给出用已知信息 计算 的递推公式。 建立这类