基于移动最小二乘法的曲线曲面拟合-Read.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 曲线曲面拟合是一种古老而常用的技术 在工 Belytschko 将其应用于无网格方法中[4] 移动最小 程 实验 统计和计算机图形等方面有着广泛的应 二乘法与传统的最小二乘法相比 有两个比较大的 用 通常 经过测量或者采集可以得到一组离散数 改进 据点(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 资助项目 G1999032805 作者简介 曾清红 1978- 男 湖南长沙人 硕士生 主要研究领域为大规模科学计算及科学计算可视化 万方数据 第1期