最小二乘法课程设计.doc

最小二乘法（1） 第 PAGE 15 页 共 NUMPAGES 15 页 1、前 言 1、1背景介绍 最小二乘法最早是由高斯提出的，这是数据处理的一种很有效的统计方法。高斯用这种方法解决了天文学方面的问题，特别是确定了某些行星和彗星的天体轨迹。这类天体的椭圆轨迹由5个参数确定，原则上，只要对它的位置做5次测量就足以确定它的整个轨迹。但由于存在测量误差，由5次测量所确定的运行轨迹极不可靠，相反，要进行多次测量，用最小二乘法消除测量误差，得到有关轨迹参数的更精确的值。最小二乘法近似将几十次甚至上百次的观察所产生的高维空间问题降到了椭圆轨迹模型的五维参数空间。 最小二乘法普遍适用于各个科学领域，它在解决实际问题中发挥了重要的作用。它在生产实践、科学实验及经济活动中均有广泛应用。比如说，我们引入等效时间的概念,根据Arrhenius函数和指数函数研究水化热化学反应速率随温度的变化,最后采用最小二乘法回归分析试验数据,确定绝热温升和等效时间的关系式。 1、2问题引入 为了更好地掌握最小二乘法，我们引入以下两个问题： (1)假设已知一组二维数据（），（i=1,2,3···n），怎样确定它的拟合曲线y=f(x)（假设为多项式形式f(x)=）,使得这些点与曲线总体来说尽量接近？ (2)若拟合模型为非多项式形式，怎样根据已知的二维数据用最小二乘线性拟合确定其系数，求出曲线拟合函数？ 怎样从给定的二维数据出发，寻找一个简单合理的函数来拟合给定的一组看上去杂乱无章的数据，正是我们要解决的问题。 1、3内容简介 本文主要介绍最小二乘法及其实现方法，通过对其原理的阐述，并以实例的形式说明了如何使用MATLAB函数求解拟合问题。利用算例中的已知二维数据，进行曲线拟合。具体内容如下： (1)最小二乘法的概念 (2)最小二乘法的原理及实现方法 (3)最小二乘法的应用 (4)具体实例的求解 2、最小二乘法的概念 在科学实验的统计方法研究中，往往要从一组实验数（）（i=1,2,3···m）中寻找自变量x与y之间的函数关系y=F(x).由于观测数据往往不准确，此时不要求y=F(x)经过所有点（），而只要求在给定上误差=F（）（i=1,2,3···m）按某种标准最小。 若记=,就是要求向量的范数最小。如果用最大范数，计算上困难较大，通常就采用Euclid范数作为误差度量的标准。 关于最小二乘法的一般提法是：对于给定的一组数据（） (i=0,1,…m)要求在函数空间Φ=span{}中找一个函数S*(x)，使加权的误差平方和 =最小，其中，是[a,b]上的权函数，它表示反应数据（）在实验中所占数据的比重。 我们说，S(x)= (nm)这就是一般的最小二乘逼近，用几何语言说就是曲线拟合的最小二乘法。 注意这里的，是线性无关的。 在研究两个变量之间的关系时，可以用回归分析的方法进行分析。当确定了描述两个变量之间的回归模型后，就可以使用最小二乘法估计模型中的参数，进而建立经验方程。 为了通过试验数据来估计参数的值，可以采用许多统计方法，而最小二乘法是目前最常用、最基本的。 3、最小二乘法的原理及实现方法 3、1最小二乘法的原理 简单地说，最小二乘的思想就是要使得观测点和估计点的距离的平方和达到最小.这里的“二乘”指的是用平方来度量观测点与估计点的远近（在古汉语中“平方”称为“二乘”），“最小”指的是参数的估计值要保证各个观测点与估计点的距离的平方和达到最小. 例如，对于回归模型y=S(x),若（）（i=1,2,3···m）为收集到的观测数据，则应该用来估计，这里是（）（i=1,2,3···m）的估计值。它们之间距离的平方和就是。 进而最小二乘估计量就是使 === (*) 达到最小值的参数。 从计算的角度看，最小二乘法与插值法类似，都是处理数据的算法.但从创设的思想看，二者却有本质的不同。前者寻求一条曲线，使其与观测数据“最接近”，目的是代表观测数据的趋势；后者则是使曲线严格通过给定的观测数据，其目的是通过来自函数模型的数据来近似刻画该函数.在观测数据带有测量误差的情况下，就会使得这些观测数据偏离函数曲线，结果使得与观测数据保持一致的插值法不如最小二乘法得到的曲线更符合客观实际. 3、2最小二乘法的实现方法 3、2、1一般的最小二乘曲线拟合 给定一组测量数据{（），i=0,1,2,…,m}，基于最小二乘原理，求得变量x和y之间的函数关系f(x,A),使它最佳地逼近已知数据。其中A=（）是一些待定参数。 为了是问题的提法更有一般性，通常把最小二乘法中的都考虑为加权平方和，即 = 其中，是[a,b]上的权函数，它表示反应数据（）在实验中所占数据的比重。 选