数学建模曲线拟合

《数学建模曲线拟合》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、曲线拟合摘要根究已有数据研究y关于x的关系，对于不同的要求得到不同的结果。问题一中目标为使的各个观察值同按直线关系所预期的值的偏差平方和为最小，利用MATLAB中lsqcurvefit函数在最小二乘法原理下拟合出所求直线。问题二目标为使绝对偏差总和为最小，使用MATLAB屮的fminsearch函数，在题冃约束条件内求的最优答案，以此方法同样求得问题三中最大偏差为最小时的直线。问题四拟合的曲线为二阶多项式，方法同前三问类似。问题五为求得最佳的曲线，将之前的一次曲线换成多次曲线进行拟合得到新的结果。经试验发现高阶多项式的阶数越高拟和效果最

2、好。关键词：函数拟合最小二乘法线性规划、问题的重述已知一个量y依赖于另一个量兀，现收集有数据如下: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.07.6&59.010.0y1.04.07.62.75.74.66.06.812.3(1)求拟合以上数据的直线y=bx-^ao目标为使y的各个观察值同按直线关系所预期的值的偏差平方和为最小。(2)求拟合以上数据的直线y=加+_日标为使y的各个观察值同按直线关系所预期的值的绝对偏差总和为最小

3、。(3)求拟合以上数据的直线，目标为使y的各个观察值同按直线关系所预期的值的最大偏差为最小。(4)求拟合以上数据的曲线y二cF+加+Q,实现(1)(2)(3)三种目标。(5)试一试其它的曲线，可否找出最好的？二、问题的分析对于问题一，利用MATLAB中的最小二乘法对数据进行拟合得到直线，目标为使各个观察值同按直线关系所预期的值的偏差平方和为最小。对于问题二、三、四均利用MATLAB中的fminscarch函数，在题日要求的约束条件下找到最佳答案。对于问题五，改变多项式最高次次数，拟合后计算残差，和二次多项式比较,再增加次数后拟合，和原多

4、项式比较残差，进而找到最好的曲线。三、基本假设1.表中数据真实可信，每个点都具有意义。四、模型的建立与求解1.问题一对给定数据点{(Xj,g)}(i=O,l,…，肋，在取定的函数类①中，求,使误差的平方和F最小，Q=》b(xJ_疔。从几何意义上讲，就是寻求与给定点{(X/)}(心0,1,…冲)的距离平方和为最小的曲线y=p{x)o函数〃(兀)称为拟合函数或最小二乘解，求拟合函数“(0的方法称为曲线拟合的最小二乘法。直接利用MATLAB屮lsqcuwefit的函数进行曲线拟合，得到目标函数如下：y=-0.0264+0.8117x函数图像如

5、下:1.问题二利用MATLAB中的fminsearch函数，在题冃要求的约束条件使y的各个观察值同按直线关系所预期的值的绝对偏差总和为最小下进行曲线拟合，得到目标函数如下：y=0.5001+0.6666•x函数图像如下:3•问题三利用MATLAB中的fminsearch函数，在题目要求的约束条件使y的各个观察值同按肓线关系所预期的值的最大偏差为最小下进行曲线拟合，得到目标函数如下：>=1.8790+1.1300讥函数图像如下:4.问题四(1)问题一同问题一相似，只是拟合的曲线为二阶多项式，利用MATLAB中Isqcurvefit的函数进

6、行曲线拟合，得到目标函数如下：y=1.3833・0.1096•兀+0.0963•x2函数图像如下：(2)问题二同问题二相似，只是拟合的曲线为二阶多项式，利用MATLAB中的fminsearch函数，在题目要求的约束条件使y的各个观察值同按直线关系所预期的值的最大偏差为最小下进行曲线拟合，得到目标函数如下：y=0.9755+0.2902•x+0.0397-x2函数图像如下：(2)问题三同问题三求解过程相似，只是拟合的曲线为二阶多项式，利用MATLAB中的fminsearch函数，在题目要求的约束条件使y的各个观察值同按直线关系所预期的值的

7、最大偏差为最小下进行曲线拟合，得到目标函数如下:y=1.7672・0.0909•x+0.0994•x2函数图像如下:5•问题五选择更高阶多项式进行曲线拟合，利用MATLAB中lsqcurvefit的函数进行曲线拟合，得到目标函数。比较方差，方差越小，得到结果越稳定，即认为曲线拟合越好，方差结果如下表所示：阶数12345方差2.68842.05771.71271.50491.4336可以得到函数阶数越高，曲线拟合越好。各阶多项式函数图像如下：三阶：四阶:五阶:五、模型的评价对于问题五屮的模型，由于我们只选择了最高为五阶的高阶多项式多项式进

8、行曲线拟合，述需要选择更多的函数进行拟合，并进行检验，找到最好的曲线。六、附录%1funl=@(a,x)a(1).*x+a(2);a=lsqcurvefit(funl,[0,0],x,y)xi=0:0.1: