构造二元切触插值函数的一种方法.pdfVIP

第29卷第2期 阜阳师范学院学报(自然科学版) Vo1．29．No．2 2012年 6月 JournalofFuyangTeachersCollege(NaturalScience) Jun．2012 构造二元切触插值函数的一种方法 剂 科 ‘ (阜阳师范学院数学与计算科学学院，安徽 阜阳 236037) 摘 要：文章利用埃米特插值基函数的方法，构造了一种矩形网格上的二元切触插值函数，并给出误差估计。最后通 过数值实例，说明该方法具有计算量低，构造过程公式化，便于编程的特点。 关键词：二元切触插值；误差估计：插值公式；二元埃米特插值 中图分类号：0241．3 文献标识码：A 文章编号20I2)O1．011-03 A method ofconstructingbivariateosculatoryinterpolation function JINGKe (Schoolo．1MathematicsandComputationalScience，FuyangTeachersCollege，FuyangAnhui236037，China) Abstract：Constructingtheosculatoryinterpolationfunctiononrectangulargrids，andgivingtheerrorestimatebymeansof themethodofhermiteinterpolationbasisfunction，throughthenumericalexample，weillustratethealgoritbm needlesscomputation， thecourseofconstructingfunctionisfommlary，anditiseasytorealizeoncomputer． Keywords：bivariateosculatoryinteprolation；errorestinmtes；interpolationformula；bivariatehermiteinteprolation 0 引言 (D|( )) = 矿 (2) yJ 二元切触插值 问题是一元切触插值问题的自 其中M=0，1，…，，n一1；=0，1，…，n，一1． 然推广，由于插值区域的复杂性和节点分布的多样 并记 ∑m=m+1，∑ =n+1。 性，已有文献多是针对矩形网格上的多项式插值问 题进行研究，如文[1．2]，下面给出二元切触插值问 文[3]利用切触拉格朗日插值算法，构造出了 题： 矩形网格上切触插值公式，但是该公式计算较为复 称 R 中的点集 杂。文[4]利用广义范德蒙矩阵的方法给出了一 兀 ={(。，Y)／(Y)∈R， 种插值公式，但求解广义范德蒙矩阵及其逆矩阵时 计算量也较大。本文利用一元埃米特插值基函数 i=0，l，…，S；=0，1，…，￡} (1) 的方法，构造 出了一种针对 171一1：1，凡，一1=1的 为矩形网格。 二元切触插值函数，并通过数值实例，说明了该方 所谓的二元切 触插值 就是求 二元多项 式