Hybrid Automation Extended Model Текст научной статьи по специальности «Математика»

ИНФОРМАТИКА

Hybrid Automation Extended Model

D. K. Andreichenko, K. P. Andreichenko, I. A. Batraeva

Dmitry K. Andreichenko, https://orcid.org/0000-0003-0525-984X, Saratov State University, 83 Astrakhanskaya St., 410012 Saratov, Russia, andre-ichenkodk@gmail.com

Konstantin P. Andreichenko, https://orcid.org/0000-0002-8729-4317, Yuri Gagarin State Technical University of Saratov, 77 Politechnicheskaya St., 410054 Saratov, Russia, kp_andreichenko@renet.ru

Inna A. Batraeva, https://orcid.org/0000-0002-6539-8473, Saratov State University, 83 Astrakhanskaya St., 410012 Saratov, Russia, batraeva-ia@info.sgu.ru

An extended model of hybrid automata for dynamic systems is considered, where, along with a discrete control subsystem and control objects with lumped parameters, there are control objects with distributed parameters (linear and stationary from the point of view of automatic control theory). The possibility of software implementation of an extended model of hybrid automata on embedded computing systems is shown.

Keywords: hybrid automata, hybrid dynamical systems, embedded computing systems.

Received: 21.10.2018 / Accepted: 22.12.2018 Published online: 28.02.2019

DOI: https://doi.org/10.18500/1816-9791 -2019-19-1 -94-104

1 . INTRODUCTION. MODEL OF THE HYBRID AUTOMATA

A hybrid automaton is a mathematical model of a hybrid system, where the switchings of a discrete control subsystem are modeled based on automata theory, and the motion of control objects with lumped parameters is based on Cauchy problems for ordinary differential equations [1-7].

xm

z(t)

y(t)

Fig. 1. Structural scheme

We consider the input-output model of a dynamic system with time t e R, a piecewise continuous input function x(t), __ . R ^ X c RN, an output function y(t), y : R ^ Y and a function z(t), z : R ^ Z, characterizing the time variation of the state of a dynamic system (Fig. 1). The system contains discrete control elements and control objects with lumped parameters interacting with each other.

x

D. K Andreichenho, K P. Andreichenho, I. A Batraewa. Hibrid Automation Extended Model Respectively

Z = S x Zc, |S| < oo, Zc G RNzc, z = (s,zc), s G S, zc G Zc, Y = Yd x Yc, |Yd| < o, Yc G rny, y = (yd,yc), yd G Yd, yc G Yc. (1)

Here S is the finite set of states of the discrete subsystem, Yd is the finite set of values of the discrete outputs, Zc is the set of states of the continuous subsystem, Yc is the set of values of the continuous outputs. The time evolution of a discrete control subsystem is characterized by a sequence of pairs of time intervals and the values of the states of the discrete subsystem that remain on them {([0, ti ],s0), ([ti, t2 ],Si),..., ([tj, tj+i],sj),...}. Therefore

T = {ti, t2, ...,tj, tj+i,...} (2)

is the sequence of switching times. The classic model of the hybrid automaton [1-7] is a tuple

