<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.0 20120330//EN" "JATS-archivearticle1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink">
  <front>
    <journal-meta />
    <article-meta>
      <title-group>
        <article-title>Research of Pareto-Optimal Schemes of Control of Availability of the Information System for Critical Use</article-title>
      </title-group>
      <pub-date>
        <year>1943</year>
      </pub-date>
      <fpage>0000</fpage>
      <lpage>0002</lpage>
      <abstract>
        <p>The relevance of the task of control the availability of information systems in general and of information systems for critical use (ISCUs) in particular is increasing simultaneously with the information needs of humanity. However given the importance of an information resources of ISCU introducing new schemes of control the availability without a simulation stage can lead to significant material and in particular human losses. These circumstances led to the expediency of research the scientific and applied results of which are presented in the article. In the article the new mathematical models for identification of the Pareto-optimal schemes of control of accessibility of ISCU with a parallel functioning set of input information ports are synthesized. The ISCU is described by the composition of single-line Markov or semi-Mark queuing systems in which unlike existing ones availability loss functions are formulated with a focus on the ability of quantitative or qualitative estimation of the scheme of control of availability and guarantee the processing of access requests with different priority degree which allows to set and solve the multicriteria optimization task for an identification of problem-oriented scheme of control of availability of ISCU taking into account its architectural features.</p>
      </abstract>
      <kwd-group>
        <kwd>information system for critical use</kwd>
        <kwd>mathematical models</kwd>
        <kwd>availability</kwd>
        <kwd>access control</kwd>
        <kwd>optimization</kwd>
        <kwd>mathematical programming</kwd>
        <kwd>operations research</kwd>
      </kwd-group>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>-</title>
      <p>Process modeling in modern information systems remains actual and evolves in sync
with the development of the information society. Such attention is driven by the
constant growth of both the material value of information and its direct impact on
people's lives, which is especially evident in information systems for critical use (ISCUs).
Applied mathematical models allow to predict the results of current and service
operations with information in such systems, to optimize information processes, to
rationally control information systems in conditions of active person-system interaction,
to design mechanisms of information security etc. However, the complex nature of
information systems makes it essential to limit the scope of mathematical models
created to prevent the complexity of the latter. However, for a complete description of
the information systems under study, the mathematical models created should be
generalized by a certain interconnected hierarchical structure, the role of which the
dependability concept [1-5] has played in critical systems for decades.</p>
      <p>Naturally, the integration of applied mathematical models of information systems
dependability attributes is simplified if the first are built in the paradigm of a single
mathematical apparatus. In particular, the authors have already explored [6-9]
mathematical models of confidentiality, accessibility, integrity, reliability of information
systems for critical use in the paradigm of Markov process theory. If we concentrate
on modeling the availability of ISCU, the access process to the information resources
of the latter, taking into account its architectural features and the technologies used to
ensure confidentiality, is described in [6, 7, 10], but the question of the synthesis of
optimal availability control schemes remained open, which led to the thematic focus
of this article. Accordingly, the object of the study will be considered the process of
synthesis of the target optimal schemes for controlling the availability of ISCU, and
the subject of the study will be methods of Markov theory, queuing theory and
operations research theory.
2</p>
    </sec>
    <sec id="sec-2">
      <title>Review of related works</title>
      <p>Mathematical modeling of optimization tasks has its own specificity and generally
includes the step of identifying a set of factors from which a subset of managed
variables is distinguished, generalizations of the latter into the objective function and a
constraint system whose expressions are compared with free members that are
determined by boundary values of the unmanaged variables. For applied optimization
tasks, which may include the synthesis and analysis of optimal ISCU availability
schemes, the optimization task statement is implemented individually, within the
chosen mathematical apparatus, the choice of which largely determines the adequacy
of the model and the method of solution. For the tasks of optimization of processes of
functioning of real technical, information, economic telecommunication network
systems, in particular, mathematical apparatus of queuing theory is used [11-17]. A
typical queuing system (QS) includes a finite number of servicing devices or service
channels. In addition to the service channels, a specific QS is identified by the values
of such characteristics as flow of requirements, function of distribution of duration of
service, queue parameters etc. QS optimization is to identify a configuration scheme,
the application of which provides the extreme value of a formalized optimality
criterion.</p>
      <p>In classical queuing theory, no intervention is assumed in the functioning of the
system under study, that is, external control effects aren’t taken into account. When
describing QSs that control the impact on the performance of which can’t be ignored,
use the derived theory of managed queuing systems. Methods of managed QS theory
allow us to set the tasks of optimal management of queue, service channels, duration
of service, input flow of requirements etc. The flexibility of the mathematical
apparatus of queuing theory allows us to describe the information flows of QS using a
powerful mathematical apparatus of theory of Markov processes [18-21] and to
identify the optimal control schemes, apply the appropriate methods of the theory of
operations research [22, 23]. This arrangement allows to provide adequate mathematical
description of large-scale real information systems [15, 17-19, 21]. However, it
should be noted that in many models of description of applied information systems in
the QS paradigm [12-17] due consideration isn’t given to the simultaneous
management of several characteristics, which is natural for ISCU in which, for example, it is
necessary to combine a high degree of confidentiality with a high characteristic
availability. Therefore, the aim of the study is to synthesize and analyze the
problemoriented optimal schemes of multicriteria accessibility management of the information
system for critical use modelled in the QS paradigm, taking into account the
architectural features of the latter.
3</p>
    </sec>
    <sec id="sec-3">
      <title>Problem statement</title>
      <p>
        Communication of multiple authorized persons (APs) I = {i} , i = 1, n , with the
information environment of ISCU is implemented through a set of open information
ports J = { j} , j = 1, m . Technologically person-system interaction is organized in
the form of closed sessions. If at a certain time moment the j th information port is
occupied by the session of the і th AP then until its completion this port is closed to
the rest of the APs and their requests addressed to this port are ordered in a queue.
We’ll consider queries sent from the і th AP to the j th information port of ISCU as
an independent Poisson process with intensity λі and characterize by the probability
of satisfaction of the request xij . The duration of the session during which the j th
port is occupied by the і th AP will be considered by stochastic value with
exponential distribution with indicator µ j . Given the described input data to model and study
the person-system interaction we use the mathematical apparatus of queuing theory
[11-17] in the paradigm of which we’ll set the problems of identification of optimal
schemes of control of availability of ISCU, which will differ in the formalization of
the personified indicator of efficiency of control of availability of ISCU – loss
function Li ( х) which characterizes the average number of denials of access to the і th AP
for a certain period of time. Taking into account the competitive nature of access to
the information environment of ISCU we formalize the process of choosing the most
optimal scheme from a set of ISCU accessibility management schemes in the
paradigm of the mathematical apparatus of game theory [22, 23]. We present a
generalized objective function of such an optimization task with a tuple
where І is a set of APs indexes, X i
of an access organization schemes for an і th AP, Li ( х) is a loss function for an і th
AP, x =x1,, ( xn ) ∈ X i . We classify the optimization task generalized by tuple (
        <xref ref-type="bibr" rid="ref37">1</xref>
        ) as
