Pith. sign in

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 →

arxiv 2602.02992 v2 pith:MAISCJBE submitted 2026-02-03 math.OC cs.SYeess.SY

classification math.OCcs.SYeess.SY MSC 93D1593B3093C05
keywords data-drivencontrolcontinuous-timesystemsdatainformativityquadraticstabilizationlinearmatrixinequalitiessynthesisoperatorsbehavioralsystemidentification
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper asks when a finite set of noisy continuous-time input-output trajectories is enough to design a stabilizing output-feedback controller for an unknown linear system in autoregressive form, and answers with an exact matrix certificate. It embeds each trajectory into finite-rank synthesis operators, which lets the set of systems consistent with the data be written as a single quadratic matrix inequality. Combining that description with a strict matrix S-lemma and a Lyapunov-style stability condition, it proves that the data are informative for quadratic stabilization if and only if one linear matrix inequality (LMI) is feasible. A feasible solution directly gives a controller that stabilizes all systems consistent with the data. The same operator machinery also characterizes when noise-free data identify the system uniquely.

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.

Watch

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

Editorial extensions of the paper, not claims the author makes directly.

  • 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.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

1 major / 4 minor

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)
  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)
  1. [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.
  2. [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.
  3. [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.
  4. [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

0 steps flagged · score 2.0 of 10

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 3 free parameters · 8 assumptions · 0 invented entities

The central theorem carries no fitted constants; its inputs are the known noise bound Theta, the known AR order L, and the surjectivity/HH* invertibility assumption. The example adds hand-chosen simulation parameters. The mathematical machinery (synthesis operators, QDF criteria, strict S-lemma) is drawn from prior literature, with two load-bearing external results cited but not reproved.

free parameters (3)
  • Noise-intensity bound Theta (example) = 10^-6 I_2
    Selected by hand in Section 6; the theorem requires Theta known. The actual embedded noise energy computed from the simulated data is <= 7.1e-7 I_2, so the choice is consistent but not derived from first principles.
  • Trajectory acceptance threshold |theta| <= 0.1 = 0.1 rad
    Data selection in Section 6 to keep the linearized model valid. Not used in the theoretical result, but shapes the example dataset.
  • Data horizon and trajectory count (example) = tau = 0.5, K = 25
    Experiment design in Section 6; these are user choices, not fitted to the claim.
assumptions (8)
  • standard math Strict matrix S-lemma (Lemma 4.2)
    Cited from [31, Theorem 4.10] and applied in Theorem 4.6 to convert inclusion of data-consistent systems into existence of a scalar alpha.
  • standard math Equivalence (54) between R in Sigma and a Z-space of N
    Delegated to [31, Theorem 3.4]; load-bearing in the proof of Theorem 4.6 because it re-encodes the data-consistent set as Z_{qL,qL}(N).
  • standard math Quadratic differential form Lyapunov criterion (Theorem A.1)
    Cited from [43]; the basis of Lemma 4.1 and therefore of the definition of quadratic stabilization used in Definition 4.2.
  • domain assumption The unknown system is exactly a continuous-time AR system of known order L (eq. (2))
    The theory assumes no unmodeled dynamics and exact knowledge of L; this is stated in Section 2 but is a real modeling limitation.
  • domain assumption Noise v_k lies in L^2 and is bounded by known Theta in synthesis-operator norm (Delta_{tau,Theta})
    Defines the data-consistent set. The simulation uses Gaussian white noise, which is not literally in the L^2 class; the paper verifies the operator bound a posteriori.
  • domain assumption Assumption 4.3: some true system R_s and noise sequence in Delta_{tau,Theta} are consistent with the data
    Used to prove nonnegativity of the Schur complement in Lemma 4.5. If the chosen Theta is too small, this assumption fails and the theorem is inapplicable.
  • domain assumption Assumption 4.4: the data-embedded operator H is surjective
    Equivalent to HH* invertible; required for N22 < 0 and the strict matrix S-lemma. It is also the identifiability condition in Proposition 5.3.
  • standard math Density of B-splines in H0^L (Lemma 5.1)
    Used in Section 5 to pass from the infinite-dimensional operator H to finite-dimensional B-spline approximations; standard spline theory.

how reviews work

0 comments
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.

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

43 extracted references · 4 linked inside Pith

  1. [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

  2. [36]

    M. Wakaiki. Data-driven control of continuous-time systems: A synthesis-operator approach. arXiv:2511.21041, 2025

  3. [1]

    D. S. Bernstein.Scalar, Vector, and Matrix Mathematics: Theory, Facts, and Formulas. Princeton Univ. Press, Princeton, NJ, 2018

  4. [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

  5. [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

  6. [4]

    Bosso, M

    A. Bosso, M. Borghesi, A. Iannelli, G. Notarstefano, and A. R. Teel. Data-driven control of continuous-time LTI systems via non-minimal realizations. arXiv:2505.22505, 2025

  7. [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

  8. [6]

    Bosso, M

    A. Bosso, M. Borghesi, A. Iannelli, B. Yi, and G. Notarstefano. Data-driven stabilization of continuous-time LTI systems from noisy input-output data. arXiv:2511.11417, 2025

Show all 43 references
  1. [7]

    Christensen.An Introduction to Frames and Riesz Bases

    O. Christensen.An Introduction to Frames and Riesz Bases. Birkh¨ auser, Cham, 2nd edition, 2016

  2. [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

  3. [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

  4. [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

  5. [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

  6. [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

  7. [13]

    Kaminaga and H

    T. Kaminaga and H. Sasahara. Data informativity for quadratic stabilization under data perturbation. arXiv:2505.01641, 2025

  8. [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

  9. [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

  10. [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

  11. [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

  12. [18]

    Version 11.0.29, 2025

    MOSEK ApS.MOSEK API for MATLAB manual. Version 11.0.29, 2025

  13. [19]

    Y. Ohta. Data informativity of continuous-time systems by sampling using linear functionals. InProc. MTNS 2024,6 pp., 2024

  14. [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

  15. [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

  16. [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

  17. [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

  18. [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

  19. [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

  20. [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

  21. [27]

    Schumaker.Spline Functions: Basic Theory

    L. Schumaker.Spline Functions: Basic Theory. Cambridge Univ. Press, Cambridge, UK, 3rd edition, 2007

  22. [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

  23. [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

  24. [30]

    Tucsnak and G

    M. Tucsnak and G. Weiss.Observation and Control of Operator Semigroups. Birkh¨ auser, Basel, 2009

  25. [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

  26. [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

  27. [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

  28. [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

  29. [37]

    M. Wakaiki. Data informativity for stabilization of discrete-time infinite-dimensional systems. arXiv:2511.16008, 2025

  30. [38]

    J. Wang, S. Baldi, and H. J. van Waarde. Experiment design for continuous-time systems using generalized filtering. arXiv:2511.09386, 2025

  31. [39]

    J. C. Willems. From time series to linear system. Part I: Finite dimensional linear time-invariant systems.Automatica, 22:561– 580, 1986

  32. [40]

    J. C. Willems. From time series to linear system. Part II: Exact modelling.Automatica, 22:675–694, 1986

  33. [41]

    J. C. Willems. From time series to linear system. Part III: Approximate modelling.Automatica, 23:87–115, 1987

  34. [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

  35. [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

Pith tools

Reviewed August 3, 2026 · model on record in the stance chip above.