插值方法.pptVIP

　　　　　第一章 插值方法 为了计算函数值或分析函数的性态，必 须首先由实验或观测数据找出函数关系 的一个近似表达式.插值与逼近就是用简 单函数为各种离散数据建立连续的数学 模型，使其既能达到精度要求，又使计 算量尽可能小.插值与逼近理论是数值计 算的最基本内容. 插值的概念 已知函数y=f(x)在n+1个互异点x0, x1, …, xn上的函数值分别为y0, y1, …, yn ，构造一个简单的函数P(x)，满足条件 　 　　P(xi) ＝yi 　(i=0,1,…n) 　　(1) 　称这类问题为插值问题，称P(x)为函数f(x)的插值函数， f(x)为被插值函数，点x0, x1, …, xn为插值节点，称(1)为插值条件. 多项式插值的概念 已知函数y=f(x)在n+1个互异点x0, x1, …, xn上的函数值分别为y0, y1, …, yn ，构造一个n次多项式?n(x)，满足条件 　 ?n(xi) ＝yi 　 (i=0,1,…n) 　　(2) 称这类问题为多项式插值问题，称?n(x)为函数f(x)的n次插值多项式，称点x0, x1, …, xn为插值节点，称(2)为插值条件. 插值多项式的唯一性 设满足条件?n(xi) ＝yi (i=0,1,…n)的插值多项式为a0+a1x+a2x2+…+anxn 由插值条件(２)知系数ai (i=0,1,…n)满足方程组 插值多项式的唯一性 ai的系数行列式是范德蒙行列式 xi当互异时方程组的解是存在且唯一的． 定理1 满足条件（2）的n次插值多项式?n(x)是存在且唯一的. Lagrange插值多项式 Lagrange插值问题 插值问题 　已知函数y=f(x)在n+1个互异点x0, x1, …, xn上的函数值分别为y0, y1, …, yn ，构造一个插值函数Ln(x)，满足条件 　1. Ln(x)为次数? n的多项式 2.Ln(xi) ＝yi　 (i=0,1,…n) 首先构造n次多项式li(x) (i=0,1,…n) ，满足 设li(x) =A(x-x0)(x-x1)…(x-xi-1)(x-xi+1)…(x-xn) 由插值条件li(xi) =1得 从而 称li(x) (i=1,2,…n)为Lagrange插值基函数． 易知 称此多项式为Lagrange插值多项式． 记?(x)=(x-x0)(x-x1) …(x-xn) 线性插值（一次插值） 　已知函数y=f(x)在两点x0, x1上的函数值分别为y0, y1，构造求一个一次式L1(x)，满足条件: L1(x0)= y0 ，L1(x1)= y1. 一次Lagrange插值多项式为 L1(x)= y0l0(x) + y1l1(x) 其中 抛物插值（二次插值） 已知函数y=f(x)在三个互异点x0, x1, x2上的函数值分别为y0, y1, y2，构造求一个二次式L2(x)，满足条件: L2(x0)= y0 ，L2(x1)= y1 ，L2(x2)= y2. 二次Lagrange插值多项式为 L２(x)= y0l0(x) + y1l1(x) + y2l2(x) 其中 Lagrange插值多项式的余项 定理2 设Ln(x) 是过点x0, x1, …, xn 的n次插值多项式，若f(x)?Cn[a,b] ， f(x)在(a, b)内存在n+1阶导数，其中 [a, b]是包含点x0, x1, …, xn的任一区间，则对任意给定的x?[a,b] ，总存在一点??(a, b) （依赖于x）使 Lagrange插值多项式的余项定理的证明 证 当x为插值节点x0, x1, …, xn中任一点时，结 论显然成立.下面设x异于x0, x1, …, xn ，由于 R(x)=f(x)-Ln(x) 满足R(xi)=0 , 故可设R(x)=K(x- x0)(x-x1)…(x-xn) ，其中K为依赖于x的待定函数. 固定x，作辅助函数 G(t)= f(t)-Ln(t) - K(t-x0)(t-x1)…(t-xn) 显然G(t)在[a, b]上有n+2个零点x ,x0, x1, …, xn ； 利用Rolly定理,知 G‘(t)在(a, b)内至少有n+1零点; 反复利用Rolly定理:G‘’(t)在(a, b)内至少有n零点; Lagrange插值多项式的余项定理的证明 …… G(n+1)(t)在(a, b)内至少有1零点;即存在一点 ?，使G(n+1)(?) =0. 由于G(n+1)(t)= f(n+1)(t) –(n+1)! K ，从而 所以 线性插值多项式的余项 抛物插值多项式的余项