一种分段曲线拟合方法研究

1、湖南人文科技学院毕业论文一种分段曲线拟合方法研究摘 要：分段曲线拟合是一种常用的数据处理方法，但在分段点处往往不能满足连续与光滑.针对这一问题，本文给出了一种能使分段点处连续的方法.该方法首先利用分段曲线拟合对数据进行处理；然后在相邻两段曲线采用两点三次 Hermite插值的方法，构造一条连结两条分段曲线的插值曲线， 从而使分段点处满足一阶连续.最后通过几个实例表明该方法简单、实用、效果较好关键词：分段曲线拟合Hermite插值 分段点 连续Study on A Method of Sub-Curve FittingAbstract: Sub-curve fitting is a common

2、ly used processing method of data, but at sub-points it often does not meet the continuation and smooth, in allusion to to solve this problem, this paper presents a way for making sub-point method continuous. Firstly, this method of sub-curve fitting deals with the data; and then uses the way of t w

3、o points ' cubic Hermite interpolation in the adjacent, structures a interpolation curve that links the two sub-curves, so the sub-point meets first-order continuation; lastly, gives several examples shows that this method is simple, practical and effective.Key words: sub-curve fitting Hermite i

4、nterpolation sub-point continuous21、入刖曰数据拟合是一种重要的数据处理方法，其中最常用的是多项式曲线拟合.然而当数 据点较多时，多项式阶数太低，拟合精度和效果不太理想，要提高拟合精度和效果就需要 提高曲线阶数，但阶数太高乂带来计算上的复杂性及其他方面的不利.因此，如果只采用一种多项式曲线函数拟合较多的数据点，难以取得较好的拟合精度和效果.为有效地解决上 述问题，一般采用分段曲线拟合.以往的分段曲线拟合方法主要是针对在自然科学领域中 测量的数据而使用的拟合方法，这些数据的变化一般都遵循一定的规律.因此，在对这些测 量数据拟合时，传统的分段曲线拟合方法一般是先根

5、据主观经验对数据分段，然后进行拟合.但是对于有些实际问题的数据，比如社会、经济生活中的大量统计数据，这些数据变 化的机理一般非常复杂，往往不像物理定律那样有着严格的规律，所以变化的不确定性很 强.因此，传统的分段曲线拟合根据主观经验对数据进行分段的做法就显现出明显地不足 针对这种不足，国内外许多文献也讨论过，文献1研究的是最小二乘法在曲线拟合中的 实现，给出了最小二乘法在多元正交基函数拟合中的计算机实现方法，以常见的二次曲 线拟合为例说明了程序编制的要点，在实验的数据处理中具有实用价值；文献2 讨论分 段最小二乘曲线拟合方法，本文在一般最小二乘的基础上提出分段最小二乘曲线拟合的 方案，讨论了连

6、接分段拟合曲线的方法，并且给出分段最小二乘多项式拟合的计算方法； 文献4主要介绍基于最小二乘原理的分段曲线拟合法，在最小二乘的基础上，运用实测数据点的分段曲线拟合法，探讨相应的模型以及用不同类型的曲线拟合同时拟合数据点 的具体应用，对一实例，应用MATLAB编程设计，完成模型的求解、显著性检验等，可 以得到拟合精度比较高的拟合曲线，该方法原理简便，其模型易用MATLAB编程求解；文献5研究的是基于最小二乘法的分段三次曲线拟合方法研究，多项式曲线拟合是一种较常用的数据处理方法，但当数据点较多时，只采用一种多项式曲线函数拟合所有数据 点难以得到较好的拟合效果，针对传统分段曲线拟合方法中对数据点分段

