---
title: Generalized Extended Dynamic Mode Decomposition (gEDMD)
url: https://www.emergentmind.com/topics/generalized-extended-dynamic-mode-decomposition-gedmd
type: topic
---

# Generalized Extended Dynamic Mode Decomposition (gEDMD)

Generalized Extended Dynamic Mode Decomposition (gEDMD) denotes a class of operator-theoretic, data-driven approximation methods that extend extended dynamic mode decomposition (EDMD) beyond a fixed, hand-chosen dictionary of observables. In the broader literature summarized here, gEDMD usually refers to EDMD with more flexible function spaces, regularization, multiple observables, learned or adaptive dictionaries, and, in an important specialization, approximation of the Koopman generator rather than only the Koopman operator [2108.01490] [2405.00539]. The unifying viewpoint is Galerkin projection: nonlinear dynamics are lifted to a finite-dimensional space of observables, and the action of a transfer operator is approximated there by a matrix learned from data.

## 1. Operator-theoretic formulation

For a discrete-time system
$$
x^+ = F(x),
$$
the Koopman operator acts on observables \(g\) by composition,
$$
Kg(x)=g(F(x)).
$$
For a continuous-time flow \(\Phi^t\), the Koopman semigroup is
$$
K^t g(x)=g(\Phi^t(x)),
$$
and, when the semigroup is strongly continuous, the generator is
$$
\mathcal{L}f(x)=\lim_{t\downarrow 0}\frac{1}{t}\big(K^t f(x)-f(x)\big).
$$
In the stochastic case
$$
dX_t=b(X_t)\,dt+\sigma(X_t)\,dW_t,
$$
the generator takes the Itô form
$$
\mathcal{L}f(x)=b(x)\cdot \nabla f(x)+\frac12\operatorname{Tr}\big(\Sigma(x)\nabla^2 f(x)\big),
$$
with \(\Sigma=\sigma\sigma^\top\) [2108.04433] [2405.00539].

EDMD begins by choosing a finite dictionary \(\Psi=\{\psi_1,\dots,\psi_N\}\), spanning a trial space \(F_N=\operatorname{span}(\Psi)\). In the classical Galerkin form, the Gram and structure matrices are
$$
[G_N]_{ij}=\langle \psi_i,\psi_j\rangle,\qquad [C_N]_{ij}=\langle A\psi_i,\psi_j\rangle,
$$
where \(A\) is the operator of interest. The finite-dimensional approximation satisfies
$$
(A_N^\Psi)^\top = C_N G_N^{-1}.
$$
In the unified Monte Carlo framework of Colbrook, Kostic, and collaborators, EDMD is the case \(A=K\) or \(K^t\), whereas gEDMD is the case \(A=\mathcal{L}\) or \(\mathcal{L}^*\) [2405.00539].

With snapshot data \(x_m\), empirical inner products produce
$$
[\hat G_{NM}]_{ij}=\frac1M\sum_{m=1}^M \psi_i(x_m)\overline{\psi_j(x_m)},\qquad
[\hat C_{NM}]_{ij}=\frac1M\sum_{m=1}^M A\psi_i(x_m)\overline{\psi_j(x_m)},
$$
and the data-driven approximation is
$$
(\hat A_{NM})^\top=\hat C_{NM}\hat G_{NM}^+.
$$
This subsumes both operator EDMD and generator EDMD within the same empirical Galerkin construction [2405.00539].

## 2. From EDMD to generalized EDMD

In its classical form, EDMD assumes that the chosen trial space is approximately Koopman-invariant. For a dictionary \(\{\psi_l\}_{l=1}^{N_o}\), one seeks a matrix \(\mathbf{K}_o\) such that
$$
\mathcal{K}^{\delta t}g \approx \sum_{l=1}^{N_o}\psi_l(\cdot)\,(\mathbf{K}_o^T a)_l
$$
for observables \(g=\sum_l a_l\psi_l\). With time-ordered data, the standard least-squares problem is
$$
\mathbf{K}_o=\arg\min_{\mathbf{K}}\|\Psi_+ - \mathbf{K}\Psi_-\|_F^2,
$$
with solution obtained from the SVD or, in Galerkin notation, by \(G^\dagger A\) [2108.04433] [2108.01490].

