数值分析21常微分方程数值解.pptVIP

*/27 */27 解一阶常微分方程欧拉法 局部截断误差与p阶精度 Range-Kutta公式 一阶常微分方程组和二阶方程 线性多步法简介 《数值分析》 23 ? ? ? ? ? * 常微分方程(ODE)的定解问题主要有初值问题(IVP)和边值问题(BVP)两大类。 急性传染病数学模型(SIR模型) 常微分方程的初值问题是描述系统发展演变的重要工具。 洛伦兹吸引子? lorenzgui 关于爱情的动力系统 ? Lover affaires and differential equations * 一阶常微分方程初值问题: 离散化 (discretization) 取定离散点 x0 x1 x2 ··· xN,步长h 其中y = y(x) 是未知函数, 右端函数 f(x, y )是已知函数, 初值 y0 是已知数据。 求未知函数 y(x) 在离散点处的近似值 y1, y2, y3, ·····, yN 考察点xn列出微分方程 Euler法与修正的Euler法 其中yn是y(xn)的近似。 足够多的点来近似连续对象 * 求解常微分方程初值问题的Euler方法 取定步长: h,记 xn = x0 + nh, ( n = 1,2, ···, N ) 称计算格式: yn+1 = yn + h f( xn, yn ) 为显式Euler公式。 对应的求初值问题数值解的方法称为Euler方法。 例1. 用Euler法求初值问题 的数值解。 解: 记 f (x, y) = y－ x y2, xn= nh (n = 0, 1, 2,···, N ) 由Euler公式得: yn+1 = yn + h( yn－ xn yn2 ) (n = 0, 1, ···,N) * 取步长 h = 2/10, 2/20, 2/30, 2/40, 用Euler法求解的数值实验结果如下. N 10 20 30 40 h 0.2 0.1 0.0667 0.05 误差 0.1059 0.0521 0.0342 0.0256 解析解: o —— 数值解 ---- —— 准确解 * y = f (x, y) 左矩形公式 用数值积分方法离散化常微分方程 ? 显式Euler法 求解常微分方程初值问题的数值方法可以由很多, 各种方法都可以看成使用不同的数值积分方法计算 * 梯形公式: 右矩阵公式: 隐式Euler法 * 预-校方法算法如下 k1 = f(xn , yn) , k2 = f( xn+1 , yn+ h k1) , 预报-校正方法(修正的Euler法): ? 预报 校正 * 预-校方法(h=0.2) 误差最大值0.0123 n 10 20 30 40 h 0.2 0.1 0.0667 0.05 误差2 0.0123 0.0026 0.0011 5.9612e-004 误差1 0.1059 0.0521 0.0342 0.0256 欧拉方法(h=0.2) 误差最大值 0.1059 * 基本概念 单步法: 即为了求得后一点xn+1上数值解yn+1,只要知道一点xn的数值解yn就可以了。 多步法: 即计算yn+1时,要用到前面多点的数值解yn, yn-1, ··· , yn-r等信息。 显式格式:数值解yn+1可以用xn和yn 解析表出。这类格式称为显式格式。反之称为隐式格式。对隐式格式而言,计算数值解yn+1形式上需要解方程(组) 。 显式方法易于计算, 但通常数值稳定性较差。而隐式方法通常稳定性较好, 但计算上显然很不方便。将显式方法和隐式方法联合使用是重要的思想方法。 * Tn+1=y(xn+1) - yn+1称为局部截断误差。设 yn= y(xn) 即 由Taylor公式 Euler公式: yn+1 = yn+ hf (xn, yn)