7、时经验成分较 多的不足，提出了一种基于最小二乘法原理的分段三次曲线拟合方法，建立三次拟合曲 线方程，通过实际数据的检验，验证了该方法的拟合效果；文献6,7,8主要研究基于分段三次曲线拟合的广州周发案量预测，随着城市化进程的不断加快，城市人口不断增多， 广州市未来治安形势预警，支持政府部门和政法部门关于治安工作的决策，首先需要对 未来时期的发案量做出比较精确的预测，由于目前广州市方案量统计数据比较少，且发 案量受农历春节影响较明显，针对传统时间序列预测方法在此情况下应用不足，提出了 基于分段三次曲线拟合的周发案量预测模型，并给出了具体的建模、计算步骤，最后通 过实际数据的检验，证明了方法预测效果

8、较好；文献9提出了分段函数的光滑方法及其 在曲线拟合中的应用，在分析复杂实验数据时，采用分段曲线拟合方法，利用此方法在 段内可以实现最佳逼近，但在段边界上却可能不满足连续性与可导性.为了克服这种现象，本文主要研究一种能使段边界连续的方法，具有一定的理论和 实际意义.在前人的基础上，本文总结分段曲线拟合的方法与步骤，介绍了分段三次曲线 的拟合方法和两点三次Hermite插值，然后讨论如何利用Hermite插值方法使得分段拟合 曲线在连接点处满足连续方法，最后通过一些实例应用，表明本文所介绍的方法具有一 定的应用价值.1最小二乘曲线拟台1.1最小二乘法1令待求的未知量为a,&,，at，它们

9、可由n(n芝t)个直接测量必必,y通过下歹0函 数关系求得：yi = fi(ai,a2,，at) y2 = f2(ai,a2, ,at) y3 = f3(ai,a2, ,at) yn = fn(ai,a2, ,at)若aj为真值，由上述已知函数求出真值y、若其测量值为y"则对应的误差为二j = yj -yj,(j =i,2, n)职小二乘法可正重表小为：(i.i.i)对不等精度的测量，应加上各测量值的权重因子 Pj ,即:(i.i.2)' Pr、= min最小二乘法是在随机误差为正态分布时，由最大似然法推出的这个结论.它可使测量误差的平方和最小，因此被视为从一组测量值中求出一

10、组未知量的最可信赖的方法.i.2最小二乘多项式曲线拟合的基本原理i.2.i线性拟合原理将拟合函数取线性函数是一种简单的数据拟合方法，将数据点(Xi, f (Xi),(X2, f(X2), ,(Xm, f(Xm)确定线性拟合函数(i.2.i.i)'(x) = a bX称为对数据的线性拟合。对丁线性拟合问题，需要求函数S(a,bT (a + bQ y(i.2.i.2)的最小值点，该问题的几何背景是寻求一条直线，使该直线与数据表所确定的平面散点 的纵向距离的平方和最小，如图i.2.i-i所示.(图 i.2.i-i)由函数对两个变量求导得:心m笠=2£ca k =心m竺=2£

11、;的k旦a+bxka +bx”k)-y，x ,其余等丁零，得正规方程组也可将其矩阵形式写出来即："mma+£ xkb kmmZ xka+E xjk 1kmm七m Z xkkz!mk：2b =ykm k，xk yk4/ mW ykJmm£ X £ x2IkAkA加mzIk；xkykJ(1.2.1.3)(1.2.1.4)(1.2.1.5)解得a,b的值，将其代入(1.2.1.1)即可得到拟合线性函数1.2.2多项式拟合原理为了确定数据拟合问题，选用籍函数1,x,x2，xn作为函数类，贝U中(x) =a°+a1x+a x2 + anxn (n+1&l

12、t;m)(1.2.2.1)这就是多项式拟合函数.为了确定拟合函数甲(x) =a°+ax+a2x2 +anxn的系数，需要求解正规方程组mmmma。寸 xk-一 x& =、' ykm z kd Lkxk3o=1m+zkT2 xk4 +k=1k=1mnH+ 2xkankmm=£km(1.2.2.2)mzn xkaomn +xka1mm+十£ x：nan=£nxk ykk=1k=1k=1也可以用矩阵形式表示为fmm/ mm£ xkV n-xkzykkmk=1'a。'k =1mmmmw xk2£ xk寸n噂&#

13、163; xkazxk ykk=1k=1 a-k=1aak*ammm )m寸 nE xk一 n-E xk一 2n£ xkzn xkyk<k=1k =1k=1)lk=1J解得a。, a1，an即可，将其代入(1.2.2.1)即可得到拟合多项式.2分段曲线拟台2.1分段曲线拟合的基本原理3.拟合函数一般然后再应用最先根据实测数据分布的特点，确定分段数目以及相应拟合曲线类型 可选为多项式函数，因为在一定范围内，连续函数可用多项式任意逼近, 小二乘法原理求得各分段拟合方程的系数.基本步骤为：第一步：将数据点分段，确定基函数 中0(x),气(x),，气(x),第二步：根据题目要求，建立正

