数值分析教学课件:6曲线拟合_第1页
数值分析教学课件:6曲线拟合_第2页
数值分析教学课件:6曲线拟合_第3页
数值分析教学课件:6曲线拟合_第4页
数值分析教学课件:6曲线拟合_第5页
已阅读5页,还剩31页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

第六章曲线拟合的最小二乘

/函数平方逼近初步NumericalValueAnalysis一.问题的提出插值法是使用插值多项式来逼近未知或复杂函数的,它要求插值函数与被插函数在插值节点上函数值相同,而在其他点上没有要求。在非插值节点上有时函数值会相差很大。若要求在被插函数的定义区间上都有较好的近似,就是最佳逼近问题。必须找到一种度量标准来衡量什么是最佳逼近.

最佳逼近是在函数空间M中选

P(x)满足

但由于绝对值函数不宜进行分析运算,常将上式化为来讨论,于是最佳逼近问题变为最佳平方逼近问题这即为连续函数的最佳平方逼近.对于离散的问题,最佳平方逼近问题为:就是常说的曲线拟合的最小二乘法.

最佳逼近二.预备知识内积:常采用的内积与范数实例:考察某种纤维的强度与其拉伸倍数的关系,下表是实际测定的24个纤维样品的强度与相应的拉伸倍数是记录:一.实例讲解6.2数据拟合(最小二乘法)纤维强度随拉伸倍数增加而增加并且24个点大致分布在一条直线附近---------(1)必须找到一种度量标准来衡量什么曲线最接近所有数据点.二、问题的提法定义平方误差(偏差平方和):我们选取的度量标准是---------(2)---------(3)三、法方程组由可知因此可假设因此求最小二乘解转化为二次函数由多元函数取极值的必要条件得即---------(4)即引入记号则由内积的概念可知---------(5)---------(6)显然内积满足交换律方程组(4)便可化为---------(7)将其表示成矩阵形式-----(8)并且其系数矩阵为对称阵.根据Cramer法则,法方程组有唯一解即是的最小值所以因此作为一种简单的情况,基函数之间的内积为平方误差例1.回到本节开始的实例,从散点图可以看出纤维强度和拉伸倍数之间近似与线性关系故可选取线性函数为拟合函数,其基函数为建立法方程组根据内积公式,可得法方程组为解得平方误差为拟合曲线与散点的关系如右图:四、加权最小二乘法各点的重要性可能是不一样的权:即权重或者密度,统称为权系数.

定义加权平方误差为-----(9)使得由多元函数取极值的必要条件得即引入记号定义加权内积-----(10)矩阵形式(法方程组)为方程组(10)式化为-----(11)---(12)平方误差为作为特殊情形,用多项式作拟合函数的法方程组为-----(13)

6.3连续函数的最佳平方逼近1.最佳平方逼近问题-----(14)2.解法(法方程)-----(15)最小二乘法方法评注曲线拟和的最小二乘法是实验数据处理的常用方法。最佳逼近可以在一个区间上比较均匀的逼近函数且具有方法简单易行,实效性大,应用广泛等特点。但当法方程组阶数较高时,往往出现病态。因此必须谨慎对待和加以巧妙处理。有效方法之一是引入正交多项式以改善其病态性。指数模型和双曲线模型-----线性化拟合超定方程组的最小二乘解Seeyounextchapter!《计

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论