计算机图形学补充04：插值与拟合.pptxVIP

插值与拟合插值法插值法是一种古老的数学方法，早在一千多年前的隋唐时期定制历法时就广泛应用了二次插值。刘焯将等距节点的二次插值应用于天文计算。插值理论却是在17世纪微积分产生后才逐步发展起来的，Newton插值公式理论是当时的重要成果。由于计算机的使用以及航空、造船、精密仪器的加工，插值法在理论和实践上都得到进一步发展，获得了广泛的应用。§4.1 引言§4.2 拉格朗日插值§4.3 均差与牛顿插值公式§4.4 差分与等距节点插值§4.5 埃尔米特(Hermite)插值与分段插值§4.6 曲线拟合§4.1 引言 问题的提出函数解析式未知,通过实验观测得到的一组数据, 即在某个区间[a, b]上给出一系列点的函数值 yi= f(xi)或者给出函数表xx0x1x2……xnyy0y1y2……yny=p(x)y=f(x)插值法的基本原理设函数 y=f(x) 定义在区间[a, b]上, 是[a, b]上取定的n+1个互异节点,且在这些点处的函数值 为已知 ,即 若存在一个f(x)的近似函数 ,满足则称 为f(x)的一个插值函数, f(x)为被插函数, 点xi为插值节点, 称(4.1)式为插值条件, 而误差函数R(x)= 称为插值余项, 区间[a, b]称为插值区间, 插值点在插值区间内的称为内插, 否则称外插 (4.1)插值函数 在n+1 个互异插值节点(i=0, 1, …, n )处与相等,在其它点x就用的值作为f(x)的近似值。这一过程称为插值，点x称为插值点。换句话说, 插值就是根据被插函数给出的函数表“插出”所要点的函数值。用 的值作为f(x)的近似值,不仅希望能较好地逼近 f(x),而且还希望它计算简单 。由于代数多项式具有数值计算和理论分析方便的优点。所以本章主要介绍代数插值。即求一个次数不超过n次的多项式。 满足 则称P(x)为f(x)的n次插值多项式。这种插值法通常称为代数插值法。其几何意义如下图所示 是函数 在区间[a, b]上的n+1个互异的节点 (i=0,1,2,…,n )上的插值多项式,则求插值多项式P(x)的问题就归结为求它的系数 (i=0,1,2,…,n )。 由插值条件: (i=0,1,2,…,n),可得 定理1 n次代数插值问题的解是存在且唯一的 证明: 设n次多项式 称为Vandermonde（范德蒙）行列式，因xi≠xj（当i≠j），故V≠0。根据解线性方程组的克莱姆（Gramer）法则，方程组的解 存在惟一，从而P(x)被惟一确定。 这是一个关于待定参数 的n+1阶线性方程组,其系数矩阵行列式为 唯一性说明，不论用何种方法来构造，也不论用何种形式来表示插值多项式,只要满足插值条件(4.1)其结果都是相互恒等的。 §4.2 拉格朗日（Lagrange）插值 为了构造满足插值条件 (i=0,1,2,…,n )的便于使用的插值多项式 L(x),先考察几种简单情形,然后再推广到一般形式。（ 线性插值与抛物插值）（1）线性插值线性插值是代数插值的最简单形式。假设给定了函数f(x)在两个互异的点的值，, , ,现要求用线性函数 近似地代替 f(x)。选择参数a和b, 使 。称这样的线性函数L1(x) 为 f(x) 的线性插值函数 。线性插值的几何意义:用通过点 和的直线近似地代替曲线 y=f(x)由解析几何知道,这条直线用点斜式表示为 为了便于推广，记 这是一次函数,且有性质 与 称为线性插值基函数。且有 于是线性插值函数可以表示为与基函数的线性组合 例4.1 已知 , , 求 解: 这里x0=100，y0=10，x1=121，y1=11, 利用线性插值 （2）抛物插值 抛物插值又称二次插值，它也是常用的代数插值之一。设已知f(x)在三个互异点x0，x1，x2的函数值y0，y1，y2，要构造次数不超过二次的多项式使满足二次插值条件：这就是二次插值问题。其几何意义是用经过3个点 的抛物线 近似代替曲线 ,如下图所示。因此也称之为抛物插值。 L2(x)的参数直接由插值条件决定，即 满足下面的代数方程组： 该三元一次方程组的系数矩阵 该行列式是范德蒙行列式，当时，方程组的解唯一。 为了与下一节的Lagrange插值公式比较,仿线性插值,用基函数的方法求解方程组。先考察一个特殊的二次插值问题： 求二次式 ,使其满足条件： 这个问题容易求解。由上式的后两个