URAL MATHEMATICAL JOURNAL, Vol. 9, No. 1, 2023, pp. 93-103
DOI: 10.15826/umj.2023.1.007
ON THE PROPERTIES OF THE SET OF TRAJECTORIES OF THE NONLINEAR CONTROL SYSTEM WITH QUADRATIC INTEGRAL CONSTRAINT ON THE CONTROL FUNCTIONS
Anar HuseyW, Nesir Huseyin^
Sivas Cumhuriyet University, 58140 Merkez, Sivas, Turkey
[email protected] tt [email protected]
Abstract: In this paper the control system described by a nonlinear differential equation is studied. It is assumed that the control functions have a quadratic integral constraint, more precisely, the admissible control functions are chosen from the ellipsoid of the space L2 ([to, Rm). Different properties of the set of trajectories are investigated. It is proved that a small perturbation of the set of control functions causes also appropriate small perturbation of the set of trajectories. It is also shown that the set of trajectories has a small change if along with the integral constraint on the control functions, a sufficiently large norm type geometric constraint on the control functions is introduced. It is established that every trajectory is robust with respect to the fast consumption of the remaining control resource, and hence every trajectory of the system can be approximated by a trajectory generated by full consumption of the total control resource.
Keywords: Nonlinear control system, Quadratic integral constraint, Set of trajectories, Robustness.
1. Introduction
The control systems described by nonlinear differential equations are investigated in a vast number of papers. Depending on the character of the control efforts the control systems are classified as a) the control systems with geometric constraint on the control functions; b) the control systems with integral constraint on the control functions; and c) the control systems with mixed constraints on the control functions which include both the geometric and the integral constraints on the control functions. The geometric constraints on the control functions appear in the case when the control resource is not exhausted by consumption. But, if the control resource is exhausted by consumption, say as energy, food, fuel, finance, etc., then the integral constraints on the control functions is inevitable (see, e.g., [1, 2, 9, 12, 15, 16]). For example, the behaviour of the flying objects with rapidly changing mass is described as a control system with integral constraint on the control functions (see, e.g., [2, 12]).
One of the important notions of the control systems theory is the set of trajectories and attainable set concepts. Attainable set of the system at the given instant of time consists of points to which arrive the trajectories of the system and can be defined as a section of the set of trajectories at the given instant of time. Different topological properties and approximate construction methods of the set of trajectories described by various types of the integral and differential equations, where the control functions have integral constraints, are considered in papers [4-8, 11, 13, 14]. In papers [4, 5, 11, 14] the compactness, closedness, path-connectedness properties and approximate construction methods of the set of trajectories and attainable sets of the control systems which are affine with respect to the control vector are discussed. In papers [6-8, 13] the same problems are investigated for nonlinear control systems. In presented paper the properties of the set of trajecto-
ries of the nonlinear control systems are studied where the admissible control functions are chosen from the ellipsoid of the space L2.
The paper is organized as follows. In Section 2, the basic conditions which have to satisfy the system's equation are formulated and preliminary properties of the system's trajectories are given. In Section 3 it is shown that introduction of the sufficiently large norm type constraint along with integral constraint and a small perturbation of the given ellipsoid, which characterizes the integral constraint, induce a small change of the set of trajectories (Theorem 2). A perturbation evaluation for the set of trajectories is presented. In Section 4 it is proved that every trajectory is robust with respect to the fast and full consumption of remaining control resource (Proposition 7). Applying this result it is proved that every trajectory can be approximated by the trajectory generated by full consumption of the total control resource (Theorem 3).
2. The system's dynamics
Consider control system described by nonlinear ordinary differential equation
£fc(t) = f (t,x(t),u(t)), x(to)= xo (2.1)
where x(t) € Rn is the phase state vector, u(t) € Rm is the control vector, t € [t0,0] is the time.
Let B(-) : [t0,0] ^ Rmxm be a continuous matrix function and B(t) be a positive definite m x m matrix for every t € [t0,0]. For given e € [0,1] and a > 0 we denote
re
)u(t),u(t))dt < 1 + e
'to '
Ua = {«(•) € Ue : ||u(t)|| < a for almost all t € [to, 0]} ,
r6
(B (t)u(t),u
t0
U = {u(0 € ¿2([io,0];Rm) : f <B(t)u(t),u(t))dt ua = MO € U : ||u(t)|| < a for almost all t € [ U* = {u(-) € L2([to,0];Rm) : jf<B(t)u(t),u(t))dt = l},
where L2 ([t0,0]; Rm) is the space of Lebesgue measurable functions u(-) : [t0,0] ^ Rm such that
r0 \ !/2
12 < Ilu(')ll2 =
||u(.)|2 < IK)Il2 = (I l|u(t)|2ds)
denotes the Euclidean norm, (■, ■) stands for inner product.
Proposition 1. The sets U£ and Ua are bounded, closed and convex subsets of the space L2 ([t0,0]; Rm). The set UQ is bounded and closed subset of the space L2 ([t0,0]; Rm).
It is not difficult to show that there exists c* > 0 such that the inequality
KOH2 < c* (2.2)
is satisfied for every u(-) € U and e € [0, l].
It is assumed that the function f (■, ■, ■) : [t0,9] xRn xRm ^ Rn satisfies the following conditions:
2.A. The function f (■, ■, ■) : [t0, 9] x Rn x Rm ^ Rn is continuous.
2.B. For every bounded set D C [t0,9] x Rn there exist 71 = 71(D) > 0, y2 = y2(D) > 0 and 73 = 73(D) > 0 such that the inequality
If (t,xi,ui) - f (t,X2,U2)|| < [71 + 72(|ui| + ||«2||)] ||x - X2 || + 73|K - U2I
is satisfied for every (t,xi,ui) € D x Rm and (t, x2,u2) € D x Rm.
2.C. There exists k > 0 such that the inequality
11/(t,x,u)|| < K (HxH +1)(HuH + 1) is held for every (t, x, u) € [t0,0] x Rn x Rm.
If the function (t, x, u) — /(t, x, u) : [t0,0] x Rn x Rm — Rn is Lipschitz continuous with respect to (x,u), then the conditions 2.B and 2.C are satisfied.
Let us define the trajectory of the system (2.1) generated by a control function u*(-) € L2 ([t0,0]; Rm). An absolutely continuous function x*(-) : [t0,0] — Rn satisfying the equation x*(t) = / (t, x* (t), u*(t)) for almost all t € [t0,0] and initial condition x*(t0) = x0 is said to be a trajectory of the system (2.1) generated by the control function u*(-) € L2 ([t0, 0]; Rm). The sets of trajectories of the system (2.1) generated by all admissible control functions u(-) € U£, u(-) € and u(-) € U* are denoted by X£(t0,x0), Xa(t0,x0) and X*(t0,x0) respectively. It is obvious that the inclusions
Xa(t0,x0) C X£(t0,x0), X*(t0, x0) C X£(t0,x0)
are verified for every e € [0,1] and a > 0.
For fixed t € [t0, 0] we set
X£(t; t0, x0) = {x(t) € Rn : x(-) € Xe(t0, x0)} , (2.3)
X£a(t; t0, x0) = {x(t) € Rn : x(-) € X£a(t0, x0)} , (2.4)
X*(t; t0,x0) = {x(t) € Rn : x(-) € X0*(t0,x0)} . (2.5)
The sets X£(t; t0,x0), Xa(t; t0,x0) and X*(t; t0,x0) are called the attainable sets of the system (2.1) at the instant of time t, generated by all admissible control functions from the sets U , Ua and U0* respectively.
It is obvious that the attainable sets consist of points to which arrive the trajectories of the system (2.1) at the instant of time t.
By symbol C ([t0,0]; Rn) we denote the space of continuous functions x(-) : [t0,0] — Rn with norm
||x(-)||e = max{||x(t)|| : t € [t0,0]} ,
hn(■, ■) and he(■, ■) stand for the Hausdorff distance between the subsets of the spaces Rn and C ([t0,0]; Rn) respectively.
Let us formulate the propositions which will be used in following arguments.
Proposition 2. Each control function u(-) € L2 ([t0, 0]; Rm) generates unique trajectory of the system (2.1).
Denote
a*
- t0) + (0 - t0)1/2c* ■ exp k (0 - t0) + c*(0 - t0)1/2 (2.6)
where c* is defined by (2.2).
The following proposition characterizes boundedness of the set of trajectories.
Proposition 3. For every e € [0,1] and x(-) € X£(t0,x0) the inequality ||x(-)||e < a* holds.
Proof. Let us choose an arbitrary e > 0 and x(-) € X£(t0,x0), generated by the control function u(-) € U£. According to the Condition 2.C, inequality (2.2) and Cauchy-Schwarz inequality we have
||x(t)||< K f t(|x(T)|| + 1)(||u(T)|| + 1) dT
J to
< K i (||u(t)|| + 1)||x(T)|| dT + K [(0 - t0) + (0 - t0)1/2c*
to
K
for every t € [t0,0]. Applying Bellman-Gronwall inequality and Cauchy-Schwarz inequality and taking into consideration (2.2) and (2.6) we conclude from the last inequality
x(t)|| < k [(9 - to)+ (9 - to)i/2c*
re
K (||u(T)|| +l) dT J to
< K
(9 - to)+ (9 - to)i/2c:
■ exp k
■ exp
(9 - to) + c*(9 - to)i/2
a*
for every t € [to,0]. The inequality (2.7) completes the proof.
(2.7)
□
Let
^(5) = K(a* + l)(5 + c*5i/2) , 5 > 0. It is obvious that ^(5) ^ 0 as 5 ^ 0+.
(2.8)
Proposition 4. For every e € [0, l], x() € X£(to,xo), ti € [to,9] and t2 € [to,9] the inequality
is verified, and hence
|x(ti) - X(t2)| < ^(|ti - t2|)
hn(X£(ti; to,xo),X£(t2; to,xo)) < ^(|ti - t21)
where is defined by (2.8).
Proof. Without loss of generality let us assume that t2 > ti. Choose an arbitrary e > 0 and x() € X£(to,xo), generated by the control function u() € U£. According to the Condition 2.C, Proposition 3, (2.2) and (2.8) we have
||x(t2) - x(ti)|| < K f t(||x(T)|| + l)(|u(T)|| + l) dT < K(a* + l) / 2(||u(t)|| + l) dT
Jt1 Jt1
< K(a* + l) [(t2 - ti) + (t2 - ti)i/2cJ = ^(|t2 - ti|).
The proposition is proved.
□
Proposition 3, Proposition 4 and Arzela-Ascoli theorem (see, e.g., [l0, p. l02]) imply the validity of the following theorem.
Theorem 1. For each e € [0, l] the set of trajectories X£(to,xo) of the system (2.1) is a precompact subset of the space C ([to,9]; Rn).
Note that in general, the set of trajectories X£(to,xo) and X,*(to,xo) are not closed subsets of the space C ([to, 9]; Rn) (see, [3, 6]). Denote
Bra(a*) = {x € Rn : ||x|| < a*} , D„(a*) = {(t, x) € [to, 9] x Rn : x € B„(a*)} ,
where a* is defined by equality (2.6).
Here and henceforth we will have in mind the cylinder Dn(a*) as the set D in Condition 2.B.
■i
3. Properties of the set of trajectories
Denote
P* = 7i(0 - to) + 272c*(0 - to)1/2, (3.1)
g* = Yac*(0 - to)1/2 ■ exp(P*), (3.2)
Be(1) = (x(-) € C ([to,0]; ) : ||x(-)||e < 1} , (3.3)
where c* is defined in (2.2).
Proposition 5. For every e € [0,1] the inequality
1
hc{X£{to,xo),Xo{to,xo)) <g* - p-holds.
Proof. Let us choose an arbitrary x(-) € X£(t0,x0) generated by the control function u(-) € U. Define new control function u0(-) : [t0, 0] ^ Rm, setting
u0(t) = -7^=u(t), t € [t0, 9]. (3.4)
\/1 + £
The equality (3.4) yields that u0(-) € U0. Now, from (2.2), (3.4) and Cauchy-Schwarz inequality it follows that
NO - M-)\\i = [ (l - Il«(r)|| dr < c*(0 - io)1/2 (l - —=) . (3.5)
Let x0(-) : [t0,0] ^ be the trajectory of the system (2.1) generated by the control function u0(-) € Uo. Then x0(-) € X0(t0,x0). From Condition 2.B, (2.1) and (3.5) it follows that
||x(t) - X0(i)|W [71 + 72(IK<r)|| + |Mr)||)] ||x(r) - X0(t)| dr
Ji0 f ^ (3.6)
+73C*(0-io)1/2 ( 1-
VT+1
for every t € [t0, 0].
Taking into consideration the inequality (2.2), Gronwall-Bellman inequality and Cauchy-Schwarz inequality, from (3.1), (3.2) and (3.6) we obtain
IMt) - ro(t)ll < 1:M6 - to)1/2 (l - ' exP [/ hi + T2(Mr)|| + IK(r)ll)] dr]
< - (l - ■ exp(A) = (l -
for every t € [t0, 0], and hence
\\x(-) - x0(-)\\c < g* (1 - 1
Since x(-) € X£(t0, x0) is an arbitrarily chosen trajectory, x0(-) € X0(t0, x0), then the last inequality implies that
X£(to,xo) C X0(t0,x0) +g*(l- y==) Bc( 1) (3-7)
where BC(l) is defined by (3.3). From inclusion Xo(to,xo) C X£(to,xo) and (3.7) we obtain the proof of the proposition. □
From Proposition 5 it follows the validity of the following corollaries.
Corollary 1. hC(X£(to, xo), Xo(to, xo)) ^ 0 as e ^ 0+.
Corollary 2. For every e € [0, l] and t € [to, 9] the inequality
l
hn(X£(t]t0, x0), X0(t]t0, x0)) <g*y 1 - j-
is verified where the sets X£(t; to,xo), e € [0, l], are defined by (2.3). Denote
r* = 27sc* ■ exp(P*) (3.8)
where P* is defined by (3.l).
Proposition 6. For every e € [0, l] and a > 0 the inequality
r
hc(X£(t0,x0), X°(t0,x0))<-
a
is satisfied where r* is defined by (3.8).
Proof. Let us choose an arbitrary e € [0,1] and y(-) € X£(t0,x0) generated by the control function v(-) € U£. Define new control function v*(■) : [to,0] ^ Rm, setting
v(t) if ||v(t)|| < a,
V*® = < V(t) -a if 11^)11 > a. (3"9)
||v(t)||
Let
A* = {t € [to,9] : ||v(t)|| > a}.
Then from (2.2) we have
a2MA*) < / ||v(T)||2 dT < ||v(T)||2 dT < c*, JAt Jt0
'At J to
and hence
c2
KA*) < ^ (3-10)
a2
where ^(A*) stands for the Lebesgue measure of the set A*.
Since v(-) € U£ and ||v(t)|| > a for every t € A*, then (3.9) implies that
[ (B(t)v*(t), v*(t)) dT J to
= [ (B(T)v(T),v(r))dT+ [ (B(T)v(T),v(r))-TrJ-r^dT
J[to,e]\At J At ||v(T )|2
< / (B(t)v(t),v(t)) dT + / (B(t)v(t),v(t)) dT
./[to,0]\At At
= f (B(t)v(t),v(t)) dT < l + e. t0
(3.ll)
Now, (3.9) and (3.11) yield that v*(-) € U£a. Let y*(-) : [t0,0] ^ Rn be the trajectory of the system (2.1) generated by the control function v*(-) € U£\ Then y*(-) € Xa(t0,x0). Now the condition 2.B, inclusions v(-) € U£, v*(-) € U£\ (2.2), (3.9) and (3.10) imply that
|y(t) - y*(t)||< / [71 + 72(|v(r)| + ||v*(r)|)] ||y(r) - y*(r)| dr + 73/ |v(r) - v*(r)| dr Jt0 J At
< T [Y1 + Y2(|v(T)| + |v,(r)|)] ||y(T) - y*(r)| dT + 73 • |^(A,)]1/2[|v(.)|2 + N^h]
J to
< r [71 + 72(INt)|| + IK(r)ll)] \\y(r) - y,(r)\\ dr + ^
J to a
for every t € [t0,0]. The last inequality, the inclusions v(-) € U£, v*(-) € U£*, Bellman-Gronwall inequality and (3.8) yield
llv(i) - V.(i)|| < ^ -exp ( f [71 +72(IKr)|| + |Mr)||)] dr) < ^ -exp(^) = £ a V./io J a a
for every t € [t0, 0], and hence
r
M-)-y*{-)\\c<^. (3.12)
Since y(-) € X£(t0,x0) is an arbitrarily chosen trajectory, y*(-) € X£*(t0,£0), then from (3.12) it follows that
r
Xe(to,xQ) C Xf(to,xQ) + - • 5C(1), (3.13)
a
where the set Be(1) is defined by (3.3). Taking into consideration that Xa(t0,x0) C X£(t0,x0), we obtain from (3.13) the proof of the proposition. □
Corollary 3. he(X£(t0,£0),X£*(t0,x0)) ^ 0 as a ^ uniformly with respect to the e € [0,1].
From Propositions 5 and 6 it follows the validity of the following theorem. Theorem 2. For every e € [0,1] and a > 0 the inequality
1 r*
hc(XQ(to,xo),X?(to,xo)) < g* 1 - -== + -
V \/1 + ez a
is satisfied where g* and r* are defined by (3.2) and (3.8) respectively. Corollary 4. For every e € [0,1], a > 0 and t € [t0,0] the inequality
hn(Xo(t-to,xo),X«(t-to,xo)) < g* (I--J==] + -
V \/1 + e) a
is satisfied where X£*(t; t0,x0) is defined by (2.4).
4. Robustness of the trajectories
Let us discuss the robustness of the trajectories with respect to the fast consumption of the remaining control resource.
Proposition 7. Let v > 0 be a given number, Q* C [a, b] be Lebesgue measurable set, z() € X0(t0,x0) be a trajectory of the system (2.1) generated by the control function w(-) € U0,
i (B(t)w(t),w(t))dT = a* < 1,
J to
re I to
the control function w*(-) € L2([io, ; Rm) be such that
t-9
[ (B(t)w*(t),w*(t))dr = 1, w*(t) = w(t), t € [io,0] \ Q*, J to
and z*() be the trajectory of the system (2.1) generated by the control function w*(). If
2
, (4.1)
MQ*) <
then
2c*Y3 exp(^*)
||z(-) - z*(-)||c < v where c* is defined by (2.2), ¿5* is defined by (3.1).
Proof. Let us underline that the equality
ro
[ (B(t) w*(t),w*(t))dr = 1 t0
implies that w*() € U* and hence z*(-) € X*(to,xo). From Condition 2.B, inclusions w(-) € Uo, w*(-) € U*, (2.2) and definition of the control function w*(-) it follows that
|z(t) - z*(i)|W [71 + 72(||w(t)|| + ||w*(t)||)] ||z(t) - z*(t)|| dT + 73/ ||w(t) - w*(t)|| dT
-'to JQ»
< T [Y1 + 72(|w(T)|| + ||w*(T)|)] ||z(T) - Z*(T)| dT + Y3 ■ [^(Q*)]1/2 [||w(.)||2 + |w*(.)|2] J to
< i [yi + 72(|w(T)| + ||w*(T)|)] ||z(T) - z*(T)| dT + 273c* ■ [^(Q*)]1/2 t0
for every t € [to,0]. The last inequality, the inclusions w(-) € U0, w*(-) € U* Bellman-Gronwall inequality, (3.1) and (4.1) imply
||z(t) - z*(t)| < 273c*HQ*)]11/2 ■ exp ^jf [71 + 72(|w(T)| + ||w*(t)|)] dT
< 273c*HQ*)]1/2 ■ exp(^*) < v
for every t € [to,0], and consequently ||z() - z*(-)||e < v.
The proof is completed. □
Theorem 3. The equality
he (Xo(to,xo),X0* (to, xo)) = 0
holds.
Proof. Let v > 0 be an arbitrary fixed number and let us choose an arbitrary trajectory x( ■ ) € X0(t0,x0) of the system (2.1) generated by the control function u( ■ ) € U0. Assume that
f (B(t)U(T),U(T))dT = c* < 1. t0
Let Q* C [t0,0] be such that
v 12
KQ*) < 7)-T^TT ,
L 2c* Y3 exp(^*)J where c* is defined by (2.2), is defined by (3.1) and let
(4.2)
/ (B(T)U(T),U(T)) dT = ci, (4.3)
B* = / B(T) dT. (4.4)
Jq*
lQ*
It is obvious that B* is positive definite m x m matrix and c1 < c* < 1. Define new control function v0( ■ ) : [i0,0] — setting
-0(t)={if t € r]\Q- <4,»
where u0 € Rm is such that
(B*u0,u0) = 1 - ci. (4.6)
From (4.3), (4.4), (4.5) and (4.6) it follows that re
f (B(t)V0(T),V0(T) >) dT = / (B(T)U(T),U(T)) dT + f (B(T)U0,U0) dT
./to J [to,0]\Q* ./Q*
= ci + (( J B(T) d^ju0,U0^ = ci + (B*u0,u0) = ci + 1 - ci = 1,
and hence v0( ■ ) € U* Let x0( ■ ) be the trajectory of the system (2.1) generated by the control function v0( ■ ) € U(*. Then x0( ■ ) € X*(t0,x0) and from (4.2) and Proposition 7 we obtain
||x(■ ) - X0(■ )||e < v.
Since x( ■ ) € X0(t0,x0) is an arbitrarily chosen trajectory, x0( ■ ) € X*(t0,x0), then the last inequality implies that
X0(t0,x0) C X0*(t0,x0) + v-Bc(1) (4.7)
where Be(1) is defined by (3.3). Taking into consideration that X*(t0,x0) C X0(t0,x0), from (4.7) we obtain that
he(X*(t0,X0),X0(t0,X0)) < v. (4.8)
Since v > 0 is an arbitrarily fixed number, then (4.8) yields the proof of the theorem. □
From Theorem 3 we obtain the validity of the following corollaries. Corollary 5. The equality
cl(X0(t0,X0)) =cl (X0*(t0,X0))
is verified where cl denotes the closure of a set.
The Corollary 5 means that every trajectory x(-) € Xo(to,xo) of the system (2.1) can be approximated by the trajectory which is generated by full consumption of the control resource.
Corollary 6. For every t € [to,0] the equality
cl (Xo(t; to,Xo)) =cl (X*(t; to,xo))
is satisfied where Xo(t; to,xo) and X,*(t; to,xo) are defined by (2.3) and (2.5) respectively.
From Theorems 2 and 3 we obtain the validity of the following theorem.
Theorem 4. For every e € [0,1] and a > 0 the inequality
hc(XS(to,xo),X?(to,xo)) <g*( 1 - -J=) + -
V V1+e/ a
is satisfied where g* and r* are defined by (3.2) and (3.8) respectively.
5. Conclusion
The results asserting that a small perturbations in the quadratic integral constraints inspire a small deviation on the set of trajectories can be applied in mathematical modelling of the control systems where the total control resource is measured with small errors. According to the obtained results, it is possible to introduce a norm type geometric constraint along with quadratic type integral constraint where upper bound of the norm type geometric constraint is sufficiently large. Since the integrally constrained control functions are not geometrically constrained, this fact simplifies the structure of the set of control functions and allows to avoid geometrical unboundedness of the admissible control functions.
Robustness of the trajectories with respect to the fast and full consumption implies that it is reasonable to spend the control resource in economical mode, i.e. it is advisable to consume the control resource on the domains with sufficiently small Lebesgue measures in small portions. This yields that if you have a superfluous control resource and you want to get rid of this resource, then by spending all of the resource on the domain with sufficiently small Lebesgue measure, you will get a small deviation from the original system's trajectory.
REFERENCES
1. Beletskii V. V. Notes on the Motion of Celestial Bodies. Moscow: Nauka, 1972. 360 p.
2. Conti R. Problem di Controllo e di Controllo Ottimale. Torino: UTET, 1974. 239 p. (in Italian)
3. Filippov A. F. Differential Equations with Discontinuous Right-Hand Sides. Dordrecht: Kluwer, 1988. 304 p.
4. Gusev M.I., Zykov I.V. On the geometry of the reachable sets of control systems with isoperimetric constraints. Tr. Inst. Mat. Mekh. UrO RAN, 2018. Vol. 24, No. 1. P. 63-75. DOI: 10.21538/0134-4889-2018-24-1-63-75 (in Russia)
5. Guseinov K.G., Ozer O., Akyar E., Ushakov V.N. The approximation of reachable sets of control systems with integral constraint on controls. Nonlin. Dif. Equat. Appl. (NoDEA), 2007. Vol. 14, No. 1-2. P. 57-73. DOI: 10.1007/s00030-006-4036-6
6. Guseinov Kh. G., Nazlipinar A. S. On the continuity properties of the attainable sets of nonlinear control systems with integral constraint on controls. Abstr. Appl. Anal., 2008. Art ID 295817, 14 pp. DOI: 10.1155/2008/295817
7. Huseyin N., Huseyin A., Guseinov Kh. G. Approximations of the set of trajectories and integral funnel of the non-linear control systems with Lp norm constraints on the control functions. IMA J. Math. Control Inform., 2022. Vol. 39, No. 4. P. 1213-1231. DOI: 10.1093/imamci/dnac028
8. Huseyin A., Huseyin N., Guseinov Kh. G. Approximations of the images and integral funnels of the Lp balls under a Urysohn-type integral operator. Funktsionalnyi Analiz i ego Prilozheniya, 2022. Vol. 56, No. 4. P. 43-58. DOI: 10.4213/faa3974
9. Ibragimov G., Ferrara M., Kuchkarov A., Pansera B. A. Simple motion evasion differential game of many pursuers and evaders with integral constraints. Dynamic Games Appl., 2018. Vol. 8, No. 2. P. 352-378. DOI: 10.1007/s13235-017-0226-6
10. Kolmogorov A. N., Fomin S. V. Introductory Real Analysis. New York: Dover Publications, Inc., 1975. 403 p.
11. Kostousova E.K. On the polyhedral estimation of reachable sets in the "extended" space for discrete-time systems with uncertain matrices and integral constraints. Tr. Inst. Mat. Mekh. UrO RAN, 2020. Vol. 26, No. 1. P. 141-155. DOI: 10.21538/0134-4889-2020-26-1-141-155 (in Russia)
12. Krasovskii N. N. Theory of Control of Motion: Linear Systems. Moscow: Nauka, 1968. 475 p. (in Russian)
13. Motta M., Sartori C. Minimum time with bounded energy, minimum energy with bounded time. SIAM J. Contr. Optimiz., 2003. Vol. 42, No. 3. P. 789-809. DOI: 10.1137/S0363012902385284
14. Rousse R., Garoche P.-L., Henrion D. Parabolic set simulation for reachability analysis of linear timeinvariant systems with integral quadratic constraint. European J. Contr., 2021. Vol. 58. P. 152-167. DOI: 10.1016/j.ejcon.2020.08.002
15. Subbotin A. I., Ushakov V. N. Alternative for an encounter-evasion differential game with integral constraints on the players controls. J. Appl. Math. Mech., 1975. Vol. 39, No. 3. P. 387-396. DOI: 10.1016/0021-8928(75)90001-5
16. Subbotina N. N., Subbotin A. I. Alternative for the encounter-evasion differential game with constraints on the momenta of the players' controls. J. Appl. Math. Mech., 1975. Vol. 39, No. 3. P. 397-406. DOI: 10.1016/0021-8928(75)90002-7