---
title: Nonlinear Fundamental Lemma in Data-Driven Control
url: https://www.emergentmind.com/topics/nonlinear-fundamental-lemma
type: topic
---

# Nonlinear Fundamental Lemma in Data-Driven Control

Searching arXiv for recent papers on the Nonlinear Fundamental Lemma and closely related formulations.
I’ll look up recent arXiv results to ground the article in the latest literature.
The **Nonlinear Fundamental Lemma** denotes a family of results that seek a nonlinear analogue of Willems’ Fundamental Lemma, namely an exact trajectory characterization for nonlinear systems from measured data. In recent data-driven control literature, the term is used primarily for nonlinear input-state-output systems that become tractable after a lifted representation, especially through Koopman embeddings and related behavioral constructions [2508.07494]. A separate usage occurs in rough-path theory, where the “non-linear sewing lemma” is a theorem about sewing almost flows into genuine flows via control of a composition defect rather than a trajectory-spanning statement [1810.11987]. In the control-theoretic sense, the central issue is that nonlinear trajectories generally do not satisfy superposition in the original coordinates, so the classical Hankel-span argument must be replaced by an exact representation in suitable lifted variables.

## 1. Problem setting and historical placement

For linear systems, Willems’ Fundamental Lemma works because a finite set of measured trajectories spans all others. For nonlinear systems, this spanning property generally fails in the original coordinates. The modern nonlinear literature therefore asks a more structured question: under what lifting, embedding, or representation can nonlinear trajectories again be characterized by a linear-algebraic data relation [2508.07494]?

Several distinct routes have been developed. One route assumes that the nonlinear system admits an exact finite-dimensional Koopman linear embedding, so that the final data-driven representation can be written directly in raw input-output coordinates without explicitly using the lifting functions online [2409.16389]. Another route establishes an LPV Fundamental Lemma and interprets it as relevant to nonlinear systems only when an LPV embedding is available [2103.16171]. A third route shows that known nonlinear extensions can often be rewritten as exact kernelized trajectory representations, equivalent to specific kernel regression problems, for nonlinear classes that become linear after a suitable lifting [2403.05368]. The most general formulation in the provided material is the product-Hilbert-space construction, which treats nonlinear controlled maps
$$
x_{t+1}=F(x_t,u_t)
$$
through an infinite-dimensional generalized Koopman operator acting on a tensor-product observable space and derives from it a nonlinear fundamental lemma in lifted bilinear coordinates [2508.07494].

| Formulation | Exact scope | Key structural requirement |
|---|---|---|
| Product-Hilbert generalized Koopman | General nonlinear maps with inputs in an infinite-dimensional lift | Product Hilbert space, complete bases, injective liftings, rank condition |
| Finite-dimensional Koopman embedding | Exact data-driven representation for a Koopman-compatible class | Exact finite-dimensional Koopman embedding |
| LPV route | Exact LPV behavioral lemma, nonlinear only via embedding | Minimal LPV representation and LPV PE |
| Kernelized route | Exact for classes linear in feature space | Correct finite-dimensional lifting or kernel |

This comparison shows that the phrase does not denote a single universal theorem. It instead denotes a class of exact nonlinear trajectory characterizations whose common strategy is to replace superposition in the original variables by linear or bilinear structure in lifted coordinates.

## 2. Product Hilbert spaces and the generalized Koopman operator

The product-Hilbert-space formulation starts from the discrete-time controlled nonlinear system
$$
x_{t+1}=F(x_t,u_t), \qquad t\in \mathbb{N},
$$
with state space \(X\subseteq \mathbb{R}^n\), input space \(U\subseteq \mathbb{R}^m\), and \(F:X\times U\to X\) assumed Lebesgue integrable. For the input-state-output setting one augments this with
$$
\begin{aligned}
x_{t+1} &= F(x_t,u_t),\\
y_t &= h(x_t),
\end{aligned}
$$
where \(h:X\to Y\), \(Y\subseteq \mathbb{R}^p\), is also Lebesgue integrable [2508.07494].

