---
title: Koopman–Nemytskii Operator in Nonlinear Control
url: https://www.emergentmind.com/topics/koopman-nemytskii-operator
type: topic
---

# Koopman–Nemytskii Operator in Nonlinear Control

The Koopman–Nemytskii operator is a bounded linear operator that generalizes the classical Koopman operator to nonlinear controlled systems with feedback. Unlike the standard Koopman operator, which provides a linear but infinite-dimensional representation for autonomous dynamical systems, the Koopman–Nemytskii operator explicitly incorporates nonlinear state-feedback laws, enabling linearization of closed-loop dynamics in reproducing kernel Hilbert spaces (RKHSs). This construction facilitates data-driven predictive modeling and controller synthesis for a broad class of smooth nonlinear control systems via operator learning methodologies such as kernel extended dynamic mode decomposition (kernel EDMD), with rigorous error guarantees [2503.18269].

## 1. Mathematical Formulation and Operator Definition

Let $X\subset \mathbb{R}^{n_x}$ represent the compact state domain and $U$ a compact family of feedback laws $u:X\to\mathbb{R}^{d_u}$. Take $\kappa$ as a Mercer (Wendland) kernel on $X$, yielding RKHS $\mathcal{N}_\kappa(X)\simeq H^s(X)$ for Sobolev index $s$, and $\varpi$ as a universal kernel on $U$ with $\mathcal{N}_\varpi(U)$. The canonical feature maps are
- $\phi:X\to\mathcal{N}_\kappa(X)$, $\phi(x) = \kappa(x,\cdot)$,
- $\psi:U\to\mathcal{N}_\varpi(U)$, $\psi(u) = \varpi(u,\cdot)$.

Define $f_u(x)=f(x,u(x))$ as the closed-loop map. The Koopman–Nemytskii operator
$$
T:\mathcal{N}_\kappa(X)\times\mathcal{N}_\varpi(U)\to\mathcal{N}_\kappa(X)
$$
is the unique bounded linear map characterized by
$$
T(\phi(x), \psi(u)) = \phi(f_u(x)),\quad\forall (x,u)\in X\times U.
$$
Alternatively, equating $\mathcal{N}_\kappa(X)\times\mathcal{N}_\varpi(U)$ with $\tilde{\mathcal{N}}$ on $X\times U$, one writes
$$
\bar{T}\,\bar{\phi}(x,u)=\phi(f_u(x)),\qquad \bar{\phi}(x,u)=\phi(x)\otimes\psi(u).
$$
This operator provides a universal, linear embedding of the family of nonlinear closed-loop maps parameterized by the feedback law $u$.

## 2. Connections to Koopman, Nemytskii, and Operator Lifting

The Koopman–Nemytskii operator unifies several classical nonlinear operator concepts:

- **Autonomous Koopman operator:** $K_f:C(X)\to C(X)$, $K_fg = g\circ f$, describing autonomous dynamics.
- **Nemytskii (substitution) operator:** $N_f:C^s(X, \mathbb{R}^{d_u})\to C^s(X, \mathbb{R}^{n_x})$, $N_f(u) = f(\cdot, u(\cdot)) = f_u$, acting on maps from state to action.
- **Koopmanizing map:** $M:C^{s+1}(X, \mathbb{R}^{n_x})\to L(H^{s+1}(X), H^s(X))$, $M(f) = K_f$.

Their composition $M N_f:u\mapsto K_{f_u}$ describes a nonlinear lifting from the feedback law space $U$ into the space of operators $L(H^{s+1}, H^s)$. By lifting $U$ to the RKHS and passing to adjoints, one constructs a linear operator $T_1:\mathcal{N}_\varpi(U)\to L(H^s(X), H^s(X))$, $T_1(\psi(u)) = K_{f_u}^*$. Identifying $L(H^s, H^s)\simeq L(H^s\otimes\mathcal{N}_\varpi(U), H^s)$ realizes the bilinear form $T:H^s(X)\times \mathcal{N}_\varpi(U)\to H^s(X)$, thereby "linearizing" the dependency on both state and policy [2503.18269].

## 3. Functional Analytic Foundations and RKHS–Sobolev Equivalence

Under appropriate regularity assumptions on the system dynamics and kernels—specifically, $f_u\in C^{s+1}(X)$ for all $u\in U$, the determinant of the Jacobian $\det D_x f_u(x)$ bounded away from zero, and choosing Wendland kernel $\kappa_{n,k}$ on $X\subset\mathbb{R}^n$ with smoothness $s=(n+1)/2+k$—there holds the isomorphism $\mathcal{N}_{\kappa_{n,k}}(X)\simeq H^{(n+1)/2+k}(X)$ [2503.18269]. This establishes well-posedness for the bounded Koopman operator $K_{f_u}:H^{s+1}(X)\to H^{s}(X)$ and, by extension, for $T\in L(\mathcal{N}_\kappa(X)\times\mathcal{N}_\varpi(U),\mathcal{N}_\kappa(X))$.

This guarantees the legitimacy of the feature space embedding, which underpins both theoretical analysis and practical computations, facilitating representation, learning, and error estimation in a functional analytic setting.

## 4. Data-driven Learning: Kernel EDMD Approximation

Empirical approximation of the Koopman–Nemytskii operator leverages kernel EDMD methodologies. Given data samples $\{(x_i, u_i, y_i)\}_{i=1}^m$ with $y_i=f(x_i,u_i(x_i))$, construct the Gram matrices:
- $G_{yy} = [\kappa(y_i,y_j)]_{i,j=1}^m$,
- $G_{xu} = [\kappa(x_i,x_j)\cdot\varpi(u_i,u_j)]_{i,j=1}^m$.

