第二章：插值.pptVIP

Applied example 某地区某年夏季时节间隔 30 天的日出日落时间为 日照时间的变化设为 y(x)=a0+ a1x + a2x2, The most useful interplation function(最常用的插值函数是 …)? Uniqueness of Solution of problem？ How to construct interplation polynomial？ Error? §2.1 Uniqueness of Solution of problem代数插值问题解的存在惟一性 Given the different n+1 points ｛xj｝nj=0 belong to [a,b］, and f(xj)，j =0，…, n， Find nth-order pn(x)∈Pn，s.t. pn(xj)=f(xj)，j=0,1,…，n. …... (1) Let pn(x)=a0+a1x+…+anxn, …... (2) (只要证明pn(x)的系数a0 ,a1,…, an存在唯一即可) Remarks: 通过解上述方程组(3)求得插值多项式pn(x)的方法并不可取. 这是因为当n较大时解方程组的计算量较大，而且方程组系数矩阵的条件数一般较大（可能是病态方程组 ill-problem）,当阶数n越高时，病态越重. l0(x0)=1 , l0(x1)=0 , l0(x2)=0. l1(x0)=0 , l1(x1)=1 , l1(x2)=0. l2(x0)=0 , l2(x1)=0 , l2(x2)=1. It is easy to find the above basis function, for example: to seek l0(x), (因x1, x2 为其零点)so So L(x)=l0(x)y0+l1(x)y1+l2(x)y2 s.t. L2(xj)=yj (j=0,1,2) x0=4,y0=2，x1=9, y1=3 ，x2=16, y2=4. Basis Functions For each pn(x)∈Pn，we have pn(x)=c0 l0(x)+c1 l1(x)+…+cn ln(x) Where c0 ,c1 ,…,cn is to seek. 当 f(x) 是n次的多项式时, Ln(x)= f(x)。即n次多项式的n次插值函数即为该n次多项式本身。 Homework P45: 1, 3, 4 §2.3 Newton’s Interpolation 牛顿插值  Lagrange 插值虽然易算，但若要增加一个节点时，全部基函数 li(x) 都需要重新计算. Divided difference 差商(亦称均差) 差商表 xk f(xk) 一阶差商 二阶差商 三阶差商 …… n 阶差商 一般地,称k阶差分的差分为k+1阶差分,如二阶向前和向后差分分别为 差分表: 差分具有如下性质 性质4(差分与差商的关系): Ex 给定f(x)在等距节点上的函数值表如下: xi 0.4 0.6 0.8 1.0 f(xi) 1.5 1.8 2.2 2.8 分别用Newton向前和向后差分公式求f(0.5)及f(0.9)的近似值. Sol 先构造向前差分表如下: xi fi △fi △2fi △3fi 0.4 1.5 0.6 1.8 0.3 0.8 2.2 0.4 0.1 1.0 2.8 0.6 0.2 0.1 x0=0.4, h=0.2, x3=1.0. 分别用差分表中对角线上的值和最后一行的值,得Newton向前和向后插值公式如下: Lab 1 Lagrange Polynomial Given a set of sample points