Optimal Warping Paths are unique for almost
every Pair of Time Series

arXiv:1705.05681v2 [cs.LG] 2 Mar 2018

Brijnesh J. Jain and David Schultz
Technische Universität Berlin, Germany
e-mail: brijnesh.jain@gmail.com

Update rules for learning in dynamic time warping spaces are based on optimal warping
paths between parameter and input time series. In general, optimal warping paths are not
unique resulting in adverse effects in theory and practice. Under the assumption of squared
error local costs, we show that no two warping paths have identical costs almost everywhere
in a measure-theoretic sense. Two direct consequences of this result are: (i) optimal warping
paths are unique almost everywhere, and (ii) the set of all pairs of time series with multiple
equal-cost warping paths coincides with the union of exponentially many zero sets of quadratic
forms. One implication of the proposed results is that typical distance-based cost functions
such as the k-means objective are differentiable almost everywhere and can be minimized by
subgradient methods.

1. Introduction
1.1. Dynamic time warping
Time series such as audio, video, and other sensory signals represent a collection of time-dependent values
that may vary in speed (see Fig. 1). Since the Euclidean distance is sensitive to such variations, its
application to time series related data mining tasks may give unsatisfactory results [6, 8, 24]. Consequently,
the preferred approaches to compare time series apply elastic transformations that filter out the variations
in speed. Among various techniques, one of the most common elastic transformation is dynamic time
warping (DTW) [20].
Dynamic time warping is based on the concept of warping path. A warping path p determines how to
stretch two given time series x and y to warped time series x0 and y 0 under certain constraints. The cost
of warping x and y along warping path p measures how dissimilar the warped time series x0 and y 0 are.
There are exponential many different warping paths [1] each of which determines the cost of warping
time series x and y. An optimal warping path of x and y is a warping path with minimum cost. Optimal
warping paths exist but are not unique in general (see Fig. 2).

1.2. The problem of non-uniqueness
Recent research is directed towards extending standard statistical concepts and machine learning methods
to time series spaces endowed with the DTW distance. Examples include time series averaging [3, 5, 16,
17, 21], k-means clustering [10, 18, 22], self-organizing maps [14], learning vector quantization [23, 12],
and warped-linear classifiers [11, 13].
The lowest common denominator of these approaches is that they repeatedly update one or more
parameter time series. In addition, update directions are based on optimal warping paths such that the
following properties hold:
• If an optimal warping path is unique, then the update direction is well-defined.

1

Figure 1: Variations in speed and their implications for comparing two time series shown in blue and
yellow. Both time series are similar in shape but differ in speed along their flat regions. Plot
(a): The Euclidean distance warps the i-th points of both time series onto one another as shown
by the gray lines. Since the Euclidean warping is sensitive to variations in speed the shapes are
not preserved and both peaks are aligned to points in flat regions. Plot (c): Euclidean warping
leaves the time series unchanged and results in a large dissimilarity score as indicated by the
gray shaded area. Plot (b): An optimal warping that preserves the shape of both time series.
Double arrows indicate segments stretched by dynamic time warping. Plot (d): Warped time
series obtained by the optimal warping in (b). The left flat segment of the yellow and the right
flat segment of the blue time series have been stretched to align both time series. The resulting
dissimilarity is zero and better reflects the similarity in shape than the Euclidean distance.
• If an optimal warping path is non-unique, then there are several update directions.
Non-uniqueness of optimal warping paths complicates the algorithmic design of learning methods in DTW
spaces and their theoretical analysis. In some situations, non-uniqueness may result in adverse effects.
For example, repulsive updating in learning vector quantization finds a theoretical justification only in
cases where the corresponding optimal warping path is unique [12].
Given the problems caused by non-uniqueness, it is desirable that optimal warping paths are unique
almost everywhere. In this case, non-unique optimal warping paths occur exceptionally and are easier to
handle as we will shortly. Therefore, we are interested in how prevalent unique optimal warping paths are.

1.3. Almost everywhere
The colloquial term “almost everywhere” has a precise measure-theoretic meaning. A measure quantifies
the size of a set. It generalizes the concepts of length, area, and volume of a solid body defined in one,
two, and three dimensions, respectively. The term “almost everywhere” finds its roots in the notion of a
“negligible set”. Negligible sets are sets contained in a set of measure zero. For example, the function

1 : x 6= 0
f (x) =
0 : x=0
is discontinuous on the negligible set {0} with measure zero. We say, function f is continuous almost
everywhere, because the set where f is not continuous is negligible. More generally, a property P is
said to be true almost everywhere if the set where P is false is negligible. The property that an optimal
warping is unique almost everywhere means that the set of all pairs of time series with non-unique optimal
warping path is negligible.
When working in a measure space, a negligible set contains the exceptional cases we can handle or
even do not care about and often ignore. For example, we do not care about the behavior of the above
function f on its negligible set {0} when computing its Lebesgue integral over [−1, 1]. Another example
is that cost functions of some machine learning methods in Euclidean spaces such as k-means or learning
vector quantization are non-differentiable on a negligible set. In such cases, it is common practice to
ignore such points or to resort to subgradient methods.

