关于曲线拟合的方法探讨

  • 格式:doc
  • 大小:155.50 KB
  • 文档页数:3

关于曲线拟合的方法探讨
平书伟
(广东白云学院机电系,广州510450)
摘 要:在现代图形造型技术中,曲线拟合是一个重要的部分,是曲面拟合的基础。

现着重对最小二乘法、移动最小二乘法、NURBS 三次曲线拟合法和基于RBF 曲线拟合法进行比较,分析拟合方法的适用场合,从而为图形造型中曲线拟合的方法选用作出更好的选择。

关键词:离散点 曲线 拟合
0 引言
随着科技的不断进步,产品造型越来越复杂,因而对于造型的方法也提出了更高的要求。

由于经常要处理数据点,并对数据点进行拟合,因此曲线拟合就成了处理离散点成线的一种常用的手段。

曲线拟合的方法有很多种,各有各的优势。

本文将对常用的几种曲线拟合方法(最小二乘法、移动最小二乘法、NURBS 三次曲线拟合和基于RBF 曲线拟合)进行论述,分析拟合方法的适用场合,以便在针对具体情况时可以采用相应的拟合方法。

1 拟合方法论述
1.1 最小二乘法
最小二乘法通过最小化误差的平方和寻找数据的最佳函数匹配,是进行曲线拟合的一种早期使用的方法。

一般最小二乘法的拟合函数是一元二次,可一元多次,也可多元多次。

该方法是通过求出数据点到拟合函数的距离和最小的拟合函数进行拟合的方法。


2()f x ax bx c =++,计算数据点到该函数所表示的曲线的距离和最小。

即:
||||δ= (1)
对(1)式求导,使其等于0,则可以求出()f x 的系数a ,b ,c ,从而求解出拟合函数。

1.2 移动最小二乘法
移动最小二乘法是在最小二乘法的基础上进行了较大的改进,通过引入紧支概念(即影响区域,数据点一定范围内的节点对该点的拟合函数值有影响),选取适合的权函数,算出拟合函数来替代最小二乘法中的拟合函数。

从而有更高的拟合精度及更好的拟合光滑度。

1.2.1 移动最小二乘法的拟合函数
设拟合函数为()f x 在求解域Ω内的n 个节点i P (i =1、2、3、……、n ),则:
()f x =1()()m
i i i x K x α=∑=()()T K x x α (2)
式中,()x α为待求系数;()K x 为线性基函数。

一般令()K x =[1,,]T
x y ,m =3;求解过程可以参照文献[1],从而可求()x α,得到()f x 。

1.2.2 移动最小二乘法的算法流程
(1)将区域进行分段。

(2)对每个分段点进行循环:1)确定网格点的影响区域大小;
2)确定包含在网格点的影响区域内的节点;3)计算型函数;4)计算网格点的节点值。

(3)连接网格点形成拟合曲线。

1.3 NURBS 三次曲线拟合
NURBS 作为定义工业产品几何形状的唯一数学方法,是现代图形学的基础,因此NURBS 曲线拟合有着重要的实际意义。

NURBS 曲线的数学模型和数学方法可以参考文献
[2]。

本文采用VC 技术[3-4],利用OpenGL 的NURBS 曲线拟合函数,即可得到NURBS 曲
线。

1.4 基于RBF 的曲线拟合
RBF (Radial Basis Function ),径向神经网络是以径向基函数(RBF )作为隐单元的“基”,构成隐含层空间,隐含层对输入矢量进行变换将低维的模式输入数据变换到高维空间内,使得在低维空间内的线性不可分问题在高维空间内线性可分。

这是一种数学分析方法,具有较快的收敛速度、强大的抗噪和修复能力。

RBF 神经网络结构图如图1所示。

图1 RBF 神经网络结构图
各算法流程如下:
最小二乘法通过建立二次函数进行拟合。

建立拟合函数2
()f x ax bx c =++,求所有数据点与二次曲线的距离和最小的二次曲线,得到a 、b 、c ,从而得到二次曲线图像。

移动最小二乘法的流程是:
(1)NURBS 曲线拟合:确定节点矢量,本文通过弦长累加来确定节点矢量。

在NURBS 曲线拟合时,设置最前4个节点矢量的值相同和最后4个节点矢量的值相同,那么拟合的曲线将通过给定型值点的第一个点和最后一个点。

由于OpenGL 有现成的NURBS 曲线拟合函数,因此本文将借助VC 进行了编程,实现NURBS 三次曲线拟合。

(2)基于RBF 曲线拟合流程:本文将采用高斯函数作为RBF 函数的核函数[5]。

1)采
用K -均值法,确定聚类中心;2)按聚类中心分组;3)计算样本均值;4)重复2)、3),直到聚类中心不再变化;5)确定半径;6)调节输出层权。

2 实例验证
为了比较上述4种方法的优劣,本文采用1组数据,用4种方法进行拟合,然后比较拟合的情况,从而进行判断。

如果型值点经比较后,在型值点变化微小的情况,拟合的曲线将趋于平稳,就难以分出其优劣,因此在给定数据的时候(如表1),特地给出一个奇异点(第
3个点)。

通过对该表的数据点进行4种的拟合方法来比较各种方法的优劣。

(a)最小二乘法拟合(b)采用移动最小二乘法拟合
(c)NURBS三次拟合曲线(d)采用RBF进行拟合曲线
图2 不同方法拟合的曲线
3 结语
通过上述4个拟合的曲线可以得出:最小二乘法的精度最差。

使用RBF进行曲线拟合的精度最高,但不易用数学表达方式去表达,而NURBS曲线易用数学表达。

能用数学方式去表达RBF拟合的曲线,将使RBF拟合方法更具发展空间。

[参考文献]
[1] 曾清红,卢德唐.基于移动最小二乘法的曲线曲面拟合.工程图学学报,2004(1):84-88
[2] 朱心雄.自由曲线曲面造型设计.北京:科学出版社,2008
[3] 吕希奎,周小平.实战OpenGL三维可视化系统开发与源码精解.北京:电子工业出版社,2009
[4] 钱能.C++程序设计教程.北京:清华大学出版社,2005
[5] 张丽艳,王宏涛,李忠文.基于RBF神经网络的三角网格曲面孔洞修补.中国机械工程,2005,16(23):2072~2079
收稿日期:2011-03-22
作者简介:平收伟(1984.12—),男,江西广丰,硕士研究生,研究方向:曲线曲面造型。