clc,clear all
x = [1 2 3 4 5];
v = [12 16 31 10 6];
xq = [0 0.5 1.5 5.5 6];
vq1 = interp1(x,v,xq,'pchip') %外插值
vq2 = interp1(x,v,xq,'linear') %内插值
vq3 = interp1(x,v,xq,'linear','extrap') %内插值和外插值结合
V1=[v,vq1];
X=[x,xq];
V2=[v,vq2];
V3=[v,vq3];
[X,ind]=sort(X);
for i=0:length(X)
V1(ind);
V2(ind);
V3(ind);
end
plot(x,v,'*:k')
hold on
plot(X,V1,'.-r')
plot(X,V2,'o:b')
plot(X,V3,'p-g')
legend('原始数据','外插值','线性插值','内外插值结合')