最小二乘法综述及举例.docx

最小二乘法综述及算例一最小二乘法的历史简介1801年，意大利天文学家朱赛普·皮亚齐发现了第一颗小行星谷神星。经过40天的跟踪观测后，由于谷神星运行至太阳背后，使得皮亚齐失去了谷神星的位置。随后全世界的科学家利用皮亚齐的观测数据开始寻找谷神星，但是根据大多数人计算的结果来寻找谷神星都没有结果。时年24岁的高斯也计算了谷神星的轨道。奥地利天文学家海因里希·奥尔伯斯根据高斯计算出来的轨道重新发现了谷神星。高斯使用的最小二乘法的方法发表于1809年他的著作《天体运动论》中。经过两百余年后，最小二乘法已广泛应用与科学实验和工程技术中，随着现代电子计算机的普及与发展，这个方法更加显示出其强大的生命力。二最小二乘法原理最小二乘法的基本原理是：成对等精度测得的一组数据，是找出一条最佳的拟合曲线，似的这条曲线上的个点的值与测量值的差的平方和在所有拟合曲线中最小。设物理量y与1个变量间的依赖关系式为:。其中是n +l 个待定参数,记其中 是测量值, 是由己求得的以及实验点得出的函数值。在设计实验时, 为了减小误差, 常进行多点测量, 使方程式个数大于待定参数的个数, 此时构成的方程组称为矛盾方程组。通过最小二乘法转化后的方程组称为正规方程组(此时方程式的个数与待定参数的个数相等) 。我们可以通过正规方程组求出a最小二乘法又称曲线拟合, 所谓“ 拟合” 即不要求所作的曲线完全通过所有的数据点, 只要求所得的曲线能反映数据的基本趋势。三曲线拟合曲线拟合的几何解释: 求一条曲线, 使数据点均在离此曲线的上方或下方不远处。（1）一元线性拟合设变量y 与x 成线性关系,先已知m个实验点，求两个未知参数。令，则应满足。即化简得从中解出（2）多元线性拟合设变量y与n 个变量的内在联系是线性的，即有下式设的第i次测量值为，对应的函数值为，则偏差平方和为使s去得最小值的方程组即。（4）将实验数据代入（4）式，即得。（3）多项式拟合 科学实验后得到一组数据时，常会遇到因变量y 与自变量x 之间根本不存在线性关系。此可以考虑用一个n 次多项式来拟合y与x 之间的函数关系。对于n次多项式，令，则可将其化为线性形式：对于i=1,2,...,m个实验点有，代入（3）式有从而得出多项式的最小二乘法拟合的方程写成矩阵的形式即为从中可以解出。（4）指数函数拟合 此时拟合函数具有形式（a，b为待定系数）。两端取自然对数有令则（*）式化为线性形式 再利用（1）式和（2）式，即可求出。从而有。故。四最小二乘法应用举例例：已知某铜棒的电阻与温度关系为：。实验测得7组数据（见表1）如下：试用最小二乘法求出参量R0、 以及确定它们的误差。表 1t / ℃19.125.130.136.040.045.150.1Rt / 76.3077.8079.7580.8082.3583.9085.10此例中只有两个待定的参量R0和，为得到它们的最佳系数，所需要的数据有n、、、、和六个累加数，为此在没有常用的科学型计算器时，通过列表计算的方式来进行，这对提高计算速度将会有极大的帮助（参见表2），并使工作有条理与不易出错。其中表内双线右边的计算是为了确定R0和的误差项用的。表 2it / ℃( xi )Rt / ( yi )t×t( x2i )Rt Rt( y2i )t×Rt( xi yi )R计算 / i / i2×10-4123456719.125.130.136.040.045.150.176.3077.8079.7580.8082.3583.9085.10364.8630.0906.01296.01600.02034.02510.05821.76052.86360.16528.66781.57039.27242.01457.31952.82400.52908.83294.03783.94263.576.2677.9979.4381.1382.2883.7585.19+0.04-0.19+0.32-0.33+0.07+0.15-0.09163611024108949225817245.5566.009340.84582620060.82845×10-4根据表2中所求得的数据，代入公式（12））则可得：把测量数据代入式（13）和（15）中可求出相关系数说明：电阻Rt与温度t的线性关系良好，所以取R0的有效数字与R对齐，即R0＝70.76；又因为t7－t1 = 31.0 ℃，R7－R1 = 8.80，取k有效数字为以上两个差值中较少的位数3位，则k = 0.288/C。由此可以得到电阻与温度的相关关系为：按补充资料中的公式计算k 和b的不确定度，可得故，则参考文献：《最小二乘法与测量平差》 郭禄光，樊功瑜著 同济大学出版社 1985《近代最小二乘法》测绘出版社