Ah = (Z,X,Y,zo,R, f, fr, fc, f0> (3)

where z0 = (s0, zc0) G Z = S x Zc is the initial state of the system; R = {(sk, Sj)} C S x S is the set of arcs of the oriented graph (S, R) which determines the possible transitions between the states of the discrete subsystem during its switching; fp : r x X x R x Zc ^ {0,1} is a predicate function that allows or prohibits switching of a discrete subsystem. For the sequence of switching times (2) is true

fp(tj, x|tj-o, (Sj-i ,Sj), zc |j-o ) = 1, tj G T, j = 1, 2, 3,... (4)

Function fr : r x X x R x Zc ^ Zc is a function of updating initial conditions at discrete switching

zc\tj +0 = fr (tj , x|ij-0, (sj-i , Sj ), zc\tj-0 ^ j = 1, ^ 3, - (5)

Function fc : r x X x S x Zc ^ Zc is a derivative function defining the right-hand sides of ordinary differential equations, which model the dynamics of a continuous subsystem over time intervals corresponding to persisting states of the discrete subsystem

zc = fc(t, x,S0, zc), zc(0)= zC0, t G [0,ti], () = d()/dt, (6)

:z

fc(t, x, Sj, zc), zc(tj) = zc|t. +0, t G [tj,tj+i], j = 1, 2, 3,...

c "cyj) "c\tj [°j

Here zc|tj+0 and zc|tj-0, j = 1,2, 3,... are connected by equation (5).

f0 : r x X x S x Zc ^ Y is an output function

y = f0(t, x, Sj, zc), j = 1, 2,3,... (7)

In order for the hybrid automaton to be deterministic, it is required that when t ^ 0 the predicate function fp can take the value 1 on no more than one arc of the graph (S, R), originating from the current state of the discrete subsystem. The hybrid automaton is effective software implemented on embedded computing systems. It is of interest to include this class of mathematical models of control objects with distributed parameters [8], and, accordingly, initial-boundary problems for partial differential equations, but with the possibility of efficient software implementation on embedded computing systems.

2 . SPECIFICATION OF THE ORIGINAL MODEL

Depending on the specific type of predicate fp and function fr, at the moment of time tj C T, j = 1,2, 3,... there can be either a single switches of the states of the discrete subsystem or a cascade of switching. However, the switching cascade is undesirable due to energy and structural constraints. We assume that switching cascades are excluded, and in (2)

0 < ti < t2 < ... < tj < tj+i < ... (8)

We assume that the initial state of a dynamical system is an equilibrium state, from which the dynamical system is derived by varying the input vector function with time t ^ 0. Those

x(t) = x0 = const, t < 0; fc(t, x0, s0, zc0) = 0, t ^ 0 (9)

Formally, the sequence of Cauchy problems (6) is equivalent to one Cauchy problem

Zc = f (t, x, Zc), f (t, x, Zc) = (zc\tj +0 - zc\tj-0)5(t - tj) +

ifc(t, x,S0, ZC), t C [0, ti], (10)

\fc(t, x, Sj, Zc), t C [tj ,tj+1 ], j = 1, 2, 3,...

Zc (0)= Zc0 (11)

where 5(t) is the Dirac impulse function. In (10) Zc|t. +0 and Zc|t.-0, j = 1, 2,3,... are related by equation (5). If the continuous dynamic subsystem contains some links with constant or variable positive lag periods, and Tmax is the maximum lag time, the initial conditions (11) take the form

Zc(t)= Zco, t C [-Tmax, 0]. (12)

The mathematical model of the hybrid automaton (1)-(7) or (1)-(5), (7), (8)-(10), (12) is the most feasible mathematical model in the MATLAB / Simulink / Stateflow dynamic system simulation software. MATLAB/ Simulink standard tools provide the generation of program code for the software implementation of this class of dynamic systems on embedded computing systems. For further, it is advisable to separate linear and nonlinear terms as a function of derivatives

fc(t, x, s, Zc) = Bx + CZc + ^1 CkT}Zc(t - Tk) + fc(1)(t, x, S, Zc),

fc(1) : r x X x S X Zc ^ Zc, B = const C r(nzc>Nx), C = const C r(nzc'nzc), Tk = const, 0 <Tk < Tmax, C^) = const C r(nzc'Nzc), k = T,NT.

Accordingly, equation (10) takes the form

Zc = Bx + CZc + ^N= 1 CkT)Zc(t - Tk) + f1 (t, x, Zc), (14)

K1),

(13)

of, ч VI I I чс/, ,ч . I fc (Ь, x, s0, Zc) Ь G J,

fl (i, x, Zc)= ¿J Zc Itj +0 - Zc Itj-0 ЖЬ - j ) + < (1)

J=^2,3,... j f (tx,Sk,zc), t G [ifc,4+1 ], (15)

k = 1, 2, 3,...

In (15) zc|tj +0 and zc|tj-0, j = 1, 2,3,... are connected by the equation (6). 96 Научный отдел

3 . CONTROL OBJECTS WITH DISTRIBUTED PARAMETERS

A dynamic system can contain not only a discrete controlling subsystem and control objects with lumped parameters, but also control objects with distributed parameters, dynamically connected with them across the boundaries. It is a hybrid dynamical system (HDS) [9], and its model equations contain a system of ordinary differential equations and partial differential equations with corresponding initial conditions connected by means of boundary conditions and constraint's conditions. In many cases that are important for applications from the point of view of the automatic control theory, objects with distributed parameters are linear and stationary. In this case, the considered non-linear HDS of a special type correspond to equations (1)-(5), (7)-(9), (15) and equations

