
clc,clear

%% 数据
x = [27.7 28 29 30];
y = [4.1 4.3 4.1 3.0];
xx = 27.7:0.01:30;

%% 三次样条插值
% coef = spline(x,[3.0, y, -4.0]); % 方法一
coef = csape(x,[3.0, y, -4.0],[1, 1]); % 方法二
Y_spline = ppval(coef,xx);


%% 画图对比
hold on;
plot(xx,Y_spline,'b-','Linewidth',1);
plot(x,y,'ro','Linewidth',2)

legend('spline','data')

axis([27.5 30 3 5])