插值与拟合教学幻灯片讲义.ppt

第 五 章 插值与拟合;5.1 插值的基本概念;5.1 插值的基本概念;插值法;　　设 f(x)为[a，b]上的函数，在互异点x0 , x1, ... , xn 处的函数值分别为 f(x0) , f (x1) , … , f (xn) ,构造一个简单函数 ?(x) 作为函数 f(x) 的近似表达式y= f(x) ? ?(x),使 ?(xi)＝f(xi) , i＝0, 1, 2, …,n 　　 (5.0) 则称?(x) 为关于节点x0 , x1, ... , xn的插值函数；称 x0 , x1, ... , xn 为插值节点；称(xi, f (xi)), i=1,2,… , n 为插值点；f(x) 称为被插值函数。 (5.0)式称为插值条件。这类问题称为插值问题。 构造出?(x)，对 f(x)在[a,b]上函数值的计算，就转化为?(x)在对应点上的计算。;5.2 Lagrange插值 选用代数多项式作为插值函数。Lagrange插值就是选用节点上的函数值作为插值条件。;解之得，;插值基函数法;显然， l0(x), l1(x) 满足：　线性插值函数可以写成节点上函数值的线性组合，即 L1(x) = l0(x) y0 + l1(x) y1 ．称l0(x), l1(x) 分别为x0, x1的插值基函数。;可设 R1(x) = k(x)(x-x0)(x-x1).;例5.1 给定sin11°=0.190809,sin12°=0.207912,求线性插 值，并计算sin11°30和sin10°30 。;例5.2 给定sin11°=0.190809,sin12°=0.207912,  sin13°=0.224951,构造二次插值，并计算  sin11°30′。;例5.3 要制作三角函数sin x的值表，已知表值有四位小数，要求用线性插值引起的截断误差不超过表值的舍入误差，试确定其最大允许的步长。;5.2.3 n次Lagrange插值多项式; 所以，(5.6)的解存在唯一。 解出(5.6)的解a0 , a1, …, an ，代入(5.6)1即得 n 次插值多项式Pn(x)。 ★ n 次插值多项式的构造 有上面的讨论知，用待定系数法要求解一个线性方程组，计算量大，实际中不可取。;插值基函数法;先构造 l0(x)。有上表知， x1 , x2, …, xn 为 l0(x) 的零点，设;于是，;例5.4 已知插值点 (-2.00,17.00), (0.00,1.00), (1.00,2.00), (2.00,17.00), 求三次插值，并计算 f (0.6)。; 三次Lagrange插值多项式为： f (0.6) ≈L3(0.6) = -0.472. ;★ n 次插值多项式的误差 定理 2 设Ln(x)是[a,b]上过插值点(xi, f(xi))的 n 次 Lagrange插值???项式，节点xi两两互异，若 f(x)∈Cn+1[a,b], 则? x ∈[a,b],存在ζ=ζ(x) ∈[a,b],使得 证明提示：因为Rn(xi)=0, i =0, 1, 2, …, n, 令 作辅助函数 显然， Ψ(t)有 n+2 个零点x0，x1，…，xn，x，反复应用Rolle 定理，由Ψ(n+1)(ζ)=0,定理得证。;★ n 次插值多项式的几点说明;实际计算中，如果 f (x)的形式未知，也就无法求出 f (n+1)(x)或估计出其较精确的上界。所以也就无法估计|Rn(x)|。 我们采用事后估计，即给定 n+2个节点，x0, x1,…, xn, xn+1,构造两个n次插值多项式Ln(x)和 ，用它们来估计|Rn(x)|。 由定理 2，得;假定 f (n+1)(x) 在 [a, b]内连续, 且变化不大, 即 f (n+1)(ζ1) ≈ f (n+1)(ζ2).; 所以我们得到 n 次Lagrange插值多项式： ;;5.3 Newton插值; 一般地，k 阶差商 f [x0,