(西南交大戴克俭版)计算方法多项式插值概述.ppt

第五章 插值法 第5章 插值法 5.1 插值法多项式存在性与唯一性 5.1 插值法多项式存在性与唯一性 5.1 插值法多项式存在性与唯一性 5.1 插值法多项式存在性与唯一性 5.2 拉格朗日插值法 5.2 拉格朗日插值法 5.2 拉格朗日插值法 5.2 拉格朗日插值法 拉格朗日插值公式(5.6)形式对称，结构紧凑，容易编写计算机程序。下图为其流程图。 5.2 拉格朗日插值法 5.2 拉格朗日插值法 5.2 拉格朗日插值法 5.2 拉格朗日插值法 5.2 拉格朗日插值法 5.3 牛顿插值公式 一、差商及其基本性质 定义1 称 为函数f(x)在x0、x1 点的一阶差商。一阶差商的差商： 称为函数f(x)在x0、x1、x2点的二阶差商。一般地， n-1阶差商的差商： 称为f(x)在x0、x1、…、xn点的n阶差商。 5.3 牛顿插值公式 性质1 差商可以表示为函数值的线性组合，即 (5.13) 其中 (5.14) 差商的这一性质称之为差商的对称性。 性质2 性质3 若f(x)在[a,b]上存在n阶导数，且节点x0、x1、 …、xn?[a,b]，则至少存在一点??[a,b]，使 5.3 牛顿插值公式 二、牛顿插值多项式 设x是[a,b]上的一点，则有一阶差商定义有 得 f(x)=f(x0)+f[x, x0](x-x0) 同理，由二阶差商定义 得 f[x, x0]=f[x0, x1]+f[x, x1, x0](x-x1) 如此继续下去，可得一系列等式： 5.3 牛顿插值公式 依次将后式代入前式，最后得： 5.3 牛顿插值公式 可见，式(5.15)所示的Pn(x)为次数不超过n的多项式，且由公式(5.17)易知Pn(x)满足插值条件(5.1)，故其为多项式插值问题的解，称之为牛顿插值多项式。其系数就是差商表5.1中横线标出的各阶差商。式(5.16)为牛顿插值多项式的余项。在实际应用中，公式(5.15)常被写成如下递推形式： 5.3 牛顿插值公式 三、等距节点的牛顿插值公式 设函数y=f(x)在等距节点xi=x0+ih(i=0,1,…,n)上的函数值为fi=f(xi)，这里步长h为常数。 定义2 分别称为f(x)在点xi处的一阶向前差分和一阶向后差分。一阶差分的差分 分别称为f(x)在点xi处的二阶向前差分和二阶向后差分。 一般地，m-1阶差分的差分为： 分别称为f(x)在点xi处的m阶向前差分和m阶向后差分 5.3 牛顿插值公式 差分有如下基本性质： 性质1 性质2 ⒈ 牛顿向前插值公式 设小x?[x0,xn]且小x在x0附近，则可令x=x0+th(0t1)，将此式及(5.18)代入牛顿插值公式(5.15)可得 5.3 牛顿插值公式 其余项公式为： ⒉ 牛顿向后插值公式及余项 令x=xn+th(-1t0)，可得 5.4 埃尔米特插值 一、问题的提出 为了使插值函数能更好地切合原来的函数，许多问题不但要求两者在节点上函数值相等，而且还要求导数值相同，甚至要求高阶导数也相等。这类插值问题称为埃尔米特(Hermite)插值。本节仅讨论两个节点且已知这两个节点上的函数值和一阶导数值的情形。 定义 设y=f(x)是区间[a,b]上的实函数，x0、x1是[a,b]上相异两点，不妨设x0x1，已知y=f(x)在xi上的函数值yi=f(xi) (i=0,1)和一阶导数值mi=f ?(xi) (i=0,1)，求三 5.4 埃尔米特插值 次多项式H3(x)，使其满