CFD2020-第10讲-常微分方程数值解法.pptVIP

Copyright by Li Xinliang * 四阶Runge-Kutta方法 四阶Adams外插方法 四阶四步方法(Taylor) 第一步计算 需要用到 。但初始条件只给出 ， 因此需要用其它方法提供另外k-1个值 。 计算量 稳定性 可能问题？ 多步法： RK方法计算4次函数值，多步法1次 RK方法优于多步法 Copyright by Li Xinliang * Admas外插方法：显式 内插方法：隐式 预测-校正法 预测公式（显式） 校正公式（隐式） Hamming 方法 迭代 Copyright by Li Xinliang * Lipschitz条件： 则称该数值方法是收敛的。 假设函数f(x,y)在区域D上连续，对区域D上的任意两点， 和 ，存在正的常数L(称为Lipschitz常数)，使得 记 则称函数f(x,y)满足Lipschitz条件。 收敛： 对 中任意固定的点 ，对数值解 ，有 相容： 步长趋于0时，数值方法的截断误差也趋于0（至少是一阶的） Copyright by Li Xinliang * 试验方程定义稳定： 稳定： 存在正常数 和 ，使得用任意步长 以及任意两个初值 和 计算得到的两组数值 和 ，成立 则称该方法是稳定的。 误差的增长有界 对初值问题 若一个数值方法在某一步由计算产生的误差在以后各步中逐步减少，则称该方法关于 是稳定的。 记 或：当 时， ，则称用步长h的这个数值积分过程是稳定的。 稳定区域： 使数值方法稳定的 的集合。 Copyright by Li Xinliang * 研究微分方程数值解的计算稳定性问题： 具有广泛的代表性，线性常系数常微分方程组通过线性变换可得到这种情形 线性非齐次微分方程组数值解的计算稳定性与齐次方程组是一样的 变系数的常微分方程组，工程上常将系数固态化后研究某个时刻邻域的性质 非线性的常微分方程组，在研究某个解邻域的误差时，可经过线化，得到误差所满足的线性方程组 试验方程： Copyright by Li Xinliang * Euler方法： 从 出发，假定在第m步按上式计算得到的值记为 ，在此基础上计算下一步的值 ， 可见，当且仅当 时，第m步产生的误差在以后各步中逐步减少，从而得到Euler方法的绝对稳定区间为： 稳定 Copyright by Li Xinliang * p阶显式Runge-Kutta的稳定性函数： 隐式Runge-Kutta的稳定性函数： RK方法的稳定区域： Euler: (-2,0) 2-RK: (-2,0) 3-RK: (-2.513,0) 4-RK: (-2.785,0) 显式RK方法： Copyright by Li Xinliang * 一般k步多步法（改变形式） 特征方程： k步多步法的稳定区域： 令 称 为线性k步法的特征多项式 线性k步法稳定(或称零稳定)的充分必要条件是 满足特征根条件： 的所有根均在单位圆内，并且在单位圆周上的根只能是单重根。 Copyright by Li Xinliang * 稳定区域的计算： 试验方程 线性k步法 特征方程 对于给定的 ，若特征方程的所有根 的模都小于1，则称线性k步法关于 绝对稳定。若对所有 ，线性多步法都绝对稳定，则称 为绝对稳定区间。 Copyright by Li Xinliang * 稳定区域的计算： 取 ，解出对应于每个 值的根 ，根的全体就是稳定区域的边界曲线。 共轭性质： Adams显式格式 1阶(Euler): (-2,0) 2阶: (-1,0) 3阶: (-6/11,0) 4阶: (-3/10,0) Euler: (-2,0) 2-RK: (-2,0) 3-RK: (-2.513,0) 4-RK: (-2.785,0) K步多步法: 令 Runge-Kutta方法: Adams隐式格式 2阶: (- ∞,0) 3阶: (-6,0) 4阶: (-3,0) 5阶: (-90