计算方法第八章.pptVIP

进一步，如果我们将恒等式中的积分区间改为 ，并在此区间上用辛普森公式，可得 第三十页，共五十八页，2022年，8月28日 1.5 绝对稳定性 一个常微分方程数值解法称为绝对稳定，是指在某一步（如第一步）产生的误差（如计算机的存储误差），在计算中会逐步减小。 称方程 为试验方程，设计算步长为 h ，设在计算开始时产生误差（存储误差），此误差在以后会逐步减弱，我们称该算法相对于 是绝对稳定的， 的全体称为该算法的绝对稳定域。 第三十一页，共五十八页，2022年，8月28日 欧拉法的绝对稳定区域 后退欧拉法的绝对稳定区域 第三十二页，共五十八页，2022年，8月28日 1.6 局部截断误差的实用估计 （1）用两种阶数相同的算法求解，计算出n+1步的近似值，从而得到局部截断误差估计。 （2）用同样的公式，用不同步长计算出n+1步的近似值，从而得到局部截断误差估计。 1.7 隐式法 隐式法具有较好的绝对稳定性！ 只不过在使用隐式法的时候，需要进行迭代，或者使用预测－校正计算格式。 第三十三页，共五十八页，2022年，8月28日 第二节 泰勒级数法与龙格－库塔法 对于方程： 取计算步长为 h ，则 ，将函数进行泰勒展开 如函数 y( x )有p+1阶导数，容易得到p阶泰勒级数展开法： 第三十四页，共五十八页，2022年，8月28日 公式中的导数用下面公式计算： 例子： 第三十五页，共五十八页，2022年，8月28日 步长0.1 步长0.01 第三十六页，共五十八页，2022年，8月28日 龙格－库塔法： 对于常微分方程的数值解法，一个关键在于选择精度高的算法计算下面公式中的积分。 要高精度的计算积分，常用的方法是适当增加计算节点，不妨 设用 m 个节点计算上面积分，节点为 则积分为 第三十七页，共五十八页，2022年，8月28日 将积分改写为： 则得公式： 这样的公式称为显式龙格－库塔公式。 第三十八页，共五十八页，2022年，8月28日 确定二阶 R－K 法: 系数满足： 此为四个未知数的三个方程，任意取 ，得 第三十九页，共五十八页，2022年，8月28日 特别，取 ，得到 通常也称为变形欧拉法，也常写为 它具有二阶精度，也称为二阶二级 R－K 方法。 第四十页，共五十八页，2022年，8月28日 例子 第四十一页，共五十八页，2022年，8月28日 第一页，共五十八页，2022年，8月28日 本章主要研究常微分方程初值问题的数值求解： 通常，假设函数 f 关于第二个变量满足李普希茨条件（L条 件），即为存在常数 L 0，使得 第二页，共五十八页，2022年，8月28日 第一节 一般概念 1.1 欧拉法及其简单改进 方法：选择适当的节点，用差分近似微分，将方程离散化， 从而求在这些节点上的函数值。 欧拉 方法 第三页，共五十八页，2022年，8月28日 例子： 下面我们分别取步长为0.1与0.01进行计算， 计算结果显示在下面的图中。 第四页，共五十八页，2022年，8月28日 步长为0.1的计算结果。 第五页，共五十八页，2022年，8月28日 步长为0.01的计算结果 第六页，共五十八页，2022年，8月28日 0.01 0.99005 0.99 0.1 0.90484 0.90438 0.2 0.81873 0.81791 0.3 0.74082 0.7397 0.4 0.67032 0.66897 0.41 0.66365 0.66228 0.59 0.55433 0.55268 0.6 0.54881 0.54716 0.9 0.40657 0.40473 0.91 0.40252 0.40068