---
title: Elastically-Scaled Zonotopic Sets
url: https://www.emergentmind.com/topics/elastically-scaled-zonotopic-sets
type: topic
---

# Elastically-Scaled Zonotopic Sets

Elastically-scaled zonotopic sets are zonotopes whose defining segment lengths or generator magnitudes are adjusted in a structured, non-uniform manner to encode geometry or uncertainty. In the geometric setting of symmetric convex random closed sets in \(\mathbb R^2\), the relevant scaling is determined from the Feret diameter process, yielding a zonotopic representation whose face lengths are directionally calibrated by sampled widths [1701.01264]. In the control setting of tube-based Model Predictive Control (MPC), the term denotes a sequence of zonotopes whose generators are scaled componentwise by time-varying nonnegative vectors along the prediction horizon, producing an “elastic tube” representation for robust constrained linear systems [2509.19824].

## 1. Geometric and algebraic foundations

A random closed set (RACS) is a random variable taking values in the family of closed subsets of \(\mathbb R^2\). The geometric framework in [1701.01264] restricts attention to random convex sets, meaning that almost every realization \(X(\omega)\) is compact and convex, and then further to symmetric sets satisfying
\[
X=\breve X:=-X,
\]
or equivalently
\[
X=\frac12(X\oplus \breve X).
\]
This symmetry is structurally decisive because the Feret diameter of a convex set determines only its symmetrized version in general; for symmetric sets it determines the set itself [1701.01264].

A zonotope is a Minkowski sum of line segments,
\[
X=\bigoplus_{i=1}^N \alpha_i S_{\theta_i}, \qquad \alpha_i\ge 0,\ \theta_i\in[0,\pi),
\]
where \(S_0=[-\tfrac12,\tfrac12]\) and \(S_\theta\) is its rotation by angle \(\theta\). In this representation, a zonotope is a compact convex symmetric polygon with at most \(2N\) faces. For regular angular grids
\[
\theta_i=\frac{(i-1)\pi}{N},\qquad i=1,\dots,N,
\]
the paper defines the class \(\mathcal C_0^{(N)}\) of 0-regular zonotopes, its rotated analogue \(\mathcal C_t^{(N)}\), and the union
\[
\mathcal C_\infty^{(N)}=\bigcup_t \mathcal C_t^{(N)},
\]
the class of regular zonotopes with \(2N\) faces [1701.01264].

In the MPC formulation of [2509.19824], the same object is expressed in generator form,
\[
Z=\langle c,G\rangle=\{x\in \mathbb R^n:\: x=c+G\xi,\ \|\xi\|_{\infty}\le 1\},
\]
with center \(c\), generator matrix \(G\), and coefficient vector \(\xi\). A scaled zonotope is obtained by componentwise scaling of the generators,
\[
G\leftarrow G\Delta,\qquad \Delta=\diag(\delta),\qquad \delta\in \mathbb R_{\ge 0}^D,
\]
so that
\[
Z(\delta)=\langle c,G\Delta\rangle.
\]
This algebraic form makes explicit that zonotopic deformation can be encoded entirely by a nonnegative scaling vector.

## 2. Feret-diameter description and directional reconstruction

The Feret diameter, or caliper diameter, of a convex set \(X\) in direction \(\theta\) is the width of \(X\) in that direction. In [1701.01264] it is written through the support function \(h_X\) as
\[
H_X(\theta)=h_X(\theta)+h_{\breve X}(\theta),
\]
and for symmetric \(X\), \(H_X\) fully characterizes \(X\). For a random convex set \(X\), the family \(\{H_X(t),\,t\in\mathbb R\}\) is the Feret diameter random process; it is \(\pi\)-periodic and continuous along each realization.

For a zonotope
\[
X=\bigoplus_{i=1}^N \alpha_i S_{\theta_i},
\]
its Feret diameter satisfies the deterministic identity
\[
H_X(\eta)=\sum_{i=1}^N \alpha_i\lvert \sin(\eta-\theta_i)\rvert.
\tag{1}
\]
This identity is the central bridge between geometry and algebra: sampling \(H_X\) at the angles \(\theta_i\) produces a linear system in the face lengths \(\alpha_i\) [1701.01264].

