讲稿 第四章插值与拟合(1～3节).ppt

Confidential, for review onlyBorland * 各阶差商可按下表排列方式逐列进行计算，称下表为差商表 Confidential, for review onlyBorland * 实用中常采取等距节点 xi+1 - xi = h，利用此特点，引入差分概念，可简化计算 1 分段线性插值 分段线性插值就是通过插值节点用折线段连接起 来逼近f(x)。 设f(x)在n+1个节点 上的函数值为 ,在每个小区间 (k=0,1,…，n）上作线性插值，得 在几何上就是用折线 替代曲线,如右图所示。 若用插值基函数表示, 则在?a,b?上 其中 显然， 是分段线性连续函数，且 称S(x)为f(x)的分段线性插值函数。 由线性插值的余项估计式知,f(x)在每个子段 上有误差估计式 其中 例17 已知f(x)在四个节点上的函数值如下表所示 30 45 60 90 1 求f(x)在区间?30,90?上的分段连续线性插值函数S(x) 解 将插值区间?30,90?分成连续的三个小区间 ?30,45?,?45,60?,?60,90? 则S(x)在区间?30,45?上的线性插值为 S(x)在区间?45,60?上的线性插值为 S(x)在区间?60,90?上的线性插值为 将各小区间的线性插值函数连接在一起，得 2 三次样条插值 高次插值不仅计算复杂,而且可能出现Runge现象. 采用分段插值虽然计算简单、也有一致收敛性,但不能保证整条曲线在连接点处的光滑性: 如分段线性插值,其图形是锯齿形的折线,虽然连续,但插值节点处都是“尖点”,因而一阶导数都不存在,这在实用上,往往不能满足某些工程技术的高精度要求。 生产中〔如船体、飞机等外形曲线的设计中〕不仅要求曲线连续,而且要有二阶光滑度,即有连续的二阶导数。这就要求分段插值函数在整个区间上具有连续的二阶导数。 定义5 .设函数定义在区间?a, b?上，给定n+1个 节点和一组与之对应的函数值，若函数 满足: （1）在每个节点上满足 S(xi)=f(xi)(i=0,1,…,n) （2）在?a, b?上有连续的二阶导数 （3）在每个小区间? xi,xi+1? (i=0,1,…,n-1) 上是一个三次多项式。 则称S(x)为三次样条插值函数。 其中四个待定系数为 ；子区间共有n个，所以要确定S(x)需要4n个待定系数。 另一方面,要求分段三次多项式S(x)及其导数 和 在整个插值区间?a,b?上连续,则要求它们在各个子区间的连接点 上连续，即满足条件 由样条函数的定义可知,三次样条插值函数S(x)是一 个分段三次多项式,要求出S(x),在每个小区间?xi,xi+1?上 要确定4个待定参数,若用Si(x)表示它在第i个子区间?xi,xi+1?上的表达式，则 （1）?????? 插值条件 （2）?????? 连接条件 (4.35) 式(4.35)共给出了4n-2个条件,而待定系数有4n个,因此还需要2个条件才能确定S(x),通常在区间端点 上 各加一个条件,称为边界条件, 常用边界条件有三种类型。 第一种类型：给定两端点f(x)的一阶导数值： 第二种类型：给定两端点f(x)的二阶导数值： 作为特例, 称为自然边界条件。 满足自然边界条件的三次样条插值函数称为自然样条插值函数。 第三种类型：当f(x)是以 为周期的函数时，则要求S(x)也是周期函数,这时边界条件应满足 当 时， 由上给定的任一种边界条件加上插值条件和连接条件，就能得出4n个方程，可以惟一确定4n个系数。从而得到三次样条插值函数S(x)在各个子区间?xi , xi+1?上的表达式S(xi）(i=1,2,…,)。