高等数值算法与应用（九） - Wenjian YU.PDFVIP

高等数值算法与应高等数值算法与应用 （九） AdvancedAdvanced NumericalNumerical Alggorithms Apppplications 计算机科学与技术系计算机科学与技术系 喻文健喻文健 内容概要内容概要 常微分方程初值问题 基本概念基本概念 ODE的稳定性的稳定性((解的稳定性解的稳定性)) 简单的数值解法与算法稳定性、截断误差分析 实用的数值解法实用的数值解法 MatlabMatlab有关命令和应用例子有关命令和应用例子 Wenjian Yu 2 OrdinaryOrdinary DifferentialDifferential EquationsEquations (I)(I) Wenjian Yu 3 本节主要内容 常微分方程(ODE)基本概念 初值问题初值问题(IVP)(IVP) 初值问题的稳定性(敏感性) 例子与Matlab演示 ODEODE-IVPIVP数值解法数值解法 欧拉方法 数值解法的稳定性 数值解法的截断误差与准确度数值解法的截断误差与准确度 数值解法的步长选择 Wenjian Yu 4 常微分方程常微分方程 要求解的是单自变量的函数(一个或多个) 自变量常为时间自变量常为时间tt, 未知函数为系统状态量未知函数为系统状态量yy 微分方程表示未知函数与自变量之间关系, 其中对 t 的最高求导阶数为常微分方程的阶数 ′′ ′′′′ ((kk)) kk阶阶ODEODE: g tt,y,y ,y , L,y 00 (( )) 一般可转化为一般可转化为显式显式常微分方程常微分方程 (k) ′ ′′ (k-1) 显式ODE: y f t, y, y , y , L, y (( )) 主要考虑显式ODE，特例是一阶、单个ODE 问题 y′′ ff t, y 例如：y′′ λy (( )) 高阶常微分方程高阶常微分方程→一阶常微分方程组阶常微分方程组 不失一般性，仅讨 论一阶常微分方程 组: y′ f t, y (( )) Wenjian Yu 6 牛顿运动第牛顿运动第二定律定律 加速度：