曲线拟合方法浅析

1、曲线拟合方法概述工业设计张静 1014201056引言：在现代图形造型技术中，曲线拟合是一个重要的部分，是曲面拟合的基础。现着重 对最小二乘法、移动最小二乘法、NURBS三次曲线拟合法和基于RBF曲线拟合法进行 比较，论述这几种方法的原理及其算法，基于实例 分析了上述几种拟合方法的特性，以 分析拟合方法的适用场合，从而为图形造型中曲线拟合的方法选用作出更好的选择。1曲线拟合的概念在许多对实验数据处理的问题中，经常需要寻找自变量和对应因变量之间的函数关 系，有的变量关系可以根据问题的物理背景，通过理论推导的方法加以求解，得到相应关 系式。但绝大多数的函数关系却很复杂，不容易通过理论推导得到相关的

2、表达式，在这种 情况下，就需要采用曲线拟合的方法 来求解变量之间的函数关系式。曲线拟合(Curve Fitting)，是用连续曲线近似地刻画或比拟平面上离散 点组所 表示的坐标之问的函数关系的一种数据处理方法。在科学实验或社会活动中，通过实验 或观测得到量x与y的一组数据对(Xjyi), i=1 , 2, 3，m，其中各Xi是彼此不同 的。人们希望用一类与数据的规律相吻合的解析表达式y=f(x)来反映量X与y之间的依 赖关系。即在一定意义下“最佳”地逼 近或拟合巳知数据。f(x)称作拟合函数，似的图 像称作拟合曲线。2 曲线拟合的方法2.1最小二乘法最小二乘法通过最小化误差的平方和寻找数据的最

3、佳函数匹配，是进行曲线拟合的 一种早期使用的方法一般最小二乘法的拟合函数是一元二次，可一元多次，也可多元多 次。该方法是通过求出数据点到拟合函数的距离和 最小的拟合函数进行拟合的方法令f(x)=ax 2+bx+c，计算数据点到该函数所表示的曲线的距离和最小即：i0对上式求导，使其等于0,则可以求出f(x)的系数a,b,c，从而求解出拟合 函数。2.2移动最小二乘法移动最小二乘法在最小二乘法的基础上进行了较大的改进，通过引入紧支概 念(即 影响区域，数据点一定范围内的节点对该点的拟合函数值有影响)，选取适合的权函数，算出拟合函数来替代最小二乘法中的拟合函数从而有更高的拟合精度及更好的拟合光滑度。

4、2.2.1移动最小二乘法的拟合函数设拟合函数为f(x)在求解域Q内的n个节点Pj (i=1、2、3、n)，贝U:mf(x)= i(x)K i(x) = Kx) (x) i 1式中，ax)为待求系数；K(x)为线性基函数。一般令K(x) = 1,x,y t, m=3 ；求 解过程 可以参照文献，从而可求x),得到f(x)。2.2.2移动最小二乘法的算法流程将区域进行分段。对每个分段点进行循环：确定网格点的影响区域大小；确定包含在网格点的影响区域内的节点；计算型函数；计算网格点的节点值。连接网格点形成拟合曲线2.3 NURBS三次曲线拟合NURBS作为定义工业产品几何形状的唯一数学方法，是现代图形

5、学的 基础，因此NURBS曲线拟合有着重要的实际意义。NURBS曲线的数学模型和数学方法可以参考文献2。本文采用VC技术，利用OpenGL的NURBS曲线拟合函数，即可得到曲线。2.4基于RBF的曲线拟合RBF( Radial Basis Function)，径向神经网络是以径向基函数(RBF)作 为隐单元 的“基”，构成隐含层空间，隐含层对输入矢量进行变换将低维的模式输入数据变换到高 维空间内，使得在低维空间内的线性不可分问题在高维空间内线性可分。这是一种数学分 析方法，具有较快的收敛速度强大的抗噪和修复能力。RBF神经网络结构图如图1所示。图1 RBF神经网给结构图各算法流程如下：最小二乘法通过建立二次函数进行拟合。建立拟合函数f(x)=ax 2+bx+c，求所有数据点与二次曲线的距离和最小的二次曲线，得到a,b,c,从而得到二次曲线图像。移动最小二乘法的流程是：(1) NURBS曲线拟合：确定节点矢量，通过弦长累加来确定节点矢量。在NURBS曲线拟合时，设置最前4个节点矢量的值相同和最后4个节点矢量的值相同，那么拟合的曲线将通过给定型值 点的第一个点和最后一个点。由于OpenGL有现成的NURBS曲线拟合函数，借助VC进 行编程，实现NURBS三次曲线拟合。（2）基于RBF曲线拟合流程：采用