REVIEW 1 major objections 4 minor 43 references
Data-driven stabilization of continuous-time systems with noisy input-output data
T0 review · 1 major / 4 minor · reviewed 2026-08-03 · deepseek-v4-flash
Pith's one-line read This paper establishes a necessary and sufficient LMI condition under which noisy continuous-time input-output data are informative for quadratic stabilization, and proves that any feasible solution yields a controller that stabilizes every
desk verdict The main LMI theorem is sound and genuinely new; referee it, and ask for code/data and a proof of the delegated set equality. 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 paper's workhorse is the synthesis operator T_ℓ, which maps a test function to a weighted integral of the ℓ-th derivative of the measured trajectory; stacking these operators gives a finite-rank data Hankel operator H and a data-and-noise matrix Π. Lemma 2.1 states that R is consistent with the noisy data exactly when (I_p, R)ᵀΠ(I_p,R) ≥ 0. The strict matrix S-lemma (Lemma 4.2) then converts the requirement 'the stability LMI holds for all such R' into the single LMI (50), with the surjectivity of H ensuring the negative-definiteness and rank conditions the S-lemma needs. The identification result uses the same H: unique recovery of the system from noise-free data is equivalent to H bein
What would settle it
Take a dataset satisfying Assumptions 4.3 and 4.4, solve LMI (50) to get C = DΦ⁻¹, then pick any R consistent with the data (i.e., with (I_p,R)ᵀΠ(I_p,R) ≥ 0) and simulate the closed-loop response from a nonzero initial condition; a bounded or non-decaying trajectory would show the controller does not stabilize all consistent systems. Conversely, if one can exhibit informative data (a stabilizing controller exists for every consistent R) for which LMI (50) is infeasible, the necessity direction fails.
Extended reading notes
Core claim
The central claim is Theorem 4.6: under the assumptions that the data were generated by some system in the noise class (Assumption 4.3) and that the data-embedded operator H is surjective (Assumption 4.4), the data are informative for quadratic stabilization if and only if the LMI (50) is feasible. The matrix N in that LMI is built from products of the synthesis operators and their adjoints, with explicit integral formulas given in Proposition 3.3. Moreover, if the LMI holds with matrices Φ and D, then the controller C = DΦ⁻¹ stabilizes every system in the consistency set Σ_{D,Θ}. The proof hinges on Lemma 2.1, which identifies Σ_{D,Θ} with the set of matrices R satisfying a quadratic inequa
Load-bearing premise
The entire equivalence rests on the data-embedded operator H being surjective, meaning the measured trajectories are rich enough that HH* is invertible; if the data are not that rich, the LMI condition in Theorem 4.6 is not established.
Editorial extensions
If this is right
- A practitioner can decide from raw input-output data alone whether stabilization is possible, and obtain a controller from the same LMI, without state or output derivatives and without filtering or intermediate identification.
- Because the condition is necessary and sufficient, failure of the LMI certifies that the available data are not informative for quadratic stabilization under the assumed noise bound.
- The resulting controller is guaranteed to stabilize every system consistent with the data, so the design is robust to the particular noise realization as long as it lies in the specified class.
- In the noise-free limit, the same data operator H determines whether the system coefficients are uniquely identifiable, with an explicit reconstruction formula.
- The operator-based embedding applies to the full continuous-time trajectories rather than sampled or filtered proxies, so no pre-processing parameters enter the stability condition.
Reading between the lines
- If H is not surjective, the paper leaves the informativity question open; in that regime the LMI may be conservative, and one could explore whether a weaker condition, perhaps after adding one more experiment, restores informativeness.
- The identification formula R_s = -Y_L H*(HH*)^-1 suggests a derivative-free parameter estimator that needs no basis choice; testing its numerical conditioning against spline-based alternatives would be a natural next step.
- The noise class bounds total energy of the embedded noise rather than pointwise size, so the results likely extend to disturbance energies across experiments; it would be interesting to compare with noise models that bound sample paths pointwise.
- A direct application of the same S-lemma logic could yield informativity certificates for other objectives, such as H2 or H∞ performance, since the stability certificate is already a quadratic differential form.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies data-driven stabilization of continuous-time LTI systems in AR form from noisy input-output trajectories. It embeds the data into synthesis operators (Section 2), characterizes the set of data-consistent systems (Lemma 2.1), derives computable integral formulas for products of synthesis operators and adjoints (Section 3), and proves a necessary and sufficient LMI condition for informativity for quadratic stabilization (Theorem 4.6), under a data-richness assumption on the operator H. A noise-free identification result is also given (Section 5), and a numerical example illustrates the approach (Section 6).
Significance. If Theorem 4.6 is correct, this is the first necessary and sufficient LMI certificate for continuous-time output-feedback stabilization from raw noisy input-output data in the behavioral/AR framework. The synthesis-operator embedding avoids filtering/sampling and yields finite-dimensional matrix conditions. The paper is technically careful: the algebra in Lemmas 2.1, 3.1-3.2, Proposition 3.3, Lemma 4.5, and Theorem 4.6 is coherent, and the use of the strict matrix S-lemma is clearly identified. The noise-free identification result also gives a clean interpretation of the surjectivity assumption on H. The main weaknesses are local: a proof typo in the converse of Theorem 4.6 and the simulation's use of white noise outside the stated L^2/H^L framework.
major comments (1)
- [Theorem 4.6, proof of converse] The converse proof defines C:=DPhi and Psi:=Phi^{-1}. With A=[-J_{q(L-1)}Phi; D; 0_{p x qL}], the (1,1) block of the LMI (50) is A^T+A = Phi(-J)^T+(-J)Phi + D^T+D. Substituting D=C Phi^{-1} gives Phi(-J)^T+(-J)Phi + Phi^{-1} C^T + C Phi^{-1}, whereas the corresponding block of M in (51) with Psi^{-1}=Phi is Phi(-J)^T+(-J)Phi + Phi C^T + C Phi. These are not equal. The intended definition is C:=D Phi^{-1}, as already stated in the theorem's controller formula; with that change, A=[-J Phi; C Phi; 0] and A^T+A equals the (1,1)-block of M. Please correct this typo in the proof.
minor comments (4)
- [Section 4.2, Eq. (54)] The set equality (54) is delegated to [31, Thm 3.4]. It follows directly from the block structure of N: the (1,1) block of [J_p; Z]^T Pi [J_p; Z] forces the first qL-p columns of any Z in Z_{qL,qL}(N) to vanish because N22=-HH*<0, after which Z=R^T J_p and Lemma 2.1 applies. A two-line derivation would make the main theorem self-contained.
- [Section 6, Example] The example applies zero-mean Gaussian white noise with E[v(t)v(s)^T]=delta(t-s)10^{-4}I_2. Such a process is almost surely not an element of L^2([0,tau];R^p), and the resulting output is not in H^L([0,tau];R^p) for L=2. This violates the standing assumptions v_k in L^2 and y_k in H^L from Section 2.1, so the numerical verification of Assumption 4.3 via Proposition 3.3 is not well-defined for the reported noise. Please replace the white noise by a band-limited or colored L^2 approximation, or state explicitly how the white noise is discretized into an L^2 representative.
- [Notation] The noise class is introduced as Delta_{tau,Theta} but written as Delta_Theta[0,tau] in Section 6. Please use consistent notation throughout.
- [Proposition 3.3] The dimensions of Lambda_ell should be stated explicitly (Lambda_ell in R^{L x n}) to avoid confusion, since tilde f_ell takes values in R^n.
Circularity Check
No significant circularity: the central LMI equivalence is derived from data and external S-lemma results; only minor non-load-bearing self-citations appear.
full rationale
Theorem 4.6 is not circular. The data-consistent set Σ_{D,Θ} is characterized in Lemma 2.1 directly from the defining equations (5) and the noise class, via the identity (I_p R^⊤)^⊤ Π (I_p R^⊤) = Θ − Σ_k (RH_k + Y_{L,k})(RH_k + Y_{L,k})^*. Lemma 4.5 verifies the strict S-lemma hypotheses using only Assumptions 4.3/4.4 and Lemma 2.1. The key set equality (54) is delegated to the external [31, Thm 3.4], not authored by Wakaiki, and the S-lemma itself is [31, Thm 4.10]; these are independent supports. The controller C = DΦ^{-1} in Theorem 4.6 is constructed from the LMI solution, not fitted to a target outcome, and the theorem is an iff statement proven in both directions. The noise-free identification result Proposition 5.3 follows from RH + Y_L = 0 and the surjectivity of H; no target result is assumed. The only self-citations ([36], [37]) concern methodology and a range approximation lemma in Lemma 5.2; they are not load-bearing for the central stabilization theorem, and the cited prior results are consistent with standard operator arguments. The simulation's use of white noise outside the stated L^2 noise class and the reliance on Assumption 4.4 are correctness/reproducibility concerns, not circularity.
Assumptions & free parameters
free parameters (3)
- Noise-intensity bound Theta (example) =
10^-6 I_2
- Trajectory acceptance threshold |theta| <= 0.1 =
0.1 rad
- Data horizon and trajectory count (example) =
tau = 0.5, K = 25
assumptions (8)
- standard math Strict matrix S-lemma (Lemma 4.2)
- standard math Equivalence (54) between R in Sigma and a Z-space of N
- standard math Quadratic differential form Lyapunov criterion (Theorem A.1)
- domain assumption The unknown system is exactly a continuous-time AR system of known order L (eq. (2))
- domain assumption Noise v_k lies in L^2 and is bounded by known Theta in synthesis-operator norm (Delta_{tau,Theta})
- domain assumption Assumption 4.3: some true system R_s and noise sequence in Delta_{tau,Theta} are consistent with the data
- domain assumption Assumption 4.4: the data-embedded operator H is surjective
- standard math Density of B-splines in H0^L (Lemma 5.1)
Cite this review
Pith. "Pith review of Data-driven stabilization of continuous-time systems with noisy input-output data." pith.science (2026). https://pith.science/paper/MAISCJBE
@misc{pith2026260202992,
author = {Pith},
title = {Pith review of: Data-driven stabilization of continuous-time systems with noisy input-output data},
year = {2026},
howpublished = {\url{https://pith.science/paper/MAISCJBE}},
note = {Machine review of arXiv:2602.02992}
}
read the original abstract
We study data-driven stabilization of continuous-time systems in autoregressive form when only noisy input-output data are available. First, we provide an operator-based characterization of the set of systems consistent with the data. Next, combining this characterization with behavioral theory, we establish a necessary and sufficient condition for the noisy data to be informative for quadratic stabilization. This condition is formulated in terms of linear matrix inequalities, whose solutions yield a stabilizing controller. Finally, we characterize data informativity for system identification in the noise-free setting.
Reference graph
Works this paper leans on
-
[33]
H. J. van Waarde, J. Eising, M. K. Camlibel, and H. L. Trentelman. A behavioral approach to data-driven control with noisy input-output data.IEEE Trans. Automat. Control, 69:813–827, 2023
2023
-
[36]
M. Wakaiki. Data-driven control of continuous-time systems: A synthesis-operator approach. arXiv:2511.21041, 2025
arXiv 2025
-
[1]
D. S. Bernstein.Scalar, Vector, and Matrix Mathematics: Theory, Facts, and Formulas. Princeton Univ. Press, Princeton, NJ, 2018
2018
-
[2]
Bisoffi, C
A. Bisoffi, C. De Persis, and P. Tesi. Data-driven control via Petersen’s lemma.Automatica, 145, Art. no. 110537, 14 pp., 2022
2022
-
[3]
Bisoffi, L
A. Bisoffi, L. Li, C. De Persis, and N. Monshizadeh. Controller synthesis for input-state data with measurement errors.IEEE Control Syst. Lett., 8:1571–1576, 2024
2024
- [4]
-
[5]
Bosso, M
A. Bosso, M. Borghesi, A. Iannelli, G. Notarstefano, and A. R. Teel. Derivative-free data-driven control of continuous-time linear time-invariant systems.Eur. J. Control, Art. no. 105965, 8 pp., 2025
2025
- [6]
Show all 43 references
-
[7]
Christensen.An Introduction to Frames and Riesz Bases
O. Christensen.An Introduction to Frames and Riesz Bases. Birkh¨ auser, Cham, 2nd edition, 2016
2016
-
[8]
Cotroneo and J
T. Cotroneo and J. C. Willems. An LMI condition for nonnegativity of a quadratic differential form along a behavior. In Proc. MTNS 2000,6 pp., 2000
2000
-
[9]
De Persis, R
C. De Persis, R. Postoyan, and P. Tesi. Event-triggered control from data.IEEE Trans. Automat. Control, 69:3780–3795, 2024
2024
-
[10]
De Persis and P
C. De Persis and P. Tesi. Formulas for data-driven control: Stabilization, optimality, and robustness.IEEE Trans. Automat. Control, 65:909–924, 2020
2020
-
[11]
Eising and J
J. Eising and J. Cort´ es. When sampling works in data-driven control: Informativity for stabilization in continuous time.IEEE Trans. Automat. Control, 70:565–572, 2025
2025
-
[12]
H. Gao, A. Bosso, L. Wang, D. Saussi´ e, and B. Yi. Input-output data-driven stabilization of continuous-time linear MIMO systems. arXiv:2511.06524, 2025
2025
-
[13]
Kaminaga and H
T. Kaminaga and H. Sasahara. Data informativity for quadratic stabilization under data perturbation. arXiv:2505.01641, 2025
2025 arXiv
-
[14]
Kaneko and T
O. Kaneko and T. Fujii. Algebraic stability criterion for high order differential equations in a behavioral framework. In Proc. 37th IEEE CDC, pages 132–137, 1998
1998
-
[15]
L¨ ofberg
J. L¨ ofberg. YALMIP: A toolbox for modeling and optimization in MATLAB. InProc. IEEE Int. Conf. Robot. Autom., pages 284–289, 2004
2004
-
[16]
V. G. Lopez and M. A. M¨ uller. On a continuous-time version of Willems’ lemma. InProc. 61st IEEE CDC, pages 2759–2764, 2022
2022
-
[17]
V. G. Lopez and M. A. M¨ uller. An input-output continuous-time version of Willems’ lemma.IEEE Control Syst. Lett., 8:916– 921, 2024. 18 MASASHI W AKAIKI
2024
-
[18]
Version 11.0.29, 2025
MOSEK ApS.MOSEK API for MATLAB manual. Version 11.0.29, 2025
2025
-
[19]
Y. Ohta. Data informativity of continuous-time systems by sampling using linear functionals. InProc. MTNS 2024,6 pp., 2024
2024
-
[20]
Ohta and P
Y. Ohta and P. Rapisarda. A sampling linear functional framework for data-driven analysis and control of continuous-time systems. InProc. 63rd IEEE CDC, pages 357–362, 2024
2024
-
[21]
Possieri
C. Possieri. A derivative free framework for data-driven stabilization of continuous-time linear systems.IEEE Control Syst. Lett., 9:1219–1224, 2025
2025
-
[22]
Possieri
C. Possieri. Derivative free data-driven stabilization of continuous-timelinear systems from input-output data. To appear in IEEE Control Syst. Lett., 2026
2026
-
[23]
H. A. Preisig and D. W. T. Rippin. Theory and application of the modulating function method–I. Review and theory of the method and theory of the spline-type modulating functions.Comput. Chem. Eng., 17:1–16, 1993
1993
-
[24]
fundamental lemma
P. Rapisarda, M. K. C ¸ amlibel, and H. J. van Waarde. A “fundamental lemma” for continuous-time systems, with applications to data-driven simulation.Syst. Control Lett., 179, Art. no. 105603, 9 pp., 2023
2023
-
[25]
Rapisarda, H
P. Rapisarda, H. J. van Waarde, and M. K. C ¸ amlibel. Orthogonal polynomial bases for data-driven analysis and control of continuous-time systems.IEEE Trans. Automat. Control, 69:4307–4319, 2024
2024
-
[26]
Schmitz, T
P. Schmitz, T. Faulwasser, P. Rapisarda, and K. Worthmann. A continuous-time fundamental lemma and its application in data-driven optimal control.Syst. Control Lett., 194, Art. no. 105950, 11 pp., 2024
2024
-
[27]
Schumaker.Spline Functions: Basic Theory
L. Schumaker.Spline Functions: Basic Theory. Cambridge Univ. Press, Cambridge, UK, 3rd edition, 2007
2007
-
[28]
Shinbrot
M. Shinbrot. On the analysis of linear and nonlinear dynamical systems from transient-response data. National Advisory Commitee for Aeronautics, Technical Note 3288, 1954
1954
-
[29]
T. R. Steentjes, M. Lazar, and P. M. Van den Hof. On data-driven control: Informativity of noisy input-output data with cross-covariance bounds.IEEE Control Syst. Lett., 6:2192–2197, 2022
2022
-
[30]
Tucsnak and G
M. Tucsnak and G. Weiss.Observation and Control of Operator Semigroups. Birkh¨ auser, Basel, 2009
2009
-
[31]
H. J. van Waarde, M. K. Camlibel, J. Eising, and H. L. Trentelman. Quadratic matrix inequalities with applications to data-based control.SIAM J. Control Optim., 61:2251–2281, 2023
2023
-
[32]
H. J. van Waarde, M. K. Camlibel, and M. Mesbahi. From noisy data to feedback controllers: Nonconservative design via a matrix S-lemma.IEEE Trans. Automat. Control, 67:162–175, 2022
2022
-
[34]
H. J. van Waarde, J. Eising, M. K. Camlibel, and H. L. Trentelman. The informativity approach: To data-driven analysis and control.IEEE Control Syst. Mag., 43:32–66, 2023
2023
-
[35]
H. J. van Waarde, J. Eising, H. L. Trentelman, and M. K. Camlibel. Data informativity: A new perspective on data-driven analysis and control.IEEE Trans. Automat. Control, 65:4753–4768, 2020
2020
-
[37]
M. Wakaiki. Data informativity for stabilization of discrete-time infinite-dimensional systems. arXiv:2511.16008, 2025
2025
-
[38]
J. Wang, S. Baldi, and H. J. van Waarde. Experiment design for continuous-time systems using generalized filtering. arXiv:2511.09386, 2025
2025 arXiv
-
[39]
J. C. Willems. From time series to linear system. Part I: Finite dimensional linear time-invariant systems.Automatica, 22:561– 580, 1986
1986
-
[40]
J. C. Willems. From time series to linear system. Part II: Exact modelling.Automatica, 22:675–694, 1986
1986
-
[41]
J. C. Willems. From time series to linear system. Part III: Approximate modelling.Automatica, 23:87–115, 1987
1987
-
[42]
J. C. Willems, P. Rapisarda, I. Markovsky, and B. L. De Moor. A note on persistency of excitation.Syst. Control Lett., 54:325–329, 2005
2005
-
[43]
J. C. Willems and H. L. Trentelman. On quadratic differential forms.SIAM J. Control Optim., 36:1703–1749, 1998. Graduate School of System Informatics, Kobe University, Nada, Kobe, Hyogo 657-8501, Japan Email address:wakaiki@ruby.kobe-u.ac.jp
1998
Reviewed August 3, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.