北京科技大学计算方法7第七章函数的插值法5v讲义.ppt

北京科技大学数理学院 卫宏儒 Weihr168@ 第7章 插值法 估计f(x)在区间[a,b]中某点 的值时，当 属于包含结点 的最小闭区间时，相应的插值称为内插，否则称为外插。 在某一逼近函数类中选取的一组线性无关的函数 ，此时对应的插值函数 为： 由插值条件确定 函数组 称为插值基函数。 Lagrange插值 例2: 已测得某地大气压强随高度变化的一组数据 Newton插值 插值型求导公式 一.问题描述 二.定义 三.定理 四.构造函数 五.例题 六.一般插值 假设函数y=f(x)是在[a,b]上有一定光滑性的函数,在xo…xn处是n+1个异点,f(x)在这些点上取值yo…...yn.求一个确定的函数p(x)在上面n+1个点上满足p(xi)=yi i=0,1,…,n.这是最简单的插值问题,如果除了知道f(x)在插值基点上的取值外,还知道f(x)在插值基点上的其他描述(如知道f(x)在插值基点上的导数值)。如何来构造插值函数呢? f(x) 是区间[ a, b] 上 n+1个互异节点a=x0x1x2……xn=b ， 定义在[a,b]上的函数f(x) 在节点上满足 f(xi) = yi f’(xi)=y i i=0,1,2……n 求一个次数不高于2n+1次的插值多项式H(x)满足2n+2个条件 H(xi) = yi H (xi)= y i i=0,1,2……n 若H(x)存在，则称为函数f(x) 的Hermite插值多项式。因为 H(x)是一个次数不高于2n+1次的多项式，常记为H2n+1(x). 定理一:满足插值条件 H(xi)= yi H(xi)= yi i=0,1,2……n 且次数不大于2n+1的多项式是唯一的。 定理二 :f(x)在区间[a,b]存在2n+2阶导数,则其Hermite插值余项为: 设Hermite插值函数 n n H2n+1(x) = ? Li(x) yi + ? hi(x) yi i=0 i=0 Li(x)，hi(x)都是不高于2n+1次的多项式，类似Lagrange插值，利用Hermite插值条件可得： Li(xj)=?ij hi(xj) = 0 Li(xj)=0 hi(xj)= ?ij i,j=0,1,2……n 从而可设 Li(x)= (aix+bi)[li(x)]2 hi(x)= (cix+di)[li(x)]2 求过0,1两点构造一个三次插值多项式,满足条件 f(0)=1, f (0)=1/2 , f(1)=2, f (1) =1/2 解: 设 H3(x)=Y0l0(x)+y1l 1(x) +y 0 h0(x) +y1h 1(x) 因为 l0(x)=(ax+b)(x-1)2 利用 l0(0)=1和l0 (0)=0, 得 b=1,a=2. 所以有l0(x)=(2x+1)(x-1)2 同理可得 l1(x)=(3-2x)x2 h0(x)=x(x-1)2 h1(x)=x2(