The basic construction introduces the observable Hilbert spaces
$$
H_x=L^2(X,\mu_x), \qquad H_u=L^2(U,\mu_u), \qquad H_y=L^2(Y,\mu_y),
$$
and then the product Hilbert space
$$
H:=H_x\otimes H_u \cong L^2(X\times U,\mu_x\otimes \mu_u).
$$
State and input observables are chosen as countable families
$$
\{\psi_{i,x}\}_{i=1}^\infty \subset H_x, \qquad \{\psi_{j,u}\}_{j=1}^\infty \subset H_u,
$$
with lifted coordinates
$$
\Psi_x(x)=\operatorname{col}(\psi_{1,x}(x),\psi_{2,x}(x),\ldots), \qquad
\Psi_u(u)=\operatorname{col}(\psi_{1,u}(u),\psi_{2,u}(u),\ldots).
$$
The product-space observable is
$$
\Psi(x,u):=\Psi_x(x)\otimes \Psi_u(u),
$$
whose basis functions are all products \(\psi_{i,x}\psi_{j,u}\). If the state and input families are orthonormal bases of \(H_x\) and \(H_u\), then their tensor products form an orthonormal basis of \(H_x\otimes H_u\); more generally, tensor products of Riesz bases again form a Riesz basis [2508.07494].

This leads to the generalized Koopman operator
$$
K:H_x\otimes H_u \to H_x
$$
defined by
$$
\Psi_x(F(x,u)) = K\big(\Psi_x(x)\otimes \Psi_u(u)\big).
$$
With \(z_t:=\Psi_x(x_t)\) and \(v_t:=\Psi_u(u_t)\), the lifted dynamics become
$$
z_{t+1}=K(z_t\otimes v_t).
$$
The operator is linear on observables, but the induced realization is bilinear in the lifted state and lifted input. This differs from stacked-observable constructions, where state and input lifts are concatenated additively and often require input-affine structure. Here the tensor product makes the state-input interaction multiplicative at the level of coordinates, and no input-affine assumption on \(F\) is needed [2508.07494].

An existence theorem is stated first for orthonormal bases and then for Riesz bases. In one key form, if \(\{\psi_{i,x}\}\) and \(\{\psi_{i,u}\}\) are independent and complete sets of Riesz basis functions spanning \(H_x\) and \(H_u\), and if
$$
\psi_{i,x}\circ F \in H_x\otimes H_u,\qquad i\in\mathbb{N},
$$
then there exists a state transition matrix \(K=Q_xR^{-1}\) such that
$$
\Psi_x(F(x,u)) = K(\Psi_x(x)\otimes \Psi_u(u)).
$$
For the output map, the same logic yields a linear operator \(C:H_x\to H_y\) satisfying
$$
\Psi_y(h(x))=C\Psi_x(x),
$$
hence the lifted input-state-output system
$$
\begin{aligned}
z_{t+1}&=K(z_t\otimes v_t),\\
w_t&=\Psi_y(y_t)=Cz_t.
\end{aligned}
$$

## 3. The nonlinear fundamental lemma in lifted bilinear coordinates

The product-space construction becomes a nonlinear fundamental lemma only after an exact equivalence between original and lifted trajectories has been established. Under independent and complete Riesz bases for \(H_x\), \(H_u\), and \(H_y\), under the composition assumptions
$$
\psi_{i,x}\circ F \in H_x\otimes H_u,\qquad
\psi_{i,y}\circ h\in H_x,
$$
and under injectivity of \(\Psi_x\) and \(\Psi_y\), trajectories of the nonlinear system are in one-to-one correspondence with trajectories of the lifted bilinear system
$$
z_{t+1}=K(z_t\otimes v_t),\qquad w_t=Cz_t
$$
via
$$
z_t=\Psi_x(x_t),\qquad v_t=\Psi_u(u_t),\qquad w_t=\Psi_y(y_t)
$$
[2508.07494].

