计算方法第八章.ppt

类似地，如果取 作为 节点，可得亚当斯(Adams)隐式公式 局部截断误差： 进一步，如果我们将恒等式中的积分区间改为 ，并在此区间上用辛普森公式，可得 1.5 绝对稳定性 一个常微分方程数值解法称为绝对稳定，是指在某一步（如第一步）产生的误差（如计算机的存储误差），在计算中会逐步减小。 称方程 为试验方程，设计算步长为 h ，设在计算开始时产生误差（存储误差），此误差在以后会逐步减弱，我们称该算法相对于 是绝对稳定的， 的全体称为该算法的绝对稳定域。 欧拉法的绝对稳定区域 后退欧拉法的绝对稳定区域 1.6 局部截断误差的实用估计 （1）用两种阶数相同的算法求解，计算出n+1步的近似值，从而得到局部截断误差估计。 （2）用同样的公式，用不同步长计算出n+1步的近似值，从而得到局部截断误差估计。 1.7 隐式法 隐式法具有较好的绝对稳定性！ 只不过在使用隐式法的时候，需要进行迭代，或者使用预测－校正计算格式。 第二节 泰勒级数法与龙格－库塔法 对于方程： 取计算步长为 h ，则 ，将函数进行泰勒展开 如函数 y( x )有p+1阶导数，容易得到p阶泰勒级数展开法： 公式中的导数用下面公式计算： 例子： 步长0.1 步长0.01 龙格－库塔法： 对于常微分方程的数值解法，一个关键在于选择精度高的算法计算下面公式中的积分。 要高精度的计算积分，常用的方法是适当增加计算节点，不妨 设用 m 个节点计算上面积分，节点为 则积分为 将积分改写为： 则得公式： 这样的公式称为显式龙格－库塔公式。 确定二阶 R－K 法: 系数满足： 此为四个未知数的三个方程，任意取 ，得 特别，取 ，得到 通常也称为变形欧拉法，也常写为 它具有二阶精度，也称为二阶二级 R－K 方法。 计算方法第八章 本章主要研究常微分方程初值问题的数值求解： 通常，假设函数 f 关于第二个变量满足李普希茨条件（L条 件），即为存在常数 L > 0，使得 第一节 一般概念 1.1 欧拉法及其简单改进 方法：选择适当的节点，用差分近似微分，将方程离散化， 从而求在这些节点上的函数值。 欧拉 方法 例子： 下面我们分别取步长为0.1与0.01进行计算， 计算结果显示在下面的图中。 步长为0.1的计算结果。 步长为0.01的计算结果 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 0.99 0.37158 0.36973 1 0.36788 0.36603 DOUBLE PRECISION h,y(0：100) OPEN(20,FILE='OUTPUT.DAT',STATUS="UNKNOWN") h=1.0/100 y(0)=1.0 do 10 i=1,100 y(i)=y(i-1)*(1.0-h) write(20,*) i*h,y(i) 10 continue END 精度更高的欧拉公式： 方法：选择计算导数值 的精度更