a multi-coalition game, the optimal solution of which is the equilibrium point
x∗ = ( x1∗ ,, xn∗ ) for which such an inequality holds:
Γ = І ,{X i }i∈І ,{Li }i∈І ,
      </p>
      <p>
=xi
</p>
      <p>m
=xi1,, ( xim ) : xij ∈ [0,1], ∑ xij
j=1</p>
      <p>
=1is a set
</p>
      <p>
        Li ( x1∗ ,, xi∗ , xn∗ ) ≤ Li ( x1,, xi ,, xn ) ,
∀i ∈ A, xi ∈ X i and functions Li ( хі ) are continuous and differentiated by хі .
Therefore, the objectives of the study are to formulate rational variants of loss functions
Li ( х) , to formalize for them the optimization tasks of the form (
        <xref ref-type="bibr" rid="ref37">1</xref>
        ), synthesis of
methods for solving sach the optimization tasks and interpretation of the obtained
optimal solutions – equilibrium points for the organization of schemes of the
availability to the target ISCU.
(
        <xref ref-type="bibr" rid="ref37">1</xref>
        )
(2)
(3)
4
      </p>
    </sec>
    <sec id="sec-4">
      <title>Models and methods</title>
      <p>The mathematical apparatus of game theory for the study of the task of finding the
optimal schemes for the availability of ISCU is chosen in particular for the potential
property of its optimal solution as Pareto optimality [24-26] which is that the
equilibrium point found for the input data indicate an optimal scheme of control of
availability of ISCU which not only improves the access conditions for the part of the APs, but
also ensures that the access conditions of the other APs don’t deteriorate which is
especially important given the industry purpose of ISCU. Therefore, we formulate
loss functions so that the solutions synthesized on the basis of their optimization tasks
are Pareto optimal. With a look at the mathematical apparatus of queuing theory we
formalize the loss function Li ( х) by an expression</p>
      <p>m
Li ( х) = λi ∑ xij Pj ( х) ,</p>
      <p>
        j=1
where х ∈ Х , i = 1, n , xi = ( xi1,, xim ) , Pj ( х) is the stationary probabilities that the
j th port is busy at the time of accessing the і th AP. Applying to the analyzed model
of the availability of ISCU with a parallel operating set of input information ports (
        <xref ref-type="bibr" rid="ref37">1</xref>
        )
the method of polynomial approximation [12-16] we decompose it into m one-line
Markov QS [17] whose stationary probabilities Pj ( x) are described by an
expression
      </p>
      <p>m
where Λ j =∑λi xij , i = 1, n . Let’s generalize expressions (3), (4) by analytically
j=1
formalizing the task of finding a Pareto-optimal scheme of availability of the ISCU
with the objective function of the form</p>
      <p>n
min ∑λi Li ( x)
x∈X i =1</p>
      <p>n m
=mx∈iXn∑i=1λi ∑j=1xij Λ j (µ j + Λ j )
−1</p>
      <p>n
=∑λiLi ( x∗ ) ,
i =1
where х∗ =x1∗,, ( xn∗ ) ∈ Х is a Pareto-optimal scheme of the control of availability of
the ISCU for which the following conditions are fulfilled:
∂Li ( x1∗ ,, xi∗−1, xi∗ , xi∗+1,, xn∗ )
∂xij
x∗j
= 0 holds;
- for all i = 1, n , j = 1, m , the equality
m
- for all xij ≥ 0 , i = 1, n , j = 1, m , the equality ∑ xij = 1 holds;
j=1
−1
Pj ( x) = Λ j (µ j + Λ j ) ,</p>
      <p>∂xij∂xik
 ∂Zi ( x1∗ ,, xi∗−1, xi∗ , xi∗+1,, xn∗ )

 ∂xij

 ∂Zi ( x1∗ ,, xi∗−1, xi∗ , xi∗+1,, xn∗ )

 ∂β i
= 0,
- the Hessian matrix of functions Li ( x1∗ ,, xi∗−1, xi∗ , xi∗+1,, xn∗ ) is positive defined
whereby the derivatives ∂Li ( x1∗ ,, xi∗−1, xi∗ , xi∗+1,, xn∗ ) , i = 1, n , j, k = 1, m , j ≠ k ,
m
that form the matrix are calculated provided that ∑ xij = 1 .
j=1</p>
      <p>We formalize the process of solving the optimization task described above with
the objective function (5) and the corresponding constraint conditions using the
Lagrange multiplier method [26].</p>
      <p>Let’s formalize the process of solving the optimization task described above with
the objective function (5) and the corresponding constraint conditions using the
Lagrange multiplier method [26]. We formulate Lagrange functions of the form
 
Zi ( x1,, xn ) =Li ( x1,, xn ) − β i  ∑ xij −1 , i = 1, n , j = 1, m ,</p>
      <p> j 
where β і are Lagrange multipliers. Then, given (6) the constraint conditions of the
optimization task (5) are generalized by the system of equations of the form</p>
      <p>For the optimization task under study with objective function (5) the system (7)
with controllable variables β і , xij is formed by nonlinear equations of the form
(Λ2j + µ j Λ j + xijλiµ j ) (µ j + Λ j )
−2
(nΛ2j + (n +1) Λ jµ j ) (Λ j + µ j )
−2</p>
      <p>n
=∑ β iλi−1 , j = 1, m .</p>
      <p>i=1</p>
      <p>
        Given that the left-hand side of the equations of form (9) doesn’t contain control
variables it can be represented as
(nΛ12 + (n +1) Λ1µ1 ) (Λ1 + µ1 )−2 = (nΛ2j + (n +1) Λ jµ j ) (Λ j + µ j )−2 , j = 2, m . (
        <xref ref-type="bibr" rid="ref21 ref26 ref29 ref3 ref7">10</xref>
        )
Developing the idea underlying expression (
        <xref ref-type="bibr" rid="ref21 ref26 ref29 ref3 ref7">10</xref>
        ) we present equation (9) like
C0 j Λ2j + C1 j Λ j + C2 j =0 , j = 2, m ,
(11)
where C0 j = nµ12 + (n −1) Λ1µ1 , C1 j = (n +1)µ12µ j − (n −1) Λ12µ j , C1 j = (n +1) Λ1µ1 −
−nΛ12µ 2j are the corresponding coefficients. The analytic solution of system (7)
formed by the equations of the form (11) for Λ j ≥ 0∀j will be an expression:
∑λi =∑ Λ j =Λ1 − 0.5∑ C0−1j (C1 j − C12j − 4С0 jC2 j ) =Λ1 + Λ1µ1−1 ∑m µ j ,
n m m
i =1 j =1 j =2 j =2
from which we get
n  m −1
Λ1 =Λ2 = =Λm =µ1 ∑λi  ∑ µ j  .
      </p>
      <p>i =1 j =1
m
Given that ∑ xij −1 =0 we express xij from expression (8) and obtain a system
j=1
of equations of the form
β iλi−2 ∑m µ −j1 (µ j + Λ j )2 − λi−1 ∑m µ −j1 (µ j + Λ j ) Λ j =1 , i = 1, n ,</p>
      <p>j =1 j =1
