Hybrid system

A hybrid system is a dynamical system that exhibits both continuous and discrete dynamic behavior – a system that can both flow (described by a differential equation) and jump (described by a state machine, automaton, or a difference equation).^[1] Often, the term "hybrid dynamical system" is used instead of "hybrid system", to distinguish from other usages of "hybrid system", such as the combination neural nets and fuzzy logic, or of electrical and mechanical drivelines. A hybrid system has the benefit of encompassing a larger class of systems within its structure, allowing for more flexibility in modeling dynamic phenomena.

In general, the state of a hybrid system is defined by the values of the continuous variables and a discrete mode. The state changes either continuously, according to a flow condition, or discretely according to a control graph. Continuous flow is permitted as long as so-called invariants hold, while discrete transitions can occur as soon as given jump conditions are satisfied. Discrete transitions may be associated with events.

Hybrid systems have been used to model several cyber-physical systems, including physical systems with impact, logic-dynamic controllers, and even Internet congestion.

A canonical example of a hybrid system is the bouncing ball, a physical system with impact. Here, the ball (thought of as a point-mass) is dropped from an initial height and bounces off the ground, dissipating its energy with each bounce. The ball exhibits continuous dynamics between each bounce; however, as the ball impacts the ground, its velocity undergoes a discrete change modeled after an inelastic collision. A mathematical description of the bouncing ball follows. Let ${\displaystyle x_{1}}$ be the height of the ball and ${\displaystyle x_{2}}$ be the velocity of the ball. A hybrid system describing the ball is as follows:

When ${\displaystyle x\in C=\{x_{1}>0\}}$, flow is governed by ${\displaystyle {\dot {x}}_{1}=x_{2},{\dot {x}}_{2}=-g}$, where ${\displaystyle g}$ is the acceleration due to gravity. These equations state that when the ball is above ground, it is being drawn to the ground by gravity.

When ${\displaystyle x\in D=\{x_{1}=0\}}$, jumps are governed by ${\displaystyle x_{1}^{+}=x_{1},x_{2}^{+}=-\gamma x_{2}}$, where ${\displaystyle 0<\gamma <1}$ is a dissipation factor. This is saying that when the height of the ball is zero (it has impacted the ground), its velocity is reversed and decreased by a factor of ${\displaystyle \gamma }$. Effectively, this describes the nature of the inelastic collision.

The bouncing ball is an especially interesting hybrid system, as it exhibits Zeno behavior. Zeno behavior has a strict mathematical definition, but can be described informally as the system making an infinite number of jumps in a finite amount of time. In this example, each time the ball bounces it loses energy, making the subsequent jumps (impacts with the ground) closer and closer together in time.

It is noteworthy that the dynamical model is complete if and only if one adds the contact force between the ground and the ball. Indeed, without forces, one cannot properly define the bouncing ball and the model is, from a mechanical point of view, meaningless. The simplest contact model that represents the interactions between the ball and the ground, is the complementarity relation between the force and the distance (the gap) between the ball and the ground. This is written as ${\displaystyle 0\leq \lambda \perp x_{1}\geq 0.}$ Such a contact model does not incorporate magnetic forces, nor gluing effects. When the complementarity relations are in, one can continue to integrate the system after the impacts have accumulated and vanished: the equilibrium of the system is well-defined as the static equilibrium of the ball on the ground, under the action of gravity compensated by the contact force ${\displaystyle \lambda }$. One also notices from basic convex analysis that the complementarity relation can equivalently be rewritten as the inclusion into a normal cone, so that the bouncing ball dynamics is a differential inclusion into a normal cone to a convex set. See Chapters 1, 2 and 3 in Acary-Brogliato's book cited below (Springer LNACM 35, 2008). See also the other references on non-smooth mechanics.

There are approaches to automatically proving properties of hybrid systems (e.g., some of the tools mentioned below). Common techniques for proving safety of hybrid systems are computation of reachable sets, abstraction refinement, and barrier certificates.

Most verification tasks are undecidable,^[2] making general verification algorithms impossible. Instead, the tools are analyzed for their capabilities on benchmark problems. A possible theoretical characterization of this is algorithms that succeed with hybrid systems verification in all robust cases^[3] implying that many problems for hybrid systems, while undecidable, are at least quasi-decidable.^[4]

