全局连续的分段最小二乘曲线拟合方法!.pdf

! #$$ 年$$ 月 重庆师范大学学报（自然科学版） %’( #$$ 第) 卷 第* 期 +,-./0 1 23.456.4 %-7/0 8.6’9-:6;（%/;,-/0 =69.9 ） ?0( ) %( * @AB ：2%CB ：D#E $$*DF %( #$$$$$# ( $D#G( ##H 全局连续的分段最小二乘曲线拟合方法! 侯超钧，曾艳姗，吴东庆，杨志伟 （仲恺农业工程学院计算科学学院，广州D$#D ） 摘要：针对传统最小二乘的曲线拟合方法不适合单一应用在具有复杂结构的试验观测数据中，本文提出一种满足全 局连续性约束的多分段区间的最小二乘数据拟合方法。通过把每个相邻分段点上要求拟合连续的约束条件转化成 一个矩阵等式! I # ，建立一个只包含线性等式约束的最小二乘模型76. J # ，最后通过应用拉格朗日的乘数 ! ! ! 方法推导出最小二乘解! 。本文的拟合方法在分段点上具有良好的拟合效果并满足全局连续，模型系数求解具有简 单的显式表达式，易于编程数值计算。 关键词：曲线拟合；最小二乘；分段拟合；拉格朗日乘数法 AK $( 文献标志码：L! ! ! 文章编号：$*ME **HG（#$$）#*E ##KKE #D 中图分类号： ! ! 在工程实践与科学实验中，常常需要从一组带噪声的试验观测数据（$ ，# ），% I $， ，…， 中找出自变量 % % $ 与因变量# 之间隐含的函数关系，一般采用数据拟合的办法来产生近似函数# I’ （$ ）。其中，线性最小二 乘方法是解决曲线拟合问题的常用方法，通过采用一组简单合适的、线性无关的基函数来逼近试验数据，可 ［$］ 以有效得出总体经验误差最小的拟合函数’ （$ ） 。 对于具有多个显著局部特征的试验观测数据，在全局定义域上应用一组基函数，将会得到性能较差的拟 合函数。一种解决手段是采用分段曲线拟合，在每段区间上进行局部最小二乘拟合。然而，拟合函数在区间 分段点上不一定连续，在相邻区间边界附近的拟合效果不理想［(K ］。其中，文献［D(* ］提出多项式基函数的 全局连续拟合方法，但只限于 个分段区间。文献［M ］给出多分段区间、全局连续的曲线拟合方法，但基函 数只限于一次多项式。文献［) ］提出分段区间重合的拟合方法，由每K 个数据点决定一个三次曲线，但分区 间太密，不适用于密集的数据拟合。目前的研究方法还没有一个可以描述全局连续的分段曲线拟合的有效 模型，以及相应的有效求解方法。 本文提出一种全局连续的多区间最小二乘曲线拟合方法，把全局连续性约束的分段拟合问题转化为带 等式约束的误差最小化模型，通过拉格朗日乘数法推导出模型的最小二乘回归系数。本文方法在实际应用 中具有以下优点： $ ）分段点上具有良好拟合效果，符合实际数据的连续性变化要求，从而避免拟合函数在分 段点上的二义性； ）模型的回归系数求解具有简单的显式表达式，易于数值计算；G ）方法对分段区间个数与 观测数据规模的限制较小，在各个分段区间上可以采用不同的基函数组，容易编程实现。 $ 两段区间的全局连续的最小二乘模型