zc = Bx + Czc + Ah + CT)zc(t - Tk) + fi(t, x, zc), h = í l(h)udS,

Js

du/dt = l1f)u + bf)x + bf)zc + bf )zc + bf)f2(t,x,zc), r G О,

iG)u + ЬГx + b^zc + bf fa(t, x, zc))

(G),

(G)

(F )

0, S = d О,

zc (t)

t G [-Tmax, 0]; u(r, 0) = uo(r).

(16)

(17)

Here h(t), h : r -

A = const G r(nn},

is the function characterizing the constraint's conditions;

f2,3 (t, x, zc ) =

c2'3) ,(2,3)

(t, x,So, zc), t G [0, ti], (t, x,Sj, zc), t G [tj ,tj+i],

j = 1, 2, 3,

(18)

(2)

: r x

S

, fc(3) : r x

Sr

Nz

3 are nonlinear func-

tions that simulate the effects of switching in the discrete subsystem; r G Q C rnr are independent spatial coordinates of individual points of control objects with distributed

parameters and their area; L(

(F ) 2

b

(F ),

(F )

b(

b(G) b2 r(G) b4

l1f )

r

Nr

c 'Nzc )

G 'Nx) , b

br = b5F)

2 (r), b2F ) (r), bf)

b3G)(r), b3G) :

N ), b

(F) 3,4

b3? (r),

r

Nr

R(n-n2), b2G) = b2G) (r),

b

2

(G)

b4G) (r),

—>

Nr

(F ) 5

, 4g) = (r), L^ : rn ^ r(ng,Nzc) g,n3) are matrices that do not depend on time t (similar to [10]); rn ) ^ (rnr ^ rn ), ) : (rn ^ rn ) ^ (rn ^ rn) are linear unbounded operators which do not depend on time t; l^g) : (rnr ^ rnu ) ^ (rnr ^ rng) is a linear bounded or unbounded operator that does not depend on time t; u = u(r, t), u : rnr x r ^ rnu. We assume that the initial conditions (17) correspond to the equilibrium state, i.e.

x(t) = x0 = const, t < 0, fkc)(t, xo, so, zco) = fkc)(0, xo, so, zco), t < 0, k = 1, 2, 3, Bxo + (c + zco + Aho + fi(c) (0, xo,so, z^ ) = 0, ho = J L(H)uodS, (19)

liF)uo + b2* 4 + b^ Jzco + ff (0, xo,so, zco) = 0, r G О,

(F)

(F)

(F) (c)

iG)uo + brxo + b^zco + bf fr (0, xo, so, zco))

(G)

(G)

(F) (c)

0, S = do.

S

h

2

c

r

zc

->

S

4 . LINEARIZED SYSTEM TRANSFER FUNCTIONS

Consider (16) and (17) under the assumption Zco = 0, u0(r) = 0. The one-sided direct and inverse integral Laplace transform in time t is

Гж ~ 1

/(A) = L [/(t)] = J / (t)e—At dt, / (t) = L—1 [/ (A)] = —

■»a+i^

/ (A)eAt dA,

(20)

a ^ aa, / : r ^ r, / : c ^ c, i = V—1.

Here aa is the abscissa of absolute convergence of the Laplace integral. Denote formally

fk(A) = L[fk(t, x(t), zc(t))], k = 1, 2,3. (21)

Under homogeneous initial conditions zc0 = 0, u0(r) = 0 the equations (16), (17) after Laplace transform take the form

: = BX + (c + ^1 ckT) e-TkA) Zc + Ah + fi,

Azc

Aü = LiF)ü + )X + (L3F) + ALf))ZC + ) f2, r G О,

й + LG)x + Lfzc + ) f3 )

From linearity (23) follows

0, S = d О; h

l(h ) üdS.

(22) (23)

h = Bu (A)x + Cu(A)Zc + C2 (A)f2 + C3 (Af, Bu (A) = Bu (A), Cu,2,3 (A) = C^A).

(24)

Let be e1N) = (1, 0,0,..., 0)T, e2N) = (0,1, 0,..., 0)T,..., e^) = (0,0, 0,..., 1)T. From linearity (23), similarly to [10], linear boundary-value problems follow for finding separate columns of matrices Bu(A), Cu(A), C2(A), C3(A) with fixed values A

