第二章插值法12节.ppt

* * * * * * * * 问题的提出 拉格朗日插值 牛顿插值 分段低次插值 埃尔米特插值 第二章 插值法 在生产实际及科学试验中，经常要研究变量之间的函数关系，但很多情况下很难找到具体的函数表达式，往往只能通过观测或测量得到一张数据表： x x0 x1 x2 …… xn y=f(x) y0 y1 y2 …… yn 表中给出某个区间[a,b]上给出一系列点的函数值yi=f(xi)。 §1问题的提出 问题：无法求出不在表中的点的函数值，也不能进一步研究函数的其他性质，如函数的积分和导数等。为了解决这些问题，需要设法通过这张表格求出一个简单函数p(x)来近似f(x) ，使得p(xi) = f(xi) (i = 0, … n) 。 y=f(x) y=p(x) ——插值问题 已知精确函数 y = f(x) 在一系列节点 x0 … xn 处测得函数值 y0 = f(x0), … yn = f(xn)，由此构造一个简单易算的近似函数 p(x) ? f(x)，满足条件p(xi) = f(xi) (i = 0, … n)。这里的 p(x) 称为f(x) 的插值函数。最常用的插值函数是 …? 多项式 x0 x1 x2 x3 x4 x p(x) ? f(x) 定义：设函数y = f(x)在区间[a,b]上有定义，且给出一系列点上的函数值yi=f(xi) (i=0,1,2,…,n)，若用函数p(x)去近似代替它，使得 p(xi)= yi (i=0,1,2,…,n) 则称函数p(x)为f(x)的插值函数；称x0,x1,… xn称为插值节点或简称节点。插值节点所界的区间[a,b]称为插值区间。p(xi)= yi 称为插值条件。 本章只讨论多项式的插值问题，即构造n次多项式 Pn(x)= a0 + a1x + a2x2+…+ anxn 使满足Pn(xi)= yi (i=0,1,2,…,n)，及利用Pn(x)进行插值计算的问题。 Pn(x)称为插值多项式。 定理 (插值多项式的存在唯一性) 满足 的 n 次插值多项式是存在且唯一的。 证明： 反证：若不唯一，则除了Pn(x) 外还有另一 n 阶多项式 Ln(x) 满足 Ln(xi) = yi 。 考察 则 Qn 的阶数 ? n 而 Qn 有 个不同的根 n + 1 x0 … xn 这与代数基本定理n次多项式只有n个根矛盾，故 ???????????? ? . 证毕. 定理 (插值多项式的存在唯一性) 满足 的 n 阶插值多项式是存在且唯一的。 插值多项式就是根据给定n+1个点 ??????????????????? ，求一个n次多项式： ????? n i y x P i i , ... , 0 , ) ( = = 使 ????这里 ?????????? 是n+1个待定系数，根据n+1个条件得到的方程组是关于参数 ?????????? 的线性方程组。 所以解是存在唯一的。但直接求解较复杂，也得不到统一的表达式。所以通常求插值多项式不用这种方法，而使用下节给出的构造基函数方法。 系数矩阵为范德蒙行列式，当节点互异时 §2 拉格朗日插值 /* Lagrange Interpolation */ n i y x P i i n , ... , 0 , ) ( = = 求 n 次多项式 使得 条件：无重合节点，即 已知 x0 , x1 ; y0 , y1 ，求 使得 1 1 1 0 0 1 ) ( , ) ( y x P y x P = = 可见 P1(x) 是过 ( x0 , y0 ) 和 ( x1, y1 ) 两点的直线。 ) ( ) ( 0 0 1 0 1 0 1 x x x x y y y x P - - - + = 1 0 1 x x x x - - 0 1 0 x x x x - - = y0 + y1 l0(x) l1(x) ? = = 1 0 ) ( i i i y x l 称为线性插值基函数 满足条件 li(xj)=1(i=j)或0( i?j ） 线性插值 n=1 直线方程的两点式： l0(x) l1(x) ? = = 1 0 ) ( i i i y x l L1