The data matrices are then built from measured lifted trajectories. For \(i\le j\),
$$
(z\otimes v)_{[i,j]} := \operatorname{col}(z_i\otimes v_i, z_{i+1}\otimes v_{i+1},\ldots,z_j\otimes v_j),
$$
and one defines the Hankel-like matrix
$$
F_N(T):=
\begin{bmatrix}
(z\otimes v)_{[0,N-1]} &
(z\otimes v)_{[1,N]} &
\cdots &
(z\otimes v)_{[T,T+N-1]}
\end{bmatrix},
$$
together with the augmented matrix
$$
F_N^w(T):=
\begin{bmatrix}
F_N(T)\;
w_{[1,N]}\;
w_{[2,N+1]}\;
\cdots\;
w_{[T+1,T+N]}
\end{bmatrix}.
$$
The paper also observes that
$$
F_N(T)=H_{N,T}(z)\odot H_{N,T}(v)
=H_{N,T}(\Psi_x(x))\odot H_{N,T}(\Psi_u(u)),
$$
where \(\odot\) is the blockwise Khatri–Rao product [2508.07494].

The resulting theorem states that, under the hypotheses above and assuming that the infinite-dimensional rectangular matrix \(F_N(T)\) satisfies the Rouché–Capelli rank condition together with the assumptions of the cited infinite-dimensional linear-algebra result, a length-\(N\) nonlinear trajectory segment
$$
\{x_{[t,t+N]},u_{[t,t+N-1]},y_{[t+1,t+N]}\}
$$
is a solution of the nonlinear system if and only if there exists a compatible real vector \(g\) such that
$$
F_N^w(T)g=
\begin{bmatrix}
(z\otimes v)_{[t,t+N-1]}\\
w_{[t+1,t+N]}
\end{bmatrix}.
$$
This is the exact nonlinear analogue of Willems’ lemma in the sense of [2508.07494]: every valid nonlinear trajectory can be represented as a linear combination of recorded lifted trajectory segments, but the representation holds in \((z\otimes v,w)\), not directly in \((x,u,y)\).

A common misconception is that such a result is merely the classical linear lemma in disguise. The cited formulation explicitly rejects that interpretation. The original map \(F(x,u)\) may be an arbitrary nonlinear map, without input-affine structure; what changes is the coordinate architecture. The linear trajectory space of the classical theory is replaced by a lifted bilinear trajectory set generated in a product Hilbert space [2508.07494].

## 4. Finite-dimensional approximation, prediction, and observables

Because the exact theory is infinite-dimensional, computation proceeds through finite-dimensional approximation. The proposed estimator is an EDMD-type regression adapted to the tensor-product lift. Given snapshots \(\{(x_t,u_t,x_{t+1})\}_{t=0}^{T-1}\), define
$$
Z_X=\begin{bmatrix}z_0&z_1&\cdots&z_{T-1}\end{bmatrix},\qquad
V_U=\begin{bmatrix}v_0&v_1&\cdots&v_{T-1}\end{bmatrix},
$$
and
$$
Z_X^+=\begin{bmatrix}z_1&z_2&\cdots&z_T\end{bmatrix}.
$$
Instead of fitting
$$
Z_X^+\approx AZ_X+BV_U,
$$
the method fits
$$
Z_X^+=\hat K\,(Z_X\odot V_U),
$$
with
$$
Z_X\odot V_U=
\begin{bmatrix}
z_0\otimes v_0 & z_1\otimes v_1 & \cdots & z_{T-1}\otimes v_{T-1}
\end{bmatrix},
$$
and, assuming full row rank,
$$
\hat K = Z_X^+(Z_X\odot V_U)^\dagger.
$$
The same lifted data also support a least-squares behavioral representation
$$
g^\ast = F_N(T)^\dagger (z\otimes v)_{[t,t+N-1]},
$$
and a multi-step output predictor
$$
w^\ast_{[t+1,t+N]}
=
\begin{bmatrix}
w_{[1,N]}&\cdots&w_{[T+1,T+N]}
\end{bmatrix}
g^\ast.
$$
Equivalently,
$$
w^\ast_{[t+1,t+N]}
=
\begin{bmatrix}
w_{[1,N]}&\cdots&w_{[T+1,T+N]}
\end{bmatrix}
F_N(T)^\dagger (z\otimes v)_{[t,t+N-1]}.
$$
The authors interpret this as a multi-step generalization of EDMD; if \(N=1\) and \(y_t=x_t\), it reduces to the one-step Koopman regression [2508.07494].

