---
title: Bilinear EDMD Control Schemes
url: https://www.emergentmind.com/topics/bilinear-edmd-based-control-schemes
type: topic
---

# Bilinear EDMD Control Schemes

Bilinear Extended Dynamic Mode Decomposition (EDMD)-based control schemes are data-driven methods for representing and controlling nonlinear, control-affine dynamical systems by lifting their evolution into a higher-dimensional space of observables where the influence of the control input is modeled in a bilinear form. This class of techniques leverages finite-dimensional approximations of the infinitesimal Koopman operator to create surrogate models that robustly capture state-input interactions, address the limitations of purely linear Koopman formulations, and facilitate tractable, high-performance control algorithms—often with certified closed-loop guarantees—even in the presence of model uncertainty and finite-data effects.

## 1. Mathematical Foundations of Bilinear Koopman Realizations

The framework is built upon the Koopman operator, a linear (infinite-dimensional) operator acting on observables of the state. For a control-affine nonlinear system
\[
\dot{x}(t) = f_0(x(t)) + \sum_{j=1}^{m} f_j(x(t)) u_j(t)
\]
the classical Extended Dynamic Mode Decomposition (EDMD) lifts $x$ to a vector of observables $z = \Psi(x)$, leading to an evolution
\[
\Psi(x_{k+1}) \approx A \Psi(x_k) + B_0 u_k + \sum_{j=1}^m u_{j,k} B_j \Psi(x_k)
\]
where $A, B_0, B_j$ are matrices estimated from data, with the bilinear terms $u_j B_j \Psi(x_k)$ encoding the interaction between lifted state and input.

This bilinear structure arises naturally when finite-dimensional projections of the control-affine Koopman generator onto an observable space are formed. Finite approximations rarely admit a purely linear structure with respect to the control, but, provided the observables are sufficiently rich, they always admit a bilinear realization for control-affine systems [2010.09961], [2202.09084], [2509.02839].

A necessary and sufficient condition for a set of observables $\{z_i\}$ to admit a bilinear Koopman realization is
\[
\frac{\partial z_i}{\partial x} \cdot F(x,u) \in \mathrm{span} \big( \{z_j\} \cup \{u_j\} \cup \{z_k u_j\} \big)
\]
for all $i$ [2010.09961]. This allows the time evolution of each observable to be written as a combination of lifted state, control, and their products, yielding the core bilinear form.

## 2. System Identification and Data-Driven Model Construction

Bilinear EDMD-based models are typically identified via regression using snapshot pairs (state and input evolution data). For each input direction, data is collected with the control held at either zero or a canonical value, and multiple autonomous Koopman operators or generators are estimated [2202.09084], [2409.18010], [2509.02839]:

- For each autonomous subsystem (i.e., at $u=0$ and $u=e_j$), the Koopman generator/operator $L^0$ and $L^{e_j}$ is estimated.
- The generator for arbitrary $u$ is constructed as a bilinear combination:
  \[
  L^u = L^0 + \sum_{j=1}^{m} u_j (L^{e_j} - L^0)
  \]
  ensuring the bilinear dependence is retained while maintaining model parsimony and mitigating the curse of dimensionality [2202.09084], [2408.12277].

This approach preserves the low input complexity of the original system and allows the construction of a global data-driven surrogate without state-control dimensional augmentation.

## 3. Error Bounds, Robustness, and Closed-Loop Guarantees

Finite data and projection onto finitely many observables introduce a structured model error, typically expressed as a residual $r(x,u)$:
\[
\Psi(x_{k+1}) = A \Psi(x_k) + B_0 u_k + \sum_{j=1}^m u_{j,k} B_j \Psi(x_k) + r(x_k, u_k)
\]
Recent advances provide proportional, uniform error bounds of the form [2202.09084], [2411.03875], [2411.10359], [2509.02839]:
\[
\| r(x,u) \|_2 \leq c_x \|\Psi(x)\|_2 + c_u \|u\|_2
\]
where $c_x, c_u$ are functions of the data richness, regularity of the dictionary, dynamics, and number of samples $m=O(M^2/(\epsilon^2\delta))$ for accuracy $\epsilon$ and confidence $1-\delta$.

Controllers are designed to account for these errors, embedding them directly into robust control designs (such as sum-of-squares (SOS) or LMI-based optimization). Feasibility of the resulting convex program establishes rigorous closed-loop certificates (e.g., exponential stability within a certified region of attraction) given explicit $c_x, c_u$ bounds [2411.03875], [2411.10359]. In min-max MPC formulations, a rational state-feedback (SOS-parametrized) law achieves robust constraint satisfaction and cost bounds [2504.04870].

## 4. Controller Synthesis: Algorithmic Strategies

Controller design is performed on the bilinear surrogate with explicit error bounds embedded. Two principal strategies appear:

