数值分析-李庆扬-第9章--常微分方程初值问题数值解法.pptVIP

2025年1月1日1例如：时，第9章常微分方程初值问题数值解法9.1引言微分方程：包含自变量、未知函数和未知函数导数或微分的方程。例如：，求定解条件：求解微分方程时，所附加的条件——定解问题。初始条件：给出积分曲线在初始时刻的值——初值问题。例如：时，边界条件：给出积分曲线在首末两端的值——边值问题。常微分方程：未知函数为一元函数。偏微分方程：未知函数为多元函数。

2025年1月1日2一阶常微分方程的初值问题：求解注意：——解函数、积分曲线；——微分函数。确定初值问题的解存在而且唯一：李普希兹条件。，

2025年1月1日3如果存在实数，使得称关于满足利普希茨条件，为的利普希茨常数。说明：条件可理解为解函数无限接近时，微分函数也无限接近。定理1设在区域上连续，且关于满足利普希茨条件，那么对任意常微分方程初值问题当时存在唯一的连续可微解。

2025年1月1日4关于方程的解对扰动的敏感性，有结论：定理2设在区域上连续，且关于满足利普希茨条件，设初值问题，，其解为，那么说明：①定理说明解对初值的敏感性，即初值不同，解也有差异；②解得敏感性与微分函数有关：当的利普希茨常数较小时，解对初值相对不敏感；当较大时，初值的扰动会引起解剧烈变化—病态问题；

2025年1月1日5数值解法：在一系列离散点上，求解近似值。“步进式”：顺着节点排列顺序，一步一步地向前推进。步长：常用等步长，节点为单步法：计算时，只用到前一点的值步法：计算时，用到前面点的值

2025年1月1日69.2简单的数值方法9.2.1欧拉法与后退欧拉法初值问题：解的形式：是通过点的一条曲线——积分曲线。特点：积分曲线上每一点的切线斜率为

2025年1月1日7尤拉方法：①将解区间离散化，选择步长，得到离散点：；②由切线，切线与交点：的近似值；③再由向前推进到，得到折线，近似。

2025年1月1日8任意折线：过点作直线，斜率，——欧拉方法假设初值，由此可逐次算出：

2025年1月1日9P281例1求解初值问题解：欧拉公式为，

2025年1月1日10局部截断误差：设前一步值准确，算下一步出现的误差假设：泰勒展开函数：局部截断误差：

2025年1月1日11后退的欧拉法：离散化：求解微分方程的关键，消除导数项，根本方法之一是用差商替代导数项。例如：——向前的欧拉公式〔显式〕

2025年1月1日12同理：——后退的欧拉公式〔隐式〕注意：①显式计算方便，隐式稳定性较好；②上式隐含，采用迭代法求解。

2025年1月1日13欧拉公式的另一种理解：将常微分方程改写对微分方程从到积分由积分左矩形公式得再以代替，以代替——向前的欧拉公式

2025年1月1日14对微分方程从到积分由积分右矩形公式得再以代替，以代替——后退的欧拉公式同理：

2025年1月1日15迭代法求解：后退的欧拉公式——逐步显示①先用尤拉格式，求