Yu.A. Bychkov, S.V. Shcherbakov
The task of studying the qualitative features of the dynamics of mathematical models of nonlinear nonautonomous systems with lumped parameters is extremely relevant. The development of the theory of nonlinear phenomena and the obtained applied results generate a variety of methods for solving this problem [1-3].
In the general case, the dynamics of mathematical models of systems of a selected class (hereinafter referred to as “systems”) are described by the following ordinary nonlinear integral-differential equation with non-stationary coefficients:
A ( D ) x ( t ) = G ( D ) f ( t ) + H ( x , f , t ) q q u a d ( 1 )
Where D - operator of generalized differentiation with respect to independent variable t ; A ( D ) - square, order L x matrix with polynomial from D and D - 1 elements where D - 1 - integration operator with a variable upper limit t ;
G (D) is a rectangular matrix of size L x t i m e s L f with polynomial from D and D - 1 elements; x ( t ) and f ( t ) - matrix columns of the coordinates of the system (the desired solutions) and external influences on it, respectively; H(x,f,t) - a column matrix with rows in the form of sums of products whose factors are non-stationary coefficients, as well as classical derivatives of any order and integrals of any multiplicity, starting with zero, from the desired solutions and external influences, in arbitrary fractional rational degrees.
Calculation of the dynamics of nonlinear nonautonomous systems under given pre-initial conditions Dnx r(0−)=Dnx r(t −0) , r=1,2,...,L x , 0−=t− 0 , n in[0;M−1] comes down to finding all existing in a given study interval [t 0;T] solutions of equation (1) and is associated with a number of algorithmic and mathematical problems. The essence of these problems determines the objective variety of qualitative features of the desired solutions.
In the general case, the nature of the dynamics of the desired solutions of equation (1) is unpredictable, this is due to the dependence of the manifestation forms of the dominant nonlinear and non-stationary properties of the system on the selected initial conditions, the ratio of parameters and the type of its characteristics. In the general case, the indicators of the desired solutions of the equation (1) irregular, reducing the dynamics of the system to the so-called deterministic chaos [2-4]. Describing its dynamics, we note that in the desired solutions of the equation (1) gaps of the first kind, including differentiable ones, and gaps of the second kind are possible. The specificity of these decisions is that they, as a rule, are associated with "rigidity", i.e. with alternating intervals of fast and slow changes in their dynamic indicators, and within such intervals transitions from local stability to instability and vice versa are possible. A direct consequence of local instability is the effect of “mixing” solutions, which greatly complicates the problem of numerical calculation of the dynamics of a nonlinear nonautonomous system and analysis of its qualitative features [5]. The possible qualitative diversity of the features of the dynamics of a nonlinear nonautonomous system, as a result of the complex synthesis of its coordinate-parametric interconnections, allows only an approximate solution of the equation (1) using one or another numerical method [6-8]. There are many numerical methods and their calculation schemes are constantly being improved. However, none of them allows you to establish in a closed form any relationship or pattern between the behavior of the desired solution of the equation (1) and information indicators of this equation. The task of identifying such relationships and organizing on their basis a subsequent analysis of the qualitative features and properties of solutions of the equation (1) important in both theoretical and applied terms. Its solution means the next step towards understanding the content of cause-effect relationships in nonlinear and non-stationary phenomena.
The article provides a solution to the problem based on the use of the nonlinear nonautonomous system of eigenvalues as information indicators of the dynamics corresponding to its equation (1) Jacobi functional matrix and decomposition by such eigenvalues of solutions of this equation. As a computational basis, an analytical-numerical method was used, which according to its computational characteristics corresponds to the task [9, 10]. An illustrative example is provided.
The composition of information indicators of the qualitative features of the desired solutions of the equation (1) It is advisable to include those that determine the interaction of the parameters of this equation and the non-stationary indicators of the solutions themselves. In a linear stationary system whose dynamics is described by the equation (1) at H(x,f,t)=const , such information indicators are the roots lambda n,n=1,2,...,N the following characteristic equation:
A(p)=0 qquad(2)
The nature and numerical values of the roots of the equation (2) , determine all the features and properties of the dynamics of a linear autonomous system in the transition process, including stability and "rigidity" [1,7,8]. Since external influences unambiguously regulate the nature of the steady state, the dynamics of a linear autonomous system is completely predictable in a semi-infinite time interval.
The directional selection of such information-intensive and fairly simple in determining the dynamics of nonlinear non-autonomous systems is not yet possible. The urgent needs of the practice, however, require already now the selection of indicators of such systems that are available for calculation and are acceptable in terms of the information content. If we follow the same path as for the linear case, then it’s logical to consider the roots as the necessary information indicators lambda n,n=1,2,...,N the equations (2) , the formation of which will be associated with the calculation of the determinant A(p) matrices A(D) selected linear part of the equation (1) . The interpretation of the semantic content of the resulting information indicators lambda n,n=1,2,...,N ambiguous because the interrelation of these essentially selected stationary indicators with the qualitative features of the dynamics of a nonlinear nonautonomous system does not lend itself to any meaningful identification, even staged. Also, when in the equation (1) matrix A(D) zero or highlighted by non-linear part H(x,f,t) this equation, then the roots of the equation (2) generally absent or their existence, and therefore the content, is conditionally subjective [9].
Thus, the target attribute of the selection of acceptable information indicators of the dynamics of a nonlinear non-autonomous system is the condition of their existence and the mathematical relationship with the numerical indicators of the desired solution of the equation (1) . This condition can be satisfied if in the equation (1) linear part matrix A(D) in a certain way to select or supplement at the expense of the matrix H(x,f,t) nonlinear part [9,10]. This selection or addition of the matrix A(D) due to the matrix H(x,f,t) perhaps in an ambiguous way. All transformations performed during this process must be equivalent, as a result of which the dynamics of the system will remain unchanged. As a result of the selection or addition of the matrix A(D) due to the matrix H(x,f,t) the necessary mathematical relationship of the roots will be achieved lambda n,n=1,2,...,N the equations (2) with the numerical values of the desired solutions of the equation (1) at discrete points in time, determined by the step of numerical calculation of the dynamics of the system. Thus, as a result of these transformations of the equation (1) the roots lambda n,n=1,2,...,N the equations (2) will acquire a functional dependence on time, both in relation to numerical indicators and qualitative characteristics.
Nonstationary root properties lambda n,n=1,2,...,N the equations (2) in the absence of a target criterion for the implementation of the corresponding matrix transformations A(D) leading to the appearance of such properties cannot by themselves guarantee the necessary informative characteristics of these roots. The desired study of the qualitative features of the dynamics of nonlinear objects dictates the choice as the target criterion for the formation or transformation of an existing matrix A(D) the equations (1) the condition of its coincidence with the Jacobi functional matrix corresponding to this equation [4–8]. The algorithm for solving this problem is described in [10]. As a result of the formation or transformation of the matrix A(D) due to the matrix H(x,f,t) to the Jacobi functional matrix of non-stationary roots lambda n,n=1,2,...,N the equations (2) acquire the status of eigenvalues of the Jacobi matrix.
Nonstationary eigenvalues lambda n,n=1,2,...,N Jacobi functional matrix of the equation (1) at discrete points in time, determined by the step of its numerical solution, have a certain well-known semantic content [4-8]. Being a local characteristic of the dynamics of a nonlinear nonautonomous system, they in the infinitely small neighborhood of a point corresponding to a discrete point in time t characterize the stability and rate of change of dynamic indicators of solutions of its dynamic equation (1) . In the neighborhood of the singular points of the equation (1) own numbers lambda n,n=1,2,...,N Jacobi’s functional matrix allows one to isolate areas of attraction and repulsion of phase trajectories on the phase plane, thus giving a forecast for the limiting states of existing solutions of this equation [6–8]. Eigenvalues corresponding to discrete points in time lambda n,n=1,2,...,N the Jacobi functional matrix, characterizing and determining the qualitative features and properties of the dynamics of nonlinear nonautonomous systems, determine the relevance of the task of identifying mathematical features and information indicators of such a relationship. The solution of this problem is associated with the establishment of a mathematical relationship between the selected non-stationary information indicators lambda n,n=1,2,...,N the dynamics of a nonlinear nonautonomous system and the parameters of the calculation scheme of the numerical method used to calculate this dynamics. A prerequisite for achieving this goal is the presence in the calculation scheme of the selected method of the corresponding analytical part. The computational algorithm of this analytical part of the method must be consistent with the peculiarities of calculating the eigenvalues of the Jacobi functional matrix, and also provide a mathematical relationship between these numbers and the dynamic coordinates of the coordinates of a nonlinear non-autonomous system. The one-step analytical-numerical method of variable order [9, 10] fully meets the formulated requirements.
The motivation for the formation of the analytical-numerical method was the desire to unify the solution of the equation (1) using only identical transformations. Hence the description of the desired solutions by Taylor series and the subsequent application of the Laplace integral for the algebraization of the problem.
Analytical-numerical method for calculating the dynamics of the selected equation (1) class of systems in a given study interval [t 0;T] consists of two parts: analytical and numerical.
The analytical part of the method is based on apparatuses of generalized functions, the Laplace transform, and function-power series. The procedure of the analytical part precedes the execution of each successive calculation step and reduces to the following. First, by expanding the functions describing external actions and non-stationary parameters of the system, into their corresponding power series, and also formally describing the regular components of the desired solutions that are included in the matrix H(x,f,t) the equations (1) corresponding to their power series, we transform the right-hand side of this equation into a matrix-column, the rows of which are now power series. The coefficients of these power series are generally unknown, since they are expressed in terms of known coefficients of power series for external influences and non-stationary parameters of the system, and unknown coefficients of power series for the regular components of some of the desired solutions. The performed operation gives the original equation (1) to the form required for the subsequent application of the integral Laplace transform. Transforming Laplace's Reformed Equation (1) and solving the resulting algebraic equation according to Cramer’s rule, for the image X l(p) the desired solution x l(t),l in[1,L x] we get the following expression:
X l(p)= fracBl l(p)A(p)= frac sum limits i=0inftyBl l.N+J l−ipN+J l−i sum limits Ni=0A ipi, qquad(3)
Where N in mathbbN;J l in mathbbZ
Fractional rational function, described by expression (3), in the general case, when J l geq0 is incorrect and allows its presentation as a sum of a whole rational X −l(p) and correct fractional rational X +l(p) functions. The subsequent decomposition of the correct fractional rational function X +l(p) the Laurent series in the neighborhood of an infinitely remote point transforms the expression (3) in the following way:
X l(p)=X −l(p)+X +l(p)= sum −J lj=0S ljp−j+ frac sum limits inftyi=1B lN−ipNi sum limits Ni=0A ipi== sum −J lj=0S ljp−j+ frac sum limits inftyi=0B l.N−i−1pN−i sum limits Ni=0A ipi frac1p== sum −J lj=0S l.jp−j+ sum inftyi=0R l.ip−(i+1). qquad(4)
S l.−J l= fracBl l.N+J lA N;
S l.−J l+j= fracBl l.N+J l+j− sum limits j−1k=0S l.−J l+kA N−j+kA N
B l.N−i=Bl l.N−i− sum N−ik=0S l.−N+i+kA k
\ begin {align} & {R \ _ {l.0} = \ frac {B \ _ {lN-1}} {A \ _ {N}};} \\ & {R \ _ {li} = \ frac {B \ _ {lN-1-i} - \ sum \ limits \ _ {k = 0} ^ {i-1} R \ _ {lk} \; A \ _ {N-i + k}} {A \ _ {N}},} \ end {align} \ qquad (5)
x l(t)=x −l(t)+x +l(t)= sum −J lj=0S lj rm delta j(t)+ sum inftyi=0R l.iti/i!, qquad(6)
h=q rm tau, qquad(7)
x +l(t k; rmI l)−| Deltax +l(t k; rmI l) rm| lex +l(t k) lex +l(t k; rmI l)+| Deltax +l(t k; rmI l) rm|. qquad(8)
R l.i= sum Nn=1r l.n lambda in,i in ge left| bfZ right|. qquad(9)
Coefficients r l.n in this case, we calculate by the formula:
r l.n= frac sum limits im=1B ∗l.N−1−m lambda −mnN+ sum limits N−1m=1(N−m)A ∗Nm lambda −mn. qquad(10)
\ begin {align} & {B \ _ {l.N-1-m} ^ {\ *} = \ frac {B \ _ {l.N-1-m}} {A \ _ {N}} \; ;} \\\\ & {A \ _ {r} ^ {\ *} = \ frac {-A \ _ {r}} {A \ _ {N}} \; .} \ end {align}
x +l(t)= sum inftyi=0 fracR l.it ii!= sum inftyi=0 sum N mn=1 fracR [n]lit ii!, Qquad(11)
x +l(t)= sum N mn=1x [n]l(t). qquad(12)
x [n]l(t)= sum inftyi=0 fracR [n]l.it ii!. qquad(13)
interrelated equality system (11)−(13) Is the result of the directed transformation of the incoming equality (6) source description of the regular component of the solution x +l(t) . According to it, at each calculation step, this component of the solution can be decomposed by poles rm lambda n, rmn=1,2,...,N rm her images X +l(p) . Decomposition, as one would expect, is not absolute, because through recurrently calculated coefficients B ∗l.N−1−m,m in left| bfZ right| included in the formula (10) , for a nonlinear nonstationary system, the relationship between all elements that always dominates in the formation of its dynamics is always preserved (13) representation (12) . Decomposition of the regular component x +l(t) the desired solution x l(t),l in[1;L x] described by a system of equalities (11)−(13) , was successfully tested when solving a number of special problems related to the formalization of the procedure for studying the existence and uniqueness of a solution, choosing a calculation step, as well as optimizing the computational costs associated with estimating the absolute local error of such a calculation [9].
Summing up, we formulate the following algorithm for studying the qualitative features of the dynamics of nonlinear nonautonomous systems. First, perform the equivalent transformation of the original equation (1) system dynamics to mind when the matrix A(D) coincides with the Jacobi functional matrix, then the poles rm lambda n, rmn=1,2,...,N rm Images X +l(p) match the eigenvalues of this matrix. Subsequent decomposition of the regular component of the solution x +l(t) on these numbers will provide the possibility of organizing a design scheme for the necessary research. In the role of non-stationary information indicators of the qualitative features of the system dynamics, we consider the functions describing the change in a given study interval. rm[t 0;T] proper numbers rm lambda n, rmn=1,2,...,N rm Jacobi functional matrix and elements (13) decomposition (12) regular component of the solution x +l(t) . The proposed decomposition scheme for the desired solutions of the equation (1) over eigenvalues of the Jacobi functional matrix corresponding to this equation is extraordinary in and of itself. This determines the novelty of the expected results, their interpretation and subsequent application prospects. We will designate only some of them, since the completeness of judgments about such a problem is a matter of time.
Signs of real parts of non-stationary eigenvalues rm lambda n, rmn=1,2,...,N rm Jacobi functional matrix at each discrete point in time t=t k determine the possibility of identifying the boundaries of the intervals of stability and instability in relation to each of the elements (13) decomposition (12) . Such information is primary to identify causal relationships between the sustainability conditions of the regular component of the solution. x +l(t),l in[1;L x] and similar conditions for the elements themselves (13) this component. Conditions of stability or instability of the regular component x +l(t) the desired solution x l(t) this will be a consequence of the existence of non-stationary eigenvalues rm lambda n, rmn=1,2,...,N rm Jacobi's functional matrix of one or several dominant ones, and also on the features of the interaction of these dominant numbers in the study interval, including the variation of their nature.
Dynamic indicators of decomposition elements (13) in conjunction with nonstationary eigenvalues rm lambda n, rmn=1,2,...,N rm Jacobi functional matrix determine the initial position for a different interpretation of the result of the calculation of the dynamics of a nonlinear non-autonomous system described by double inequality (8) . Mathematical interrelation and correlation in the course of calculation of indicators of double inequality (8) and double inequalities like inequalities (8) but for corresponding elements of decomposition (13) with the irregularity of the dynamic properties of the system, serve as the basis for the selection of characteristic ratios and combinations between the elements of the decomposition (13) for which the desired solutions of the equation (1) most sensitive to changes in parameters and external influences of the system. The analysis of the results obtained in this case makes it possible to isolate the conditions that determine the occurrence of bifurcation, or the ranges of initial conditions leading to the “mixing” of phase trajectories [5].
Qualitative features and characteristic properties of the functions describing the change in a given study interval rm[t 0;T] nonstationary eigenvalues rm lambda n, rmn=1,2,...,N rm Jacobi’s functional matrix, including the continuity and monotony of these functions, the presence of extremum points, the variability of the nature of such numbers and the sign of their real parts, undoubtedly reflect the essence of the internal causal relationships of the equation (1) . As a result, the features and properties of these functions are highly adaptive indicators of the qualitative features of the dynamics of a nonlinear nonautonomous system. The diversity of the content spectrum and the variability of the manifestations of such indicators of the system’s dynamics are still to be determined and this constitutes an independent task, the solution of which may serve as a starting point for a more complete understanding of the complex and ambiguous behavior of the nonlinear object of study. Separating and classifying the relationships between these functional indicators and the qualitative features of a nonlinear nonautonomous system should lead to the emergence of new computational algorithms for analyzing nonlinear and non-stationary phenomena, including determining the conditions for the existence of regular dynamics or the transition to deterministic chaos [2-5].
Described by the system of equalities (11)−(13) decomposition of the regular components of the desired solutions of the equation (1) introduces into consideration and activates as an independent computational toolkit of analysis of the dynamics of a nonlinear nonautonomous system, the eigenvalues of the Jacobi functional matrix. They generate features and dynamic properties of decomposition elements. (13) . As shown in the following example, the analysis of such features and dynamic properties of the decomposition elements (13) leads to the selection of new conditions for the existence of one of the characteristic properties of nonlinear dynamics - "stiffness" [1,4,5].
Summarizing what has been said, we note that there are many variants of manifestations of qualitative features and characteristic properties of nonlinear dynamics, they are diverse in terms of the conditions of occurrence and the forms of manifestations. Irregularity, as a property of such dynamics, is rather the norm than the exception [2]. The reasons for such a complex and unpredictable dynamics of a nonlinear non-autonomous system are laid down in the structure and composition of the internal parametric relationships of the dynamic equation. (1) . One of the components of the complex indicator of such interrelations is undoubtedly eigenvalues rm lambda n, rmn=1,2,...,N rm Jacobi functional matrix. Analysis of the qualitative features of the dynamics of a nonlinear nonautonomous system based on the decomposition of the desired solutions of its dynamic equation (1) into components, each of which corresponds to one of the eigenvalues rm lambda n, rmn=1,2,...,N rm Jacobi’s functional matrix, while maintaining the logical and mathematical relationship between these components, reflects the essence of the proposed approach to conducting relevant research and a new methodology.
Example
Investigation of the qualitative features of the dynamics of a non-linear autonomous system based on the decomposition of solutions of its dynamic equation by eigenvalues of the Jacobi functional matrix corresponding to this equation
In the shape of (1) the dynamic equation of the considered nonlinear autonomous system has the following form:
$$ display $$ \ begin {equation} \ left \ | \ begin {array} {cc} {a \ _ {1.1} ^ {[1]} D + a \ _ {1.1} ^ {[0]}} & {0} \\\\ {a \ _ {2.1 } ^ {[0]}} & {a \ _ {2.2} ^ {[1]} D} \ end {array} \, \ right \ | \, \ left \ | \, \ begin {array} {c} {x \ _ {1} (t)} \\\\ {x \ _2 (t)} \ end {array} \ right \ | = \ left \ | \ begin {array} {c} {g \ _ {1.1} ^ {[0]}} \\\\ {0} \ end {array} \ right \ | f (t) + \ left \ | \ begin { array} {c} {h \ _ {1.1} x \ _ {1} ^ 2 (t) x \ _2 (t)} \\\\ {h \ _ {2.1} x \ _ {1} ^ {2 } (t) x \ _ {2} (t)} \ end {array} \ right \ |, \ qquad (14) \ end {equation} $$ display $$
Where a left[1 right]1.1=a left[1 right]2.2=h kern1pt1.1=1;h kern1pt2.1=−1;f left(t right)= rm delta kern1pt1 left(t right);g kern1pt left[0 right]1.1= textitA;a kern1pt left[0 right]1.1=B+1;a kern1pt left[0 right]2.1=−B.
After substitution of the values of the parameters specified in the explication to the equation (14) , obtained the following unified form of this equation:
$$ display $$ \ begin {equation} \ left \ | \ begin {array} {cc} {D {\ rm + B + 1}} & {0} \\\\ {-B} & {D} \ end {array} \ right \ | \, \ left \ | \ begin {array} {c} {x \ _ {1} (t)} \\\\ {x \ _ {2} (t)} \ end {array} \ right \ | = \ left \ | \ begin {array} {c} {A} \\\\ {0} \ end {array} \ right \ | \ delta \ _ {1} (t) + \ left \ | \ begin {array} {c} {x \ _ {1} (t) ^ {2} x \ _ {2} (t)} \\\\ {-x \ _ {1} (t) ^ {2 } x \ _ {2} (t)} \ end {array} \ right \ | . \ end {equation} \ qquad (15) $$ display $$
Where R_1.0 and R_2.0 - x_r(t),r=1,2 . (16) - , (6) :
x_r(t)=x_r+(t)=∑_i=0∞R_r.itii ! ,r=1,2,( 17 )
Formed for the desired solutions x_r(t),r=1,2 (16) (17) , . - [0; 10] . ε_r(h)=1⋅10−5 x_1(0)=1,x_2(0)=2,5 , A=2,B=6 .1. , [4]. «», .. (2) . , .2 , , «» , . [3-4], [8-9], , , , , 100 . t=4 , , , [4-8] , , . «», , . , , . λ_n,n =1,2 , , , R_1.0 and R_2.0 . , (9) , (10) , R_l.i . , , «» - (16) . R_l.i x_r(t),r=1,2 , , , λ_n,n=1,2 , , «» . , , «», , , .
Fig.1. (16) A = 2, B = 6 (16) .
«» . 2-4. (11) - (13) , N_m=N=2 and l=1,2 x_l(t)=x_l+(t) (16) λ_n,n=1,2 . , , .2 λ_n,n=1,2 (16) . , «», ,
()
()
Fig.2. () () λ_n,n=1,2 A= 2, B = 6
Fig.3. () () (13), l=1,n=1,2 x_1(t) (16) A = 2, B = 6
Fig.4. () () (13), l=2,n=1,2 x_2(t) (16) A = 2, B = 6
[0-10] λ_n,n=1,2 , λ_n(t) . , (16) . , .2, λ_n(t) , . , (16) , . , , , , λ_n(t),n=1,2 . (16) : =A2+1 .
λ_n,n=1,2 , , (9) , (10) (11)−(13) , x_l[n](t) , n=1,2,l=1,2 x_l(t) (16) . [0 -2,5] [4 -6,5], , , .1,3,4 , , x_l[n](t) , n=1,2,l=1,2 x_l(t) , (12) , N=2,l=1,2 , x_l+(t_k;I_l) . x_l[n](t) , n=1,2,l=1,2 x_l(t) , x_l+(t_k;I_l) . , (16) x_l+(t_k;I_l) - .
[3 — 4] [7 — 8], , .1,3,4 , (16) , , x_l[n](t) , n=1,2,l=1,2 , , λ_n,n=1,2 , . x_l[n](t) , n=1,2,l=1,2 x_l(t) λ_n,n=1,2 , .2 . x_l[n](t) , n=1,2,l=1,2 x_l(t) , . -. , x_1(t) , .3 , x_l[n](t) , n=1,2,l=1,2 , .4, , x_2[1](t) , x_2[2](t) . , (16) , x_l[n](t) , n=1,2,l=1,2 x_l(t),l=1,2 .
(16) , , , .5-8. , (16)A=2,B=5.
Fig.5. (16) A = 2, B = 5
, , A=2,B=6.
.5 (16) at A=2,B=5 , =A2+1 , >(A2+1) , : . , .6-8, λ_n,n=1,2 (13) at l=1,2 and n=1,2 x_l[n](t) x_l(t) .
.6 λ_n,n=1,2 A=2 and B=5 . .6 , t=3,5 λ_n,n=1,2 . .2 , , . , , x_l(t),l=1,2 (16) (A=2,B=6) (A=2,B=5) , , , λ_n,n=1,2 , .
.5 , (16) «». (A=2,B=6) t , λ_n,n=1,2 , . (16) , A=2,B=5 , - . , «» (16) , , λ_n,n=1,2 , .
x_l[n](t) , n=1,2,l=1,2 x_l(t) (16) , λ_n,n=1,2 . , .3,4 .7,8 , , (A=2,B=6) (A=2,B=5) x_l[n](t) , n=1,2,l=1,2 x_l(t) . , x_l[n](t) , n=1,2,l=1,2 x_l(t),l=1,2 . .7,8 . , , x_l(t),l=1,2 (16) (A=2,B=6) ( A = 2 , B = 5 ) , is the exclusion of any dominance among the componentsx _l [ n ] ( t ) , n = 1 , 2 , l = 1 , 2 in the formation of dynamic indicators of the desired solutionsx _ l ( t ) .
Source: https://habr.com/ru/post/321212/
All Articles