用最小二乘法拟合数据并求均方偏差

《用最小二乘法拟合数据并求均方偏差》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、用最小二乘法拟合数据并求均方偏差摘要：数值分析是研究分析用计算机求解数学计算问题的数值计算方法及其理论的学科，是数学的一个分支，它以数字计算机求解数学问题的理论和方法为研究对象。本文通过对一组数据用origin先得到散点图，然后根据散点图来预测函数。接着用matlab软件采用最小二乘法拟合，得到两个不同的函数，并计算它们的均方偏差，以便比较这两个拟合函数的优劣。关键字：数值分析;origin;matlab;最小二乘法;均方偏差Abstraction:Numericalanalysisisacomputationalmethodanditsnumericalcompu

2、tationproblemissolvedbycomputeranalysisofmathematicalresearchsubject,isabranchofmathematics,itisbasedonthetheoryandmethodologyofdigitalcomputertosolvemathematicalproblemsastheresearchobject.Thisarticlethroughtoasetofdatawithorigintogetscatterplot,thenpredictfunctionaccordingtothescatte

3、rplot.Thenusingleastsquaresfittingwiththematlabsoftware,gettwodifferentfunction,andthemeansquaredeviation,theycalculatedtocomparetheadvantagesanddisadvantagesofthetwofittingfunction.Keyword:Numericalanalysis;origin;matlab;theleastsquaremethod;MeanSquareError1、引言数值分析主要介绍现代科学计算中常用的数值计算方法

4、及其基本原理，研究并解决数值问题的近似解，是数学理论与计算机和实际问题的有机结合。随着科学技术的迅速发展，运用数学方法解决科学研究和工程计算领域中的实际问题，已经得到普遍重视。数学建模是数值分析联系实际的桥梁。在数学建模过程中，无论是模型的建立还是模型的求解都要用到数值分析课程中所涉及的算法,如插值方法、最小二乘法、拟合法等，那么如何在数学建模中正确的应用数值分析内容，就成了解决实际问题的关键。其中最小二乘法（又称最小平方法）是一种数学优化技术。它通过最小化误差的平方和寻找数据的最佳函数匹配。利用最小二乘法可以简便地求得未知的数据，并使得这些求得的数据与实际数据之间

5、误差的平方和为最小。最小二乘法还可用于曲线拟合。其他一些优化问题也可通过最小化能量或最大化熵用最小二乘法来表达。本文主要运用的数值分析算法就是最小二乘法。2、用最小二乘法解决问题2.1、最小二乘法的历史简介1801年，意大利天文学家朱赛普·皮亚齐发现了第一颗小行星谷神星。经过40天的跟踪观测后，由于谷神星运行至太阳背后，使得皮亚齐失去了谷神星的位置。随后全世界的科学家利用皮亚齐的观测数据开始寻找谷神星，但是根据大多数人计算的结果来寻找谷神星都没有结果。时年24岁的高斯也计算了谷神星的轨道。奥地利天文学家海因里希·奥尔伯斯根据高斯计算出来的轨道重新发现了谷神星。高斯使

6、用的最小二乘法的方法发表于1809年他的著作《天体运动论》中。法国科学家勒让德于1806年独立发现“最小二乘法”，但因不为世人所知而默默无闻。勒让德曾与高斯为谁最早创立最小二乘法原理发生争执。1829年，高斯提供了最小二乘法的优化效果强于其他方法的证明，因此被称为高斯-马尔可夫定理。2.2、最小二乘法及其计算在函数的最佳平方逼近中f(x)∈C[a,b],如果f(x)只在一组离散点{xi,i=0,1,……,m}上给出，这就是科学实验中常见到的实验数据{（xi，yi）,i=0,1,……,m}的曲线拟合，这里yi=f(xi)(i=0,1,……,m),要求一个函数y=s*(

7、x)与所给的数据{（xi，yi）,i=0,1,……,m}拟合，若记误差设φ0(x)，φ1(x)，……φn(x)是C[a,b]上线性无关函数族，在φ=span{φ0(x)，φ1(x)，……φn(x)}中找一个函数s*(x)，使误差平方和这就是一般的最小二乘逼近，用几何语言说，就称为曲线拟合的最小二乘法。已知数据如表2-1x1.22.84.35.46.87.9y2.111.528.141.972.391.4根据表2-1所给的数据，用0rigin软件得出它的散点图2-1图2-1根据散点图2-1可看出该图形可能为多项式函数，所以不妨设该方程式为y=a*x2+b*x+c。