The following dynamic systems are implemented in CORA. More information can be found in Section 4 in the CORA manual.

• Linear Systems

  In CORA, linear systems have the form

  $$ \dot{x}(t) = A x(t) + B u(t) + c + w(t),$$ $$ y(t) = C x(t) + D u(t) + k + v(t), $$

  where \(x(t) \in \mathbb{R}^n\) is the system state, \(u(t)\in \mathbb{R}^m\) is the system input, \(w(t) \in \mathbb{R}^n\) is the disturbance, \(y(t) \in \mathbb{R}^p\) is the system output, \(v(t) \in \mathbb{R}^p\) is the sensor noise, and \(A \in \mathbb{R}^{n \times n}\), \(B \in \mathbb{R}^{n \times m}\), \(c \in \mathbb{R}^n\), \(C \in \mathbb{R}^{p \times n}\), \(D \in \mathbb{R}^{p \times m}\), \(k \in \mathbb{R}^p\).

  More information can be found in Section 4 in the CORA manual.

  Example

  The code below gives a reachability analysis example for linear systems in CORA. All dynamic systems share a set of identical operations, such as reach to compute the reachable set and simulateRandom for random simulations. Further examples can be found at ./examples/contDynamics on GitHub.

• Linear Systems with Uncertain Parameters

  This class extends linear systems by uncertain parameters. We provide two implementations, one for uncertain parameters that are fixed over time and one for parameters that can arbitrarily vary over time. For the case with fixed parameters, a linear parametric system is defined as

  $$\dot{x}(t) = A(p)~x(t) + B(p)~u(t), ~~ p \in \mathcal{P},$$

  which can be equivalently formulated as

  $$\dot{x}(t) = A x(t) + B u(t), ~~ A \in \mathcal{A},~B \in \mathcal{B},$$ $$\text{with} ~~ \mathcal{A} = \{ A(p) ~|~ p \in \mathcal{P} \},~ \mathcal{B} = \{ B(p) ~|~ p \in \mathcal{P} \},$$

  where \(x(t) \in \mathbb{R}^n\) is the system state, \(u(t)\in \mathbb{R}^m\) is the system input, \(p \in \mathbb{R}^p\) is the parameter vector, and \(\mathcal{P} \subset \mathbb{R}^p\) is the set of parameters. For the case with fixed parameters, a linear parametric system is defined as

  $$ \dot{x}(t) = A(t)~x(t) + B(t)~u(t), ~~ A(t) \in \mathcal{A},~B(t) \in \mathcal{B}.$$

  More information can be found in Section 4 in the CORA manual.

• Linear Discrete-Time Systems

  In addition to continuous-time linear systems, CORA also supports discrete-time linear systems defined as

  $$x[i+1] = A x[i] + B u[i] + c + w[i],$$ $$y[i] = C x[i] + D u[i] + k + v[i],$$

  where \(x[i] \in \mathbb{R}^n\) is the system state, \(u[i] \in \mathbb{R}^m\) is the system input, \(w[i] \in \mathbb{R}^n\) is the disturbance, \(y[i] \in \mathbb{R}^p\) is the system output, \(v[i] \in \mathbb{R}^p\) is the sensor noise, and \(A \in \mathbb{R}^{n \times n}\), \(B \in \mathbb{R}^{n \times m}\), \(c \in \mathbb{R}^n\), \(C \in \mathbb{R}^{p \times n}\), \(D \in \mathbb{R}^{p \times m}\), \(k \in \mathbb{R}^p\).

  More information can be found in Section 4 in the CORA manual.

• Linear Probabilistic Systems

  In contrast to all other systems, we consider stochastic properties here. The system under consideration is defined by the following linear stochastic differential equation (SDE), which is also known as the multivariate Ornstein-Uhlenbeck process:

  $$\dot{x} = A x(t) + u(t) + C \xi(t),$$ $$x(0) \in \mathbb{R}^{n}, \, u(t)\in \mathcal{U} \subset \mathbb{R}^n, \, \xi \in \mathbb{R}^{m},$$

  where \(A\) and \(C\) are matrices of proper dimension and \(A\) has full rank. There are two kinds of inputs: the first input \(u\) is Lipschitz continuous and can take any value in \(\mathcal{U} \subset \mathbb{R}^n\) for which no probability distribution is known. The second input \(\xi\in\mathbb{R}^m\) is white Gaussian noise. The combination of both inputs can be seen as a white Gaussian noise input, where the mean value is unknown within the set \(\mathcal{U}\).

  More information can be found in Section 4 in the CORA manual.

• Linear ARMAX

  Discrete-time Autoregressive Moving Average models with exogenous Input (ARMAX models) are considered in the \(linearARMAX\) class. They can be described by the following equation:

  $$y[i] = \sum_{j=1}^p \bar{A}_{j} y[i-j] + \sum_{j=0}^p \bar{B}_{j} {u}[i-j]$$

  where \([i] \in \mathbb{R}^o\) is the system output at time step \(i\) and \(u[i] \in \mathbb{R}^m\) is the system input which can also include disturbances and measurements noise.

  More information can be found in Section 4 in the CORA manual.

