计算机数值方法第三章 插值法与最小二乘法.pptx

计算机数值方法第三章 插值法与最小二乘法.pptx

  1. 1、本文档共145页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
计算机数值方法第三章 插值法与最小二乘法.pptx

延安大学 数学与计算机科学学院;本章要点:;第3章 插值法 ;§1 插值问题;其插值函数的图象如下图;;线性插值与二次插值 ;问题 求作一次式 ,使满足条件 从几何图形上看, 表示过两点 的直线。;过两点的直线 ,可表示为如下点斜式:;其中,;例1: 已知 , ,求 代入点斜式插值多项式 得 y=10.71428 精确值为 10.723805, 故这个结果有3位有效数字。; 问题 求作二次式 ,使满足条件 二次插值的几何解释是用通过三个点 的抛物线来近似考察曲线,故又称为拋物插值。;类似于线性插值,构造基函数,要求满足下式:;称;代数多项式插值的存在唯一性 ; 这样的多项式是否存在并且唯一呢? 根据插值条件,代数多项式(3―1)中的各个系数a0,a1,…,an应满足下列n+1阶线性方程组;线性方程组的系数行列式为范德蒙特行列式;代数多项式的余项 ; 则对插值区间上的任何x,都存在ξ∈(a,b),使得; (2)若f(x)为次数不超过n的多项式,那么以n+1个点为基点的插值多项式就一定是其本身,即Pn(x)≡f(x)。这是因为此时 Rn(x)=0。 (3)从余项Rn(x)中的ω(n+1)(x)知,当点x位于x0,x1,…,xn的中部时,|ωn+1(x)|比较小,精度要高一些,而位于两端时,精度要差一些;若x位于x0,x1,…,xn的外部,一般称为外插(或外推),此时精度一般不理想,使用时必须注意。;且在不同的基底下有不同的形式;;n+1次多项式;令;其中;Lagrange线性插值基函数为;例1:;且;例2.;Lagrange线性插值多项式为;拉格朗日插值多项式的缺点:; 3 分段低次插值法;;也称折线插值,如下图,;;二、分段二次Lagrange插值;;一般;;2. 分段二次插值的误差估计;分段低次Lagrange插值的特点;因为任何一个n次多项式都可以表示成;线性无关,;有;再继续下去待定系数的形式将更复杂;定义1.;显然;差商的主要性质:;(2) 对称性,即任意调换节点的次序,差商的值不变;;称;由插值多项式的唯一性,Newton基本插值公式的余项为;则Newton均差插值公式;定义2.;依此类推;可以证明;;在等距节点的前提下,差商与差分有如下关系;依此类推;由??商与向前差分的关系;则均差插值公式;插值余项为;插值余项为;例 已知数据表; 例2 分别作出 f(x)=x2+x+1 在区间[1,4]的前差、后差表,其中h=1。 解:前差表、后差表如下: ;前 差 表; 例3 给出正弦函数sinx由x=0.4到0.7的值(h=0.1),试分别 用牛顿前差和后差公式计算sin0.57891的近似值。 解 作差分表:; 利用牛顿前差公式; 利用牛顿后差公式 ;;;;;希望插值系数与Lagrange插值一样简单;;可得;;;二、两点三次Hermite插值的余项;构造辅助函数;即;例1.; 作为多项式插值,三次已是较高的次数,次数再 高就有可能不稳定,因此,对有n+1个节点的插值问题, 我们可以使用分段两点三次Hermite插值;三、分段两点三次Hermite插值;其中;我们称; 6 三次样条插值;一、三次样条插值函数;------(1);二、三次样条插值多项式;------(3);------(4);少两个条件;也就是所谓的边界条件:;;------(9);;由条件;;;;由(11)式,可知;;(19)式与(14)一样,都是三对角方程组,并且都严格对角占优;例1. 对于给定的节点及函数值;由(19)式得基本方程组;定理 . ;插值法及主要特性:;问题提出;第七节 数据拟合的最小二乘法;;;; 已知函数y=f(x)的数据点( xi , yi )(i=0,…,m),

文档评论(0)

czy2014 + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档