基于移动最小二乘法的曲线曲面拟合_曾清红.pdfVIP

2004年 工 程 图 学 学 报 2004 第1期 JOURNAL OF ENGINEERING GRAPHICS No.1 基于移动最小二乘法的曲线曲面拟合 曾清红， 卢德唐 （中国科学技术大学力学和机械工程系，合肥 230027 ） 摘 要：建立了一种基于移动最小二乘(Moving Least-Squares， MLS) 法的曲线曲 面拟合方法。这种方法对传统的最小二乘(LS) 法的作了比较大的改进，使生成的曲线曲面具 有精度高、光滑性好等许多优点。详细介绍了移动最小二乘法的原理、应用和特点，并且给 出了使用移动最小二乘法进行曲线曲面拟合的程序设计流程。最后给出了曲线拟合和空间散 乱数据曲面拟合算例，将拟合结果与最小二乘拟合结果作了比较，分析了MLS 拟合曲线曲 面的光滑性和拟合质量，表明了该方法的优越性和有效性。 关 键 词：计算机应用；移动最小二乘法；曲线曲面拟合；权函数 中图分类号：TP 391 文献标识码：A 文 章 编 号：1003-0158(2004)0 1-0084-06 曲线曲面拟合是一种古老而常用的技术，在工 [4] Belytschko 将其应用于无网格方法中 。移动最小 程、实验、统计和计算机图形等方面有着广泛的应 二乘法与传统的最小二乘法相比，有两个比较大的 用。通常，经过测量或者采集可以得到一组离散数 改进： 据点(x , y ), i=1, 2, …, n ，这里x 为坐标值。由于这 （1）拟合函数的建立不同。这种方法建立拟 i i i 些数据点并非完全精确，而且函数y =f (x) 的表达式 合函数不是采用传统的多项式或其它函数，而是由 预先无法知道，需要在给定的函数类f 上根据这些 一个系数向量a(x)和基函数p (x)构成，这里a(x)不 离散数据作出逼近曲线( 曲面) 。因为离散数据有误 是常数，而是坐标x 的函数。 差，并不要求逼近曲线( 曲面)经过数据点，而只是 （2 ）引入紧支（Compact Support ）概念，认 要求逼近曲线( 曲面) f (x) 的误差的某个指标达到最 为点x 处的值y 只受x 附近子域内节点影响，这个 小。传统的曲线( 曲面)拟合方法一般使用最小二乘 子域称作点x 的影响区域，影响区域外的节点对x [1], [2] 法 ，通过使误差的平方和最小，得到一个线性 的取值没有影响。在影响区域上定义一个权函数 方程组，求解线性方程组就可以得到拟合曲线( 曲 w(x) ，如果权函数在整个区域取为常数，就得到传 面) 。如果离散数据量比较大、形状复杂，还需要 统的最小二乘法。 进行分段(分块)拟合和平滑化，这在实际中往往带 这些改进能够带来许多优点，减缓或解决传统 来一定的困难。 曲线曲面拟合过程中存在的困难。可以取不同阶的 笔者使用移动最小二乘(MLS)法建立一种新 基函数以获得不同的精度，取不同的权函数以改变 的曲线( 曲面)拟合方法，这种方法能够克服以上困 拟合曲线( 曲面) 的光滑度，这是其它拟合方法无法 难，并且还具有许多其它优点。Lancaster 和 做到的。 Salkauskas 最先在曲面生成中使用了MLS[3] ，后来 收稿日期：2002-10-21 基金项目：国家自然科学基金资助项目；国家973 资助项目（G199903280