---
title: Finite-Dimensional Koopman Representation
url: https://www.emergentmind.com/topics/finite-dimensional-koopman-representation
type: topic
---

# Finite-Dimensional Koopman Representation

Searching arXiv for recent and foundational papers on finite-dimensional Koopman representations.
Searching arXiv for exact finite-dimensional Koopman embeddings and controlled/input-driven extensions.
Searching arXiv for delay-coordinate, RKHS, and learning-based finite-dimensional Koopman approximations.
A finite-dimensional Koopman representation is a lifted model of nonlinear dynamics in which a finite vector of observables evolves under a linear—or, for input-driven variants, structured linear parameter-varying or bilinear—state transition. For a discrete-time system \(x^+=F(x)\), the Koopman operator acts on observables \(g\) by \(\mathcal K g=g\circ F\); for a continuous-time system \(\dot x=f(x)\), the Koopman semigroup is \(\mathcal K^t g=g\circ x(t;\cdot)\) with generator \(\mathcal L g=\langle \nabla g,f\rangle\). The finite-dimensional problem is to find a lifting \(z=\Phi(x)\) and matrices such that \(z^+=Az\) or \(\dot z=Az\), often together with a reconstruction \(x=Cz\). Exact realizations require a finite Koopman-invariant subspace; otherwise one obtains projected, truncated, or learned approximations of an intrinsically infinite-dimensional operator [2607.01819].

## 1. Operator-theoretic formulation

The classical finite-dimensional setting starts with a family of observables \(\{\varphi_i\}_{i=1}^n\) spanning a subspace \(\mathbb V\). If \(\mathbb V\) is Koopman-invariant, \(\mathcal K\mathbb V\subseteq\mathbb V\), then the restriction of the Koopman operator to \(\mathbb V\) is represented by a finite matrix. Writing \(\varphi=[\varphi_1,\dots,\varphi_n]^\top\), one has \(\mathcal K\varphi=K\varphi\) in discrete time or \(\dot \Phi(x)=A\Phi(x)\) in continuous time. When the state itself is included among the observables, the original nonlinear dynamics can be recovered by a linear projection \(x=Cz\), so the lifted linear model is simultaneously predictive and state-reconstructive [2607.01819].

Koopman eigenfunctions furnish the canonical coordinates for such representations. If \(U\phi_j=\lambda_j\phi_j\), then a diagonal representation arises immediately; if generalized eigenfunctions are required, Jordan blocks appear instead. In the diagonalizable case, any finite-dimensional linear Koopman representation is conjugate to one generated by eigenfunctions. This connects the existence of a finite-dimensional representation directly to the spectral structure of the composition operator rather than to local tangent linearization [2010.05377].

The geometry of the state space is encoded by level sets of eigenfunctions. If \(\phi\) is an eigenfunction, its level sets are mapped into one another by the flow, and joint level sets of several eigenfunctions define invariant partitions. A representation is faithful precisely when the joint level sets of the chosen observables separate points. This geometric viewpoint clarifies why finite-dimensional Koopman coordinates are more than a feature map: they organize invariant sets, stable and unstable structures, and equivalence classes under the dynamics [2010.05377].

A broader usage of the term includes finite-dimensional nonlinear autonomous closure maps of the form \(U(\mathbf f)=\Phi(\mathbf f)\), sometimes described as a representation eigenproblem. The linear invariant-subspace picture remains the central case in Koopman-based prediction and control, but this broader formulation is useful for learned latent dynamics in which closure is finite-dimensional without being strictly linear [2010.05377].

## 2. Existence, obstructions, and structural conditions

Exact finite-dimensional Koopman representations are exceptional rather than generic. A state-inclusive invariant subspace is especially restrictive. One precise obstruction is topological: systems with multiple fixed points or more complicated attractors are not globally topologically conjugate to a finite-dimensional linear system, and therefore cannot admit a finite-dimensional linear Koopman subspace that includes the state globally. By contrast, systems with a single isolated fixed point may admit state-inclusive closures, particularly in special polynomial classes [1510.03007].

The Koopman spectrum imposes a second obstruction. Point spectrum supports finite-dimensional linear closure, whereas continuous spectral components obstruct conjugacy to linear finite-dimensional dynamics. In the terminology of representation theory for Koopman systems, a finite-dimensional representation is not conjugate to a linear one when the chosen observables have nontrivial overlap with the continuous-spectrum component. This identifies a precise spectral boundary between exact finite linear closure and intrinsically non-finite or nonlinearly representable behavior [2010.05377].

