Behavior of DEVS


Behavior of DEVS

Behaviors of a given DEVS model is a set of sequences of timed events including null events, called event segments which make the model move one state to another within a set of legal states. To define this way, the concept of a set of illegal state as well a set of legal states are needed to be introduced.

In addition, since the behaviors of a given DEVS model needs to define how the state transition change both when time is passed by and when an event occurs, it has been described by a much general formalism, called general system [ZPK00]. In this article, we use a sub-class of General System formalism, called timed event system instead.

Depending on how to define the total state and its external state transition function of DEVS, two ways to define the behavior of a DEVS model using Timed Event System. Since behavior of a coupled DEVS model is defined as an atomic DEVS model, behavior of coupled DEVS class is defined by timed event system.

Contents

View 1: total states = states * elapsed times

Suppose that a DEVS model, \mathcal{M}=<X,Y,S,s_0,ta,\delta_{ext},\delta_{int},\lambda> has

  1. the total state set Q=\{(s,t_e)| s \in S, t_e \in (\mathbb{T} \cap [0, ta(s)])\} where te denotes elapsed time since last event and  \mathbb{T}=[0,\infty) denotes the set of non-negative real numbers, and
  2. the external state transition  \delta_{ext}:Q \times X \rightarrow S.

Then the DEVS model, \mathcal{M} is a Timed Event System \mathcal{G}=<Z,Q,q_0,Q_A,\Delta> where

  • The event set Z=X \cup Y^\phi.
  • The state set Q=\mathcal{M}.Q \cup \{(illegal,t_e)| illegal \not \in S, t_e \in \mathbb{T} \}.
  • The initial state  \,q_0 = (s_0,0).
  • The set of acceptance states  Q_A = \mathcal{M}.Q.
  • The state trajectory function  \Delta: Q \times \Omega_{Z,[t_l,t_u]}
\rightarrow Q is defined for an total state  q=(s,t_e) \in Q at time  t \in \mathbb{T} and an event segment  \omega \in \Omega_{Z,[t_l,t_u]} as follows.
If unit event segment ω is the null event segment, i.e.  \omega=\epsilon_{[t, t+dt]}
\, \Delta(q, \omega)=(s, t_e+dt).\,

If unit event segment ω is a timed event ω = (x,t) where the event is an input event  x \in X,

 
\Delta(q, \omega)= 
\begin{cases}
(\delta_{ext}(s,t_e,x), 0)& \textrm{if } ~s \in S\\
(illegal, t_e)& \textrm{otherwise}.
\end{cases}

If unit event segment ω is a timed event ω = (y,t) where the event is an output event or the unobservable event  y \in Y^\phi,

 
\Delta(q, \omega)= 
\begin{cases}
(\delta_{int}(s), 0)& \textrm{if } ~s \in S, t_e = ta(s), y = \lambda(s)\\
(illegal, t_e)& \textrm{otherwise}.
\end{cases}

Computer algorithms to simulate this view of behavior are available at Simulation Algorithms for Atomic DEVS.

View 2: total states = states * lifespans * elapsed times

Suppose that a DEVS model, \mathcal{M}=<X,Y,S,s_0,ta,\delta_{ext},\delta_{int},\lambda> has

  1. the total state set Q=\{(s,t_s, t_e)| s \in S, t_s\in \mathbb{T}^\infty, t_e \in (\mathbb{T} \cap [0, t_s])\} where ts denotes lifespan of state s, te denotes elapsed time since last tsupdate, and  \mathbb{T}^\infty=[0,\infty) \cup \{ \infty \} denotes the set of non-negative real numbers plus infinity,
  2. the external state transition is  \delta_{ext}:Q \times X \rightarrow S \times \{0,1\} .

Then the DEVS Q=\mathcal{D} is a timed event system \mathcal{G}=<Z,Q,q_0,Q_A,\Delta> where

  • The event set Z=X \cup Y^\phi.
  • The state set Q=\mathcal{M}.Q \cup
