Next Article in Journal
Improving the Return Loading Rate Problem in Northwest China Based on the Theory of Constraints
Next Article in Special Issue
A New Approach of Soft Joint Based on a Cable-Driven Parallel Mechanism for Robotic Applications
Previous Article in Journal
Bitcoin and Fiat Currency Interactions: Surprising Results from Asian Giants
Previous Article in Special Issue
Modular and Self-Scalable Origami Robot: A First Approach
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Study of Synergistic Effects in Complex Stochastic Systems

by
Gurami Tsitsiashvili
Institute for Applied Mathematics, Far Eastern Branch of Russian Academy Sciences, Radio Str. 7, IAM FEB RAS, 690041 Vladivostok, Russia
Mathematics 2021, 9(12), 1396; https://doi.org/10.3390/math9121396
Submission received: 23 March 2021 / Revised: 19 May 2021 / Accepted: 13 June 2021 / Published: 16 June 2021
(This article belongs to the Special Issue Control, Optimization, and Mathematical Modeling of Complex Systems)

Abstract

:
In this paper, a method for detecting synergistic effects of the interaction of elements in multi-element stochastic systems of separate redundancy, multi-server queuing, and statistical estimates of nonlinear recurrent relations parameters has been developed. The detected effects are quite strong and manifest themselves even with rough estimates. This allows studying them with mathematical methods of relatively low complexity and thereby expand the set of possible applications. These methods are based on special techniques of the structural analysis of multi-element stochastic models in combination with majorant asymptotic estimates of their performance indicators. They allow moving to more accurate and rather economical numerical calculations, as they indicate in which direction it is most convenient to perform these calculations.

1. Introduction

