曲线拟合问题讲解

1、曲线拟合问题摘要本文首先对给定数据根据不同要求进行多次直线拟合，分别求得使所拟直线预期值的偏差平方和、绝对偏差总和和最大偏差最小的三类拟合直线，然后再求得二次曲线条件下满足三类要求的二次拟合曲线，最后运用其他曲线对给定数据进行拟合，得到吻合度最高的曲线。针对问题一，构建线性回归方程，运用最小二乘法及lingo软件使得目标函数预期值的即拟合偏差平方和达到最小，从而得到拟合曲线y二0.8031048x-0.0123077。i针对问题二，构建给定数据的线性回归方程，使得目标函数即预期值的绝对偏差综合最小，但由于绝对偏差较难处理，采用转化的思想将对绝对偏差的求解转化为对偏差平方和开方的求解，从而得到拟

2、合曲线y二0.65x+0.575。i针对问题三，构建给定数据的线性回归方程，运用lingo软件使得目标函数即预期值的最大偏差最小，从而得到拟合曲线y二1.13x-1.879。i针对问题四，构建给定数据的二次方程，运用lingo软件分别求得三类不同条件下的最优拟合曲线，偏差平方和达到最小：y=0.09703011x2-0.138534x+1.425301，绝1ii对偏差总和达到最小：y二0.04148148x2+0.2711111x+1，观测值与预测值最大偏差为1ii最小：y二0.02556818x2+0.7659091x-0.6923295。1ii针对问题五，本文做出给定数据散点图，构建不同曲

3、线类型进行拟合，得到R2即吻合度最高的曲线类型，运用Matlab软件求得该曲线类型的方程。本文的特色在于利用图标直观表达拟合曲线，增强文章可靠性及真实性，并构建不同的曲线类型，得到吻合度最高的拟合曲线。关键词：曲线拟合、线性回归、lingo1问题的重述已知一个量y依赖于另一个量x，现收集有数据如下：x0.00.51.01.51.92.53.03.54.04.5y1.00.90.71.52.02.43.22.02.73.5x5.05.56.06.67.67.68.59.010.0y1.04.07.62.75.74.66.06.812.3（1）求拟合以上数据的直线y=bx+a。目标为使y的各个观察

4、值同按直线关系所预期的值的偏差平方和为最小。（2）求拟合以上数据的直线y=bx+a，目标为使y的各个观察值同按直线关系所预期的值的绝对偏差总和为最小。（3）求拟合以上数据的直线，目标为使y的各个观察值同按直线关系所预期的值的最大偏差为最小。（4）求拟合以上数据的曲线y=cx2+bx+a，实现（1）（2）（3）三种目标。（5）试一试其它的曲线，可否找出最好的？2模型假设（1）假设数据可靠，并且数据量足够充分来反映y与x的依赖关系;2）假设原始数据精度较大，拟合求得的曲线足够精确；3.符号说明y预测值y观测值观测值与预测值的差值4问题一的分析与求解4.1 问题的分析对于给定的点(x,y)(i=1,

5、2,n)，假定x和y之间满足线性模型y=bx+a+b,iiiiN(Oq2),covQq)=0,(i丰j;i,j=1,2,n)，据此，建立线性回归方程$=bx+a,i使得拟合偏差平方和S=2=工y-y|y-bx-a|2达到最小，即利用最小二乘iiiii=1i=1i=1法求解。4.2 问题的求解将数据代入lingo软件(程序代码见附录1),约束条件为:可解得：minS二工y-y二工|y-bx-a|2iiiii=1i=1Ay=0.8031048x0.0123077,S=48.55993imin偏差平方和最小时的函数与数据散点图y线性5问题二的分析与求解5.1 问题的分析对于给定的点(x,y)(i=1

6、,2,n)，假定x和y之间满足线性模型y=bx+a+，iiiiN(0q2),cov(qq)=0，(i丰j;i,j=1,2,n)，据此，建立线性回归方程y=bx+a，i使得拟合绝对偏差总和S二i=1i=1以计算，因此，可将其看作S=工22=工i=1i=1=Y|y-bx-a达到最小。但是，绝对值难ii2厂i=1y-bx-aii55.2 问题的求解将数据代入lingo软件(程序代码见附录2)，约束条件为:S丄i=1i=1y-bx-aii可解得：Ay=0.65x+0.575，S=18.455imin观测值与预测值最大偏差为最小时的函数与数据散点图yy=0.65x+0.575线性(y=0.65x+0.5

7、75)6问题三的分析与求解6.1 问题的分析对于给定的点(x,y)(i=1,2,n)，假定x和y之间满足线性模型y=bx+a+b，iiiiN(0q2),cov(qq)=0，(i丰j;i,j=1,2,n)，据此，建立线性回归方程$=bx+a,i使y的各个观察值同按直线关系所预期的值的最大偏差为最小，即求各线性方程中maxS=y-y=|y-bx-a达到最小时，对应的线性方程。iii6.2 问题的求解将数据代入lingo软件(程序代码见附录3)，约束条件为：min(maxS)二min(maxy-y)二min(max|y-bx-a|)iii可解得：Ay二1.13x-1.879，maxS=2.879im

8、iny=1.13x-1.879线性(y=1.13x-1.879)7问题四的分析与求解7.1 问题的分析该问题分析方法与上述三个问题相同，对于给定的点(x,y)(i=1,2,n)，假定X和iiy之间满足模型y=cx2+bx+a+b，cN(0,c2),cov©,c)=0，(i丰j;i,j=1,2,n)，iii据此，建立方程y=cx2+bx+a，分别求iS=EC2=2Li=1i=LS=i=1i=1Ay-yiiA-yii=工卜-bx-a|iii=1=Qy一bx-a|iii=1使得=|y一bx一a|ii达到最小。7.2 问题的求解将数据代入lingo软件(程序代码见附录4),约束条件分别对应为

