---
title: Koopman Control Factorization
url: https://www.emergentmind.com/topics/koopman-control-factorization
type: topic
---

# Koopman Control Factorization

Koopman control factorization is the separation of autonomous dynamics and control channels within Koopman-operator descriptions of nonlinear systems with inputs. For control-affine continuous-time systems it appears at the generator level as \(L(u)=L_0+\sum_{j=1}^m u_j L_j\), or equivalently \(\frac{d}{dt}\mathcal U(t)=\big(\mathcal L_f+\sum_{i=1}^m u_i(t)\mathcal L_{g_i}\big)\mathcal U(t)\); in lifted coordinates this induces bilinear evolution, while in sampled-data and weak-input regimes it is often approximated by affine input models such as \(z_{k+1}=A z_k+B u_k\). The factorization organizes controlled Koopman modeling around fixed autonomous and input-channel operators, and underlies structured identification, LPV and bilinear realizations, model predictive control, reinforcement learning, and several physics-informed or robust extensions [2102.02522] [2211.07112] [2207.12132].

## 1. Operator-theoretic foundation

In the autonomous discrete-time setting, for a map \(F:x\mapsto f(x)\) and an observable \(g\), the Koopman operator acts by composition, \((Kg)(x)=g(f(x))\). In autonomous continuous time, for \(\dot x=f(x)\), the infinitesimal generator acts as the Lie derivative \((Lg)(x)=\nabla g(x)\cdot f(x)\). For controlled systems, one either considers a family of controlled Koopman operators \(K(u)\) in discrete time or a controlled generator \(L(u)\) in continuous time. In the control-affine case \(\dot x=f(x)+\sum_{i=1}^m g_i(x)u_i\), the generator splits as \(L(u)=L_0+\sum_{i=1}^m u_iL_i\), with \(L_0 g=\nabla g\cdot f\) and \(L_i g=\nabla g\cdot g_i\) [2102.02522].

A more geometric formulation writes the Koopman evolution itself as an operator-valued differential equation:
\[
\frac{d}{dt}\mathcal U(t)=\Big(\mathcal L_f+\sum_{i=1}^m u_i(t)\mathcal L_{g_i}\Big)\mathcal U(t),\qquad \mathcal U(0)=\mathrm{Id}.
\]
Here \(\mathcal L_f\) and \(\mathcal L_{g_i}\) are Lie derivatives on \(C^\infty(M)\), and their commutators satisfy \([\mathcal L_a,\mathcal L_b]=\mathcal L_{[a,b]}\). This gives a global, coordinate-free bilinearization on an infinite-dimensional Lie group of Koopman operators, rather than a local state-space linearization. In this formulation, Cartan’s formula \(\mathcal L_v=d\,\iota_v+\iota_v\,d\) connects Koopman generators to de Rham operators, which in turn yields controllability statements and an infinite-dimensional extension of the Lie algebra rank condition [2211.07112].

This operator-level viewpoint is the conceptual core of the subject. It identifies control not as an external add-on to an autonomous Koopman model, but as a structured perturbation through fixed generators or operator slices. That distinction is decisive for both theory and computation.

## 2. Canonical lifted forms

The generator factorization induces specific lifted state-space forms. In eigenfunction or observable coordinates \(\phi\), continuous-time control-affine systems admit bilinear evolution
\[
\dot \phi=\Lambda \phi+\sum_{j=1}^m u_j B_j \phi,
\]
where \(\Lambda\) contains autonomous eigenvalues and \(B_j\) encodes the action of the control generators. For sampled-data systems with piecewise-constant inputs, integrating the generator yields
\[
z_{k+1}=A z_k+\sum_{j=1}^m u_{k,j}N_j z_k.
\]
A simpler exogenous-input approximation,
\[
z_{k+1}=A z_k+B u_k,
\]
is often used when input nonlinearity is weak, but the review literature emphasizes that assuming \(\dot z=A z+B u\) in lifted coordinates is generally only local unless bilinear terms are included; control-affine systems generically admit bilinear realizations [2102.02522].

Earlier generalizations of Koopman theory to inputs and control cast the lifted predictor as
\[
\psi(x_{k+1})\approx K_{xx}\psi(x_k)+K_{xu}\phi(u_k),
\]
with a larger observable space \(H_X\otimes H_U\otimes H_{XU}\) to accommodate state-only, input-only, and mixed observables. This framework recovers DMDc when identity observables are used and extends it by allowing mixed terms such as \(x_i u_j\) or nonlinear products to appear in the regressor while keeping the output restricted to state observables [1602.07647].

A central correction to the naïve lifted-LTI picture is that exact controlled Koopman forms are typically LPV or bilinear rather than LTI. For continuous-time systems with inputs, the lifted dynamics take the form
\[
\dot \Phi(x)=A\Phi(x)+B(x,u)u,
\]
and for control-affine continuous-time systems the dependence simplifies to \(B(x)u\). In discrete time, the input matrix is generally state and input dependent. The resulting representation has linear state transition in the lifted coordinates but a specially structured LPV input map. Error bounds are available for approximating this structure by a constant \(\hat B\), and these bounds quantify when an LTI Koopman approximation is acceptable and when LPV or bilinear control should be preferred [2207.12132].

