数值第六章_数值插值方法教材.pptVIP

第六章 插值方法;6.1 代数插值问题;6.1 代数插值问题;日照时间的变化设为 y(x)=a0+ a1x + a2x2,;定义;设 Pn (x)=a0 + a1x + a2x2 + ……+ anxn 是y=f(x)在[a,b]上的n+1个互异节点x0,x1,…,xn的插值多项式，则求Pn (x)问题归结为求系数a0,a1,…,an。;故Pn (x)存在且唯一。虽然直接求解上述方程组可求得插值多项式，但繁琐复杂，一般不用。;给定插值节点 x0,x1， y0=f(x0)，y1=f(x1). 求线性插值多项式L1 (x)=a0+ a1x，使满足: L1(x0)=y0 , L1(x1)=y1.;由两点式可以看出， L1 (x)是由两个线性函数;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.;满足上式的插值基函数很容易求出。如求l0(x), 因x1, x2 为其零点，故可表为;显然 L(x)=l0(x)y0+l1(x)y1+l2(x)y2 满足条件 L2(xj)=yj (j=0,1,2); 取x0=4,y0=2，x1=9, y1=3 ，x2=16, y2=4.;设有n+1个互异节点x0 x1…xn，且 yi=f(xi) (i=0,1,2…,n) 构造Ln (x)，使 Ln (xj)= yj (j = 0,1,2,…,n);由n=1,2时的讨论可得 ;定理 设 f(x)在[a,b]上具有n阶连续导数, 且 f (n+1)(x) 存在,在节点a ≤ x0 x1…xn≤b， Ln (x)是满足条件Ln (xj)= yj (j = 0,1,2,…,n)的插值多项式，则对任何x?[a,b]，插值余项;证明：;;根据Rolle定理,;所以;当 f(x) 是n次的多项式时, Ln(x)= f(x)。即n次多项式的n次插值函数即为该n次多项式本身。;例:;6.3 Newton插值;例：;xk;这说明截断误差很小。;此例中，五阶均差f[x,x0,x1,……,x4]是用f[x0,x1,……,x5]来近似的。;6.4 分段低次插值;Runge反例:;n;从表中可以看出，随着n的增加，R(xn-1/2)的绝对值几乎成倍地增加，这说明当n－∞时Ln在[-5,5]上不收敛。;一、分段线性Lagrange插值;;内插;故也称折线插值，如右图：;由前述余项定理可知，n次Lagrange插值多项式的余项为：;二、分段二次Lagrange插值;2. 分段二次插值的误差估计;例:;同理;(2) 分段二次Lagrange插值的公式为;样条：是 指飞机或轮船等的制造过程中为描绘出光滑的外形曲线(放样)所用的工具。;一、三次样条插值函数的定义;(1)S(x)在每一个子区间[xj-1 , xj ] ( j= 0,1,2,···,n)上是一个三次多项式； (2) S(x)在每一个内接点xj ( j= 1,2,···,n-1)上具有直到二阶的连续导数；;三次样条插值多项式的确定：;可得3n-3个方程，又由条件(3);例 已知f(x)：f(-1)=1,f(0)=0,f(1)=1，求f(x)在[-1,1]上的三次自然样条插值函数， f(x) 在-1,1处的二阶导数为0。;由一阶及二阶导数连续得：; 用二阶导数值构造三次样条插值函数 （又称三弯矩算法）;其中c1，c2为积分常数。;上式称M表达式，只需确定Mi，即可确定三次样条插值函数。;故;记;得;则求Mi的方程的矩阵形式为;上述方程组是关于Mj(j=1,…,n-1)的三对角方程组。;6.5 埃尔米特插值（Hermite）;例1.;作为多项式插值，三次已是较高的次数，次数再高就有可能发生Runge现象，因此，对有n+1个节点的插值问题，我们需要使用分段插值等其他方法。;二. 添加“修正项”的算法;P137 习题五：5 ，11，13，16