样条插值的算例三次样条的概念用一阶导数表示的样条三次样条的极性.pptVIP

* 样条插值的算例 三次样条的概念 用一阶导数表示的样条 三次样条的极性 《数值分析》 16 ? ? ? ? 例1. 飞机机翼剖面图 1.数据采集 2/17 2. 数据插值处理 X 0 -0.4552 -0.6913 -0.8640 -0.9689 -0.9996 Y 0 0.3285 0.3467 0.2716 0.1408 -0.0160 T=1:6;t=1:.2:6; x=spline(T,X,t); y=spline(T,Y,t); 例2:龙格函数的插值逼近 3/17 7结点等距插值 7结点切比雪夫插值 7结点样条插值 7结点埃尔米特插值 利用龙格函数的数据表做样条插值第一步 4/17 x -5.0 -3.33 -1.66 0 1.66 3.33 5.0 y 0.038 0.082 0.264 1.0 0.264 0.082 0.038 m 0.014 -0.0054 0.4142 0. -0.4142 0.0054 -0.014 y’ 0.014 0.045 0.233 0 -0.233 -0.045 -0.014 估算结点处导数值 mk, 由三对角方程组求解得出 (k=1，2，3，4，5，6)， h=10/6 定义 5.4 给定区间[a , b]上的一个分划: a = x0 x1 … xn = b 已知 f(xj) = yj (j = 0,1,···,n), 如果 满足: (1) S(x)在 [xj，xj+ 1]上为三次多项式; (2) S”(x)在区间[a，b]上连续; (3) S(xj) = yj ( j = 0,1,···,n). 则称 S(x)为三次样条插值函数. 5/17 当x∈[xj , xj+ 1] ( j= 0,1,…n-1 )时 Sj(x)= aj + bj x + cj x2 + dj x3 插值条件: S(xj) = yj ( j = 0,1,···,n) 连续性条件: S(xj+0) =S(xj-0) ( j = 1,···,n-1) S’(xj+0) =S’(xj-0) ( j = 1,···,n-1) S”(xj+0) =S”(xj-0) ( j = 1,···,n-1) 由样条定义,可建立方程(4n-2)个！！ n个三次多项式, 待定系数共4n个！！ 方程数少于未知数个数 ？？ 6/17 (1)自然边界条件: S”(x0)=0, S”(xn)=0 例 5.7 已知f(–1) = 1, f(0) = 0, f(1) = 1.求[–1，1] 上的三次自然样条(满足自然边界条件). 解 设 则有: – a1+b1–c1+d1=1,d1=0, a2+b2+c2+d2=1 d1=d2, c1=c2, b1=b2 (2)周期边界条件: S’(x0)=S’(xn), S”(x0)=S”(xn) (3)固定边界条件: S’(x0)=f ’(x0), S’(xn)=f ’(xn) 7/17 由自然边界条件: – 6a1+2b1=0, 6a2+2b2=0 解方程组,得 a1=-a2=1/2, b1=b2=3/2, c1=c2=d1=d2=0 问题的解 y= x2 8/17 y= S (x) 分段Hermite插值公式导出的样条方法 已知函数表 x x0 x1 ······ xn f(x) y0 y1 ······ yn 设 f(x) 在各插值节点 xj 处的一阶导数为 mj 取 xj+1 – xj = h，( j = 0,1,2,···,n).当 x∈[xj, xj+ 1]时, 分段Hermite插值 9/17 由S”(x)连续,有等式: S”(xj + 0)=S”(xj – 0) 考虑 S”(x) 在区间[xj , xj+1]和[xj-1 , xj]上表达式. 当 x∈[xj , xj+1]时, S(