The observable design remains central. The paper discusses universal kernels, especially inverse multiquadrics,
$$
\psi_{x,i}(x)=\left(1+\frac{\|\xi_{i,x}-x\|_2^2}{\sigma_x^2}\right)^{-\beta},\qquad
\psi_{u,i}(u)=\left(1+\frac{\|\xi_{i,u}-u\|_2^2}{\sigma_u^2}\right)^{-\beta},
$$
deep neural networks with separate state and input encoders combined through a Kronecker product of last hidden layers, and Takens/HAVOK delay embeddings
$$
\Psi_x(x_t)=
\begin{bmatrix}
x_t\\
x_{t+1}\\
\vdots\\
x_{t+n_z-1}
\end{bmatrix},
\qquad
\Psi_u(u_t)=
\begin{bmatrix}
u_t\\
u_{t+1}\\
\vdots\\
u_{t+n_v-1}
\end{bmatrix}.
$$
In all cases, the tensor-product architecture is essential [2508.07494].

The controlled Van der Pol oscillator serves as the main illustration:
$$
\dot x_{1,t}=x_{2,t},\qquad
\dot x_{2,t}=\mu(1-x_{1,t}^2)x_{2,t}-x_{1,t}+u_t,
$$
with \(\mu=1.2\). The comparison is between a benchmark Koopman-with-inputs method based on observables on the stacked variable \(\operatorname{col}(x,u)\) and the proposed generalized Koopman method with separate state and input kernels and their tensor product. The reported outcome is that, as the number of observables increases, the proposed method’s prediction error consistently improves, whereas the benchmark stacked-observable method does not improve despite a much larger feature dimension [2508.07494].

## 5. Related formulations and terminological variants

A second exact control-theoretic formulation appears in the finite-dimensional Koopman-embedding setting. There, the nonlinear system
$$
x_{k+1}=f(x_k,u_k),\qquad y_k=g(x_k,u_k)
$$
is assumed to admit a Koopman linear embedding
$$
z_{k+1}=Az_k+Bu_k,\qquad y_k=Cz_k+Du_k,
$$
with \(z_k=\Phi(x_k)\). Under lifted excitation of order \(L\), and with \(T_{\mathrm{ini}}\ge n_z\), every valid length-\(L=T_{\mathrm{ini}}+N\) nonlinear trajectory is characterized exactly by
$$
\operatorname{col}(U_P,Y_P,U_F,Y_F)g
=
\operatorname{col}(u_{\mathrm{ini}},y_{\mathrm{ini}},u_F,y_F).
$$
The distinctive point is that the final representation uses only measured input-output data and bypasses explicit lifting functions in the online equation, although the proof depends on the existence of \(\Phi\) [2409.16389].

The LPV literature provides a different, conditional route. The main theorem in that line is an LPV Fundamental Lemma for discrete-time LPV systems in a behavioral framework, with scheduling-dependent coefficients in \(\mathcal R(\mathbb P)\). Its relevance to nonlinear systems is indirect: it applies only to the extent that the nonlinear dynamics admit an LPV representation or embedding. The result is therefore a pathway toward nonlinear data-driven analysis, not a direct theorem for arbitrary nonlinear systems [2103.16171].