The main misconception addressed across this literature is therefore straightforward: a controlled nonlinear system does not, in general, become an exact lifted LTI system merely by adding a constant \(B\) matrix.

## 3. Structured finite-dimensional approximation

Finite-dimensional Koopman control factorization is learned from data by choosing dictionaries or kernels that preserve the autonomous/input-channel split. In EDMDc and related regressions, one lifts \(x_k\) to \(z_k=\Psi(x_k)\), forms snapshot matrices \(Z_-=[z_0\ \cdots\ z_{T-1}]\), \(Z_+=[z_1\ \cdots\ z_T]\), and inputs \(U=[u_0\ \cdots\ u_{T-1}]\), then solves
\[
[A\ B]=Z_+[\,Z_-\ \ U\,]^\dagger.
\]
To learn bilinear models, the regressor is augmented with cross-terms \(u_{k,j}z_k\), so that \(A\), \(B\), and \(N_j\) are estimated simultaneously. The same literature distinguishes structured approximations, which impose factorization and spectral or invariance constraints, from unstructured large-dictionary fits, which are described as data-inefficient, locally predictive, and spectrally misleading [2102.02522].

The original Koopman-with-inputs-and-control formulation makes the same point in operator language. The regressor may include \(\psi(x)\), \(\phi(u)\), and mixed observables \(\chi(x,u)\), while the output space is restricted to state observables. This allows closure without forcing future inputs to be predicted when they are exogenous. In the SIR-with-vaccination example, including the mixed term \(SI\) in the input space but not in the output space avoids an otherwise unclosed lifted system [1602.07647].

