6讲：ch2-5样条插值(录课).ppt

内容小结 故所求的三次样条插值函数 三次样条插值函数的构造 三次样条插值函数的构造 Maths程序如下： Clear[x,y,a,b,c,n,M] x[i_]:=2i; y[1]=3; y[2]=7; y[3]=13; B=Table[{x[i],y[i]},{i,1,3}]; y[1]=1; y[3]=-1; h[j_]:=2; a[j_]:=h[j-1]/(h[j-1]+h[j]); a[3]=1; 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]); 三次样条插值函数的构造 不能单独用c 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[%] 注意它的含义 到此求出三弯矩方程的系数矩阵 解方程求出M[I] 构造三次样条插值函数 三次样条插值函数的构造 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]] 注：三对角矩阵的构造 A=Table[Switch[i-j,-1,a,0,b,1,c,_,0],{i,1,m},{j,1,n}]; MatrixForm[%] 行标减列标= -1时为对角线上方元素。 作图观察效果 三次样条插值函数的构造 三次样条插值函数的构造 程序运行结果如下： 系数矩阵 常数项 求得的弯矩值M 样条函数 样条插值图形 三次样条插值函数的构造 2、若满足S″(a)=M0 = f″(a) , S″(b)=Mn= f″(b）， 则事实上只有n-1个未知数，其矩阵形式为： 三次样条插值函数的构造 三次样条插值函数的构造 总结以上论述，可得求三次样条的步骤为： （1）确定边界条件，判定是第几类插值问题； （2）根据所确定的条件计算各值，形成方程组(**)； （3）解三对角方程组(**)，求得M0, M1, M2, ???, Mn； （4）将求得的Mi值代回S(x)的表达式中，从而可求 得函数y=f(x)在任一点的近似值S(x)。 三次样条插值函数的构造 例2．已知f(x)在若干点处的值为f(0)=0, f(1)=1, f(2)=1, f(3)=0,试求f(x)满足条件（1）f′(0)=1,f′(3)=2 （2） f”(0)=1,f”(3)=2的三次样条插值函数s(x)以及f(2.5)的近似值。 解：构造一阶均差表 例5.3 三次样条插值函数的构造 对条件(1)有 三次样条插值函数的构造 方程组为： 将数据代入可得样条插值函数 三次样条插值函数的构造 对条件(2)有 三次样条插值函数的构造 三次样条插值函数的构造 1.一般插值函数的不足； 内 容 小 结 2.三次样条函数的概念； 3.三次样条插值函数的构造. 二 元 样 条 曲 面 第五节 样条插值法 插值法 样条插值的研究背景 样条函数的力学意义 三次样条插值多项式的构造 一般的插值问题 拉格朗日(Lagrange) 插值多项式，是插值基函数的 线性组合，形式简单，结构对称，便于理论分析。 一、拉格朗日(Lagrange) 插值 样条插值的研究背景 牛顿(Newton)插值引入了差商的概念，在增加新的 节点时，只是增加一项，前面结果可再利用。 二、牛顿 (Newton) 插值 样条插值的研究背景 称为龙格Runge现象。 三、分段插值 样条插值的研究背景