第二章 函数基本逼近(一).ppt

* ? 构造三次样条插值函数的三弯矩法 在 上，记 ] , [ for ) ( ) ( 1 ] [ j j j x x x x S x S - ? = 对每个j, 此为3次多项式 则 S[j]”(x) 为 次多项式，需 个点的值确定之。 1 2 设 S[j]”(xj?1) = Mj?1， S[j]”(xj) = Mj 对应力学中的梁弯矩，故名 对于x ?[xj?1, xj ] 可得到 S[j]”(x) = j j j j j j h x x M h x x M 1 1 - - - + - 积分2次，可得 S[j]’(x) 和 S[j](x) ： j j j j j j j A h x x M h x x M + - + - - - - 2 ) ( 2 ) ( 2 1 2 1 S[j]’(x) = j j j j j j j j B x A h x x M h x x M + + - + - - - 6 ) ( 6 ) ( 3 1 3 1 S[j](x) = 利用已知S[j](xj?1) = yj?1 S[j](xj) = yj 可解 * j j j j j j j h M M h y y A 6 1 1 - - - - - = j j j j j j j j j j j j h x x h M y h x x h M y B x A 1 2 2 1 1 ) 6 ( ) 6 ( - - - - - + - - = + 利用已知S[j](xj?1) = yj?1 S[j](xj) = yj 可解 * 下面解决 Mj ： 利用S’ 在 xj 的连续性 [xj?1, xj ]: S[j]’(x) = j j j j j j j j j j j h M M x x f h x x M h x x M 6 ] , [ 2 ) ( 2 ) ( 1 1 2 1 2 1 - - - - - - + - + - - 1 1 1 1 2 1 1 2 1 6 ] , [ 2 ) ( 2 ) ( + + + + + + + - - + - + - - j j j j j j j j j j j h M M x x f h x x M h x x M [xj , xj+1]: S[j+1]’(x) = 利用S[j]’(xj) = S[j+1]’(xj)，合并关于Mj?1、 Mj、 Mj+1的同类项，并记 , , , 整理后得到： 1 1 j j j j h h h + + + = l 1 j j - = l m ]) , [ ] , [ ( 6 1 1 1 j j j j j j j x x f x x f h h g - + + - + = 2 1 1 g M M M j j j j j j = + + + - l m ? j ? 1 n?1 即：有 个未知数， 个方程。 n?1 n+1 还需2个边界条件 /* boundary conditions */ * ? 第1类边界条件 ： S’(a) = y0’， S’(b) = yn’ [a , x1 ]: S[1]’(x) = 1 0 1 1 0 1 2 1 1 2 1 0 6 ] , [ 2 ) ( 2 ) ( h M M x x f h a x M h x x M - - + - + - - 0 1 0 1 1 0 ) ] , [ ( 6 2 g y0’ x x f h M M = - = + n n n n n n g x x f yn’ h M M = - = + - - ]) , [ ( 6 2 1 1 类似地利用[ xn?1, b ] 上的 S[n]’(x) * ? 第2类边界条件： S”(a) = y0” = M0， S”(b) = yn” = Mn 这时： * ? 第3类边界条件： 当 f 为周期函数时， yn = y0 ， S’(a) = S’(b) ， S” (a) = S”(b) ?M0 = Mn * 例 给定插值条件 0 1 2 3 0 0 0 0 附加边界条件为： .试求 满足上述条件的三次插值样条函数的分段表达式. 解：（三弯矩方法） 利用二阶导数 为待定参数,得 * 边界条件方程为: * 利用公式 求得三次插值样条函数： * 若利用三转角方法,选取 作待定 参数，也可获