Analytical solution of equation (13) with respect to expression (12) is formalized as
(8)
(9)
(12)
(13)</p>
      <p>The values of the objective functions Li ( x ) at the point xi∗j are calculated by the
expression</p>
      <p>β i = λi2µ (λ + µ )−2 + λiλ (λ + µ ) , i = 1, n ,
n m
where λ = ∑λi , µ = ∑ µ j . Expressions (12) and (14) can be generalized by
substii=1 j=1
tuting them into expression (13). As a result we obtain an expression to identify the
equilibrium point for the optimization task with the objective function (5):
(14)
(15)
(16)</p>
      <p>Assume that the constraint conditions of the optimization task with the objective
function (5) for the equilibrium point (15) are satisfied. Substitute in the second
derivative of the functions Li ( x ) the coordinates of the equilibrium point xi∗j calculated
by
the
expression
(15):
∂2 Li ( x )
∂xi2j
= µ j Λ j (µ j + Λ j − λi xij ) (µ j + Λ j )
−3
=
=λµ −1 (µ + λ − λi ) (µ + λ )−1 &gt; 0 . Accordingly for other points:
j, k = 1, m , j ≠ k . The Hessian matrix H ( Li ( x )) is positive defined and the rest of
the constraint conditions formulated in the setting of the optimization task with the
objective function (5) are satisfied. We also prove that the optimal solution obtained
in the form (15) is also Pareto optimal. According to the statement of task (5):
m n  m  n
xi1 = 1 − ∑ xij , i = 1, n , then Λ1 =∑ 1 − ∑ xij λi , Λ j =∑ xijλi , j = 2, m ;
j=2 n = 1 j =2 i=1</p>
      <p> m  m −1 
Li ( x1,, xn ) =λi  Λ1 (µ1 + Λ1 )−1 1 − ∑ xij  + ∑ Λ j (µ j + Λ j ) xij  , i = 1, n ,
respec  j =2 j =2 
∂2 Li ( x )
∂xij ∂xik
= 0 , i = 1, n ,
tively,
∂L ( x )
∂xij</p>
      <p>n  m −1 
+β ∑ λi ∑ Λ j (µ j + Λ j ) xij  . We identify the extremum of the function L ( x ) :
i =1 j =2 
= −βλi ( Λ1 ( 2µ1 + Λ1 ) (µ1 + Λ1 )−2 − Λ j ( 2µ j + Λ j ) (µ j + Λ j )−2 ) = 0 ,
i = 1, n ,
j = 2, m , why simplify the latter: µ12Λ12 + 2µ12µ j Λ j − µ 2jΛ1 ( 2µ1 + Λ ) =0 and solve
1
the
resulting
expression
with
respect
to
=−1 µ1λµ . On the basis of the expression for Λ1 we determine
n
Λ j : Λ j = µ j Λ1µ1−1 = µ jλµ −1 and form a system of linear equations ∑λi xij = Λ j , an
i=1
analytic view of the solution whose, xij = µ jµ −1 , i = 1, n , j = 1, m , is identical to the
expression (15). Therefore formalized by expression (15) the optimal solution of the
optimization task with the objective function (5) is Pareto-optimal.</p>
      <p>
        Applying to the analyzed model of the availability of ISCU with a parallel
operating set of input information ports (
        <xref ref-type="bibr" rid="ref37">1</xref>
        ) the method of polynomial approximation we
decompose it into m one-line semi-Markov QS. In this case the stationary
probabilities Pj ( x) defined by expression (4) will be redefined as:
−1
      </p>
      <p>Pj ( x) = Λ jτ j (1+ Λ jτ j ) ,
m
where Λ j =∑λi xij are the intensity of incoming requests from the APs, i = 1, n , τ j
j=1
is the instantaneous value of the duration of the incoming request service for the j th
open port of ISCU - a stochastic variable with arbitrary distribution. Accordingly the
analytic appearance of the objective function (5) will change:</p>
      <p>n n
mx∈iXn ∑i1 =Λjτ λi Li ( x) =mx∈iXn j (1+ Λ jτ j )−1 =∑i1 =.(18) λi Li ( x∗∗ )
(17)
(19)</p>
      <p>The constraint conditions formulated for the optimization task with the objective
function (5) remain relevant also for the optimization task with the objective function
(18). Let’s carry out the analytical formalization of the process of identification of the
equilibrium point х∗∗ =x1∗∗,, ( xn∗∗ ) ∈ Х which is the Pareto-optimal scheme for the
control of the availability of ISCU whose loss function Li ( х) given (17) is described
by the expression</p>
      <p>m −1
Li ( х) =λi ∑ xij Λ jτ j (1+ Λ jτ j ) , i = 1, n .</p>
      <p>j=1</p>
      <p>As in the previous study we use the Lagrange multiplier method to describe the
process of solving the optimization task with the objective function (18). On the basis
of the Lagrange function (6) we generalize the constraint conditions for the
optimization task with the objective function (18) with respect to the controlled variables β і ,
xij by a system of nonlinear equations of the form
(Λ2jτ 2j + Λ jτ j + λiτ j xij ) (1+ Λ jτ j )
−2
=λi−1β i
m
at ∑ xij −1 =0 , i = 1, n . Simplifying the system formed by the equations of the form
j=1
(20) we consistently obtain analytical expressions to calculate the values of Λ1 , xi∗j :
 m −1 n
Λ1 =τ 1 ∑τ −j1  ∑λi ,</p>
      <p> j =1 i =1
xi∗j = τ j ∑τ −1 
 m</p>
      <p>j  .
 j=1 </p>
      <p>−1
∂2 Li ( x∗∗ )
(21)
(22)
(23)
matrix H ( Li ( x)) at the point xi∗j∗ is positive defined. Therefore, the point xi∗j∗
calculated by expression (22) is an equilibrium point and the values of the loss functions at
n  n m 
this point are determined by the expression Li ( x∗∗ ) = =∑λi λi  ∑λi + ∑τ j  ,
i =1 i =1 j =1
i = 1, n . The procedure for proving the Pareto-optimality of the equilibrium point xi∗j
calculated by expression (15) is also valid for the point xi∗j∗ calculated by expression
(22).</p>
      <p>In addition to the variants formalized by expressions (3), (19) the formulation of
loss functions Li ( х) with an orientation on a quantitative estimation of the optimal
scheme for the control of the availability of ISCU is of practical importance. We
formalize the loss functions in the form of such an expression:</p>
      <p>m
Li ( х) = λi ∑ xij K jν j ,</p>
      <p>
        j=1
