参数估值理论大作业.docVIP

最小二乘法及其在matlab中的应用 摘要：物理量之间的函数关系在实际研究工作中有很重要的作用。本文首先介绍了最小二乘法拟合的基本原理， 其次介绍了用Matlab实现曲线拟合以得到函数关系的方法和步骤，最后举例详细介绍了该方法的应用。 关键词：最小二乘法；Matlab;曲线拟合 引言 在现代科学研究中，物理量之间的相互关系通常是用函数来描述的。有些函数关系是由经典理论分析推到得出的，这些函数关系不仅为我们进一步的分析研究工作提供了物理的理论基础，也是我们可以十分方便的运用丰富的数学知识来解决物理问题。在现实的物理研究过程中，有一些问题很难由经典物理理论推导出物理量的函数表达式，或者此推导出的表达式也十分复杂，不利于进一步的分析。但由于研究需要，又很希望能得到这些量之间的函数关系，这时就可以利用曲线拟合的方法，用实验数据结合数学方法得到物理量之间的曲线函数表达式。 在处理实验数据时，成用最小二乘法，通常使用计算器的统计功能来间接计算，计算量大，不但极易出现按键错误，而且不方便核查校对，耗时费力，同时也无法找出误差较大甚至错误的数据点。Matlab是Math Works公司推出的一种科学计算软件，是集数值计算、符号运算及出色的图形处理、程序语言设计等强大功能于一体的科学计算语言。用Matlab处理实验数据仅需编写十几行几乎像通常笔算式的简练程序，运行后就可得到所需的结果。应用它既克服了最小二乘法计算量大等缺点，又使繁琐、枯燥的数值计算变成种简单、直观的可视化操作过程，且能较准确地标记实验数据点和绘出拟合曲线。 最小二乘法拟合的基本原理 曲线拟合又称为函数逼近，是指对一个复杂函数f(x),求出一个简单的便于计算的函数p(X),它要求使f(x)与p(x)的误差在某种度量意义下最小。 我们把近似值yi*和实测值yi之差称为残差，即： =yi-yi* 显然，残差的大小是衡量拟合好坏的重要标志。我们经常采用的三种衡量的准则为： ⑴使残差的最大绝对值最小，即： Max| |=min ⑵使残差的绝对值之和最小，即： ||=min ⑶使残差的平方和最小，即： ||2=min 分析上面的三种准则，准则⑴、⑵的提法都比较自然，但是由于含有绝对值，所以不利 于实际计算，而按照准则⑶来确定参数，得到拟合曲线的方法称作曲线拟合的最小二乘法，它的计算比较简单，是工程实际当中常用的一种函数逼近的方法。 曲线拟合的最小二乘法问题可以描述为：根据已知的数据(xi,yi) (i=1,2,﹍,n),选取一个近似函数，使得||2=[yi-2 最小。这种求近似函数的方法称为曲线拟合的最小二乘法，函数称为这组数据的最小二乘法。 用Matlab实现曲线拟合 Matlab软件是是集数值计算、符号运算及出色的图形处理、程序语言设计等强大功能于一体的科学计算语言。用Matlab处理实验数据仅需编写十几行几乎像通常笔算式的简练程序，运行后就可得到所需的结果。应用它既克服了最小二乘法计算量大等缺点，又使繁琐、枯燥的数值计算变成种简单、直观的可视化操作过程，且能较准确地标记实验数据点和绘出拟合曲线。 本文通过举例详细介绍用Matlab实现最小二乘法处理实验数的方法。 算法： ⑴输入各参量x,y 的测量值； ⑵用Matlab语言中的plot函数x,y的曲线散点趋势图，以此图对比典型曲线图，选择合适的经验公式； ⑶编写Matlab的函数，用来完成经验公式中待定系数的计算。 Matlab程序： 将上面的算法编写成Matlab程序计算数据最小二乘拟合系数。 Function S=squar_least(x,y,n,w) 数据的最小二乘拟合，其中x,y为数据（x，y）坐标； n为数据拟合的次数,缺省是n=1; w为权值，缺省值w=1;s为数据拟合的系数。 举例说明 例：已知在测量小车位移（yi）和时间（xi）关系时，测得的数据如下： 表1 小车时间（xi）和位移（yi）关系 xi 0 1 2 3 4 5 6 7 8 9 yi 0 2 4 7 8 9 12 14 15 18 解：⑴首先画出数据的散点图，输入 x=[0 1 2 3 4 5 6 7 8 9];y=[0 2 4 7 8 9 12 14 15 18]; subplot(1,2,1);plot(x,y,’o’) grid on ⑵选择经验公式： 我们从图中可以看出其基本为直线趋势，所以拟合的曲线应该是一次直线，我们分别用三次样条插值法和最小二乘法一次拟合进行比较，下面是借助Matlab工具进行作图，其程序如下： x=[0 1 2 3 4 5 6 7 8 9];y=[0 2