最小二乘法概述

  • 格式:doc
  • 大小:448.43 KB
  • 文档页数:6

最小二乘法

一、简介

最小二乘法,又称最小平方法,是一种数学技术。它通过最小误差的平方和寻找数据函数的最佳匹配。最小二乘法是提供“观测组合”的主要工具之一,它依据对某事件的大量观测而获得“最佳”结果或“最可能”表现形式。如已知两变量为线性关系bxay,对其进行)2(nn次观测而获得n对数据。若将这n对数据代入方程求解a,b之值则无确定解。最小二乘法提供了一个求解方法,其基本思想就是寻找“最接近”这n个观测点的直线。最小二乘法不仅是19世纪最重要的统计方法,而且还可以称为数理统计学之灵魂。相关回归分析、方差分析和线性模型理论等数理统计学的几大分支都以最小二乘法为理论基础。作为其进一步发展或纠正其不足而采取的对策,不少近现代的数理统计学分支也是在最小二乘法基础上衍生出来的。

最小二乘法之于数理统计学,有如微积分之于数学,这并非夸张之辞。统计学应用的几个分支如相关分析、回归分析、方差分析和线性模型理论等,其关键都在于最小二乘法的应用不少现代的统计学研究是在此法的基础上衍生出来,作为其进一步发展或纠正其不足之处而采取的对策,如回归分析中一系列修正最小二乘法而产生的估计方法等就是最好的例子。

二、创立思想

勒让德在先驱者解线性方程组的基础上,以整体的思想方法创立了最小二乘法;高斯由寻找随机误差函数为突破,以独特的概率思想导出了正态分布,详尽地阐述了最小二乘法的理论依据。

最小二乘法(OLSE)的思想就是要使得观测点和估计点的距离平方和达到最小,在各方程的误差之间建立一种平衡,从而防止某一极端误差,对决定参数的估计值取得支配地位,有助于揭示系统的更接近真实的状态。这里的“二乘”指的是用平方来度量观测点与估计点的远近,“最小”指的是参数的估计值要保证各个观测点与估计点的距离的平方和达到最小。

三、原理

设一组数据(,)iixy(1,2,,)in,现用近似曲线)(xy拟合这组数据,“拟合得最好”的标准是所选择的x在ix处的函数值ix(1,2,,)in与iy(1,2,,)in相差很小,即偏差(也称残差)iixy(1,2,,)in都很小.一种方法是使偏差之和1niiixy-很小来保证每个偏差都很小.但偏差有正有负,在求和的时候可能相互抵消.为了避免这种情况,还可使偏差的绝对值之和1||niiixy为最小.但这个式子中有绝对值符号,不便于分析讨论.由于任何实数的平方都是正数或零,因而我们可选择使“偏差平方和21niiixy[()]最小”的原则来保证每个偏差的绝对值都很小,从而得到最佳拟合曲线y=x.这种“偏差平方和最小”的原则称为最小二乘原则,而按最小二乘法原则拟合曲线的方法称为最小二乘法或称最小二乘曲线拟合法.

一般而言,所求得的拟合函数可以使不同的函数类,拟合曲线x都是由m个线性无关函数1x,2x ,„, mx的线性组合而成,即

1122mmxaxaxax…)1(nm,

其中1a,2a,„,ma为待定系数.线性无关函数1x,2x ,„mx,称为基函数,常用的基函数有:

多项式:1,x, 2x,„,mx;

三角函数: sinx,sin2x,„,sinmx;

指数函数:xxxmeee,,,21,x2e,„,xme.

最小二乘法又称曲线拟合,所谓“ 拟合” ,即不要求所作的曲线完全通过所有的数据点,只要求所得的近似曲线能反映数据的基本趋势,它的实质是离散情况下的最小平方逼近.

四、运用曲线拟合做最小二乘法

1 一元线性拟合

已知实测到的一组数据(,)iixy(1,2,,)in,求作这组数据所成的一元线性关系式.设线性关系式为yabx,求出a和b即可.

法一:即要满足则)(令,0,0,,12bsasbabxaysniii,则,ab要满足sa=0,sb=0.即

11()()niiiniiiisyabxasyabxxb=-2=0=-2=0

化简得

112111nniiiinnniiiiiiibaxynnaxbxxy1+=+=

从中解出

1112211111nnniiiiiiinniiiinniiiinxyxybnxxbayxnn-=-=- (1)

法二:将ix,iy代入yabx得矛盾方程组

1122nyabxyabxyabxn (2)

令A=12111nxxx,B=12nyyy,则(2)式可写成bBAa,则对应的正规方程组为TTabABAA,所以ab=1()TTAAAB,其中A称为结构矩阵,B称为数据矩阵,TAA称为信息矩阵,TAB 称为常数矩阵.

2 多元线性拟合

设变量y与n个变量1x,2x,„,nx(1n)内在联系是线性的,即有如下关系式

njjjxaay10,设jx的第i次测量值为ijx,对应的函数值为iy(1,2,,)im,则偏差平方和为s='220111()()mmniiiiijiijyyyaax,为了使s取最小值得正规方程组

011001111011202020mnijijijmnijijiijmnijijinijnsyaaxasyaaxxasyaaxxa (3)

011101111nmmijjijiimnmmikijikjikiijiimaxayxaxxaxy1,2,,kn. (4)

将实验数据(,)iixy代入(4)式,即得maaa,,,10.

3 指数函数拟合

科学实验得到一组数据(,)iixy(1,2,,)in时,还可以考虑用指数函数为基函数来拟合,此时设拟合函数具有形式bxyae(,ab为待定系数).对上式两端取自然对数可得:

lnlnyabx (9)

令Y=lny,0lnba,则(9)式可转化为一元线性函数形式0Ybbx,此时将指数函数拟合转化成了一元线性拟合,利用一元线性拟合中的两种方法均可求出0b和b,继而根据0bae可求出a,从而得出因变量y与自变量x之间的函数关系式0bbxbxyaee

4 对数函数拟合

科学实验得到一组数据(,)iixy(1,2,,)in时,还可以考虑用对数函数为基函数来拟合,此时设拟合函数具有形式lnyabx(0)x(,ab为待定系数).0b时,y随x增大而增大,先快后慢;0b时,y随x增大而减小,先快后慢.当以y和lnx绘制的散点图呈直线趋势时,可考虑采用对数函数描述y与x之间的非线性关系,式中的b和a分别为斜率和截距.这时令X=lnx,就可以利用一元线性拟合的方法来求解.更一般的对数函数还可设为ylnabxk,式中k为一常量.

五 举例

例1

使电流通过2的电阻,用伏特表测量电阻两端的电压V.测得数据如下表:

tI/A 1 2 4 6 8 10

tV/V 1.8 3.7 8.2 12.0 15.8 20.2

试用最小二乘法建立I与V之间的一元经验公式(有效数字保留到小数点后第3位).

解:可取一次线性关系式VabI作为I与V之间的一元经验公式.

将数据代入得矛盾方程组

1.823.748.2612.0815.81020.2abababababab

令1112141618110A,1.83.78.212.015.820.2B,则上述矛盾方程组可写成矩阵形式0aABb

由此得出其正规方程组0TTaAAABb,将数据代入即得63161.7031221442.4ab,解之得0.212.032ab,故所求经验公式为0.2152.VI.

例2 在在开发一种抗过敏性的新药时,要对不同剂量的药效进行实验.10名患者各服用了该新药的一个特定的剂量.药物消失时立即纪录.观测值列于下表中.x是剂量,y是症状消除持续的日数.用7个不同的剂量, 其中3个剂量重复给两名患者.试给出y与x之间的一元经验公式(保留3位有效数字).

1 2 3 4 5 6 7 8 9 10 ∑

/ixmg 3 3 4 5 6 6 7 8 8 9 59

/iyd 9 5 12 9 14 16 22 18 24 22 151 2ix 9 9 16 25 36 36 49 64 64 81 389

iixy 27 15 48 45 84 96 154 144 192 198

1003

解:可设y与x之间的经验公式为yabx.

由上表可知,

101iix59,101iiy151,101iiixy1003,1021iix389,2101iix3481

再由(1)式可求得,

1010101112101021110101003591512.7410389348110iiiiiiiiiiixyxybxx

101011112.74151591.0710101010iiiibayx

所以y与x之间的经验公式为1.072.74yx.

最小二乘法能将从实验中得出的一大堆看上去杂乱无章的数据中找出一定的规律,拟合成一条曲线来反映所给数据特点。