数值分析4祥解.pptVIP

} 例 给定区间[0,3]上3个点的函数值f(0)=0,f(1)=2,f(3)=4, 试求数 a,b,c,d,使函数 S(x)为给定点上的三次样条插值函数。 解 设 根据定义，由 得 d=0，故 则 由 得 由 得 由 得 由 得 求得 } 给定n+1个样点(xi , yi )(i=0,1,…,n) ，确定一个三次样条插值函数需要4n个独立条件。在定义中，已指定了4n-2个条件，即 所以，一般需补充指定2个边界条件。 } 4.4.2 三次样条插值函数的求法 三弯矩法：以各基点处的二阶导数作为待定参数。 在区间[xk , xk+1]上，对二阶导函数S″(x)利用拉格朗日线性插值公式得 由泰勒公式有 其中 用x=xk+1 代入上式可解得 令 } 得到 S(x) 在[xk , xk+1 ]上的表达式(记为Sk(x)) S(x) 在[xk-1 , xk ]上的表达式为 根据 经过一些推导整理可得到 其中 } 这是含有n+1个未知量m0 , m1 ,…,mn 的n-1个方程的方程组，要确定未知量还需补充两个边界条件。 常用的边界条件有三种: (1)第一种边界条件：给定端点处的一阶导数，即 (2)第二种边界条件：给定端点处的二阶导数m0 ,mn。特别取m0 = mn=0时称为自然边界条件，求得的S(x)称为自然样条函