Av = liF)v + L2F)ef*), r G О, (l1G)v + l2G)ef*)) Bu(A)ef x) = J l(h)vdS, j = 1Ж,

= 0,

Av = LiF ) v + (L3" + AL^

(f )

(F h(Nzc )

r G О,

,1G) v + L3G) ef^ ) =0, Cu (A)ef-): s j

Av = liF »v + l5f »e<N2 >.

1 5 j

j l(h) vdS, j = 1,Nz

rS

r G О,

(25)

(llG)v) = 0, C2(A)ef2) = f l(h)vdS, j = 1,N, S «/s

Av = l1F)v, r G О,

) = 0, Cs(A)ef3) = f l(h)vdS, j = 17Ж. S Js

<G) v + Lf ejN3)

For small and medium values of A linear boundary value problems (25) can be solved numerically. When |A| > 1 they are asymptotic integration. Substituting (24) into (22) leads to the result

AE - C - ^^ 1 CkT)e-TkA - ACu(A)] Zc = (B + ABU(A))x+

a—«œ

S

S

S

+f1 + AC2(A)f2 + AC3(A)f3, E = diag(1,1,..., 1),

(26)

where it comes from

3

f = $(A)x + J2 k=1 Фк (A)fk.

(27)

In (27) the transfer function matrices are represented as quasi-rational fractions [9,10]

N n — 1

Ф1 (A) = |AE - C - ^k=i ckT)e—Tkл - ACu(A)J = Qi(A)/D(A),

Ф(А) = Ф1 (A)(B + ABu (A)) = Q(A)/D(A), Ф2Й) = Ф^^ (A) = Q2 (A)/D(A), Ф3 (A) = Ф1 (A)AC3 (A) = Q3 (A)/D(A):

D(A) = det ÍAE - C - ^N= 1 Ck )e—TkЛ - ACu(A)

(28)

Here D : c ^ c is the characteristic quasi-polynomial of HDS, Q : c ^ c(nzcN Q1 : c ^ c(Nzc-Nzc), q2 : c ^ c(n-N), Q3 : c ^ c(n-N) are matrices whose elements are perturbing quasi-polynomials of HDS. In [10] the conditions for the ana-lyticity of the functions (A) and Cu(A), as well as the characteristic and perturbing quasi-polynomials D(A) and Q(A) are formulated and proved in the domain |A| > 1, Re A > — | A| sin a, 0 < a < n/2 (conditions with respect to the functions C2(A) and C3(A), as well as perturbing quasi-polynomials Q1 (A), Q2(A) and Q3(A) can be formulated and proved similarly). In particular, with sufficiently general assumptions, this is true for a number of typical mathematical models of control objects with distributed parameters (heat conduction and diffusion processes, damping and supporting layers of a viscous incompressible fluid, the dynamics of an elastic medium with a small but finite energy dissipation).

As a rule, for the most common mathematical models of control objects with distributed parameters

|Bu(A)|| ==o(|A|), ||Cu (A) У ==o(|A|), ||C2 (A) У ==o(|A|), ЦС3 (A)|| ==o(|A|), A —> то, Re A > — |A| sin a, 0 < a < n/2.

(29)

5 . HYBRID AUTOMATA MODEL EXTENSION

Definition 1. By the extended model of a hybrid automaton, we mean a mathematical model in which the switching dynamics of a discrete subsystem is modeled by equations (1)-(5), (7), (12), and the dynamics of a continuous subsystem is modeled by some integral or integro-differential equations.

Lemma 1. If the functions (A), Cu(A), C2(A), C3(A) are analytic with |A| > 1, Re A > — | A| sin a, 0 < a < n/2, and conditions (29) are satisfied, then under homogeneous initial conditions zc0 = 0, u0 (r) = 0 the nonlinear HDS (16), (17) reduces to a nonlinear integral convolution type equation.

Proof. In this case

Ф^) — 0, Фк (A) — 0, k = 1,2, 3, A — то, Re A> — |A| sin a, 0 < a < n/2 (30)

and transfer functions are analytical in domain (30). In this case, for t = 0 and sufficiently large a the Mellin integral (20) for inversion of the Laplace transform converges, and there exist matrices of impulse transition functions

F (t) = L-1 [Ф(Л)], Fk (t) = L-1 [Фк (Л)], ||F(t)|| < ro, IIFfe(t)|| < ro, 0 <t< ro (31)

