16_数值分析5_3_1样条插值.pptVIP

分段低次插值函数都有一致收敛性，但光滑性较差，对于像高速飞机的机翼形线，船体放样等型值线往往要求有二阶光滑度，即有二阶连续导数。早期工程师制图时，把富有弹性的细长木条（所谓样条）用压铁固定在样点上，在其它地方让它自由弯曲，然后画下长条的曲线，称为样条曲线。从数学上看，它实际上是由分段三次曲线拼接而成，在连接点处即样点上不仅函数自身是连续的，而且它的一阶导数和二阶导数也是连续的，所以样条曲线具有非常好的光滑性，从数学上加以概括就得到数学样条这一概念。 1946年舒恩柏格（I.J.Schoenberg) 首次提出了样条函数这一名称，但一开始样条函数并没有得到应有的重视，只是在电子计算机的迅速发展和样条函数也找到了很多重要的应用以后，它才得到了蓬勃发展，建立了完备的理论。自上世纪60年代以来，样条函数越来越流行，它不仅是现代函数逼近的一个活跃的分支，而且也是现代数学计算中一个十分重要的数学工具。它已经以各种方式应用到逼近论、数据拟合、数值微分、数值积分、微分方程和积分方程等数值求解中。尽管样条函数内容十分丰富，应用非常广泛，但最常用的是“三次样条插值函数”。 一、函数组在区间上线性无关 ?定义 设函数组 ??0?x?, ?1?x?, ?, ?n?x?? 在区间 [a, b]上连续,如果 a0?0?x? ? a1?1?x? ??? an?n?x? ? 0 (a ? x ? b) ? a0 ? a1 ? ? ? an ? 0 . 则称 ??k?x?? 在 [a, b] 上是线性无关的; 否则, 叫线性相关. 二、在区间[a, b]上线性无关的函数系 ?定义 若 ??i ?x?| i= 0,1, ? ? 中任何有限个函数, 在区间[a, b] 上都是线性无关的, 则称 ??i ?x?| i= 0,1,?? 是线性无关的函数系. 三、k 次半截单项式 ?定义 记 为k次半截单项式. 并规定 四、k 次半截单项式的性质 (1) 在区间(? ?, ?) 上有 k ?1阶连续的导数； (2) 当 k ? 2 时, 有 (xk+)(r) = k (k ?1)… (k ?r ?1) xk-r+；(r=1,2,…,k-1) (3) 在 x = 0处, xk+ 的 k 阶导数不存在。 对任意的实数 a, 有 五、k 次多项式样条函数 ?定义 对于区间[a, b]上的一个分划 ? : a ? x0 ? x1 ? ? ? xn ? b 如果函数 s(x) 满足条件 (1) s(x) 在每个子区间 [xi , xi+1] (i=0,1,…,n-1)上是次数不高于 k 的多项式; (2) s(x) 在 (a, b) 上有 k ?1 阶连续导数. 则称 s(x) 是定义在 [a, b] 上, 对应分划? 的 k 次多项式样条函数 (简称 k 次样条); x0, x1, ?, xn称为样条节点; x0, xn叫做边界节点, 其它点叫内节点. ? k 次样条函数 s(x) 是次数不超过 k 的多项式的线性组合, 它所在的线性空间, 记为 Dk,? . 可以证明 Dk,? ?Span{1, x,…, x k, (x ? x1) k+ ,…, (x ? xn ?1) k+}. 定义在 [a,b] 上, 对应分划? 的样条函数s(x)，可以表示为 ?注 样条函数定义中涉及到的一些问题 k 次样条函数 s(x) 在每个子区间上都是一个 k 次多项式， 它有 k ?1 个系数； 分划? 共有 n 个子区间， 所以 共有 n(k ?1 ) 个待定参数。 样条函数定义中仅提供了 k (n-1) 个约束条件： s(x) 及其导数 (阶数 k) 在内节点处连续。即， s(r)(xi-) ? s(r)(xi+) (i ? 1, 2,…, n?1; r ?0,1,…,k?1). 所以对应分划? 的 k 次样条函数组成了一个函数集， 而且它是一个线性空间， 记为Dk,? ，维数不超过 n+k ? n(k ?1) ? k(n?1) . ? 根据样条函数定义, 易知 x j?Dk,? , j ?0,1,…,k; (x ? xj )k+ ?Dk,? , j ?1,2,…, n?1; 函数组 {x j ，(x ? xj )k+}是Dk,? 的一组基。 即线性空间 Dk,? ?Span{1, x,…, x k, (x ? x1) k+ ,…, (x ? xn ?1) k+ }, 于是，定义在[a,b]上对应