The kernel-regression viewpoint gives yet another reformulation. In that perspective, meaningful nonlinear extensions are exact kernelized trajectory representations for nonlinear classes that become linear after a suitable lifting. The implicit kernel equation is equivalent to a specific kernel regression problem, and the nonlinear analogue of persistency of excitation becomes a rank condition on the kernel Gramian. Exact results are described for classes such as Hammerstein systems and SISO differentially flat systems in the cited formulations [2403.05368].

The term also has a distinct meaning outside data-driven control. Brault and Lejay’s “non-linear sewing lemma” concerns a two-parameter family of maps \(\phi_{t,s}:V\to V\) with controlled composition defect
$$
\phi_{t,s}\circ\phi_{s,r}\approx \phi_{t,r},
$$
and proves that such an almost flow can be sewn into a genuine flow, measurably under weak assumptions and uniquely as a Lipschitz flow under stronger assumptions [1810.11987]. This is a nonlinear analogue of additive and multiplicative sewing lemmas, not a nonlinear trajectory-spanning theorem.

## 6. Assumptions, limitations, and open directions

The exact product-Hilbert-space theory is broad in the sense that it does not require input-affine structure and is formulated for general square-integrable nonlinear maps, but it is exact only under nontrivial assumptions. These include Lebesgue integrability of \(F\) and \(h\), independent and complete Riesz bases for the observable spaces, closure of compositions in the relevant Hilbert spaces,
$$
\psi_{i,x}\circ F\in H_x\otimes H_u,\qquad
\psi_{i,y}\circ h\in H_x,
$$
and injectivity of \(\Psi_x\) and \(\Psi_y\) when one wants exact correspondence back to the original state and output trajectories [2508.07494].

The nonlinear fundamental lemma itself requires more. In addition to the embedding assumptions, the infinite-dimensional data matrix \(F_N(T)\) must satisfy the Rouché–Capelli rank condition and the summability and convergence hypotheses imported from the cited infinite-dimensional bilinear-systems result. The paper explicitly identifies this as the most challenging assumption. It also notes that even if the lifted Hankel matrices \(H_{N,T}(z)\) and \(H_{N,T}(v)\) have full row rank, their blockwise Khatri–Rao product
$$
F_N(T)=H_{N,T}(z)\odot H_{N,T}(v)
$$
need not. The excitation problem is therefore more subtle than in the linear case, and the conditions under which such infinite-dimensional data matrices arise from trajectories “require further investigation” [2508.07494].

The finite-dimensional Koopman-embedding route has a different limitation profile. Its theorem is exact only for nonlinear systems that admit an exact finite-dimensional Koopman embedding, and the embedding dimension \(n_z\) is generally unknown from finite data. The associated data requirements can become large because width must satisfy \(l\ge mL+n_z\), while depth must satisfy \(T_{\mathrm{ini}}\ge n_z\) unless a smaller observability index is known [2409.16389].

The LPV route is similarly conditional. It depends on existence of an LPV embedding, measurability of the scheduling variable, minimality of the LPV realization, and an LPV notion of persistency of excitation that is substantially more demanding than the LTI version because it depends on the pair \((u,p)\), the representation order, and the scheduling-dependency class [2103.16171]. The kernelized route is exact only when the chosen kernel matches the true finite-dimensional lifting. The cited analysis explicitly notes that for most popular kernels the RKHS is infinite-dimensional, and then exact rank conditions of the form used in the theorem cannot be satisfied with finite data [2403.05368].

Taken together, these results suggest a precise synthesis. The nonlinear fundamental lemma is not a direct restoration of linear superposition in the original coordinates. It is an exact trajectory characterization obtained only after the nonlinear dynamics are recast through an embedding—bilinear in a product Hilbert space, linear in a finite Koopman lift, linear in an LPV representation, or linear in a finite-dimensional feature space. In that sense, the contemporary theory replaces the classical trajectory span by an exact lifted behavioral representation.

Source: https://www.emergentmind.com/topics/nonlinear-fundamental-lemma