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

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