dtw mathematica

gistfile1.mathematica

dtwMatrix2[s_List, t_List] :=
 With[{n = Length@s, m = Length@t},
  Fold[
   {acc, pos} \[Function] pos /. {i_, j_} :>
      ReplacePart[acc,
       pos -> (s[[i - 1]] - t[[j - 1]])^2 + 
         Min[Part[acc, ##] & @@@ {{i - 1, j}, {i, j - 1}, {i - 1, 
             j - 1}}]],
   SparseArray[{{1, 1} -> 
      0, {1, _} -> \[Infinity], {_, 1} -> \[Infinity]}, {n, m} + 1],
   Tuples[{Range[n], Range[m]} + 1]]]