最小二乘法和拉格朗日插值法在数据测试处理中的综合应用
- 格式:pdf
- 大小:590.30 KB
- 文档页数:12
最小二乘法在数据处理中的应用嘿,朋友!想象一下这样一个场景:你正在为一个科学实验收集数据,一堆数字摆在你面前,就像一群调皮的小精灵,让你眼花缭乱,不知所措。
这时候,“救星”出现了,那就是最小二乘法!比如说,有个叫小李的科研工作者,正为他的实验数据愁眉苦脸。
他的实验是研究植物在不同光照条件下的生长速度。
经过一段时间的辛苦观察和记录,他得到了一堆光照时长和植物生长高度的数据。
可这些数据杂乱无章,怎么从中找出规律呢?这时候,最小二乘法就大显身手啦!它就像一个神奇的魔法棒,能把这些看似混乱的数据变得有条有理。
最小二乘法到底是怎么施展魔法的呢?简单来说,它就是要找到一条最合适的线或者曲线,来尽可能地靠近这些数据点。
这就好比你要穿过一片树林,找到一条最顺畅的小路,让你能以最省力的方式通过。
假设小李的数据点分布得比较接近一条直线,那最小二乘法就能算出这条直线的方程。
它会考虑每个数据点与这条假设直线的距离,然后通过一系列巧妙的计算,让这些距离的平方和最小。
这是不是很神奇?想象一下,如果没有最小二乘法,小李就得靠自己的眼睛和感觉去估摸数据的规律,那得多不靠谱啊!就像闭着眼睛在黑屋子里找东西,全凭运气。
在实际生活中,最小二乘法的应用可广泛啦!不只是科研领域,经济领域也少不了它。
比如说,预测股票价格的走势,分析市场的需求和供应关系等等。
它就像一个聪明的参谋,为决策者提供可靠的依据。
再比如,在工程领域,测量建筑物的变形、评估机器的性能,最小二乘法都能发挥巨大的作用。
它能帮助工程师们更准确地了解物体的状态,提前发现潜在的问题,避免出现大的失误。
你可能会想,这么厉害的方法,是不是很难掌握呢?其实不然!只要你有一些基本的数学知识,再加上一点耐心和细心,就能理解和运用它。
总之,最小二乘法在数据处理中简直就是一把“万能钥匙”,能打开数据背后隐藏的秘密之门,让我们更加清晰地看到事物的本质和规律。
它就像一位默默无闻的英雄,在幕后为我们的科学研究、经济决策和工程建设等众多领域提供着强大的支持和帮助。
函数逼近的几种算法及其应用汇总函数逼近是数值计算中非常重要的技术之一,它主要用于用已知函数逼近未知函数,从而得到未知函数的一些近似值。
在实际应用中,函数逼近广泛用于数据拟合、插值、信号处理、图像处理等领域。
下面将介绍几种常用的函数逼近算法及其应用。
1. 最小二乘法(Least Square Method)最小二乘法将函数逼近问题转化为最小化离散数据与拟合函数之间的残差平方和的问题。
它在数据拟合和插值中应用广泛。
例如,最小二乘法可以用于拟合数据点,找出最佳拟合曲线;也可以用于信号处理中的滤波器设计。
2. 插值法(Interpolation)插值法旨在通过已知数据点之间的连线或曲线,来逼近未知函数在这些数据点上的取值。
常见的插值方法有拉格朗日插值、牛顿插值和分段线性插值等。
插值法在图像处理中广泛应用,例如可以通过已知的像素点来重构图像,提高图像的质量和分辨率。
3. 最小二乘曲线拟合(Least Square Curve Fitting)最小二乘曲线拟合是一种将渐近函数与离散数据拟合的方法,常见的函数包括多项式、指数函数、对数函数等。
最小二乘曲线拟合可以在一定程度上逼近原始数据,从而得到曲线的一些参数。
这种方法在数据分析和统计学中经常使用,在实际应用中可以拟合出模型参数,从而做出预测。
4. 正交多项式逼近(Orthogonal Polynomial Approximation)正交多项式逼近是一种通过正交多项式来逼近未知函数的方法。
正交多项式具有良好的性质,例如正交性和递推关系,因此可以用于高效地逼近函数。
常见的正交多项式包括勒让德多项式、拉盖尔多项式和切比雪夫多项式等。
正交多项式逼近广泛应用于数值计算和信号处理中,例如用于图像压缩和数据压缩。
5. 插值样条曲线(Interpolating Spline)插值样条曲线是将多个局部的多项式插值片段拼接在一起,从而逼近未知函数的方法。
插值样条曲线在实现光滑拟合的同时,还能逼近离散数据点。
最小二乘法及其在数据拟合中的应用在现代科学和工程领域,数据拟合是一项重要的任务。
通过拟合数据,我们可以找到数据背后的规律,并用数学模型来描述这些规律。
而最小二乘法是一种常用的数据拟合方法,它可以帮助我们找到最佳的拟合曲线或者函数。
最小二乘法的基本原理是通过最小化误差的平方和来拟合数据。
在数据拟合中,我们通常会有一组离散的数据点,我们的目标是找到一条曲线或者函数,使得这些数据点到曲线的距离最小。
而这个距离可以通过计算每个数据点到曲线的垂直距离来表示。
假设我们有一组数据点(x1, y1), (x2, y2), ..., (xn, yn),我们要找到一个函数f(x)来拟合这些数据点。
最小二乘法的思想是,我们要找到一个函数f(x),使得数据点到函数的垂直距离的平方和最小。
换句话说,我们要找到一个函数f(x),使得Σ(yi - f(xi))^2最小。
为了实现最小二乘法,我们需要选择一个合适的函数形式来拟合数据。
常见的函数形式包括线性函数、多项式函数、指数函数等。
以线性函数为例,我们要找到一个直线y = ax + b来拟合数据。
通过最小二乘法,我们可以得到最佳的a和b的取值,使得数据点到直线的垂直距离的平方和最小。
最小二乘法的求解过程可以通过数学推导得到闭式解,也可以通过数值优化算法来求解。
在实际应用中,我们通常会使用计算机来进行求解。
计算机可以通过迭代的方式,逐步调整函数的参数,使得误差平方和不断减小,最终找到最佳的拟合曲线或者函数。
最小二乘法在数据拟合中有着广泛的应用。
它可以用于拟合实验数据,找到实验结果背后的数学模型。
例如,科学家可以通过最小二乘法来拟合实验数据,找到物理定律的数学表达式。
最小二乘法还可以用于拟合观测数据,找到数据背后的规律。
例如,经济学家可以通过最小二乘法来拟合经济数据,找到经济模型的参数。
除了数据拟合,最小二乘法还有其他的应用。
例如,在信号处理中,最小二乘法可以用于滤波和降噪。
通过最小二乘法,我们可以找到一个滤波器或者降噪算法,使得信号的噪声被最小化。
最小二乘法及其应用1. 引言最小二乘法在19世纪初发明后,很快得到欧洲一些国家的天文学家和测地学家的广泛关注。
据不完全统计,自1805年至1864年的60年间,有关最小二乘法的研究论文达256篇,一些百科全书包括1837年出版的大不列颠百科全书第7版,亦收入有关方法的介绍。
同时,误差的分布是“正态”的,也立刻得到天文学家的关注及大量经验的支持。
如贝塞尔( F. W. Bessel, 1784—1846)对几百颗星球作了三组观测,并比较了按照正态规律在给定范围内的理论误差值和实际值,对比表明它们非常接近一致。
拉普拉斯在1810年也给出了正态规律的一个新的理论推导并写入其《分析概论》中。
正态分布作为一种统计模型,在19世纪极为流行,一些学者甚至把19世纪的数理统计学称为正态分布的统治时代。
在其影响下,最小二乘法也脱出测量数据意义之外而发展成为一个包罗极大,应用及其广泛的统计模型。
到20世纪正态小样本理论充分发展后,高斯研究成果的影响更加显著。
最小二乘法不仅是19世纪最重要的统计方法,而且还可以称为数理统计学之灵魂。
相关回归分析、方差分析和线性模型理论等数理统计学的几大分支都以最小二乘法为理论基础。
正如美国统计学家斯蒂格勒( S. M. Stigler)所说,“最小二乘法之于数理统计学犹如微积分之于数学”。
最小二乘法是参数回归的最基本得方法所以研究最小二乘法原理及其应用对于统计的学习有很重要的意义。
2. 最小二乘法所谓最小二乘法就是:选择参数10,b b ,使得全部观测的残差平方和最小. 用数学公式表示为:21022)()(m ini i i i ix b b Y Y Y e--=-=∑∑∑∧为了说明这个方法,先解释一下最小二乘原理,以一元线性回归方程为例.i i i x B B Y μ++=10 (一元线性回归方程)由于总体回归方程不能进行参数估计,我们只能对样本回归函数来估计即:i i i e x b b Y ++=10)...2,1(n i =从上面的公式可以看出:残差i e 是i Y 的真实值与估计值之差,估计总体回归函数最优方法是,选择10,B B 的估计量10,b b ,使得残差i e 尽可能的小.总之,最小二乘原理就是选择样本回归函数使得所有Y 的估计值与真实值差的平方和为最小,这种确定10,b b 的方法叫做最小二乘法。
函数逼近的几种算法及其应用函数逼近是数值计算中的一种重要技术,用于在给定的函数空间中找到与目标函数最相近的函数。
函数逼近算法可以在不知道目标函数解析表达式的情况下,通过对给定数据进行处理来逼近目标函数的结果。
这篇文章将介绍几种常见的函数逼近算法及其应用。
1.多项式逼近:多项式逼近是一种利用多项式函数逼近目标函数的方法。
多项式逼近算法有很多种,常见的有最小二乘法、拉格朗日插值法和牛顿插值法等。
多项式逼近广泛应用于数据拟合、信号处理和图像处理等领域。
最小二乘法是一种通过最小化实际观测值与多项式模型之间的差异来确定多项式系数的方法。
最小二乘法可以用于拟合非线性和线性函数。
拉格朗日插值法和牛顿插值法是通过插值多项式来逼近目标函数的方法,可以用于填充缺失数据或者生成曲线过程中的中间点。
2.三角函数逼近:三角函数逼近是一种利用三角函数来逼近目标函数的方法。
三角函数逼近算法有傅里叶级数逼近和小波变换等。
傅里叶级数逼近是一种利用三角函数的线性组合来逼近目标函数的方法。
这种方法广泛应用于信号处理、图像处理和数学建模等领域。
小波变换是一种通过特定的基函数来逼近目标函数的方法。
小波变换可以用于信号去噪、图像压缩和模式识别等应用。
3.插值逼近:插值逼近是一种通过已知数据点在给定区间内的函数值来确定目标函数的方法。
常见的插值逼近方法有拉格朗日插值法、牛顿插值法和差值多项式法等。
插值逼近广泛应用于任何需要通过已知数据点来逼近目标函数的领域。
在实际应用中,函数逼近常用于数据分析和模型构建。
例如,在金融领域,函数逼近可以用于确定股票价格走势的模型和预测。
在工程领域,函数逼近可以用于建立复杂系统的模型和优化控制。
在计算机图形学领域,函数逼近可以用于生成真实感图像和动画。
总结起来,函数逼近是一种重要的数值计算技术,有多种算法可供选择。
多项式逼近、三角函数逼近和插值逼近是常见的函数逼近算法。
函数逼近广泛应用于数据分析、模型构建和优化控制等领域,对于解决实际问题具有重要作用。
最小二乘法在数据拟合中的应用最小二乘法是数学中的一种常见方法,用于在一组数据中找到最符合数据特征的函数模型。
在数据分析和拟合中,使用最小二乘法可以对实验数据进行比较准确的模型推导和预测。
最小二乘法的原理最小二乘法的核心思想是通过对目标函数中的平方误差求和,并将误差平方和最小化来确定函数参数值。
简言之,就是用一个函数去拟合一些数据点,找出最能够符合这些数据点的函数方程,从而得到预测或分析的标准。
具体而言,最小二乘法会先提供一组有n个坐标的点(x1, y1), (x2, y2), ..., (xn, yn),然后根据这些数据点来求出一个一定形式的函数y = f(x),使得y值与每个点的目标值yr之间的误差平方之和最小。
求这个函数就是求最小二乘的函数方程。
应用最小二乘法的过程用最小二乘法对数据进行拟合的步骤如下:1. 收集实验数据,并把数据图表显示出来;2. 根据数据情况选择函数模型;3. 对选择的函数模型变量进行求解;4. 通过计算每组实验数据与模型曲线之间的距离平方和,得到拟合函数的误差;5. 对误差函数求导取极小值,从而确定拟合函数的系数和截距;6. 最后将得到的拟合函数与实验数据绘制到同一张图表上,检验拟合效果。
实际应用在实际的科学研究和工程应用中,最小二乘法在数据分析和拟合中被广泛应用。
例如,最小二乘法可以用于分析物理实验数据,以推导出实验工作曲线;在经济学中,最小二乘法可以用于分析价格和销售数据之间的关系,以预测市场走势;在金融学中,最小二乘法可以应用于证券交易中,以实现资产组合优化。
最小二乘法还可以应用于数字信号处理和机器学习等领域。
例如,在数字信号处理中,最小二乘法可以用于降噪和滤波信号;在机器学习中,最小二乘法可以用于模型训练和预测。
总结从原理到实际应用,最小二乘法在科学研究和工程领域中具有广泛的应用。
这种方法可以对实验数据进行准确的模型推导和预测,帮助科学家和工程师更好地理解数据,并从中获得更多信息。
最小二乘法在计算机算法中的应用分析随着计算机科学的发展,越来越多的数学算法被应用于计算机编程中,提高了编程的准确性和效率。
其中,最小二乘法是一种常用的数学算法,可以在多个领域中被应用。
本文将分析最小二乘法在计算机算法中的具体应用,并探讨其优缺点。
1. 最小二乘法的基本概念最小二乘法是一种数学优化方法,用于通过最小化误差平方和来拟合数据。
在最小二乘法中,误差是指观测值和拟合值之间的差距。
其基本公式为:min Σ(y - f(x))^2其中,y为观测值,f(x)为拟合值。
最小二乘法可以求出最优的拟合函数,使得误差平方和最小。
2. 最小二乘法在曲线拟合中的应用最小二乘法在计算机算法中最常见的应用是曲线拟合。
在曲线拟合中,需要找到一条最能代表观测数据的曲线,这就需要用到最小二乘法。
最小二乘法可以拟合多项式、正弦曲线、指数曲线等多种类型的曲线。
例如,想要通过一组x和y的观测值来拟合一条多项式曲线,就可以用最小二乘法。
首先,需要选择多项式的阶数,比如2、3、4等。
然后,通过最小二乘法求得多项式系数,即可得到拟合曲线。
3. 最小二乘法在回归分析中的应用回归分析是一种统计学方法,用于分析两个或多个变量之间的关系。
最小二乘法在回归分析中是一种常用的方法,用于对变量之间的关系进行建模。
例如,考虑一个简单的线性回归模型:y = a + bx,其中y是被解释变量,x是解释变量,a和b是常数。
可以用最小二乘法计算出最优的a和b值,使得拟合函数最能代表数据。
最小二乘法可以拟合不同类型的回归模型,比如一个单一的解释变量、多个解释变量、定性变量、非线性关系等。
在实际应用中,回归分析可以用于预测和控制因素,比如销售量、股票价格等。
4. 最小二乘法的优缺点最小二乘法作为一种常用的数学算法,具有一些优点和缺点。
优点:最小二乘法易于使用,且可以用于拟合不同类型的数据,包括线性和非线性数据。
其算法简单易懂,而且具有广泛的应用领域,比如机器学习、图像处理、信号处理等。
最小二乘法在数据拟合中的应用数据拟合是科学研究、工程设计等领域中常见的问题,它是指根据已知的数据,通过一定的方法,建立一个能够描述数据特征的数学模型。
通常情况下,这些数据点之间存在着一些误差,因此拟合出的模型不可能完全精确地描述实际情况。
在这种情况下,最小二乘法就成了处理数据拟合问题的重要工具之一。
最小二乘法是一种数学优化方法,它的主要思想是寻求一种数学函数,使得该函数与一组给定的数据点之间的误差平方和最小。
这个方法的应用范围非常广泛,不仅可以用于数据拟合问题,还可以用于解决图像处理、信号处理、回归分析等许多实际问题。
需要注意的是,最小二乘法只适用于某些类型的函数拟合问题。
具体来说,当我们希望拟合的函数可以表示为线性组合形式时,最小二乘法就成为了一种有效的求解方法。
在最小二乘法中,我们通常会选择一种函数形式,然后调整函数中的参数,使得该函数与所有数据点之间的距离最小。
为了达到这个目标,我们通常会计算每个数据点到拟合函数的垂直距离,即残差,然后将所有的残差平方相加,得到误差平方和。
最终我们需要调整函数中的参数,使得误差平方和最小。
例如,我们可以尝试使用一个一阶多项式函数拟合一组数据点:$y = ax + b$对于每一个数据点 $i$,我们可以计算该数据点到拟合的直线上的垂直距离为:$e_i = y_i - (ax_i + b)$我们可以将所有的残差平方相加,得到误差平方和:$S(a,b) = \sum\limits_{i}^n e_i^2 = \sum\limits_{i}^n (y_i - ax_i - b)^2$现在的问题是如何求解 $a$ 和 $b$,使得误差平方和$S(a,b)$ 最小。
这时候,最小二乘法就派上用场了。
我们可以使用偏导数的方法,求得误差平方和对 $a$ 和 $b$ 的导数,然后令它们等于零,解出最小二乘估计值。
具体来说,我们可以得到以下两个方程:$\dfrac{\partial S}{\partial a} = -2\sum\limits_{i}^n (y_i - ax_i - b)x_i = 0$$\dfrac{\partial S}{\partial b} = -2\sum\limits_{i}^n (y_i - ax_i - b) = 0$解这两个方程得到:$a = \dfrac{\sum\limits_{i}^n x_i y_i -\dfrac{1}{n}(\sum\limits_{i}^n x_i)(\sum\limits_{i}^ny_i)}{\sum\limits_{i}^n x_i^2 - \dfrac{1}{n}(\sum\limits_{i}^nx_i)^2}$$b = \bar{y} - a\bar{x}$其中,$\bar{x}$ 和 $\bar{y}$ 分别是 $x$ 和 $y$ 的平均值。
最小二乘法在实验数据处理中的应用最小二乘法,简称最小二乘,是一种广泛用于拟合的数学方法。
它通过最小化误差的平方和来得到最优的参数估计值。
在实验数据处理中,最小二乘法可以用来拟合实验数据。
例如,在物理实验中,我们常常需要测量一些物理量,比如电阻、电动势、势能等。
这些物理量往往是经过实验测量得到的,并且存在一定的误差。
为了得到更加精确的结果,我们可以使用最小二乘法来拟合这些数据。
那么,最小二乘法是如何进行拟合的呢?我们可以用一个函数来描述实验数据的规律,比如说:y = ax + b其中,a和b是我们要求的参数,y是测量得到的数据,x是一些已知的量。
我们的目标就是通过最小二乘法来求出a和b的最优估计值。
为了求出a和b的最优估计值,我们需要定义一个误差函数,用来衡量预测值与真实值之间的差异。
这个误差函数通常是平方误差函数,即:E = Σ(y - y')^2其中,y是测量得到的数据,y'是通过拟合函数得到的预测值,Σ表示将所有的误差平方值加起来。
我们的目标就是让这个误差函数的值最小。
为了使误差函数的值最小,我们需要求出a和b的最优估计值。
这就需要使用一种叫做梯度下降法的优化算法。
梯度下降法的原理是:我们首先随机给出a和b的初始值,然后不断调整a和b的值,使得误差函数的值不断减小。
这样,我们就能求出a和b的最优估计值。
最小二乘法在实验数据处理中的应用非常广泛。
它可以用来拟合各种各样的数据,包括物理量、化学量、生物量等。
它的优势在于,它能够提供精确的参数估计值,并且易于实现。
总之,最小二乘法是一种非常有用的数学方法,在实验数据处理中有着广泛的应用。
通过对最小二乘法的了解和掌握,我们能够更好地处理实验数据,为科学研究提供有力的支持。
在结束这篇文章之前,我们再来总结一下最小二乘法在实验数据处理中的应用。
首先,最小二乘法是一种用于拟合数据的数学方法。
它通过最小化误差的平方和来得到最优的参数估计值。
其次,最小二乘法在实验数据处理中有着广泛的应用。
第一题:曲线拟合最小二乘法和切比雪夫的相同和不同,以及适用的场合背景及意义:在很多日常生活以及科研活动中,我们需要对一些离散的点集进行拟合,使得拟合的曲线尽量多的穿过所给出的离散点,并且误差小。
从而通过拟合的函数,找出离散点的规律,以此进行进一步的研究。
下面,就最小二乘法和切比雪夫两种拟合方法进行研究和分析。
1、最小二乘法它的标准是,所求得的拟合函数*()y S x =与给出的实际离散点{(,),0,1,,i i x y i m = 之间的误差平方和最小。
公式为:*()(1,2,,)i i i S x y i m δ=-=01(,,,)T m σδδδ=22222()0||||[*()]min[()]m mmii i iiS x i i i S x y S x y ϕσδ∈=====-=-∑∑∑其中ϕ是规定区间上的线性无关函数族,01{,,,}m span ϕϕϕϕ= 。
为了使问题提法更具一般性,在各自的离散点的区间中添加权函数()0x ω≥以表示各个离散点数据的比重不同。
要想求出函数*()S x ,就要求出其各阶系数,转而变成求多元函数极小点***01,,,n a a a2010(,,,)()[()()]mnn i j j i i i j I a a a x a x f x ωϕ===-∑∑其中:()j x ϕ取21,,,,n x x x的问题。
为了求取极值,其必要条件为2002()[()()]()0,0,1,,.m ni j j i i k i i j k Ix a x f x x k n a ωϕϕ==∂=-==∂∑∑ 简化上式可得到矩阵形式Ga d =其中01(,,,)T n a a a a = ,01(,,,)T n d d d d = ,00010(,)(,)(,)(,)n n n G ϕϕϕϕϕϕϕϕ⎛⎫⎪= ⎪ ⎪⎝⎭要想使所求极值有唯一解,就要求G 非奇异。
又因()j x ϕ的组所组成向量为非奇异,则G 为非奇异,故而存在唯一的解*,0,1,.k k a a k n == 使得*()S x 为所求最优解。
计算方法
题目最小二乘法和拉格朗日插值法在数据测试处理中的综合应用
学院机电工程学院
专业控制工程
姓名徐进
学号 1504122120
最小二乘法和拉格朗日插值法在数据测试处理中的综合应用
一、应用背景概述
在工程实践和科学研究中,常常需要对系统进行某些性能或现象的测试研究,来探求其中的某些内在规律。
一般而言,采用实验的方法获得实验数据,然后经过处理实验数据来表征出目标变量间的相互关系。
在实际的数据测试处理中,不仅需要利用数值计算方法来探求出最终的表征函数关系,而且对于测试过程中的粗大误差数据,可以通过分析观察被测数据与表征目标函数曲线的吻合情况,可以判断出测试过程中含有粗大误差的数据,并加以剔除,以免对后续系统的研究产生干扰和不必要的麻烦。
数据处理和分析中应用较多的是最小二乘法和拉格朗日插值法,这两种数值分析方法是运算简便且应用广泛的数据测试处理表示方法。
本次应用实践,充分考虑两种方法的优势,通过两种方法得到两种拟合曲线,然后通过对两种拟合曲线函数多项式系数相应进行平均值计算,得到整合后的新的多项式。
二、应用方法原理概述
①粗大误差
粗大误差又称疏忽误差或过失误差,它是由于技术不成熟,测量时不小心或外界的突然干扰(例如突然振动、仪器电源电压的突然变化)等原因造成的。
含有粗大误差的测量数据,常比正常数据相差较大(过大或过小)。
当对某一量值作多次独立的等精度重复测量,如其中个别或少数数据明显地偏大或偏小时,则可怀疑数据中含有粗大
误差。
本次应用中,通过被测数据与拟合函数图形进行对比的方法,剔除偏离较大的粗大误差。
②最小二乘法和拉格朗日插值法综合应用原理
关于最小二乘法和拉格朗日插值法各自的详细介绍,在这里就不再赘述。
本次应用,主要通过利用最小二乘法和拉格朗日插值法在数值分析中的优点,对其进行融合,得到一种改进的探求拟合目标函数的方法。
其应用原理如下:
例如,利用最小二乘法拟合得到的目标函数的多项式为: y=a0+a1x+a2x2+…+a n x n
利用拉格朗日插值法得到的目标函数多项式为
Y=b0+b1x+b2x2+…+b n x n
则结合方法得到的多项式为:
简单来说,就是将最小二乘法得到的拟合多项式各自的系数和拉格朗日插值法得到的拟合多项式各自的系数分别相应进行算术平均值计算,得到的各项新系数,构成新的最终拟合目标函数多项式。
三、MATLAB实践应用
本次是以某控制系统中的光电传感装置的光强度——电压为测试参数对,记录了9组测试数据,如下表所示:(数据已归一化处理)
一般在应用数值分析方法探求测试数据拟合曲线之前,应该对测试数据的走势有一个定性的判断,以便后续准确的设置拟合函数多项式的次数,因此先对测试样本数据进行离散图显示观察,其MATLAB 程序及结果图为:
原始测试数据坐标绘点显示程序:
x=[7.378,7.496,7.601,7.696,7.783,7.836,7.937,8.006,8.134];
y=[2.05,2.65,3.16,3.60,3.99,4.38,4.72,5.01,5.83];
plot(x,y,'*')
图一原始样本数据图形显示结果
通过对原始测试数据的直观显示,可以发现数据是具有一次性函数的特点,因此后续的拟合函数都应设置为一次函数型多项式。
1.最小二乘法剔除测试数据中的粗大误差数据
由于测试数据中的粗大误差对后续的分析研究干扰较大,容易导致错误结论,因此首先应剔除掉粗大误差的数据
MATLAB源程序如下:
最小二乘法函数定义m文件为:
function c=lspoly(x,y,m)
n=length(x);
b=zeros(1:m+1);
f=zeros(n:m+1);
for k=1:m+1
f(:,k)=x'.^(k-1);
end
a=f'*f;
b=f'*y';
c=a\b;
c=flipud(c);
最小二乘法算法运行程序:
lspoly([7.378,7.496,7.601,7.696,7.783,7.836,7.937,8.006,8.134],[2. 05,2.65,3.16,3.60,3.99,4.38,4.72,5.01,5.83],1)
ans = 4.8686
-33.8629
图形对比结果为;
图二最小二乘法处理结果图(含粗大误差数据)可以看出最后一组数据偏离较大,可判断为粗大误差,予以剔除
2.最小二乘法和拉格朗日插值法综合数据处理
由前面的粗大误差处理后,得到8组正常的测试数据样本,新的表格数据如下:
同理应用最小二乘法MATLAB程序算出8组数据下的拟合函数表达式,程序如下:
lspoly([7.378,7.496,7.601,7.696,7.783,7.836,7.937,8.006],[2.05,2.6 5,3.16,3.60,3.99,4.38,4.72,5.01],1)
ans =4.7522
-32.9759
因此可得由最小二乘法算出的拟合函数表达式为;y=4.7522x-32.9759。
再者,通过拉格朗日插值法对8组数据进行拟合函数计算,由于测试数据实际较符合一次函数模型,因此拟合函数设定为一次函数。
通过对8组数据两两为一组进行拉格朗日插值法计算出28个拟合函数表达式,此过程也可通过MATLAB计算出,拉格朗日MATLAB 函数M文件程序为:
function [c,l]=lagran(x,y)
w=length(x);
n=w-1;
l=zeros(w,w);
for k=1:n+1
v=1;
for j=1:n+1
if k~=j
v=conv(v,poly(x(j)))/(x(k)-x(j));
end
end
l(k,:)=v;
end
c=y*l;
计算结果如下:
Y=5.08x-35.47 y=4.98x-34.69 y=4.87x-33.91 y=4.8x-33.39 y=4.77x-33.19 y=4.79x-33.29 y=4.71x-32.73 y=4.86x-33.76
y=4.75x-32.96 y=4.67x-32.30 y=4.71x-32.69 y=4.69x-32.54
y=4.63x-32.04 y=4.63x-32.04 y=4.56x-31.5 y=4.66x-32.23
y=4.7x-32.17 y=4.57x-31.56 y=4.48x-33.95 y=4.67x-32.34
y=4.64x-32.17 y=4.54x-31.4 y=4.88x-33.95 y=4.74x-32.9
y=4.57x-31.61 y=4.59x-31.75 y=4.41x-30.26 y=4.2x-28.64 为了保证数据处理精度,对以上28个拟合函数各项系数进行均值化处理,得出一个总的均值化后的拟合函数,处理如下:
=4.68
K=5.08+4.98+⋯+4.41+4.2
28
b=−(32.47+34.69+⋯+30.26+28.64)
=-32.47
28
因此拉格朗日插值法拟合函数为y=4.68x-32.47
其matlab程序及图形为
x=[7.378,7.496,7.601,7.696,7.783,7.836,7.937,8.006];
y=[2.05,2.65,3.16,3.60,3.99,4.38,4.72,5.01];
t=7:0.1:9;
z=4.68*t-32.47;
plot(x,y,'*',t,z);grid
图三拉格朗日插值法处理结果图
最后将最小二乘法和拉格朗日插值法算得的拟合函数各项系数进行均值化后,得到本次数据处理最终的拟合函数,处理如下:
最小二乘法拟合函数y=4.7522x-32.9759
拉格朗日插值法拟合函数y=4.68x-32.47
K’=4.7522+4.68
=4.7161
2
b’=−32.9759−32.47
=-32.72295
2
因此最终函数多项式系数均值化处理后数据拟合函数为
y=4.7161x-32.72295
MATLAB程序及图形为:
x=[7.378,7.496,7.601,7.696,7.783,7.836,7.937,8.006];
y=[2.05,2.65,3.16,3.60,3.99,4.38,4.72,5.01];
t=7:0.1:9;
z=4.7161*t-32.72295;
plot(x,y,'*',t,z);grid
图四最小二乘法和拉格朗日插值法综合处理结果图
从结果可以看出,通过最小二乘法和拉格朗日插值法综合处理得到的拟合函数与实测数据吻合度很高,比单独使用一种方法处理测试数据,其效果和准确性要提高很多,具有较大的优越性。
四、小结
本次实践应用主要用于数据的测试处理,通过最小二乘法这种数值分析方法来剔除数据中较明显的粗大误差,将最小二乘法和拉格朗日插值法两者结合,优势叠加,形成一种改进的探求拟合函数的方法并加以应用于光电传感系统的数据测试处理,其效果明显优于单独使用一种数值分析的方法。
数值分析方法在工程上具有很大的应用价值,不仅在数据测试处理中,也在我正在学习的图像处理方面也有很多应用,比如中值滤波和边沿检测等处理环节,其用到的数值分析方法更为复杂。