REVIEW 3 major objections 6 minor 1 cited by
Unsupervised Physics-Informed Neural Network-based Nonlinear Observer design for autonomous systems using contraction analysis
T0 review · 3 major / 6 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read An unsupervised physics-informed neural network can learn contraction-based observer gains that guarantee global exponential convergence of the estimation error for autonomous nonlinear systems.
desk verdict The PINN idea is new, but the training target is infeasible for the paper's own examples, and the loss function's sign error makes the central claim unsupported. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the symmetric matrix $D(\hat{x},y)$ from the contraction condition (7), whose negative semi-definiteness is enforced via Sylvester's criterion using its leading principal minors. The loss function $L = \mu_1 L_{\text{MPDI}} + \mu_2 L_{\text{BC}}$ penalizes the signed principal minors (positive for odd indices, negative for even indices) and the boundary condition $\hat{k}_\theta(\hat{x}, h(\hat{x}))=0$, with collocation points sampled uniformly over the domain of interest. Automatic differentiation provides the Jacobian $\partial \hat{k}_\theta / \partial \hat{x}$ needed to form $D$.
What would settle it
Take the trained gain and evaluate the largest eigenvalue of $D(\hat{x},y) = \operatorname{He}\{\frac{\partial f}{\partial \hat{x}} + \frac{\partial \hat{k}_\theta}{\partial \hat{x}}\} + 2\lambda I$ on a dense grid over the domain $\mathcal{X}\times\mathcal{Y}$; if any eigenvalue is positive at a grid point, the MPDI (7) is violated and the claimed exponential convergence guarantee is not applicable to that gain.
Extended reading notes
Core claim
The paper's central claim is that the loss function (6)-(8), built from the leading principal minors of the matrix $D(\hat{x},y)=\operatorname{He}\{\frac{\partial f}{\partial \hat{x}} + \frac{\partial \hat{k}_\theta}{\partial \hat{x}}\} + 2\lambda I$, enforces the contraction MPDI (7) on the training domain, so the learned gain $\hat{k}_\theta$ yields a globally exponentially convergent observer for the autonomous nonlinear system (1). Theorem 3 further states that under a Lipschitz gain and bounded noise and approximation error, the estimation error satisfies the exponential input-to-state bound (14) whenever the contraction rate $\lambda > 2$. The claim is supported by simulation on the Van der Pol and reverse Duffing oscillators, where increasing $\lambda$ reduces the mean squared estimation error under 15% measurement noise.
Load-bearing premise
The whole guarantee rests on the unverified Assumption 3(A2) that the trained network's gain is uniformly close to some ideal gain that satisfies the strict contraction inequality with $P=I$ and $\lambda>2$; if no such ideal gain exists or the training is not $\varepsilon$-close to it, the exponential bounds in Theorem 3 do not follow.
Editorial extensions
If this is right
- The approach removes the need to solve the MPDI analytically; any smooth autonomous system with a known Jacobian can be plugged into the training loop directly.
- The learned observer is robust to bounded measurement noise and network approximation error, with the steady-state error bound shrinking as the contraction rate $\lambda$ increases.
- Because the loss is unsupervised (no state trajectories required), the method needs only collocation points, not data, to design an observer.
- The exponential input-to-state bound (14) gives a quantitative guarantee that the estimation error remains within a neighborhood determined by noise level and approximation error.
- Higher $\lambda$ improves noise rejection but makes it harder to satisfy the MPDI, revealing a design trade-off in choosing the contraction rate.
Reading between the lines
- A natural extension is to verify post-training whether the learned gain actually satisfies the MPDI on a dense grid; if not, the theoretical guarantee is void even if simulations look good. This check is not performed in the paper.
- The method could be extended to non-autonomous systems by making the gain depend explicitly on time, though the paper leaves this as future work.
- The choice $P=I$ simplifies the MPDI, but allowing a non-constant metric $P(x)$ could widen the class of acceptable systems; the same PINN framework might be adapted to learn both $P$ and $k$ together.
- The ISS bound depends on the Lipschitz constant $L$ of the gain and the bounds of noise and error; a post-training computation of these constants could turn the qualitative guarantee into a computable error certificate.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes an unsupervised physics-informed neural network (PINN) approach for designing the correction term of a nonlinear observer for autonomous systems. The observer gain is learned by minimizing a loss that penalizes violation of the contraction matrix partial differential inequality (MPDI) from Theorem 1, with the metric fixed to the identity, plus a boundary-condition loss. The authors also state a robustness result (Theorem 3) bounding the estimation error under measurement noise and neural-network approximation error, and they report simulations on a Van der Pol oscillator and a reverse Duffing oscillator. The central claim is that the learned gain enforces the contraction conditions and therefore yields a globally exponentially convergent observer whose robustness is quantified by Theorem 3.
Significance. The problem is relevant: numerically solving the contraction-based observer MPDI is genuinely challenging, and a learning-based unsupervised formulation is a reasonable idea to explore. I credit the authors for avoiding dependence on labeled trajectory data, for using automatic differentiation to compute the required Jacobians, and for attempting an explicit robustness analysis with noise and approximation error. However, the formal core of the paper is not sound as written. The loss function does not correctly encode negative semidefiniteness, and, more seriously, for the two numerical examples there exists no C1 gain satisfying the MPDI with P=I, so the assumptions underlying Theorem 3 are vacuous in the very systems used for validation. These are load-bearing issues that invalidate the main claims. The paper would need substantial reworking—correcting the loss, choosing a feasible metric and/or examples, and repairing Theorem 3—before its contribution could be credited.
major comments (3)
- [Section IV.B, Eq. (8)] The loss for even principal minors has the wrong sign. Since (7) requires D to be negative semidefinite, every even principal minor must be nonnegative, so the penalty should grow when Δ_i is negative. Writing l_i = min(0, Δ_i) for even i makes the loss decrease as Δ_i becomes more negative, thereby actively encouraging violation of the contraction condition. Independently, Theorem 2 as stated is false: for symmetric matrices, alternation of the leading principal minors is necessary but not sufficient for negative semidefiniteness; for example, diag(0,1) has leading minors Δ1=0 and Δ2=0, satisfying (−1)^1Δ1≥0 and (−1)^2Δ2≥0, yet the matrix is not negative semidefinite. A correct criterion requires all principal minors, not only the leading ones. Thus minimizing (6) does not enforce (7).
- [Sections IV.B and VI, Eqs. (7), (19), (20)] With P=I and h(x)=x1, no C1 gain satisfying (7) exists for either test system. The boundary condition k(x,h(x))=0 in (4), together with C1 regularity, forces ∂k₂/∂x̂₂=0 at every point of the form (x̂₁,x̂₂,x̂₁). For the reverse Duffing oscillator (20), ∂f₂/∂x₂=0, so the (2,2) entry of D in (7) is 2λ>0. For the Van der Pol oscillator (19), at x₁=0 that entry is 1+2λ>0. A symmetric matrix with a strictly positive diagonal entry cannot be negative semidefinite, and the condition (7) is required to hold for all (x̂,y), including y=h(x̂). Therefore Assumption 3(A2) is vacuous in the examples, and the learned observer cannot inherit the guarantees of Theorem 1 or Theorem 3. The boundary loss (9) does not repair the obstruction, because it penalizes the value of k at the diagonal, not its derivative along x̂₂.
- [Section V, Theorem 3 and Eq. (14)] The bound in (14) does not follow from the proof. From (17) and the claimed inequality (18), the differential inequality is ˙V ≤ −(λ−2)V + (1/2)∥ε∥² + (L²/2)∥v∥². The solution of this inequality yields an estimation bound with decay rate (λ−2)/2 and asymptotic gain proportional to (∥ε̄∥² + L²∥v̄∥²)^{1/2}/(λ−2)^{1/2}, not the stated e^{−ηt} factor and the term (ε̄+Lv̄)/(2√η). The step from (17) to (18) is also insufficiently justified: the contraction condition gives, at best, e^T[f(x)−f(x̂)−k(x̂,y)] ≤ −λ∥e∥², and the additional ∥e∥² term from Young's inequality changes the effective rate; the citation to [31, Theorem 4.3, p. 231] is too terse. Since Theorem 3 is presented as the formal robustness guarantee, this gap is significant.
minor comments (6)
- [Throughout] The manuscript inconsistently uses both “MPDI” and “MDPI” for the matrix partial differential inequality; please standardize.
- [Section IV.B] “Principle minors” should be “principal minors” in the text describing Eq. (8).
- [Section IV.C] The training algorithm is said to be “provided upon request,” but for a methods paper the pseudocode should be fully self-contained, including batch size and the exact stopping criteria, or the code should be released.
- [Section VI] Figure 3 is captioned “Inverse Duffing” while the text and Eq. (20) call the system the “reverse Duffing” oscillator; the nomenclature should be consistent.
- [Section VI] The dataset region is written as “(X , Y) = ([ −1, 1]², [−1, 1])”; this notation is ambiguous and should be clarified, since X is a subset of R² and Y is a subset of R.
- [Section V, proof of Theorem 3] The “identity” from [31] that converts the path term into the contraction-rate term should be stated explicitly, including its assumptions, because the proof depends on it.
Circularity Check
No significant circularity: the contraction loss enforces an imported sufficient condition, and the stability/ISS claims follow from external contraction-observer theorems, not from fitted outputs.
full rationale
The derivation chain is not circular. The learned gain k̂θ is trained by minimizing a physics loss (6)-(9) that directly penalizes violation of the matrix partial differential inequality (7), i.e., the contraction condition from Theorem 1 with P = I. This is a standard PDE-constrained learning design: the learned object is constructed to satisfy a sufficient condition, and the exponential-convergence claim (Section IV and Theorem 3) is a consequence of that enforced inequality plus the external theorem [31], not an output inferred from fitted data. The parameter λ is a user-chosen design constant, and Table II is an empirical consistency check of Theorem 3, not a fit used to produce the bound (14). The paper does cite works by the same authors ([7], [26], [27]), but only as related work on algebraic estimators and KKL observer design; none of those citations carries the central existence or stability argument. The robustness proof invokes the 'identity in [31, see Theorem 4.3 p. 231]', an external survey result, so the paper is not relying on a self-citation chain. The known concerns—that the leading-principal-minor criterion in Theorem 2 is not by itself sufficient for negative semidefiniteness, and that Assumption 3(A2)'s existence of an ideal gain with P = I may fail for the example systems—are correctness/falsifiability issues, not circular reductions: they question whether the imported sufficient condition is satisfied, not whether the paper's derivation re-labels its own inputs.
Assumptions & free parameters
free parameters (3)
- lambda (contraction rate) =
2.5, 4, 5
- loss weights mu1, mu2, rho1, rho2 =
VdP: mu1=1e-3, mu2=1, rho1=1, rho2=1e-1; Duffing: mu1=1, mu2=1, rho1=1, rho2=1
- Training region X, Y =
VdP: X=[-2,2]x[-3,3], Y=[-2,2]; Duffing: X=[-1,1]^2, Y=[-1,1]
assumptions (6)
- standard math Theorem 1 (observer contraction condition)
- standard math Theorem 2 (Sylvester's criterion) applied to negative semidefiniteness
- domain assumption Assumption 1: forward invariance on compact sets X, Y
- domain assumption Assumption 2: learned gain is Lipschitz in its second argument
- domain assumption Assumption 3: bounded noise and bounded approximation error
- ad hoc to paper Existence of an ideal gain k satisfying (4) with P=I and lambda>2 for the considered systems
Cite this review
Pith. "Pith review of Unsupervised Physics-Informed Neural Network-based Nonlinear Observer design for autonomous systems using contraction analysis." pith.science (2026). https://pith.science/paper/QFDJJVQM
@misc{pith2026241109237,
author = {Pith},
title = {Pith review of: Unsupervised Physics-Informed Neural Network-based Nonlinear Observer design for autonomous systems using contraction analysis},
year = {2026},
howpublished = {\url{https://pith.science/paper/QFDJJVQM}},
note = {Machine review of arXiv:2411.09237}
}
read the original abstract
Contraction analysis offers, through elegant mathematical developments, a unified way of designing observers for a general class of nonlinear systems, where the observer correction term is obtained by solving an infinite dimensional inequality that guarantees global exponential convergence. However, solving the matrix partial differential inequality involved in contraction analysis design is both analytically and numerically challenging and represents a long-lasting challenge that prevented its wide use. Therefore, the present paper proposes a novel approach that relies on an unsupervised Physics Informed Neural Network (PINN) to design the observer's correction term by enforcing the partial differential inequality in the loss function. The performance of the proposed PINN-based nonlinear observer is assessed in numerical simulation as well as its robustness to measurement noise and neural network approximation error.
Figures
Forward citations
Cited by 1 Pith paper
-
Designing Robust Software Sensors for Nonlinear Systems via Neural Networks and Adaptive Sliding Mode Control
A neural network that learns time-varying observer gains, combined with adaptive sliding mode control, estimates states of general nonlinear systems directly from sensor measurements.
Reference graph
Works this paper leans on
-
[1]
D. G. Luenberger. Observing the state of a linear system. IEEE Transactions on Military Electronics , 8(2):74–80, 1964
work page 1964
-
[2]
R. E. Kalman. A New Approach to Linear Filtering and Prediction Problems. Journal of Basic Engineering , 82(1):35–45, 03 1960
work page 1960
-
[3]
Hassan K. Khalil. High-gain observers in nonlinear feedback control. In 2008 International Conference on Control, Automation and Systems, pages xlvii–lvii, 2008
work page 2008
-
[4]
Invariant manifold based reduced-order observer design for nonlinear systems
Dimitrios Karagiannis, Daniele Carnevale, and Alessandro Astolfi. Invariant manifold based reduced-order observer design for nonlinear systems. IEEE Transactions on Automatic Control, 53(11):2602–2614, 2008
work page 2008
-
[5]
Linearization by output injection and nonlinear observers
Arthur J Krener and Alberto Isidori. Linearization by output injection and nonlinear observers. Systems & Control Letters, 3(1):47–52, 1983
work page 1983
- [6]
-
[7]
Yasmine Marani, Ibrahima N’Doye, and Taous Meriem Laleg-Kirati. Non-asymptotic neural network-based state and disturbance estimation for a class of nonlinear systems using modulating functions. In 2023 American Control Conference, 2023
work page 2023
-
[8]
N. Kazantzis and C. Kravaris. Nonlinear observer design using Lyapunov’s auxiliary theorem. In Proceedings of the 36th IEEE Conference on Decision and Control , volume 5, pages 4802–4807 vol.5, 1997
work page 1997
Show all 38 references
-
[9]
Anderson and J.B
B.D.O. Anderson and J.B. Moore. Optimal Filtering. Information and system sciences series. Prentice-Hall, 1979
1979
-
[10]
Lohmiller and J.-J.E
W. Lohmiller and J.-J.E. Slotine. On metric observers for nonlinear systems. In Proceeding of the 1996 IEEE International Conference on Control Applications , pages 320–326, 1996
1996
-
[11]
Lohmiller and Slotine
W. Lohmiller and Slotine. On contraction analysis for non-linear systems. Automatica, 34(6):683–696, 1998
1998
-
[12]
Nonlinear dynamical control systems
Henk Nijmeijer and Arjan van der Schaft. Nonlinear dynamical control systems. Springer-Verlag, Berlin, Heidelberg, 1990
1990
-
[13]
A. Isidori. Nonlinear Control Systems . Communications and Control Engineering. Springer London, 1995
1995
-
[14]
H.K. Khalil. Nonlinear Systems . Pearson Education. Prentice Hall, 2002
2002
-
[15]
Contracting nonlinear observers: Convex opti- mization and learning from data
Ian R Manchester. Contracting nonlinear observers: Convex opti- mization and learning from data. In 2018 annual American Control Conference (ACC), pages 1873–1880. IEEE, 2018
2018
-
[16]
Lohmiller and J.-J.E
W. Lohmiller and J.-J.E. Slotine. Applications of metric observers for nonlinear systems. In Proceeding of the 1996 IEEE International Conference on Control Applications , pages 367–372, 1996
1996
-
[17]
Lohmiller and J.-J.E
W. Lohmiller and J.-J.E. Slotine. Simple observers for hamiltonian systems. In Proceedings of the 1997 American Control Conference (Cat. No.97CH36041), volume 5, pages 2748–2753 vol.5, 1997
1997
-
[18]
Convergence of nonlinear observers on Rn with a riemannian metric (part i)
Ricardo G Sanfelice and Laurent Praly. Convergence of nonlinear observers on Rn with a riemannian metric (part i). IEEE Transactions on Automatic Control , 57(7):1709–1722, 2011
2011
-
[19]
Observer de- sign for stochastic nonlinear systems via contraction-based incremental stability
Ashwin P Dani, Soon-Jo Chung, and Seth Hutchinson. Observer de- sign for stochastic nonlinear systems via contraction-based incremental stability. IEEE Transactions on Automatic Control , 60(3):700–714, 2014
2014
-
[20]
Review on contraction analysis and computation of contraction metrics
Peter Giesl, Sigurdur Hafstein, and Christoph Kawan. Review on contraction analysis and computation of contraction metrics. Journal of Computational Dynamics , 10(1):1–47, 2023
2023
-
[21]
Physics- informed neural networks: A deep learning framework for solving forward and inverse problems involving nonlinear partial differential equations
Maziar Raissi, Paris Perdikaris, and George E Karniadakis. Physics- informed neural networks: A deep learning framework for solving forward and inverse problems involving nonlinear partial differential equations. Journal of Computational physics , 378:686–707, 2019
2019
-
[22]
Parameter es- timation and modeling of nonlinear dynamical systems based on runge–kutta physics-informed neural network
Weida Zhai, Dongwang Tao, and Yuequan Bao. Parameter es- timation and modeling of nonlinear dynamical systems based on runge–kutta physics-informed neural network. Nonlinear Dynamics , 111(22):21117–21130, 2023
2023
-
[23]
Nonlinear discrete-time observers with physics-informed neural networks
Hector Vargas Alvarez, Gianluca Fabiani, Nikolaos Kazantzis, Ioan- nis G Kevrekidis, and Constantinos Siettos. Nonlinear discrete-time observers with physics-informed neural networks. Chaos, Solitons & Fractals, 186:115215, 2024
2024
-
[24]
Learning-based design of luenberger observers for autonomous nonlinear systems
Muhammad Umar B Niazi, John Cao, Xudong Sun, Amritam Das, and Karl Henrik Johansson. Learning-based design of luenberger observers for autonomous nonlinear systems. In 2023 American Control Conference (ACC), pages 3048–3055. IEEE, 2023
2023
-
[25]
Deep learning-based luenberger observer design for discrete-time nonlinear systems
Johan Peralez and Madiha Nadri. Deep learning-based luenberger observer design for discrete-time nonlinear systems. In 2021 60th IEEE Conference on Decision and Control (CDC) , pages 4370–4375. IEEE, 2021
2021
-
[26]
Deep-learning based design of cascade observers for discrete-time nonlinear systems with output delay
Yasmine Marani, Ibrahima N’Doye, and Taous Meriem Laleg-Kirati. Deep-learning based design of cascade observers for discrete-time nonlinear systems with output delay. IFAC-PapersOnLine, 56(2):9869– 9874, 2023. 22nd IFAC World Congress
2023
-
[27]
Deep-learning based kkl chain observer for discrete-time nonlinear systems with time-varying output delay
Yasmine Marani, Ibrahima N’Doye, and Taous Meriem Laleg-Kirati. Deep-learning based kkl chain observer for discrete-time nonlinear systems with time-varying output delay. Automatica, 171:111955, 2025
2025
-
[28]
Physics-informed neural nets for control of dynamical systems
Eric Aislan Antonelo, Eduardo Camponogara, Laio Oriel Seman, Jean Panaioti Jordanou, Eduardo Rehbein de Souza, and Jomi Fred H¨ubner. Physics-informed neural nets for control of dynamical systems. Neurocomputing, 579:127419, 2024
2024
-
[29]
Physics-informed neural networks with skip connections for modeling and control of gas-lifted oil wells
Jonas Ekeland Kittelsen, Eric Aislan Antonelo, Eduardo Camponog- ara, and Lars Struen Imsland. Physics-informed neural networks with skip connections for modeling and control of gas-lifted oil wells. Applied Soft Computing , 158:111603, 2024
2024
-
[30]
Stability of Motion
Wolfgang Hahn. Stability of Motion . Springer, 1967
1967
-
[31]
Observer design for continuous-time dynamical systems
Pauline Bernard, Vincent Andrieu, and Daniele Astolfi. Observer design for continuous-time dynamical systems. Annual Reviews in Control, 53:224–248, 2022
2022
-
[32]
Positive definite matrices and sylvester’s criterion
George T Gilbert. Positive definite matrices and sylvester’s criterion. The American Mathematical Monthly , 98(1):44–46, 1991
1991
-
[33]
Adam: A method for stochastic optimization
Diederik P Kingma. Adam: A method for stochastic optimization. arXiv preprint arXiv:1412.6980 , 2014
2014 arXiv
-
[34]
A limited memory algorithm for bound constrained optimization
Richard H Byrd, Peihuang Lu, Jorge Nocedal, and Ciyou Zhu. A limited memory algorithm for bound constrained optimization. SIAM Journal on scientific computing , 16(5):1190–1208, 1995
1995
-
[35]
Challenges in training pinns: A loss landscape perspective
Pratik Rathore, Weimu Lei, Zachary Frangella, Lu Lu, and Madeleine Udell. Challenges in training pinns: A loss landscape perspective. arXiv preprint arXiv:2402.01868 , 2024
2024 arXiv
-
[36]
Physics-informed neural networks for solving forward and inverse flow problems via the boltzmann-bgk formulation
Qin Lou, Xuhui Meng, and George Em Karniadakis. Physics-informed neural networks for solving forward and inverse flow problems via the boltzmann-bgk formulation. Journal of Computational Physics , 447:110676, 2021
2021
-
[37]
A mixed pressure-velocity formulation to model flow in heterogeneous porous media with physics-informed neural networks
Franc ¸ois Lehmann, Marwan Fahs, Ali Alhubail, and Hussein Hoteit. A mixed pressure-velocity formulation to model flow in heterogeneous porous media with physics-informed neural networks. Advances in Water Resources, 181:104564, 2023
2023
-
[38]
Improved training of physics-informed neural networks with model ensembles
Katsiaryna Haitsiukevich and Alexander Ilin. Improved training of physics-informed neural networks with model ensembles. In 2023 International Joint Conference on Neural Networks (IJCNN) , pages 1–8. IEEE, 2023
2023
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.