14、规方程组，第三步：解正规方程组，求出待定系数，第四步：写出拟合函数.下面以分段线性拟合与分段三次曲线拟合为例讨论分段拟合的基本过程 2.1.1分段线性拟合我们把给出的数据点分成k组N1,N2，Nk,即*、(X11, yi1)，( X12，yi2 )，( X1N1，yiN1 ).* . .* .*.(X21，y21 )，( X22，y22)，( X2N2，y2N2 ).* . .* .*.(Xk1，yk1 )，( Xk2，yk2)，( XkNk，ykNk )其中N1，N2,，Nk为每组数据的个数.首先考虑线性拟合这种简单的情形，对k组数据点分别应用最小二乘线性拟合，得到各组数据点所对应的近似线性

15、函数，g(X) =a+bX(命小小卬)g2(X) =a2 b X (*卬 £ x ： X2n2)gk(X)=ak bkX(Xjn % x £、)而在整个考虑的拟合区间上就得到了 k -1条直线段:现在就这k -1条直线段所在各 区间的左端点定义g (xm ) = g+3凡)，该函数就成为整个区间上的数据拟合函数.这就是 分段最小二乘线性拟合问题.然而有些数据组并不是每段都呈线性关系，如数据 (xyji =1,2，n，根据其散点 图却发现其前m个点较接近直线，后n-m个点呈现非线性关系，则可分两段拟合.分别 以一次多项式Y和n次多项式进行拟合，即Y=kx + b(2.1.1.

16、1)为了说明具体的方法，不妨选E的阶数为2,即Y2 =a0x2 +%x + a2(2.1.1.2)要保证在边界点(xym)连续光滑，所以存在两个约束条件 2kXm+b =aoXm+aXm+a2和k=2a°Xm+a，因此，式(2.1.1.1)和(2.1.1.2)的系数是相关的.解得 aa0xm，故式(2.1.1.1)为2Y =(2a°Xm +a)x aoXm +a2令S为最小二乘估计量，则mn2、2，2、2S = '(2 a0xm，a)Xia2'a0Xm yi 一 .二.(a0X|3ixi3 yi )(2.1.1.1)i 4im通过模型 =0 ; i=0,1,

17、2，可求得最小方差S的a0,a1,a 2的值，从而确定出式 ：ai与(2.1.1.2)中的回归系数.最后，通过Il mnZ (Yi -V)2 +£ (Ei 3)2i 4i -mr ：n -:L (" - y)2和F检验值F =(二)22 ,对回归方程进行显著性检验，式中22Yi =( 2血 xm + a) x a mx+ ;236i=a0X +ax + a2.当然，根据不同的数据，可分三段进行拟合，或根据不同的数据特点，采用多次曲 线拟合方式.2.1.2分段三次曲线拟合4,5设有N个数据乙,Z2,Z3,Zn .因为四个数据点可确定一条三次曲线，但在选取分段 点时，必须考虑分

18、段后相邻曲线必须连续，即边界点连续，因此用五个数据点拟合一条三次曲线.拟合方法：首先对数据进行一定的分段，将第一到第五数据分为第一段，再将第五到第九个数据分为第二段，将第九到第十三个数据分为第三段，依次类推进行分组，即前一段末尾的数据为下一段数据的首位，这样便保证了数据分段的连续性.然后再对个分段数据进行三次曲线拟合即可.令某段数据的三次拟合曲线函数为：Wt =a+bt+ct2 +dt3(t = -2,-1,0,1,2)可以将此曲线函数分解为奇偶两个函数：奇函数 v=bt+dt3和偶函数ut=a + ct2.下面应用最小二 乘法的基本原理求三次拟合曲线的系数6,由丁在每段数据中第一点和最后一点

19、均两次 参与拟合，因此，在求一段曲线的拟合方差时需要加权.按照平均分配的原则口，求方差1的权值史2=2= ,，=九0=气=1,得到该段曲线拟合的方差2(2.1.2.1)(2.1.2.2)(2.1.2.3)(2.1.2.4)S2=' "(W-Zt)2 17曲线表小为奇偶函数的形式如下Wt =Ut Vt,u_t =Ut,v_t = -Vt由(2.1.2.2)可以推导出下式.、1,、ut = 2伽 w_t), vt = 2 (W - w) 令 Zt =x +yt，x_t =x，y_t = % M1_1 _Xt"(Zt *曹=2(乙一*)因此拟合方差为2./_ 、2，、2S

20、t(Wt-Zt)t(Ut VtXt - yt)t瑚t瑚22 ,、2 一 ，、2=' ，t (ut xt ) - ._ - .t (vt - yt)(2.1.2.5)t /t/=S2即wt对乙的平滑可以看作是奇函数和偶函数分别平滑的叠加.从（2.1.2.5）式中可知奇函数拟合的方差.2S*t(xt -yt)2t / 2=，2人(btdt3 yt)2(2.1.2.6)t a=2(b d - y1)2 (2b 8d - y2)2令b d - y1 = 02b 8d - y2 = 0，解出b=(y2-2y)6 d =(8y y2>.6.因此S奇=0,即奇函数的拟合方差为0,达到最佳逼近.

21、同样，从（2.1.2.5）式中可知偶函数拟合方差为2S2&（Ut X）2= （ a °x 户 2（a -c x） +（ a -4 c2 &）17由（2.1.2.3）式得知在边界点上1 ,、u2 = (w2 W/) = a 4c .2考虑到边界点连续这一约束条件,令(2.1.2.7)(2.3.2.8)因此由式(2.3.2.7)可令S2 =篱 -(a 4c-x2)2 =(a-x。)2 2(a c*)2*-为)2423= (a-X0)2(a4人一含31解令 W=0,有 2(a-x0)+3(a+-e-x1)=0，得fa44a=(12x 8x0 -3e).17从(2.1.2.1

22、0)式可知三次曲线函数的系数a,c的取值与边界点值有关， 式代入(2.1.2.9)式中可得S12 =编-(a 4c-x2)2 =编-(e-x2)2 =1(3x° e-4x)217(2.1.2.9)(2.1.2.10)将(2.1.2.10)e = a 4c所以得出 S2 =(ex2)2+1(3x0+e4x1)2,再令 旦=0,有17；e、2小,、-2(ex2) (3x0 e4x1) = 0,174x1 - 3< 01 x 2e =.(2.1.2.11)解得18联立式(2.1.2.8)、式(2.1.2.10)、式(2.1.2.11),解得 a = (3后 4毛 一 x2) 6 c

23、=(-3-2x1 5x2)18最后得到三次拟合曲线表达式为(3x° 4x * . (y2 -2y)t (-3x° -2x(8y - y2)1863基于两点三次Hermite插值的分段曲线拟台3.1插值的定义定义3.1.19设函数y = f(x)在区间a,b上有定义，且已知在点a _ 尚：为：x2 ：: xn _ b处的函数值yj= f(xj),( j =0,1,2，n),若存在n次多项式pn(x) =a0 +%x +a2x2 + +anxn(3.1.1)使得Pn(xj)=yj,(j =0,1,2，,n)(3.1.2)成立，则称Pn(x)为f (x)的插值多项式,为,为法2，

24、X为插值结点,f ( x)为插值函数.3.2 Hermite插值方法Hermite插值方法可以处理插值条件中合导数值的插值问题，即知道插值结点处的 函数值以及导数值，求插值多项式的插值问题.3.2.1三次Hermite插值考虑两个插值结点的情形,设a'xocxb，函数f(x)在ca,b 且已知f (尚)=y°, f(x)= y，f (均)= m°, f'(x)=m ,在区间a,b上求三次插值函数(3.2.1.1)23H (x) =a° ax a2xa3x使其满足插值条件H(x)=y,H(x3 甲用 0,(3.2.1.2)定理3.2.1.19满足插值

25、条件(3.2.1.2)的三次Hermite插值多项式是存在且唯一的.证明：由插值条件得线性方程组1100x。x112 x°2 X2x02x,a0aya2 lm0八a3 )<m1 J3x°3 x3x23x2(3.2.1.3)(3.2.1.4)考虑系数矩阵行列式，利用行列式的拉普拉斯展开定理，可得故系数矩阵非奇异，线性方程组(3.1.2.3)有唯一解，从而三次多项式存在且唯例1求满足插值条件 "(为)=1,0(00) =0,口(冷)=0,口3) = 0的三次插值多项式 a(x)，以及满足插值条件 E(A)= 0,E xi * 0P' x0习1, X(=)

26、"次插值多项式 -(x).解：由丁为是三次多项式口(x)的二重零点，故可设2一 (x) =(gx C2)(x-x)口(尚)=1,心。)=0得2(Gx0 +0)(x0 x)=1 ,2G(x0-x)2(c1x0 C2)(x0-x)T0由插值条件求解得21G = _,C2 =3(x°xi)(x°xi)(x0x1)2x0代入2-(x) =(Gx C2)(x-x)整理得由插值条件x - x°x 一 x 2:(x) =(1 2 )(一)x - xa x - x°由丁 x1是三次多项式P(x)的二重零点，x0是一重零点，故可设-(x) =c(x x0)(x

27、x1)2，3) =1得c(x° 次)2 2(x° -处)(><0 -*) =1求解得1c= J(x -x。)所以：(x) = (x-x°)d2X -x°注：例题中的两个特殊的插值函数实际上是两点Hermite插值的基函数.定理3.2.1.29两点Hermite插值函数可以用基函数的方法表示为H(x) =)0财*) + 辩1江 +mA(x) +mR(x) , x" x°,为(3.2.1.5)其中1x02x03x0123Xxx012x03x2012x13x2=一(x1 - xo):o(x) =(1 2°)()2x x

