实验5插值.pptVIP

插值法(Interpolation Method) 邓国强 数学与计算科学学院 实验内容 Newton’s Interpolation Lagrange Interpolation Spline Interpolation 构造基函数 差商表 . xk f(xk) 一阶差商 二阶差商 三阶差商 …… n 阶差商 (1) 输入: xi, fi; di=fi (i=0,1,…,n); 计算差商 对i=1,2,…,n 做 (3.1) 对j=i,i+1,…,n 做 fj=(dj-dj-1)/(xj-xj-i); (3.2) 对j=i,i+1,…,n 做 dj=fj; (4) 计算插值N(u) (4.1)输入插值点u; (4.2) v=0; (4.3) 对i=n,n-1,…,1,0 做 v=v(u-xi)+fi; (5) 输出u,v. Cubic Spline Interpolation Lagrangr Interpolation 为了利用多项式插值而克服高次插值的缺陷，便引入分段插值，Cubic Spline Interpolation 也是一种分段插值的方法 因 是分段3次多项式 ,故在每个区间 上 都是3次多项式 ,从而 共须 个独立条件确定 . ① 和 在 个内结点连续,即满足条件(4.4),因而 (4.4)给出了 个条件； ②(4.2)提供了 个独立条件; ③还差2个条件,有多种给法.最常见的给法是: (i) （简支边界，导致三弯矩关系式, 关系式）, 特别地, (自然边界,三次自然样条); (ii) (固支边界,导致三转角关系式, 关系式). * 当精确函数 y = f(x) 非常复杂或未知时，在区间[a,b]上一系列节点 x0 … xm 处测得函数值 y0 = f(x0), …, ym = f(xm)，由此构造一个简单易算的 近似函数 g(x) ? f(x)，满足条件 g(xj) = f(xj) (j = 0, … m) (*) 这个问题称为“插值问题” 插值问题的定义 这里的 g(x) 称为f(x) 的插值函数。 节点 x0 … xm称为插值节点, 条件（*)称为插值条件，区间[a,b]称为插值区间 x0 x1 x2 x3 x4 x f(x) g(x) 精确函数 y = f(x) 非常复杂或未知 g(x) 称为f(x) 的插值函数。 节点 x0 … x4称为插值节点 基本思想：在n次多项式空间Pn中找一组合适的基函数 ?0(x),?1(x),…, ?3（x),使 pn(x)=a0 ?0(x) +a1 ?1(x) +…+an ?3（x) 不同的基函数的选取导致不同的插值方法 Lagrange插值 Newton插值 （2） 与 节点有关，而与f 无关 这里每个lj(x)都是n次多项式，且由(1)式容易验证 lj(x）满足 j=0,1,…，n （1） Ln(x)是插值基函数 ? Lagrange Interpolation基本原理 因此得到插值多项式 Ln(x) = f(x0)l0(x)+f(x1) l1(x)+…+ f(xn) ln(x) 记为Ln(x)＝? f(xj)lj(x) 称Ln(x)为n次Lagrange插值多项式 ? 计算实习： Lagrange Polynomial 1 输入n,x,y,x* 2 赋初始值Pai=1.0，Slag=0.0 3for i=0，1，…，n 3.1Pai*(x*-xi)→Pai end for(i)  4for j=0,1,…，n 4.1Tai=1.0 4.2for i=0,1，…，n 4.2.1if　i≠j then Tai*(xj-xi)→Tai end if end for(i)  4.3Slag+((yj