matlab求解三次条样插值

函数原型 csape(x,y,conds,valconds)

须要安装Curve Fitting Toolbox这个工具箱

conds 表示是哪一种边界条件

　　'not-a-knot'  非扭结条件,不用给边界值.

        'periodic'  周期性边界条件,不用给边界值.

        'second'  给定边界二阶导数.

　　 'variational'  天然样条(边界二阶导数为0)，默认

　　 'complete'  给定导数状况

valconds表示给定的条件

eg：

 

x	27.7	28	29	30
y	4.1	4.3	4.1	3.0

 

 

 

 

 

边界条件：

 

S'(27.7) = 3.0 ，S'(30) = -4.0 

 

 

x = [27.7 28 29 30];
y = [4.1 4.3 4.1 3.0]
pp=csape(x,y,'complete',[3,-4]);
disp(pp.coefs);
xi=27.7:0.05:30;
yi=ppval(pp,xi);
plot(x,y,'o',xi,yi);

 

运行结果：

y =

　　4.1000 4.3000 4.1000 3.0000

　　13.2930 -11.7657 3.0000 4.1000
　　0.0723 0.1980 -0.4703 4.3000
　　-1.6574 0.4149 0.1426 4.1000

 

补充：除了csape还有spline函数能实现天然边界条件和第二边界条件及能够实现一维或者高维的曲线插值。