ν j = Λ j (µ j (µ j − Λ j ))−1
where K j is the cost per unit of waiting time for the release of j th information port
of the ISCU and a parameter
describes the average waiting time for the release of this information port. For
functions Li ( х) defined by expression (23) expression (8) will look like this:
( K j (µ j (Λ j + λi xij ) − Λ2j )) (µ j (µ j − Λ j ) )
2 −1
we sum up equation (25) by i = 1, n and obtain quadratic equation with respect to Λ j :
(n + Cm ) Λ2j − µ j (n + 2Cm +1) Λ j + Cmµ 2j =0 ,
where Cm =(µ mΛm (n +1) − Λ2mn) (µ m − Λm )−2 . Solving equation (27) we obtain the
value of Λ j substituting which in equation (25) we get the expression to calculate the
coordinates of the equilibrium point xi∗j :
(25)
(
        <xref ref-type="bibr" rid="ref36">26</xref>
        )
(27)
(28)
      </p>
      <p>The value of the loss functions at the equilibrium point xi∗j is analytically
formalized by the expression</p>
      <p>Li ( xi∗j ) =λiλµ −1 ∑m K j (µ − λ )−1 =λiλ (µ (µ − λ ))−1 ∑m K j ,</p>
      <p>
        j =1 j =1
n m
where λ = ∑λi , µ = ∑ µ j . The procedure for proving the Pareto-optimality of the
i=1 j=1
equilibrium point for the task of finding the optimal scheme for the control of
accessibility of ISCU with the loss functions of the species (23) and the objective function
of the species (5) is methodologically identical to that described for expression (15)
and allows us to claim that the equilibrium point (28) is a Pareto-optimal if condition
(
        <xref ref-type="bibr" rid="ref36">26</xref>
        ) is satisfied. Condition (
        <xref ref-type="bibr" rid="ref36">26</xref>
        ) is rational since it governs the proportional
relationship between the cost of waiting for the AP to gain access and the speed of service of
incoming requests of the ISCU.
      </p>
      <p>Finally in order to handle emergencies we mathematically describe the precedent
that a і th AU with a probability ωі or 1−ωі can qualify its request as a priority or
ordinary respectively. Consider this information in the corresponding loss function:
Li (ω1,,ωn ) = Ki1λi M {γ 1}ωi + Ki2λi M {γ 2}(1−ωi ) ,
(29)
where M {γ 1} = M {γ 1 (ω1,,ωn )} and M {γ 2} = M {γ 2 (ω1,,ωn )} are the average
delay in satisfaction of priority and ordinary requests for access, arranged in the
corresponding queues γ 1 and γ 2 the unit cost of the request in which for the і th AP are
Ki1
and</p>
      <p>Ki2
respectively.</p>
      <p>Therefore, two
