实验数据曲线拟合方法研究

本科毕业设计论文题 目 实验数据曲线拟合方法研究 专业名称学生姓名指导教师毕业时间毕业 任务书一、题目实验数据曲线拟合方法研究二、指导思想和目的要求通过毕业设计，使学生对所学自动控制原理、现代控制原理、控制系统仿真、电子技术等的基本理论和基本知识加深理解和应用；培养学生设计计算、数据处理、文件编辑、文字表达、文献查阅、计算机应用、工具书使用等基本事件能力以及外文资料的阅读和翻译技能；掌握常用的实验数据曲线拟合方法，培养创新意识，增强动手能力，为今后的工作打下一定的理论和实践基础。要求认真复习有关基础理论和技术知识，认真对待每一个设计环节，全身心投入，认真查阅资料，仔细分析被控对象的工作原理、特性和控制要求，按计划完成毕业设计各阶段的任务，重视理论联系实际，写好毕业论文。三、主要技术指标设计系统满足以下要求：数据拟合误差要尽量的小的同时保证曲线的线形形状最佳。四、进度和要求1、搜集中、英文资料，完成相关英文文献的翻译工作，明确本课题的国内外研究现状及研究意义； （第 1、2 周）2、撰写开题报告； （第 3、4 周）3、应用最小二乘法进行曲线拟合； （第 5、6 周）4、应用 Matlab 命令曲线拟合； （第 7、8周）5、应用 Matlab 图形用户界面曲线拟合； （第 9、10设计论文周）6、研究其他曲线拟合方法； （第 11 周）7、整理资料撰写毕业论文；（1）初稿 ； （第 12、13 周）（2）二稿 ； （第 14 周）8、准备答辩和答辩。 （第 15 周）五、主要参考书及参考资料1卢京潮， 自动控制原理 ，西北工业大学出版社，2010.62胡寿松， 自动控制原理 ，科学出版社，2008，63薛定宇，陈阳泉， 系统仿真技术与应用 ，清华大学出版社，2004.44王正林， Matlab/Simulink 与控制系统仿真 ，电子工业出版社，2009.75李桂成， 计算方法 ，电子工业出版社，2013.86蒋建飞，胡良剑，唐俭.数值分析及其 Matlab 实验【M】.北京：科学出版社，2008西北工业大学明德学院本科毕业设计论文学生 指导教师 系主任 西北工业大学明德学院本科毕业设计论文I摘 要在我们实际的实验和勘探中，都会产生大量的数据。为了解释这些数据或者根据这些数据作出预测、判断，给决策者提供重要的依据。需要对测量数据进行拟合，寻找一个反映数据变化规律的函数。本文介绍了几种常用的数据拟合方法，线性拟合、二次函数拟合、数据的n 次多项式拟合等。并着重对曲线拟合进行了研究，介绍了线性与非线性模型的曲线拟合方法，最小二乘法、牛顿迭代法等。在传统的曲线拟合基础上，为了提高曲线拟合精度，本文还研究了多项式的摆动问题，从实践的角度分析了产生这些摆动及偏差的因素和特点，总结了在实践中减小这些偏差的处理方法。采用最小二乘法使变量转换后所得新变量离均差平方和最小，并不一定能使原响应变量的离均差平方和最小，所以其模型的拟合精度仍有提高的空间。随着计算机技术的发展，实验数据处理越来越方便。但也提出了新的课题，就是在选择数据处理方法时应该比以往更为慎重。因为稍有不慎，就会非常方便地根据正确的实验数据得出不确切的乃至错误的结论。所以提高拟合的准确度是非常有必要的。关键词：数据拟合、最小二乘法、曲线拟合、多项式摆动西北工业大学明德学院本科毕业设计论文IIABSTRACTIn our experiments and exploration, it will produce large amounts of data. In order to explain these data to make predictions based on these data to determine, provide an important basis for policy makers . Need to fit the measured data to find a function to reflect data changes in the law. This article describes several commonly used data fitting methods, and focused on a nonlinear curve fitting of the model.This paper introduces some commonly used data fitting method, linear fitting, secondary function fitting, data n times polynomial fitting etc. T And focuses on the curve fitting, introduced the linear and nonlinear model of curve fitting method, the least square method, Newton iterative method, etc. In the traditional curve fitting basis, in order to improve the curve fitting precision, this paper also studies the polynomial swing, from the perspective of the practice the oscillation and deviation of factors and characteristics, and summarizes the decrease in practice the treatment method of these deviations. The least square method to variable after converting from new variables are the sum of squared residuals minimum, not necessarily make the original response from all the variables of the sum of squared residuals minimum, so the model fitting precision still has room to improve. With the development of computer technology, the experiment data processing more and more convenient. But also put forward the new subject, which is in the data processing method of choice should be more careful than ever before. Because carelessly a bit, it can be very easily according to the correct experimental data that not the exact and even the wrong conclusion. Therefore, to raise the fitting accuracy is very necessary.KEY WORDS： Data Fitting , Least square method , Curve fitting , Polynomial swing西北工业大学明德学院本科毕业设计论文III目录摘 要 .IABSTRACT .II第 1 章 绪论 .51.1 引言 .51.2 研究背景 .51.3 研究意义 .71.4 本论文主要内容 .7第 2 章 曲线拟合及最小二乘法 .92.1 线性模型的曲线拟合 .92.2 最小二乘法基本原理 .92.3 用正交多项式作最小二乘拟合 .112.4 非线性模型的曲线拟合 .132.4.1 常见非线性模型 .132.4.2 牛顿迭代 .15第 3 章 基于 MALTAB 实现最小二乘法 .183.1 Matlab 简介 .183.1.1 Matlab 的概况 .183.1.2 Matlab 的语言特点 .183.1.3 Matlab 工作界面 .193.1.4 优势特点 .203.2 用 MALTAB 实现曲线拟合 .203.2.1 最小二乘法 .203.2.2 非线性曲线拟合 .243.2.3 多项式曲线 .25第 4 章 多项式的摆动 .274.1 多项式摆动介绍 .274.2 影响多项式拟合偏差的因素 .294.2.1 实验数据的不均匀性 .294.2.2 数据的密度 .304.2.3 拟合曲线的适用区间 .304.3 使用多项式拟合的注意事项 .304.3.1 尽量避免高阶多项式的拟合 .30西北工业大学明德学院本科毕业设计论文IV4.3.2 保持密度 .314.3.3 其它的非线性拟合方法 .31第 5 章 全文总结 .32参考文献 .33致 谢 .34毕业设计小结 .35西北工业大学明德学院本科毕业设计论文5第 1 章 绪论1.1 引言在我们实际的实验和勘探中，都会产生大量的数据。为了解释这些数据或者根据这些数据作出预测、判断，给决策者提供重要的依据。需要对测量数据进行拟合，寻找一个反映数据变化规律的函数。在解决实际工程问题和科学实验的过程中，经常需要通过研究某些变量之间的函数关系，帮我们去认识事物内在的规律和本质属性，这些变量间的未知的关系一般隐含在从观测、试验而得到的一组离散的数据之中。所以，是否能够根据一组试验观测数据来找到变量之间的相对准确的函数关系成为了解决工程实际问题的关键。归纳总结数据拟合理论在工程中实际应用，发掘各个数据拟合算法的在实际应用中的应用范围适用性。通过对本项目的研究和分析，使得实际中的工程问题根据不同的需求使用最合适的拟合算法，从而提高拟合的精确度。研究和发展数据拟合理论，发掘各种数据拟合的优化方案。根据离散的数据，我们想要得到连续的函数或更加密集的离散方程与已知数据相吻合。如何选择数学模型，如何减小误差，如何使得逼近函数图像最靠近那些数据点，使得优化拟合算法变得十分重要。1.2 研究背景在实验中，实验和戡测常常会产生大量的数据。为了解释这些数据或者根据这些数据作出预测、判断，给决策者提供重要的依据。需要对测量数据进行拟合，寻找一个反映数据变化规律的函数。数据拟合方法与数据插值方法不同，它所处理的数据量大而且不能保证每一个数据没有误差，所以要求一个函数严格通过每一个数据点是不合理的。数据拟合方法求拟合函数，插值方法求插值函数。这两类函数最大的不同之处是，对拟合函数不要求它通过所给的数据点，西北工业大学明德学院本科毕业设计论文6而插值函数则必须通过每一个数据点。例如，在某化学反应中，测得生成物的质量浓度 y (10 3 g/cm3)与时间 t （min）的关系如表(1-1)所示 表 1-1 某化学反应数据t 1 2 3 4 6 8 10 12 14 16y 4.00 6.41 8.01 8.79 9.53 9.86 10.33 10.42 10.53 10.61显然，连续函数关系y(t)是客观存在的。但是通过表中的数据不可能确切地得到这种关系。何况，由于仪器和环境的影响，测量数据难免有误差。因此只能寻求一个近拟表达式 )(ty（1-1）寻求合理的近拟表达式，以反映数据变化的规律，这种方法就是数据拟合方法。数据拟合需要解决两个问题：第一，选择什么类型的函数 作为拟合函数)(t（数学模型） ；第二，对于选定的拟合函数，如何确定拟合函数中的参数。数学模型应建立在合理假设的基础上，假设的合理性首先体现在选择某种类型的拟合函数使之符合数据变化的趋势（总体的变化规律） 。拟合函数的选择比较灵活，可以选择线性函数、多项式函数、指数函数、三角函数或其它函数，这应根据数据分布的趋势作出选择。为了问题叙述的方便，将例 1 的数据表写成一般的形式表 1-2 函数一般形式t x1 x2 x3 x4 x5 x6 x7 x8 x9 x10y y1 y2 y3 y4 y5 y6 y7 y8 y9 y100 2 6 8 45689 图 1-1 )(ty西北工业大学明德学院本科毕业设计论文71.3 研究意义曲线拟合，俗称拉曲线，是一种把现有数据透过数学方法来代入一条数式的表示方式。科学和工程问题可以通过诸如采样、实验等方法获得若干离散的数据，根据这些数据，我们往往希望得到一个连续的函数（也就是曲线）或者更加密集的离散方程与已知数据相吻合，这过程就叫做拟合。在科学实验或社会活动中，人们常常需要观测很多数据的规律， 通过实验或者观测得到量 x 与 y 的一组数据对（ ）(i=1，2， ，N)，其中 是iyx, ix彼此不同的。人们希望用一类与数据本质规律相适应的解析表达式 来反映量 x 与 y 之间的依赖关系，即在一定意义下“最佳” 地逼近或拟),(cfy合已知数据。 常称作拟合模型，当 c 在 中线性出现时，称为线性模型，),(f f否者称为非线性模型。线性模型是回归模型中最常见的一种，但在实际中，许多现象之间的关系往往并不是线性的，而是呈现某种曲线关系。如服药后血药浓度与时间的关系，病毒剂量与致死率的关系，化学反应的反应物浓度与反应速度的关系。这就产生的曲线拟合，用连续曲线近似地刻画或比拟平面上离散点组所表示的坐标之间的函数关系。用解析表达式逼近离散数据的一种方法。为了提高曲线拟合精度，本文还研究了多项式的摆动问题，从实践的角度分析了产生这些摆动及偏差的因素和特点，总结了在实践中减小这些偏差的处理方法。采用最小二乘法使变量转换后所得新变量离均差平方和最小，并不一定能使原响应变量的离均差平方和最小，所以其模型的拟合精度仍有提高的空间。随着计算机技术的发展，实验数据处理越来越方便。但也提出了新的课题，就是在选择数据处理方法时应该比以往更为慎重。因为稍有不慎，就会非常方便地根据正确的实验数据得出不确切的乃至错误的结论。所以提高拟合的准确度是非常有必要的。1.4 本论文主要内容本文介绍了几种常用的数据拟合方法，并着重对曲线拟合进行了研究，介绍了线性与非线性模型的曲线拟合方法，最小二乘法、牛顿迭代法等。介绍一种科学计算软件 Matlab， Matlab 是集数值计算、符号运算及出色的图形处理、西北工业大学明德学院本科毕业设计论文8程序语言设计等强大功能于一体的科学计算语言。应用 Matlab 处理既克服了最小二乘法计算量大等缺点，又使繁琐、枯燥的数值计算变成种简单、直观的可视化操作过程，且能较准确地标记实验数据点和绘出拟合曲线。并用 Matlab 逐一实现最小二乘法，非线性曲线，多项式曲线的仿真，并着重对最小二乘法进行研究。在传统的曲线拟合基础上，为了提高曲线拟合精度，本文还研究了多项式的摆动问题，从实践的角度分析了产生这些摆动及偏差的因素和特点，总结了在实践中减小这些偏差的处理方法。西北工业大学明德学院本科毕业设计论文9第 2 章 曲线拟合及最小二乘法在科学实验或社会活动中，人们常常需要观测很多数据的规律，通过实验或者观测得到量 x 与 y 的一组数据对（ ）(i=1，2， ，N)，其中 是彼iyx, ix此不同的。人们希望用一类与数据本质规律相适应的解析表达式， 来),(cfy反映量 x 与 y 之间的依赖关系，即在一定意义下“最佳” 地逼近或拟合已知数据。常称作拟合模型，当 c 在 中线性出现时，称为线性模型，否者称为非),(cf f线性模型。2.1 线性模型的曲线拟合已知某函数的若干离散函数值f1，f2，fn ，通过调整该函数中若干待定系数 f(1， 2，m)， 使得该函数与已知点集的差别 (最小二乘意义)最小。如果待定函数是线性，就叫线性拟合 1。下面介绍计算线性拟合的基本方法。2.2 最小二乘法基本原理曲线拟合又称函数逼近，是指对一个复杂函数 f(x)，求出一个简单的便于计算的函数 p(x)，要求使 f(x)与 p(x)的误差在某种度量意义下最小。我们把近似值 和测定值 的插值称为残余误差 ，即 显然，残差的大小)(ixsiy1)(iiyxs事衡量拟合好坏的重要标志。经常采用的三种衡量的准则为：(1)是残差的最大绝对值最小： ；minaxi(2)使残差的绝对值之和最小： ；i(3)使残差的平方和最小： 。2i分析上面的三种准则，准则(1)，(2) 的提法都比较自然，但是由于含有绝对西北工业大学明德学院本科毕业设计论文10值，所以不利于实际计算，而按照准则(3)来确定参数。得到拟合曲线的方法称作曲线拟合的最小二乘法，它的计算比较简单，是工程实际当中常用的一种函数逼近的方法。设给定的一组实验数据 ，及各点的权系数 ，求出自变,.21),(miyxii量 x 与因变量 y 的函数关系 ，最小二乘法不要求 y=s(x)(;0xasn通过测量点，而只要求残余误差最小。设逼近函数为：(2-1)(.)()()(10 xaxaxs n假设给出一组数据 以及对应的一组权：,.21,miyi求 使 .0,i )(,.;(10mxsyn(2-2)mi iinxsap02102 )(),.(的值最小，则 是多元函数 取得极小值的点，这就是最,.*1*n ),.(10nap小二乘法曲线拟合。求由极值的必要条件可得：(2-3) mi ikiiniik nxyxaap010 ,.10),()(.)()(2 根据内积定义引入相应的带权内积记号：(2-4)miiki ikijkj xyy0)(),(,则式（2-3 ）可改写成：(2-5)nkyaankkk ,.1),(),(.),(),(10这是关于参数 的线性方程组，用矩阵的形式表示为：n,.西北工业大学明德学院本科毕业设计论文11(2-6) ),(,)(),(),(),( , )()()( 10101011 0100 nnnnn nya由于 线性无关，故方程组存在唯一的解),.(10xx nka,.0,*于是有： (2-7)nkxaxp0*)()(2.3 用正交多项式作最小二乘拟合如果 是关于点集 带权)(),.(10xnx),.10(mix正交的函数族，即,.)(mi(2-8)kjAxxKikijmiikj ,0)()(),(0(2-9)nxfmiikii ijikka.,1,)(),(02*且平方误差为 2*022)(| knkaAf现在我们根据给定节点 及权函数 ，造出带权 正交mx,.10 x)(x的多项式 ，注意 ，用递推公式表示 ，即)(xpnn)(pn(2-10) .1,.0),()(), 11010 kxpaxppkkk 这里 是首项系数为 1 的 k 次多项式，根据 的正交性，得)(xpk西北工业大学明德学院本科毕业设计论文12(2-11) 1,.0,)()(1,.0,)( )(,()(1021021 nkpxpnkxxpxpkimikiiiik kmiiii ikik下面用归纳法证明这样给出的 是正交的，由（2-11）式第二次及(2-k11)式中 的表达式，有1 0),(,(),(),(),(),( 0000100 pxppaxp现假定 对 及 均成立，要证,slsl,.l ),.1nk对 均成立。由（2-10）式有0)(1skpk,.1(2-12),(),(),(),()( 111 sksksksksk ppapaxp 由归纳法假定 时，20.0),(0),(1sksk另外， 是首项系数为 1 的 s+1 次多项式，它可由)(xps 1,.sp的线性组合表示，而 ，故由归纳法假定又有ks0),(),(sksk于是由（2-12）式，当 时 。2k0),(1skp再看(2-13),(),(),(),( 1111 kkkkkk pxp由假定有 0,1 k pk，0),(1kp ),(,),(),(1011 kkjjkkk pcxpx 西北工业大学明德学院本科毕业设计论文13利用（2-11 ）式中 表达式及以上结果，得k 0),(),(),(),(),( 111 kkkkk pppxpp最后，由（2-7）式有 ),(),(),(1 kkkpx至此已证明了由(2-10)式及(2-11) 式确定的多项式 组成),.0mn一个关于点集 的正交系。ix用正交多项式 的线性组合作最小二乘曲线拟合，只要根据公式(2-10)及)(px(2-11)逐步求 的同时，相应计算出系数Pk,.1,0,)()(),(02* nkxppfami ikii iiikk 并逐步把 累加到 )中去，最后就可得到所求的拟合曲线)*xkxS)(.)( *1*0* xpaxay n这里 n 可事先给定或在计算过程中根据误差确定。用这种方法编程序不用解线性方程组，只用递推公式，并且当逼近次数增加一次时，只要把程序中循环数加 1，其余不用改变。这就是目前用多项式作曲线拟合最后的计算方法。2.4 非线性模型的曲线拟合当前研究的非线性模型主要是指参数或自变量是非线性的，形式复杂多样，常见的有多项式形式、双曲线形式、对数形式、幂函数形式等等，更复杂的有修正指数曲线 Compterz 曲线以及 Logistic 曲线等。如何根据数据的大致规律来选择合适的模型，是拟合的关键。总的来说有两中可参考的方法：一是根据散点图来确定类型，即由散点图的形状大体确定模型类型；二是根据专业知识和经验，判断研究的数据曲线属于什么类型。现在研究非线性模型的方法用得最多的就是最小二乘法。2.4.1 常见非线性模型对于解释变量是非线性的，但参数之间是线性的模型，可以利用变量直接西北工业大学明德学院本科毕业设计论文14代换的方法将模型线性化，通过线性拟合来计算。1.多项式函数模型多项式函数形式 uxxyk.210令 原模型可化为线性形式kzxz,21即可利用多元线性回归分析的方法处理了。这类模型广泛地用于生产和成本函数。例如总成本函数可表示为： 3210iiii xxy其中，y 表示总成本，表示产出。2双曲线模型双曲线函数形式 uxy103.双对数函数模型函数形式 ylnln10Exyxdx/)(ln*1所以弹性为一常数。它表示 x 变动 1%，y 变动 了。由于这个特%1殊的性质，双对数模型又称为不变弹性模型。4.半对数函数模型函数形式 uxy10lnuxyln对于线性-对数模型 lxydxy/)(ln*1 x1它表示 x 变动 1%，y 将变动 个单位的绝对量。即 y 的绝对变化量等1于 乘以 x 的相对变化量。15.逻辑斯蒂（Logistic）曲线西北工业大学明德学院本科毕业设计论文15函数形式 xbeay1令 y x则有 bxa6.指数曲线函数形式 xy两边取对数得： balogllog令 ly 则有 xba7.幂函数曲线函数形式 bdy两边取对数得： xlogllog令 ylx da则有 bxa8.龚伯兹（Gompertz ）曲线函数形式 bxdey两边取对数得： ln令 yl da则有 bx2.4.2 牛顿迭代无论采取什么方式变换都不可能实现线性化，这样的模型称为不可线性化模型 2。对于不可线性化模型，一般采用高斯一牛顿迭代法进行参数估计，即西北工业大学明德学院本科毕业设计论文16借助于泰勒级数展开式进行逐次的线性近似估计。第一步：做 Logit-Ln 线性回归，求 ， ， x 和 p 的初值。此时 x 不1A0能为 0 值，若输入的 x 有 0 值，则将其设为一小值（例如：0.00001） 。首选将原方程变形为如下线性形式：xpyAlnlln010将 初值设为输入的 y 值的最大值加 1， 的初值设为输入的 y 值的最小0A值减 0.1。通过简单的直线拟合即可求出 p 和 的初值。0x第二步：对 Logistic 方程四个参数求偏微分，得到 y 对给定系数的增量（， ， ， ）的泰勒级数展开式。 1A2xppxAy01ppxpxy02010 200100lnppxAy泰勒级数展开式为： )( 0010 pyxAyy 由此，将曲线回归转化为多元线性回归，通过迭代计算，得到四个参数的变量 ， ， x， p，逐步修正四参数的值。多元线性回归与多项式1A2拟合方法相同。每一次迭代可计算出参数变量值，新的参数值为原参数值与变量值的叠加。西北工业大学明德学院本科毕业设计论文17第三步：为保证迭代收敛，在计算相关系数时，引入一系数 a，初值设为2，将 a 与参数的变量矩阵相乘，计算相关系数。 a=a/2，循环 10 次，每次 a 的值减半。取循环中得到的相关系数最大的变量矩阵 ， ， x， p。1A2第四步：默认总的迭代次数为 1000 次，或者当相关系数不再减小时，则迭代停止。返回得到的四参数值。西北工业大学明德学院本科毕业设计论文18第 3 章 基于 MALTAB 实现最小二乘法3.1 Matlab 简介 3.1.1 Matlab 的概况 Matlab 是矩阵实验室（ Matrix Laboratory）之意。除具备卓越的数值计算能力外，它还提供了专业水平的符号计算，文字处理，可视化建模仿真和实时控制等功能。 Matlab 的基本数据单位是矩阵，它的指令表达式与数学，工程中常用的形式十分相似，故用 Matlab 来解算问题要比用 C，FORTRAN 等语言完相同的事情简捷得多。当前流行的 Matlab 5.3/Simulink 3.0包括拥有数百个内部函数的主包和三十几种工具包(Toolbox).工具包又可以分为功能性工具包和学科工具包。功能工具包用来扩充 Matlab 的符号计算，可视化建模仿真，文字处理及实时控制等功能。学科工具包是专业性比较强的工具包，控制工具包，信号处理工具包，通信工具包等都属于此类。开放性使 Matlab 广受用户欢迎，除内部函数外，所有 Matlab 主包文件和各种工具包都是可读可修改的文件，用户通过对源程序的修改或加入自己编写程序构造新的专用工具包。 3.1.2 Matlab 的语言特点 一种语言之所以能如此迅速地普及，显示出如此旺盛的生命力，是由于它有着不同于其他语言的特点，正如同 FORTRAN 和 C 等高级语言使人们摆脱了西北工业大学明德学院本科毕业设计论文19需要直接对计算机硬件资源进行操作一样，被称作为第四代计算机语言的Matlab，利用其丰富的函数资源，使编程人员从繁琐的程序代码中解放出来。Matlab 最突出的特点就是简洁。Matlab 用更直观的，符合人们思维习惯的代码，代替了 C 和 FORTRAN 语言的冗长代码。Matlab 给用户带来的是最直观，最简洁的程序开发环境。以下简单介绍一下 Matlab 的主要特点。 语言简洁紧凑，使用方便灵活，库函数极其丰富。Matlab 程序书写形式自由，利用起丰富的库函数避开繁杂的子程序编程任务，压缩了一切不必要的编程工作。由于库函数都由本领域的专家编写，用户不必担心函数的可靠性。可以说，用 Matlab 进行科技开发是站在专家的肩膀上。 具有 FORTRAN 和 C 等高级语言知识的读者可能已经注意到，如果用FORTRAN 或 C 语言去编写程序，尤其当涉及矩阵运算和画图时，编程会很麻烦。例如，如果用户想求解一个线性代数方程，就得编写一个程序块读入数据，然后再使用一种求解线性方程的算法（例如追赶法）编写一个程序块来求解方程，最后再输出计算结果。在求解过程中，最麻烦的要算第二部分。解线性方程的麻烦在于要对矩阵的元素作循环，选择稳定的算法以及代码的调试动不容易。即使有部分源代码，用户也会感到麻烦，且不能保证运算的稳定性。解线性方程的程序用 FORTRAN 和 C 这样的高级语言编写，至少需要四百多行，调试这种几百行的计算程序可以说很困难。3.1.3 Matlab 工作界面图 3-1 Matlab 工作界面Matlab 是 matrix&laboratory 两个词的组合，意为矩阵工厂（矩阵实验室）是由美国 mathworks 公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为西北工业大学明德学院本科毕业设计论文20科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、FORTRAN）的编辑模式，代表了当今国际科学计算软件的先进水平。Matlab 和 Mathematica、Maple 并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。Matlab 可以进行 矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。Matlab 的基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，故用 Matlab 来解算问题要比用 C，FORTRAN 等语言完成相同的事情简捷得多，并且 Matlab 也吸收了像 Maple 等软件的优点，使 Matlab 成为一个强大的数学软件。在新的版本中也加入了对 C，FORTRAN，C+，JAVA的支持。可以直接调用，用户也可以将自己编写的实用程序导入到 Matlab 函数库中方便自己以后调用，此外许多的 Matlab 爱好者都编写了一些经典的程序，用户直接进行下载就可以用。3.1.4 优势特点(1) 高效的数值计算及符号计算功能，能使用户从繁杂的数学运算分析中解脱出来；(2) 具有完备的图形处理功能，实现计算结果和编程的可视化；(3) 友好的用户界面及接近数学表达式的自然化语言，使学者易于学习和掌握；(4) 功能丰富的应用工具箱(如信号处理工具箱、通信工具箱等) ，为用户提供了大量方便实用的处理工具。3.2 用 MALTAB 实现曲线拟合3.2.1 最小二乘法Matlab 是一种功能强大的系统分析和仿真工具，我们选用它作为实现曲线拟合的软件工具。用 Matlab 语言编程实现最小二乘法的思路 3:(1)输入各参量 x、y 的测量值(以数组形式输入，这样便于在计算过程中引西北工业大学明德学院本科毕业设计论文21用)；(2)用 Matlab 语言中的 plot 函数 x、y 的曲线关系图，以此图对比典型曲线图，选择合适的经验公式；(3)按照上例中的方法，选一个系数 a，求 Q(a，b)对它的偏导数，求出其计算表达式；(4)编写 Matlab 的 M 函数，用来完成经验公式中待定系数 a 的计算，该函数输入量为 x、y、b，输出量为 a、Q，按照由最小二乘法推导出的公式代入数值由 x、y、b 计算 a、Q;(5)改变 b 的取值，多次调用该 M 函数，比较结果中的 Q 值，最小的 Q 值所对应的 a、b 值即为所求。改变 b 的取值 #这部分工作也可编一个循环函数，输入 b 可能取的区间，计算不同 b 对应的 Q，再进行比较，保留使 Q 最小的 b 及对应的 a。但通常 b的改变对 Q 的影响不是线性的，为方便观察结果并选择适当的 b， 改变 b 的取值#这部分工作最好还是编程者自己完成。最后，只要将得到的函数图像和 x、y 的曲线关系图进行对比，就可以直观的看到拟合的效果了。另外 Matlab 语言提供了一个函数，可以完成线性曲线拟合，这就是函数polyfit。函数 polyfit 的输入量为 x、y、n，其中 x、y 即为需要建立相互关系的两个量的测量值，以数组的形式输入，n 为多项式的次数;输出的是多项式系数的行向量，而得到的多项式是降幂的 4。如某一消失模滤波器中螺钉归一化插入深度与其对应的电纳值的数据表 3-1 螺钉归一化插入深度与其对应的电纳值h CB h CB0.2411 4.029 0.1829 6.2190.2328 4.226 0.1746 6.8290.2245 4.439 0.1663 7.6280.2162 4.690 0.1580 8.646西北工业大学明德学院本科毕业设计论文220.2079 4.984 0.1497 10.3330.1995 5.296 0.1413 13.6680.1912 5.724 0.1330 23.038现在按照上面介绍的方法用 Matlab 语言求Bc= f (h)。(1)输入数据:h= 0. 2411 0. 2328 0. 2245 0. 2162 0. 2078 0. 1995 0