曲线拟合问题_第1页
曲线拟合问题_第2页
曲线拟合问题_第3页
曲线拟合问题_第4页
曲线拟合问题_第5页
已阅读5页,还剩9页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

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

2、程,运用lingo软件使得目标函数即预期值的最大偏差最小,从而得到拟合曲线。针对问题四,构建给定数据的二次方程,运用lingo软件分别求得三类不同条件下的最优拟合曲线,偏差平方和达到最小:,绝对偏差总和达到最小:,观测值与预测值最大偏差为最小:。针对问题五,本文做出给定数据散点图,构建不同曲线类型进行拟合,得到即吻合度最高的曲线类型,运用Matlab软件求得该曲线类型的方程。本文的特色在于利用图标直观表达拟合曲线,增强文章可靠性及真实性,并构建不同的曲线类型,得到吻合度最高的拟合曲线。关键词:曲线拟合、线性回归、lingo 1问题的重述已知一个量依赖于另一个量,现收集有数据如下:0.00.51

3、.01.51.92.53.03.54.04.51.00.90.71.52.02.43.22.02.73.55.05.56.06.67.67.68.59.010.01.04.07.62.75.74.66.06.812.3(1)求拟合以上数据的直线。目标为使的各个观察值同按直线关系所预期的值的偏差平方和为最小。(2)求拟合以上数据的直线,目标为使的各个观察值同按直线关系所预期的值的绝对偏差总和为最小。(3)求拟合以上数据的直线,目标为使的各个观察值同按直线关系所预期的值的最大偏差为最小。(4)求拟合以上数据的曲线,实现(1)(2)(3)三种目标。(5)试一试其它的曲线,可否找出最好的?2模型假设(

4、1)假设数据可靠,并且数据量足够充分来反映与的依赖关系;(2)假设原始数据精度较大,拟合求得的曲线足够精确;3.符号说明预测值观测值观测值与预测值的差值4问题一的分析与求解4.1 问题的分析对于给定的点,假定和之间满足线性模型,,,据此,建立线性回归方程,使得拟合偏差平方和达到最小,即利用最小二乘法求解。4.2 问题的求解将数据代入lingo软件(程序代码见附录1),约束条件为:可解得:,图 15问题二的分析与求解5.1 问题的分析对于给定的点,假定和之间满足线性模型,,,据此,建立线性回归方程,使得拟合绝对偏差总和达到最小。但是,绝对值难以计算,因此,可将其看作。5.2 问题的求解将数据代入

5、lingo软件(程序代码见附录2),约束条件为:可解得:, 图 2 6问题三的分析与求解6.1 问题的分析对于给定的点,假定和之间满足线性模型,,,据此,建立线性回归方程,使的各个观察值同按直线关系所预期的值的最大偏差为最小,即求各线性方程中达到最小时,对应的线性方程。6.2 问题的求解将数据代入lingo软件(程序代码见附录3),约束条件为:可解得:, 图 37问题四的分析与求解7.1 问题的分析该问题分析方法与上述三个问题相同,对于给定的点,假定和之间满足模型,,,据此,建立方程,分别求使得达到最小。7.2 问题的求解将数据代入lingo软件(程序代码见附录4),约束条件分别对应为:可解得

6、:偏差平方和达到最小:,图 4绝对偏差总和达到最小:,图 5观测值与预测值最大偏差为最小:,图 68问题五的分析与求解8.1 问题的分析该问题分析方法主要是采用最小二乘法拟合函数图像,利用观测值与预测值的平方和达到最小,即达到最小。采用最小二乘法进行参数估计时,表示回归平方和与总离差平方和的比值,这一比例越大越好,模型越精确,回归效果越显著。介于01之间,越接近1,回归拟合效果越好,一般认为超过0.8的模型拟合度比较高。8.2 问题的求解将原始数据录入Matlab拟合工具箱,观察其散点图分布情况,发现指数函数、三次多项式和三角函数能够较好的反映散点分布,分别得拟合图像(图 1、图 2、图 3)

7、。图 1图 2图 3通过对值的比较,发现相较于其他来说最大,所以采用三次多项式作为目标函数。9模型的优缺点分析与改进方向9.1优点:(1)本文首先画出所给数据散点图,然后运用图表直观表达拟合曲线,避免了只有函数方程,增强了本文的可读性和理解性。(2)本文在构建绝对偏差最小目标函数时,由于绝对偏差较难处理,采用转化的思想,将对绝对偏差的求解转化为对偏差平方和开方的求解,从而巧妙地解决这一难题,使得文章更加通俗易懂。 (3)本文所用程序增加注释,增强了程序的可读性和规范性。9.2缺点:本文虽建立不同拟合曲线类型,但也仅限于有限的曲线类型,无法保证此曲线即是吻合度最高的曲线。9.3改进方向:针对模型

8、特点,在此基础上,找到吻合度更高的拟合曲线类型,从而得到最优的拟合曲线。参考文献1 周鹏,许钢. 基于LabVIEW的广义线性拟合在成本预测中的应用J,安徽工程大学学报,2013,(3)2 王礼想,刘利姣,黄光明. 基于EM算法的线性拟合问题研究J. 廊坊师范学院学报,自然科学版,2013,(4)3 狄晓敏,谢红薇. 多疾病共同危险因素挖掘与MARS预测模型研究J.电子学报,2009,(6)4 郝海燕,郝春蕾,康荣雷等. 我国普通高等学校生源规模预测及高职院校发展趋势分析. 承德石油高等专科学校学报,2013,(4)附录:1.第一问:model:sets:!定义集math,属性x,y,共119

9、个成员数据;math/1.19/:x,y;endsetsdata:x=00.511.51.92.533.544.555.566.67.67.68.5910;y=10.90.71.522.43.222.73.5147.62.75.74.666.812.3;Enddata!约束条件;min=sum(math:(a*x+b-y)2);!解除变量a,b的非负限制;free(a);free(b);end2.第二问:model:sets:!定义集math,属性x,y,共119个成员数据;math/1.19/:x,y;endsetsdata:x=00.511.51.92.533.544.555.566.67

10、.67.68.5910;y=10.90.71.522.43.222.73.5147.62.75.74.666.812.3;Enddata!约束条件;min=sum(math(i):(y(i)-z(i)2)(1/2);for(math(i):z(i)=a*x(i)+b);!解除变量a,b的非负限制;free(a);free(b);end3.第三问model:sets:!定义集math,属性x,y,共119个成员数据;math/1.19/:x,y;endsetsdata:x=00.511.51.92.533.544.555.566.67.67.68.5910;y=10.90.71.522.43.2

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

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

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论