Questions of composition and decomposition in multi-element stochastic systems are relevant for solving a number of problems. These include paralleling of algorithms and programs, modeling of supercomputers, the Internet, computer networks, mobile telephone communication systems, development of software packages for modeling catastrophic events in complex systems, design and improvement of technological and economic processes, and so on.
The term “synergy” (the result of the interaction of many elements of the system) originated in statistical physics, but recently it has been used by specialists from other fields: economics, biology, engineering, etc. Furthermore, research in these areas no longer leads to microscopic, but to phenomenological considerations. Here are some examples, taken from science history and devoting the detection of synergistic effects in complex systems, which have been obtained by famous researchers in their objective areas, using observation and mathematical intuition.
The economist A. Smith investigated the transition from shop production to manufacturing on the example of the production of safety pins. In the workshop method, the pin was made entirely by one master, performing all the operations sequentially. In manufacturing, each operation was performed by a separate master, which significantly increased labor productivity.
The physiologist I.P. Pavlov discovered the conditioned reflex by detecting feedback in the nervous system of the body. Stochastic feedback theory was developed by N. Wiener. A detailed study of the conditioned reflex led to the creation by P.K. Anokhin of the concept of a functional system that is urgently formed in the body when it is necessary to achieve the desired result and quickly disintegrates after it is achieved. N. Wiener and P.K. Anokhin collaborated in the development of this scientific direction, actively discussing the possibilities of mathematical methods in this area.
The physicist E. Rutherford discovered the atomic nucleus and proposed to P.L. Kapitsa to create an installation for the effect of a strong magnetic field on the atomic nucleus. Long-running installations were melted under the influence of a strong magnetic field. P.L. Kapitsa constructed an installation that creates a strong magnetic field for a short time, which turned out to be long enough for the processes occurring in the atomic nucleus.
Synergistic effects are a source of explicit dependencies between the characteristics of the system against the background of sufficiently large random perturbations. To study them, it was necessary to develop special techniques based on the structural analysis of multi-element stochastic models in combination with majorant asymptotic estimates of their performance indicators. This, in turn, required new techniques for working with statistical data, as well as skills in using the limit theorems of probability theory and the accompanying asymptotic expansions and estimates.
According to the author, many works on the analysis of complex multi-element systems require at the initial stage the construction of simpler models that allow you to determine the main performance indicators and the main parameters with which you can influence these indicators. For this purpose, it is very convenient to build procedures for comparing systems with different (alternative) structures to study their effectiveness with a large number of elements, with a large load, etc. For this purpose, schemes and/or modes of complex systems, computational algorithms, etc. can be used as objects of comparison. At the same time, at the initial stage of the study, a reasonable proportion should be observed between the accuracy of the calculations, which may be relatively small; the complexity of the calculations, which also should not be large; and the significance of the obtained results, which should be sufficiently large. Comparing systems with an alternative structures allows us to take these requirements into account, as with an increase in the number of elements, the differences between systems with an alternative structures are quite large.
The first author’s works, devoted to the study of synergistic effects, are analytical generalizations of the results of numerical and field experiments conducted by his colleagues in the modeling of telecommunications systems, container terminals, etc. In this connection, it should be remembered that in hydrodynamics, nonlinear soliton waves were also first discovered in the course of numerical experiments, and then their analytical theory was constructed. The use of computational experiments allows to obtain more accurate estimates of the synergistic effects. This can be used when working with models used in the programs “digital economy”, “smart city”, when modeling remote modes of operation that have become popular, on-line conferences, when using smart phones, etc. Currently, new information technologies are rapidly entering our life and their research helps us to adapt to them and to adapt these technologies themselves to the needs of potential users (for example, the use of smart phones by aged users). Nevertheless, analytical research helps to determine the direction of such research and to carry them out. In a sense, this avoids very complex structural optimization problems, a significant part of which are NP problems. Along with this, it becomes possible to use observations of complex systems, which also contributes to the study of synergistic effects in them.
In this paper, the synergistic effect is understood as a significant change in the performance indicators of a complex system when its structure changes, i.e., the connections between its elements. The complexity issues play an important role in modern systems analysis [1,2]. To reduce complexity, various techniques are used, among which the structural transformation of the system plays an important role [3,4]. This methodological technique is closely linked to the issues of the stability of a complex system [5].
Such a statement of the problem can be the comparing the reliability of separate and block reserving elements of a two-pole with unreliable edges [6]. This result is a classic in the mathematical theory of reliability and its refinement or amplification can be significant in itself. Note that the study of the reliability of two poles is widely used in various theoretical and applied studies (see, for example, in [7,8,9,10], etc. However, such a comparative analysis made in the monograph [6] of Barlow and Proshan did not develop in subsequent works, while the synergistic effects identified in this paper were very significant.
The peculiarity of this task is the use of probabilistic models, which, at first glance, complicate the task. This transition requires the selection of a new performance indicator—the required amount of reserve, which reflects the content of the reservation procedure and the novelty of the proposed approach. On the other hand, it is also necessary to construct sufficiently weak (logarithmic) dependencies of the so-introduced indicator on a number of the scheme elements. When selecting a new reserve efficiency indicator, its analogy with probabilistic metrics [11] is used. Therefore, the results obtained in this section are new, original, and significant.
Another problem considered in this paper and connected with synergistic effects appeared at the ITMM 2018 and ITMM 2020 conferences, when discussing the multi-server model of the RQ queuing system [12,13,14]. The Conference ITMM 2018 was held in conjunction with 12th International Workshop on Retrial Queues (RQ) and Related Topics with a wide representation of queuing theory researchers from Russia, India, Bulgaria, the Netherlands, and other countries. A good source of recent works on this topic is a collection of articles on RQ systems [15].
The RQ queuing system is a system in which a customer received in the presence of busy servers is not rejected, but is sent to the so-called orbit, from where it is extracted in accordance with some protocol for queuing, when one of the servers is released. A.N. Dudin remarked that it is most often assumed Poisson input flow to such a system. However, this requirement is not met because there is a dependence and even a long range dependence between the random variables, characterizing numbers of arrived customers in disjoint time intervals. Moreover, despite the large number of analytical results in which the distributions in RQ systems are calculated in formula form, their use for numerical calculations is difficult due to the high complexity of such calculations, especially for multi-server RQ systems.
The novelty of the proposed approach is that instead of a stationary distribution of the process describing a multi-server RQ system, the probability of customers appearing in the orbit of this system for a fixed period of time is investigated and the convergence of this indicator to zero is established when the number of channels proportional to the intensity of the input flow tends to infinity. Secondly, with an increase in the number of servers, even in a system with a Poisson input flow and exponentially distributed service times, the computational complexity of the problem of calculating the limit distribution in an RQ system increases quite strongly.
Thus, a new problem arises for calculating the RQ queuing system with a large number of servers and a non-Poisson input flow. Using asymptotic theorems for multichannel queuing systems, it is possible, on the contrary, to simplify the problem of analyzing RQ systems. For this purpose, it is convenient to use limit theorems based on topological concepts of convergence in the space of random processes defined on a finite time interval [16]. The significance of this approach lies in the broad scope of its application and in the ability to circumvent the computational complexity of the problem by reducing it to the limit theorems of probability theory.
A continuation of the study of multi-server queuing systems in this paper is the analysis of a system with failures. This system arises when modeling modern data transmission networks (of fifth generation), formulated by leading Russian specialists in the mathematical theory of communication [17]. This task is quite important and leading Russian specialists in modeling of transmission networks Samouylov K.E. and Gaidamaka Yu.V. even organized a seminar with the author’s participation to find alternative approaches to this problem solution with publishing of obtained results [18]. The solution to this problem is based on the recently installed a synergistic effect in multi-server queuing system with failures, when the stationary probability of failure tends to zero as the number of servers and proportionally the input flow intensity tend to infinity. Moreover, the obtained asymptotic results were quite accurate.
This study is based on the classical Erlangian model of loss multi-server system (see, for example, in [19,20]). Asymptotic behaviors of the blocking probability and parameters of the Equivalent Random Theory method was analyzed in [21] for the case when both the number of servers and the input flow intensity tend to infinity.
However, the inclusion in this model of the assumption that load factor equals one, the intensity of the input flow is proportional to the number of servers n, and the tendency of n to infinity allowed us to establish that the probability of failure tends to zero also. The exact asymptotic rate of this convergence is established. Moreover, when the load factor is less than one, it is possible to construct an upper estimate of the rate of convergence to zero in the form of a geometric progression. Therefore, the synergistic effect found in this paper is very strong and so can be used in the design of data transmission systems of the fifth generation.
The features of probabilistic models of complex systems discovered in this way can also be used in the estimation of their parameters. In particular, in the deterministic model of logistic growth [22] (which is very important and classical in mathematical biology), the problem of estimating the growth parameter from inaccurate observations arises and attracts specialists again and again. The solution of this problem by the method of least squares leads to quite large errors. In this paper, the unknown growth parameter is expressed in terms of the trajectory averages of the deterministic sequence of the model. In turn, the trajectory averages are estimated from observations over a sufficiently long period of time, which leads to the leveling of observation errors. These estimates are based on the use of probabilistic metrics developed in [11] and are new.
Thus, the solution of the above problems of system analysis required a combination of probabilistic and deterministic methods of system analysis, among which the methods of studying the synergistic effects arising from the structural restructuring of a complex system play a decisive role. The benefit of received results is to establish sufficiently strong dependencies of performance indicators on changes in the system structure. This approach opens up new opportunities in solving problems of structural optimization of stochastic systems: queuing, reliability, etc.

2. Separate Redundancy in a Two-Pole System

Consider m sequentially connected and independently operating elements with a failure-free probability of p , 0 < p < 1 . The probability of failure-free operation of such a chain is p m . Let us focus on two alternative ways to reserve this network. In the first method, n independently functioning duplicates are connected in parallel (see Figure 1).
Reliability of the network obtained in this way is H n ( m ) = 1 ( 1 p m ) n . In the second method, each element of the original chain is n-multiple reserved separately (see Figure 1). Reliability of the newly formed network H n ( m ) = ( 1 q n ) m , q = 1 p .
From the results of the monograph [6], it follows that H n ( m ) H n ( m ) (this inequality is valid for any bipolar). However, this inequality gives only a qualitative idea of the possibilities of separate reservation. To give a quantitative assessment, it is convenient to move from the reliability function to the required amount of reserve.
For δ > 0 , denote n * ( m , δ ) = min ( n : H n ( m ) 1 δ ) , n * ( m , δ ) = min ( n : H n ( m ) 1 δ ) the required volume of reserve, in which the reliability exceeds 1 δ . To calculate the reliability of general type two-pole, it is necessary to solve an NP-complex problem. However, to compare the different ways of reserving a chain of the length m , it is necessary only to solve a few simple inequalities. Moreover, the results of this comparison are very contrasting and the most interesting consideration of a separate reserve may be applied to general type two-pole also. Let us denote [ a ] the integer part of the real number a .
Proposition 1.
The following inequalities are met:
n * ( m , δ ) 1 δ p m + 1 , n * ( m , δ ) ln ( δ / m ) ln q + 1 .
Proof. 
Indeed, for all a , 0 < a < 1 , the inequality holds
( 1 a ) m 1 m a , m = 1 , 2 , , H n ( m ) = 1 ( 1 p m ) n n p m ,
n * ( m , δ ) min ( n : n p m 1 δ ) = min n : n 1 δ p m 1 δ p m + 1 ,
so the first relation in Formula (1) takes place. In turn,
H n ( m , δ ) ( 1 q n ) m 1 m q n , n * ( m , δ ) min ( n : m q n δ )
n * ( m , δ ) min ( n : ln m + n ln q ln δ ) ln ( δ / m ) ln q + 1 .
Therefore, the second relation in Formula (1) is valid.
Table 1 demonstrates how much n * ( m , δ ) is greater than n * ( m , δ ) .
A comparison of these relations shows that for a large chain length of m, the split-reservation scheme provides special advantages, as the lower bound, which grows as a geometric progression, is replaced by the upper logarithmic bound. Note that the upper estimate of the required reserve in the scheme of separate reservation of a sequential chain is logarithmic in m and can be easily extended to the general case.
Indeed, consider a two-pole consisting of m independently operating edges with probabilities of operation p 1 , , p m 1 q , 0 < q < 1 . Let us construct a two-pole in which each edge of the original two-pole is a reserve of n identical elements and denote H n ( p 1 , , p m ) the probability of the existence of a working path from the initial to the final vertex in this two-pole.
Proposition 2.
For the value n * ( p 1 , , p m , δ ) , the relation is valid
n * ( p 1 , , p m , δ ) ln ( δ / m ) ln q + 1 .
Proof. 
The proof of this statement is based on the inequality H n ( p 1 , , p n ) ( 1 ( 1 p 1 ) n ) · · ( 1 ( 1 p m ) n ) ( 1 q n ) m 1 m q n . This inequality follows from the fact that the reliability of an arbitrary two-pole with m independently functioning elements is not less than the reliability of a chain of m elements connected in series. Therefore, the second inequality in the formula (1) is true also. □
For an arbitrary two-pole, a logarithmic by m upper estimate of the value of the required reserve in the separate reservation scheme is performed. Note that this result is obtained using trivial inequalities and does not require calculating the reliability of H n ( p 1 , , p m ) , which in general is an NP-problem.
Indeed, if α 1 , , α m are independent boolean random variables which describe states of two-pole elements and boolean function A ( α 1 , , α m ) describes a workability of two-pole dependently on states of its elements then its reliability
H A ( p 1 , , p m ) = α 1 , , α m = 0 1 A ( α 1 , , α m ) j = 1 m p j α j , w i t h p j 1 = p j , p j 0 = 1 p j , j = 1 , , m .
Calculation of the reliability H ( p 1 , , p m ) formally requires performing of 2 m arithmetical operations.
Thus, a convenient choice of the reserve efficiency indicator in the form of the required reserve volume solves two problems. It allows us to obtain a strong (logarithmic) dependence of the chosen efficiency indicator on the number of edges of the two-pole m and makes it possible to abandon the solution of the NP-problem of calculating the reliability of the two-pole.

3. RQ-Queuing Systems with a Large Number of Servers

Consider an RQ-system, i.e., a queuing system, in which, if there is a free server, the customer that has come to the system immediately begins to be served on it. If there are no free servers, then the customer is sent to the orbit, from where it can return to the newly released server in accordance with some protocol [12,14,23]. A good source for RQ systems in recent years has been Conference ITMM 2018 in Tomsk, which was held in conjunction with 12th International Workshop on Retrial Queues and Related Topics (WRQ 2018).
To solve this problem, we propose to use the theorem on the asymptotic behaviour of an n-server queuing system for n . In this theorem, we prove that at T > 0 for n , the probability P n ( T ) that on the segment [ 0 , T ] in the system there will be customers going into orbit tends to zero. Already in this result, the transition from the limit distribution to the above probability is made. Moreover, this characteristic becomes a new indicator of efficiency, which is convenient to use when analyzing a multichannel RQ queuing system.
Consider the series scheme in which the characteristics of n-server queuing systems are defined by the parameter n , which characterizes an intensity of input flow tending to infinity. Denote e n ( t ) a number of input flow customers arriving before the moment t, e n ( 0 ) = 0 . Assume that q n ( t ) is a number of busy servers in this system at the moment t , q n ( 0 ) = 0 , τ j is the service time of j input flow customer and τ j , j 1 , is a sequence of independent and identically distributed random variables (s.i.i.d.r.v.’s) with the distribution function (d.f.) F ( t ) ( F ¯ = 1 F ) , which has continuous and bounded by f ¯ > 0 density f ( t ) . All results of this section are based on ([16], Chapter II, § 1, Theorem 1):
Theorem 1.
Assume that the following conditions are true.
(1)
For some a > 0 the equality E e n ( t ) = n a t , t 0 , takes place.
(2)
There is the function B ( n ) such that for A ( n ) = max ( n 1 / 2 , B ( n ) ) the limit relations take place for n
B ( n ) A ( n ) B 0 , n A ( n ) K 0 , n A ( n ) .
so that max ( B , K ) = 1 ).
(3)
The sequence of random processes x n ( t ) = e n ( t ) E e n ( t ) B ( n ) for n C-converges to the centred Gaussian process z ( t ) .
(4)
Random process ζ ( t ) = 0 t F ¯ ( t u ) d z ( u ) + K Θ ( t ) , 0 t T , where Θ ( t ) is centred Gaussian process independent with z ( t ) , which has the covariance function R ( t , t + u ) = 0 t F ¯ ( v + u ) F ( v ) a d v and satisfies the relation P ( sup 0 t T ζ ( t ) > L ) 0 , L .
(5)
If the inequality ρ = a E τ j < 1 is true then for any T > 0 we have the relation P sup 0 t T q n ( t ) n 0 , n .
Here, the concept of C-convergence used in Theorem 1 is defined as follows. Denote by F 1 the space of deterministic functions on the segment [ 0 , T ] with uniform metric ρ . Designate by F the set of bounded functionals f defined on F 1 and continuous in the metric ρ : if z = z ( t ) , z 1 = z 1 ( t ) , z 2 = z 2 ( t ) , F 1 and ρ ( z , z n ) 0 , n , then f ( z n ) f ( z ) , n . Say that the sequence of random processes z n = z n ( t ) , n 1 , C-converges to the random process z = z ( t ) if for any functional f F we have that E f ( z n ) E f ( z ) , n .
Deterministic input flow of customer groups. Let at times 1 , 2 , in n-server RQ-queuing system come groups of customers of the size of η 1 0 , η 2 0 , , where η 1 , η 2 , – i.i.d.r.v.‘s with integer values, E η 1 = a , V a r η 1 < . Define the input flow by the equality e n ( t ) = k = 1 [ n t + ψ ] η k , t 0 , where ψ – independent of η k , k 1 , τ j , j 1 , a random variable with a uniform distribution on the segment [ 0 , 1 ] ( [ d ] is the integer part of the real number d). Here and in two next models random variable ψ has uniform distribution to ensure the proportionality t of the mathematical expectation E e n ( t ) .
Theorem 2.
Suppose that, for some D > 0 , almost certainly η 1 < D and the inequality a E τ 1 < 1 is true. Then for any T > 0 the relation P n ( T ) 0 , n , is valid.
Proof. 
In [23] it is proved that under the conditions of this theorem,
P sup 0 t T q n ( t ) n 0 , n .
Connecting this relation with the inequality P n ( T ) P sup 0 t T q n ( t ) n , n 1 , one obtains the proof of the theorem. □
Alternating input flow. Consider a n-server RQ-queuing system, assuming n = n ( N ) , N . Let us define the input flow to this system using the following construction. Following the works in [24,25], we define a continuous random flow defined by ON and OFF periods. Let a sequence of i.i.d.r.v’s X 0 0 , X 1 0 , X 2 0 , consists of lengths of ON-periods, the sequence of i.i.d.r.v.‘s Y 0 0 , Y 1 0 , Y 2 0 , consists of the lengths of OFF-periods and these random sequences are independent. Denote F 1 ( t ) = P ( X 1 < t ) , F 2 ( t ) = P ( Y 1 < t ) , t 0 , and suppose that
F ¯ 1 ( t ) = t α 1 L 1 ( t ) , F ¯ 2 ( t ) = t α 2 L 2 ( t ) , 1 < α 1 < α 2 < 2 ,
where the function L 1 ( t ) l 1 > 0 , t , and L 2 ( t ) is a slowly varying function. Let b ( t ) is the inverse of the function 1 / F ¯ 1 ( t ) : b ( 1 / F ¯ 1 ( t ) ) = t .
We introduce independent r.v.‘s B , X , Y , which are independent of X n , Y n , n 1 , and Y 0 with distributions
P ( B = 1 ) = μ 1 μ , P ( B = 0 ) = μ 2 μ , μ = μ 1 + μ 2 , μ 1 = E X 1 , μ 2 = E Y 1 ,
P ( X x ) = 1 μ 1 0 x F ¯ 1 ( s ) d s , P ( Y x ) = 1 μ 2 0 x F ¯ 2 ( s ) d s .
Then, a random sequence
T 0 = B ( X + Y 0 ) + ( 1 B ) Y , T n = T 0 + i = 1 n ( X i + Y i ) , n 1 ,
generates an ON–OFF process
W ( t ) = B I [ 0 , X ) ( t ) + n = 0 I [ T n , T n + X n + 1 ) ( t ) , t 0
(here I A ( t ) is the indicator function of a random event t A ). The random process W ( t ) is binary: W ( t ) = 1 , if t is contained in an ON-period, W ( t ) = 0 , if t is contained in the OFF-period and stationary, and E W ( t ) = μ 1 / μ = α .
Denote A ( t ) = 0 t W ( s ) d s , then E A ( t ) = α t , t 0 . Let n = n ( N ) = N M ( N ) , M = M ( N ) = [ N γ ] , γ > 0 , and random functions A m ( t ) , m = 1 , . . . , M , are independent copies of a random function A ( t ) . We introduce the function e n ( t ) = m = 1 M A m ( N t ) + ψ , specifying the alternating input flow.
Theorem 3.
If γ > α 1 1 and α E τ j < 1 , then for any T > 0 the relation P n ( T ) 0 , n , is true.
The proof of Theorem 3 repeats the proof of Theorem 2 verbatim.
Erlangian input flow. Let E n ( t ) a Poisson flow of customers with intensity n α . Define the input flow to the n-server system described above by the equality
e n ( t ) = E n ( t ) r + ψ , t 0 ,
where ψ is a random variable independent of η k , k 1 , τ j , j 1 , with a uniform distribution on the segment [ 0 , 1 ] , and the fixed r takes natural values. It is obvious that for any fixed ψ , 0 ψ 1 , the moments of single jumps of the process e n ( t ) form an Erlangian flow. Here, the Erlangian flow is obtained from E n ( t ) by allocation of moments with numbers that are multiples of r .
Theorem 4.
If α E τ j < 1 , then for any T > 0 the relation P n ( T ) 0 , n , holds.
Proof. 
In [26] it is proved that Formula (2) is valid under the conditions of the theorem. Connecting it with inequality P n ( T ) P sup 0 t T q n ( t ) n , one obtains the proof of the theorem. □
The choice of the probability P n ( T ) as an efficiency indicator allows us to apply the known theorems to the analysis of a multi server RQ-system with a fairly general protocol for the transfer of customers from orbit to the vacated server almost without additional consideration.

