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

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

ID:20799738

大小:120.70 KB

页数:5页

时间:2018-10-16

最小二乘法综述及举例.doc_第1页
最小二乘法综述及举例.doc_第2页
最小二乘法综述及举例.doc_第3页
最小二乘法综述及举例.doc_第4页
最小二乘法综述及举例.doc_第5页
资源描述:

《最小二乘法综述及举例.doc》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库

1、最小二乘法综述及算例一最小二乘法的历史简介1801年,意大利天文学家朱赛普·皮亚齐发现了第一颗小行星谷神星。经过40天的跟踪观测后,由于谷神星运行至太阳背后,使得皮亚齐失去了谷神星的位置。随后全世界的科学家利用皮亚齐的观测数据开始寻找谷神星,但是根据大多数人计算的结果来寻找谷神星都没有结果。时年24岁的高斯也计算了谷神星的轨道。奥地利天文学家海因里希·奥尔伯斯根据高斯计算出来的轨道重新发现了谷神星。高斯使用的最小二乘法的方法发表于1809年他的著作《天体运动论》中。经过两百余年后,最小二乘法已广泛应用与科学实验和工程技术中,随着现代电子计算机的普及与发展,这个方

2、法更加显示出其强大的生命力。二最小二乘法原理最小二乘法的基本原理是:成对等精度测得的一组数据,是找出一条最佳的拟合曲线,似的这条曲线上的个点的值与测量值的差的平方和在所有拟合曲线中最小。设物理量y与1个变量间的依赖关系式为:。其中是n+l个待定参数,记其中是测量值,是由己求得的以及实验点得出的函数值。在设计实验时,为了减小误差,常进行多点测量,使方程式个数大于待定参数的个数,此时构成的方程组称为矛盾方程组。通过最小二乘法转化后的方程组称为正规方程组(此时方程式的个数与待定参数的个数相等)。我们可以通过正规方程组求出a最小二乘法又称曲线拟合,所谓“拟合”即不要求所

3、作的曲线完全通过所有的数据点,只要求所得的曲线能反映数据的基本趋势。三曲线拟合曲线拟合的几何解释:求一条曲线,使数据点均在离此曲线的上方或下方不远处。(1)一元线性拟合设变量y与x成线性关系,先已知m个实验点,求两个未知参数。令,则应满足。即化简得从中解出(2)多元线性拟合设变量y与n个变量的内在联系是线性的,即有下式设的第i次测量值为,对应的函数值为,则偏差平方和为使s去得最小值的方程组即。(4)将实验数据代入(4)式,即得。(3)多项式拟合科学实验后得到一组数据时,常会遇到因变量y与自变量x之间根本不存在线性关系。此可以考虑用一个n次多项式来拟合y与x之间的

4、函数关系。对于n次多项式,令,则可将其化为线性形式:对于i=1,2,...,m个实验点有,代入(3)式有从而得出多项式的最小二乘法拟合的方程写成矩阵的形式即为从中可以解出。(4)指数函数拟合此时拟合函数具有形式(a,b为待定系数)。两端取自然对数有令则(*)式化为线性形式再利用(1)式和(2)式,即可求出。从而有。故。四最小二乘法应用举例例:已知某铜棒的电阻与温度关系为:。实验测得7组数据(见表1)如下:试用最小二乘法求出参量R0、a以及确定它们的误差。表1t/℃19.125.130.136.040.045.150.1Rt/W76.3077.8079.7580.

5、8082.3583.9085.10此例中只有两个待定的参量R0和a,为得到它们的最佳系数,所需要的数据有n、、、、和六个累加数,为此在没有常用的科学型计算器时,通过列表计算的方式来进行,这对提高计算速度将会有极大的帮助(参见表2),并使工作有条理与不易出错。其中表内双线右边的计算是为了确定R0和a的误差项用的。表2it/℃(xi)Rt/W(yi)t×t(x2i)Rt´Rt(y2i)t×Rt(xiyi)R计算/Wui/Wui2×10-4123456719.125.130.136.040.045.150.176.3077.8079.7580.8082.3583.90

6、85.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))则可得:把测

7、量数据代入式(13)和(15)中可求出相关系数说明:电阻Rt与温度t的线性关系良好,所以取R0的有效数字与R对齐,即R0=70.76W;又因为t7-t1=31.0℃,R7-R1=8.80W,取k有效数字为以上两个差值中较少的位数3位,则k=0.288W/°C。由此可以得到电阻与温度的相关关系为:按补充资料中的公式计算k和b的不确定度,可得故,则参考文献:1.《最小二乘法与测量平差》郭禄光,樊功瑜著同济大学出版社19852.《近代最小二乘法》测绘出版社19803.《最小二乘法的拟合及其应用》邓亮章兰州教育学院学报2012.114.《最小二乘法的创立及其思想方法》贾

8、小勇,徐传胜,白欣西北大

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。