For controlled nonlinear systems, recent work gives a sharp characterization of exact finite-dimensional Koopman linear embeddings. For the discrete-time system
\[
x^+=f(x,u),
\]
an exact embedding requires a lifting \(\Psi(x)\) and constant matrices \(A_K,B_K,C_K\) such that
\[
\Psi(f(x,u))=A_K\Psi(x)+B_Ku,\qquad x=C_K\Psi(x).
\]
This exists if and only if the system can be transformed into a control-affine preserved (CAP) structure and the autonomous nonlinear subsystem admits a finite-dimensional Koopman closure that contains the state and all nonlinear terms needed for reconstruction. In this sense, the source of approximation error can be separated into an intrinsic structural limitation—absence of CAP—and a lifting limitation—failure to choose a sufficiently rich invariant subspace when CAP does exist [2602.14537].

Even when the autonomous part admits a finite-dimensional lift, systems with inputs generally do not yield an exact lifted LTI model. For continuous-time systems the lifted input term is typically state-dependent; for discrete-time systems it is generally state- and input-dependent. The resulting exact finite-dimensional form is therefore an LPV Koopman model,
\[
\dot z=Az+B_z(p_t)u_t
\quad\text{or}\quad
z_{k+1}=Az_k+B_z(p_k)u_k,
\]
rather than \(z^+=Az+Bu\) with constant \(B\). Constant-input-matrix LTI Koopman models remain approximations whose error depends on the deviation of the true input matrix from the chosen constant approximation [2207.12132].

## 3. Constructive exact embeddings for special system classes

Although general existence is rare, several structured nonlinear classes admit constructive exact finite-dimensional embeddings. A basic example is the polynomial slow-manifold family studied through state-inclusive observables. For systems such as
\[
\dot x_1=\mu x_1,\qquad \dot x_2=\lambda(x_2-P(x_1)),
\]
with polynomial \(P\), the observables formed by the state and the monomials appearing in \(P\) close exactly under differentiation. This produces a finite-dimensional linear lifted system whose eigenfunctions recover intrinsic coordinates on the invariant manifold [1510.03007].

A more systematic constructive result is available for autonomous continuous-time lower-triangular polynomial systems,
\[
\dot x_1=a_1x_1,\quad
\dot x_2=a_2x_2+f_2(x_1),\quad
\dots,\quad
\dot x_n=a_nx_n+f_n(x_1,\dots,x_{n-1}),
\]
where each \(f_i\) is polynomial. For this class there exists an exact finite-dimensional lifting \(\Phi\) containing the state such that
\[
\dot z=Az,\qquad z=\Phi(x),\qquad x=Cz.
\]
The construction is recursive: one augments the state with all monomials generated under differentiation until closure is reached, and finite termination follows from the lower-triangular structure and finite polynomial degree. In the worked 4-state example, the resulting lifted state has dimension \(19\), and the trajectories of the nonlinear and lifted linear systems overlap exactly up to numerical error on the order of \(10^{-15}\) [2301.06557].

An even broader exact theory covers block-oriented polynomial systems assembled from series and parallel interconnections of LTI blocks and static polynomial nonlinearities. For these systems one can build a finite lifted state from Kronecker powers and obtain an exact polynomial input time-invariant Koopman form,
\[
\dot z = A z + L(z)R(u)u,
\qquad
y = C z + L(z)R(u)u,
\]
with constant \(A\) and \(C\), \(L(z)\) linear in \(z\), and \(R(u)\) polynomial in \(u\). If the LTI blocks have no feedthrough, this simplifies to a bilinear time-invariant model. In the reported examples, duplicate monomials can be removed by projection, reducing a 17-dimensional raw lift to 12 dimensions in one case and a 931-dimensional raw lift to 103 dimensions in another, with output errors around \(10^{-13}\), i.e. machine precision [2507.15093].

These constructive results identify a common mechanism behind exact finite-dimensionality: closure is achieved when the nonlinearities belong to a finitely generated algebra that remains closed under the Koopman generator. Triangularity, polynomial degree bounds, and block-oriented acyclic interconnection are specific structural realizations of that mechanism.

## 4. Approximate finite-dimensional models from data and learning

Outside such special classes, finite-dimensional Koopman representations are typically approximate. The standard approximation is EDMD. Given snapshot pairs \(\{(x^{(i)},y^{(i)})\}_{i=1}^N\) with \(y^{(i)}=F(x^{(i)})\) and dictionary \(\psi(x)=[\psi_1(x),\dots,\psi_n(x)]^\top\), EDMD computes
\[
K_{\Psi}=\Psi_Y\Psi_X^\top(\Psi_X\Psi_X^\top)^\dagger
\]
as a least-squares approximation of the projected Koopman operator on the dictionary span. In continuous time, generator EDMD uses
\[
L_\Psi=\dot\Psi_X\Psi_X^\top(\Psi_X\Psi_X^\top)^\dagger.
\]
The accuracy of such a model depends on both finite-sample estimation and the non-invariance of the chosen dictionary. This latter effect is quantified by the invariance proximity
\[
\mathcal I_{\mathcal K(\mathbb V)}
=
\sup_{g\in\mathbb V,\ \|\mathcal Kg\|\neq 0}
\frac{\|\mathcal K g-\mathcal K_{\mathrm{approx}}g\|}{\|\mathcal K g\|},
\]
which vanishes only for an exactly invariant subspace [2607.01819].