The paper’s “elastic scaling” interpretation follows from this reconstruction formula. In the approximation
\[
\alpha = F^{(N)-1}H_X^{(N)},
\]
the face lengths are not chosen freely; they are the coefficients required to match the sampled caliper widths exactly. Under a global homothety \(Y=kX\), the Feret diameter scales linearly,
\[
H_{kX}=kH_X,
\]
and therefore
\[
\alpha(kX)=k\,\alpha(X).
\]
This suggests that the elastic aspect, in the geometric sense, is a directional calibration mechanism rather than a uniform rescaling: each face length responds to the directional content of the observed width process [1701.01264].

## 3. Approximation of symmetric convex random closed sets by zonotopes

The main deterministic approximation result in [1701.01264] states that any symmetric convex set \(X\) can be approximated as closely as desired by a 0-regular zonotope constructed from Feret diameter samples. Define
\[
F^{(N)}=\bigl(|\sin(\theta_i-\theta_j)|\bigr)_{1\le i,j\le N},
\]
and
\[
H_X^{(N)} = {}^t\bigl(H_X(\theta_1),\dots,H_X(\theta_N)\bigr).
\]
Then the \(\mathcal C_0^{(N)}\)-approximation is
\[
X_0^{(N)}=\bigoplus_{i=1}^N \bigl(F^{(N)-1}H_X^{(N)}\bigr)_i\, S_{\theta_i}.
\tag{2}
\]

This approximation has three stated properties. First, \(X_0^{(N)}\in\mathcal C_0^{(N)}\). Second, it is the unique element of \(\mathcal C_0^{(N)}\) satisfying
\[
H_{X_0^{(N)}}(\theta_i)=H_X(\theta_i),\qquad i=1,\dots,N.
\]
Third, it converges to \(X\) in Hausdorff distance with the explicit bound
\[
d_H(X,X_0^{(N)})\le (6+2\sqrt2)\sin\!\Bigl(\frac{\pi}{2N}\Bigr)\operatorname{diam}(X),
\tag{3}
\]
hence
\[
d_H(X,X_0^{(N)})\to 0 \quad \text{as } N\to\infty.
\]

In the random setting, the same construction is applied realization-wise, producing a random zonotope approximation with consistency as \(N\to\infty\). The paper also defines a rotationally optimized approximation \(X_\infty^{(N)}\in\mathcal C_\infty^{(N)}\), obtained by selecting the best orientation among all rotations. This removes the orientation sensitivity of the class \(\mathcal C_0^{(N)}\) [1701.01264].

A common misconception is that zonotopic approximation in this framework is merely polygonal curve fitting. The deterministic statement is stronger: the approximation is explicitly tied to the Feret diameter samples, is unique within the chosen regular class, and comes with a Hausdorff-distance estimate. The role of symmetry is equally non-optional; without symmetry, the Feret diameter determines the symmetrized body rather than the original set.

## 4. Tube parameterization and elastic scaling in robust MPC

In [2509.19824], elastic scaling is formalized for robust tube MPC. The scaling vector is allowed to vary along the prediction horizon,
\[
\delta_k\in\mathbb R_{\ge 0}^D,\qquad k=0,\ldots,N,
\]
so the tube cross-section at time \(k\) is a zonotope whose generators are stretched or shrunk over time. This differs from rigid tubes, which use fixed-shape translation only, from homothetic tubes, which use one scalar scaling factor for all generators, and from elastic polyhedral tubes, where half-space offsets vary but inclusion conditions are generally more complex and can become bilinear [2509.19824].

The tube is parameterized as
\[
x_k \in \bar x_k \oplus \langle c_k, G\Delta_k\rangle, \qquad \Delta_k=\diag(\delta_k),
\]
where \(\bar x_k\) is the nominal trajectory, \(c_k\) is the tube center offset, \(G\) is a fixed seed generator matrix, and \(\delta_k\) is the scaling vector. The tuple \((\bar x_k,c_k,\delta_k)\) is admissible if there exists \(u_k\in\mathcal U\) such that one-step reachability and state admissibility hold, namely
\[
x_{k+1}\in \bar x_{k+1}\oplus \langle c_{k+1},G\Delta_{k+1}\rangle,
\]
and
\[
x_k \in \bar x_k \oplus \langle c_k,G\Delta_k\rangle \subseteq \langle c_X,G_X\rangle.
\]
Input admissibility is reduced to
\[
u_k\in \langle c_U,G_U\rangle.
\]

