计算方法电子教案第二章插值法与数值微分课件.pptVIP

第二章 插值法与数值微分 Newton插值 Newton插值 把直线方程利用点斜式表示： 因此, 是f(x) 在 x1 , x0处的一阶 均差， 利用均差的对称性， （＊） 式可以表示为： 这种形式的插值叫做牛顿（Newton）插值． 定理2.1：设给定 x x0 x1 y y0 y1 是过x0,x1的线性插值函数，[a,b]是包含(x0,x1)的任一区间，并设 在[a,b]上存在，则对任意给定 , 总存在一点 使 且可以证明： roll定理 如果函数f(x)满足： 在闭区间[a,b]上连续； 在开区间(a,b)内可导； 在区间端点处的函数值相等，即f(a)=f(b)， 那么在(a,b)内至少有一点ξ(aξb)，使得 f(ξ)=0 2.2 二次插值 线性插值是用两点 和 来构造y=f(x)的插值函数，下面用三个点 来构造y=f(x)过三点的插值函数。（过三点可以作一条抛物线）。 为了确定系数将三点代入方程得： 构造l0(x):由于l0(x)有x1和x2二个零点，因此有因子(x-x1)(x-x2)，又因有l0(x）是一个次数不高于二次的多项式，所以，还可能相差一个常数因子，于是把l0(x）写成： 利用条件 ，可求得： 那么x0点的二次插值基函数为： 同理构造 为： 插值函数： 牛顿二次插值多项式： 假设过 的二次插值多项式具有下面的形式： 确定系数A,B,C： 利用 我们有： 再利用 有： B = 最后 确定C： （为一阶均差的均差） 其中 f [xi, xj , xk ]为f (x)在点xi, xj, xk处的二阶差商 一般的称 为f (x)在点x0, x1,…, xn处的n阶差商。 于是得到二次牛顿插值多项式： 例2: 已测得某地大气压强随高度变化的一组数据 例3.取节点x0=0,x1=1和　　　　　　　对　　　　　　　　　　　　建立线性插值多项式和二次插值多项式。 解：先构造x0=0,x1=1两点的线性插值多项式。 x 0 1 y 1 e-1 （1）拉格朗日插值多项式 先选过(0,1)和(1, e-1)的一次插值函数 这样： = (2)牛顿型插值多项式： 因为 ，所以 因此 = 因此 2.3 n次插值 设给定函数表 x x0 x1 ……xn y y0 y1 ……yn 要求构造一个插值多项式 满足条件： （1） 是次数不高于n的多项式 （2） 写成方程组形式： 拉格朗日型n次插值多项式： （1）先构造n+1个插值节点x0,x1,……xn上的n次插值基函数li(x) （2） li(x)的数值表： x0 x1 ,…… xn l0(x) 1 0 ,…… 0 l1(x) 0 1 ,…… 0 ………… ln(x) 0 0 ,…… 1 （3）确定li(x)的零点，构造li(x): （4）利用 得： 为了得到n次牛顿型插值多项式：（1）构造均差表： 这里 N次插值表见教