数值分析课件第六章_数值插值方法选编.pptVIP

第六章 插值方法 插值的基本概念 Lagrange插值 Newton插值 分段插值 埃尔米特（Hermite）插值 6.1 代数插值问题 6.1 代数插值问题 例. 某地区某年夏季时节间隔 30 天的日出日落时间为 插值：研究用简单函数为各种离散数据建立连续数学模型的方法。 日照时间的变化设为 y(x)=a0+ a1x + a2x2, 求出a0，a1，a2，即可得到5、6月份的日照时 间的变化规律。 根据三组数据: (1, 13.2167)， (31, 14.35)， (61, 14.6667) 导出关于a0，a1，a2的线性方程组 定义 已知函数y=f(x)在[a,b]有定义，且已知它在n+1个互异节点 a ≤ x0 x1…xn≤b 上的函数值 y0=f(x0)，y1=f(x1) ,…,yn=f(xn)， 若存在一个次数不超过n次的多项式 Pn (x)=a0 + a1x + a2x2 + ……+ anxn 满足条件 Pn (xk)= yk (k = 0,1,…,n) 则称Pn (x)为f(x)的n次插值多项式。（分段、三角插值） 点x0,x1,…,xn称插值节点， f(x)为被插值函数。[a,b]称插值区间，点x称插值点。 插值点在插值区间内的叫内插，否则叫外插。 设 Pn (x)=a0 + a1x + a2x2 + ……+ anxn 是y=f(x)在[a,b]上的n+1个互异节点x0,x1,…,xn的插值多项式，则求Pn (x)问题归结为求系数a0,a1,…,an。 定理 n次多项式插值问题的解是存在且唯一的。 证明： 由插值条件： Pn (xk)= yk (k = 0,1,…,n) 得关于a0,a1,…,an的n+1阶线性方程组 故Pn (x)存在且唯一。虽然直接求解上述方程组可求得插值多项式，但繁琐复杂，一般不用。 因 故上式不为0。 据Cramer法则，方程组解存在且唯一。 其系数行列式是Vandermonde行列式 给定插值节点 x0,x1， y0=f(x0)，y1=f(x1). 求线性插值多项式L1 (x)=a0+ a1x，使满足: L1(x0)=y0 , L1(x1)=y1. 6.2 Lagrange插值 一、线性插值与抛物插值 1. 线性插值：n=1情形 y= L1 (x)的几何意义就是过点(x0, y0)，(x1, y1)的直线。 L1 (x)的表达式： 点斜式: 两点式: 由两点式可以看出， L1 (x)是由两个线性函数 的线性组合得到，其系数分别为y0， y1。即 显然， l0 (x)及l1 (x)也是线性插值多项式，在节点x0,x1上满足条件： l0(x0)=1 , l0(x1)=0. l1(x0)=0 , l1(x1)=1. 称l0 (x)及l1 (x)为线性插值基函数。 (j,k=0,1) 即 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. 2. 抛物插值：n=2情形 假定插值节点为x0, x1, x2 ，求二次插值多项式 L2 (x)，使 L2(xj)=yj (j=0,1,2) y= L2 (x)的几何意义就是过 (x0, y0)，(x1, y1) ，(x2, y2)三点的抛物线。 采用基函数方法，设 L2 (x)=l0(x)y0+l1(x)y1+l2(x)y2 此时基函数l0(x), l1(x), l2(x)是二次函数，且在节点上满足： 满足上式的插值基函数很容易求出。如求l0(x), 因x1, x2 为其零点，故可表为 故 即 (j,k=0,1,2) 其中A为待定系数，由l0(x0)=1 , 得 显然 L(x)=l0(x)y0+l1(x)y1+l2(x)y2 满足条件 L2(xj)=yj (j=0,1,2) 同理 将l0(x), l1(x), l2(x)代入得 取x0=4,y0=2，x1=9, y1=3 ，x2=16, y2=4. 取x0=4, x1=9, x2=16 例 已知 求 解 （1）线性插值： 取x0=4, x1=9 （2）抛物插值： 设有n+1个互异节点x0 x1…xn，且 yi=f(xi) (i=0,1,2…,n) 构造Ln (x)，使 Ln (xj)= yj (j = 0,1,2,…,n) 二、Lagrange插值多项式 定义 若n次多项式lj(x) (j = 0,1,…,n)在n+1个节点x0 x1…xn上满足条件 (j,k=0,1,…,n) 则称这n+1个n次多项式l0(x), l1(x),…, ln(x)为节点x0 ,