\{(illegal,\infty, t_e)| illegal \not \in S, t_e \in \mathbb{T} \}.
  • The initial state  \,q_0 = (s_0,ta(s_0),0).
  • The set of acceptance states  Q_A = \mathcal{M}.Q.
  • The state trajectory function  \Delta: Q \times \Omega_{Z,[t_l,t_u]}
\rightarrow Q is defined for a total state  q=(s,t_s, t_e) \in Q at time  t \in \mathbb{T} and an event segment  \omega \in \Omega_{Z,[t_l,t_u]} as follows.
If unit event segment ω is the null event segment, i.e.  \omega=\epsilon_{[t, t+dt]}
 \,\Delta(q, \omega)=(s, t_s, t_e+dt).\,

If unit event segment ω is a timed event ω = (x,t) where the event is an input event  x \in X,

 
\Delta(q, \omega)= 
\begin{cases}
(s', ta(s'), 0)& \textrm{if } ~s \in S, \delta_{ext}(s,t_s,t_e,x)=(s',1)\\
(s', t_s, t_e)& \textrm{if } ~s \in S, \delta_{ext}(s,t_s,t_e,x)=(s',0)\\
(s, t_s, t_e)& \textrm{otherwise}.
\end{cases}

If unit event segment ω is a timed event ω = (y,t) where the event is an output event or the unobservable event  y \in Y^\phi,

 
\Delta(q, \omega)= 
\begin{cases}
(s', ta(s'),0)& \textrm{if } ~s \in S, t_e = t_s, y = \lambda(s), \delta_{int}(s)=s'\\
(illegal, \infty, t_e)& \textrm{otherwise}.
\end{cases}

Computer algorithms to simulate this view of behavior are available at Simulation Algorithms for Atomic DEVS.

Comparison of View1 and View2

Features of View1

View1 has been introduced by Zeigler [Zeigler84] in which given a total state  q=(s,t_e) \in Q and

\, ta(s)=\sigma

where σ is the remaining time [Zeigler84] [ZPK00]. In other words, the set of partial states is indeed S=\{(d,\sigma)| d \in S', \sigma \in \mathbb{T}^\infty \} where S' is a state set.

When a DEVS model receives an input event  x \in X, View1 resets the elapsed time te by zero, if the DEVS model needs to ignore x in terms of the lifepan control, modellers have to update the remaining time

\, \sigma = \sigma - t_e

in the external state transition function δext that is the responsibility of the modellers.

Since the number of possible values of σ is the same as the number of possible input events coming to the DEVS model, that is unlimited. As a result, the number of states  s=(d, \sigma) \in S is also unlimited that is the reason why View2 has been proposed.

If we don't care the finite-vertex reachability graph of a DEVS model, View1 has an advantage of simplicity for treating the elapsed time te = 0 every time any input event arrives into the DEVS model. But disadvantage might be modelers of DEVS should know how to manage σ as above, which is not explicitly explained in δext itself but in Δ.

Features of View2

View2 has been introduced by Hwang and Zeigler[HZ06][HZ07] in which given a total state  q=(s, t_s, t_e) \in Q , the remaining time, σ is computed as

\, \sigma = t_s - t_e.

When a DEVS model receives an input event  x \in X, View2 resets the elapsed time te by zero only if δext(q,x) = (s',1). If the DEVS model needs to ignore x in terms of the lifepan control, modellers can use δext(q,x) = (s',0).

Unlike View1, since the remaining time σ is not component of S in nature, if the number of states, i.e. | S | is finite, we can draw a finite-vertex (as well as edge) state-transition diagram [HZ06][HZ07]. As a result, we can abstract behavior of such a DEVS-class network, for example SP-DEVS and FD-DEVS, as a finite-vertex graph, called reachability graph [HZ06][HZ07].

See also

  • Behavior of Coupled DEVS
  • Simulation Algorithms for Atomic DEVS
  • Simulation Algorithms for Coupled DEVS

References


Wikimedia Foundation. 2010.