File Exchange

image thumbnail

Dynamic Time Warping (DTW)

version 1.4 (3.18 KB) by

This package implements Dynamic Time Warping (DTW).

15 Ratings



View License

DTW is implemented in both Matlab and C/MEX. The C/MEX function is very fast. The algorithm we use is the standard one from wiki:
There is a demo in the package. To use, just call:
d=dtw(s,t,w) or d=dtw_c(s,t,w)

Comments and Ratings (28)


Peng (view profile)

Yongjun Wang

x = -10:0.1:10;
y = sin(x);
len = length(y);
y2 = interp1(1:len,y,linspace(1,len,2*len));


Zheng (view profile)

Hi Quan,

I'm very new to dynamic time warping. Does "D" indicate the similarity between two vectors?



Zheng (view profile)

Beil Zheng


What is a common value for the w parameter?

yabing hu

good job

Yijun Yan

Hi Quan,

Just wonder, which is the corresponding paper for this code?

Hi Quan,
It is nice. It would be useful if you can add an option to force the alighment of start/end points of two signals. Or force the alignment of selected points (so it becomes interpolation in a sense).

Chanmann Lim

Hi Quan, This is a great piece of work and I have made slight changes to normalize the DTW distance by its warping path for both matlab and c versions for my project. Please let me know if I can contribute that here.


Grace Pan

Hi,Quan Wang
Could i ask you a question? The code works perfectly when two input data have the same length. But if two input data have different length, such as [10 12 14 18 14 12 10 10 12 14 18 14 12 10] and [NaN 12 14 18 14 12 10 10 12 14 18 NaN NaN NaN], the output is NaN. Could you think how to deal with it?

nguyen anh

Hi Quan Wang.
I have a question:
What the paper does algorithm DTW using adapt window size belong to ?

Thank you.

Haoyi Xu

Tipu Rahaman

I am getting error as
'Maximum variable size allowed by the program is exceeded.'

for below line.

D=zeros(ns+1,nt+1)+Inf; % cache matrix

Quan Wang

Quan Wang (view profile)

Hi Yongqiang,

The window parameter w in function dtw(s,t,w) is the parameter that prevents big gradients.


Hi Quan,

Could I ask you a question? I just want to confirm with you, so that I correctly understand what is written in your file. You didn't implement the slope constraint (the one that prevents too much or too little gradient, namely the fifth condition in Sakoe & Chiba 1978), did you?

Sakoe & Chiba 1978: Dynamic programming algorithm optimization for spoken word recognition.

Chang hsiung


Quan Wang

Quan Wang (view profile)

Hi Daniel,

You need to learn what is a MATLAB function first.


Daniel (view profile)

What does MATLAB Compiler do? When I try to run the code, I get an error message for Line 1 (function d=dtw(a,b,w)).


Quan Wang

Quan Wang (view profile)

Hi Steven du,

Yes, you are right. I forgot to update it when I extended it for multi-dimensional signals.

Steven du

I thinks line 19 of dtw.m
Ought to be ns=size(s,1)

As length(s)=max(size(s))


great job!


Blaise (view profile)

works perfectly!

Yan Ou

Yan Ou (view profile)


Siqi (view profile)

Quan Wang

Quan Wang (view profile)

Hi Leonel Rozo, it now supports multidimensional signals.

Leonel Rozo

Does it consider multidimensional signals?



Minor bug fixed according to the suggestion from Steven du, NTU. Thanks.


added Matlab compiler as required products


removed mex file

MATLAB Release
MATLAB 8.0 (R2012b)

Download apps, toolboxes, and other File Exchange content using Add-On Explorer in MATLAB.

» Watch video