MATLAB曲线拟合的应用.doc

《MATLAB曲线拟合的应用.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、MATLAB曲线拟合的应用王磊品吴东新疆泒犨泰克石油科技有限公司新疆油田公司准东采油厂信息所摘要：1.阐述MATLAB数学分析软件的基本功能；2.对MATLAB在生产数据分析中的应用进行了研究,指出曲线拟合的基本方法；3.以实例阐明MATLAB与行业生产数据结合对生产数据进行分析的原理。关键词：MATLAB；曲线拟合；插值1.引言在生产开发过程中,复杂的生产数据之间或多或少的存在着这样或者那样的联系,如何利用现今普及的计算机以及网络资源在最短的时间内找到这个联系,以指导我们的生产开发,这对于行业科研人员来说无疑是一个最为关心的问题。MA

2、TLAB矩阵分析软件,自推出以来，已成为国际公认的最优秀的数学软件之一，其范围涵盖了工业、电子、医疗以及建筑等各个领域，以其强大的科学计算功能使众多科研机构纷纷采用。为此，本文从介绍MATLAB软件开始，以实例讲述如何使用MATLAB对生产开发数据进行计算与分析，从而达到高效、科学指导生产的目的。2.MATLAB简介MATLAB是MathWorks公司于1982年推出的一套高性能的数值计算和可视化数学软件。由于使用编程运算与人进行科学计算的思路和表达方式完全一致，所以不象学习其它高级语言那样难于掌握，用Matlab编写程序犹如在演算纸上

3、排列出公式与求解问题，所以又被称为演算纸式科学算法语言。在这个环境下，对所要求解的问题，用户只需简单地列出数学表达式，其结果便以数值或图形方式显示出来。MATLAB的含义是矩阵实验室（MATRIXLABORATORY），主要用于方便矩阵的存取，其基本元素是无须定义维数的矩阵。自问世以来,就是以数值计算称雄。MATLAB进行数值计算的基本单位是复数数组（或称阵列），这使得MATLAB高度“向量化”。经过十几年的完善和扩充，现已发展成为线性代数课程的标准工具。由于它不需定义数组的维数，并给出矩阵函数、特殊矩阵专门的库函数，使之在求解诸如信号

4、处理、建模、系统识别、控制、优化等领域的问题时，显得大为简捷、高效、方便，这是其它高级语言所不能比拟的。美国许多大学的实验室都安装有供学习和研究之用。MATLAB中包括了被称作工具箱（TOOLBOX）的各类应用问题的求解工具。工具箱实际上是对MATLAB进行扩展应用的一系列MATLAB函数（称为M文件），它可用来求解各类学科的问题，包括信号处理、图象处理、控制系统辨识、神经网络等。随着MATLAB版本的不断升级，其所含的工具箱的功能也越来越丰富，因此，应用范围也越来越广泛，成为涉及数值分析的各类工程师不可不用的工具。MATLAB中包括了

5、图形界面编辑GUI，改变了以前单一的“在指令窗通过文本形的指令进行各种操作的状况。这可让使用者也可以像VB、VC、VJ、DELPHI等那样进行一般的可视化的程序编辑。在命令窗口（matlabcommandwindow）键入simulink，就出现(SIMULINK)窗口。以往十分困难的系统仿真问题，用SIMULINK只需拖动鼠标即可轻而易举地解决问题，这也是近来受到重视原因所在。1.曲线拟合与插值在大量的应用领域中，人们经常面临用一个解析函数描述数据(通常是测量值)的任务。对这个问题有两种方法。在插值法里，数据假定是正确的，要求以某种方

6、法描述数据点之间所发生的情况。这种方法在下一节讨论。这里讨论的方法是曲线拟合或回归。人们设法找出某条光滑曲线，它最佳地拟合数据，但不必要经过任何数据点。1.1.曲线拟合曲线拟合涉及回答两个基本问题：最佳拟合意味着什么？应该用什么样的曲线？当最佳拟合被解释为在数据点的最小误差平方和，且所用的曲线限定为多项式时，那么曲线拟合是相当简捷的。数学上，称为多项式的最小二乘曲线拟合。对各数据点距离求平方，并把平方距离全加起来，就是误差平方和。这条虚线是使误差平方和尽可能小的曲线，即是最佳拟合。最小二乘这个术语仅仅是使误差平方和最小的省略说法。在MA

7、TLAB中，函数polyfit求解最小二乘曲线拟合问题。»x=[0.1.2.3.4.5.6.7.8.91]；»y=[-.4471.9783.286.167.087.347.669.569.489.3011.2]；为了用polyfit，我们必须给函数赋予上面的数据和我们希望最佳拟合数据的多项式的阶次或度。如果我们选择n=1作为阶次，得到最简单的线性近似。通常称为线性回归。相反，如果我们选择n=2作为阶次，得到一个2阶多项式。现在，我们选择一个2阶多项式。»n=2;%给定阶数»p=polyfit(x,y,n)%输出多项式向量polyfit输

8、出一个多项式系数的行向量。其解是y=-9.8108x2+20.1293x-0.0317。为了将曲线拟合解与数据点比较，让我们把二者都绘成图。»xi=linspace(0,1,100);%生成100个点的行向