三次样条插值2.pptx

;三弯矩算法; 只要能求出所有的{M i}，就能求出样条插值函数S(x).;; 第一型边界条件：;化简得;;说明; 例1 已知函数f(x)的数值表如下： x 2 4 6 f(x) 3 7 13 1 -1 试求f(x) 在[2,6]上的三次样条插值函数 ;由公式知α2 =β0 =1 ，; 得方程组 2 M0 + M1 = 3 0.5 M0 + 2M1 +0.5 M2 = 1.5 M1 +2 M2 = -12;Maths程序如下：;b[1]=1; b[j_]:=1-a[j]; c[1]=6/h[1]((y[2]-y[1])/h[1]-y[1]); c[j_]:=6((y[j+1]-y[j])/h[j]-(y[j]-y[j-1])/h[j-1])/(h[j-1]+h[j]); c[3]=6/h[3-1](y[3]-(y[3]-y[3-1])/h[3-1]); A=Table[Switch[i-j,-1,b[j-1],0,2,1,a[j+1],_,0],{i,1,3},{j,1,3}]; MatrixForm[%] CC=Table[c[j],{j,1,3}]; MatrixForm[%] LinearSolve[A,CC]; MatrixForm[%];;M[j_]:=LinearSolve[A,CC][[j]] Table[M[j],{j,1,3}] S[j_]:=M[j+1](x-x[j])^3/(6h[j])-M[j](x-x[j+1])^3/(6h[j])+ (y[j+1]-M[j+1]h[j]^2/6)(x-x[j])/h[j]- (y[j]-M[j]h[j]^2/6)(x-x[j+1])/h[j] Table[S[j],{j,1,2}]; Expand[%]; MatrixForm[%];g1=Plot[%[[1]],{x,2,4}] g2=Plot[%%[[2]],{x,4,6}] g3=ListPlot[B,Prolog-AbsolutePointSize[15]] Show[g1,g2,g3,Prolog-AbsolutePointSize[15]];这种情况下只有n-1个未知数，其矩阵形式为：;第三型边界条件：已知f(x)是以b -a为周期的周期函数 ，要求S(x)满足周期条件; （1）确定边界条件，判定是第几型插值问题；;;对条件(1)有;方程组为：;对条件(2)有;　　构造一阶导数值 表示的三次样条插值函数。 在力学上解释为细梁在 截面处的转角，并且得到的转角与相邻两个转角有关，故称用 表示 的算法为三转角算法。;　　根据Hermite插值函数的唯一性和表达式 可设 S(x)在区间[xi , xi+1](i=0,1,…n-1)的表达式为 ;于是有;其中;由此可解得m1,m2,…, mn-1 ，从而得 S(x)的表达式.; (2)对于边界条件(a),可导出两个方程: ;若令;可解出 方程组的矩阵形式为;　　在实际应用中，如果不需要规定内节点处的一阶导数值，那么使用三次样条插值函数会得到很好的效果。三次样条插值函数s(x)不仅在内节点处的二阶导数是连续的，而且s(x)逼近f(x)具有很好的???敛性，也是数值稳定的。由于误差估计与收敛性定理的证明比较复杂，下面只给出误差估计的结论。;