Two basic hybrid system modeling approaches can be classified, an implicit and an explicit one. The explicit approach is often represented by a hybrid automaton, a hybrid program or a hybrid Petri net. The implicit approach is often represented by guarded equations to result in systems of differential algebraic equations (DAEs) where the active equations may change, for example by means of a hybrid bond graph.

As a unified simulation approach for hybrid system analysis, there is a method based on DEVS formalism in which integrators for differential equations are quantized into atomic DEVS models. These methods generate traces of system behaviors in discrete event system manner which are different from discrete time systems. Detailed of this approach can be found in references [Kofman2004] [CF2006] [Nutaro2010] and the software tool PowerDEVS.

• HyEQ Toolbox: Hybrid system solver for MATLAB and Simulink
• PowerDEVS: General-purpose tool for DEVS (Discrete Event System) modeling and simulation oriented to the simulation of hybrid systems

• Ariadne: C++ library for (numerically rigorous) reachability analysis of nonlinear hybrid systems
• CORA: A MATLAB Toolbox for reachability analysis of cyber-physical systems, including hybrid systems
• Flow*: A tool for reachability analysis of nonlinear hybrid systems
• HyCreate: A tool for overapproximating reachability of hybrid automata
• HyPro: C++ library for state set representations for hybrid systems reachability analysis
• JuliaReach: A toolbox for set-based reachability

• C2E2: Nonlinear hybrid system verifier
• HyTech: Model checker for hybrid systems
• HSolver: Verification tool for hybrid systems
• KeYmaera: Theorem prover for hybrid systems
• PHAVer: Polyhedral hybrid automaton verifier
• S-TaLiRo: MATLAB toolbox for verification of hybrid systems with respect to temporal logic specifications

• SCOTS: Tool for the synthesis of correct-by-construction controllers for hybrid systems
• SpaceEx: State-space explorer

• Hybrid automaton
• Sliding mode control
• Variable structure system
• Variable structure control
• Joint spectral radius
• Cyber-physical system
• Behavior trees (artificial intelligence, robotics and control)
• Jump process (in the context of probability), an example of a (stochastic) hybrid system with zero flow component
• Piecewise-deterministic Markov process (PDMP), an example of a (stochastic) hybrid system and a generalization of the jump process
• Jump diffusion, an example of a (stochastic) hybrid system and a generalization of the PDMP

• Henzinger, Thomas A. (1996), "The Theory of Hybrid Automata", 11th Annual Symposium on Logic in Computer Science (LICS), IEEE Computer Society Press, pp. 278–292, archived from the original on 2010-01-27
• Alur, Rajeev; Courcoubetis, Costas; Halbwachs, Nicolas; Henzinger, Thomas A.; Ho, Pei-Hsin; Nicollin, Xavier; Olivero, Alfredo; Sifakis, Joseph; Yovine, Sergio (1995), "The algorithmic analysis of hybrid systems", Theoretical Computer Science, 138 (1): 3–34, doi:10.1016/0304-3975(94)00202-T, hdl:1813/6241, archived from the original on 2010-01-27
• Goebel, Rafal; Sanfelice, Ricardo G.; Teel, Andrew R. (2009), "Hybrid dynamical systems", IEEE Control Systems Magazine, 29 (2): 28–93, doi:10.1109/MCS.2008.931718, S2CID 46488751
• Acary, Vincent; Brogliato, Bernard (2008), Numerical Methods for Nonsmooth Dynamical Systems, Lecture Notes in Applied and Computational Mechanics, vol. 35, doi:10.1007/978-3-540-75392-6, ISBN 978-3-540-75391-9
• [Kofman2004] Kofman, E (2004), "Discrete Event Simulation of Hybrid Systems", SIAM Journal on Scientific Computing, 25 (5): 1771–1797, Bibcode:2004SJSC...25.1771K, CiteSeerX 10.1.1.72.2475, doi:10.1137/S1064827502418379
• [CF2006] Francois E. Cellier and Ernesto Kofman (2006), Continuous System Simulation (first ed.), Springer, ISBN 978-0-387-26102-7
• [Nutaro2010] James Nutaro (2010), Building Software for Simulation: Theory, Algorithms, and Applications in C++ (first ed.), Wiley
• Brogliato, Bernard; Tanwani, Aneel (2020), "Dynamical systems coupled with monotone set-valued operators: Formalisms, Applications, well-posedness, and stability" (PDF), SIAM Review, 62 (1): 3–129, doi:10.1137/18M1234795, S2CID 212727046

• IEEE CSS Committee on Hybrid Systems