4. Multiserver Loss Systems

Consider n-server queuing system M | M | n | 0 with a Poisson input flow of intensity n λ and exponentially distributed service times having intensity μ on all n servers, ρ = λ / μ . This system can be considered as combining n single-server systems with input flow intensities λ (see Figure 2).
The number of customers in the system M | M | n | 0 describes the process of death and birth with the intensities of birth and death λ n ( k ) = n λ , 0 k < n , μ n ( k ) = k μ , 0 < k n .
Let us denote P n ( ρ ) the stationary probability of failure in the system A n for a given ρ . It is not difficult to establish that P 1 ( 1 ) = 1 / 2 . However, the combined system A n satisfies new relation, which characterizes the synergistic effect of such a combination.
Theorem 5.
The following limit ratio is true: P n ( 1 ) 2 π n , n .
Proof. 
Let δ > 0 , consider the function f ( x ) = 1 x exp ( ( 1 + δ ) x ) . The f ( x ) function satisfying the condition: f ( 0 ) = 0 , f ( 1 ) < 0 , and such that the inequalities
f ( x ) > 0 , 0 < x < ln ( 1 + δ ) 1 + δ , f ( x ) < 0 , ln ( 1 + δ ) 1 + δ < x 1
hold. Therefore, on the half interval [ 0 , 1 ) there exists a single x ( δ ) , satisfying the condition f ( x ( δ ) ) = 0 and such that the inequalities 1 x exp ( ( 1 + δ ) x ) , 0 x x ( δ ) < 1 hold. Let p n ( k ) = lim t P ( x n ( t ) = k ) , 0 k n , then in force [16] [Chapter 2, § 1]
p n ( n 1 ) = p n ( n ) μ λ n n , p n ( n 2 ) = p n ( n ) μ λ 2 n ( n 1 ) n 2 ,
Therefore, the stationary blocking probability in virtue of the integral theorems of recovery and the law of large numbers for the recovery process [1] [Chapter 9, § 4, 5] satisfies the equality
P n ( ρ ) = p n ( n ) = k = 0 n ρ k j = 0 k 1 1 j n 1 ,
where j = 0 1 equals 1 . From Formula (3), we obtain the inequality
P n 1 ( 1 ) 0 k n x ( δ ) j = 0 k 1 1 j n 0 k n x ( δ ) j = 0 k 1 exp ( ( 1 + δ ) j / n )
1 k n x ( δ ) exp ( ( 1 + δ ) k 2 / 2 n ) .
This implies that
P n 1 ( 1 ) 1 n x ( δ ) e ( 1 + δ ) x 2 / 2 n d x = n 1 + δ 1 + δ n x ( δ ) n ( 1 + δ ) e y 2 / 2 d y ,
consequently
P n ( 1 ) n ( 1 + δ ) 1 + δ n x ( δ ) n ( 1 + δ ) e y 2 / 2 d y 1 ( 1 + δ ) 2 π , n .
and so lim sup n P n ( 1 ) π n 2 1 + δ .
Using Formula (3) and the inequality 1 x exp ( x ) , 0 x 1 , we obtain
P n 1 ( 1 ) 1 k n e k ( k 1 ) / 2 n 1 k n e ( k 1 ) 2 / 2 n 0 e x 2 / 2 n d x ,
thus it follows that 1 lim inf n P n ( 1 ) π n 2 . Obtained above inequalities for upper and lower limits lead to the statement of Theorem 5. □
Remark 1.
In aggregated M | M | n | 0 system at ρ < 1 following relations are valid [18]:
e n ln 2 ρ / 2 2 π n ρ 8 P n ( ρ ) ( e n ln 2 ρ / 2 ) ( ρ 1 ) / ln ρ 2 π n ln ρ ρ 1 .
And if ρ = ρ ( n ) = 1 n γ , γ > 0 , then Theorem 5 gives
1 2 1 π n P n ( ρ ) 2 π n , γ 1 2 ,
1 2 1 π n P n ( ρ ) exp n 1 2 γ 2 2 π n , γ < 1 2 .
Similar results were obtained for Erlang‘s loss function in [27,28] but in a more complex way.
Remark 2.
In aggregated M | M | n | system following relations are valid [29] for A n —stationary mean waiting time and B n —stationary mean queue length:
(1)
If ρ < 1 , then for some c < , q < 1 the relation holds A n c q n , n 1 .
(2)
If ρ = 1 n α , 0 < α < , then for n
A n 0 , α < 1 , 1 / μ , α = 1 , , α > 1 . B n 0 , α < 1 / 2 , , α 1 / 2 .
Suppose that we have m independently functioning n k -server queuing systems with Poisson input flows of intensity λ k , k = 1 , , m . In the k-th system, the customer of the input flow is served exponentially distributed time simultaneously on c k channels with intensity μ k . Let l k = n k / c k be a natural number and the equality ρ k = λ k / ( l k μ k ) = 1 holds.
We combine n copies of each of the n k -server systems under consideration, denoting P n k stationary probability of failure in each of the combined systems, k = 1 , , m . Using Theorem 5, it is not difficult to obtain the following limit relations
P n 1 2 π n l 1 , , P n m 2 π n l m , n .
This solution allows us to distribute the total number of n ( n 1 + + n m ) servers between flows so that the failure probabilities of customers of different flows tend to zero with the growth of a large parameter n . To solve this problem, one could use the exact multiplicative formula obtained in [17], but this would lead to significantly more complex calculations.

