计算方法_04.ppt

* ⒊ Ⅲ型样条函数 对于Ⅲ型插值问题，由插值条件(4-42)知： S?(x0+0)=S?(xn-0) S??(x0+0)=S??(xn-0) 上述第一式表明： m0=mn 而由第二式以及(4-45)、(4-46)可得： 化简得 其中： 5.6 三次样条插值 * 将式(4-57)与(4-50)合并，即得： 5.6 三次样条插值 * 四、求样条插值函数的三弯矩法 三次样条插值函数S(x)还可用节点上的二阶导数值S ??(xi)=Mi (i=0,1,…,n)为参数进行表示，这种通过确定Mi来求S(x)的方法称为三弯矩法。 因为S(x)在[xi,xi+1]上是三次多项式，所以S ??(x)在[xi,xi+1]上是一次多项式，且可表示为： 对S ??(x)积分两次并利用S(xi)=yi和S(xi+1)=yi+1定出积分常数得： 5.6 三次样条插值 * 为了使S(x)成为所求的三次样条，同时为了确定参数Mi，对S(x)求导得： 所以 5.6 三次样条插值 * 由S?(xi+0)=S?(xi-0) (i=1,2,…,n-1)得： 其中?i、?i (i=1,2,…,n-1)由式(4-48)表示，而di为： 方程组(4-63)与方程组(4-50)完全类似，它是关于n+1个未知数Mi (i=0,1,…,n)的n-1个方程，为确定Mi必须再加2个条件。 5.6 三次样条插值 * ⒈ Ⅰ型样条函数 对Ⅰ型插值问题，已知S ?(x0)=m0， S ?(xn)=mn，由公式(4-61)、(4-62)得： 即： 其中： 5.6 三次样条插值 * 合并(4-63)、(4-65)有： 从中解出Mi (i=0,1,…,n)，代入式(4-60)得到Ⅰ型三次样条S(x)。 5.6 三次样条插值 * ⒉ Ⅱ型样条函数 对于Ⅱ型样条插值问题，M0，Mn已知，故方程组(4-63)为： 从中解出Mi (i=1,2,…,n-1)，代入式(4-60)得到Ⅱ型三次样条S(x)。 5.6 三次样条插值 * ⒊ Ⅲ型样条函数 对于Ⅲ型插值问题，已知S?(x0+0)=S?(xn-0)，则由公式(4-61)、(4-62)得： 整理得： 其中： 5.6 三次样条插值 * 综合式(4-63)、(4-69)有： 从中解出Mi (i=1,2,…,n)，代入式(4-60)得到Ⅲ型三次样条S(x)。 5.6 三次样条插值 * 五、余项估计及收敛性、稳定性 设f(x)在[a,b]上有直到四阶的连续导数，则其Ⅰ型和Ⅱ型三次样条插值函数以及导数的误差有如下估计式： 其中， 5.6 三次样条插值 * * 二、 * * * 次多项式H3(x)，使其满足： 这样的H3(x)称为三次埃尔米特插值多项式。 二、埃尔米特插值多项式 解决上述问题，类似于拉格朗日插值多项式的构造过程，仍采用基函数的思想，首先求分别满足下列插值条件的四个最基本的三次埃尔米特插值多项式，即基函数?0(x)、?1(x)、?0(x)、?1(x) 5.4 埃尔米特插值 * 5.4 埃尔米特插值 * 5.4 埃尔米特插值 * 现仅以求?0(x)为例进行推导。由式(4-18)可知， ?0(x1)=? ?0(x1)=0 即x1为三次多项式?0(x)的二重根，所以可将?0(x)表示为： ?0(x)=[c(x-x0)+d](x-x1)2 式中常数c、d待定。在由式(4-18)知 解得 5.4 埃尔米特插值 * 类似地，可以构造出其他基函数 5.4 埃尔米特插值 * 根据条件式(5.21) ~ (5.24)，可得满足条件式(5.20)的三次埃尔米特插值多项式为： H3(x)=y0?0(x)+y1?1(x)+m0?0(x)+m1?1(x) (5.29) 实际计算中经常用到下面两个函数： ?0(x)=(1+2x)(1-x)2 ?1(x)=x(1-x)2 (5.30) 它们分别满足： ?0(0)=1 ?0(1)=?0?(0)=?0?(1)=0 ?1?(0)=1 ?1(0)=?1(1)=?1?(1)=0