基于散乱数据的球面自然样条插值法.pdfVIP

Vd．27 第27卷第5期 成都信息工程学院学报 No．5 OFCHENGDUUNIVERSITYOFINFORMATIONTECHNOLOGY Oct．2012 2012年10月 JOURNAL 文章编号：1671．1742{2012}05．0520．05 基于散乱数据的球面自然样条插值法 王丽丽， 徐应祥 (中山大学新华学院，广州510520) 摘要：针对球面散乱数据插值问题，结合三元多项式自然样条方法的基本原理，提出一种新的逼近球面上散乱 数据的自然样条插值方法，利用新的方法对给定区域内的三元函数进行数值逼近实验。将数值实验误差与理论误 差对比(5％)，结果验证孑这种方法在球面散乱数据逼近问题上能够有效的使最大误差控制在2％以内。研究结果 表明，这是一种比较有效的球面散乱数据逼近方法。 关键词：计算数学；数值逼近；球面散乱数据；自然样条；插值 中图分类号：0241．3 文献标志码：A O 引言 随着信息数字化技术的发展，信息数据在数量和维度上都发生了显著地变化。越来越多专家和学者的研究 于此，3D空间数据模型中球面数据逼近就成了最基本的研究问题【1,2J。 球面数据逼近的研究工作，不仅有重要的理论意义，而且有相当地实际应用背景。近年来，关于球面数据逼 近问题的研究，主要表现在以下5个方面：地球物理及地震学研究方面的应用；物理大地测量学中的重力场的测 量；气象学、海洋学等相关环境资源的数据逼近；地球气温、气压、臭氧或磁场等方面的研究；医疗方面对某种脑电 图的分析研究。至今，专家与学者们针对球面数据逼近问题已经提出了多种不同的有效方法∞．7J。但球面散乱 数据插值问题仍没有非常有效的解决方法。 针对球面上散乱数据插值问题，基于三元多项式自然样条[8—10,1213l，构造了一种新的三元自然样条插值方 法，并进行了数值实验．数值实验结果表明：这是一种比较有效的方法。 1球面散乱数据的三元自然样条插值 给定单位球面B(O，1)={(z1，z2，z3)l =1，…，N}，其值为伽i，净1，…，N．设X=H(声’(n)为Sobolev空间，其中n=Hai，bi]为包含给定散乱数 据的方型域，P=P1+P2+夕3，p1，户2，夕3均为给定的正整数． 边界条件 f“(‘1，户2，o’(口1，z2，z3)=0 (1) 心’_、(¨z，z，宁 、。 l“(il，i2，≯3’(口1，n2，z3)=0 【“(pl’Pz，i3)(z1，z2，口3)=0 其中i女=0，1，…，么一1；k=1，2，3。 记Z=RN为N维欧氏空间，W=(W1，…，wN)。A：x—z是一个插值算子，满足： Au=(U(z11，X21，X31)，…，“(z1N，z2．N，X3N))． 收稿日期：2012．04—16 第5期 王丽丽等：基于散乱数据的球面自然样条插值法 521 题．该问题的解称之为三元球面散乱自然样条，简称球面自然样条。 问题P寻找一个函数a(x1，X2，X3)∈X，满足 1|乃lI2=rain…而||2} (2) Ut^．m2Ⅲ 这里，l|TuIl 2=…(“(户t，pz，户s’)2dxldx2dx3，带形如(1)的自然边界条件。 《 定理1球面自然样条叮(zl，z2，X3)具有如下的显式及紧凑格式的表达式： N