2

Figure 2: Example of non-unique warping paths (see Section 2.1 for definitions). Section (a) Time series
x and y. (b) Local cost matrix with elements (xi − yj )2 . (c)√Accumulated cost obtained by
dynamic programming [20] gives the DTW distance δ(x, y) = 15. The gray shaded cells show
two different optimal warping paths. (d) Alignments of x and y by the two different optimal
warping paths.

1.4. Contributions
Consider the following property P on the set F m × F n of pairs of time series of length m and n: The pair
of time series x ∈ F m and y ∈ F n satisfies P if there are two different (not necessarily optimal) warping
paths between x and y with identical costs. Under the assumption of a squared error local cost function,
the main result of this article is Theorem 1:
Property P is negligible on F m × F n .
Direct consequences of Theorem 1 are (i) optimal warping paths are unique almost everywhere, and (ii)
property P holds on the union of exponentially many zero sets of quadratic forms. The results hold for
uni- as well as multivariate time series.
An implication of non-unique optimal warping paths is that adverse effects in learning are exceptional
cases that can be safely handled. For example, learning amounts in (stochastic) gradient descent update
rules almost everywhere.

2. Background
This section first introduces warping paths and then defines the notions of negligible and almost everywhere
from measure theory.

2.1. Time Series and Warping Paths
We first define time series. Let F = Rd denote the d-dimensional Euclidean space. A d-variate time series
of length m is a sequence x = (x1 , . . . , xm ) consisting of elements xi ∈ F. By F m we denote the set of all
time series of length m with elements from F.
Next, we describe warping paths. Let [n] = {1, . . . , n}, where n ∈ N. An (m × n)-lattice is a set of
the form Lm,n = [m] × [n]. A warping path in lattice Lm,n is a sequence p = (p1 , . . . , pL ) of L points
pl = (il , jl ) ∈ Lm,n such that
1. p1 = (1, 1) and pL = (m, n)
2. pl+1 − pl ∈ {(1, 0), (0, 1), (1, 1)} for all l ∈ [L − 1].
The first condition is called boundary condition and the second one is the step condition.

3

By Pm,n we denote the set of all warping paths in Lm,n . A warping path departs at the upper left
corner (1, 1) and ends at the lower right corner (m, n) of the lattice. Only east (0, 1), south (1, 0), and
southeast (1, 1) steps are allowed to move from a given point pl to the next point pl+1 for all 1 ≤ l < L.
Finally, we introduce optimal warping paths. A warping path p ∈ Pm,n defines an alignment (warping)
between time series x ∈ F m and y ∈ F n by relating elements xi and yj if (i, j) ∈ p. The cost of aligning
time series x and y along warping path p is defined by
X
2
Cp (x, y) =
kxi − yj k ,
(i,j)∈p

where k·k denotes the Euclidean norm on F. A warping path p∗ ∈ Pm,n between x and y is optimal if
p∗ ∈ argmin Cp (x, y).
p∈Pm,n

By P∗ (x, y) we denote the set of all optimal warping paths between time series x and y. The DTW
distance is defined by
q
δ(x, y) = min
Cp (x, y).
p ∈ Pm,n
p
is the DTW distance. Observe that δ(x, y) = Cp (x, y) for all p ∈ P∗ (x, y).