The central limitation is the dictionary. Classical EDMD requires the observables to be chosen a priori, and poorly chosen dictionaries can produce misleading or spurious spectral information when the span of the dictionary is not close to a Koopman-invariant subspace [2108.04433]. This is the primary motivation for gEDMD: the generalization lies not in abandoning the Galerkin structure, but in enlarging the admissible observable spaces, changing the projection norm or operator, and introducing regularization or data-adaptive feature maps [2108.01490] [2405.00539].

A second axis of generalization is output structure. In multi-output EDMD, one distinguishes the operator approximation from the representation of observables themselves. If \(\bm g=[g_1,\dots,g_p]^T\) is a vector observable and \(\bm W\) collects its coefficients in the dictionary basis, then one approximates
$$
\bm g^T = \bm\psi^T \bm W,
$$
and evolves all outputs by \(K^k\bm W\). This leads naturally to multi-task regression and to regularized estimation of both \(K\) and \(\bm W\) [2108.01490].

Regularization is not peripheral in gEDMD. For the Koopman matrix, ridge regularization gives
$$
K=(G+\beta I)^{-1}A.
$$
For multiple outputs, Tikhonov regularization with prior structure yields
$$
\bm W=(G+\bm Q)^{-1}(\bm B+\bm Q\bm W^0),
$$
allowing known observables, such as state coordinates, to be regularized differently from unknown sensor outputs [2108.01490]. This makes gEDMD a framework for constrained inverse problems on observable spaces, not merely an eigen-analysis routine.

## 3. Principal methodological branches

The modern literature uses the label gEDMD for a heterogeneous but structurally coherent family of methods. Some papers use the term explicitly; others fall into the same class because they retain EDMD’s projection-and-regression core while generalizing the dictionary, function space, or operator.

| Branch | Core generalization | Representative paper |
|---|---|---|
| Generator gEDMD | Approximate \(\mathcal{L}\) or \(\mathcal{L}^*\) instead of \(K\) | [2405.00539] |
| Multi-output EDMD | Jointly regress multiple observables and their iterates | [2108.01490] |
| Dictionary learning | Learn \(\Psi_\theta\) from data with neural networks | [1707.00225], [2110.01450] |
| Autoencoder / latent EDMD | Learn an embedding and perform EDMD in latent coordinates | [2108.04433] |
| Eigenfunction-based lifting | Use learned Koopman eigenfunctions as lifting functions | [1911.08751] |
| Analytic RKHS EDMD | Project in analytic RKHSs via Taylor or polynomial subspaces | [2405.15945] |
| Symmetry-constrained EDMD | Impose equivariance and group-convolution structure | [2411.00905] |
| Rigged DMD | Approximate generalized eigenfunctions and continuous spectrum | [2405.00782] |

Neural dictionary learning provides one of the most visible gEDMD instantiations. Li, Dietrich, Bollt, and Kevrekidis formulate a joint optimization
$$
(K,\theta)=\arg\min \sum_n \|\Psi(y_n;\theta)-K\Psi(x_n;\theta)\|^2+\lambda\|K\|_F^2,
$$
in which the dictionary \(\Psi(\cdot;\theta)\) is represented by a neural network and alternates with a regularized EDMD step [1707.00225]. Miyatake and collaborators replace the multilayer perceptron dictionary by a neural ODE, retaining the same EDMD-DL loss while reducing parameter count and memory footprint [2110.01450].

A closely related branch uses latent-coordinate learning. The DLDMD construction employs an encoder \(\mathcal{E}\) and decoder \(\mathcal{D}\), performs EDMD on encoded snapshots, and couples reconstruction, one-step EDMD residual, and multi-step prediction losses:
$$
\mathcal{L}=\alpha_1\mathcal{L}_{\text{recon}}+\alpha_2\mathcal{L}_{\text{dmd}}+\alpha_3\mathcal{L}_{\text{pred}}+\alpha_4\|\mathbf{W}_g\|_2^2.
$$
In the paper’s own interpretation, this searches for a learned feature map whose span is approximately Koopman-invariant and therefore sits squarely in the conceptual class of generalized EDMD [2108.04433].

