REVIEW 3 major objections 6 minor 1 cited by
Issues with Input-Space Representation in Nonlinear Data-Based Dissipativity Estimation
T0 review · 3 major / 6 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read For nonlinear data-based dissipativity estimation, dense input grids need super-exponential samples, while PAC generalization bounds make the trade-off tractable.
desk verdict A genuinely useful tutorial with a solid negative result on δ-covering; the ML positive claim needs a sharper statement about the sampling distribution before it supports worst-case dissipativity. 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 argument runs on two objects. One is the δ-cover of the input space $U_{A1234}$—the set of inputs that are bounded in amplitude, bounded below in norm, have a long enough horizon, and are represented by $b$ orthonormal basis functions—whose sample count $K=(\rho b)^b$ is what makes the robustness guarantee cost super-exponential. The other is the PAC generalization error $|L_D - L_S|$ for the one-class SVM that fits the half-plane $\langle \Pi, \Gamma(u) \rangle \ge 0$, bounded by the cited one-class SVM generalization theory; this bound does not grow with $b$, so it decouples the sample complexity from how finely the input space is approximated. Also load-bearing are the conic-sector formulas in the appendix, which compute cone$(a,b)$ with $O(1)$ space, turning sampled trajectories into concrete dissipativity bounds.
What would settle it
Run the δ-covering estimator on the system $\hat{G}(\omega)=1/(j\omega+1)+1/4$ with the first four Legendre basis functions and $T=10$, sending $K$ to infinity; the paper reports the estimated cone converges to cone(0.68, 1.25) instead of the true cone(1/4, 5/4), directly demonstrating the failure of the robustness claim for practical basis counts. For the ML side, train the OC-SVM on samples from one distribution $D$ and test on inputs drawn from a different distribution; if the error exceeds the PAC bound, then the distributional premise is doing the work and the robustness guarantee is conditional rather than worst-case.
Extended reading notes
Core claim
The central claim is that δ-covering methods for input-space representation in nonlinear data-based dissipativity estimation face an intractable trade-off: either the sample complexity is too large to implement, or the robustness guarantee does not extend to the true input space $L_{2e}$. For a covering radius $\delta$ and sampling density $\rho$, the number of samples is $K=(\rho b)^b$, so $K$ is $O(b^b)$ for constant density; and because an $L_\infty$ cover requires at least $\rho^b$ points, every δ-covering scheme has super-exponential complexity in $b$. Existing case studies use $b \le 5$, and the paper shows that with small $b$ the estimated dissipativity can converge to wrong conic bounds even for a known linear system: with the first four Legendre bases and $T=10$, the estimate for the system in Equation (5) converges to cone(0.68, 1.25) rather than the true cone(1/4, 5/4). As an alternative, the paper proposes quantifying the robustness of ML-based dissipativity estimation using a PAC generalization bound for the one-class SVM problem; the bound does not scale with $b$ or input/output dimensions, so a large number of basis functions can be used while keeping the sample complexity fixed.
Load-bearing premise
The load-bearing premise is Assumption 6: the input trajectories used in training are drawn from the same probability distribution $D$ that describes the inputs the system will actually see, and the PAC guarantee says nothing about inputs that are not drawn from that distribution.
Editorial extensions
If this is right
- For any δ-covering scheme, the sample count is at least $\rho^b$ in the $L_\infty$ norm and $O(b^b)$ at constant density, so robust nonlinear dissipativity estimation by grid search is out of reach for input spaces large enough to approximate $L_{2e}$.
- Keeping $b$ small to make δ-covering feasible changes which inputs are represented, and the estimated dissipativity can converge to wrong values even in the infinite-data limit; data collection and basis choice must be aligned with the true input space.
- Using the OC-SVM PAC bound, the number of samples needed to certify dissipativity with high probability does not grow with the number of basis functions, so input representations with $b=100$ or more become practical.
- Uniformly sampling bounded basis coefficients and using Wiener-process inputs with a range of time horizons both yield good conic-sector estimates in the case studies, the latter without requiring amplitude or norm bounds.
- A single Wiener process trimmed at multiple horizons can estimate dissipativity directly on truncated signals, removing the assumption that a long enough horizon guarantees dissipativity.
Reading between the lines
- The PAC certificate depends on the distribution $D$, so if the inputs that threaten closed-loop stability are rare or adversarial, the probabilistic guarantee can be arbitrary; a natural extension is to shape $D$ to cover the frequency ranges that matter for stability, or to combine the certificate with adversarial search over basis coefficients.
- The super-exponential sample count in $K=(\rho b)^b$ suggests that any input-space covering approach, not just δ-covering, will fail for high-dimensional bases; the ML route reframes the problem from covering to coverage, which is why it scales.
- The $O(1)$-space conic formulas in the appendix could be used as streaming estimators, updating dissipativity bounds incrementally as data arrive, which would fit online data-based control loops.
- A testable extension: for a nonlinear system whose dissipativity is determined by a narrow frequency band, compare the PAC certificate trained on uniform basis-coefficient sampling against one trained on a distribution that emphasizes that band; the certificate should tighten, suggesting $D$ should encode prior physical knowledge.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This manuscript is a tutorial-style study of data-based dissipativity estimation for nonlinear systems. It formalizes the assumptions commonly used to restrict the input space L2e to a finite-dimensional set UA1234, then analyzes two families of methods. For δ-covering methods (Section IV), it derives the sample count K=(ρb)^b, argues that this is super-exponential in the number of basis functions b, and gives numerical examples showing that the estimated conic bounds depend strongly on the choice of basis and time horizon. For machine-learning methods (Section V), it proposes applying a PAC generalization bound from Schölkopf et al. to an OC-SVM formulation of dissipativity, claiming that this yields a probabilistic robustness guarantee whose sample complexity is independent of b. Numerical experiments with an LTI system and a pendulum illustrate both the basis-dependence of δ-covering and the behavior of randomly sampled inputs.
Significance. The negative result on δ-covering methods is a useful tutorial contribution: the paper makes explicit that the super-exponential sample complexity, not just the input dimension, is the binding limitation, and the Legendre-basis examples in Figure 1 concretely demonstrate basis-dependent failure even for b=4. The proposed PAC-based approach is a natural and promising direction, and the observation that the generalization bound does not explicitly grow with b is valuable. However, the positive claim is not established as stated: the PAC guarantee is distribution-relative rather than worst-case, and it does not yet provide the certificate required by the Dissipativity Theorem. The manuscript also contains a technical error in the lower-bound argument of Section IV.A, although that conclusion appears repairable. Overall, the paper is likely to be of interest to the data-driven control community, but the central robustness claim needs substantial revision before the manuscript can be accepted.
major comments (3)
- [Section V, Assumption 6] The PAC bound displayed in Section V controls P_{u∼D}(⟨Π,Γ(u)⟩<ρ−γ) for a distribution D, but Definition 1 and the Dissipativity Theorem require Eq. (1) to hold for every u∈U_e. The paper concedes that the meaning of D 'is not immediately clear.' Neither proposed interpretation resolves the gap: uniform sampling over b basis coefficients has no well-defined limit in L2e as b→∞, and the Wiener-process sampling demonstrably under-represents high-frequency signals for large T (Figure 3). Therefore, the claimed 'tractable trade-off between robustness and sample complexity' is at present a statement about a chosen sampling distribution, not a certificate of robustness on the input space relevant to stability. The paper should either restrict the claim to a clearly defined distributional robustness notion and show that this suffices for the intended closed-loop stability guarantee, or construct a D whose support equals UA1234 (with a well-defined limit as b,T→∞) and prove the corresponding guarantee.
- [Section IV.A, Eq. (4)] The lower-bound argument uses the statement B_L2_δ ⊆ B_L∞_δ. This inclusion is false for continuous-time signals: an L2 ball of radius δ can contain a narrow spike of amplitude M with L2 norm δ and L∞ norm M, where M can be much larger than δ (take support length δ^2/M^2). The correct inclusion is B_L∞_δ ⊆ B_L2_{√T δ} up to a horizon-dependent factor, so the derivation as written is invalid. The Ω(ρ^b) lower bound can likely be recovered by an L2 covering-number or volume argument in the b-dimensional coefficient space, and the paper's negative conclusion may survive, but the proof needs to be corrected because the current inclusion is load-bearing for the claim that all δ-covering methods are intractable.
- [Section V, Eq. (8)] The cited generalization bound from [22] is applied to the solution of the hard OC-SVM problem (8), but the manuscript does not verify that the hypotheses of the theorem in [22] hold for this learned solution. In particular, the feature map Γ(u) is not normalized, the kernel associated with the linear product ⟨Π,Γ(u)⟩ is not specified, and the constant c=10^3 is imported without justification. Because this bound is the central evidence for the claimed b-independent sample complexity, the authors should either verify that the solution of (8) satisfies the conditions of the theorem, or derive a suitable bound directly for the hypothesis class defined by (8).
minor comments (6)
- [Section II] There are several typos: 'cardonality' should be 'cardinality', and Assumption 5 uses 'Lipschiz' instead of 'Lipschitz'.
- [Section IV.A, Eq. (3)] The set-builder notation for U′ is unclear; the sampling grid and the definition of the coefficients α_i should be stated explicitly.
- [Section V and Figures 3-4] The Wiener process is consistently misspelled as 'Weiner'; please correct throughout.
- [Figure 2 caption] The caption does not explain how uniform sampling over UA1234 is implemented for different b, nor what exactly is plotted as the 'upper bound' and 'lower bound' in each panel.
- [Appendix, Eqs. (10)-(12)] The ± notation inside M(u_i) and in the subsequent formulas is ambiguous; the text should clarify that both sign choices are evaluated and the extremal value is selected.
- [References] Reference [12] lists 'G. Allgöwer' while the other Allgöwer references use 'F. Allgöwer'; this should be corrected.
Circularity Check
No significant circularity: the PAC certificate is an external, distribution-relative bound; the sole same-author citation is a non-load-bearing appendix reference.
full rationale
The paper's central negative claim (Sec. IV) is a counting/volume bound on the size of delta-covers, derived from the definition of a delta-cover and the sampling scheme in Eqs. (3)-(4); it is not an output fitted into the claim. The positive claim (Sec. V) applies an external PAC generalization bound from Schoelkopf et al. [22] to the OC-SVM formulation from [13]. The bound is a standard conditional guarantee: if training data are drawn from D, then with high probability future D-draws are correctly classified. The paper explicitly concedes the limitation: 'This generalization error is only achieved if the data is sampled according to distribution D, the meaning of which is not immediately clear.' That is a distributional premise and a mismatch with the worst-case Dissipativity Theorem, but not a circular reduction: the guarantee is not defined in terms of the conclusion it is used to support. The only same-author citation is [23] (Bridgeman and Forbes) in the appendix, used for two known ways to compute extremal conic bounds; it is not load-bearing for the intractability or PAC results. Technical concerns such as the L2/L-infinity ball inclusion direction in Sec. IV and the average-case versus worst-case gap in Sec. V are correctness risks, not circularity. No step reduces to its inputs by construction.
Assumptions & free parameters
free parameters (4)
- Covering radius δ
- Lipschitz constant L
- PAC margin γ =
γ=ρ
- Input-space parameters b, T, ū, ε
assumptions (9)
- domain assumption Assumption 1: the amplitude of permissible inputs is absolutely bounded above by ū.
- domain assumption Assumption 2: the induced norm of permissible inputs is bounded below by ε.
- ad hoc to paper Assumption 3: if Eq. 1 holds for one sufficiently large T, it holds for all T.
- domain assumption Assumption 4: the input signal can be represented by a finite number b of orthonormal basis functions.
- domain assumption Assumption 5: the unknown operator G is Lipschitz continuous with known constant L.
- domain assumption Assumption 6: data u_i are sampled i.i.d. from a probability distribution D over Ue.
- ad hoc to paper The inclusion B_L2_δ ⊆ B_L∞_δ used to lower-bound sample complexity.
- standard math Parseval's theorem and Fourier series representation of orthonormal basis functions.
- standard math The one-class SVM generalization bound of [22].
Cite this review
Pith. "Pith review of Issues with Input-Space Representation in Nonlinear Data-Based Dissipativity Estimation." pith.science (2026). https://pith.science/paper/GL5WIDH6
@misc{pith2026241113404,
author = {Pith},
title = {Pith review of: Issues with Input-Space Representation in Nonlinear Data-Based Dissipativity Estimation},
year = {2026},
howpublished = {\url{https://pith.science/paper/GL5WIDH6}},
note = {Machine review of arXiv:2411.13404}
}
read the original abstract
In data-based control, dissipativity can be a powerful tool for attaining stability guarantees for nonlinear systems if that dissipativity can be inferred from data. This work provides a tutorial on several existing methods for data-based dissipativity estimation of nonlinear systems. The interplay between the underlying assumptions of these methods and their sample complexity is investigated. It is shown that methods based on delta-covering result in an intractable trade-off between sample complexity and robustness. A new method is proposed to quantify the robustness of machine learning-based dissipativity estimation. It is shown that this method achieves a more tractable trade-off between robustness and sample complexity. Several numerical case studies demonstrate the results.
Figures
Forward citations
Cited by 1 Pith paper
-
Learning the Integral Quadratic Constraints on Plant-Model Mismatch
A one-class SVM can learn frequency-domain IQC bounds on nonlinear plant-model mismatch from sampled trajectories, demonstrated on a time-delay mismatch and a two-phase reactor.
Reference graph
Works this paper leans on
-
[10]
Some problems aris ing in controller design from big data via input-output methods,
J. M. Montenbruck and F. Allg¨ ower, “Some problems aris ing in controller design from big data via input-output methods,” 55th IEEE Conf. Decis. Ctrl. , pp. 6525–6530, Dec. 2017
work page 2017
-
[22]
Estimating the support of a high-dimensional distribu- tion,
B. Sch¨ olkopf, J. C. Platt, J. Shawe-Taylor, A. J. Smola , and R. C. Williamson, “Estimating the support of a high-dimensional distribu- tion,” Neural Comp. , vol. 13, pp. 1443–1471, 2001
work page 2001
-
[1]
Data-enabled policy optimization for direct adaptive learning of the LQR,
F. Zhao, F. D¨ orfler, A. Chiuso, and K. Y ou, “Data-enabled policy optimization for direct adaptive learning of the LQR,” Arχ iv, pp. 1– 16, Apr. 2024
work page 2024
-
[2]
Guarantees fo r data- driven control of nonlinear systems using semidefinite prog ramming: A survey,
T. Martin, T. B. Sch¨ on, and F. Allg¨ ower, “Guarantees fo r data- driven control of nonlinear systems using semidefinite prog ramming: A survey,” Ann. Rev. Ctrl. , vol. 56, pp. 2–23, 2023
work page 2023
-
[3]
Vidyasagar, Input-output analysis of large-scale interconnected systems
M. Vidyasagar, Input-output analysis of large-scale interconnected systems. Berlin, Germany: Springer-V erlag, 1981
work page 1981
-
[4]
Synthesis of positive rea l H2 controllers,
J. C. Geromel and P . B. Gapski, “Synthesis of positive rea l H2 controllers,” IEEE Tran. Aut. Ctrl. , vol. 42, no. 7, pp. 988–992, 1997
work page 1997
-
[5]
Dissipative dynamical systems: b asic input- output and state properties,
D. Hill and P . Moylan, “Dissipative dynamical systems: b asic input- output and state properties,” J. Franklin Inst. , vol. 309, no. 5, pp. 327–357, May 1980
work page 1980
-
[6]
A note on persistancy of excitation,
J. C. Willems, R. Rapisarda, I. Markovsky, and B. L. De Moo r, “A note on persistancy of excitation,” Sys. Ctrl. Lett. , vol. 54, pp. 325– 329, 2005
work page 2005
Show all 23 references
-
[7]
Provably robus t verification of dissipativity properties from data,
A. Koch, J. Berberich, and F. Allg¨ ower, “Provably robus t verification of dissipativity properties from data,” IEEE Tran. Aut. Ctrl. , vol. 67, no. 8, pp. 4248–4255, 2022
2022
-
[8]
Dissipativity verification with guarantees for polynomial systems from noisy input-state data,
T. Martin and F. Allg¨ ower, “Dissipativity verification with guarantees for polynomial systems from noisy input-state data,” IEEE Ctrl. Sys. Lett, vol. 5, no. 4, pp. 1399–1404, Oct. 2021
2021
-
[9]
Data-driven system analysis of nonlinear systems u sing polyno- mial approximation,
——, “Data-driven system analysis of nonlinear systems u sing polyno- mial approximation,” IEEE Tran. Aut. Ctrl. , vol. 69, no. 7, pp. 4261– 4274, 2024
2024
-
[11]
Determi ning dis- sipation inequalities from input-output samples,
A. Romer, J. M. Montenbruck, and F. Allg¨ ower, “Determi ning dis- sipation inequalities from input-output samples,” IF AC-PapersOnline, vol. 50, no. 1, pp. 7789–7794, 2017
2017
-
[12]
Data-driven inf erence of passivity properties via Gaussian process optimization,
A. Romer, S. Trimpe, and G. Allg¨ ower, “Data-driven inf erence of passivity properties via Gaussian process optimization,” Euro. Ctrl. Conf., pp. 29–35, Jun. 2019
2019
-
[13]
Input-output data-driven co ntrol through dissipativity learning,
W. Tang and P . Daoutidis, “Input-output data-driven co ntrol through dissipativity learning,” Amer . Ctrl. Conf., pp. 4217–4222, Jul. 2019
2019
-
[14]
Dissipativity learning control (dlc): A framewor k of input- output data-driven control,
——, “Dissipativity learning control (dlc): A framewor k of input- output data-driven control,” Amer . Ctrl. Conf. , vol. 130, no. 2, pp. 1–16, 2019
2019
-
[15]
Dissipativity learning control (DLC): Theoretic al foundations of input-output data-driven model-free control,
——, “Dissipativity learning control (DLC): Theoretic al foundations of input-output data-driven model-free control,” Sys. Ctrl. Lett. , vol. 147, pp. 1–9, 2021
2021
-
[16]
Shalev-Shwartz and S
S. Shalev-Shwartz and S. Ben-David, Understanding machine learn- ing: From theory to algorithms . New Y ork, NY , USA: Cambridge University Press, 2014
2014
-
[17]
On the input-output stability of time-varyi ng nonlinear feedback systems parts I & II,
G. Zames, “On the input-output stability of time-varyi ng nonlinear feedback systems parts I & II,” IEEE Tran. Aut. Ctrl. , vol. AC-11, no. 2-3, pp. 228–238,465–476, Apr./Jul. 1966
1966
-
[18]
Regeneration theory,
H. Nyquist, “Regeneration theory,” Bell Sys. Tech. J. , vol. 11, no. 1, pp. 126–147, 1932
1932
-
[19]
Computing controlled invariant sets from da ta using convex optimization,
M. Korda, “Computing controlled invariant sets from da ta using convex optimization,” SIAM J. Ctrl. Opt. , vol. 58, no. 5, pp. 2871– 2899, 2020
2020
-
[20]
A pivoting algorithm for convex h ulls and vertex enumeration of arrangements and polyhedra,
D. Avis and K. Fukuda, “A pivoting algorithm for convex h ulls and vertex enumeration of arrangements and polyhedra,” Discrete Comput. Geom., vol. 8, pp. 297–313, 1992
1992
-
[21]
L2e-gain analysis of nonlinear systems and nonlinear state feedback H∞ control,
A. J. V an Der Schaft, “ L2e-gain analysis of nonlinear systems and nonlinear state feedback H∞ control,” IEEE Tran. Aut. Ctrl. , vol. 37, no. 6, pp. 770–784, 1992
1992
-
[23]
Conic-sector-based c ontrol to circumvent passivity violations,
L. J. Bridgeman and J. R. Forbes, “Conic-sector-based c ontrol to circumvent passivity violations,” Int. J. Ctrl. , vol. 87, no. 8, pp. 1467– 1477, 2014
2014
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.