cubic与spline都是Matlab的三次样条插值法,但是它们在插值点处仍然有着很微妙的区别,这个区别说明不了两种方法的好坏,只能根据实际情况进行合理筛选。以一维插值为例:
clc
clear
%% load data
load Fig1m
x1=Fig1m(:,1);y1=Fig1m(:,2);
m=x1;n=y1;
m(57)=0.0226;m(83)=0.0326;
m(95)=0.0374;m(97)=0.0378;m(99)=0.0382;
mx=0.0:0.0001:0.065;
ny=interp1(m,n,mx,‘cubic‘); % or ny=interp1(m,n,mx,‘pchip‘);
m1=log(1+mx);n1=ny.*(1+mx);
%% Cubic spline
figure(1)
H=plot(m1,n1,‘-‘,mx,ny,‘-‘,‘LineWidth‘,1.5,‘MarkerSize‘,5);
text(0.005,2.7,‘(a) cubic or pchip‘,‘interpreter‘,‘latex‘,‘fontsize‘,18);
xlabel(‘Strain‘,‘Interpreter‘,‘LaTex‘,‘fontsize‘,18);
ylabel(‘Stress (GPa) ‘,‘interpreter‘,‘latex‘,‘fontsize‘,18);
h=legend(H,‘in1‘,‘in2‘,‘Location‘,‘South‘);
po=get(h,‘Position‘ );