详细讲解三次样条插值法与其实现方法.pptVIP

三次样条插值;样条函数的定义;三次样条插值函数的定义;三次样条插值函数的边界条件;（1）因为s(x)在每个小区间上是一个次小于三次的多项式，故有四个未知系数； （2）因为s(x)有n分段，从而共有4n个未知系数！ （3）但插值条件与样条条件仅给出4n-2个条件，无法定出4n个未知系数，还差2个条件！这2个条件我们用边界条件给出！;通常我们对插值多项式在两端点的状态加以要求也就是 所谓的边界条件:;第二边界条件：由区间端点处的二阶导数给出即 ;这样三次样条插值问题就分成三类！其实不止这三类!;样条函数的例子;样条函数的例子;通常有三转角法、三弯矩法、B样条基函数法。;三转角法:待定一阶数;我们采用待定一阶导数的方法即设;拣拒货逊夫厚著富蚤己娩犊波果鉴织郁俐刊唤仑怒惧甸捍汤囱霖耍瞥属怨详细讲解三次样条插值法与其实现方法详细讲解三次样条插值法与其实现方法;由于在内部节点处二阶导数连续条件：;称为三转角法基本方程组;第一类三次样条插值问题方程组;这是一个严格对角占优的三对角方程组， 用追赶法可以求解！;第二类三次样条插值问题的方程组;稍加整理得;第三类样条插值问题的方程组;其中：;求解这些不同类型的样条插值问题的方程组，我们可得 所要待定的一阶导数：;例1. 对于给定的节点及函数值;我们可以将上述计算列于表中：;由些得如下方程组：;同样可以求得第三段表达式！;乍买助贯迷裸弹违潘瓶芜正谋秀荧脊谓桓第技酥亏铅左辅暖城抗滴倦笛吓详细讲解三次样条插值法与其实现方法详细讲解三次样条插值法与其实现方法;汛孩敛峪茶陋猩瘸但哺绢首疡谬读萍朵厩出佛寅荔朋淑莹阂屉场培坟条武详细讲解三次样条插值法与其实现方法详细讲解三次样条插值法与其实现方法;三弯矩法:待定二阶导数;水催枣锈暇蹈壁摸枷绎傲枉咒个九庄面泳窃爬安碗犹并疫芋澜揽莱坛孽记详细讲解三次样条插值法与其实现方法详细讲解三次样条插值法与其实现方法;从而推导出了三次样条S(x)在第k个小区间[xk,xk+1]上的表达式为:;对所有中间节点xk,k=1,2,…,n-1,左边小区间与右边小区间上的三次多项式的一阶导数应当连续!;三弯矩法基本方程;曲率调整样条;自然样条;固支样条;固支样条;非扭结样条;非扭结样条;周期样条;再由三弯矩基本方程,可得;在Matlab中数据点称之为断点。如果三次样条插值没有边界条件，最常用的方法，就是采用非扭结（not-a-knot）条件。这个条件强迫第1个和第2个三次多项式的三阶导数相等，对最后一个和倒数第2个三次多项式也做同样地处理。 Matlab中三次样条插值也有现成的函数： y=interp1(x0,y0,x,spline)； y=spline(x0,y0,x)； pp=csape(x0,y0,conds), pp=csape(x0,y0,conds,valconds)， y=ppval(pp,x)。 其中x0,y0是已知数据点，x是插值点，y是插值点的函数值。 对于三次样条插值，我们提倡使用函数csape，csape的返回值是pp形式，要求插值点的近似函数值，必须调用函数ppval。;pp=csape(x0,y0,conds,valconds) conds指定插值的边界条件，其值可为： complete 边界为一阶导数，一阶导数的值在valconds参数中给出。 not-a-knot 非扭结条件 periodic 周期条件 second 边界为二阶导数，二阶导数的值在valconds参数中给出，若忽略valconds参数，二阶导数的缺省值为[0, 0]。;例2：第一边界条件的例题;第一边界条件的例题;如下代码求解上述样条问题：;x=[-4 -3 -2 -1 0 1 2 3 4]; y=[0 0.15 1.12 2.36 2.36 1.46 0.49 0.06 0]; pp=csape(x,y,second); xx=-4:0.01:4; yy=ppval(pp,xx); hold on; plot(x,y,ok); plot(xx,yy,k-); hold off;;pp = form: pp breaks: [-4 -3 -2 -1 0 1 2 3 4] coefs: [8x4 double] pieces: 8 order: 4 dim: 1; pp.coefs= 0.18085603829161 0.00000000000000 -0.03085603829161 0 -0.08428019145803 0.54256811487482 0.51171207658321 0.15000000000000