Kernel and RKHS formulations replace a fixed dictionary by the finite span of kernel sections centered at the data. In this setting, DMD and EDMD become dual value-based and function-based approximations, and the learned finite-dimensional subspace is shaped by the sample locations. The regularized kernel Koopman estimator can be written entirely in terms of Gram matrices, for example
\[
U_f^{(r,K)}=
\bigl(K(\bar x,\bar x)+\sigma^2I_M\bigr)^{-1}K(\bar y,\bar x),
\]
so the infinite-dimensional operator is approximated on a data-adaptive finite span without prescribing explicit basis functions a priori [2103.13752].

A closely related but conceptually sharper result is the representer theorem for Koopman learning. In the discrete-time autonomous setting, regularized empirical risk minimization over the infinite-dimensional operator space \(\mathcal L(\mathcal H,\mathcal H)\) admits an exact finite-rank optimizer of the form
\[
\hat K=\sum_{k=1}^N\sum_{l=1}^n a_{kl}\,v_k\otimes g_l,
\]
under broad assumptions on the observables and regularizer. This yields a finite-dimensional optimization problem with no approximation or loss of precision at the optimization level, even though the learned operator still models nonlinear dynamics through finite data and finite observables [2208.01681].

Recent machine-learning approaches focus on learning nearly invariant subspaces rather than fixing them. One line learns locally supported, approximately orthogonal basis functions \(\{\phi_j\}_{j=1}^N\) and a latent linear map \(\mathcal G\), producing
\[
A=\mathcal D\circ\mathcal G\circ\mathcal P
\]
as a finite-dimensional approximation of the transfer or Koopman operator. Reconstruction of both \(g\) and \(\mathcal L(g)\), together with an \(L^1\)-based local-support penalty, encourages the learned span to be nearly invariant and spectrally informative [2505.05085]. Another line, Koopmanizing Flows, assumes a nonlinear system is diffeomorphically related to a latent stable linear system \(\dot y=\underline A y\), lifts \(y\) through a truncated monomial basis of order \(\bar p\), and learns the diffeomorphism, linear dynamics, and decoder jointly; the finite feature dimension is
\[
D=\binom{d+\bar p}{d}-1.
\]
The use of a Hurwitz parameterization guarantees asymptotic stability regardless of operator approximation accuracy, and on the LASA handwriting benchmark the method outperforms SKEL on DTWD, RMSE, and PCM [2112.04085].

Approximate closure can also be engineered through specialized dictionaries. State Inclusive Logistic Lifting augments the state with multivariate conjunctive logistic functions,
\[
\psi(x)=\begin{bmatrix}1 & x & \Lambda\end{bmatrix},
\]
to obtain finite approximate closure under differentiation. Its error is governed by the steepness parameter \(\alpha\) and the mesh parameter \(\epsilon\), with improved closure as the logistic functions become sharper and the dictionary resolves the vector field more densely [1712.03132]. Information-theoretic latent models add a different perspective: they treat Koopman learning as an information bottleneck problem and use latent mutual information together with von Neumann entropy to avoid collapse onto a few dominant modes, thereby preserving a finite-dimensional latent Koopman subspace that remains predictive over long horizons [2510.13025].

## 5. Delay coordinates, infinite-dimensional systems, and related operator variants

Time-delay observables provide a distinct route to finite-dimensional Koopman coordinates. Instead of selecting instantaneous observables, one embeds trajectory segments and projects them onto an orthonormal basis on a delay window. The resulting convolutional coordinates
\[
w_j(x)=\langle \phi_j,\mathcal K_{(\cdot)}g(x)\rangle
\]
obey a universal generator representation
\[
\mathcal K w_j(x)=\sum_{k=0}^{\infty}K_{jk}w_k(x),
\qquad
K_{jk}=\langle\phi_j,\phi_k'\rangle,
\]
where the matrix depends only on the chosen delay basis, not on the underlying dynamics. When the measured observable has a finite Koopman mode decomposition of order \(r\), the first \(r\) SVD delay coordinates span an invariant subspace and form an optimal finite-dimensional basis; the induced linear model coincides with DMD computed in those coordinates [1810.01479].

