pith. sign in

arxiv: 2407.11256 · v2 · submitted 2024-07-15 · 📡 eess.SY · cs.SY

Controlled Invariant Sets for Gaussian Process State Space Models

Pith reviewed 2026-05-23 22:33 UTC · model grok-4.3

classification 📡 eess.SY cs.SY
keywords gaussian processcontrolled invariant setssemidefinite programmingstate feedbacknonlinear systemsprobabilistic safetyquadrotor
0
0 comments X

The pith

A semidefinite programming scheme designs state-feedback controllers that maximize the probability trajectories remain inside probabilistic controlled invariant sets for Gaussian process state-space models while respecting input constraints

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

The paper computes probabilistic controlled invariant sets for nonlinear systems whose unknown dynamics are captured by data-driven Gaussian process state-space models. It formulates a semidefinite program that synthesizes state-feedback controllers to maximize the probability that trajectories stay inside these sets. The same program enforces input constraints. The approach is tested on a quadrotor both in simulation and on physical hardware.

Core claim

We compute probabilistic controlled invariant sets for nonlinear systems using Gaussian process state space models, which are data-driven models that account for unmodeled and unknown nonlinear dynamics. We propose a semidefinite programming scheme for designing state-feedback controllers that maximize the probability of the trajectories staying within a probabilistic controlled invariant set while satisfying input constraints. The results are validated on a quadrotor, both in simulation and on a physical platform.

What carries the argument

semidefinite programming scheme that maximizes the probability of remaining inside a probabilistic controlled invariant set defined from a Gaussian process state-space model

Load-bearing premise

The Gaussian process state-space model learned from data sufficiently captures the unmodeled nonlinear dynamics so that the computed probabilistic invariant sets remain valid for the true system.

What would settle it

Run the designed controller on the physical quadrotor and measure the empirical fraction of time trajectories exit the computed set; a large mismatch with the designed probability would indicate the GP model does not support the invariance claim.

Figures

Figures reproduced from arXiv: 2407.11256 by Bingzhuo Zhong, Majid Zamani, Marco Caccamo, Murat Arcak, Paul Griffioen.

