Jointly Lipschitz plant-conditioned operator
Implementation & benchmark of arXiv:2607.04362 — Approximate Feedback Linearization for a Nonlinear Hyperbolic PDE Class -- Part II: Neural Operator
Source paper: Approximate Feedback Linearization for a Nonlinear Hyperbolic PDE Class -- Part II: Neural Operator arXiv:2607.04362 ⓘ · analyzed Aug 29, 2026
AI-generated research hypothesis, automatically tested. Not peer-reviewed.
Idea description
Build a neural operator that maps both a function-valued plant description and the current state trajectory to an action, while explicitly controlling sensitivity to each input. This should improve interpolation to unseen dynamics and prevent large action changes caused by small state or plant-estimation errors.
Formulas
Mathematical statement
Theorem 1 states that the truncated controller operator \(\mathcal{U}_{N}\) satisfies \(|\mathcal{U}_{N}(\mathbf{f}_{N},u)-\mathcal{U}_{N}(\tilde{\mathbf{f}}_{N},\tilde u)|\le M_{N}(R)\|\mathbf{f}_{N}-\tilde{\mathbf{f}}_{N}\|_{\star}+M'_N(R)\|u-\tilde u\|_{L^{\infty}(0,1)}\). Here \(\mathbf f_N=(f_2,\ldots,f_N)\) is the collection of plant nonlinearity kernels, \(u\) is the state function on \([0,1]\), \(R\) bounds the state magnitude, \(\|\cdot\|_{\star}\) is the paper's kernel-coefficient norm, and \(M_N(R)=\sum_{n=2}^{N}R^n\lambda_n/n!\), \(M'_N(R)=\sum_{n=2}^{N}nR^{n-1}\kappa_n/n!\). The constants \(\lambda_n\) and \(\kappa_n\) bound the nth-order controller contribution with respect to plant kernels and state, respectively. The transferable property is separate, explicit sensitivity accounting for plant uncertainty and state perturbations.
Implementation notes
Integrate this at the policy or transition-model head. Represent each coefficient function \(f_n\) on a fixed spatial grid or basis, encode the concatenated coefficient fields with a branch network, encode the sampled state trajectory \(u(x_j)\) with a second branch network or Fourier neural operator, and combine the two embeddings through a DeepONet-style bilinear head to produce the scalar action or boundary control. During each minibatch, create pairs \((\mathbf f,u)\) and \((\tilde{\mathbf f},\tilde u)\) by perturbing coefficients and states within known radii. Compute \(d_f=\|\mathbf f-\tilde{\mathbf f}\|_{\star}\) and \(d_u=\|u-\tilde u\|_\infty\), then add the penalty \(\max(0,|\hat U-\hat U'|-L_fd_f-L_ud_u)^2\) to the ordinary imitation or prediction loss. Use spectral normalization in the branch and trunk layers, or estimate local gains with Jacobian-vector products. The paper's exact quantities are \(M_N(R)\) and \(M'_N(R)\); in a generic network, estimate \(\lambda_n\) and \(\kappa_n\) from perturbation slopes over training samples. First test on Burgers-like 1D systems with quadratic and cubic kernel nonlinearities, comparing an unconstrained MLP or FNO against this conditioned operator at equal parameter count. Measure action sensitivity under 1%, 5%, and 10% plant/state perturbations, held-out-plant prediction error, and closed-loop rollout failure rate. Success means lower worst-case action amplification and better stability without sacrificing nominal control cost.
Verification
This idea has not been verified yet.
Verification happens in two stages: Stage 1 — a mechanism check on a toy system confirms the claimed mathematical phenomenon reproduces; Stage 2 — a benchmark implements the idea on a real (small) neural network task and compares it against a tuned baseline over 8 paired seeds with a permutation test.
Artifacts
Artifacts unavailable.