For nonlinear infinite-dimensional systems such as PDEs, the Koopman semigroup can be defined on bounded continuous functionals, then projected onto a finite-dimensional space of functionals \(\mathcal E_n\). With basis \(\{\zeta_i\}_{i=1}^n\), one obtains a compressed semigroup \(K_n^t\) and Lie generator \(L_n\), and the data-driven approximation remains EDMD-like:
\[
\mathbf K=\mathbf\Theta_1^\dagger\mathbf\Theta_2.
\]
The associated matrix logarithm
\[
\widetilde{\mathbf L}^{(t_s)}
=
\frac{1}{t_s}\log(\mathbf\Theta_1^\dagger\mathbf\Theta_2)
\]
can be used both for spectral approximation and for identification of PDE coefficients through the first column of the generator estimate [2103.12458].

A further extension replaces observables or densities by wavefunctions. In the Koopman–von Neumann framework, the generator is
\[
\mathcal Q=\tfrac12(\mathcal L^*-\mathcal L),
\]
and Galerkin projection onto a finite basis yields
\[
Q=\tfrac12 G^{-1}(A^\top-A).
\]
After whitening the basis, the finite-dimensional generator becomes skew-symmetric, so its exponential is unitary. This produces finite matrices suitable for spectral computation and, in principle, quantum-circuit realization, even for non-Hamiltonian dynamics [2604.08414].

The phrase “Koopman representation” also appears in operator-algebraic settings unrelated to state lifting. For self-similar groupoid actions on \(L^2(E^\infty,\nu)\), the Koopman representation is generally infinite-dimensional, but it preserves a filtration by finite-dimensional cylinder-function spaces
\[
H_n=\operatorname{span}\{\chi_{Z(a)}:|a|=n\},
\qquad
\dim(H_n)=|E^0|p^n.
\]
The resulting \(C^*\)-algebra is residually finite-dimensional, meaning that finite-dimensional representations separate points even though the original Koopman representation is not itself finite-dimensional [2112.14341].

## 6. Exactness, approximation error, and unresolved issues

A recurrent misconception is that finite-dimensional Koopman modeling is principally a matter of enlarging the dictionary. The literature instead shows a structural dichotomy. If an exact invariant subspace exists, enlarging the observable set can help recover it. If the system lacks the requisite structure, no observable design can remove the intrinsic approximation error. In controlled systems this distinction is explicit: absence of CAP rules out exact finite-dimensional Koopman linear embeddings, whereas existence of CAP shifts the difficulty to finding a sufficiently rich lifting [2602.14537].

A second misconception is that a richer dictionary is automatically superior. In projection-based methods, larger spaces can reduce representation error while worsening non-invariance error and long-horizon drift. The core tradeoff is not only expressivity but also closure. This is why several recent methods optimize invariance explicitly, either by learning dynamically adapted bases, enforcing latent linearity with stable parameterizations, or regularizing against mode collapse in the latent covariance spectrum [2607.01819].

The exact-versus-approximate distinction is especially important for systems with inputs. Much of the applied literature assumes a lifted LTI model \(z^+=Az+Bu\), but the exact theory shows that even when the autonomous part is finitely representable, the input channel typically produces LPV or bilinear structure rather than a constant \(B\). A plausible implication is that model structure selection—LTI versus LPV versus bilinear—is not a secondary identification choice but a consequence of the nonlinear input mechanism [2207.12132].

Finally, the literature contains claims of very strong general exactness that require careful scrutiny. The abstract of “Functional Dimensionality of Koopman Eigenfunction Space” states that the general form solution of the Koopman Partial Differential Equation is obtained, that the functional dimensionality of the Koopman eigenfunction space is finite and equal to the dimensionality of the dynamics, and that nonlinear dynamics therefore admit an error-free finite set of Koopman eigenfunctions together with a simple numerical method [2401.02272]. However, the supplied document text contains only template structure—title, metadata, theorem environments, and empty section headings—and does not include definitions, propositions, equations, derivations, flowbox statements, or algorithms. From the available text, the claimed theorem, its assumptions, and its numerical consequences cannot be assessed [2401.02272].

In current usage, finite-dimensional Koopman representation therefore denotes a family of closely related but technically distinct constructions: exact invariant-subspace realizations for special structural classes, exact optimization-level reductions for certain learning problems, and approximate projected or learned latent models for generic nonlinear systems. The central question is not whether the Koopman operator is linear—it always is—but whether the nonlinear dynamics admit a finite observable algebra that closes under that linear action, and if not, how the residual non-closure should be quantified and controlled.

Source: https://www.emergentmind.com/topics/finite-dimensional-koopman-representation