版本 (18.3 KB) 作者: Jan
Fast linear interpolation of equally spaced data (C-MEX and M)

2.8K 次下载

更新时间 2020/10/20


ScaleTime - fast linear matrix interpolation

Yi = ScaleTime(Y, T)
Where T is a vector with values between 1 and size(Y,1). This is equivalent to
Yi = interp1(1:size(Y, 1), Y, T, 'linear')
If T is defined by uniform steps: T = linspace(Ti, Tf, Tn) this is faster:
Yi = ScaleTime(Y, Ti, Tf, Tn)

For nun-uniform input use the 2nd function:
Yi = ScaleTimeNU(X, Y, Xi)

I have recorded timeseries of e.g. size [100'000 x 10] and cut out chunks of about 100 frames according to some events and normalize them to 101 frames.
The direct approach with Matlab's INTERP1 works, but it is slow:
interp1(1:size(Y, 1), Y, Ti, '*linear')

Therefore I've created the function ScaleTime for:
- linear interpolation
- equally spaced input data
- equally or non-equally spaced interpolation steps
- matrix input
- no handling of NaNs, no extrapolation.

To my surprise, ScaleTime is remarkably faster than INTERP1 and even the modern griddedInterpolant.
The functions lininterp1f and qinterp1 from the FEX do not accept matrices and they are significantly slower than the MEX ScaleTime:
If you have installed these tools, their speed is compared in uTest_ScaleTime.

Included: M-source, C-source.
Compilation: Automatically at running ScaleTime the first time, or:
mex -O ScaleTime.c
Pre-compiled mex function coming soon on:


Jan (2023). ScaleTime (, MATLAB Central File Exchange. 检索来源 .

MATLAB 版本兼容性
创建方式 R2018b
Windows macOS Linux
Help CenterMATLAB Answers 中查找有关 Time Series Objects 的更多信息

Community Treasure Hunt

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

Start Hunting!
版本 已发布 发行说明

Bug fix in test function, double and single input. A 2nd function for treating non-uniform inputs is added.