舉個例子
clear;clc;
x=1:4;
y=sin(x);
p=spline(x,y)
運行後可以看到結構體p如下
p =
form: 'pp'
breaks: [1 2 3 4]
coefs: [3x4 double]
pieces: 3
order: 4
dim: 1
它是壹個3段多項式分段函數(pp),斷點為[1 2 3 4](即1~2、2~3、3~4三段),每個分段函數為壹個多項式,多項式的系數保存在coefs裏,每壹行代表壹個分段函數。
輸入p.coefs
ans =
0.1177 -0.7711 0.7212 0.8415
0.1177 -0.4180 -0.4679 0.9093
0.1177 -0.0649 -0.9508 0.1411
就可以知道多項式的表達式了,比如第二行表示函數:0.1177*(x-2)^3-0.4180*(x-2)^2 -0.4679*(x-2)+0.9093(這個2就是斷點)
這個分段函數的定義域為2~3,如果妳要算2.5的插值,只要把他代入這個函數,就得到了這個多項式即可,當然用matlab只要輸入命令ppval(p,2.5)就可以了