18. 常微分方程数值解法初步

数值分析



周维祺

一般理论

设定

例子

  • 人口增长:; 解:

  • 成对增长:; 解:
    在有限时间内就能增长到无穷


  • logistic曲线:;解:

常用来描述随时间变化的过程

解的局部存在性

  • 相对于连续,则

  • 一般可证满足一定条件时,解局部存在(在的某个邻域上存在)

  • 相关定理:Picard-Lindelöf, Peano, Caratheodory 等

  • 拓展阅读:Gronwall不等式、解对初值的敏感性

研究方法

  • 解析法:讨论解的存在性、唯一性、各种性质(光滑性、奇异性等)

  • 数值法:找出近似解,解法的收敛性、稳定性

  • 找出上的近似解

  • 计算只用到上一步的信息:单步法,否则称为多步法

欧拉法

设定

  • ,在上求解

  • 将区间等分,记



策略

计算公式

例子

  • ,精确解:



练习

,借助计算器,分别计算以下两个方程在处的解

解答

  • 精确的通解分别是

  • 满足的特解都是,因此

  • 计算得分别约等于

  • 第二个方程实际解得的是时的近似解

常微分方程数值解法的收敛性

,则称该方法收敛

欧拉法的收敛性

若在存在且连续,且存在常数使

对任意都成立,则欧拉法收敛

证明的主要步骤

  • 为第步的误差

  • 泰勒展开:
    欧拉法:



证明的主要步骤







  • 时,误差

局部截断误差和精度

设有数值方法

则称

为局部截断误差,若误差为阶,则称该法具有阶精度

局部截断误差

以真实值代入数值公式计算所得的误差

即相应数值积分的误差

欧拉法的局部截断误差和精度

欧拉法:

泰勒展开:

局部截断误差:,一阶精度

小结

  • 常微分方程数值解法的一般概念

  • 欧拉法的公式、计算

  • 欧拉法的收敛性和精度

  • 局部截断误差和精度