F(t) = Fk(t) = 0, t< 0, k = 1,2, 3

wherein

lim Fi(t) = lim ЛФ1 (Л) = E, lim i F(e)de = lim ЛФ(Л) = 0,

t^+o л^оо Л (32)

lim i Fk(e)d£ = lim ЛФк(Л) = 0, k = 2, 3.

t^+0 Jo Л^о

Reversing the integral Laplace transform in (27), we find

zc (t) = Г [f (t - e )x(c) + Fk (t - e )fk (e, x(e), zc (e ))] de. (33)

Jo k=1 J

Equation (33) is a non-linear integral equation of convolution type for the unknown function zc(t). □

When 0 < t ^ 1 and when t > 1 for the inversion of the integral Laplace transform in (31) approximate formulas can be obtained on the basis of its asymptotic properties. For medium values of t the inversion of the Laplace transform in (31) can be performed numerically.

Lemma 2. If the initial conditions (17) correspond to the equilibrium state, i.e. conditions (19) are fulfilled, functions (Л), Cu(Л), C2 (Л), C3(Л) are analytic with | Л| > 1, Re Л > — | Л| sin a, 0 < a < n/2, and conditions (29) are satisfied, then nonlinear HDS (16), (17) reduces to a nonlinear integral convolution type equation.

Proof. Putting x*(t) = x(t) — x0, z*(t) = zc(t) — zc0, u*(r,t) = u(r,t) — u0(r), h*(t) = = h(t) — h0, fk*(t,x*,z*) = fk(t,x0 + x*(t), zc0 + z*(t)) — fk(0,x0, zc0), we obtain a system of equations for quantities ()*, similar to (16) with homogeneous initial conditions, whence

