第二章__插值法.pptVIP

第二章 插值法 §0 插值问题 §1 线性插值 §2 二次插值 §3 n次插值 §4 分段线性插值 §5 Hermite插值 §0 插值问题 实际问题中，函数 y = f (x) 往往是通过实验观测数据给出的，即在某个区间［a, b］上给出一系列点的函数值 yi = f (xi), i=0, 1, 2, …, n (2.1) 或者给出一张函数表： 目标：根据给定的函数表，构造出一个函数? (x)，使得 ? (xi)=yi , i=0,1,2,…,n 作为函数 y =f(x)的近似表达式。 用途： 1）利用 ? (x) 近似计算出函数表中没有列出的点上的函数值； 2）利用 ? (x) 大致分析原函数 f(x) 的性质。 构造方法：1）插值法；2）数据拟合法 插值法定义： 给定函数表 ， 要构造一个多项式?n (x)，使其满足条件： （1） ?n(x)是一个不超过n次的多项式； （2） ?n(xi) = yi = f(xi)，i=0,1,2,…,n (2.2) 称 ?n(x) 为 f(x) 的 n 次插值多项式， f(x)为被插函数，x0，x1，x2，?，xn为插值节点。 §1 线性插值 已知： (x0, x1互异)，求?1(x)。 由定义，?1(x)不超过1次，且与f(x)取相同的两点。 最直接的方法： 将 (x0, y0)，(x1, y1)代入 ?1(x)＝a1x+a0 一、?1(x)的构造 Lagrange（拉格朗日）线性插值 用两点式表示直线?1(x) 插值基函数性质 插值基函数个数与节点个数相同 例：计算 的值 Newton（牛顿）线性插值 用点斜式表示直线?1(x) 用行列式形式表示直线?1(x) 二、线性插值的误差 定理2.1 设给定 (x0, x1互异)， ?1(x) 是过 x0, x1 的线性插值，[a, b] 是包含 (x0, x1) 的任一区间，并设 f(x)?C1[a,b]，f”(x) 在 [a, b]上存在，则对任意给定的 x?[a, b] 总存在一点 ??(a, b)（依赖于x）使 且，若 x 在 x0 和 x1 之间，则 证明： 说明： 1）[a,b]区间的取定：一般书上采用的方法是 a=min(x0, x1），b=max(x0, x1） 2）对不同的x，?不同，即? 的选取依赖于x §2 二次插值 已知： (x0, x1, x2互异)，求?2(x)。 一、?2(x)的构造 Lagrange二次插值 构造二次插值基函数l0(x)，l1(x)，l2(x) 推出 Newton二次插值 定义 f(x) 过点( xi, yi )，( xj, yj )，( xk, yk )的二阶均差为： 设 ?2(x)＝A + B (x-x0) + C (x-x0) (x-x1) 由插值条件 逐次线性插值 例：已知函数 y=f(x) 的观测数据如下，试求拉格朗日插值多项式。 二、?2(x)的误差 定理2.2 设?2(x)是过x0, x1 , x2的二次插值多项式，f(x)?C2[a,b]，f”’(x)在[a,b]上存在，其中[a,b]是包含点x0, x1, x2的任一区间，则对任意给定的x?[a,b]，在(a,b)上总存在一点?（依赖于x）使 例：取节点x0＝0, x1＝1, x2＝1/2，对函数 y=e-x 做二次插值，并分析误差。 Lagrange插值 3.行列式插值 误差: 例：已知函数 y=f(x) 的观测数据如下，试求拉格朗日插值多项式。 §3 n次插值 已知： (x0, x1,…, xn互异)， 构造?n(x)，满足条件： 1）?n(x)是次数不高于n的代数多项式 2）?n(xi)= f(xi)＝ yi，i=0,1,…,n 一、?n(x)的存在性和唯一性 由于?n(x)满足插值条件，故可设： ?n(x)＝a0+a1x+a2x2+…+anxn 又∵?n(xi)= f(xi)＝ yi，i=0,1,…,n ∴