MATLAB一维插值加快运行速度的方法

经常使用MATLAB的人也许会注意到MATLAB自带的一维插值函数vq = interp1(Xvec,Vvec,xq,method)调用时间过长的问题。了解他的算法可以发现，它是对xq对于Xvec向量从头到尾依次比较大小而定位到特定index的，即

for i=1:length(Xvec)-1
	if xq<Xvec(i)&&xq>Xvec(i+1)
		在[Xvec(i),Xvec(i+1)]区间内线性插值;
	else
		继续比较大小，以定位xq，确定要用哪一段的线性插值;
	end
end

如果Xvec向量是均匀排布的，这样比较大小来定位index就太慢了，可以利用xq在Xvec中的大小反算出index，即

lx=length(Xvec);% Xvec向量的长度
index=(Xvec(lx)-Xvec(1))/xq*lx；

然后在index和index+1之间进行线性插值即可。