Fast linear matrix interpolation for time data NOT equally spaced

3 次查看(过去 30 天)
I desperately need a fast linear matrix interpolation function for NOT equally spaced time data. Matlab's interp1 is slow and your contribution ScaleTime, albeit very fast, only accepts equally spaced time data (right?). Any suggestions?
Thanks,
Fredrik

采纳的回答

Jan
Jan 2013-8-15
编辑:Jan 2013-8-15
I assume, you mean FEX: ScaleTime. The new version is still in the pipeline. But this works efficiently for not equally spaced data as long as X is sorted:
function Yi = myLinearInterp(X, Y, Xi)
% X and Xi are column vectors, Y a matrix with data along the columns
[dummy, Bin] = histc(Xi, X); %#ok<ASGLU>
H = diff(X); % Original step size
% Extra treatment if last element is on the boundary:
sizeY = size(Y);
if Bin(length(Bin)) >= sizeY(1)
Bin(length(Bin)) = sizeY(1) - 1;
end
Xj = Bin + (Xi - X(Bin)) ./ H(Bin);
Yi = ScaleTime(Y, Xj);
A full M-version can be found here.
  1 个评论
Spigge
Spigge 2013-8-16
It works great, 5 x faster than interp1 using my matrix data. This means I can do the calculations I need without reducing the framerate in my application. Many thanks!
Fredrik

请先登录,再进行评论。

更多回答(0 个)

类别

Help CenterFile Exchange 中查找有关 Interpolation 的更多信息

产品

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by