The applied control law is
\[
u_k=\bar u_k + K(x_k-\bar x_k),
\]
with stabilizing feedback gain \(K\) and \(A_K=A+BK\). Under this feedback, one-step reachability is guaranteed if
\[
\left\langle A\bar x_k+B\bar u_k+A_Kc_k,\ A_KG\Delta_k\right\rangle\oplus\langle c_W,G_W\rangle \subseteq \langle \bar x_{k+1}+c_{k+1},G\Delta_{k+1}\rangle.
\]
State and input admissibility are ensured by
\[
\langle \bar x_k+c_k,G\Delta_k\rangle \subseteq \langle c_X,G_X\rangle,
\]
and
\[
\langle \bar u_k+Kc_k,KG\Delta_k\rangle \subseteq \langle c_U,G_U\rangle.
\]
The essential point is that the tube shape is not fixed a priori; it is parameterized by the horizon-indexed vectors \(\delta_k\), which act as elastic degrees of freedom [2509.19824].

## 5. Linear inclusion conditions for scaled zonotopes

A major contribution of [2509.19824] is a linear inclusion test for zonotopes. For zonotopes \(Z_i=\langle c_i,G_i\rangle\subset\mathbb R^d\), a sufficient condition for
\[
Z_1\oplus\cdots\oplus Z_{n-1}\subseteq Z_n
\]
is the existence of matrices \(\Gamma_i\in \mathbb R^{D_n\times D_i}\) and a vector \(\gamma\in\mathbb R^{D_n}\) such that
\[
G_i = G_n\Gamma_i,\qquad \sum_{i=1}^{n-1}c_i = c_n + G_n\gamma,
\]
and
\[
\sum_{i=1}^{n-1} |\Gamma_i|\mathbf 1_{D_i} + |\gamma| \le \mathbf 1_{D_n}.
\]
The paper then rewrites this containment criterion as a linear condition using an auxiliary matrix \(\Phi\):
\[
G_n\Phi^\top V^\top = \begin{bmatrix} -c_n+\sum_{i=1}^{n-1}c_i & G_1 & \ldots & G_{n-1} \end{bmatrix},
\]
\[
\Phi \ge \mathbf 0_{2\bar D\times D_n},\qquad \Phi^\top \mathbf 1_{2\bar D}\le \mathbf 1_{D_n},
\]
with
\[
V=\begin{bmatrix}I_{\bar D} & -I_{\bar D}\end{bmatrix},\qquad \bar D=1+\sum_{i=1}^{n-1}D_i.
\]

For scaled zonotopes
\[
Z_i(\delta_i)=\langle c_i,G_i\diag(\delta_i)\rangle,
\]
the inclusion condition becomes
\[
G_n\Phi^\top V^\top = \begin{bmatrix} -c_n+\sum_{i=1}^{n-1}c_i & G_1\Delta_1 & \ldots & G_{n-1}\Delta_{n-1} \end{bmatrix},
\]
\[
\Phi \ge \mathbf 0_{2\bar D\times D_n},\qquad \Phi^\top \mathbf 1_{2\bar D}\le \delta_n.
\]
The paper emphasizes that this encoding is linear in \(\Phi\) and linear in the scaling factors \(\delta_i\), so it can be embedded directly into optimization problems without bilinearities.

To reduce online complexity, the paper proposes a precomputed matrix \(\Phi_0\) satisfying the unscaled inclusion condition, and then imposes the sufficient condition
\[
\Phi_0^\top \left( \begin{bmatrix}1\end{bmatrix}\otimes \begin{bmatrix} 1\ \delta_1\ \vdots \ \delta_{n-1} \end{bmatrix} \right)\le \delta_n.
\]
The corresponding precomputation problem is
\[
\Phi_0 = \arg\min_\Phi \|\Phi \cdot \mathbf 1_{\bar D}\|_\infty \quad \text{s.t. }\eqref{eq:zon_inclusion_test2}.
\]
This yields the stated trade-off: \(\Phi\) is more flexible and more expensive, whereas \(\Phi_0\) is cheaper but more restrictive [2509.19824].

## 6. Moment formulas, inversion, and optimization structure