queues
Λ2 =∑λі (1−ωі ) are formed on the open information port of ISCU for priority and
ordinary requests respectively. We extend this concept to m open ports of ISCU the
parameter of session duration for which is a stochastic quantity with exponential
distribution with indicator µ and the queues γ 1 and γ 2 are characterized by the
inequalities (Λ1 + Λ2 ) (mµ )−1 &lt; 1 and Λ1 (mµ )−1 &lt; 1 respectively. As a result we define
Λ1 =∑λіωі
and
the parameters M {γ 1} and M {γ 2} by expressions M {γ 1}
and</p>
      <p>M {γ 2} = С (m) ((1− ρ m−1 ) ×</p>
      <p>  −1
×1− m−1 ∑n ρiωi   ,</p>
      <p> i=1  
  −1
C (m) = mµ 1+ (m −1)!(m − ρ ) ρ −m ∑m−1 ((k !)−1 ρ k )   ,</p>
      <p>  k =0  
Let’s redefine the loss function (29) taking into account the above:
n
ρ =∑ ρi =µ −1 (Λ1 + Λ2 ) .</p>
      <p>i=1
Li (ω1,,ωn ) =λiС (m) ( Ki1ωi + Ki2 (1−ωi )) ×
×1− m−1 ∑n ρiωi  (1− m−1ρ )−1 , i = 1, n .</p>
      <p> i=1 </p>
      <p>The optimization task is to identify the point ω ∗ = (ω1∗ ,,ωi∗ ,,ωn∗ ) for which
condition Li (ω ∗ ) ≤ Li (ω ) ∀i ∈ I ,ω ∈ Ω is satisfied. To find the coordinates of a point
ω ∗ we use the necessary conditions for finding the extremum of a convex function:
 −1
=С(m) 1− m−1 ∑n ρiωi </p>
      <p> i=1 
where
ρi = λiµ −1 ,
(30)
(31)
 ∂Li (ω1,,ωi ,,ωn ) = 0,

 ∂ωi
 ∂2 Li (ω1,,ωi ,,ωn ) &gt; 0,
 ∂ωi2
i = 1, n .</p>
      <p>The first condition (31) for the task of identifying the optimal scheme for control
of availability of ISCU with loss functions (30) and the objective function (5) is
generalized by a system of linear equations with respect to ωі :
m−1 ∑n ρ kωk = ( Ki1 (1− m−1ρ ) − Ki2 (1− m−1ρi )) ( Ki1 (1− m−1ρ ) − Ki2 ) , i = 1, n , (32)
k =1
k ≠i
solving which we obtain the set of admissible schemes for control of availability of
ISCU, investigating which using the second condition of (31) we identify the optimal
scheme:
xi∗ =m ((n −1) ρi )
−1 n
∑ (rk − (n − 2) ri ) =m ((n −1) ρi )
kk ≠=1i,
−1 n
∑ (rk − (n −1) ri ) , i = 1, n , (33)
k =1
where ri = ( Ki1 (1− m−1ρ ) − Ki2 (1− m−1ρi )) ( Ki1 (1− m−1ρ ) − Ki2 ) and the values of
−1
the coefficients Ki1 and Ki2 satisfy the inequalities Ki1Ki−21 &gt; (1− m−1ρi ) (1− m−1ρ ) .</p>
      <p>Thus the section synthesizes mathematical models for the identification of
Paretooptimal schemes for control of accessibility of ISCU with a parallel functioning set of
input information ports described by a composition of single-line Markov or
semiMarkov QS. The optimization tasks formulated differ in the view of the loss functions
of APs which allow one to calculate a quantitative or cost estimate of the identified
Pareto-optimal scheme for control of availability of ISCU which takes into account
the processing of emergency access requests indicated by the corresponding priority.
5</p>
    </sec>
    <sec id="sec-5">
      <title>Experiment and results</title>
      <p>Approbation of the models proposed in the article for Pareto-optimal control of
availability of ISCU was carried out on a real object - ISCU of the Situation Center of the
Information Technology Department of Vinnytsia City Council. The Situation Center
has been operational since 2018. Its main function is to capture and store video events
in the city, generalizing the incoming video stream from more than 550 camcorders.
Employees of the Vinnitsa City Council, the National Police of Ukraine, the Security
Service of Ukraine, and municipal enterprises have access to a distributed database of
video recordings. In the investigated ISCU B open information ports are allocated
for ordinary requests servicing and С − В open information ports are allocated for
priority requests, where С is the total number of open information ports available for
the APs. Ordinary and priority requests are assumed to arrive at a set of open ISCU
ports with intensities λО and λР respectively. New requests from the APs at the time
of receipt of which less than В open information ports are free in ISCU are
forwarded to queue Q1 with capacity NO . Priority requests at the time of receipt of which all
the С open information ports of ISCU are occupied are forwarded to the queue Q2
with capacity NР . Both queues are served by information environment of ISCU in
accordance with the FCFS algorithm. If a priority access request is in the queue Q2 it
will be satisfied with the ISCU as soon as at least one of the С open information
ports is released. If an ordinary access request is registered in queue Q1 then it will be
satisfied with the ISCU information environment if queue Q2 is empty and the
number of free open information ports exceeds В . If at some time moment both queues of
ISCU are filled then when a new request is received it is rejected. An arbitrarily
accepted by ISCU input request is satisfied with intensity µ + or rejected by a system
with intensity µ − . In general an arbitrary query is analyzed by ISCU over a time
interval whose value is a stochastic variable with exponential distribution and
mathematical expectation µ −1 =(µ + + µ − )−1 . The AP whose request is pending in the queue
Q1 can cancel it with the intensity µ − = µО− . Similarly the AP whose request is
pending in queue Q2 can cancel it with intensity µ Р− . The processes that determine the
availability of the investigated ISCU are described in the structural scheme shown in
Fig. 1.</p>
      <p>=(t) {( L1 , L2 (t )) , t ≥ 0} which is realized at a continuous time and can be in
discrete
states
defined
by
the
set</p>
      <p>X
={(i, j ) : B ≤ i ≤ C + NP , 0 ≤ j ≤ NO }
{(i, 0) : 0 ≤ i ≤ B −1} . In this case the functions L1 (t ) and L2 (t ) will describe the
total number of incoming requests received by the ISCU and at the time moment
t ≥ 0 waiting for satisfaction in the queues Q1 and Q2 respectively. We classify the
formalized random process X (t ) as a Markov process over the state space Х and
visualize by the transition intensity diagram shown in the form of the UML state
diagram in Fig. 2.</p>
      <p>Applying the mathematical apparatus described above to numerically solve the
optimization task of control the availability of the investigated ISCU described by QS
we find the Pareto-optimal equilibrium point х∗ = { р (i, j )} which identifies the
stationary probabilities p (i, j ) , (i, j ) ∈ X of states of a Markov process X (t ) .</p>
      <p>At the identified equilibrium point х∗ it is possible to analytically formalize such
applied probabilistic characteristics of the investigated ISCU as:
- the probability of an incoming request being rejected as a result of the queue Q1
overflow: RO =</p>
      <p>C + NP
∑ p (i, NO ) ;
i=B
the probability of an incoming request being rejected as a result of the queue Q2
overflow: RP
- the probability that a priority request will be received after a priority input request
to the ISCU: PPP = A− AP−1 ;
- the probability of premature termination of an active session sanctioned by the
investigated ISCU: Pf = RP + (1− RP ) PPP .</p>
      <p>In Fig. 3 shows the results of a numerical experiment to identify the dependence
of the probability of rejection of new incoming requests from the intensity of their
receipt calculated for the investigated ISCU with parameters B = 80 , C − B =20 ,
NO = 3⋅103 , NP = 2 ⋅103 , µ + = 0.5 , µ − = 0.5 , µ P− = 8 .</p>
      <p>Parameters the probability of rejection of new incoming requests and the intensity
of their receipt are the main ones that determine the availability of ISCU. Studies have
shown that with increasing queue Q1 capacity the probability of rejection of new
incoming requests is significantly reduced. By varying the value of the queue Q2
capacity with unlimited increase of the intensity of new input requests it is possible to
decompose the model of availability of the investigated ISCU into two independent
models of QS of the form M / M / B / QO and M / M / C − B / QP because the high
value of the probability of rejection of new input requests due to the increase of λО
the value of intensity λР stabilizes and ceases to depend on the value of λО .
However, increasing the capacity of queue Q2 contributes to the positive dynamics of the
probability of rejection of received priority requests but increases the probability of
rejection of new incoming requests.</p>
      <p>Further investigations of the effectiveness of the optimal schemes of control the
availability of ISCU was performed while varying the number of open information
ports in the range С = {30, 50,100} and the value of the reduced intensity of the input
requests which for queue Q1 would be defined as ρ1 =µ(λO + λP ) and for queue Q2
as ρ 2 = µλР . In further investigations we will take ρ = ρ1 at ρ 2 = const which will
allow us to represent the results of experiments in the two-dimensional space.</p>
      <p>Determine the empirical dependence of PW – the probability that the incoming
request received by the investigated ISCU will be allocated into the queue, from the
values of parameters ρ and С . Taking into account the model of the investigated
ISCU shown in Fig. 1 and 2 we describe analytically functional dependence
PW = f (ρ , C ) by the expression PW = ρ1С ((С −1)!(С − ρ 2 ))−1 Р0 where the
probability P0 = (1+ ρ1 (1!)−1 + ρ12 (2!)−1 + + ρ1C−1 ((C −1)!)−1 +
+ρ1C ((C −1)!(C − ρ 2 ))−1 )−1 .</p>
      <p>The empirical values of PW for the optimal scheme of control the availability of
investigated ISCU for С = {30, 50,100} and ρ = 0, 0.1,,10 are shown in Fig. 4.</p>
      <p>From the results shown in Fig. 4 it can be seen that when the value of ρ = 1 when
the intensity of the analysis of the incoming requests by ISCU and the intensity of the
receipt of new requests coincide there are enough 30 &lt; С &lt; 50 &lt; 100 active input
information ports to prevent the incoming requests from being rejected by the
investigated system. It is also obvious that the dependence of PW = f (ρ , C ) is nonlinear but
as the value of parameter С increases the steepness of the value of parameter PW
decreases.</p>
      <p>Determine the empirical dependence of LQ – the average length of the queue of
received incoming requests in the investigated ISCU on the parameters ρ and С .
Taking into account the model of the investigated ISCU shown in Fig. 1 and 2 we
describe analytically the functional dependence of LQ = f (ρ , C ) by the expression
LQ =(С − ρ 2 )−1 ρ P . The empirical values of LQ for the optimal scheme of control
2 W
the availability of investigated ISCU for С = {30, 50,100} and ρ = 0, 0.1,,10 are
shown in Fig. 5.</p>
      <p>From the results shown in Fig. 5 it can be seen that for optimal scheme of control
the availability of investigated ISCU the queue begins to form when the value of the
intensity of the new requests is exceeded the intensity of the analysis of the incoming
requests by the ISCU is more than: - 3 times when С = 30 ; - 5 times when С = 50 ;
9 times when С = 100 . Therefore the Pareto-optimal scheme of control the
availability of investigated ISCU is effective both with a limited number of input information
ports and demonstrates a nonlinear increase in the efficiency of use of the input of the
information system as it expands.</p>
      <p>Determine the empirical dependence of twq – the mean time of the received input
request staying in the queue waiting for access on the parameters ρ and С . Taking
into account the model of the investigated ISCU shown in Fig. 1 and 2 we describe
analytically the functional dependence of twq = f (ρ , C ) by the expression
twq
=(µ (C − ρ 2 ))−1 P . The empirical values of twq for the optimal scheme of control</p>
      <p>W
the availability of investigated ISCU for С = {30, 50,100} and ρ = 0, 0.1,,10 are
shown in Fig. 6.</p>
      <p>From the results shown in Fig. 6 it can be seen that for optimal scheme of control
the availability of investigated ISCU the mean time of the received input request
staying in the queue waiting for access is non-zero at: - ρ &gt; 0.3 when С = 30 ; - ρ &gt; 1.9
when С = 50 ; - ρ &gt; 5 when С = 100 . The steep rise in the value of the parameter
twq with the increase of the value of the parameter ρ decreases markedly with the
increase of the value of the parameter С which confirms the effectiveness of the
Pareto-optimal scheme of control the availability of investigated ISCU.
The relevance of the task of control the availability of information systems in general
and of ISCU in particular is increasing simultaneously with the information needs of
humanity. However given the importance of an information resources of ISCU
introducing new schemes of control the availability without a simulation stage can lead to
significant material and in particular human losses. These circumstances led to the
expediency of research the scientific and applied results of which are presented in the
article.</p>
      <p>In the article the new mathematical models for identification of the Pareto-optimal
schemes of control of accessibility of ISCU with a parallel functioning set of input
information ports are synthesized. The ISCU is described by the composition of
single-line Markov or semi-Mark queuing systems in which unlike existing ones
availability loss functions are formulated with a focus on the ability of quantitative or
qualitative estimation of the scheme of control of availability and guarantee the processing
of access requests with different priority degree which allows to set and solve the
multicriteria optimization task for an identification of problem-oriented scheme of
control of availability of ISCU taking into account architectural features of the latter.
In the article on the basis of the identified Pareto-optimal scheme of control of
availability of ISCU which represented by the equilibrium point analytically formalized
such applied probabilistic characteristics of an investigated system as: the probability
of an incoming request being rejected as a result of the queue of ordinary requests
overflow; the probability of an incoming request being rejected as a result of the
queue of priority requests overflow; the probability that an incoming request received
by the system will be queued; the intensity of the flow of new incoming requests
rejected by the ISCU; the intensity of the flow of priority requests received by the
ISCU; the probability that a priority request will be received after a priority input
request to the ISCU; the probability of premature termination of an active session
sanctioned by the investigated ISCU.</p>
      <p>To verify the models synthesized in the article the Pareto-optimal scheme of
control of availability of information system of the Situation Center of the Information
Technology Department of the Vinnytsia City Council was implemented. It turned out
that parameters the probability of rejection of new incoming requests and the intensity
of their receipt are the main ones that determine the availability of ISCU.</p>
      <p>Further research is planned to be devoted to the analysis of the conformity of the
obtained Pareto-optimal scheme of control of availability of ISCU with current
information security standards.
7</p>
    </sec>
    <sec id="sec-6">
      <title>Acknowledgements</title>
      <p>The research was carried out within the framework of the departmental scientific
research work number 46K4 "Methods of modeling and optimization of complex
systems on the basis of intellectual technologies" at the Department of Computer
Control Systems (CCS) of Vinnytsia National Technical University (VNTU) with the
staff Department of Automation and Intelligent Information Technologies.</p>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          <source>lenges of Critical and Extreme Dependability and Security</source>
          .
          <source>2017 IEEE 22nd Pacific Rim</source>
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          <source>International Symposium on Dependable Computing (PRDC)</source>
          ,
          <year>Christchurch</year>
          .
          <fpage>92</fpage>
          -
          <lpage>97</lpage>
          (
          <year>2017</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          <source>doi: 10</source>
          .1109/PRDC.
          <year>2017</year>
          .
          <volume>21</volume>
          . 1.
          <string-name>
            <surname>Shang</surname>
            ,
            <given-names>Z.:</given-names>
          </string-name>
          <article-title>The Performance and Dependability of Control Plane in OpenFlow Networks</article-title>
          .
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          <source>2019 IEEE International Symposium on Software Reliability Engineering Workshops</source>
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          <source>(ISSREW)</source>
          , Berlin, Germany.
          <fpage>108</fpage>
          -
          <lpage>111</lpage>
          (
          <year>2019</year>
          ). doi:
          <volume>10</volume>
          .1109/ISSREW.
          <year>2019</year>
          .
          <volume>00054</volume>
          . 2.
          <string-name>
            <surname>Zhao</surname>
            ,
            <given-names>F.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Jin</surname>
            ,
            <given-names>H.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Zou</surname>
            ,
            <given-names>D.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Qin</surname>
            ,
            <given-names>P.</given-names>
          </string-name>
          :
          <article-title>Dependability analysis for fault-tolerant computer sys-</article-title>
        </mixed-citation>
      </ref>
      <ref id="ref6">
        <mixed-citation>
          <article-title>tems using dynamic fault graphs,"</article-title>
          <source>in China Communications</source>
          .
          <volume>9</volume>
          (
          <issue>11</issue>
          ),
          <fpage>16</fpage>
          -
          <lpage>30</lpage>
          (
          <year>2014</year>
          ). doi:
        </mixed-citation>
      </ref>
      <ref id="ref7">
        <mixed-citation>
          10.1109/CC.
          <year>2014</year>
          .
          <volume>6969708</volume>
          . 3.
          <string-name>
            <surname>Abidi</surname>
            ,
            <given-names>S. S. A.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Farooqui</surname>
            ,
            <given-names>M. F.</given-names>
          </string-name>
          :
          <article-title>Software Dependability: Classification and Categoriza-</article-title>
        </mixed-citation>
      </ref>
      <ref id="ref8">
        <mixed-citation>
          tion.
          <source>2019 6th International Conference on Computing for Sustainable Global Develop-</source>
        </mixed-citation>
      </ref>
      <ref id="ref9">
        <mixed-citation>
          <string-name>
            <surname>ment (INDIACom)</surname>
          </string-name>
          , New Delhi, India.
          <fpage>1124</fpage>
          -
          <lpage>1128</lpage>
          (
          <year>2019</year>
          ). 4.
          <string-name>
            <surname>Abidi</surname>
            ,
            <given-names>S. S. A.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Farooqui</surname>
            ,
            <given-names>M. F.</given-names>
          </string-name>
          :
          <article-title>Analytical Hierarchical Process for Software Dependabil-</article-title>
        </mixed-citation>
      </ref>
      <ref id="ref10">
        <mixed-citation>
          ity.
          <source>2019 6th International Conference on Computing for Sustainable Global Development</source>
        </mixed-citation>
      </ref>
      <ref id="ref11">
        <mixed-citation>
          <source>(INDIACom)</source>
          , New Delhi, India.
          <fpage>1184</fpage>
          -
          <lpage>1188</lpage>
          (
          <year>2019</year>
          ). 5.
          <string-name>
            <surname>Bisikalo</surname>
            ,
            <given-names>O. V.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Grischuk</surname>
            ,
            <given-names>T. V.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Kovtun</surname>
            ,
            <given-names>V. V.</given-names>
          </string-name>
          :
          <article-title>The automatic speaker recognition system</article-title>
        </mixed-citation>
      </ref>
      <ref id="ref12">
        <mixed-citation>
          <article-title>of critical use classifier optimization</article-title>
          .
          <source>Radio Electronics</source>
          , Computer Science, Control.
          <volume>4</volume>
          ,
          <fpage>71</fpage>
          -
        </mixed-citation>
      </ref>
      <ref id="ref13">
        <mixed-citation>
          <volume>84</volume>
          (
          <year>2018</year>
          ). doi:
          <volume>10</volume>
          .15588/
          <fpage>1607</fpage>
          -3274-2018-4
          <issue>-7</issue>
          . 6.
          <string-name>
            <surname>Bisikalo</surname>
            ,
            <given-names>O. V.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Kovtun</surname>
            ,
            <given-names>V. V.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Yukhimchuk</surname>
            ,
            <given-names>M. S.</given-names>
          </string-name>
          :
          <article-title>Modeling the security policy of the</article-title>
        </mixed-citation>
      </ref>
      <ref id="ref14">
        <mixed-citation>
          <article-title>information system for critical use</article-title>
          .
          <source>Radio Electronics</source>
          , Computer Science, Control,
          <volume>1</volume>
          ,
          <fpage>132</fpage>
          -
        </mixed-citation>
      </ref>
      <ref id="ref15">
        <mixed-citation>
          <volume>149</volume>
          (
          <year>2019</year>
          ). doi:
          <volume>10</volume>
          .15588/
          <fpage>1607</fpage>
          -3274-2019-1-
          <issue>13</issue>
          . 7.
          <string-name>
            <surname>Bykov</surname>
            ,
            <given-names>M. M.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Kovtun</surname>
            ,
            <given-names>V. V.</given-names>
          </string-name>
          , et al:
          <article-title>Research of neural network classifier in speaker</article-title>
        </mixed-citation>
      </ref>
      <ref id="ref16">
        <mixed-citation>
          <article-title>recognition module for automated system of critical use</article-title>
          .
          <source>Proc. SPIE</source>
          .
          <volume>10445</volume>
          ,
          <fpage>1044521</fpage>
        </mixed-citation>
      </ref>
      <ref id="ref17">
        <mixed-citation>
          (
          <year>2017</year>
          ).
          <source>doi:10.1117/12.2280930</source>
          . 8.
          <string-name>
            <surname>Bykov</surname>
            ,
            <given-names>M. M.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Kovtun</surname>
            <given-names>V. V.</given-names>
          </string-name>
          , et al:
          <article-title>The automated speaker recognition system of critical</article-title>
        </mixed-citation>
      </ref>
      <ref id="ref18">
        <mixed-citation>
          use.
          <source>Proc. SPIE</source>
          .
          <volume>10808</volume>
          ,
          <issue>108082V</issue>
          (
          <year>2018</year>
          ).
          <source>doi: 10.1117/12.2501688</source>
          . 9.
          <string-name>
            <surname>Bisikalo</surname>
            ,
            <given-names>O. V.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Kovtun</surname>
            <given-names>V. V.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Sholota</surname>
            ,
            <given-names>V. V.</given-names>
          </string-name>
          :
          <article-title>The Information System for Critical Use</article-title>
        </mixed-citation>
      </ref>
      <ref id="ref19">
        <mixed-citation>
          <string-name>
            <given-names>Access</given-names>
            <surname>Process Dependability Modeling</surname>
          </string-name>
          . 9th International Conference on Advanced Com-
        </mixed-citation>
      </ref>
      <ref id="ref20">
        <mixed-citation>
          <string-name>
            <surname>puter Information</surname>
          </string-name>
          <article-title>Technologies (ACIT), Ceske Budejovice</article-title>
          ,
          <source>Czech Republic. 5-8</source>
          (
          <year>2019</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref21">
        <mixed-citation>
          <source>DOI: 10.1109/ACITT</source>
          .
          <year>2019</year>
          .
          <volume>8780013</volume>
          . 10.
          <string-name>
            <surname>Zhou</surname>
            ,
            <given-names>X.</given-names>
          </string-name>
          :
          <source>Research on Optimization of Queuing System Based on Computer Simulation.</source>
        </mixed-citation>
      </ref>
      <ref id="ref22">
        <mixed-citation>
          <source>(ICCTEC)</source>
          , Dalian, China.
          <fpage>9</fpage>
          -
          <lpage>12</lpage>
          (
          <year>2017</year>
          ). doi:
          <volume>10</volume>
          .1109/ICCTEC.
          <year>2017</year>
          .
          <volume>00011</volume>
          . 11.
          <string-name>
            <surname>Balla</surname>
            ,
            <given-names>H. A. M. N.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Sheng</surname>
            ,
            <given-names>C. G.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Weipeng</surname>
          </string-name>
          , J.: Reliability Enhancement in Cloud Compu-
        </mixed-citation>
      </ref>
      <ref id="ref23">
        <mixed-citation>
          <source>Queuing Theory. 1st International Conference on Data Intelligence and Security</source>
          (ICDIS),
        </mixed-citation>
      </ref>
      <ref id="ref24">
        <mixed-citation>
          <string-name>
            <surname>South Padre</surname>
            <given-names>Island</given-names>
          </string-name>
          , TX.
          <fpage>127</fpage>
          -
          <lpage>130</lpage>
          (
          <year>2018</year>
          ).
          <volume>12</volume>
          .
          <string-name>
            <surname>Li</surname>
            ,
            <given-names>M.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Fang</surname>
          </string-name>
          , Y.:
          <article-title>Research and Application of the Queuing Theory of the Virtualization</article-title>
        </mixed-citation>
      </ref>
      <ref id="ref25">
        <mixed-citation>
          <string-name>
            <surname>tion</surname>
          </string-name>
          ,
          <source>Big Data &amp; Smart City (ICITBS)</source>
          ,
          <year>Changsha</year>
          .
          <fpage>245</fpage>
          -
          <lpage>248</lpage>
          . (
          <year>2016</year>
          ). doi:
        </mixed-citation>
      </ref>
      <ref id="ref26">
        <mixed-citation>
          10.1109/ICITBS.
          <year>2016</year>
          .
          <volume>123</volume>
          . 13.
          <string-name>
            <surname>Yuping</surname>
            ,
            <given-names>Z.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Chen</surname>
            ,
            <given-names>L.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Bing</surname>
            ,
            <given-names>L.:</given-names>
          </string-name>
          <article-title>The real-time information scheduling method of industrial</article-title>
        </mixed-citation>
      </ref>
      <ref id="ref27">
        <mixed-citation>
          <article-title>control network based on queuing theory</article-title>
          .
          <source>Proc. of 2nd International Conference on Meas-</source>
        </mixed-citation>
      </ref>
      <ref id="ref28">
        <mixed-citation>
          <string-name>
            <surname>urement</surname>
          </string-name>
          , Information and Control, Harbin.
          <fpage>225</fpage>
          -
          <lpage>228</lpage>
          (
          <year>2013</year>
          ). doi:
        </mixed-citation>
      </ref>
      <ref id="ref29">
        <mixed-citation>
          10.1109/MIC.
          <year>2013</year>
          .
          <volume>6757953</volume>
          . 14.
          <string-name>
            <surname>Dudin</surname>
            ,
            <given-names>A.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Nazarov</surname>
            ,
            <given-names>A.</given-names>
          </string-name>
          :
          <source>Information Technologies and Mathematical Modelling</source>
          , Queuing
        </mixed-citation>
      </ref>
      <ref id="ref30">
        <mixed-citation>
          <source>Theory and Applications</source>
          . Springer.
          <volume>443</volume>
          (
          <year>2015</year>
          ). ISBN-
          <volume>10</volume>
          :
          <fpage>3319258605</fpage>
          . 15.
          <string-name>
            <surname>Rahmandad</surname>
            ,
            <given-names>H.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Oliva</surname>
            ,
            <given-names>R.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Osgood</surname>
            ,
            <given-names>N. D.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Richardson</surname>
          </string-name>
          , G.:
          <article-title>Optimal Control for Complex</article-title>
        </mixed-citation>
      </ref>
      <ref id="ref31">
        <mixed-citation>
          <string-name>
            <surname>Systems</surname>
          </string-name>
          .
          <article-title>Analytical Methods for Dynamic Modelers</article-title>
          . MIT Press.
          <volume>337</volume>
          -
          <fpage>372</fpage>
          (
          <year>2015</year>
          ).
          <fpage>16</fpage>
          .
          <string-name>
            <surname>Lv</surname>
            ,
            <given-names>X.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Fan</surname>
            ,
            <given-names>B.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Wang</surname>
            ,
            <given-names>X.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Zhao</surname>
            ,
            <given-names>X.:</given-names>
          </string-name>
          <article-title>A modeling and simulation method of queue systems</article-title>
        </mixed-citation>
      </ref>
      <ref id="ref32">
        <mixed-citation>
          <source>based on stateflow. 3rd International Conference on Systems and Informatics (ICSAI),</source>
        </mixed-citation>
      </ref>
      <ref id="ref33">
        <mixed-citation>
          Shanghai. 541-
          <fpage>545</fpage>
          (
          <year>2016</year>
          ). doi:
          <volume>10</volume>
          .1109/ICSAI.
          <year>2016</year>
          .
          <volume>7811014</volume>
          . 17.
          <string-name>
            <surname>Sericola</surname>
            ,
            <given-names>B.</given-names>
          </string-name>
          :
          <article-title>Markov Chains: Theory and Applications</article-title>
          . Wiley.
          <volume>416</volume>
          (
          <year>2014</year>
          ).
          <fpage>18</fpage>
          .
          <string-name>
            <surname>Ching</surname>
          </string-name>
          , W.-K., et al.:
          <source>Markov Chains: Models, Algorithms and Applications</source>
          . Springer. 258
        </mixed-citation>
      </ref>
      <ref id="ref34">
        <mixed-citation>
          (
          <year>2013</year>
          ). ISBN:
          <volume>1461463114</volume>
          ,
          <fpage>9781461463115</fpage>
          . 19.
          <string-name>
            <surname>Huang</surname>
            ,
            <given-names>J.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Lin</surname>
            <given-names>C.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Wan</surname>
          </string-name>
          , J.:
          <article-title>Modeling, Analysis and Optimization of Dependability-Aware</article-title>
        </mixed-citation>
      </ref>
      <ref id="ref35">
        <mixed-citation>
          <string-name>
            <surname>vices Computing</surname>
          </string-name>
          , Santa Clara, CA.
          <fpage>683</fpage>
          -
          <lpage>690</lpage>
          (
          <year>2013</year>
          ). doi:
          <volume>10</volume>
          .1109/SCC.
          <year>2013</year>
          .
          <volume>63</volume>
          . 20.
          <string-name>
            <surname>Kohlik</surname>
            ,
            <given-names>M.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Kubatova</surname>
          </string-name>
          , H.:
          <source>Hierarchical Dependability Models Based on Markov Chains.</source>
        </mixed-citation>
      </ref>
      <ref id="ref36">
        <mixed-citation>
          <source>26th International Conference on Architecture of Computing Systems</source>
          , Prague, Czech Re-
        </mixed-citation>
      </ref>
      <ref id="ref37">
        <mixed-citation>
          public. 1-
          <fpage>8</fpage>
          (
          <year>2013</year>
          ).
          <fpage>21</fpage>
          .
          <string-name>
            <surname>Herrmann</surname>
          </string-name>
          , J.:
          <source>Handbook of Operations Research for Homeland Security</source>
          . Springer. 232
        </mixed-citation>
      </ref>
      <ref id="ref38">
        <mixed-citation>
          (
          <year>2013</year>
          ). ISBN:
          <volume>9781461452782</volume>
          . 22.
          <string-name>
            <surname>Bhunia</surname>
            ,
            <given-names>A. K.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Laxminarayan</surname>
            ,
            <given-names>S.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Shaikh</surname>
            ,
            <given-names>A. A.</given-names>
          </string-name>
          :
          <article-title>Advanced Optimization</article-title>
          and Operations
        </mixed-citation>
      </ref>
      <ref id="ref39">
        <mixed-citation>
          Research. Springer.
          <volume>626</volume>
          (
          <year>2019</year>
          ).
          <source>ISBN 9789813299665</source>
          . 23.
          <string-name>
            <surname>Bauso</surname>
            ,
            <given-names>D.</given-names>
          </string-name>
          : Game Theory: Models,
          <string-name>
            <given-names>Numerical</given-names>
            <surname>Methods</surname>
          </string-name>
          and Applications. Game Theory:
        </mixed-citation>
      </ref>
      <ref id="ref40">
        <mixed-citation>
          <string-name>
            <surname>Models</surname>
            ,
            <given-names>Numerical</given-names>
          </string-name>
          <string-name>
            <surname>Methods</surname>
          </string-name>
          and Applications.
          <source>Now Foundations and Trends</source>
          .
          <volume>166</volume>
          (
          <year>2014</year>
          ).
          <fpage>24</fpage>
          .
          <string-name>
            <surname>Pavel</surname>
            ,
            <given-names>L.</given-names>
          </string-name>
          :
          <article-title>Game Theory for Control of Optical Networks</article-title>
          . Springer.
          <volume>260</volume>
          (
          <year>2012</year>
          ). 25.
          <string-name>
            <surname>McEachern</surname>
            ,
            <given-names>A.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Ashlock</surname>
            ,
            <given-names>D.</given-names>
          </string-name>
          :
          <article-title>Game Theory: A Classical Introduction</article-title>
          , Mathematical
        </mixed-citation>
      </ref>
      <ref id="ref41">
        <mixed-citation>
          <string-name>
            <surname>Games</surname>
          </string-name>
          , and the Tournament. Morgan &amp; Claypool.
          <volume>117</volume>
          (
          <year>2017</year>
          ).
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>