5. Parameter Estimation in the Logistics Growth Model

The recurrent model of logistic growth
x 0 = a , x n + 1 = b x n ( 1 x n ) , n = 0 , 1 , . . . ,
where the parameters a , b satisfy the conditions 0 < a < 1 , 1 < b < 4 , attracts increased attention from biologists and physicists. For this model, both practically and theoretically, it is important to evaluate the parameter b based on inaccurate observations. Due to the nonlinearity of the recurrence relation (5), the least squares method applied to the estimation of the parameter b seems somewhat unnatural, which is confirmed by numerous computational experiments that give quite large errors. It seems more natural to apply such qualitative properties of the sequence, as the existence of its limit cycle or limit distribution [30] depending on the value of b in combination with the method of probability metrics [11].
Consider an additive model for introducing errors in observations y n = x n + ε n , n = 1 , Here, ε n , n = 1 , , is a sequence of i.i.d.r.v.’s having a distribution with mean zero and variance σ 2 . We introduce the following notation
X n = i = 0 n 1 x i n , Y n = i = 0 n 1 y i n , X n = i = 0 n 1 x i 2 n , Y n = i = 0 n 1 y i 2 n .
Using the results of [30], it is possible to establish that for the deterministic sequence x n , n = 1 , , with a given b there are limits
lim n X n = x ¯ lim n X n = x 2 ¯ .
Indeed, say that the sequence x n , n = 1 , . . . , has a limit cycle x ( 1 ) , . . . , x ( q ) of length q 1 , if lim k x q k + j = x ( j ) , j = 1 , . . . , q . Denote x ¯ = 1 q j = 1 q x ( j ) , x 2 ¯ = 1 q j = 1 q [ x ( j ) ] 2 , then we have
X N q = 1 N q i = 1 N q x i x ¯ , X N q = 1 N q i = 1 N q x i 2 x 2 ¯ , N ,
so Formula (6) is true in the case, when the sequence x n , n = 1 , , has limit cycle.
Let p ( d x ) be a probability measure on the σ -algebra of Lebesgue-measurable subsets of the segment [ 0 , 1 ] . Let us say that p ( d x ) is the limiting distribution of the sequence x n , n = 1 , . . . , if for any Lebesgue-measurable set C [ 0 , 1 ] the equality holds lim n k ( C , n ) n = C p ( d x ) = p ( C ) , where k ( C , n ) is the number of x i satisfying the inclusion x i C , i = 1 , . . . , n . Then, we define x ¯ = 0 1 x p ( d x ) , x 2 ¯ = 0 1 x 2 p ( d x ) and prove Formula (6) as follows.
Let us take an arbitrary δ > 0 and put m = 2 δ + 1 , γ = δ 2 + 2 δ 1 . Divide the half interval [ 0 , 1 ) into disjoint segments
C 1 = 0 , 1 m , C 2 = 1 m , 2 m , , C m = m 1 m , 1 .
Choose N ( δ ) so that for any n > N ( δ ) we have k ( C j , n ) n γ , j = 1 , , m . It is sufficiently simple to prove for n N ( δ ) the following inequalities
1 n · x i C j , i = 1 , , n x i j m · k ( C j , n ) n j m ( p ( C j ) + γ )
γ j m + C j x + 1 m p ( d x ) = γ j m + p ( C j ) m + C j x p ( d x ) .
Summing these inequalities by j = 1 , . . . , m , and using the equality for m we get for n N ( δ ) the inequality
X n γ ( m + 1 ) 2 + 1 m + x ¯ = x ¯ + δ .
Analogously it is possible to obtain
1 n · x i C j , i = 1 , , n x i j 1 m · k ( C j , n ) n j 1 m ( p ( C j ) γ )
C j x 1 m p ( d x ) γ ( j 1 ) m = C j x p ( d x ) 1 m γ ( m 1 ) 2 x ¯ δ ,
consequently X n x ¯ , n . Similarly we have the relation X n x 2 ¯ , n , so Formula (6) is true in the case, when the sequence x n , n = 1 , , has limit distribution also.
Note that formally the limits x ¯ , x 2 ¯ may depend on the initial state x 0 . However, in the logistics growth model there is no such dependence.
We will evaluate the parameter b in two stages. First, we express b in terms of the path averages: b = x ¯ / ( x ¯ x 2 ¯ ) . Using the ratio
E Y n = 1 n i = 0 n 1 E ( x i + ε i ) = X n x ¯ ,
E Y n = 1 n i = 0 n 1 E ( x i + ε i ) 2 = 1 n i = 0 n 1 ( x i 2 + σ 2 ) = X n + σ 2 x 2 ¯ + σ 2 , n ,
let us estimate the parameter b by the formula
b n = E Y n E Y n ( E Y n σ 2 ) b , n .
As a result, the parameter b is evaluated by the formula b ^ n = Y n Y n ( Y n σ 2 ) . The convergence in probability b ^ n b , n , follows from the relations
V a r Y n = 1 n 2 i = 0 n 1 V a r ( x i + ε i ) = 1 n 2 i = 0 n 1 V a r ε i = σ 2 n 0 ;
V a r Y n = 1 n 2 i = 0 n 1 V a r ( x i + ε i ) 2 = 1 n 2 i = 0 n 1 V a r ( 2 x i ϵ i + ε i 2 ) 4 n ( 4 σ 2 + σ 4 ) 0 , n .
The following is an illustrative example of estimating parameter b for a logistic growth model. Calculations of b ^ n were performed for the case x 0 = 0.75 ; a = 0.5 ; b = 3 at n = 1000 (see Figure 3). An additive model of introducing errors was considered under the assumption that ε i , i = 0 , . . , n 1 , have a uniform distribution on the segment [ 1 / 4 , 1 / 4 ] .
This method can be applied to the estimation of the parameters of the Rikker model (see, for example, in [22]). Here, the Rikker model is described by recurrent relation
x 0 = 1 , x n + 1 = a x n exp ( b x n ) , a , b > 0 ,
and observations are following: y n = x n exp ( ε n ) , where ε n has normal distribution with zero mean and known variation, n 0 . Another application of described method is the finite-difference approximation of the system of Lorentz differential equations (see, for example, in [31]), etc.

