
an important necessary condition for flatness. Even though
for discrete-time systems the property of forward-flatness
can be checked efficiently by a generalization of the test
for static feedback linearizability (see Kolar et al. (2022))
which is based on a certain decomposition property derived
in Kolar et al. (2021), for the more general case including
both forward- and backward-shifts of the system variables
a computationally feasible test does not yet exist. 1Hence,
as we shall illustrate by our second example, the derived
necessary condition is a useful possibility to prove that a
given discrete-time system is not flat.
The paper is organized as follows: First, Section 2 deals
with the concept of discrete-time flatness for time-varying
systems. The core of the paper is then contained in Section
3, which studies the relation between a flat system and
the linear time-varying system obtained by a linearization
along a trajectory. The presented results are illustrated by
two examples in Section 4.
Notation Since we apply differential-geometric concepts,
we use index notation and the Einstein summation conven-
tion to keep formulas short and readable. However, to high-
light the summation range especially for double sums, we
also frequently indicate the summation explicitly. For coor-
dinates that represent forward- or backward-shifts of sys-
tem variables, we use a notation with subscripts in brack-
ets. For instance, the α-th forward- or backward-shift of a
component yj,j∈ {1, . . . , m}of a flat output ywith α∈Z
is denoted by yj
[α], and y[α]= (y1
[α], . . . , ym
[α]). Furthermore,
to facilitate the handling of expressions which depend on
different numbers of shifts of different components of a flat
output, we use multi-indices. If A= (a1, . . . , am) is some
multi-index, then y[A]= (y1
[a1], . . . , ym
[am]).
2. FLATNESS OF TIME-VARYING DISCRETE-TIME
SYSTEMS
In this contribution, we consider nonlinear time-varying
discrete-time systems
xi,+=fi(k, x, u), i = 1, . . . , n (1)
with dim(x) = n, dim(u) = m, and smooth functions
fi(k, x, u). In addition, we assume that the system (1)
meets the submersivity condition
rank(∂(x,u)f) = n , (2)
which is quite common in the discrete-time literature, for
all time-steps k.
As proposed in Diwold et al. (2022b), where only time-
invariant systems are considered, we call a time-varying
discrete-time system (1) flat if there exists a one-to-one
correspondence between its trajectories (x(k), u(k)) and
the trajectories y(k) of a trivial system with dim(y) =
dim(u). The trajectories of a trivial system are not re-
stricted by any difference equation and hence completely
free. By one-to-one correspondence, we mean that the
values of x(k) and u(k) at a time-step kare determined by
an arbitrary but finite number of future and past values of
y(k), i.e., by the trajectory y(k) in an arbitrarily large but
1An interesting approach can be found in Kaldm¨ae (2022) but
requires the solution of partial differential equations.
Fig. 1. One-to-one correspondence of the trajectories.
finite time window. Conversely, the value of y(k) at a time-
step kis determined by an arbitrary but finite number of
future and past values of x(k) and u(k). Consequently,
the one-to-one correspondence of the trajectories can be
expressed by maps of the form
(x(k), u(k)) = F(k, y(k−r1), . . . , y(k+r2)) (3)
and
y(k) = Φ(k, x(k−q1), u(k−q1), . . . , x(k+q2), u(k+q2))
(4)
with suitable integers r1, r2, q1, q2that describe the length
of the corresponding finite time windows, cf. Fig. 1. Since
the number of forward- and backward-shifts in (3) and (4)
can of course be different for the individual components
of y,x, and u, we will later use appropriate multi-indices
where it is important.
In the remainder of this section, the framework used in
Diwold et al. (2022b) for the analysis of flat time-invariant
discrete-time systems is adapted to the time-varying case.
First, it is important to note that the representation of a
trajectory of the system (1) by both sequences x(k) and
u(k) contains redundancy, as these sequences are coupled
by the system equations (1). By a repeated application of
(1), all forward-shifts x(k+α), α≥1 of the state variables
are determined by x(k) and the input trajectory u(k+α)
for α≥0:
x(k+ 1) = f(k, x(k), u(k))
x(k+ 2) = f(k+ 1, x(k+ 1), u(k+ 1))
.
.
.
In the case rank(∂xf) = n, the same is also true for the
backward-direction. However, even if the system meets
only the weaker submersivity condition (2), there exist
mfunctions g(k, x, u) such that the (n+m)×(n+m)
Jacobian matrix ∂xf ∂uf
∂xg ∂ug(5)
is regular for all k. With such functions, the map
x+=f(k, x, u)
ζ=g(k, x, u)(6)
is locally invertible for all k, and by a repeated application
of its inverse
(x, u) = ψ(k, x+, ζ) (7)
all backward-shifts x(k−β), u(k−β), β≥1 of the state-
and input variables are determined by x(k) and backward-
shifts ζ(k−β), β≥1 of the system variables ζdefined by
(6):
(x(k−1), u(k−1)) = ψ(k−1, x(k), ζ(k−1))
(x(k−2), u(k−2)) = ψ(k−2, x(k−1), ζ(k−2))
.
.
.
Hence, every trajectory of the system (1) is uniquely
determined both in forward- and backward-direction by