高维李查逊外推算法研究和应用.pdfVIP

高维李查逊外推算法的研究和应用 林柏洪 1 张权 2 杨小远 3 （1，2 北京航空航天大学物理科学与核能工程学院，北京，100191） （3 北京航空航天大学数学与系统科学学院，数学、信息、行为教育部重点实验室北京，100191） 摘要：本文在李查逊外推法算法的基础上，提出了广义李查逊外推法算法，在此基础上构造 了高精度高阶导数的数值计算法，实验结果表明，本文提出的算法有效性。 关键字 外推，Tayolor公式，导数的数值计算 中图分类号：0177.2 文献标示码:A 1. 李查逊外推法的简介 李查逊外推法是科学计算领域提高算法的精度重要方法，广泛应用于数值积分，有限元 和偏微分方程数值解等领域。下面介绍理查逊外推法 定理 1 李查逊(Richardson)外推法 设步长为 h 的算法Ι(h) 去逼近量 ，若 和Ι(h)之间的截断误差有渐近展开式： Ι Ι 1 1 Ι−Ι (h) C hp 1 =+C hp 2 ++C hp k + O(hp 1 ) 1 1 2 k 其中p k p k −1 p 1 0 Ι (rh) −rp m Ι (h) 定义Ι (h) m m ，则Ι (h) 逼近Ι的截断误差阶为hp m+1 。 m+1 1−rp m m+1 从以上算法中，我们可以看到，此算法局限于单变量函数的计算。那么对于多元函数而言， 理查逊外推法是否一样成立，下面我们将理查逊外推法推广为广义理查逊外推。 2. 广义李查逊外推 这节我们给出广义广义李查逊外推算法以及理论分析。首先我们给出多元函数的 多元函数的 Taylor 定理 定理 1 多元函数的 Taylor 公式 假设 n m+1 D ⊂R 是一凸域，f ∈C （D）。x (x , x ,x ) 和 1 2 n → x +h (x =+h , x +h ,x +h ) 是 D 中两点，则必存在θ ∈（0，1），使得： 1 1 2 2 n n → → → m Dα f (x ) α Dα f (x +θh ) α f (x +h ) ∑∑ h =+Rm ，其中Rm ∑ h ，这里 α! α! k 0 |α| k |α| m=+1 α +α +α α → ∂ 1 2 n f → α α α α D f (x ) α α α (x ),h h1 1h2 2 ...hn n 其中 1 2 n ∂x ∂x ∂x 1 2 n |α| α =+α +α ,α! α !α !α ! 1 2 n 1 2 n 基金项目：北京市精品课程建设项目资助 作者简介： 林柏洪，男，1989 年 11 月出生，现为北京航空航天大学物理科学与核能工程学院学生。张权， 男，1991 年 1 月出生，现为北京航空航天大学物理科学与核能工程学院学生。杨小