6. Discussion

All the problems of system analysis considered in this paper are based on the choice of changes in the structure of the system, the efficiency indicator, and the computational algorithm with an assessment of its complexity. In some cases, it is possible to replace the NP-problem with a fairly simple computational procedure, abandoning the high accuracy of the resulting solution in favor of a significant change in the performance indicator. Apparently, such problems require a certain proportion between the accuracy and efficiency of the resulting solution.
The proposed approach to the study of synergistic effects in complex systems can be applied to the construction of queuing systems with a large load and a small queue, to backup systems with recovery, to insurance models and other stochastic systems. It allows you to explore and find the main parameters in such popular technologies in applications as powder metallurgy, 3-D printing, fast mixing of fuel in engines, etc. The emphasis on economical, but not highly accurate calculations, makes it possible at the initial stage to correctly select the main parameters of the analyzed systems before performing more detailed and accurate calculations. This property of the proposed approach to the analysis of complex systems can be used in programs of digital economy, smart city, etc., when at the initial stage of the study it is important to determine the main indicators of the effectiveness of a complex system.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing is not applicable to this article.

Acknowledgments

The author thanks Marina Osipova for her help in the design of the work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Skrimizea, E.; Haniotou, H.; Parra, C. On the complexity turn in planning: An adaptive rationale to navigate spaces and times of uncertainty. Plan. Theory 2019, 18, 122–142. [Google Scholar] [CrossRef]
  2. Battiston, S.; Caldarelli, G.; May, R.; Roukny, T.; Stiglitz, J. The price of complexity in financial networks. Proc. Natl. Acad. Sci. USA 2016, 113, 10031–10036. [Google Scholar] [CrossRef] [PubMed]
  3. Majdandzic, A.; Braunstein, L.; Curme, C.; Vodenska, I.; Levy-Carciente, S.; Eugene, S.; Havlin, S. Multiple tipping points and optimal repairing in interacting networks. Nat. Commun. 2016, 7, 1–10. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Lever, J.; Leemput, I.; Weinans, E.; Quax, R.; Dakos, V.; Nes, E.; Bascompte, J.; Scheffer, M. Foreseeing the future of mutualistic communities beyond collapse. Ecol. Lett. 2020, 23, 2–15. [Google Scholar] [CrossRef] [PubMed]
  5. Limiao, Z.; Guanwen, Z.; Daqing, L.; Hai-Jun, H.; Eugene, S.; Shlomo, H. Scale-free resilience of real traffic jams. Proc. Natl. Acad. Sci. USA 2019, 116, 8673–8678. [Google Scholar]
  6. Barlow, R.; Proshan, F. Mathematical Theory of Reliability; John Wiley and Sons: New York, NY, USA, 1965. [Google Scholar]
  7. Ryabinin, I.A. Reliability and Safety of Structurally Complex Systems; St. Petersburg University Press: Saint-Petersburg, Russia, 2007. (In Russian) [Google Scholar]
  8. Solojentsev, E.D. Scenario Logic and Probabilistic Management of Risk in Business and Engeneering; Springer: Berlin/Heidelberg, Germany, 2004. [Google Scholar]
  9. Gertsbakh, I. Statistical Reliability Theory; Marcel Dekker: New York, NY, USA, 1989. [Google Scholar]
  10. Rocchi, P. Reliability Is A New Science. Gnedenko Was Right; Springer: Berlin/Heidelberg, Germany, 2017. [Google Scholar]
  11. Zolotarev, V. Modern Theory of Summation of Random Variables; VSP: Utrecht, The Netherlands, 1997. [Google Scholar]
  12. Dudin, A.; Nazarov, A. On a tandem queue with retrials and losses and state dependent arrival, service and retrial rates. Int. J. Oper. Res. 2017, 29, 170–182. [Google Scholar] [CrossRef]
  13. Artalejo, J.; Gomez-Corral, A. Retrial Queueing Systems. A Computational Approach; Springer: Berlin/Heidelberg, Germany, 2008. [Google Scholar]
  14. Nazarov, A.; Moiseeva, E. Investigation of the RQ-system MMPP|M|1 by the method of asymptotic analysis in the condition of a large load. Izv. Tomsk. Polytech. Univ. 2013, 322, 19–23. (In Russian) [Google Scholar]
  15. Information Technologies and Mathematical Modelling-Queueing Theory and Applications. In Proceedings of the 17th International Conference, ITMM 2018, Named After A.F. Terpugov, and 12th Workshop on Retrial Queues and Related Topics, WRQ 2018, Tomsk, Russia, 10–15 September 2018; Volume 912.
  16. Borovkov, A. Asymptotic Methods in Queueing Theory; John Wiley and Sons: New York, NY, USA, 1984. [Google Scholar]
  17. Basharin, G.P.; Gaidamaka, Y.V.; Samouylov, K.E. Mathematical Theory of Teletraffic and Its Applications to the Analysis of Multiservice Communication of Next Generation Networks. Autom. Control. Comput. Sci. 2013, 47, 62–69. [Google Scholar] [CrossRef]
  18. Tsitsiashvili, G.S.; Osipova, M.A.; Samoulov, K.E.; Gaidamaka, Y.V. Synergetic effects in multiserver loss systems. In Proceedings of the VIII Moscow International Conference on Operations Research for the Centenary of Yu. B. Hermeyer at the Moscow State University and the VC RAS, Moscow, Russia, 17–21 October 2016; Volume I, pp. 350–355. [Google Scholar]
  19. Gnedenko, B.V.; Kovalenko, I.N. Introduction to Queuing Theory; Nauka: Moscow, Russia, 1966. (In Russian) [Google Scholar]
  20. Ivchenko, G.I.; Kashtanov, V.A.; Kovalenko, I.N. Queuing Theory; Visshaya Shkola: Moscow, Russia, 1982. (In Russian) [Google Scholar]
  21. Naumov, V.A. On the behavior of the parameters of the Equivalent Random Theory method at low load. In Numerical Methods and Informatics; RUDN Publisher: Moscow, Russia, 1988. (In Russian) [Google Scholar]
  22. Geritz, S.; Kisdi, E. On the mechanistic underpinning of discrete-time population models with complex dynamics. J. Theor. Biol. 2004, 228, 261–269. [Google Scholar] [CrossRef] [PubMed]
  23. Tsitsiashvili, G.; Osipova, M. Synergetic effects for number of busy servers in multiserver queuing systems. Commun. Comput. Inf. Sci. Ser. 2015, 564, 404–414. [Google Scholar]
  24. Heath, D.; Resnick, S.; Samorodnitsky, G. Heavy tails and long range dependence in on/off processes and associated fluid models. Math. Oper. Res. 1998, 23, 145–165. [Google Scholar] [CrossRef] [Green Version]
  25. Mikosch, T.; Resnick, S.; Rootzen, H.; Stegeman, A. Is network traffic approximated by stable Levy motion or fractional Brownian motion? Ann. Appl. Probab. 2002, 12, 23–68. [Google Scholar] [CrossRef]
  26. Tsitsiashvili, G.; Markova, N. Synergistic effects in a multi-channel queuing system with an Erlangian input flow. Bull. Pac. State Univ. 2015, 4, 17–22. (In Russian) [Google Scholar]
  27. Jagerman, D.L. Some Properties of the Erlang Loss Function. Bell Syst. Tech. J. 1974, 53, 525–551. [Google Scholar] [CrossRef]
  28. Mitra, D.; Weiss, A. The Transient Behavior in Erlang’s Model for Large Trunk Groups and Various Traffic Conditions. Proc. 1988 Int. Teletraffic Congr. 1988, 26, 223. [Google Scholar] [CrossRef] [Green Version]
  29. Tsitsiashvili, G.S.; Osipova, M.A. Phase Transitions in Multiserver Queuing Systems. Inf. Technol. Math. Model. Queueing Theory Appl. 2016, 638, 341–353. [Google Scholar]
  30. Sharkovskiy, A.; Sharkovskiy, A.N. Difference Equations and Population Dynamics. Preprint 82.18; Institute of mathematics of the Academy of Sciences of the Ukrainian SSR: Kiev, Ukraine, 1982. (In Russian) [Google Scholar]
  31. Leonov, G.; Kuznetsov, N.; Korzhemanova, N.; Kusakin, D. Lyapunov dimension formula for the global attractor of the Lorenz system. Commun. Nonlinear Sci. Numer. Simul. 2016, 41, 84–103. [Google Scholar] [CrossRef] [Green Version]
