《地球物理计算方法》第三章 微分方程差分法——Adams方法.pptVIP

《地球物理计算方法》

第3章微分方程差分法——Adams方法

设计思想由于龙格-库塔法在计算yn+1时,需要预报[xn,xn+1]区间上的若干斜率值。比如四阶龙格-库塔法，每一步都需要预报4个K值，即求4次函数f的值，计算量大。如果考虑利用yn+1之前的一系列“老”节点xn,xn-1，…上的斜率值，则可以大大减少计算量。亚当斯法属于线性多步法龙格-库塔法(4阶)亚当斯法(4阶)

特别地，格式与隐式是一阶方法.利用的信息，得到n阶的显式及隐式基本思想：利用xn，xn-1，xn-2…上的斜率值减少计算yn+1的计算量或提高精度。

4取

方法：取合理的λ，使上述格式具有二阶精度

——二阶Adams格式1、显式Admas方法

5选取参数λ=-1/2，该格式具有二阶精度.对上式进行展开,得到：将上式与泰勒展开公式进行对比：

6这样导出的计算格式?

7?其中他们都是显式的，算法简单。缺点：利用了n,n-1,…近似n+1点导数，精度差，加入n+1的导数信息，即采用隐式算法。

8?

9选取参数λ=1/2，该格式具有二阶精度隐式Admas格式.对上式进行展开，得到：将上式与泰勒展开公式进行对比：梯形公式

10同样，可以得到三阶隐式Admas格式.同样，可以得到四阶隐式Ad