- 1、本文档共51页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
数值分析(广东工业大学课件).ppt
* * * * * * * * * * * * * * * * * * * * * * 分段低次插值只能保证各段曲线在节点上的 连续性,而不能保证整条曲线在节点上的充 分光滑性(即一阶二阶导数存在)。 §2.7 样条插值 一、样条插值的概念 注:三次样条与分段Hermite 插值的根本区别在于S(x)自身光滑,不需要知道f 的导数值(除了在2个端点可能需要);而Hermite插值依赖于f 在所有插值点的导数值。 二、三次样条插值函数的建立 三、误差界与收敛性 * * * * * * * * * * * * 第二章 插值 §2.4 差分与等距节点的Newton插值 上节讨论任意分布节点的插值公式,应用时常碰到等距节点的情形,此时插值公式可简化,为此先介绍差分. 一、差分及其性质 差分的基本性质: * 可以证明 如 * 差分表 §3 Newton’s Interpolation 牛顿公式 Nn(x)= f(x0)+ f[x0,x1](x-x0)+...+ f[x0,...,xn](x-x0)...(x-xn-1) 牛顿前插公式 /* Newton’s forward-difference formula */ 注:一般当 x 靠近 x0时用前插,靠近 xn时用后插,故两种 公式亦称为表初公式和表末公式。 牛顿后插公式 /* Newton’s backward-difference formula */ 将节点顺序倒置: Nn(x)= f(xn)+ f[xn,xn-1](x-xn)+...+ f[xn,...,x0](x-xn)...(x-x1) §3 Newton’s Interpolation 使用Newton前插或后插公式,先构造差分表如下: x 0.4 0.5 0.6 0.7 sinx 0.38942 0.47943 0.56464 0.64422 例 已知数值表如下,分别用向前、向后Newton插值公式 求sin0.57891的近似值。 解:作差分表 §3 Newton’s Interpolation -0.00563 2 = 0.47943+ 0.7891×0.08521+ 0.7891×(0.7894-1) 误差 = 0.54714 若用向后插公式,则可取x0=0.6,x-1=0.5, x-2=0.4, x=x0+th, t=-0.2109,于是 §3 Newton’s Interpolation 使用向前插公式,取x0=0.5,x1=0.6, x2=0.7, x=x0+th, h=0.1, t=(x-x0)/h=0.7891,于是 N2(0.57891) = §3 Newton’s Interpolation * Newton插值法的优点是计算较简单,尤其是增加节点时, 计算只要增加一项,这是Lagrange插值无法比的. §2.5 埃尔米特(Hermite)插值 Lagrange插值公式所求得L(x)保证了节点处的函数值相等,也就是保证了函数的连续性,但 埃尔米特插值又称为带导数的插值,这类插值要求在节点处的导数值也要相等,导数的阶数越高则光滑度越高。 2.5.1 低次埃尔米特插值多项式 1.二点三次埃尔米特(Hermite)插值多项式 由于给出了四个条件,故可以唯一定出一个三次多项式,不妨设 由于给出了四个条件,故可以唯一定出一个三次多项式,不妨设 2.5.2 一般埃尔米特(Hermite)插值多项式 2.5.3 利用重节点的差商构造埃尔米特(Hermite)插值多项式 书上P36 例2.6 问:这种方法和前面方法得到的插值多项式是否相同? 2.三点三次带一个导数值的插值多项式 P34 §2.6 分段低次插值法 一、高次插值的病态性质 龙格(Runge)现象…. 二、分段线性插值 所谓分段线性插值就是用通过插值点的折线段逼近f(x). 二、分段三次埃尔米特插值 分段线性插值函数导数间断,若已知节点上函数值和导数,可构造一个导数连续的插值函数Ih(x),满足 * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
文档评论(0)