9-10.09级8.2euler方法.pptVIP

* §2 Euler方法 考虑问题 特点 简单，精度低. 先以Euler方法为例介绍离散化途径、 数值解法中的基本概念、 术语、加速方法等. 2.1 显式Euler方法 (折线法) 1 显式欧拉公式 设节点为: 则(1.1),(1.2) 的显式Euler方法为 步长 其中 推导公式 (1) Taylor展开法 将 在 点进行Taylor展开 忽略 这一高阶项, 由于 因此 分别用 近似 得 结合初值条件 即得(2.1). (2) 向前差分近似微分法 向前差商 近似 得 将近似号改为等号， 近似 并结合初始条件即得(2.1). (3) 左矩数值积分法 将(1.1)两端从 到 积分,得 数值积分采用左矩形公式 , 即 由初始条件亦得(2.1). o 用 近似 得 2 几何意义 (1.1),(1.2)的解曲线 过点 具有斜率 从 出发 以 为斜率作直线段,交 于点 此时 从 出发, 以 为斜率 于点 此时 以次类推.最终到达点 这样得到了一条折线 用折线 作为(1.1),(1.2)解曲线 的近似曲线 ， 折线法. o 过 的解曲线如图示, 具有斜率 作直线段,交 过 的解曲线如图. 它在点 的右侧具有斜率 与(1.1)过 的解曲线相切. 折线 称为欧拉折线，所以欧拉方法又称为 2.1 隐式Euler方法和梯形方法 (Euler方法的改进) 将 在 点进行Taylor展开 忽略 项, 用 分别近似 可得隐式欧拉方法: 1 隐式 Euler方法 在 点展开 显式Euler方法 或 说明 (1)隐式 Euler方法也可用 近似微分 或者用右矩数值求积公式来建立. 向后差商,即 (自己推) (2)隐式 Euler方法(2.5)是关于 的方程，若求 需要 解方程(2.5). 2 梯形方法 取平均,得 当 三次连续可微时, 忽略 项,用 得梯形方法: ? 与 因为(2.6)式也可通过梯形数值求积公式 其中 忽略数值求积余项来建立, 表示关于 的全微分. 说明 梯形方法也是隐式方法,要求 ,需解方程(2.6). 由 分别近似 ? 3 Euler方法、隐式Euler方法、梯形方法与单步法计算公式的 显式单步法 对应关系 隐式单步法 显式 Euler方法 隐式Euler方法 梯形方法(隐式) 4 迭代函数解的唯一性及收敛性 (以(2.6)式为例) 当 时,以y 为变量的函数 在 上关于y 满足Lipschitz条件,且L－常数为 , 从而由第七章 定理3,压缩不动点定理得方程 若 在 上满足基本条件, 关于 的Lipschitz常数为 ；同理只要 (2.6)确定了唯一的 L时,只要 由 Th.3, (2.5)确定了唯一的 且从任意初始值出发,迭代(2.5),(2.6) 都收敛到 . 有唯一不动点 事实上， 在实际计算中希望有较好的 , 用较少的迭代步(次数), 取得 称 为 的m次迭代 最常用的方法之一是先用显式Euler方法所得的 为 量较大 ,往往取 作为 来用. 有唯一不动点 , 都收敛到 . 而且从 出发,迭代 更精确的 (足够精度的 ). 而在实际计算中 的计算 改进. (初始值) ,再由梯形方法改进一次. 即是预估－校正Euler方法 (或称为改进Euler方法) . 2.3 预估－校正Euler方法 (显式Euler公式) (梯形公式) 或写成 或写成 例1 分别用显式Euler方法，梯形方法和预估－校正Euler方法解初值问题 解 取 h =0.1， Euler方法为 梯形方法为 预估－校正Euler方法 预估－校正方法 梯形方法 Euler方法 6.6×10－4 1.368514 3.1×10－4 1.367573 1.9×10－2 1.348678 1.0 6.6×10－4 1.307228 3.1×10－4 1.306264 1.9×10－2 1.287420 0.9 6.5×10－4 1.249975 3.0×10－4 1.249019 1.9×1