第7章2三次样条插值.pptVIP

第四节 三次样条插值 样条插值的提出： 拉格朗日高次插值多项式可保证插值曲线光滑, 但计算不稳定。 分段低次插值虽然计算稳定, 但又不保证曲线光滑。 希望的插值函数：采用分段光滑插值，又要保证在插值节点处导数连续，插值函数的曲线是光滑的。 样条插值 样条插值：采用分段光滑插值。在每两个插值节点的小区间内用低次插值，但又要保证在插值节点处导数连续，也就是在整个插值区间内，插值函数的曲线是光滑的。 三次样条插值函数 定义：已知在[a,b]上的各点xi (i=0,1,2,…,n) 处的 yi (i=0,1,2,…,n)值, 则要求构造函数S(x)使其满足 1) S (xi)=yi (i=0,1,2,…,n)， 2)在[a,b]上，S (x)具有连续的二阶导数，即 ??????????? 3)在子区间[xi-1, xi]( i=0,1,2,…,n)上，Si (x)为不高于三次的多项式。 则称S (x)为关于结点 (xi, yi) (i=0,1,2,…,n)的三次样条插值函数, S (x)则为分段三次多项式。 三次样条插值函数的计算 三次样条插值函数S(x)，即n个Si(x)的确定 Si(x) 在子区间[xi-1, xi]( i=1,2,…,n)上为不高于三次的多项式， 在 n个子区间，每个三次多项式有四个待定系数， 共有 4n个待定系数。 三次样条插值的边界条件 求解方程（4n-2)个： 边界条件 第一类边界条件： 给定端点处的一阶导数， 第二类边界条件：给定端点处的二阶导数, 特别是 ，为自然边界条件 第三类边界条件：周期性条件 三次样条插值函数的构造 根据不同的边界条件，利用(n+1)个节点作不同的构造 节点处的二阶导数Mi表示的三次样条插值函数 节点处的一阶导数mi表示的三次样条插值函数 构造过程：假设已知Mi 或mi ，将 表示为Mi或mi的线性方程，利用样条函数的要求，求解线性方程组得到Mi或mi ， 得到 的表达式 例：以二阶导数表示三次样条 假设在各节点处的二阶导数值为 例：以二阶导数表示三次样条 对 积分，分别得到： Mi的求解方程 第一类边界条件， Mi的求解 三次样条函数的收敛误差 一般分别收敛到 MATLAB三次样条插值 MATLAB提供的三次样条插值函数有spline与 interp1两个， spline函数的调用格式 y0=spline(x,y,x0) interp1函数的调用格式： y0=interp1(x,y,x0,’spline’) 约束条件： MATLAB样条工具箱 MATLAB样条工具箱提供的函数，可用于构造样条函数以及对样条函数进行各种操作其调用格式为 pp=csape(x,y,?conds?,[valconds]) 其中(x,y) 为插值点的序列，pp为指定条件下以(x,y)为插值点所返回的pp形式的三次样条函数。conds为边界条件 “conds”边界条件 complete —— 给定端点的斜率，斜率大小在valconds参数中给出； not-a-knot —— 两个端点存在三阶连续导数； periodic —— 给定周期特性； second —— 给定端点的二阶导数，大小在valconds参数中给出； variational —— 给定端点的二阶导数，且大小均为零。 ppval函数可以计算在节点处x样条函数pp的值，调用格式为ppval(pp,x) 常用语句pp,ppval,unmkpp pp=spline(x,y); yi=ppval(pp,xi) yi=spline(x,y,xi),计算插值 [breaks,coefs,l,k]=unmkpp(pp)，调样条函数 pp=mkpp(breaks,coefs), 生成多项式 第五节 数值微分 数值微分是用离散方法近似地计算y= f (x)在某节点x=a的导数值。 根据导数定义，可以用差商近似微商（导数），有 其中h为步长。式7-24和7-25分别称为前差公式和后差公式。如果将两式平均，得中心差商。 差商的几何意义 数值微分的计算误差 数值微分的计算误差 为了提高精度可以用二次插值函数代替曲线f (x)， 三点公式，其误差为 数值微分的其他问题 二阶数值微分， 高阶多项式的数值微分 非节点处的数值微分，误差难以估计 利用三次样条插值函数求导数 三次样条插值函数S(x)作为表格函数f (x)的近似表达式，函数值非常接近，而且导数值也很接近