2012数值分析 -lec16--样条插值.pptVIP

数值分析 朱立永 北京航空航天大学 数学与系统科学学院 Email: numerical_analysis@ Password:beihang2012 答疑时间：星期三下午2:00－5:00 答疑地点：主216 第十六讲样条插值 第五章插值与逼近 代数插值 一元函数插值（一元Lagrange插值） 二元函数插值（二元Lagrange插值） Hermite插值 分段低次插值 样条插值 样条插值 分段低次插值，收敛性好，但光滑性不够理想。为了得到光滑度更高的插值函数，引入样条插值函数。 “样条”名词来源于工程中船体和汽车等的外形设计：给出外形曲线上的一组离散点(样点),(xi , yi),i = 0, 1, 2, …, n, 将有弹性的细长木条或钢条(样条)在样点上固定,使其在其它地方自由弯曲,这种样条所表示的曲线,称为样条曲线(函数). 这样，整个曲线不仅通过样点，并且在整个区间上其一阶导数，二阶导数是连续的。 样条函数的定义：对于区间[a, b]上的一个分划 π：a=x0x1…xn-1xn=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 称为样条节点，其中x1, …, xn-1 称为内节点，x0 和xn 称为边界节点。 样条插值是分段低次插值函数； 实际应用中常取k=3，即三次样条函数。 三次样条插值 余下的n+3个条件的确定： （1）n+1个插值节点条件，即s3(xk)=f(xk)=yk; （2）两个边界条件！ 三次样条插值的边界 （3）基于B样条的三次样条插值函数 （2）三转角方法 三弯矩法求三次样条插值函数 我们接下来的任务是确定s(x) ，也就是要确定它的待定系数，如何来确定，有多种方法。在这里，我们选取s(x)在插值节点的二阶导数为待定参数来列方程组（称为三弯矩法），也可以用其他方法，如用一阶导数作为待定参数。 1. 对于第一种边值条件有： 2. 对于第二种边值条件有：（同上） 3. 对于第三种边值条件有： * * * 例 已知离散点: (1.1, 0.4000), (1.2, 0.8000), (1.4, 1.6500), (1.5, 1.8000), 取自然边界条件 M0 = Mn = 0, 构造三次样条插值函数，并计算 f(1.25). 解 n = 3. ∵ h0= x1- x0 = 0.1, h1= 0.2, h2=0.1, 因此，分段的三次样条插值函数为 例 已知离散点: (1.1, 0.4000), (1.2, 0.8000), (1.4, 1.6500), (1.5, 1.8000), 取自然边界条件 M0 = Mn = 0, 构造三次样条插值函数，并计算 f(1.25). 解 n = 3. ∵ h0= x1- x0 = 0.1, h1= 0.2, h2=0.1, 因此，分段的三次样条插值函数为 B（皮埃尔·贝塞尔（Pierre Bézier））样条 样条函数插值 定义：记 称为 k 次半截单项式，并规定 性质： x+k 有 k-1阶连续导数； x+k 在x=0处k阶导数不存在； 小结 多项式插值法，其目的是利用节点上的值，构造通过这些节点的多项式，从原则上说，利用n+1个节点的值，可以构造n次多项式，而且这种构造是唯一的！利用待定系数法，将节点值带入后，得到一个n+1阶线性方程，即可求出多项式。 为便于在计算机上实现，本章引入了Lagrange插值公式和牛顿插值公式，各有千秋， Lagrange插值公式便于理解和记忆，牛顿公式便于计算机计算。但必须指出的是，不管用什么方法插值，所得到的插值公式实际上是完全一样的，包括上面所说的待定系数法，这就是插值公式唯一性。唯一性的一个直接推论就是各种插值公式的余项完全一样，都是Lagrange余项。 Lagrange公式的构造思想非常重要！后面的埃尔米特公式出发点也是利用了这一点。 如果节点上不仅已知函数值，同时还已知函数的导数值，即要求插值多项式在节点上与函数具有相同的函数值和导数值，这时要用埃尔米特公式。 高次插值有时会引起较大误差，对此，解决的方案是分段低次插值。如果对插值多项式要求较好的光滑性，这时就需要进行样条函数插值。 作业 教材第146页习题：20、22、25、26