KEEDMD generalizes EDMD in a different direction. Instead of arbitrary lifting functions, it learns a diffeomorphism \(h\) that pulls back eigenfunctions of a nominal linear model to approximate Koopman eigenfunctions of the nonlinear system, then uses those eigenfunctions as the EDMD dictionary. The lifted state is
$$
\mathbf{z}=\begin{bmatrix}\mathbf{x}\\ \boldsymbol{\phi}(\mathbf{x})\end{bmatrix},
$$
and the identified linear model has explicit structural blocks reflecting kinematics and eigenfunction evolution [1911.08751].

Two structurally specialized branches are especially notable. Analytic EDMD works in an RKHS of analytic functions, projects orthogonally onto polynomial subspaces, and exploits the triangular structure of the Koopman operator near a hyperbolic equilibrium. Because the projection is a Taylor projection in an analytic RKHS, the method does not suffer from spectral pollution and can achieve arbitrary spectral accuracy with fixed finite dimension by increasing data [2405.15945]. Group-convolutional EDMD, by contrast, generalizes EDMD under finite-group equivariance assumptions. Under those assumptions, the optimal EDMD matrix is equivariant, can be represented by a group convolution kernel, and becomes block-diagonal under the generalized Fourier transform, which yields data-efficient learning and fast eigenfunction approximation [2411.00905].

## 4. Spectral interpretation and function-space dependence

A defining issue in gEDMD is that the meaning of the computed spectrum depends on the function space. For analytic expanding circle maps, Slipantschuk, Wormell, and collaborators show that EDMD with Fourier-type dictionaries approximates compact Perron–Frobenius and Koopman operators on Hardy–Hilbert spaces of analytic functions, not merely a naive \(L^2\) Koopman operator. In that setting, the finite-dimensional matrices are genuine Galerkin approximations of compact operators, and eigenvalues converge exponentially with dictionary size [1905.09266].

A parallel result for chaotic circle maps clarifies the least-squares side. Wormell proves that, in the infinite-data limit and for trigonometric polynomial dictionaries, the least-squares projection error is exponentially small even under non-uniform analytic sampling measures, by means of a new approximation theorem in orthogonal polynomials on the unit circle. Forecasts and Koopman spectral data therefore converge exponentially fast with dictionary size to physically meaningful limits in this analytic setting [2305.08074].

These results sharpen a recurring lesson: gEDMD is not meaningful independently of the chosen observable space. If the function space is too large or poorly matched to the dynamics, continuous spectrum and spectral pollution can dominate. If it is chosen to reflect analyticity, hyperbolicity, or symmetry, the finite-dimensional approximation can inherit compactness, triangularity, or equivariance and thus acquire a legitimate spectral interpretation [1905.09266] [2405.15945] [2411.00905].

Continuous spectrum exposes the limitation of ordinary eigenpair extraction most starkly. Rigged DMD addresses this by combining measure-preserving EDMD with high-order kernels and a rigged Hilbert space \(\mathcal{S}\subset L^2\subset \mathcal{S}^*\). Instead of approximating only \(L^2\) eigenfunctions, it computes wave-packet approximations of generalized eigenfunctions through the resolvent and smoothed spectral projections, thereby covering both discrete and continuous spectral components [2405.00782]. A plausible implication is that gEDMD, in its most general form, includes not only richer dictionaries but also richer spectral objects.

## 5. Convergence, error, and robustness

The 2024 convergence theory of Colbrook, Kostic, and collaborators places EDMD and gEDMD inside a unified Monte Carlo Galerkin framework and proves convergence of the approximating operator and its spectrum under non-restrictive conditions [2405.00539]. For fixed dictionary size \(N\), the empirical matrix \(\hat A_{NM}\) converges almost surely to the Galerkin matrix \(A_N\) as the number of samples \(M\to\infty\). Under approximation assumptions on the trial spaces \(F_N\) and, for unbounded operators, on the graph-norm spaces \(D_N\), the Galerkin approximation converges strongly to the target operator as \(N\to\infty\) [2405.00539].