By the representer theorem, any Hilbert–Schmidt estimator $\hat{T}$ admits a finite-rank expansion:
$$
\hat{T} = \sum_{i,j=1}^m \theta_{ij}\phi(y_i)\otimes\bar{\phi}(x_j,u_j),
$$
where $\Theta=(\theta_{ij})$ minimizes the empirical loss with (or without) regularization.

Approximation modalities include:
- **Kernel EDMD (no regularization):** $\Theta=G_{xu}^{-1}$; $\hat{T}=TS$, $S$ projects onto $\operatorname{span} \{\bar{\phi}(x_i,u_i)\}$.
- **Reduced-rank regression (regularized, rank $\leq r$):** Solve the generalized eigenproblem
  $$
  \frac{1}{m^2}G_{yy}G_{xu}v = \sigma^2\left(\frac{1}{m}G_{xu}+\beta I\right)v,
  $$
  extract the $r$ leading modes $V=[v_1\ldots v_r]$, and set $\Theta = \frac{1}{m}VV^TG_{xu}$.

This operator learning framework exploits the high expressiveness of RKHSs while maintaining computational tractability through matrix representations.

## 5. Error Bounds and Predictive Guarantees

The Koopman–Nemytskii operator framework yields explicit, nonasymptotic error bounds for state prediction and cost approximation:

- **Generalization for reduced-rank regression (Theorem 4.1):** With probability at least $1-\delta$,
  $$
  \ell(\hat{T}) - \hat{\ell}_0(\hat{T}) \leq \frac{1}{m}\log\frac{6}{\delta} + \sqrt{\frac{8}{m}\log\frac{6}{\delta}} + \beta(\beta+2\sqrt{r})\left[\frac{6}{m}\log\frac{12m^2}{\delta} + \sqrt{\frac{9}{m}\log\frac{12m^2}{\delta}}\right].
  $$
- **Single-step prediction (Thm 4.3):** Under an interior-cone condition and mesh fill-distance $\eta_{X\times U}$,
  $$
  \|T - \hat{T}\| \leq \|T\|c\eta_{X\times U}^{k+1/2}.
  $$
- **Multi-step prediction (Thm 4.4):** Writing $\phi_x^0 = \phi(x)$, $\phi_x^{t+1}=T(\phi_x^t\otimes\psi(u))$ and similarly for $\hat{T}$, then with $\beta=\|T\|+c\eta^{k+1/2}$,
  $$
  \|\hat{\phi}_x^t - \phi_x^t\| \leq \frac{1-\beta^t}{1-\beta}c\eta^{k+1/2},\quad \text{if }\beta<1,
  $$
  which gives a uniform bound for all $t$.
- **Accumulated cost (Thm 4.5):** Let $Q$ (on $\mathcal{N}_\kappa(X)$) and $R$ (on $\tilde{\mathcal{N}}(X\times U)$) be kernel-quadratic costs, discount $\gamma$, with $c_Q=\|Q\|$, $c_R=\|R\|$. Then
  $$
  |\hat{\psi}(x,u)-\psi(x,u)| \leq \sum_{t=0}^\tau \gamma^t(c_Q+c_R)\left[2c\eta^{k+1/2}\beta^t\frac{1-\beta^t}{1-\beta} + \left(c\eta^{k+1/2}\frac{1-\beta^t}{1-\beta}\right)^2\right].
  $$

These results provide quantitative control over the performance of data-driven operator approximations for both prediction and cost evaluation tasks.

## 6. Applications in Nonlinear Closed-Loop Control

Upon learning the Koopman–Nemytskii operator $\hat{T}$, several functionalities become available:

- Predict closed-loop trajectories in feature space, then reconstruct state-space quantities via $\langle\phi, \mathrm{Id}\rangle$.
- Approximate cost-to-go $\psi(x,u)$ for any $u\in U$.
- Optimize control policies using grid search or gradient-based methods in $U$.

Two application domains are detailed in [2503.18269]:
- **Single-tank level control:** For $d_x=1$, $U$ parameterized by feedback gain $k$ in $u(x)=-k x$. Utilizing Wendland kernels on $X$, Gaussian kernel on $U$, kernel EDMD accurately predicts single-/multi-step states as well as cumulative quadratic cost functions on test grids.
- **Williams–Otto reactor:** For $d_x=6$, $U$ parameterized by feedback $u=k_1 x_3 + k_2 x_6$ over a specified range. Using $m\approx 1000$ closed-loop samples, Wendland kernels on $X$, Gaussian on $U$, the learned $\hat{T}$ achieves low error in state prediction ($x_6(t)$ for $t=4,8,16,32$) and accumulated cost approximation across policies.

This demonstrates the utility of Koopman–Nemytskii operator learning for both trajectory synthesis and policy optimization in smooth nonlinear systems with state-dependent feedback.

## 7. Significance and Theoretical Implications

The Koopman–Nemytskii operator framework provides an infinite-dimensional, linear representation of closed-loop, nonlinear controlled dynamics with explicit feedback law dependence. This framework overcomes limitations of prior extensions of the Koopman operator, such as the restriction to autonomous systems or open-loop input dependence only. By employing Sobolev-compatible kernels and leveraging RKHS–Sobolev equivalence, the theory ensures rigorous operator well-posedness and favorable generalization error properties.

A plausible implication is the potential for principled, data-driven control synthesis and evaluation for nonlinear systems in regimes where classical linear or hybrid representations are inadequate. The approach accommodates rich policy classes, exploits regularity of nonlinearities, and provides explicit, nonasymptotic error quantification, setting a foundation for robust nonlinear control in data-centric settings [2503.18269].

Source: https://www.emergentmind.com/topics/koopman-nemytskii-operator