• Nonlinear Systems

  Although a fairly large group of dynamic systems can be described by linear systems, the extension to nonlinear systems is an important step towards the analysis of more complex systems. We consider general nonlinear continuous systems defined by the differential equation

  $$\dot x(t) = f(x(t),u(t)), \label{eq:nonlinearSystem}$$ $$y(t) = g(x(t),u(t)), \label{eq:nonlinearSystem_output}$$

  where \(x(t) \in \mathbb{R}^n\) is the system state, \(u(t) \in \mathbb{R}^m\) is the system input, \(y(t) \in \mathbb{R}^o\) is the system output, and \(f:\mathbb{R}^n \times \mathbb{R}^m \to \mathbb{R}^n\) and \(g:\mathbb{R}^n \times \mathbb{R}^m \to \mathbb{R}^o\) are sufficiently smooth.

  More information can be found in Section 4 in the CORA manual.

  Example

  The code below gives a reachability analysis example for nonlinear systems in CORA. Further examples can be found at ./examples/contDynamics on GitHub.

• Nonlinear Systems with Uncertain Parameters

  Nonlinear parametric systems extend nonlinear systems by additionally considering uncertain parameters \(p\):

  $$\dot{x}(t) = f(x(t),u(t),p), ~~ p\in \mathcal{P} \subset \mathbb{R}^{p}, \label{eq:nonlinParamSystem} $$ $$y(t) = g(x(t),u(t),p) \label{eq:nonlinParamSystem_output},$$

  where \(x(t) \in \mathbb{R}^n\) is the system state, \(u(t) \in \mathbb{R}^m\) is the system input, \(p \in \mathbb{R}^p\) is the parameter vector, \(y(t) \in \mathbb{R}^o\) is the system output, and \(f: \mathbb{R}^n \times \mathbb{R}^m \times \mathbb{R}^p \to \mathbb{R}^n\) and \(g: \mathbb{R}^n \times \mathbb{R}^m \times \mathbb{R}^p \to \mathbb{R}^o\) are sufficiently smooth. As for linear parametric systems, the parameters \(p \in \mathcal{P}\) can be constant over time or time-varying.

  More information can be found in Section 4 in the CORA manual.

• Nonlinear Discrete-Time Systems

  We consider nonlinear discrete-time systems defined as

  $$x[i+1] = f \big( x[i],u[i] \big),$$ $$y[i] = g \big( x[i],u[i] \big),$$

  where \(x[i] \in \mathbb{R}^n\) is the system state, \(u[i] \in \mathbb{R}^m\) is the system input, \(y[i]\) is the system output, and \(f:\mathbb{R}^n \times \mathbb{R}^m \to \mathbb{R}^n\) and \(g:\mathbb{R}^n \times \mathbb{R}^m \to \mathbb{R}^o\) are continuous functions.

  More information can be found in Section 4 in the CORA manual.

• Nonlinear Differential-Algebraic Systems

  We consider time-invariant, semi-explicit, index-1 differential-algebraic systems defined as

  $$\dot{x} = f(x(t),y(t),u(t)), $$ $$0 = g(x(t),y(t),u(t)), $$ $$z = h(x(t),y(t),u(t)), $$

  where \(x(t) \in \mathbb{R}^n\) is the vector of differential variables, \(y(t) \in \mathbb{R}^q\) is the vector of algebraic variables, \(u(t) \in \mathbb{R}^m\) is the vector of inputs, \(z(t) \in \mathbb{R}^o\) is the system output, and \(f: \mathbb{R}^n \times \mathbb{R}^q \times \mathbb{R}^m \to \mathbb{R}^n\), \(g: \mathbb{R}^n \times \mathbb{R}^q \times \mathbb{R}^m \to \mathbb{R}^q\), and \(h: \mathbb{R}^n \times \mathbb{R}^q \times \mathbb{R}^m \to \mathbb{R}^o\) are sufficiently smooth continuous functions. The initial state is consistent when \(g(x(0),y(0),u(0)) = 0\), while for DAEs with an index greater than 1, further hidden algebraic constraints have to be considered. For an implicit DAE, the index-1 property holds if and only if \(\forall t: \, \det(\frac{\partial g(x(t),y(t),u(t))}{\partial y})\neq 0\), i.e., the Jacobian of the algebraic equations is non-singular. Loosely speaking, the index specifies the distance to an ODE (which has index 0) by the number of required time differentiations of the general form \(0 = F(\dot{\tilde{x}}, \tilde{x}, u, t)\) along a solution \(\tilde{x}(t)\), in order to express \(\dot{\tilde{x}}\) as a continuous function of \(\tilde{x}\) and \(t\).

  More information can be found in Section 4 in the CORA manual.

• Neural Network Controlled Systems

  Due to the current trend towards artificial intelligence, the system class of neural network controlled system is gaining more and more importance. In CORA, we consider neural network controlled systems with discrete feedback, where the neural network controller updates the control input at the end of the time steps defined by the sampling time \(\Delta t\). The dynamic system \(\dot{x} = f(x,u)\) can be any of the above dynamic systems.

  More information can be found in Section 4 in the CORA manual or read more here.