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函數能實現天然邊界條件和第二邊界條件及能夠實現一維或者高維的曲線插值。