2012数值分析 -lec23--常微分方程数值解法-Euler.pptVIP

数值分析 * 朱立永 北京航空航天大学 数学与系统科学学院 Email: numerical_analysis@ Password:beihang2012 答疑时间：星期四下午2:30－5:30 答疑地点：主216 * 第二十二讲欧拉方法 第七章 常微分方程初值问题的数值解法 * 常微分方程及其求解概述 初值问题解法 1.Euler方法 2.线性多步法 3.Runge--Kutta法 4.方程组及刚性问题的Gear方法 * 常微分方程及求解概述 基本概念 描述自由落体的ODE: * ? 只有一个自变量的微分方程为 ODE, 否则称为偏微分方程 PDE。 ? 方程中未知函数导数的最高阶数称为方程的阶。 ? 方程中关于未知函数及其各阶导数均是一次的，则称为线性微 分方程。 是线性二阶ODE。 ? 初始条件亦称定解条件、初值问题 * ? 在没有给定解条件时，方程一般有一族解曲线 y(x,c) ，如： ? 对任意的n阶ODE，如果能写成： 则称该方程为显式的。前面提到的方程是显式的。 而下面方程是隐式的。 * ? 对于高阶显式方程。通过定义n-1个新变量，可以写成 n维一阶方程组。即令： * ? 在讨论初值问题时，我们从一阶方程开始： ? 当 f(x,y)与y无关时，f(x,y)=g(x) 然后毫不费力地套用来解方程组。 * 数值解及其重要性 ? ? * ODE数值解的基本思想和方法特点 基本思想有两点 1. 离散化 用Taylor级数，数值积分和差商逼近导数等手 段，把ODE转化为离散的代数方程(称差分方程)。 2. 递推化 在具有唯一解的条件下，通过步进法逐步计 算出解在一系列离散点上的值。从而得到原ODE 的数值近似解。 * 初值问题解法 我们讨论一阶ODE，而高阶可能化为一阶 ODEs。一阶初值问题可以一般地写成： 欧拉（Euler）方法 Euler方法是求解一阶初值问题最简单方法， 但精度差，故不实用。然而对理论分析很有用。 * 欧拉方法原理及推导 设一阶初值问题满足： * 欧拉方法的几何意义： 用分段的折线逼近函数,此为“折线法”而非“切线法”，除第一个点是曲线上的切线，其它都不是。 * * * 多 步 法 单步法 自动起步 显式 隐式 半隐式 图 3 ODE求解方法的类型 * Euler公式的误差 Euler方法的误差估计 一般其它方法的误差估计也类似。 这里误差是指截断误差(算法理论误 差) 而不是舍入误差。后者由计算机字长等决定，属于稳 定性问题。 i) 几何分析 * * 整体截断误差：设ym是Euler公式精确解，而 y(x) 是 初值问题的解, 则整体截断误差定义为 它是局部截断误差的积累。 定理：若f(x,y)关于y满足Lipschitz条件, 则有估计式： * * * * 注意： 稳定性： * * * * * Thanks for your attention! Any questions?!