Figure 2
Figure 2. Figure 2: to validate the probabilistic guarantee proposed in this paper via Monte Carlo analysis. We consider state constraints xx, xy ∈ [−5, 5] m and vx, vy ∈ [−7, 7] m/s, and input constraints ux, uy ∈ [−5, 5] m/s2 [PITH_FULL_IMAGE:figures/full_fig_p008_2.png] view at source ↗
Figure 1
Figure 1. Figure 1: A quadrotor moving on an x-y plane. model the quadrotor dynamics as a GPSSM accord￾ing to (1) with states and inputs described according to [13]. The system states and control inputs are given by x ≜ h xx vx xy vy iT and u ≜ h ux uy iT , where xi , vi , and ui are the position, velocity, and acceleration of the quadrotor on the i th axis, respectively, with i ∈ {x, y}. 5.1 High Fidelity SIMULINK Model Firs… view at source ↗
Figure 3
Figure 3. Figure 3: Simulation results for the high fidelity [PITH_FULL_IMAGE:figures/full_fig_p009_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Left: Physical quadrotor for the experiment. Right: [PITH_FULL_IMAGE:figures/full_fig_p009_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Top: Evolution of the quadrotor’s positions in the [PITH_FULL_IMAGE:figures/full_fig_p010_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Sequences of the quadrotor’s velocity and control [PITH_FULL_IMAGE:figures/full_fig_p010_6.png] view at source ↗
read the original abstract

We compute probabilistic controlled invariant sets for nonlinear systems using Gaussian process state space models, which are data-driven models that account for unmodeled and unknown nonlinear dynamics. We propose a semidefinite programming scheme for designing state-feedback controllers that maximize the probability of the trajectories staying within a probabilistic controlled invariant set while satisfying input constraints. The results are validated on a quadrotor, both in simulation and on a physical platform.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 0 minor

Summary. The manuscript proposes computing probabilistic controlled invariant sets for nonlinear systems via Gaussian process (GP) state-space models that capture unmodeled dynamics from data. It introduces a semidefinite programming (SDP) formulation to synthesize state-feedback controllers maximizing the probability that closed-loop trajectories remain inside the probabilistic invariant set while respecting input constraints. Validation is reported on a quadrotor both in simulation and on hardware.

Significance. If the probabilistic invariance guarantees transfer from the GP model to the true nonlinear plant, the SDP-based design would offer a principled data-driven route to safe control under uncertainty, with direct relevance to robotics applications. The combination of GP uncertainty quantification with SDP optimization and the inclusion of physical-platform experiments are strengths that could support broader adoption if the model-mismatch issue is resolved.

major comments (2)
  1. [Abstract (and implied theoretical sections)] The SDP designs invariance and probability bounds explicitly for the learned GP state-space model, yet the central claim asserts validity for the underlying nonlinear system. No section supplies a rigorous bound on the GP approximation error, a robustness margin, or a mismatch analysis that would guarantee the computed probabilities remain valid on the true dynamics; this gap is load-bearing for the transfer of the invariance result.
  2. [Abstract (validation paragraph)] The physical quadrotor validation is cited, but the manuscript does not report quantitative comparison between observed trajectory frequencies and the GP-predicted probabilities, nor does it include an error analysis or statistical test confirming that the empirical invariance rates align with the computed values. This leaves the hardware claim unsupported for the probabilistic guarantee.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive comments. We address each major comment below, with proposed revisions to improve clarity and validation.

read point-by-point responses
  1. Referee: [Abstract (and implied theoretical sections)] The SDP designs invariance and probability bounds explicitly for the learned GP state-space model, yet the central claim asserts validity for the underlying nonlinear system. No section supplies a rigorous bound on the GP approximation error, a robustness margin, or a mismatch analysis that would guarantee the computed probabilities remain valid on the true dynamics; this gap is load-bearing for the transfer of the invariance result.

    Authors: The probabilistic invariance guarantees and SDP formulation are derived explicitly for the GP state-space model. The manuscript does not provide a rigorous bound on GP approximation error or a formal transfer result to the true nonlinear dynamics. We will revise the abstract and introduction to state that the theoretical results hold for the learned GP model, with the physical experiments serving as empirical validation rather than a formal guarantee. revision: yes

  2. Referee: [Abstract (validation paragraph)] The physical quadrotor validation is cited, but the manuscript does not report quantitative comparison between observed trajectory frequencies and the GP-predicted probabilities, nor does it include an error analysis or statistical test confirming that the empirical invariance rates align with the computed values. This leaves the hardware claim unsupported for the probabilistic guarantee.

    Authors: We agree that a direct quantitative comparison would strengthen the hardware results. In the revised manuscript we will add an analysis comparing the empirical frequency of invariant trajectories from the hardware experiments to the probabilities obtained from the SDP, including basic error measures. revision: yes

Circularity Check

0 steps flagged

No circularity: SDP scheme computes invariants directly on the GP model without self-referential reduction

full rationale

The derivation consists of formulating an SDP that optimizes a state-feedback law to maximize an invariance probability under the posterior of a learned GP state-space model, subject to input constraints. The probabilistic sets and controller are defined and solved explicitly from the GP mean and covariance; the result is not obtained by fitting a parameter to a subset and relabeling it as a prediction, nor does any load-bearing step invoke a self-citation chain or uniqueness theorem that collapses to the authors' prior ansatz. External quadrotor validation supplies an independent check. The central claim therefore remains a constructive computational procedure rather than a tautology.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract-only review yields no identifiable free parameters, axioms, or invented entities.

pith-pipeline@v0.9.0 · 5598 in / 844 out tokens · 17114 ms · 2026-05-23T22:33:38.405626+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

37 extracted references · 37 canonical work pages

  1. [1]

    Probabilistic Reachability and Safety for Controlled Discrete Time Stochastic Hybrid Systems

    Alessandro Abate, Maria Prandini, John Lygeros, and Shankar Sastry. Probabilistic Reachability and Safety for Controlled Discrete Time Stochastic Hybrid Systems. Automatica, 44(11):2724–2734, 2008

  2. [2]

    On Estimation Error Bounds for Receding-Horizon Filters Using Quadratic Boundedness

    Angelo Alessandri, Marco Baglietto, and Giorgio Battistelli. On Estimation Error Bounds for Receding-Horizon Filters Using Quadratic Boundedness. IEEE Transactions on Automatic Control, 49(8):1350–1355, 2004

  3. [3]

    Set Invariance in Control

    Franco Blanchini. Set Invariance in Control. Automatica, 35(11):1747–1767, 1999

  4. [4]

    Linear Matrix Inequalities in System and Control Theory

    Stephen Boyd, Laurent El Ghaoui, Eric Feron, and Venkataramanan Balakrishnan. Linear Matrix Inequalities in System and Control Theory . SIAM, 1994

  5. [5]

    Quadratic Boundedness of Nonlinear Dynamical Systems

    Mark L Brockman and Martin Corless. Quadratic Boundedness of Nonlinear Dynamical Systems. In Proceedings of 1995 34th IEEE Conference on Decision and Control, volume 1, pages 504–509. IEEE, 1995

  6. [6]

    Quadratic Boundedness of Nominally Linear Systems

    Mark L Brockman and Martin Corless. Quadratic Boundedness of Nominally Linear Systems. International Journal of Control , 71(6):1105–1117, 1998

  7. [7]

    Path Planning Using Positive Invariant Sets

    Claus Danielson, Avishai Weiss, Karl Berntorp, and Stefano Di Cairano. Path Planning Using Positive Invariant Sets. In 2016 IEEE 55th Conference on Decision and Control (CDC), pages 5986–5991. IEEE, 2016

  8. [8]

    Gaussian Processes for Data-Efficient Learning in Robotics and Control

    Marc Peter Deisenroth, Dieter Fox, and Carl Edward Rasmussen. Gaussian Processes for Data-Efficient Learning in Robotics and Control. IEEE Transactions on Pattern Analysis and Machine Intelligence , 37(2):408–423, 2013

  9. [9]

    A Stochastic Games Framework for Verification and Control of Discrete Time Stochastic Hybrid Systems

    Jerry Ding, Maryam Kamgarpour, Sean Summers, Alessandro Abate, John Lygeros, and Claire Tomlin. A Stochastic Games Framework for Verification and Control of Discrete Time Stochastic Hybrid Systems. Automatica, 49(9):2665–2674, 2013

  10. [10]

    Identification of Gaussian Process State Space Models

    Stefanos Eleftheriadis, Tom Nicholson, Marc Deisenroth, and James Hensman. Identification of Gaussian Process State Space Models. Advances in Neural Information Processing Systems, 30, 2017

  11. [11]

    Variational Gaussian Process State-Space Models

    Roger Frigola, Yutian Chen, and Carl Edward Rasmussen. Variational Gaussian Process State-Space Models. Advances in Neural Information Processing Systems , 27, 2014

  12. [12]

    Bayesian Inference and Learning in Gaussian Process State-Space Models with Particle MCMC

    Roger Frigola, Fredrik Lindsten, Thomas B Sch¨ on, and Carl Edward Rasmussen. Bayesian Inference and Learning in Gaussian Process State-Space Models with Particle MCMC. Advances in Neural Information Processing Systems , 26, 2013

  13. [13]

    Analytical Design and Experimental Verification of Geofencing Control for Aerial Applications

    Azad Ghaffari. Analytical Design and Experimental Verification of Geofencing Control for Aerial Applications. IEEE/ASME Transactions on Mechatronics , 26:1106–1117, 2021. 10

  14. [14]

    Probabilistic Invariance for Gaussian Process State Space Models

    Paul Griffioen, Alex Devonport, and Murat Arcak. Probabilistic Invariance for Gaussian Process State Space Models. In Learning for Dynamics and Control Conference , pages 458–468. PMLR, 2023

  15. [15]

    On a Correspondence Between Probabilistic and Robust Invariant Sets for Linear Systems

    Lukas Hewing, Andrea Carron, Kim P Wabersich, and Melanie N Zeilinger. On a Correspondence Between Probabilistic and Robust Invariant Sets for Linear Systems. In 2018 European Control Conference (ECC) , pages 1642–

  16. [16]

    Discrete Time Stochastic Hybrid Dynamical Games: Verification & Controller Synthesis

    Maryam Kamgarpour, Jerry Ding, Sean Summers, Alessandro Abate, John Lygeros, and Claire Tomlin. Discrete Time Stochastic Hybrid Dynamical Games: Verification & Controller Synthesis. In 2011 50th IEEE Conference on Decision and Control and European Control Conference , pages 6122–6127. IEEE, 2011

  17. [17]

    Probabilistic Set Invariance and Ultimate Boundedness

    Ernesto Kofman, Jos´ e A De Don´ a, and Maria M Seron. Probabilistic Set Invariance and Ultimate Boundedness. Automatica, 48(10):2670–2676, 2012

  18. [18]

    Continuous-Time Probabilistic Ultimate Bounds and Invariant Sets: Computation and Assignment

    Ernesto Kofman, Jos´ e A De Don´ a, Maria M Seron, and Noelia Pizzi. Continuous-Time Probabilistic Ultimate Bounds and Invariant Sets: Computation and Assignment. Automatica, 71:98–105, 2016

  19. [19]

    YALMIP: A Toolbox for Modeling and Optimization in MATLAB

    Johan Lofberg. YALMIP: A Toolbox for Modeling and Optimization in MATLAB. In 2004 IEEE International Conference on Robotics and Automation , pages 284–289. IEEE, 2004

  20. [20]

    The MOSEK Optimization Toolbox for MATLAB Manual

    MOSEK ApS. The MOSEK Optimization Toolbox for MATLAB Manual. Version 9.3.6 , 2019

  21. [21]

    Data-Driven Synthesis of Robust Invariant Sets and Controllers

    Sampath Kumar Mulagaleti, Alberto Bemporad, and Mario Zanon. Data-Driven Synthesis of Robust Invariant Sets and Controllers. IEEE Control Systems Letters , 6:1676–1681, 2021

  22. [22]

    Probabilistic ultimate bounds and invariant sets in nonlinear systems

    Noelia Pizzi, Ernesto Kofman, Dami´ an Edgardo Marelli, Jos´ e Adrian De Don´ a, and Maria M Seron. Probabilistic ultimate bounds and invariant sets in nonlinear systems. Automatica, 133:109853, 2021

  23. [23]

    Gaussian Processes for Machine Learning , volume 2

    Carl Edward Rasmussen and Christopher KI Williams. Gaussian Processes for Machine Learning , volume 2. MIT Press Cambridge, MA, 2006

  24. [24]

    Exploiting Model Uncertainty Estimates for Safe Dynamic Control Learning

    Jeff Schneider. Exploiting Model Uncertainty Estimates for Safe Dynamic Control Learning. Advances in Neural Information Processing Systems, 9, 1996

  25. [25]

    Direct and Inverse Addition in Convex Analysis and Applications

    Alberto Seeger. Direct and Inverse Addition in Convex Analysis and Applications. Journal of Mathematical Analysis and Applications, 148(2):317–349, 1990

  26. [26]

    An Engineering Method for Safety Region Development

    Danbing Seto and Lui Sha. An Engineering Method for Safety Region Development. Citeseer, 1999

  27. [27]

    Computationally Efficient Bayesian Learning of Gaussian Process State Space Models

    Andreas Svensson, Arno Solin, Simo S¨ arkk¨ a, and Thomas Sch¨ on. Computationally Efficient Bayesian Learning of Gaussian Process State Space Models. In Artificial Intelligence and Statistics , pages 213–221. PMLR, 2016

  28. [28]

    Learning for Safety-Critical Control with Control Barrier Functions

    Andrew Taylor, Andrew Singletary, Yisong Yue, and Aaron Ames. Learning for Safety-Critical Control with Control Barrier Functions. In Learning for Dynamics and Control , pages 708–717. PMLR, 2020

  29. [29]

    State- Space Inference and Learning with Gaussian Processes

    Ryan Turner, Marc Deisenroth, and Carl Rasmussen. State- Space Inference and Learning with Gaussian Processes. In Proceedings of the Thirteenth International Conference on Artificial Intelligence and Statistics , pages 868–875. JMLR Workshop and Conference Proceedings, 2010

  30. [30]

    Learning Stochastically Stable Gaussian Process State–Space Models

    Jonas Umlauft and Sandra Hirche. Learning Stochastically Stable Gaussian Process State–Space Models. IFAC Journal of Systems and Control , 12:100079, 2020

  31. [31]

    Safe Learning of Quadrotor Dynamics Using Barrier Certificates

    Li Wang, Evangelos A Theodorou, and Magnus Egerstedt. Safe Learning of Quadrotor Dynamics Using Barrier Certificates. In 2018 IEEE International Conference on Robotics and Automation (ICRA) , pages 2460–2465. IEEE, 2018

  32. [32]

    Learning While Tracking: A Practical System Based on Variational Gaussian Process State-Space Model and Smartphone Sensory Data

    Ang Xie, Feng Yin, Bo Ai, Sha Zhang, and Shuguang Cui. Learning While Tracking: A Practical System Based on Variational Gaussian Process State-Space Model and Smartphone Sensory Data. In 2020 IEEE 23rd International Conference on Information Fusion (FUSION) , pages 1–7. IEEE, 2020

  33. [33]

    A Gaussian Process State Space Model Fusion Physical Model and Residual Analysis for Fatigue Evaluation

    Aijun Yin, Junlin Zhou, and Tianyou Liang. A Gaussian Process State Space Model Fusion Physical Model and Residual Analysis for Fatigue Evaluation. Sensors, 22(7):2540, 2022

  34. [34]

    Gaussian Process-Integrated State Space Model for Continuous Joint Angle Prediction from EMG and Interactive Force in a Human-Exoskeleton System

    Yan Zeng, Jiantao Yang, and Yuehong Yin. Gaussian Process-Integrated State Space Model for Continuous Joint Angle Prediction from EMG and Interactive Force in a Human-Exoskeleton System. Applied Sciences , 9(8):1711, 2019

  35. [35]

    The Schur Complement and Its Applications, volume 4

    Fuzhen Zhang. The Schur Complement and Its Applications, volume 4. Springer Science & Business Media, 2006

  36. [36]

    Automata-Based Controller Synthesis for Stochastic Systems: A Game Framework via Approximate Probabilistic Relations

    Bingzhuo Zhong, Abolfazl Lavaei, Majid Zamani, and Marco Caccamo. Automata-Based Controller Synthesis for Stochastic Systems: A Game Framework via Approximate Probabilistic Relations. Automatica, 147:110696, 2023

  37. [37]

    Synthesizing Safety Controllers for Uncertain Linear Systems: A Direct Data-Driven Approach

    Bingzhuo Zhong, Majid Zamani, and Marco Caccamo. Synthesizing Safety Controllers for Uncertain Linear Systems: A Direct Data-Driven Approach. In 2022 IEEE Conference on Control Technology and Applications (CCTA), pages 1278–1284. IEEE, 2022. 11 Paul Griffioen is currently a post- doctoral researcher in the Department of Electrical Engineering and Com- put...