Figure 1. n-multiple block redundancy (top), split redundancy (bottom) of a chain of length m .
Figure 1. n-multiple block redundancy (top), split redundancy (bottom) of a chain of length m .
Mathematics 09 01396 g001
Figure 2. n isolated M | M | 1 | 0 systems (left), aggregated M | M | n | 0 system (right). b ^ n .
Figure 2. n isolated M | M | 1 | 0 systems (left), aggregated M | M | n | 0 system (right). b ^ n .
Mathematics 09 01396 g002
Figure 3. Frequency histogram for b ^ n .
Figure 3. Frequency histogram for b ^ n .
Mathematics 09 01396 g003
Table 1. Meanings of n * ( m , δ ) , n * ( m , δ ) for p = 0 , 7 , δ = 0 , 1 .
Table 1. Meanings of n * ( m , δ ) , n * ( m , δ ) for p = 0 , 7 , δ = 0 , 1 .
m1234536789101112131415
n * ( m , δ ) 2469131927395681116166237339484
n * ( m , δ ) 233444444444555
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tsitsiashvili, G. Study of Synergistic Effects in Complex Stochastic Systems. Mathematics 2021, 9, 1396. https://doi.org/10.3390/math9121396

AMA Style

Tsitsiashvili G. Study of Synergistic Effects in Complex Stochastic Systems. Mathematics. 2021; 9(12):1396. https://doi.org/10.3390/math9121396

Chicago/Turabian Style

Tsitsiashvili, Gurami. 2021. "Study of Synergistic Effects in Complex Stochastic Systems" Mathematics 9, no. 12: 1396. https://doi.org/10.3390/math9121396

APA Style

Tsitsiashvili, G. (2021). Study of Synergistic Effects in Complex Stochastic Systems. Mathematics, 9(12), 1396. https://doi.org/10.3390/math9121396

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop