东北大学数值分析课件第六部分插值与逼近学习课件.pptVIP

当x?[xi-1,xi]时,由于 所以 于是有 由连续性条件S??(xi-0)=S??(xi+0)可得 两侧同除以 3(?i?[xi-1,xi]+?i?[xi,xi+1])=gi,则有 ?imi-1+2mi+?imi+1=gi , i=1,2,…,n-1. (6.9) 再结合不同的边界条件,可得关于mi的方程. 若边界条件为:m0=y?0 ,mn=y?n ,带入(6.9)式可得 若边界条件为:S??(x0)=y??0 ,S??(xn)=y??n ,则有 连同(6.9)式一起,可得 若边界条件为周期性边界条件，由S?(x0+0)=S?(xn-0) ,和S??(x0+0)=S??(xn-0),有 m0=mn ?nmn-1+2mn+?nm1=gn 和 其中 于是有 对应不同的边界条件，只要求出相应的线性方程组的解，便得到三次样条函数在各区间[xi-1,xi]上的表达式. 由于三个方程组的系数矩阵都是严格对角占优矩阵,所以都有唯一解,前两个方程组均可用追赶法求解,第三个方程组可用LU分解法或Gauss消元法求解. 例7 设?(0)=1,?(1)=0 ,?(2)=-1 ,?(3)=0 ,??(0)=1, ??(3)