In the geometric setting, zonotopic representation is constructive. Once the Feret diameters \(H_X(\theta_i)\) are known at regular angles, the face-length vector is recovered from
\[
F^{(N)}\alpha = H_X^{(N)},
\qquad
\alpha = F^{(N)-1}H_X^{(N)}.
\]
For the 0-regular approximation, the paper gives the direct inversion formulas
\[
\alpha(\omega)=F^{(N)-1}H_X^{(N)}(\omega),
\tag{10}
\]
\[
\mathbb{E}[\alpha]=F^{(N)-1}\mathbb{E}[H_X^{(N)}],
\tag{11}
\]
\[
C[\alpha]=F^{(N)-1}C[H_X^{(N)}]\,{}^t\!F^{(N)-1}.
\tag{12}
\]
For a zonotope
\[
X=\bigoplus_{i=1}^N \alpha_i S_{\theta_i},
\]
its perimeter and area are
\[
U(X)=2\sum_{i=1}^N \alpha_i,
\tag{4}
\]
and
\[
A(X)=\frac12\sum_{i=1}^N\sum_{j=1}^N \alpha_i\alpha_j \,|\sin(\theta_i-\theta_j)|.
\tag{5}
\]

For isotropic random regular zonotopes, [1701.01264] introduces the kernel
\[
k_S(t)=\frac1\pi\int_0^\pi |\sin(t+z)\sin z|\,dz,
\]
with explicit form on \([0,\pi]\),
\[
k_S(t)=\frac{1}{2\pi}\Bigl(2\sin^3(t)+\cos(t)\bigl(\pi-2t+\sin(2t)\bigr)\Bigr),
\tag{6}
\]
and derives the second-order Feret relation
\[
\mathbb{E}\!\left[H_X(t)H_X(t+t')\right] = \sum_{i=1}^N\sum_{j=1}^N \mathbb{E}[\alpha_i\alpha_j]\, k_S(t'+\theta_i-\theta_j).
\tag{7}
\]
For a central face-length vector \(\alpha\) of an isotropic regular random zonotope,
\[
\mathbb{E}[\alpha_i]=\frac{\pi}{2N}\,\mathbb{E}[H_X(x)], \qquad \forall i,
\tag{8}
\]
and
\[
V[\alpha]=\frac{1}{N}K(0)^{-1}V[H_X^{(N)}].
\tag{9}
\]
If noisy data are available, the paper suggests the least-squares estimation
\[
\tilde V[\alpha] =\arg\min_{V\in\mathbb R_+^N}\bigl\|\hat V[H_X^{(2(N'-1))}] - QV\bigr\|^2.
\]

In the MPC setting, the full elastic-tube optimization problem is
\[
\bar{\mathbf u}^\star_{\mathbf N} = \arg\min_{\bar{\mathbf u}_{\mathbf N}, \Phi_{\mathbf N}, \delta_{\mathbf N},\bm c_{\mathbf N}, \bm \lambda_{\mathbf N}} V_T(\bar x_N, \delta_N) +\sum_{k=0}^{N-1}\ell(\bar x_k,\bar u_k,\delta_k),
\]
subject to constraints including
\[
\bar x_0\oplus \langle c_0,G\Delta_0\rangle \subseteq \mathcal X_0,
\]
\[
c_k=G\lambda_k,\qquad |\lambda_k|\le 1,
\]
\[
F_X(c_k+\bar x_k)+|F_XG|\delta_k\le \theta_X,
\]
\[
F_U(Kc_k+\bar u_k)+|F_UKG|\delta_k\le \theta_U,
\]
\[
G\Phi_k^\top V^\top= \begin{bmatrix} -c_{k+1}-\bar x_{k+1}+A\bar x_k+B\bar u_k+A_Kc_k+c_\omega \mid A_KG\Delta_k\mid G_\omega \end{bmatrix},
\]
\[
\Phi_k\ge \mathbf 0_{2\bar D\times D_n},\qquad \Phi_k^\top\mathbf 1_{2\bar D}\le \delta_{k+1},
\]
and terminal admissibility
\[
F_T(c_N+\bar x_N)+|F_TG|\delta_N\le \theta_T.
\]
The stage and terminal costs are
\[
\ell(\bar x,\bar u,\delta) = \bar x^\top Q_x \bar x + \bar u^\top Q_u \bar u + (\delta-\mathbf 1)^\top Q_\delta(\delta-\mathbf 1),
\]
and
\[
V_T(\bar x,\delta) = \bar x^\top P_x \bar x + (\delta-\mathbf 1)^\top P_\delta(\delta-\mathbf 1).
\]
These formulas show that elastic scaling enters not only the feasible-set description but also the optimization objective [2509.19824].

## 7. Computational trade-offs, invariance, and applications

The control-theoretic treatment in [2509.19824] makes the computational trade-offs explicit. The paper compares polyhedral and zonotopic tubes, together with rigid, homothetic, and elastic parameterizations. Its stated conclusions are that polyhedral elastic tubes are expensive because the polyhedral robust positively invariant representation can have
\[
q = 2{D\choose n-1} \gg D
\]
constraints; zonotopic formulations use only \(D\) generators, so the number of scaling variables is much smaller; \(\mathcal P+\mathcal E\) becomes impractical for \(n\ge 3\); \(\mathcal Z+\mathcal E(\Phi)\) has the most flexibility but highest online complexity; \(\mathcal Z+\mathcal E(\Gamma)\) and \(\mathcal Z+\mathcal E(\Phi_0)\) are more compact; and \(\mathcal Z+\mathcal E(\Phi_0)\) removes the need for online auxiliary variables. The resulting three-way trade-off is among size of the domain of attraction, computational burden, and flexibility of the tube shape.

The same paper gives a zonotopic approximation of the minimal robust positively invariant (RPI) set for the error dynamics
\[
z_{k+1}=A_K z_k+\omega_k.
\]
With
\[
\bar G = \begin{bmatrix}I & A_K & \ldots & A_K^s\end{bmatrix}G_\omega,
\]
the associated program is
\[
\max_{c,\Phi,\delta} \quad \delta
\]
subject to
\[
\bar G\Phi^\top V^\top = \begin{bmatrix} -c+A_Kc+c_\omega & A_K\bar G\Delta & G_\omega \end{bmatrix},
\]
\[
\Phi\ge 0,\qquad \Phi^\top\mathbf 1_{2\bar D}\le \delta.
\]
This returns a zonotope \(\langle c,\bar G\Delta\rangle\) that is robust positively invariant. The paper states that this LP uses
\[
2sD_\omega(1+D_\omega+sD_\omega)
\]
decision variables, fewer than the earlier CDC’24 formulation with
\[
3sD_\omega(1+D_\omega+sD_\omega)
\]
variables.

For terminal conditions, the nominal terminal control law is
\[
\bar u_T(\bar x)=K_T\bar x,
\]
with terminal closed-loop dynamics
\[
\bar x_{N+i+1}=A_{K_T}\bar x_{N+i},
\]
and
\[
\delta_{N+i+1} = \Phi_0^\top \left( \mathbf 1_2\otimes \begin{bmatrix} \mathbf 1_{1+D_w}\ \delta_{N+i} \end{bmatrix} \right) = L\delta_{N+i}+d.
\]
Terminal admissibility is encoded by
\[
F_X(\bar x_{N+i}+c_{N+i})+|F_XG|\delta_{N+i}\le \mathbf 1,
\]
and
\[
F_U(K_T\bar x_{N+i}+Kc_{N+i})+|F_UKG|\delta_{N+i}\le \mathbf 1,
\]
while the Lyapunov-like conditions are guaranteed by the LMIs
\[
A_{K_T}^\top P_xA_{K_T}-P_x+Q_x+K_T^\top Q_uK_T\prec 0,
\]
and
\[
L^\top P_\delta L-P_\delta+Q_\delta\prec 0.
\]

The practical implications differ across the two literatures but are conceptually aligned. In stochastic geometry and image-based shape analysis, the geometric construction replaces an arbitrary symmetric convex random set by a finite-dimensional random vector of face lengths linked to observable Feret-diameter data [1701.01264]. In robust MPC, elastic generator scaling yields a tube description that is more scalable than polyhedral elastic tubes and, in the numerical examples reported, maintains feasibility in higher-dimensional settings where polyhedral methods fail [2509.19824]. A plausible implication is that elastically-scaled zonotopic sets are best understood not as a single narrowly defined object, but as a common zonotopic principle: fixed combinatorial structure combined with adaptable directional magnitudes, calibrated either by geometric width data or by horizon-dependent robustness requirements.

Source: https://www.emergentmind.com/topics/elastically-scaled-zonotopic-sets