28、o xi xo：i(x) =(1 2)(°)2x1 - xo xi - x0r(x) =(x x)( x x0 )2F(x)=(x*(%)2xi - x01:i(xj) = c0(勺=0注：定理3.2.1.2中的cc0(x)1(x)J0(x)J1(x)为Hermite插值基函数，其中'=.")=0；i = ji = j.:j"。i = j例 2 给定 f (1) =0, f (1) = 4, "(1) = 2, f (1) = 0 ,求 Hermite 插值多项式. 解：H3(x) =0h0(x) 4h1(x) 2g°(x) 0g(x)

29、.显然本题不必计算ho(x), g1 (x).(x-1) 丫一1 一1 人1 _ (一1)¥h(x) 112= (2x)(x+1)2/4)i'x-1 f2 /g°(x) =(x-(-1)(x 1)(x-1) 4-1 -1H3(x) =4h(x) 2g°(x)H3(x) =(2x)(x 1)2 1 (x 1)(x1)2 23.3基于Hermite插值的分段曲线拟合基本原理的主要步骤第一步：根据给出的数据做出其散点图，第二步：分析散点图的特点，通过拟合试验确定分段拟合函数，第三步：采用MATLAB编程求得分段拟合函数的表达式，第四步：利用Hermite插值求出

30、分段边界点的插值多项式，第五步：将插值多项式与分段拟合函数连接成连续的拟合曲线4实例应用例3在农业生产试验研究中，对某地区土豆的产量与化肥的关系做了一实验, 得到磷肥的施肥量与土豆产量的对应关系如表 4-1所示.表4-1磷肥的施肥量与土豆产量的对应关系磷施肥(公 斤/公顷)024497398147196245294342土豆产量(公斤)33.4632.4736.0637.9641.0440.0941.2642.1740.3642.73根据上表的数据给出土豆产量与磷肥的关系做出其散点图，如图 4-1所示.土豆产量（公斤）图4-1磷肥的施肥量与土豆产量对应关系的散点图从图可看出从0到98、从98到

31、342之间分别呈明显的线性关系， 由此 可选取所求拟合函数为一分段的线性函数作拟合试验， 换言之，用前5点作一线 性拟合函数，再用后5个点也作一线性拟合函数.采用MATLAB编程（见附录1）求得，对磷肥的分段拟合函数0.0844x 32.07710.0090X 39.1303分段拟合图如图4-2所示.图4-2磷肥的施肥量与土豆产量分段拟合曲线图考虑到边界点不连续，采用两点三次Hermite插值使边界点连续的方法，由丁y(98) =40.3483, y(147) =40.4533, y'(98) = 0.0844, y'(147) = 0.0090，故可以设其 Hermite 插

32、值多项式为H3(x) =40.3483h0(x) 40.4533hi(x) 0.0844g。(x) 0.0090gi(x)40.3483(147 -x)2(2x-147) 40.4533(x - 98)2(343 - 2x) H3(x):经计算得117649-2 .20.0844(147 -x) (x-98) 0.0090(x-98) (x-147)2401即H3（x） =0.0000371155<3 -0.0144093294x2 1.8392571429-36.444499972将插值多项式与分段边界点连接便可以得到连续的拟合曲线图，达到较好的拟合效果.拟合曲线图如图4-3所示（程序

33、见附录2）.图4-3磷肥的施肥量与土豆产量的Hermite插值分段拟合曲线图例4弹簧受力F的作用伸长x, F与x在一定范围内服从虎克定律：F = kx （ x为 弹性系数），呈线性关系；但当F增加到一定值后，不再服从虎克定律.一次试验测得的 数据如表4-2所示，其散点图如图4-4所示.表4-2弹簧受力与伸长量的关系x(cm)135791112141618F(N)1.95.38.612.115.716.819.220.721.421.8(N)图4-4弹簧受力与伸长量的关系的散点图通过散点图先拟合试验，得出前5个点可用线性拟合，后5个点可作二次函数拟合； 同样采用分段拟合的方法，方法同例3 (可设

34、y1 =职+b, y2 = a2+a1x2+a0x22).运行程序(见附录 3)可得 ao =-0.1350向=4.5518;a2 =-16.5508;k=1.7200;b = 0.1200.同样将 拟合函数的边界点采用两点三次 Hermite插值.由 y(9) =15.6000, y(11)=17.1840, y'(9) =1.7200, y(11) = 1.5818，采用 MATLAB 编程 (见附录4)求得插值多项式为 32H3 =0.4294499999999997 x -12.91805000000000 x 129.8885499999999x -420.103949999

35、9994再用插值多项式连接分段拟合曲线的边界点便可得到较好的拟合图形，拟合曲线如图 4-5所示(程序见附录5).图4-5弹簧受力与伸长量的Hermite插值分段拟合曲线图例5在油贞高温分解的过程中，一种苯有机分解成沥宵及其他物质，要了解沥宵在一定温度下随时间t (分钟)变化的相对浓度y (%)之间的关系.试验如表4-3所示,散点图如图4-6所示.表4-3沥宵的相对浓度与时间变化的关系t51520506580100120160180y08.015.120.120.522.020.918.211.55.5图4-6沥宵的相对浓度与时间变化的关系的散点图同样通过散点图先作拟合试验，得出前5个点可采用三

36、次多项式拟合，后5个点可 采用二次多项式拟合，可设分段拟合函数为32.2y1 = a1x a2x a3x1 a4, y2 = b|X2 b2x2 b3运行程序(见附录6)得出a1 =0.0002;a2 = -0.0285;a3 =1.5456;a4 = -7.4491§ = -0.0013；b =0.1747;b3 =16.2750再将拟合函数的边界点采用两点三次 Hermite插值，由y(65) =20.4974, y(80) =22.2531, y (65) =0.3756, y (80) =0.3827采用MATLAB编程(见附录7)计算求得插值多项式为H3 = 0.00232

37、9807407407407c3 - 0.5064964444444440k2 36.68982888888888<-864.2173592592588将插值多项式连接分段拟合曲线的边界点后得到的拟合曲线图，拟合曲线见图4-7(程序见附录8).图4-7沥宵的相对浓度与时间变化的关系的Hermite插值分段拟合图5结束语本文介绍最小二乘多项式曲线拟合的基本原理，在具体介绍线性拟合、多项式拟合 的基本及方法的基础上，给出了分段曲线拟合的方法与步骤.分段曲线拟合是一种常用 的数据处理方法，但是在分段点处往往不能满足连续与光滑，针对这一问题，本文进一 步给出了 Hermite插值的基本原理，并采用

38、两点三次 Hermite插值连接分段曲线，从而 使分段点处满足一阶连续，最后通过三个实例表明该方法的拟合效果较好.另外，本文仅讨论了使用Hermite插值连接分段线性、分段多项式曲线拟合的方法， 对其他种类的曲线未作讨论.事实上，两点三次Hermite插值的方法连接其他种类的拟合 曲线同样适用.参考文献1 聂翔，张瑞林.最小二乘法在曲线拟合中的实现J.陕西工学院学报，2000, 3: 79-82.2 张东林.分段最小二乘曲线拟合J.沈阳大学学报(自然科学版),1994, 2: 80-83.3 文U晓莉，陈春梅.基于最小二乘原理的分段曲线拟合法J.伊犁教育学院学报，2004, 17(3):131

39、-136.4 蔡山，张浩，陈洪辉，等.基于最小二乘法的分段三次曲线拟合方法研究J.科学技术与工程，2007, 7(3): 352-355.5 张浩，任义广，沙基昌.基于分段三次曲线拟合的广州周发案量预测J.计算机仿真，2008, 25(6):257-260.6 Roychowdhury S. Fuzzy curve fitting using least square principlesJ. IEEE International Conference on Systems, Man and Cybemetics, 1998, 4: 4022-4027.7 高伟，姜水生.分段曲线拟合与离散度加

40、权的数据误差处理方法J.中国测试技术，2005, 11:55-56.8 张兴元.分段函数的光滑方法及其在曲线拟合中的应用J.西南民族大学学报(自然科学版),2007, 33(3): 486-490.9 钟尔杰，黄延祝.数值分析(第四版)M.北京：高等教育出版社，2004.10 韩中庚.数学建模方法及其应用M.北京：高等教育出版社，2005.11 文U卫国.MATLAB程序设计与应用(第二版)M.北京：高等教育出版社，2006.致谢经过几个月的努力和忙碌，本次毕业论文即将完成，对为一个本科生的毕业论文， 由丁经验不足，难免有许多地方考虑不全面，如果没有指导老师的督促与辛勤的指导， 以及一起学习的

41、同学们的帮助与支持，想顺利的完成这篇论文比较难 .值此论文完成之际，首先对指导老师李军成老师表示最诚挚的感谢与崇高的敬意 李老师严谨的治学态度，深厚渊博的学术素养，敏锐的思维，积极进去的精神，严以律 己，宽以待人的崇高品质，乐观向上的人生态度，谦逊和蔼的为人品德，平等的师生关 系，尤其是认真负责的工作态度均给我留下了不可磨灭的印象，相信对我今后的学习、 工作以及生活都会有着深远的影响.感谢陈国华主任、杨笃庆书记、谭本远主任等数学系领导们，你们认真负责的治学 态度和高速度、高效率的办事方式深深的感染了我们， 让我们能够时时刻刻提醒自己要 认真负责对待每件事情、每一个环节，感谢梁经珑老师、杨涤尘老

42、师、余星老师、李军 成老师、邓华老师、钟月娥老师、孙红果老师、李兵老师、龙承星老师等数学系的老师 们，你们的授课方式与渊博的知识深化了我们的知识面，拓广了我们的视野，使我们对数学有了更浓厚的兴趣与体会.感谢杜鹃老师、郑丽峰老师，你们热忱的帮助使我们有 一个很好的学习余围来完成论文.在本文的写作过程中，李军成老师，石小芳、彭迪、方其斌等同学提出了许多宝贵 的意见，此论文的完成离不开他们的指导，特别是李军成老师；他们渊博的学识与敏锐的头脑让我受益匪浅.再次对帮助我的人表示衷心的感谢.附录1磷肥的施肥量与土豆产量的分段拟合函数程序x=0,24,49,73,98,147,196,245,294,342

43、;y=33.46,32.47,36.06,37.96,41.04,40.09,41.26,42.17,40.36,42.73;plot(x,y,'r+')a1=polyfit(x(1:5),y(1:5),1)a2=polyfit(x(6:10),y(6:10),1)xx1=0:98; yy1=a1 (1) *xx1+a1(2);xx2=147:342; yy2=a2(1)*xx2+a2(2);plot(x,y,'r+',xx1,yy1,xx2,yy2)附录2磷肥的施肥量与土豆产量的Hermite插值的分段拟合曲线图程序x=0,24,49,73,98,147,19

44、6,245,294,342;y=33.46,32.47,36.06,37.96,41.04,40.09,41.26,42.17,40.36,42.73;plot(x,y,'r+')a1=polyfit(x(1:5),y(1:5),1)a2=polyfit(x(6:10),y(6:10),1)xx1=0:0.01:98; yy1=a1 (1) *xx1+a1(2);xx2=147:0.01:342; yy2=a2 (1) *xx2+a2 (2);xx3=98:0.01:147;yy3=0.0000371155*xx3.A3-0.0144093294*xx3.A2+1.839257

45、1429*xx3-36.444499972;plot(x,y,'r+',xx1,yy1,xx2,yy2,xx3,yy3)附录3弹簧受力与伸长量的关系的分段拟合函数程序x=1,3,5,7,9,11,12,14,16,18;y=1.9,5.3,8.6,12.1,15.7,16.8,19.2,20.7,21.4,21.8;plot(x,y,'r+')a1=polyfit(x(1:5),y(1:5),1)a2=polyfit(x(6:10),y(6:10),2)xx1=1:0.01:9; yy1=a1 (1) *xx1+a1(2);xx2=11:0.01:18; yy2

46、=a2(1)*xx2.A2+a 2(2) *xx2+a2(3);plot(x,y,'r+',xx1,yy1,xx2,yy2)附录4弹簧受力与伸长量的关系的两点三次Hermite插值多项式程序format long eclf,clear,x0=9;x1=11;y0=15.6000;y1=17.1840;m0=1.7200;m1=1.5818;x=linspace(9,11,100);y=linspace(15.6000,17.1840,100);m=2*(y0-y1)+(m0+m1)*(x1-x0);n=3*(x0+x1)*y1-3*(x0+x1)*y0-(2*x1+x0)*m0

47、*(x1-x0)-(2*x0+x1)*(x1-x0)*m1;k=6*x0*x1*(y0-y1)+(x1-x0)*(m0*x1.A2+m1*x0.A2)+2*x1*x0*(m0+m1);q=x1.A2*(x1-3*x0)*y0+x0.A2*(3*x1-x0)*y1-x0*x1*(x1-x0)*(x1*m0+x0*m1);p=(x1-x0).A3;a=m/pb=n/pc=k/pd=q/py=a*x.A3+b*x.A2+c*x+d;plot(x,y,'r-')附录5弹簧受力与伸长量的Hermite插值的分段拟合曲线图程序x=1,3,5,7,9,11,12,14,16,18;y=1.9,5.3,8.6,12.1,15.7,16.8,19.2,20.7,21.4,21.8;plot(x,y,'r+')a1=p