it follows z*(t) = /0 [f(t — e)x*(e) + sk=1 Fk(t — ef*(e,x*(e), z*(e»]de. we get

-t

zc(t)= zC0 + / [F(t - f)(x(0 - xo) +

J0

+ Ek=1 Fk(t - 6[ffc(?,x(e), Zc(?)) - fk(0,xo, Zco)] d£ (34)

after returning to the original variables. □

Theorem 1. If the initial conditions (17) correspond to the equilibrium state, i.e. conditions (19) are fulfilled, functions (A), Cu(A), C2 (A), C3(A) are analytic with |A| > 1, Re A > — | A| sin a, 0 < a < n/2, and conditions (29) are fulfilled, then nonlinear special type HDS is isomorphic to an extended model of a hybrid automaton in which the dynamics of a continuous subsystem is modeled by a nonlinear integral convolution type equation.

Proof. On the basis of Lemmas 1 and 2 and using (15) and (33), we find

zc(t)= zco + / [F(t - f)(x(0 - xo) +

Jo

+ £1 Fk(t - 0[fc№)(i,x(?),Sdo, zc(0) - fc(k)(i,xo,Sdo, Zco)]] t e [0, ti],

Er t

m /

. (ZcItj+0 - ZcItj-0)Fi(t - tj) + / [F(t - e)(x(i) - xo)+

j «/ 0

+EL Fk (t - e )fc(k) (e, x(e), ^, zc(e)) - fc(k) (0, xo, ^, zCo)]] de,

(35)

t G [tm,tm+1 ], m = 1, 2, 3,..., tm G T.

The extended model of the hybrid automaton corresponds to equations (1)-(5), (7), (12) and (35), where zc|tj +0 and zc|tj_0 are related by equation (5). □

The generalized degree n G r of the characteristic quasi-polynomial D(A) is determined from the condition [9,10]

lim A_nD(A) = ca, 0 < |ca| < oo, A ^ Re A > ao, ao G (-to, 0) (36)

When conditions (29) are fulfilled, n = Nzc, and condition (36) is fulfilled in a wider domain A ^ to, Re A > — | A| sin a, 0 < a < n/2.

The numerical solution of auxiliary linear boundary value problems (25) with |A| < R, R > 1, Re A > —e, 0 < e ^ 1 is performed on the basis of the Galerkin projection method [10]. In the considered domain, the absence of singularities of functions

(A), Cu(A), C2(A), C3(A) and, consequently, of the characteristic and perturbing quasi-polynomials D(A), Q(A), Q1 (A), Q2(A), Q3(A) is verified numerically on the basis of principle of argument. The presence of singularities of the characteristic and perturbing quasi-polynomials when Re A > 0 means the instability of control objects with distributed parameters as such and requires a change in the parameters of the designed system. Further, we assume that the functions (A), Cu(A), C2(A), C3(A), and also D(A), Q(A), Q1 (A), Q2(A), Q3(A) are analytic in domain Re A > a0, a0 G (—to, 0).

If the quasi-polynomial D(A) is analytic with Re A > a0, a0 G (—to, 0) and is stable, that is, all its roots (zeros) lie strictly to the left of the imaginary axis Re A = 0, then according to the theorem on the stable quasi-polynomial [9]

A arg D(iu) = 1 nn. (37)

2

When optimizing the frequency change step u (37) this is a "fast" algorithm for checking the stability of linearized HDS [9].

Theorem 2. If the initial conditions (17) correspond to the equilibrium state, the functions (A), Cu(A), C2(A), C3(A) are analytic with |A| > 1, Re A > — | A| sin a, 0 < a < n/2, conditions (29) are satisfied, and the characteristic quasi-polynomial D(A) is stable, then the kernel of nonlinear integral equations of convolution type (35) decays exponentially.

Proof. By the conditions of the theorem, the characteristic quasi-polynomial D(A) is analytic in domain Re A > a0, a0 G (—to, 0).

ReX

lHf ß Fig. 2. Integration contour

All the roots (zeros) of D(A) lie strictly to the left of the imaginary axis Re A = 0. Consequently, all the singular points of the transfer functions (28) lie strictly to the left of the imaginary axis Re A = 0. If conditions (29) are true, then from (30) it follows that the singular points of the transfer functions cannot be in the domain Re A > - | A| sin a, 0 < a < n/2, | A| > 1. In formulas (31) the path of integration in the Mellin integral (20) for inversion of the Laplace transform can be along a straight line Re A = —a, a > 0. Further, it can be deformed into the contour U l U (Fig. 2), completely lying in the domain where the transfer functions are analytical. In particular,

F (t) = i. L

Ф(А)ел*dA, l = {A = -a + iu, и G [-R, R]},

-) ulul(+)

l(± = {A = -a ± iR ± ie±ie f, f G [0, то)}, R > 1, 0< в < a < n/2, a > 0,

i Ф(А)ел dA = 1

Ji t

Ï Ф(А)ел^ dA

(Ф(А)е < Me"ai

л*

- JФ'(A)e^dA

sin e df

Л=а-iR Л=-а-гЯ то

M

t sin в

"o(e), t ^то, ==o(e-at), t ^то,

M = max |ф(-а ± iR ± ie±ie£)l < то. Those F(t) = =o(e-at), t ^ то, and similarly Fk(t) = o(e-at), t ^ то, k = 1,2, 3,

a > 0.

□

Corollary 2. Under the conditions of Theorem 2, the asymptotic complexity of the numerical implementation of the extended model of a hybrid automaton based on integral equations (35) has the same order as the asymptotic complexity of the numerical implementation of the model of a hybrid automaton based on Cauchy problems for ordinary differential equations (6).

Remark 1. The numerical integration of linear boundary value problems (25) and the numerical inversion of the integral Laplace transform in (31) can be performed on a sufficiently powerful computing system. Further numerical solution of integral equations (35) can be performed on the embedded computer system.

References

1. Kashevnik A. M., Ponomarev A. V., Savosin S. V. Hybrid Systems Control Based on Smart Space Technology. SPIIRAS Proceedings, 2014, iss. 4(35), pp. 212-226 (in Russian).

2. Meslem N., Ramdani M., Candau Y. Guaranteed Parameter Set Estimation for Monotone Dynamical Systems Using Hybrid Automata. In: Reliable Computing, Springer Verlag, 2010, pp. 88-104.

3. Karoui M. F., Alla H., Chatti A. Monitoring of dynamic processes by rectangular hybrid automata. Nonlinear Analysis: Hybrid Systems, 2010, vol. 4, iss. 4, pp. 766-774. DOI: https://doi.org/10.1016/j-.nahs.2010.05.004

0

4. Thiagarajan P. S., Yang S. Modular discrete time representation of distributed hybrid automata. Theoretical Computer Science, 2012, vol. 429, pp. 292-304. DOI: https://doi.org/10.1016/j.tcs.2011.12.050

5. Konecny M., Taha W., Bartha F. A., Duracz J., Duracz A., Ames A. D. Enclosing the behavior of a hybrid automation up to and beyond a Zeno point. Nonlinear Analysis: Hybrid Systems, 2016, vol. 20, pp. 1-20. DOI: https://doi.org/10.1016/j-.nahs.2015.10.004

6. Elmetennani S., Laleg-Kirati T.M., Djemai M., Tadjine M. New MPPT algorithm for PV applications based on hybrid dynamical approach. Journal of Process Control, 2016, vol. 48, pp. 14-24. DOI: https://doi.org/10.1016/jjprocont.2016.10.001

7. Iovine A., Valentini F., De Santis E., Di Benedetto M. D., Pratesi M. Safe human-inspired mesoscopic hybrid automation for autonomous vehicles. Nonlinear Analysis: Hybrid Systems, 2017, vol. 25, pp. 192-210. DOI: https://doi.org/10.1016/j-.nahs.2016.08.008

8. Shornikov Yu., Bessonov A., Dostovalov D. Specification and instrumental analysis of hybrid systems. Science Bulletin of the Novosibirsk State Technical University, 2015, no. 4(61), pp. 101-117 (in Russian). DOI: http://dx.doi.org/10.17212/1814-1196-2015-4-101-117

9. Andreichenko D. K., Andreichenko K. P. On the theory of hybrid dynamical systems. Journal of Computer and Systems Sciences International, 2000, vol. 39, no. 3, pp. 383398.

10. Portenko M. S., Melnichuk D. V., Andreichenko D. K. Analyticity conditions of characteristic and disturbing quasipolynomials of hybrid dynamical systems. Izv. Saratov Univ. (N. S.), Ser. Math. Mech. Inform., 2016, vol. 16, no. 2, pp. 208-217 (in Russian). DOI: https://doi.org/10.18500/1816-9791-2016-16-2-208-217

Cite this article as:

Andreichenko D. K., Andreichenko K. P., Batraeva I. A. Hibrid Automation Extended Model. Izv. Saratov Univ. (N.S.), Ser. Math. Mech. Inform., 2019, vol. 19, iss. 1, pp. 94-104. DOI: https://doi.org/10.18500/1816-9791-2019-19-1-94-104

УДК 519.713.8: 517.935.2

Расширенная модель гибридного автомата Д. К. Андрейченко, К. П. Андрейченко, И. А. Батраева

Андрейченко Дмитрий Констатинович, доктор физико-математических наук, заведующий кафедрой математического обеспечения вычислительных комплексов и информационных систем, Саратовский национальный исследовательский государственный университет имени Н. Г. Чернышевского, Россия, 410012, Саратов, ул. Астраханская, д. 83, andreichenkodk@gmail.com

Андрейченко Константин Петрович, доктор технических наук, профессор кафедры прикладной математики и системного анализа, Саратовский государственный технический университет имени Гагарина Ю. А, Россия, 410054, Саратов, ул. Политихническая, д. 77, kp_andreichenko@renet.ru

Батраева Инна Александровна, кандидат физико-математических наук, и.о. заведующего кафедрой технологий программирования, Саратовский национальный исследовательский государственный университет имени Н. Г. Чернышевского, Россия, 410012, Саратов, ул. Астраханская, д. 83, batraevaia@info.sgu.ru

Рассмотрена расширенная модель гибридных автоматов для динамических систем, где наряду с дискретной управляющей подсистемой и объектами управления с сосредоточенными по пространству параметрами имеются объекты управления с распределенными по пространству параметрами (линейные и стационарные с точки зрения теории автоматического управления). Показана возможность программной реализации расширенной модели гибридных автоматов на встроенных вычислительных системах.

Ключевые слова: гибридные автоматы, гибридные динамические системы, встроенные вычислительные системы.

Образец для цитирования:

Andreichenko D. K., Andreichenko K. P., Batraeva I. A. Hibrid Automation Extended Model [Андрейченко Д. К., Андрейченко К. П., Батраева И. А. Расширенная модель гибридного автомата] // Изв. Сарат. ун-та. Нов. сер. Сер. Математика. Механика. Информатика. 2019. Т. 19, вып. 1. С. 94-104. DOI: https://doi.org/10.18500/1816-9791-2019-19-1-94-104