The same framework yields spectral convergence. For empirical eigenpairs \((\lambda_M,f_M)\) of \(\hat A_{NM}\), subsequences converge to eigenpairs of the Galerkin operator; in the dictionary limit, subsequences of Galerkin eigenpairs converge to eigenpairs of the full operator. In the joint limit of growing data and growing dictionary, subsequences of empirical eigenpairs converge to eigenpairs of the target operator \(A\), including the generator in the gEDMD case [2405.00539].

Quantitatively, the sampling error is governed by matrix concentration estimates for the empirical Gram and structure matrices. Up to logarithmic factors and basis-dependent quantities such as \(\kappa(G_N)\), the data complexity scales like \(\epsilon^{-2}\), and the empirical operator error decays like \(M^{-1/2}\) [2405.00539]. Noise changes constants but not the convergence order: the noisy-data estimates retain the same \(\epsilon^{-2}\) scaling, provided the noise is mean-zero and bounded in probability [2405.00539].

Conditioning remains decisive. The error bounds depend explicitly on \(\|G_N^{-1}\|\) and \(\kappa(G_N)\), and numerical experiments show that Gaussian and finite-element dictionaries are markedly more robust to noise than monomials, largely because of Gram-matrix conditioning [2405.00539]. This is consistent with the warning, already emphasized in the EDMD literature, that rich dictionaries can be numerically harmful when they are redundant or nearly linearly dependent [2108.01490] [1707.00225].

## 6. Applications, limitations, and recurrent misconceptions

The application range of gEDMD is broad. Learned-dictionary variants have been tested on the Duffing oscillator, Van der Pol oscillator, Lorenz-63, and the Kuramoto–Sivashinsky equation, where they outperform standard DMD or reduce parameter count substantially relative to fixed-dictionary baselines [2108.04433] [2110.01450]. KEEDMD improves state prediction and closed-loop trajectory tracking for a simulated cart-pole system by using learned Koopman eigenfunctions as lifting functions in an EDMD-style model [1911.08751]. Multi-output EDMD is aimed at settings where observables are sensor outputs rather than full states, including cases where some outputs are known and others are only observed through data [2108.01490]. Symmetry-aware group-convolutional EDMD has been demonstrated on the Kuramoto–Sivashinsky equation and a \(\lambda\)-\(\omega\) spiraling wave system, both nonlinear PDEs [2411.00905]. Rigged DMD has been applied to systems with Lebesgue spectrum, integrable Hamiltonian systems, the Lorenz system, and a high-Reynolds-number lid-driven flow [2405.00782].

Several misconceptions recur. First, gEDMD is not synonymous with deep learning. Neural dictionaries and autoencoder embeddings are prominent, but analytic RKHS projections, regularized multi-output formulations, symmetry-constrained EDMD, generator approximations, and rigged-Hilbert-space methods are equally part of the generalization [2405.15945] [2108.01490] [2405.00782]. Second, generalization does not remove the problem of function-space design. It redistributes it across basis choice, kernel choice, latent dimension, regularization, symmetry assumptions, and operator choice [2108.04433] [2405.00539]. Third, finite-dimensional Koopman models remain limited for strongly chaotic systems. In DLDMD, for example, Lorenz-63 attractor geometry is recovered reasonably well, but long-horizon pointwise prediction remains poor [2108.04433]. Rigged DMD can be read as a response to this limitation: when continuous spectrum dominates, generalized eigenfunction decompositions become more appropriate than ordinary eigendecompositions [2405.00782].

Open problems in the cited literature are correspondingly structural. They include automated selection or regularization of latent dimension in learned EDMD [2108.04433], extension of rigorous rates from i.i.d. sampling to time-correlated trajectory data [2405.00539], treatment of non-analytic or singular sampling measures in chaotic dynamics [2305.08074], extension of symmetry-constrained EDMD beyond finite groups or exact equivariance [2411.00905], and generalization of rigged spectral methods to non-unitary Koopman operators [2405.00782]. Across these branches, the central theme remains unchanged: gEDMD seeks finite, data-driven representations of transfer operators whose spectral content is both computationally accessible and dynamically meaningful.

Source: https://www.emergentmind.com/topics/generalized-extended-dynamic-mode-decomposition-gedmd