第三章 插值方法与数据拟合.docVIP

第三章 插值方法与数据拟合 所讨论的问题给复杂的函数找一简单的函数如多项式、三角函数等，并让其满足一定的条件，让其近似的取代原函数。 或 有一数据表格，我们需要找一函数取近似的表征该表数据。 §1 拉格朗日（Lagrange）插值 在函数类中多项式具有最简单的性质。 设在区间[a，b]连续的实函数已知在该区间上n+1个不同点的函数值 或 有数据表有对数据 插值节点 我们需要找一个n次多项式 使得在这些点上函数值等于插值节点的值。 1、线性插值 已知两个点的函数值： 做一线性函数使得在两个节点上函数值为节点值。 函数为： 基函数为一次函数，且在节点上 几何意义：过两点做直线。按变化量平均。 2、抛物线插值 已知三个点的函数值： 做二次函数使得在三个节点上函数值为节点值。 函数为： 基函数为二次函数，且在节点上 3、拉格朗日插值 已知n+1个点的函数值： 做n次函数使得在n+1个节点上函数值为节点值。 函数为： 基函数为n次函数，且在节点上 为了程序设计我们作如下推导： 程序结构是什么样子？ DIMENSION X(N),Y(N) 读入节点数据 X=..... T1=0.0 DO 100 I=0,N T2=1.0 DO 200 J=0,N IF(J.EQ.I) GOTO 200 T2=T2* 200 CONTINUE T1=T1+T2* 100 CONTINUE WRITE(*,*)X,T1 STOP END 拉格朗日插值的唯一性 已知n+1个点的函数值： 做n次函数 使得在n+1个节点上函数值为节点值。 即： 求得系数则多项式确定 线性方程组，系数矩阵对应的行列式： 范德蒙（Vandermonde）行列式 所以方程有唯一解。所以多项式唯一。 §2、牛顿插值公式 差分、差商 定义 4.1 称为函数的的一阶差商，为的二阶差商。 一般地， 为的n阶差商。为统一起见，补充定义为零阶差商。 由差商定义，显然有 如此类推，可以证明 由此看出，差商的值与节点的排列次序无关，即 这种性质称为差商对称性。 定义4.2 设函数在等距节点上的函数值为，其中为常数，称做步长。称 向前差分 向后差分 中心差分 分别为在处以为步长的一阶向前差分，一阶向后差分和一阶中心差分。 符号、、分别称为向前差分算子，向后差分算子和中心差分算子。 由一阶差分的定义出发，可定义二阶差分 ， ， 。 一般地刻定义阶差分为 如果令，，则称为不变算子，为移位算子。由于 ， 于是 。 同理可得 ， 牛顿插值公式 有了差商的概念，前面介绍的线性插值与抛物插值可表示为 ， 抛物插值 事实上，按差商定义 ， ， 反复用后一个式子代入前面的式子，可得 记 （*） （**） 于是 由于 ， 故必有 所以式（*）为次插值多项式，由插值多项式的唯一性讨论可知 ， 且有 。 式（*）称作牛顿基本插值公式，它的各项系数就是函数的各阶差商，每增加一个插值节点，只需在原来的基础上多计算一项，这一性质称作牛顿插值公式的承继性。 利用牛顿插值公式时，一般分两步： 第一步：造差商表 节点 函数值 一阶差商 二阶差商 三阶差商 。。。。 X0 X1 X2 X3 X4 X5 … Xn F(x0) F(x1) F(x2) F(x3) F(x4) F(x5) … F(xn) F(x0,x1) F(x1,x2) F(x2,x3) F(x3,x4) F(x4,x5) … F(xn-1,xn) F(x0,x1,x2) F(x1,x2,x3) F(x2,x3,x4) F(x3,x4,x5) … F(xn-2,xn-1,xn) F(x0,x1,x2,x3) F(x1,x2,x3,x4) F(x2,x3,x4,x5) … F(xn-3,xn-2,xn-1,xn) 第二步：利用 例如：节点表 i 0 1 2 3 4 5 Xi 0.0 1.0 2.0 3.0 4.0 5.0 F(xi) 0.0 2.0 12.0 42.0 116.0 282.0 差商表 Xi F（xi) 一阶 二阶 三阶 四阶 五阶 1 2 3 4 5 6 0 2 12 42 116 282 2 10 30 74 166 4 10 22 46 2 4 8 0.5