2.2. Measure-Theoretic Concepts
We introduce the necessary measure-theoretic concepts to define the notions of negligible and almost
everywhere. For details, we refer to [9].
One issue in measure theory is that not every subset of a given set X is measurable. A family A
of measurable subsets of a set X is called σ-algebra in X . A measure is a function µ : A → R+ that
assigns a non-negative value to every measurable subset of X such that certain conditions are satisfied.
To introduce these concepts formally, we assume that P(X ) denotes the power set of a set X , that is the
set of all subsets of X . A system A ⊂ P(X ) is called a σ-algebra in X if it has the following properties:
1. X ∈ A
2. U ∈ A implies X \ U ∈ A
S
3. (Ui )i∈N ∈ A implies i∈N Ui ∈ A.
A measure on A is a function µ : A → [0, +∞] that satisfies the following properties:
1. µ(U) ≥ 0 for all U ∈ A
2. µ(∅) = 0
3. For a countable collection of disjoint sets (Ui )i∈N ∈ A, we have
!
[
X
µ
Ui =
µ (Ui ) .
i∈N

i∈N

A triple (X , A, µ) consisting of a set X , a σ-algebra A in X and a measure µ on A is called a measure space.
The Borel-algebra B in Rd is the σ-algebra generated by the open sets of Rd . The Lebesgue-measure µ
on B generalizes the concept of d-volume of a box in Rd . The triple Rd , B, µ is called Borel-Lebesgue
measure space.
Let (X , A, µ) be a measure space, where X is a set, A is a σ-algebra in X , and µ is a measure defined
on A. A set N ⊂ X is µ-negligible if there is a set N 0 ∈ A such that µ(N 0 ) = 0 and N ⊆ N 0 . A property
of X is said to hold µ-almost everywhere if the set of points in X where this property fails is µ-negligible.

4

3. Results
We first show that optimal warping paths are unique almost everywhere. Then we geometrically describe
the location of the non-unique set. Finally, we discuss the implications of the proposed results on learning
in DTW spaces.
Let X = F m × F n bet the set of all pairs (x, y) of time series, where x has length m ∈ N and y has
length n ∈ N. We regard the set X as a Euclidean space and assume the Lebesgue-Borel measure space
(X , B, µ).1 The multi optimal-path set of X is defined by
NX∗ = {(x, y) ∈ X : |P∗ (x, y)| > 1} .
This set consists of all pairs (x, y) ∈ X with non-unique optimal warping path. To assert that NX∗ is
µ-negligible, we show that NX∗ is a subset of a set of measure zero. For this, consider the multi-path set
[
NX =
{(x, y) ∈ X : Cp (x, y) = Cq (x, y)},
p,q∈Pm,n
p6=q

The set NX consists of all pairs (x, y) that can be aligned along different warping paths with identical
cost. Obviously, the set NX∗ is a subset of NX . The next theorem states that NX is a set of measure zero.
Theorem 1. Let (X , B, µ) be the Lebesgue-Borel measure space. Then µ(NX ) = 0.
From Theorem 1 and NX∗ ⊆ NX immediately follows that NX∗ is µ-negligible.
Corollary 2. Under the assumptions of Theorem 1 the set NX∗ is µ-negligible.
Thus, optimal warping paths are unique µ-almost everywhere in X . Even more generally: The property
that all warping paths have different cost holds µ-almost everywhere in X .
We describe the geometric form of the multi-path set NX . For this we identify F m × F n with F k ,
where k = m + n. Thus, pairs (x, y) ∈ X of time series are summarized to z ∈ F k , henceforth denoted as
z ∈ X . By X 2 = F k×k we denote the set of all (k × k)-matrices with elements from F. Finally, the zero
set of a function f : X → R is of the form
Z(f ) = {z ∈ X : f (z) = 0} .
From the proof of Theorem 1 directly follows that the set NX is the union of zero sets of quadratic forms.
Corollary 3. Under the assumptions of Theorem 1, there is an integer D ∈ N and symmetric matrices
A1 , . . . , AD ∈ X 2 such that
D
[

NX =
Z zT Ai z .
i=1

The number D of zero sets in Corollary 3 grows exponentially in m and n. From the proof of Theorem
1 follows that D = Dm,n (Dm,n − 1)/2, where
min {m,n}

Dm,n =

X

  
m n
2
.
i
i
i

i=0

is the Delannoy number [1]. The Delannoy number Dm,n = |Pm,n | counts the number of all warping
paths in lattice Lm,n . Table 1 presents the first Delannoy numbers up to m = 10 and n = 9. We see
that there are more than half a million warping paths in a (10 × 9)-lattice showing that NX is the union
of more than 178 billion zero sets. For two time series of length 20, the number of warping paths is
D20,20 = 260, 543, 813, 797, 441, which is more than 260 trillions. Thus, the multi-path set NX of two
time series of length 20 is the union of more than 33 octillion zero sets, that is d > 33 · 1027 . An open
question is the number D∗ ≤ D of zero sets that form the multi optimal-path set NX∗ . The example in
Figure 3 indicates that the multi optimal-path set NX∗ can be much smaller than the multi-path set NX .
1 See

Remark 7 for an explanation of why we regard X as a Euclidean space.

5

m\n

1

2

3

4

5

6

7

8

9

1
2
3
4
5
6
7
8
9
10

1
1
1
1
1
1
1
1
1
1

1
3
5
7
9
11
13
15
17
19

1
5
13
25
41
61
85
113
145
181

1
7
25
63
129
231
377
575
833
1,159

1
9
41
129
321
681
1,289
2,241
3,649
5,641

1
11
61
231
681
1,683
3,653
7,183
13,073
22,363

1
13
85
377
1,289
3,653
8,989
19,825
40,081
75,517

1
15
113
575
2,241
7,183
19,825
48,639
108,545
224,143

1
17
145
833
3,649
13,073
40,081
108,545
265,729
598,417

Table 1: Delannoy numbers Dm,n for m ∈ [10] (rows) and n ∈ [9] (columns).

Figure 3: Illustration of the multi-path set NX and multi optimal-path sets NX∗ for the pair of time series
x = (x1 , x2 , 0) and y = (4, 8, 7) with varying values for x1 and x2 . The Delannoy number D3,3 is
13. Hence, the multi-path set NX consists of D = 78 zero sets of quadratic forms as indicated by
colored curves. The subset of multi optimal-path set NX∗ is highlighted by fat curve segments.

3.1. Discussion
We discuss the implications of Theorem 1 for learning in DTW spaces.
3.1.1. Learning
This section shows that almost-everywhere uniqueness implies almost-everywhere differentiability of the
underlying cost function.
To convey the line of argument, it is sufficient to restrict to the problem of averaging time series
as representative for other, more complex learning problems. In contrast to computing the average in
Euclidean spaces, time series averaging is a non-trivial task for which the complexity class is currently
unknown [3].
Let x(1) , . . . , x(N ) be a sample of N time series, possibly of varying length. Consider the cost function
J : F n → R of the form
N
1 X   (k) 
J(z) =
` δ x ,z ,
N
k=1

where ` : R → R is a loss function. Common loss functions for averaging time series are the identity
`(a) = a and the squared loss `(a) = a2 . The goal of time series averaging is to find a time series z∗ ∈ F n
of length n that minimizes the cost J(z).

6

The challenge of time series averaging is to minimize the non-differentiable cost function J(z). We show
that almost-everywhere uniqueness of optimal warping paths implies almost-everywhere differentiability
of the cost function J(z) and provides a stochastic (incremental) update rule.
We express the DTW distance δ(x, z) as a parametrized function. Suppose that x ∈ F m is a time series.
Then the parametrized DTW function δx : F n → R restricted to the set F n is of the form δx (z) = δ(x, z).
We have the following result:
Proposition 4. Suppose that x ∈ F m and z ∈ F n are two time series with unique optimal warping path.
Then the function δx (z) = δ(x, z) is differentiable at z and its gradient is a time series of the form
∇z δx (z) = ∇z Cp∗ (x, z) ∈ F n ,
where p∗ ∈ P∗ (x, z) is an optimal warping path.
The proof follows from [19] after reducing δx (z) to a piecewise smooth function. By construction, the
cost Cp (x, z) of warping x and z along warping path p ∈ Pm,n is differentiable as a function of the second
argument z ∈ F n . Non-differentiability of δx (z) is caused by non-uniqueness of an optimal warping path
between x and z. In this case we have
δx (z) = Cp∗ (x, z) = Cq∗ (x, z),
where p∗ , q∗ ∈ P∗ (x, z) are two distinct optimal warping paths. Then it can happen that
∇z Cp∗ (x, z) 6= ∇y Cq∗ (x, z)
showing that δx is non-differentiable at z.
Next, suppose that the loss function ` : R → R is differentiable and an optimal warping path p∗ between
time series x and z is unique. Then the individual cost Jx (z) = ` (δ(x, z)) is also differentiable at z with
gradient
∇z Jx (z) = `0 (Cp∗ (x, z))∇z Cp∗ (x, z).
Differentiability of Jx (z) gives rise to a stochastic (incremental) update rule of the form
z ← z − η `0 (Cp∗ (x, z))∇z Cp∗ (x, z),

(1)

where η is the step size and p∗ is an optimal warping path between x and z. We can also apply the
stochastic update rule (1) in cases where an optimal warping path between x and z is not unique. In this
case, we first (randomly) select an optimal warping path p∗ from the set P∗ (x, z). Then we update z by
applying update rule (1).
Updating at non-differentiable points according to the rule (1) is not well-defined. In addition, it is
unclear whether the update directions are always directions of descent, for learning problems in general.
The next result confines both issues to a negligible set.
Corollary 5. Suppose that X = F m × F n and ` : R → R is a differentiable loss function. Then the
functions
• δ(x, z) = δx (z)
• Jx (z) = ` (δ(x, z))
are differentiable µ-almost everywhere on X .
Corollary 5 directly follows from Prop. 4 together with Theorem 1. In summary, almost-everywhere
uniqueness of optimal warping paths implies almost-everywhere differentiability of the individual cost
Jx (z). The latter in turn implies that update rule (1) is a well-defined stochastic gradient step almost
everywhere.
The arguments in this section essentially carry over to other learning problems in DTW spaces such
as k-means, self-organizing maps, and learning vector quantization. We assume that it should not be
a problem to transfer the proposed results to learning based on DTW similarity scores as applied in
warped-linar classifiers [11, 13].

7

3.1.2. Learning Vector Quantization in DTW Spaces
Learning vector quantization (LVQ) is a supervised classification scheme introduced by Kohonen [15]. A
basic principle shared by most LVQ variants is the margin-growth principle [12]. This principle justifies
the different learning rules and corresponds to stochastic gradient update rules if a differentiable cost
function exists. As the k-means algorithm, the LVQ scheme has been generalized to DTW spaces [12, 23].
In this section, we illustrate that a unique optimal-warping path is a necessary condition to satisfy the
margin-growth principle in DTW spaces as proved in [12], Theorem 12.
As a representative example, we describe LVQ1, the simplest of all LVQ algorithms [15]. Let X = Rd be
the d-dimensional Euclidean space and let Y = {1, . . . , C} be a set consisting of C class labels. The LVQ1
algorithm assumes a codebook C = {(p1 , z1 ), . . . , (pK , zk )} of K prototypes pk ∈ X with corresponding
class labels zk ∈ Y. As a classifier, LVQ1 assigns an input point x ∈ X to the class zc of its closest
prototype pc ∈ C, where
2
c ∈ argmin kpk − xk .
k

LVQ1 learns a codebook C on the basis of a training set D = {(x1 , y1 ), . . . , (xN , yN )} ⊆ X × Y. After
initialization of C, the algorithm repeats the following steps until termination: (i) Randomly select a
training example (xi , yi ) ∈ D; (ii) determine the prototype pc ∈ C closest to xi ; and (iii) attract pc to xi
if their class labels agree and repel pc from xi otherwise. Step (iii) adjusts pc according to the rule
pc ← pc ± η (xi − pc ),
where η is the learning rate, the sign ± is positive if the labels of pc and xi agree (zc = yi ), and negative
otherwise (zc 6= yi ). The update rule guarantees that adjusting pc makes an incorrect classification of xi
more insecure. Formally, if the learning rate η < θ is bounded by some threshold θ, the LVQ1 update
rule guarantees to increase the hypothesis margin
(
2
2
kxi − p− k − kxi − pc k : yi = zc
,
µc (xi ) =
2
2
kxi − pc k − kxi − p+ k : yi 6= zc
where p+ (p− ) is the closest prototypes of x with the same (different) class label.
The different variants of LVQ have been extended to DTW spaces by replacing the squared Euclidean
distance with the squared DTW distance [12, 23]. The update rule is based on an optimal warping
path between the current input time series and its closest prototpye. In asymmetric learning [12], the
margin-growth principle always holds for the attractive force and for the repulsive force only when the
optimal warping path is unique (as a necessary condition).
3.1.3. Comments
We conclude this sections with two remarks.
Remark 6. Proposition 4 states that uniqueness of an optimal warping path between x and y implies
differentiability of δx at y. The converse statement does not hold, in general. A more general approach to
arrive at Prop. 4 and Corollary 5 is as follows: First, show that a function f is locally Lipschitz continuous
(llc). Then invoke Rademacher’s Theorem [7] to assert almost-everywhere differentiability of f . By the
rule of calculus of llc functions, we have:
• δx (y) is llc on F n , because the minimum of continuously differentiable functions is llc.
• If the loss ` is llc, then Jx (z) is llc, because the composition of llc functions is llc.
Remark 7. Measure-theoretic, geometric, and analytical concepts are all based on Euclidean spaces
rather than DTW spaces. The reason is that contemporary learning algorithms are formulated in such
a way that the current solution and input time series are first projected into the Euclidean space via
optimally warping to the same length. Then an update step is performed and finally the updated solution
is projected back to the DTW space. Therefore, to understand this form of learning under warping, we
study the DTW distance δ(x, y) as a function restricted to the Euclidean space F m × F n , where m is the
length of x and n is the length of y.

8

4. Conclusion
The multi-path set is negligible and corresponds to the union of zero sets of exponentially many quadratic
forms. As a subset of the multi-path set, the multi optimal-path set is also negligible. Therefore optimal
warping paths are unique almost everywhere. The implications of the proposed results are that adverse
effects on learning in DTW spaces caused by non-unique optimal warping paths can be controlled and
learning in DTW spaces amounts in minimizing the respective cost function by (stochastic) gradient
descent almost everywhere.
Acknowledgements.

B. Jain was funded by the DFG Sachbeihilfe JA 2109/4-1.

A. Proofs
The appendix presents the proof of Theorem 1 and Proposition 4. We first consider the univariate case (d = 1) in
Sections A.1 and A.2. Section A.1 introduces a more useful representation for proving the main results of this
contribution and derives some auxiliary results. Section A.2 proves the proposed results for the univariate case.
Finally, Section A.3 generalizes the proofs to the multivariate case.

A.1. Preliminaries
We assume that elements are from F = R, that is d = 1. We write Rm instead of F m to denote a time series of
length m. By ek ∈ Rm we denote the k-th standard basis vector of Rm with elements

1 : i=k
k
ei =
.
0 : i 6= k
Definition 8. Let p = (p1 , . . . , pL ) ∈ Pm,n be a warping path with points pl = (il , jl ). Then

T
Φ = ei1 , . . . , eiL ∈ RL×m

T
Ψ = ej1 , . . . , ejL ∈ RL×n
is the pair of embedding matrices induced by warping path p.
The embedding matrices have full column rank n due to the boundary and step condition of the warping
path. Thus, we can regard the embedding matrices of warping path p as injective linear maps Φ : Rm → RL and
Ψ : Rn → RL that embed time series x ∈ Rm and y ∈ Rn into RL by matrix multiplication Φx and Ψy. We can
express the cost Cp (x, y) of aligning time series x and y along warping path p by the squared Euclidean distance
between their induced embeddings.
Proposition 9. Let Φ and Ψ be the embeddings induced by warping path p ∈ Pm,n . Then
Cp (x, y) = kΦx − Ψyk2
for all x ∈ Rm and all y ∈ Rn .
Proof. [21], Proposition A.2.
Next, we define the warping and valence matrix of a warping path.
Definition 10. Let Φ and Ψ be the pair of embedding matrices induced by warping path p ∈ Pm,n . Then the
valence matrix V ∈ Rm×m and warping matrix W ∈ Rm×n of warping path p are defined by
V = ΦT Φ
W = ΦT Ψ.
The definition of valence and warping matrix are oriented in the following sense: The warping matrix W ∈ Rm×n
aligns a time series y ∈ Rn to the time axis of time series x ∈ Rm . The diagonal elements vii of the valence matrix

9

V ∈ Rm×m count the number of elements of y warped onto the same element xi of x. Alternatively, we can define
the complementary valence and warping matrix of w by
V = ΨT Ψ
W = ΨT Φ = W T .
The complementary warping matrix W ∈ Rn×m warps time series x ∈ Rm to the time axis of time series y ∈ Rn .
The diagonal elements v ii of the complementary valence matrix V ∈ Rn×n counts the number of elements of x
warped onto the same element yi of y.
Let p ∈ Pm,n be a warping path of length L with induced embedding matrices Φ ∈ RL×m and Ψ ∈ RL×n . The
aggregated embedding matrix Θ induced by warping path p is defined by
Θ = (Φ, −Ψ) ∈ RL×k ,
where k = m + n. Then the symmetric matrix ΘT Θ is of the form


V
−W
.
ΘT Θ =
V
−W
We use the following notations:
X = Rm × Rn = Rk
X 2 = Rk×k .
The next result expresses the cost Cp (x, y) by the matrix ΘT Θ.
Lemma 11. Let Θ be the aggregated embedding matrix induced by warping path p ∈ Pm,n . Then we have
Cp (z) = zT ΘT Θz
for all z ∈ X .
Proof. Suppose that Θ = (Φ, −Ψ), where Φ and Ψ are the embedding matrices induced by p. Let z = (x, y) ∈ X .
Then we have
Cp (z) = kΦx − Ψyk2
= xT ΦT Φx − xT ΦT Ψy − yT ΨT Φx + yT ΨT Ψy
= xT V x − xT W y − yT W x + yT V y
 

 V
−W
x
= xT , yT
y
−W
V
= zT ΘT Θz.

The last auxiliary result shows that the zero set of a non-zero quadratic form has measure zero.
Lemma 12. Let matrix A ∈ Rn×n be non-zero and symmetric. Then


µ x ∈ Rn : xT Ax = 0 = 0,
where µ is the Lebesgue measure on Rn .
Proof. Since A is symmetric, there is an orthogonal matrix Q ∈ Rn×n such that Λ = QT AQ is a diagonal matrix.
Consider the function
n
X
f (x) = xT Λx =
λii x2i ,
i=1

where the λii are the diagonal elements of Λ. Since A is non-zero, there is at least one λii 6= 0. Hence, f (x) is a
non-zero polynomial on Rn . Then
the set U = {x ∈ Rn : f (x) = 0} is measurable and has measure zero [4].

We show that the set Ũ = x̃ ∈ Rn : x̃T Ax̃ = 0 is also a set of measure zero. Consider the linear map
φ : U → Rn ,
First, we show that φ(U) = Ũ.

10

x 7→ QT x.

• Ũ ⊆ φ(U): Suppose that x̃ ∈ Ũ. With Qx̃ = x we have
0 = x̃T Ax̃ = x̃T QT ΛQx̃ = xT Λx.
This shows that x ∈ U . From φ(x) = QT x = x̃ follows that x̃ ∈ φ(U).
• φ(U) ⊆ Ũ: Let x̃ ∈ φ(U). Then x̃ = QT x = φ(x) for some x ∈ U. Hence, x = Qx̃ and we have
0 = xT Λx = x̃T QT ΛQx̃ = x̃T Ax̃.
This shows that x̃ ∈ Ũ.
Next, we show that µ(Ũ) = 0. Observe that the linear map φ is continuously differentiable on a measurable set
U with Jacobian Jφ (x) = QT . Applying [2], Prop. 3.7.3 gives
Z
µ (φ(U)) ≤
det QT dx.
U

Since QT is orthogonal, we have det QT = 1. Thus, we find that
Z
µ (φ(U)) ≤
dx = µ(U) = 0.
U

Finally, the assertion µ(Ũ) = 0 follows from Ũ = φ(U).

A.2. Proof of Theorem 1 and Proposition 4
This section assumes the univariate case (d = 1).
Proof of Theorem 1:
Suppose that Pm,n = {p1 , . . . , pD }. We use the following notations for all i ∈ [D]:
1.
2.
3.
4.

Θi denotes the aggregated embedding matrix induced by warping path pi .
Vi and Wi are the valence and warping matrices of pi .
V i and W i are the complementary valence and warping matrices of pi .
Ci (z) with z = (x, y) denotes the cost Cpi (x, y) of aligning x and y along warping path pi .

For every i, j ∈ [D] with i 6= j and for every z = (x, y) ∈ X , we have
Ci (z) − Cj (z) = zT ΘTi Θi z − zT ΘTj Θj z = zT A(ij) z,
where A(ij) ∈ X 2 is a symmetric matrix of the form
(ij)

A

=

ΘTi Θi

−

ΘTj Θj

=

Vi − Vj
−W i + W j

−Wi + Wj
Vi−Vj

!
.

For i 6= j the warping paths pi and pj are different implying that the warping
n matrices Wi and Wj ,oresp., are also
different. Hence, A(ij) is non-zero and from Lemma 12 follows that Uij = z ∈ X : zT A(ij) z = 0 has measure
zero. Then the union
[
U=
Uij
i<j

of finitely many measure zero sets also has measure zero. It remains to show that NX = U.
• NX ⊆ U: Suppose that z = (x, y) ∈ NX . Then there are indices i, j ∈ [D] with i < j such that the costs of
aligning x and y along warping paths pi and pj are identical, that is Ci (x, y) = Cj (x, y). Setting z = (x, y)
gives
0 = Ci (z) − Cj (z) = zT A(ij) z.
Hence, z ∈ Uij ⊆ U and therefore NX ⊆ U.
• U ⊆ NX : Let z = (x, y) ∈ U. Then there is a set Uij containing z. From Ci (z) − Cj (z) = 0 follows that pi
and pj are two warping paths between x and y with identical costs. Hence, (x, y) ∈ NX . This proves the
assertion.

11

Proof of Proposition 4:
To show the proposition, we first define the notion of piecewise smooth function. A function f : Rn → R is
piecewise smooth if it is continuous on Rn and for each x0 ∈ Rn there is a neighborhood N (x0 ) of x0 and a finite
collection (fi )i∈I of continuously differentiable functions fi : N (x0 ) → R such that
f (x) ∈ {fi (x) : i ∈ I}
for all x ∈ N (x0 ).
Proof. We show that the function δx (y) = minp Cp (x, y) is piecewise smooth. The function δx (y) is continuous,
because all Cp are continuous and continuity is closed under the min-operation. In addition, the functions Cp
are continuously differentiable as functions in the second argument. Let y0 ∈ F n and let N (y0 ) ⊆ F n be a
neighborhood of y0 . Consider the index set
I = {p ∈ Pm,n : ∃y ∈ N (y0 ) s.t. δx (y) = Cp (x, y)} .
By construction, we have δx (y) ∈ {Cp (x, y) : p ∈ I} for all y ∈ N (y0 ). This shows that δx (y) is piecewise smooth.
Then the assertion follows from [19], Lemma 2.

A.3. Generalization to the Multivariate Time Series
We briefly sketch how to generalize the results from the univariate to the multivariate case. The basic idea is to
reduce the multivariate case to the univariate case. In the following, we assume that x ∈ F m and y ∈ F n are two
d-variate time series and p = (p1 , . . . , pL ) ∈ Pm,n is a warping path between x and y with elements pl = (il , jl ).
First observe that a d-variate time series x ∈ F m consists of d individual component time series x(1) , . . . , x(d) ∈
Rm . Next, we construct the embeddings of a warping path. The d-variate time warping embeddings Φd : F n → F L
and Ψd : F m → F L induced by p are maps of the form


 
yj1
xi1

 

Φd (x) =  ...  ,
Ψd (y) =  ...  .
xiL
yjL
The maps Φd and Ψd can be written as


Φd (x) = Φx(1) , . . . , Φx(d)


Ψd (x) = Ψy (1) , . . . , Ψy (d) ,
where Φ and Ψ are the embedding matrices induced by p. Since Φ and Ψ are linear, the maps Φd and Ψd are also
linear maps. We show the multivariate formulation of Prop. 9.
Proposition 13. Let Φd and Ψd be the d-variate embeddings induced by warping path p ∈ Pm,n . Then
Cp (x, y) = kΦd (x) − Ψd (y)k2 .
for all x ∈ F m and all y ∈ F n .
Proof. The assertion follows from
kΦd (x) − Ψd (y)k2 =

d
X

2

Φx(k) − Ψy (k)

k=1

=

d
X
X 

(k)

2

(k)

2

(k)

− yj

(k)

− yj

xi

k=1 (i,j)∈p

=

d 
X X

xi

(i,j)∈p k=1

=

X

kxi − yj k 2 = Cp (x, y).

(i,j)∈p

Due to the properties of product spaces and product measures, the proofs of all other results can be carried out
componentwise.

12

References
[1] C. Banderier and S. Schwer. Why Delannoy numbers?
135(1):40–54, 2005.

Journal of Statistical Planning and Inference,

[2] V. Bogachev. Measure Theory. Springer-Verlag Berlin Heidelberg, 2007.
[3] M. Brill, T. Fluschnik, V. Froese, B. Jain, R. Niedermeier, D. Schultz, Exact Mean Computation in Dynamic
Time Warping Spaces. SIAM International Conference on Data Mining, 2018 (accepted).
[4] R. Caron and T. Traynor. The zero set of a polynomial. WSMR Report, University of Windsor, 2005.
[5] M. Cuturi and M. Blondel. Soft-DTW: a differentiable loss function for time-series. International Conference
on Machine Learning, 2017.
[6] P. Esling and C. Agon. Time-series data mining. ACM Computing Surveys, 45:1–34, 2012.
[7] L.C. Evans and R.F. Gariepy. Measure theory and fine properties of functions. CRC Press, 1992.
[8] T. Fu. A review on time series data mining. Engineering Applications of Artificial Intelligence, 24(1):164–181,
2011.
[9] P.R. Halmos. Measure Theory. Springer-Verlag, 2013.
[10] V. Hautamaki, P. Nykanen, P. Franti. Time-series clustering by approximate prototypes. International
Conference on Pattern Recognition, 2008.
[11] B. Jain. Generalized gradient learning on time series. Machine Learning 100(2-3):587–608, 2015.
[12] B. Jain and D. Schultz. Asymmetric learning vector quantization for efficient nearest neighbor classification
in dynamic time warping spaces. Pattern Recognition, 76:349-366, 2018.
[13] B. Jain. Warped-Linear Models for Time Series Classification arXiv preprint, arXiv:1711.09156, 2017.
[14] T. Kohonen and P. Somervuo. Self-organizing maps of symbol strings. Neurocomputing, 21(1-3):19–30, 1998.
[15] T. Kohonen. Self-Organizing Maps Springer-Verlag Berlin Heidelberg, 2001.
[16] J.B. Kruskal and M. Liberman. The symmetric time-warping problem: From continuous to discrete. Time
warps, string edits and macromolecules: The theory and practice of sequence comparison, 1983.
[17] F. Petitjean, A. Ketterlin, and P. Gancarski. A global averaging method for dynamic time warping, with
applications to clustering. Pattern Recognition 44(3):678–693, 2011.
[18] F. Petitjean, G. Forestier, G.I. Webb, A.E. Nicholson, Y. Chen, and E. Keogh. Faster and more accurate
classification of time series by exploiting a novel dynamic time warping averaging algorithm. Knowledge and
Information Systems, 47(1):1–26, 2016.
[19] R.T. Rockafellar. A Property of Piecewise Smooth Functions. Computational Optimization and Applications,
25:247–250, 2003.
[20] H. Sakoe and S. Chiba. Dynamic programming algorithm optimization for spoken word recognition. IEEE
Transactions on Acoustics, Speech, and Signal Processing, 26(1):43–49, 1978.
[21] D. Schultz and B. Jain. Nonsmooth analysis and subgradient methods for averaging in dynamic time warping
spaces. Pattern Recognition, 74:340–358, 2018.
[22] S. Soheily-Khah, A. Douzal-Chouakria, and E. Gaussier. Progressive and Iterative Approaches for Time
Series Averaging. Workshop on Advanced Analytics and Learning on Temporal Data, 2015.
[23] P. Somervuo and T. Kohonen, Self-organizing maps and learning vector quantization for feature sequences.
Neural Processing Letters, 10(2):151–159, 1999.
[24] Z. Xing, J. Pei, and E. Keogh. A brief survey on sequence classification. ACM SIGKDD Explorations
Newsletter, 12(1):40–48, 2010.

13

