5.5 样条插值.ppt

* 24 * 24 * 24 第五节、样条插值 问题概述 三次样条 构造三次样条 算法设计 样条类型 回顾前面几节讲过的各种代数插值,它们有一个共同的弱点,那就是: 它们都是相当刚性(stiff)的.也就是说, 局部数据误差易向远处传播、放大. 如果采用分段多项式插值，则由于插值基函数只是局部活跃，结点上的误差可以被控制在小的范围内，这是分段插值的一大优势。 一、样条函数的物理背景 许多实际问题希望插值函数具有较高阶的整体光滑性，此时高次埃尔米特插值或分段高次埃尔米特插值可以利用，但提高埃尔米特插值多项式的次数就要增加约束条件—给出插值点被插函数及其直到足够高阶导数之值。 作为约束条件的所有数据都是通过观测得到的。观测总难免有误差。高次插值不仅增添了数据准备与计算困难，也将导致更大的误差。 高次埃尔米特 还有许多应用不仅要求插值函数具有足够高阶的整体光滑性，还要求在某些结点处转折灵活。例如若赶不上点处加载集中力的杆、梁或板弯曲。这就导致本节要讨论的样条函数插值。 数学里的样条一词来源于它的直观几何背景：绘图员或板金工人常用弹性木条或金属条加压铁来绘制或者放样成曲线、曲面。 序 样条：是指飞机或轮船等的制造过程中为描绘出光滑的外形曲线(放样)所用的工具。 对数据集进行多项式曲线拟合在CAD,CAM和计算机图形系统中有许多应用，操作者希望画出经过数据点的无误差光滑曲线。 从传统上讲，一般使用曲线板或设计师的样条主观地画出曲线，只要通过眼睛看着是光滑的就可以了。 二、样条曲线 从数学角度上分析，在每个区间[xk,xk+1]上可以构造一个三次函数sk(x)。使得分段曲线y=s(x)和它的一、二阶导数在更大区间[x0,xn]内连续。 一阶导函数连续意味着曲线没有急弯。 二阶导函数连续意味着足够光滑且每点的曲率半径有定义。 数学样条 设 有n+1个点，其中 若存n个三次多项式 满足如下性质 k=0,1,…,n k=0,1,…,n-2 k=0,1,…,n-2 k=0,1,…,n-2 k=0,1,…,n-2 则称函数s(x)为三次样条函数 三次样条函数 是否可能构造一个三次样条满足以上五条性质？ 每个三次多项式Sk(x)有4个未知常数，需要求解4n个系数，即要确定4n个自由度。 为了保证样条插值问题求解的唯一性，至少需要4n个条件。 分析性质2-5，插值条件共有n＋1个，内点原函数和一、二阶导函数连续条件各有n-1个，总共4n-2个条件。 还剩两个自由度，实际应用中常采用区间端点处满足某种边界条件。 三次样条的存在性 由于S(x)是分段三次多项式，它的二阶导数 区间[x0,xn]内是分段线性的，由线性拉格朗日插值有 构造三次样条 令 代入，有 将上次积分两次，引入两个积分常数，得如下形式 将xk和xk+1代入，则有 从而得 代回方程 表达式化为只含{mk}形式，再将xk点代入方可求出未知{mk}，求导 而一阶导函数在xk点连续,即 从而得 方程组未知数为{mk}，而 {hk}、{dk}、{uk}通数据点集简单计算即可得到。方程组包含n+1个未知数mk ，具有n-1个线性方程的欠定方程组。 需在端点处找两个方程才能消去第一个方程中的m0和最后一个方程中的mn 求解出mk后，各段三次多项式的系数为 端点处确实两方程的不同方法对应不同的插值类型 各段系数 k=1,2,…，n-1 样条是柔软有弹性的木杆经过所有数据点形成的曲线，但让端点的斜率自由的地某一位置保持平衡，使得曲线的摇摆最小。 m0=0 mn=0 它在对多位有效数字精度的试验数据进行曲线拟合时很有用。 自然样条 例 已知(0,0),(1,0.5),(2,2.0),(3,1.5), 求三次自然样条样条。 解： 自然样条，m0=mn=0,代入方程组有 即： 解之得 紧压样条在端点有斜率。 S‘ (a),S’ ( b) 已知 紧压样条可想像为，用外力使柔软而有弹性的木杆经过数据点，并在端点处使其具有固定斜率 紧压样条 例 已知(0,0),(1,0.5),(2,2.0),(3