经常使用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之间进行线性插值即可。