9、:7minS=工=Qy-bx-a|iiS=工i=1y-yii2i=1Ay-bx-aiimin(maxS)=min(maxy一y)=min(max|y一bx一a|)ii可解得：偏差平方和达到最小：y=0.09703011x2-0.138534x+1.425301,S=36.934921iimin偏差平方和最小时的函数与数据散点图xy=0.09703011x辽-0.1385340x+1.425301多项式(y=0.09703011x辽-0.1385340x+1.425301)图49观测值与预测值最大偏差为最小：y=a-ebx+c+d，maxS=2.77642min观测值与预测值观测值与预测值最大偏

10、差为最小时的函数与数据散多项式(y二0.02556818x2+0.765y=0.02556818X2+0.7659091X-0.69232958问题五的分析与求解8.1 问题的分析该问题分析方法主要是采用最小二乘法拟合函数图像，利用观测值与预测值的平方和达到最小，即S=2二工y-y二工卜-bx-a|2达到最小。米用最小二乘法进行iiiii=1i=1i=1参数估计时，R2表示回归平方和与总离差平方和的比值，这一比例越大越好，模型越精确，回归效果越显著。R2介于01之间，越接近1,回归拟合效果越好，一般认为超过0.8的模型拟合度比较高。8.2 问题的求解将原始数据录入Matlab拟合工具箱，观察其

11、散点图分布情况，发现指数函数、三次多项式和三角函数能够较好的反映散点分布，分别得拟合图像（图1、图2、图3）。11y二1.088eo.3728x-i.47i+1.0281y=0.02869X3-0.327x2+1.485x+0.2639x图2y二-0.1515sinx+0.7992x+0.02906图3通过对R值的比较R=0.7927,R=0.8088,R=0.688发现R相较于其他来说最1232大，所以采用三次多项式作为目标函数。9模型的优缺点分析与改进方向9.1 优点：（1）本文首先画出所给数据散点图，然后运用图表直观表达拟合曲线，避免了只有函数方程，增强了本文的可读性和理解性。（2）本文

12、在构建绝对偏差最小目标函数时，由于绝对偏差较难处理，采用转化的思想，将对绝对偏差的求解转化为对偏差平方和开方的求解，从而巧妙地解决这一难题，使得文章更加通俗易懂。（3）本文所用程序增加注释，增强了程序的可读性和规范性。9.2 缺点：本文虽建立不同拟合曲线类型，但也仅限于有限的曲线类型，无法保证此曲线即是吻合度最高的曲线。9.3 改进方向：针对模型特点，在此基础上，找到吻合度更高的拟合曲线类型，从而得到最优的拟合曲线。参考文献1 周鹏，许钢.基于LabVIEW的广义线性拟合在成本预测中的应用J,安徽工程大学学报，2013，(3)2 王礼想，刘利姣，黄光明.基于EM算法的线性拟合问题研究J.廊坊师

13、范学院学报,自然科学版,2013,(4)3 狄晓敏，谢红薇.多疾病共同危险因素挖掘与MARS预测模型研究J.电子学报，2009，(6)4 郝海燕，郝春蕾，康荣雷等.我国普通高等学校生源规模预测及高职院校发展趋势分析.承德石油高等专科学校学报，2013，(4)17附录：1.第一问：model:sets:!定义集math,属性x,y,共1一19个成员数据;math/1.19/:x,y;endsetsdata:x=00.511.51.92.533.544.555.566.67.67.68.5910y=10.90.71.522.43.222.73.5147.62.75.74.666.812.3Endd

14、ata!约束条件;min=sum(math:(a*x+b-y)入2);!解除变量a,b的非负限制；free(a);free(b);end2.第二问：model:sets:!定义集math，属性x,y,共1一19个成员数据;math/1.19/:x,y;endsetsdata:x=00.511.51.92.533.544.555.566.67.67.68.5910y=10.90.71.522.43.222.73.5147.62.75.74.666.812.3Enddata!约束条件;min=sum(math(i):(y(i)-z(i)八2)八(1/2);for(math(i):z(i)=a*x(

15、i)+b);!解除变量a,b的非负限制；free(a);free(b);end3.第三问model:sets:!定义集math,属性x,y,共1一19个成员数据；math/1.19/:x,y;endsetsdata:x=00.511.51.92.533.544.555.566.67.67.68.5910y=10.90.71.522.43.222.73.5147.62.75.74.666.812.3Enddata!约束条件;min=max(math:abs(a*x+b-y);!解除变量a,b的非负限制；free(a);free(b);End4.第四问4.1model:sets:!定义集math,

16、属性x,y,共1一19个成员数据;math/1.19/:x,y;endsetsdata:x=00.511.51.92.533.544.555.566.67.67.68.5910y=10.90.71.522.43.222.73.5147.62.75.74.666.812.3Enddata!约束条件;min=sum(math:(a*x2+b*x+c-y)入2);!解除变量a,b,c的非负限制；free(a);free(b);free(c);End4.2model:sets:!定义集math,属性x,y,共1一19个成员数据;math/1.19/:x,y;endsetsdata:x=00.511.51.92.533.544.555.566.67.67.68.5910y=10.90.71.522.43.222.73.5147.62.75.74.666.812.3Enddata!约束条件;min=sum(math:abs(a*x2+b*x+c-y);!解除变量a,b,c的非负限制；free(a);free(