A nonparametric extension replaces fixed dictionaries by RKHS machinery. In the control-Koopman operator regression framework, the joint state-input kernel
\[
k((x,u),(x',u'))=k(x,x')+\sum_{j=1}^{n_u}u_j u'_j k(x,x')
\]
induces a single operator regression \(G:H_u\to H_x\). The learned lifted predictor is bilinear,
\[
z_{t+1}=A^T z_t+\sum_{j=1}^{n_u} B_j^T u_j(t)\,z_t,
\]
but the input dependence is pushed into the domain kernel rather than parameterized by an explicit tensor-product basis. The resulting kernel ridge estimator yields \(A_\gamma\) and \(B_{\gamma,k}\) from one stabilized Gram inversion, and Nyström sketching reduces the complexity to \(O(m^3+m^2N)\), independent of input dimension \(n_u\) [2405.07312].

Across these constructions, dictionary choice, persistence of excitation, and closure under generator action remain the decisive practical issues. The factorization is useful only if the chosen lifted coordinates actually support it.

## 4. Control synthesis and optimization

Once a factorized lifted model is available, standard and nonstandard control methods become available in the lifted coordinates. For affine lifted models \(z_{k+1}=A z_k+B u_k\), the literature uses discrete LQR with the usual Riccati equation and feedback \(u_k=-K_c z_k\). For bilinear models \(\dot z=A z+\sum_i u_i B_i z\) or \(z_{k+1}=A z_k+\sum_i u_{k,i}B_i z_k\), nonlinear MPC or bilinear optimal control is used instead [2110.08442] [2105.08036].

A canonical example is Koopman NMPC for control-affine robotics. In a planar quadrotor study, the open-loop prediction MSEs were \(8.71\times 10^{-2}\) for DMD, \(5.60\times 10^{-2}\) for EDMD, and \(7.53\times 10^{-3}\) for a bilinear EDMD model. In closed loop, realized cost normalized by NMPC was \(1.04\pm0.05\) for DMD, \(1.02\pm0.03\) for EDMD, and \(1.00\pm0.00\) for Koopman NMPC. The key computational advantage was that SQP linearizations of the bilinear lifted model use only matrix additions and multiplications with the learned \(A\) and \(B_i\) blocks [2105.08036].

The geometric control line yields a different synthesis route: Koopman feedback linearization. If an output has a relative degree on the controllable submanifold, then observables built from the output and its successive Lie derivatives produce a lifted linear closed-loop system \(\frac{d}{dt}(\mathcal U(t)\phi)=A\,\mathcal U(t)\phi+B\,v(t)\). The feedback law is localized to the controllable submanifold and does not require global controllability of the original system on the full manifold [2211.07112].

Koopman factorization has also entered reinforcement learning and predictive control. Koopman-assisted reinforcement learning parameterizes an action-conditioned operator \(K^u\) by a Koopman tensor, so that the Bellman continuation term becomes \(w^\top K^u\phi(x)\), enabling soft value iteration and a Koopman-assisted SAC critic [2403.02290]. In data-driven predictive control, the lifted state/control separation appears as \(z_{k+1}=\mathfrak A z_k+\mathfrak B u_k\), but the actual controller avoids explicit operator identification and instead combines Koopman lifting with Willems’ fundamental lemma and a bilevel prediction-control scheme [2102.05122]. Deep Koopman models with control introduce a learned latent control variable \(v_t=g_\psi(x_t,u_t)\) and then use
\[
z_{t+1}=K z_t+B v_t,
\]
so that state-dependent or nonlinear input effects are absorbed into the control encoding rather than forced into a fixed \(B u_t\) term [2202.08004].

## 5. Specialized variants

Several later developments refine what is meant by factorization. One direction uses set-valued analysis. For controlled systems, the set-valued Koopman semigroup is
\[
K_{(\tau,t)}(\varphi)=\{\varphi\circ\Phi_{(\tau,t)}^u:\ u(\cdot)\in\mathcal U\},
\]
and its generator is the set-valued Liouville operator
\[
L(\varphi)=\{\nabla \varphi\cdot f_u:\ u\in U\}.
\]
For control-affine systems this yields \(L=\mathcal L_0+\mathrm{cone}\{\mathcal L_1,\dots,\mathcal L_m\}\), which formalizes the practice of bundling Liouville operators across admissible controls and leads to set-valued spectral mapping results [2401.11569].

Another direction focuses on physical coherence of the input matrix. In control-coherent Koopman modeling, actuator states are explicitly included in the observables, the autonomous Koopman operator is learned with control removed from the actuator update, and the exact lifted input matrix is then imposed as
\[
B=\begin{bmatrix}B_p\\0\end{bmatrix}.
\]
This construction is designed so that inputs act only on actuator observables in one step, while non-actuator coordinates respond only through autonomous propagation. On multi-degree-of-freedom robotic arms, reported tracking errors on circular trajectories were markedly smaller for the control-coherent model than for DMDc with a least-squares \( \hat B \) [2403.16306].

A distinct recent use of the term refers to closed-loop controller parameterization. For discrete-time control-affine lifts, if the algebraic compatibility condition \((S\psi_x(x))\otimes\psi_u(x)=H\psi_x(x)\) holds and the plant admits a bilinear Koopman lift, then feedback \(u_k=K_u\psi_u(x_k)\) produces a closed-loop lifted operator
\[
\tilde K=K_{xx}+K_{xu}(I_{d_S}\otimes K_u)H.
\]
For fixed system matrices, this is affine in the controller matrix \(K_u\), which allows Lyapunov-stable controller synthesis by semidefinite programming with a fixed \(P\) [2510.05359].

The same factorization philosophy has also been extended to multiscale and stochastic settings. For hierarchical control with time-scale separation, observables are partitioned into slow, fast, actuator, and supervisory components, and fast equilibria are analytically collapsed into reduced slow operators such as
\[
K_{\mathrm{comb},xx}=K_{xx}+(I-K_{xx})K_{xy}K_{yx},
\]
which are then used for cross-scale stability analysis and optimal control at the slow time scale [2506.15586]. For control-affine stochastic systems, the Itô generator splits as \(L=L_0+\sum_i u_iL_i\), with \(L_0\) including the diffusion term; training on constant inputs \(u\in\{0,e_1,\dots,e_m\}\) then yields \(L_i=L_{e_i}-L_0\), which supports bilinear lifted prediction and optimal control of metastable stochastic dynamics [2410.09452].

## 6. Limitations, misconceptions, and open problems

The literature is unusually explicit about limitations. First, the affine lifted model \(z_{k+1}=A z_k+B u_k\) is not a universal controlled Koopman form. For control-affine systems, bilinear or LPV realizations are generally the faithful finite-dimensional models; in continuous time the input matrix is state dependent, and in discrete time it is generally state and input dependent [2207.12132] [2102.02522].

Second, finite-dimensional closure is fragile. Large fixed dictionaries are described as data-inefficient and prone to ill-conditioning; fixed unstructured bases often fail to deliver Koopman-invariant subspaces; measurement noise distorts spectra; and regression estimates for \(B\) or \(N_j\) require persistently exciting inputs and careful experimental design [2102.02522]. Even in nonparametric RKHS formulations, extrapolation outside the training domain remains limited by kernel generalization, and exact bilinearization is proved only for control-affine systems on non-recurrent domains [2405.07312].

Third, not all variants resolve the same issue. Set-valued Koopman theory provides existential spectra over control sets rather than robust “for all \(u\)” spectra [2401.11569]. Deep latent-control factorizations improve prediction when actuation is strongly nonlinear, but general nonlinear control encodings may not be invertible, which prevents direct recovery of physical inputs from latent control variables [2202.08004]. Control-coherent models remove one major structural error by fixing the lifted \(B\) from actuator physics, but they rely on the availability of explicit actuator states and a linear local actuator model [2403.16306].

The open problems identified across these works are consistent: better dictionary or kernel design, stronger invariance and stability constraints, robust spectra and uncertainty handling across control sets, adaptive or online identification, and controller synthesis that retains the algebraic advantages of factorization without sacrificing fidelity. Taken together, these works place Koopman control factorization not as a single algorithm, but as a structural principle: control is most effective in Koopman coordinates when autonomous evolution and input channels are represented separately, and when that separation is respected by learning, reduction, and synthesis.

Source: https://www.emergentmind.com/topics/koopman-control-factorization