第一章插值法(下)探析.ppt

课堂小测验 给出 及 的推导方法 多点Hermite多项式插值误差估计 Thank you! * * * * * * * 在数学分析中我们用y=f(x)描述一条平面曲线，但在实际问题中，函数y=f(x)往往是通过实验观测得到的一组数据来给出的，即在某个区间[a,b]上给出一系列点的函数值yk=f(xk)，或者给出一个函数表，如何通过这些对应关系去找函数f(x)的近似表达式呢？就可以利用插值． 简单的说，插值的目的就是根据给定的数据表，寻找一个解析形式的函数p(x)，近似代替f(x) 例： 已知x=1,4,9的平方根为1,2,3，利用牛顿基本差商 公式求 的近似值。 解： 从而得二阶牛顿基本差商公式为 因此计算得 的近似值为 Newton插值公式余项 1 2 … … … … n+1 Nn(x) Rn(x) ai = f [ x0, …, xi ] Newton插值公式余项 反复的把后一项带入前一项： Newton插值公式及其余项 差商具有如下性质: Warning: my head is exploding… What is the point of this formula? 差商的值与 xi 的顺序无关！ 定理4 P26 练习 牛顿插值 牛顿插值方法的引出 牛顿插值方法的推导 等距牛顿插值 上面我们讨论了节点任意分布的插值公式，但实际应用时经常会遇到等距节点的情形，这时插值公式可以进一步简化，计算也简单多了，为了给出等距节点的插值公式，我们先来看一个新概念； 利用一阶差分可以定义二阶差分 差分表 差分与函数值之间的关系 归纳可知，k阶差商可表示为 在等距节点的前提下,差商与差分有如下关系 依此类推 由差商与差分的关系 Newton插值基本公式为 如果假设 Newton差分插值公式 则插值公式 化为 其余项 化为 称 为Newton差分插值公式 插值余项为 Newton插值法的优点是计算较简单,尤其是增加节点时, 计算只要增加一项,这点是Lagrange插值无法比的. 但是Newton插值仍然没有改变Lagrange插值的插值曲线 在节点处有尖点，不光滑，插值多项式在节点处不可导 等缺点. 例子 xi yi Δyi Δ2yi Δ3yi Δ4yi 1 2.71828 1.5 4.48169 1.76341 2 7.38906 2.90737 1.14396 2.5 12.18249 4.79343 1.88606 0.74210 3 20.08554 7.90305 3.10962 1.22356 0.48146 例子（解） 例子（解） 给出函数表 x x0 x1 x2 …… xn y y0 y1 y2 …… yn 其中 求 满足 插值方法小结 由插值条件: ,可得 待 定 系 数 法 拉格朗日（Lagrange）插值法 设 用已知的函数值作为组合系数去待定相应的基函数．即 其中 基函数方法 Newton插值法 其基函数为 其中待定系数为 余 项 校 正 法 提纲 埃特金插值 牛顿插值 埃尔米特插值 Hermite插值 Hermite插值 Hermite插值 第一章 插值方法 李书杰 合肥工业大学 计算机学院 Lagrange插值多项式的缺点 我们知道,Lagrange插值多项式的插值基函数为 理论分析中很方便，但是当插值节点增减时全部插值基函数就要随之变化，整个公式也将发生变化，这在实际计算中是很不方便的； Lagrange 插值虽然易算，但若要增加一个节点时，全部基函数 li(x) 都需重新算过。 提纲 埃特金插值 牛顿插值 埃尔米特插值 埃特金算法(Aitken) 算法的承袭性在计算方法中有很重要的作用，这样可以节省计算量。这里介绍一种具有这种性质的算法：埃特金算法 . 一、符号说明 对于给定的插值点 ，如果除顺序排列的 个节点 外，再增加一个节点 进行 次插值，则插值结果依赖于所给定的次数 与节点 ，记为 。 例如， 表示取 进行线性插值，即 特别地 二、例子分析 又例如， 表示取 进行抛物插值的结果，特别有 或（改变记号） 线性插值 一般地有 约定： 例5 利用下表左部所给数据求正弦积分 在x=0.462的值。 已知 线性