东北大学数值分析 第6章 插值与逼近.ppt

当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)=0,试求?(x)在区间[0,3]的三次样条插值函数S(x). 解 这里h1=h2=h3=1 ,y?0=1, y?3=0,计算参数有 ?1=?2=?1=?2=1/2, g1=-3,g2=0 于是有 解得 故有 三次样条函数S(x)也可以利用在节点处的二阶导数为参数来表示,设S??(xi)=Mi,i=0,1,…,n,则对x?[xi-1,xi]有 连续积分两次,并利用S(xi-1)=yi-1 ,S(xi)=yi ,确定积分常数,可得 其中hi=xi-xi-1 .为了确定S(x),只需确定Mi ,i=0,1,…,n.可利用S?(xi-0)=S?(xi+0)来求出Mi . 对上式求导易得: 于是有 因此 若记 则有 再结合不同的边界条件,可得关于Mi的方程. 若边界条件为:M0=y??0 ,Mn=y??n ,可得 若边界条件为:S?(x0)=y?0 ,S?(xn)=y?n ,则有 可得 若边界条件为周期性边界条件，由S?(x0+0)=S?(xn-0) ,和S??(x0+0)=S??(xn-0),有 M0=Mn ?nM1+?nMn-1+2Mn=dn 和 其中 于是有 而且M0=Mn . 例8 设?(0)=0,?(1)=1 ,?(2)=0 ,?(3)=1 ,???(0)=1, ???(3)=0,试求?(x)在区间[0,3]的三次样条插值函数S(x). 解 这里h1=h2=h3=1 ,y??0=1, y??3=0,计算参数有 ?1=?2=?1=?2=1/2, d1=-6,d2=6 于是有 解得 故有 §6 正交多项式 记区间[a,b]上所有连续函数的全体为C[a,b],可以证明C[a,b]是一个线性空间,把所有次数不超过n的多项式全体记为Pn,则Pn是C[a,b]的子空间. 若?(x),g(x)?C[a,b], 则称 为?(x)与g(x)的内积,记为(?,g),满足 (1)(?,g)=(g,?); (2)(c?,g)=c(?,g); (3)(?1+?2,g)=(?1,g)+(?2,g); 若(?,g)=0,称?(x)与g(x)正交,记为??g . 利用内积可以定义函数的平方模 函数的平方模满足 (1) ?????2?0,而且?????2=0??(x)=0; (2) ??c???2=|c|?????2; (3) ???+g??2??????2+??g??2 (4) ?(?,g)???????2 ??g??2 考虑到?(x)在区间[a,b]上各点的函数值比重不同,常引进加权形式的定义 这里函数?(x)是非负连续函数,称为[a,b]上的权函数.它的物理意义可以解释为密度函数.