- **Convex/robust control synthesis:** LMI-based (or SOS-based) programs are constructed to robustly stabilize
  \[
  \Psi(x_{k+1}) = A\Psi(x_k) + B_0 u_k + \sum_j u_{j,k} B_j \Psi(x_k) + r(x_k, u_k)
  \]
  with $r(x,u)$ bounded as above. The classic synthesis yields controllers of the form [2411.03875], [2411.10359]:
  \[
  u(x) = [I - L_w (\Lambda^{-1} \otimes \Psi(x))]^{-1} L P^{-1} \Psi(x)
  \]
  or, in rational-SOS form,
  \[
  u(x) = \frac{L_n(\Psi(x)) P^{-1} \Psi(x)}{u_d(\Psi(x))}
  \]
  with $L_n$ and $u_d$ SOS polynomials ensuring global or semi-global robust stability.

- **Predictive control (MPC/DeePC):** Bilinear Koopman surrogates enable efficient implementation of predictive controllers. In model predictive control, the surrogate model replaces the true dynamics in the horizon updates, allowing for convex (when locally linearized) or bilinear constrained optimization [2010.09961], [2505.03346]. In DeePC frameworks, the bilinear structure is exploited directly in the data, removing explicit model identification and using Hankel data and bilinear consistency constraints [2505.03346].

When properly accounting for $r(x,u)$, such controllers guarantee constraint satisfaction and closed-loop stability within robustness margins set by the data-derived error bound [2411.03875], [2504.04870].

## 5. Scalability, Curse of Dimensionality, and Modularization

Bilinear models circumvent the bulk of the curse of dimensionality: unlike state-augmented EDMDc approaches, the model's dimension depends linearly on the number of state variables and affine in the number of control channels [2202.09084]. Modularized schemes for interconnected systems further decompose the problem: local Koopman generators for each subsystem are learned individually on lower-dimensional state spaces, then globally coupled through known network topology [2408.12277]. This supports transfer learning, network adaptation, and error propagation analysis via graph-theoretic tools.

Empirical illustrations confirm that modularized EDMD-based surrogates require $\mathcal{O}(s \cdot e^{n_0})$ data (with $s$ subsystems each of dimension $n_0$) versus $\mathcal{O}(e^{s n_0})$ for a monolithic EDMD approach [2408.12277].

## 6. Applications and Limitations

Bilinear EDMD-based schemes have been applied to:

- **Robotics and mechatronic systems:** Real-time trajectory-following and MPC of robots; high-accuracy predictive control with computational efficiency [2010.09961], [2112.07983], [2411.07192].
- **Quantum control:** Optimal control pulse synthesis and system identification leveraging known Hamiltonian structure [2010.14577].
- **Modular large-scale networks:** Coupled oscillators and multi-component physical systems, where scalability and transfer learning are critical [2408.12277].
- **Data-driven inverse optimal control:** Estimation of cost and dynamics via bilinear Koopman representations, enabling a tractable extension of the classical inverse LQR solution to unknown nonlinear systems [2501.18318].

However, precise dictionary selection remains crucial; improper choice reduces model accuracy or increases conservatism in robust control [2509.02839]. For systems with geometric structure (e.g., nonholonomic robots), domain-specific knowledge must inform the selection of observables and cost functions to achieve high-performance closed-loop behavior [2411.07192]. The data requirements implied by theory may be conservative. Scaling computational tools (especially SOS or LMI-based synthesis) to high-dimensional lifted spaces remains an area for future progress.

## 7. Theoretical Implications and Future Directions

Bilinear Koopman realizations provide a systematic, universally applicable surrogate for control-affine systems, bridging data-driven system identification and robust controller synthesis. This structurally enforces that, as dictionary richness and data volume increase, model fidelity and control performance improve monotonically (in contrast to the stagnation of purely linear lifts) [2010.09961], [2202.09084], [2509.02839].

Future directions include:

- **Automated/learned observable selection**—potentially using neural networks or adaptive dictionaries—to optimize tradeoffs between expressivity and computational tractability.
- **Extension to output-only (partial observation) settings** and the development of observer architectures with guaranteed error bounds.
- **Integrated data-driven robust performance synthesis** with joint consideration of state/input constraints, performance indices, and model uncertainty.
- **Scalable algorithms** for high-dimensional, real-time applications—especially regarding the efficient solution of large-scale SOS or LMI programs.
- **Hybrid pipelines** combining deep learning with operator-theoretic error quantification and geometric insights for systems with inherent structure.

Bilinear EDMD-based control thus provides an algorithmically flexible, theoretically grounded platform for robust, data-driven nonlinear control with quantifiable performance under realistic data and computational constraints.

Source: https://www.emergentmind.com/topics/bilinear-edmd-based-control-schemes