Next Article in Journal
On Some Examples of Trajectories in R7
Next Article in Special Issue
Dynamic Path Planning for the Differential Drive Mobile Robot Based on Online Metaheuristic Optimization
Previous Article in Journal
Prime Spectrum of the Ring of Adeles of a Number Field
Previous Article in Special Issue
Indirect Neural-Enhanced Integral Sliding Mode Control for Finite-Time Fault-Tolerant Attitude Tracking of Spacecraft
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Collision Avoidance Problem of Ellipsoid Motion

1
School of Business, Qingdao University, Qingdao 266071, China
2
School of Mathematics and Statistics, Qingdao University, Qingdao 266071, China
3
College of Science, Jiamusi University, Jiamusi 154007, China
*
Authors to whom correspondence should be addressed.
Mathematics 2022, 10(19), 3478; https://doi.org/10.3390/math10193478
Submission received: 22 August 2022 / Revised: 18 September 2022 / Accepted: 20 September 2022 / Published: 23 September 2022

Abstract

:
This paper studies the problem of target control and how a virtual ellipsoid can avoid the static obstacle. During the motion to the target set, the virtual ellipsoid can achieve a motion under collision avoidance by keeping the distance between the ellipsoid and obstacle. We present solutions to this problem in the class of closed-loop (feedback) controls based on Hamilton–Jacobi–Bellman (HJB) equation. Simulation results verify the validity and effectiveness of our algorithm.

1. Introduction

Control problems play an important role in modern scientific research. Data-driven control and fuzzy control problems for both linear and nonlinear system remain a popular task, and many theoretical researches have been carried out, such as Virtual Reference Feedback Tuning (VRFT) [1] and indirect adaptive iterative learning control (iAILC) scheme [2]. Formation control is one of the significant interest topics in the research of multi-agent system (MAS). Formation control can be realized by using various types of vehicles, such as unmanned aerial vehicles, mobile robots, underwater vehicles and so on. It has found a wide range of applications including military, aerospace, industry, etc.
The research of formation control mainly includes the following aspects: formation generation, formation maintenance, formation switching, formation avoidance and formation adaptive problems. The approaches to formation control are roughly categorized as leader-follower [3,4,5,6], behavioral [7,8], virtual structure [9,10] and artificial potential field [11,12].
Generally speaking, the collision avoidance problem in formation control includes two aspects. On the one hand, the members of MAS should avoid collision with each other in the process of motion. On the other hand, MAS should avoid collision with obstacles. A lot of research has been done in the field of obstacle avoidance algorithm. Many of them are based on artificial potential field [13,14]. Sometimes the control is realized with genetic algorithms (GA) [15], which is a powerful tool based on models of natural selection and evolution and allow an exhaustive search over large spaces. In [16] a visibility graph approach is used, while in [17] collision avoidance is achieved using A-Star algorithm.
In [18], Kurzhanski presents a theoretical framework of formation control approach for MAS based on the virtual ellipsoidal motions. This motion is achieved by means of a virtual ellipsoidal container inside which the members are assumed to be a sphere with safety volume. Kurzhanski considers the motions of MAS within the container so that the motions reach the target set together with the virtual container containing them. Therefore, his work considers about moving to the target set while avoiding collisions. In light of the above discussion, the motions of virtual ellipsoids can be fully portrayed by the linear differential equation satisfied with the coordinates of the center and the configuration matrix. Kurzhanski proposes the motion trajectory of virtual ellipsoid and MAS based on a suitable value function and dynamic programming method, thereby giving a solution scheme for MAS. Kurzhanski’s related work can be found in the literature [19,20,21,22,23,24,25,26].
In fact, we can regard Kurzhanski’s approach to formation control of MAS based on virtual ellipsoid as an improvement of virtual structure method. Traditional virtual structure views the formation of MAS as virtual rigid structure, in which each intelligent tracks a fixed virtual point motion on a rigid structure. This rigid structure limits on the application range of the method. While the formation control approach based on virtual ellipsoid motion not only inherits the advantage that it is easy to describe system behavior as a whole and can obtain higher-cell control accuracy from traditional virtual structures, but also fully compensates for the defects of rigid motion in traditional virtual structures.
Based on the work of Kurzhanski, this paper achieves the optimal target control problem based on virtual ellipsoid. This means that we should obtain optimal controls for the ellipsoid center and configuration matrix, which is significantly different from the traditional virtual structure method. Then, taking the actual situation of obstacle avoidance for formation control into account, we present the integration of virtual ellipsoid volume constraint, ellipsoid center constraint and configuration constraint into value function. Thereby we achieve optimal control to ellipsoid. Finally, we verify the effectiveness of our method by numerical simulations.

2. Statement of the Problem

We start this section with the definition of an ellipsoid in n .
Definition 1
([22]). A nondegenerate ellipsoid ε ( q , Q ) in n with center q n and configuration matrix Q n × n is the set
ε ( q , Q ) = { p n : p q , Q 1 ( p q ) 1 } ,
where Q is positive definite. Here , denotes inner product of two vector and Q 1 is the inverse matrix of Q .
On the time interval [ t 0 , θ ] , consider the virtual ellipsoidal motions of type
E c [ t ] = ε ( q ( t ) , Q ( t ) ) ,
with q t n and positive definite matrix function Q ( t ) n × n which are continuous at time t . E c [ t ] ( t [ t 0 , θ ] ) is referred as an ellipsoidal tube.
Now consider an ellipsoidal tube E c [ t ] ( t [ t 0 , θ ] ) with the following systems [23]:
q ˙ ( t ) = A q ( t ) q ( t ) + B q ( t ) u ( t , q ) , q ( t 0 ) = q 0 ,
Q ˙ ( t ) = T ( t ) Q ( t ) + Q ( t ) T ( t ) + B Q ( t ) U ( t , Q ) B Q ( t ) , Q ( t 0 ) = Q 0 ,
where u ( t , q ) m 1 controls the trajectory of q ( t ) and U ( t , Q ) m 2 × m 2 controls the configuration matrix Q ( t ) . The matrix parameters A q ( t ) n × n , B q ( t ) n × m 1 , B Q ( t ) n × m 2 , T ( t ) n × n are assumed to be continuously differentiable at time t , and T ( t ) is symmetric matrix. Here the prime of a matrix stands for the transposition of the matrix. The solvable conditions of system (1) and (2) have been discussed in [22] and Appendix B.
Next the used control constraint, geometric constraints, barrier constraint for the systems (1) and (2) are stated and the main problem is explained.

2.1. Control Constraint

Assume the admissible control sets to be u ( t , q ) m 1 and U ( t , Q ) m 2 × m 2 , where u ( t , q ) and U ( t , Q ) are given compact convex sets. For any control u ( t , q ) u ( t , q ) and U ( t , Q ) U ( t , Q ) , we have
t 0 θ ( u ( t , q ) , u ( t , q ) ) + [ U ( t , Q ) , U ( t , Q ) ] d t μ 2 ,
where the constant μ > 0 . Here [ U , U ] = tr ( U U ) denotes the inner product of two matrices.

2.2. State Constraints

Systems (1) and (2) are also considered under the additional state constraints
  [ Q ( t ) , Q ( t ) ] δ 2 ,
σ 2 [ Q ( t ) , Q ( t ) ] ,
where constants δ > σ > 0 . Note that inequality (4) defines a convex constraint, and inequality (5) defines a constraint complementary to a convex one. These inequalities restrict the possible size of an ellipsoid with configuration matrix Q ( t ) . Condition ( 4 ) implies that E c [ t ] is contained in a ball with a radius of δ / n and condition ( 5 ) implies that E c [ t ] contains a ball with a radius of σ / n .

2.3. Barrier Constraint

In this paper, obstacles are also regarded as virtual ellipsoids. Let the static ellipsoid B = ε ( q , Q ) be fixed obstacle on the path of the ellipsoid motion, where q n is the center and Q n × n is the configuration matrix which is positive definite. Let the eigenvalues of Q be λ 1 , λ 2 , , λ n , where λ i > 0 ,   i = 1 , 2 , , n . Denote
max { λ i | i = 1 , , n } = d 2 2 ,
where d / 2 is the radius of the smallest n-dimensional ball that includes B . While moving, E c [ t ] has to avoid obstacle B . We realize a motion under collision avoidance by keeping the distance between the E c [ t ] and obstacle B . Actually, we only limit the distance between q ( t ) and q to be greater than δ / n + d / 2 . Letting r = δ / n + d / 2 , we have
q ( t ) q , q ( t ) q r 2 .

2.4. Main Problem

With the ellipsoid E M = ε ( m , M ) being the target state, where m n is center, M n × n is configuration matrix and positive definite. Let the target set
M ω = { ( q , Q ) n × n × n | q m , q m + [ Q M , Q M ] ω 2 , ω > 0 }
be given in the form of a neighborhood of the ellipsoid E M .
This paper deals with the control problem of ellipsoid motion. The control objective in this paper is to drive the system to arrive at the target set, with the collision avoidance. By summarizing the above descriptions, we formulate the solution for systems ( 1 ) and ( 2 ) as follows.
Problem 1.
Given equation systems ( 1 ) and ( 2 ) on an interval [ t 0 , θ ] , let an initial state be { t 0 , q ( t 0 ) , Q ( t 0 ) } , therefore we get the initial ellipsoid E c [ t 0 ] = ε ( q ( t 0 ) , Q ( t 0 ) ) . Find feedback controls u t , q and U ( t , Q ) that transfer the ellipsoid E c [ t ] from the state E c [ t 0 ] into terminal set M ω under the constraints (3)–(6).
In order to find the optimal feedback controls u t , q and U ( t , Q ) of problem, we define the objective function
Ψ ( u ( ) , U ( ) ) = t 0 θ ( α 1 ( u ( t , q ) , u ( t , q + [ U ( t , Q ) , U ( t , Q ) ] ) + α 2 ( δ 2 [ Q ( t ) , Q ( t ) ] ) + α 3 ( [ Q ( t ) , Q ( t ) ] σ 2 ) β q ( t ) q , q ( t ) q ) d t + γ ( q ( θ ) m , q ( θ ) m + [ Q ( θ ) M , Q ( θ ) M ] ) ,
where constants α 1 , α 2 , α 3 , β , γ are referred to as weight coefficients.
In this paper, we use the distance between the ellipsoid and the obstacle to constrain the objective function, in which the movement of the ellipsoid is restricted by term q ( t ) q , q ( t ) q to realize collision avoidance.

3. Solutions Developed

In this section, we use dynamic programming method to solve the above-mentioned optimal control problem. The value function defined on the trajectories of systems (1) and (2) is
V E ( t , q , Q ) = min u , U { Ψ ( u ( ) , U ( ) ) | q ( t ) = q , Q ( t ) = Q } .
For fixed α 1 , α 2 , α 3 , β , γ , we obtain ([24])
W [ t ] = q ( t ) , Q ( t ) | V E ( t , q , Q ) 1 ,
where W [ t ] is the backward reach set relative to M ω under the constraints (3)–(6), i.e., the set of points { q ( t ) , Q ( t ) } for which there exist controls u ( t , q ) and U ( t , Q ) that steers E c [ t ] to E ω under the constraints (3)–(6). The value function for this problem is
V E ( t , q , Q ) = max α 1 , α 2 , α 3 , γ min β { V E ( t , q , Q ) | α 1 , α 2 , α 3 , β , γ }
over all { α 1 , α 2 , α 3 , β , γ 0 :   α 1 μ 2 + α 2 δ 2 + α 3 σ 2 + β r 2 + γ ω 2 = 1 } . The meanings of the parameters have been listed in the Table 1 below.
We can substitute the weight coefficients of this problem into value function to obtain the corresponding control. The obtained solutions will depend on μ , δ , σ , r , ω .
Solving this problem is equivalent to minimizing this value function
V E ( t , q , Q ) = min u , U t θ ( α 1 ( u ( τ , q ) , u ( τ , q ) + [ U ( τ , Q ) , U ( τ , Q ) ] ) + α 2 ( δ 2 [ Q ( τ ) , Q ( τ ) ] ) + α 3 ( [ Q ( τ ) , Q ( τ ) ] σ 2 ) β q ( τ ) q , q ( τ ) q ) d τ + γ ( q ( θ ) m , q ( θ ) m + [ Q ( θ ) M , Q ( θ ) M ] ) .
For the convenience of calculation, we rewrite the system ( 2 ) in vector form [25]. Let A = ( a i j ) n × n , introducing the notation
A ¯ = [ a 11 , a 12 , , a 1 n , a 21 , a 22 , , a 2 n , , a n 1 , a n 2 , , a n n ] .
Let B = ( b i j ) n × n , the Kronecker product of matrices is denoted as
A B = a 11 B a 1 n B a n 1 B a n n B n 2 × n 2 .
Using the identity A X B ¯ = ( A B ) X ¯ , then system ( 2 ) can be rewritten as
Q ¯ ˙ = ( T ( t ) I ) Q ¯ + ( I T ( t ) ) Q ¯ + ( B Q ( t ) B Q ( t ) ) U ¯ ( t , Q ) .
Denoting A ( t ) = T ( t ) I + I T ( t ) , B ( t ) = B Q ( t ) B Q ( t ) , we obtain the relation
Q ¯ ˙ = A ( t ) Q ¯ + B ( t ) U ¯ ( t , Q ) .
Then the value function can be replaced by the form
V E ( t , q , Q ) = min u , U ¯ t θ ( α 1 ( u ( τ , q ) , u ( τ , q ) + U ¯ ( τ , Q ) , U ¯ ( τ , Q ) ) + α 2 ( δ 2 Q ¯ ( τ ) , Q ¯ ( τ ) ) + α 3 ( Q ¯ ( τ ) , Q ¯ ( τ ) σ 2 ) β q ( τ ) q , q ( τ ) q ) d τ + γ ( q ( θ ) m , q ( θ ) m + Q ¯ ( θ ) M ¯ , Q ¯ ( θ ) M ¯ ) .
Then function ( 9 ) derives a solution to the HJB equation ([27])
V E t + min u , U { V E q , A q ( t ) q + B q ( t ) u ( t , q ) + V E Q ¯ , A ( t ) Q ¯ + B ( t ) U ¯ ( t , Q ) + α 1 ( u ( t , q ) , u ( t , q ) + U ¯ ( t , Q ) , U ¯ ( t , Q ) ) + α 2 ( δ 2 Q ¯ , Q ¯ ) + α 3 ( Q ¯ , Q ¯ σ 2 ) β q q , q q } = 0 ,
with the boundary condition
V E ( θ , q ( θ ) , Q ( θ ) ) = γ ( q ( θ ) m , q ( θ ) m + Q ¯ ( θ ) M ¯ , Q ¯ ( θ ) M ¯ ) .
The solution to Equation ( 10 ) is given as Equation ( 9 ) . Denoting
H ( t , u ( t ) , U ¯ ( t ) ) = V E q , A q ( t ) q + B q ( t ) u ( t , q ) + V E Q ¯ , A ( t ) Q ¯ + B ( t ) U ¯ ( t , Q ) + α 1 ( u ( t , q ) , u ( t , q ) + U ¯ ( t , Q ) , U ¯ ( t , Q ) ) + α 2 ( δ 2 Q ¯ , Q ¯ ) + α 3 ( Q ¯ , Q ¯ σ 2 ) β q q , q q } = 0 .
Then Equation ( 10 ) is equivalent to
V E t + H ( t , u , U ¯ ) = 0 , H ( t , u , U ¯ ) u = 0 , H ( t , u , U ¯ ) U ¯ = 0 .
First, we have the partial derivative of H ( t , u , U ¯ ) with respect to u and U
u ( t , q ) = 1 2 α 1 B q ( t ) V E q ,
U ¯ ( t , Q ) = 1 2 α 1 B ( t ) V E Q ¯ .
Substituting ( 12 ) , ( 13 ) into V E t + H ( t , u , U ¯ ) = 0 , Equation ( 10 ) can be rewritten as
V E t + V E q , A q ( t ) q 1 4 α 1 V E q , B q ( t ) B q ( t ) V E q + V E Q ¯ , A ( t ) Q ¯ 1 4 α 1 V E Q ¯ , B ( t ) B ( t ) V E Q ¯ + α 2 ( δ 2 Q ¯ , Q ¯ ) + α 3 ( Q ¯ , Q ¯ σ 2 ) β q q , q q = 0 .
The resulting expression is a quadratic form in the state coordinates and the spatial derivatives of the value function. The latter permits one to seek the value function in quadratic form as well,
V E ( t , q , Q ) = q , p ( t ) q + [ q , k ( t ) ] + [ Q , P ( t ) Q ] + [ Q , K ( t ) ] + s ( t ) ,  
the above equation is equivalent to (16)
V E ( t , q , Q ) = q , p ( t ) q + q , k ( t ) + Q ¯ , P ( t ) Q ¯ + Q ¯ , K ( t ) + s ( t ) ,
where P ( t ) and p ( t ) are symmetric positive definite, P ( t ) = P ( t ) I n × n , K ( t ) = K ¯ ( t ) . We have the partial derivative of function V E ( t , q , Q ) with respect to t , q and Q ¯
V E t = q , p ˙ ( t ) q + q , k ˙ ( t ) + Q ¯ , P ˙ ( t ) Q ¯ + Q ¯ , K ˙ ( t ) + s ˙ ( t ) , V E q = 2 p ( t ) q + k ( t ) , V E Q ¯ = 2 P ( t ) Q ¯ + K ( t ) .
Equations ( 12 ) and ( 13 ) can be rewritten as
u ( t , q ) = 1 2 α 1 B q ( t ) ( 2 p ( t ) q + k ( t ) ) ,
U ¯ ( t , Q ) = 1 2 α 1 B ( t ) ( 2 P ( t ) Q ¯ + K ( t ) ) .
By substituting Equation ( 15 ) into Equations ( 11 ) and ( 14 ) , we thereby obtain equations for the parameters of the form ( 15 ) , such as
p ˙ ( t ) + 2 p ( t ) A q ( t ) 1 α 1 p ( t ) B q ( t ) B q ( t ) p ( t ) β I n × n = 0 , p ( θ ) = γ I n × n , k ˙ ( t ) + A q ( t ) k ( t ) 1 α 1 p ( t ) B q ( t ) B q ( t ) k ( t ) + 2 β q = 0 , k ( θ ) = 2 γ m , P ˙ ( t ) + 2 P ( t ) A ( t ) 1 α 1 P ( t ) B ( t ) B ( t ) P ( t ) ( α 2 α 3 ) I n 2 × n 2 = 0 , P ( θ ) = γ I n 2 × n 2 , K ˙ ( t ) + A ( t ) K ( t ) 1 α 1 P ( t ) B ( t ) B ( t ) K ( t ) = 0 , K ( θ ) = 2 γ M ¯ , s ˙ ( t ) 1 4 α 1 k ( t ) , k ( t ) 1 4 α 1 K ( t ) , K ( t ) + α 2 δ 2 α 3 σ 2 β q , q = 0 , s ( θ ) = γ ( m , m + M ¯ , M ¯ ) .
Under the condition of continuous differentiability of the matrix parameters of the systems ( 1 ) and ( 2 ) , the functions ( 10 ) and ( 11 ) have a unique classical solution. We have thereby proved the following assertion.
Theorem 1.
The value function ( 15 ) in which the parameters are determined by systems ( 18 ) specifies a solution to Problem 1. In this case, the optimal controls u ( t , q ) and U t , Q ¯ are given by ( 16 ) and ( 17 ) .
Since the above formulas involve nonlinear terms, it is difficult for us to obtain analytical expressions even for this relatively simple model, so we focus on numerical solutions. In a similar way as [26], forward Explicit Euler method (in reverse time) is introduced to perform numerical discretization, and the ellipsoid trajectory tube is obtained by MATLAB software. The corresponding implemented algorithm has been attached in Appendix A.

4. Numerical Simulation

In order to verify the control method of collision avoidance proposed in this paper, numerical simulation results are presented in this section. The simulation includes two cases: the presence of obstacle constraint and the absence of obstacle constraint.
Consider the following parameters of the dynamics ( 1 ) and ( 2 ) of the center and the configuration matrix for a solution that uses Equation ( 8 ) .
q 0 = [ 0 , 0 ] , Q 0 = 0.005 0 0 0.005 , t 0 = 0 , θ = 1 , A q = 9 t 6 20 6 20 9 t , B q = t 3 t 3 t 3 t 3 , T = t 2 2 t , B Q = cos π t sin π t sin π t cos π t .
The center and the configuration matrix of the target ellipse are given as
m = [ 1 , 1 ] , M = 0.005 0 0 0.005 .
The center and the configuration matrix of the obstacle are described as
q = [ 0.5 , 0.5 ] , Q = 0.005 0 0 0.005 .
Now we choose the parameters of Equations (3)–(7) as follows
μ 2 = ω 2 = r 2 = 0.2 ,   δ 2 = 0.28 ,   σ 2 = 0.1 .
Omitting the presence of obstacle constraint, we consider Equation ( 8 ) with the parameters specified as
α 1 = 0.019 ,   α 2 = 0.075 ,   α 3 = 0.035 ,   β = 0 ,   γ = 0.057 .
Considering the presence of obstacle constraint, the parameters in Equation ( 8 ) are defined as follows
α 1 = 0.012 ,   α 2 = 0.057 ,   α 3 = 0.029 ,   β = 0.048 ,   γ = 0.04 .
The elliptical tubes of trajectories without regard of the obstacle constraints are shown in Figure 1, and the elliptical tubes of trajectories with regard of the obstacle constraints are shown in Figure 2, where the vertical trajectories tube is the obstacle, and the curved trajectories tube is E c [ t ] . It can be seen from the simulation results in Figure 1 and Figure 2 that when obstacle constraint is not considered, E c [ t ] cannot avoid the obstacle. When considering obstacle constraint, E c [ t ] can avoid the obstacle.
The comparison diagram of the terminal ellipse and the target ellipse without obstacle constraint and with obstacle constraint are given. The simulation results are provided in Supplementary Materials and shown in Figure 3, Figure 4 and Figure 5. In Figure 3, the center of E c [ t ] reaches the target set without obstacle constraint. In Figure 4, the center of E c [ t ] reaches the target set with obstacle constraint. In Figure 5, the configuration matrix E c [ t ] reaches the target set with or without obstacle constraint.

5. Conclusions

This paper studies the collision avoidance problem in formation control and presents a solution to realize collision avoidance based on dynamic programming. By introducing a barrier constraint into the value function, we obtain the optimal controls on the virtual ellipsoid to pass through the obstacle. For the three-dimensional case of a virtual ellipsoid, the method using this article can obtain similar results.
Our work focuses on a situation with static obstacles. This framework is flexible and extensible, and the results of collision avoidance are stable. Meanwhile, several limitations remain to be studied in the future. For example, the current algorithm has not yet achieved parallel computing and still needs to be improved. Additionally, when the obstacle is movable or when actively interfering movement towards the virtual ellipsoid is considered, the target control problem will be converted to a kind of complex game problem which will be studied in the future.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/math10193478/s1, Txt S1: Value of center q on time interval [0, 1] step 0.01.txt, Txt S2: Value of configuration matrix Q on time interval [0, 1] step 0.01.txt, Txt S3: Value of control u on time interval [0, 1] step 0.01.txt, Txt S4: Value of control bar{U} on time interval [0, 1] step 0.01.txt.

Author Contributions

Conceptualization, S.G.; Methodology, S.G., L.J., Z.D. (Zhaopeng Dai) and Y.Y.; Software, A.S. and Z.D. (Zhiqing Dang); Formal analysis, Z.Y., J.G. and Y.S.; Validation, S.G. and H.G.; Writing—original draft preparation, S.G.; Writing—review and editing, S.G., L.J., Z.D. (Zhaopeng Dai) and Y.Y.; Supervision, H.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Natural Science Foundation of China (No. 72171126, 11872220).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Algorithm A1. Numerical Solutions for Theorem 1.
Require: The total time interval t 0 , θ , the number of time step n t , the initial center q 0 and configuration matrix Q 0 , the target center q and configuration matrix Q , and the obstacle center m and configuration matrix M for the ellipsoid. The coefficients of dynamic formulae A q , B q , T , B Q , μ , δ , σ , r , ω depending on the system, and Practitioner designed coefficients α 1 , α 2 , α 3 , β , γ .
 
Ensure: The optimal controls u t , q for the center and U t , Q for the configuration matrix of the ellipsoid ε q , Q which allow ε q , Q to move to the target ε q , Q without collision with ε m , M , and the center q t , configuration matrix Q t of ε q , Q . The above solution was given at the following time node: t 0 , t 0 + Δ t ,   t 0 + 2 Δ t , , t 0 + n t Δ t = θ , where Δ t = θ t 0 n t .
 
%Solving p t , k t , P t , K t , s t in Equation (18) using forward Explicit Euler Method in reverse time.
 
 1: p θ = γ I n × n , k θ = 2 γ m , P θ = γ I n 2 × n 2 ,
   K θ = 2 γ M ¯ , s θ = γ m , m + M ¯ , M ¯ ;
   A t = T t I n × n , B t = B Q B Q ;
 
 2: for t = θ ; t t 0 ; t = Δ t {
    p t Δ t = p t 1 Δ t 2 A q t p t + 1 α p t B q t B q t p t + β I n × n ;
    P t Δ t = P t 1 Δ t 2 A t P t + 1 α P t B t B t P t + α 2 α 3 I n 2 × n 2 ;
    k t Δ t = k t 1 Δ t A q t k t + 1 α p t B q t B q t k t 2 β m ;
    K t Δ t = K t 1 Δ t A t K t + 1 α P t B t B t P t + α 2 α 3 I n 2 × n 2 ;
 }
 
%Solving u t , q , U t , q using equation (16) and (17), and solving q t , Q t using Equations (1) and (2) in the same loop.
 
 3: u t 0 = 0 , U t = Z e r o   M a t r i x   w i t h   s i z e   n × n ,   q t 0 = q 0 , Q t 0 = Q 0 ;
 
 4: for t = t 0 ; t θ ; t + = Δ t {
    u t + Δ t = 1 2 α B q t 2 p t q t + k t ;
    U t + Δ t = 1 2 α B t 2 P t Q t + K t ;
    q t + Δ t = q t + A q t q t + B q t u t Δ t ;
    Q t + Δ t = Q t + A t Q t + B t U t Δ t ;
 }

Appendix B. Solvability Analysis

Discussion about the solvability of system provided by Equations (1) and (2) are provided in this appendix. In the view of formation control, Equations (1) and (2) form a coupled system that should not be analyzed separately. Thus, we introduce the concept of Ellipsoidal Dynamics with notations in reference [22] that make our demonstration theoretical and concise. Several definitions below have already been given in [22] and without special notation we omit the name of this book and only mark the corresponding content with page number.
First, we rewrite Equations (1) and (2) into an Attainable Domain (Definition 1.2.1 Page 9 in [22]) Equation:
X t = ε q 0 , Q 0 + t 0 t ε 0 , P t d t , t t 0 , θ ,
where X t = { x t | c o n t r o l   u t , q , U t , q ,   s . t .   x t ε q t , Q t   w i t h   u , U } , ε 0 , P t denotes the ellipsoid constructed by an eligible control, and by Equation (3) it forms an ellipsoid (actually it is a sphere with isotropy) with the center as origin. The + and integral sign in Equation (A1) denotes Ellipsoid Sum and Integral, defined on Pages (128) and (161) in [22], respectively.
Next, by carrying out the similar procedure from Section 3.1 (on Page 178) to 3.4 (on Page 194) in [22], Corresponding Theorem 3.4.1 (on Page 195) in [22] holds for our system, which implies a condition for solvability. We transform this conclusion with our notation as below.
Theorem A1.
ε q t , Q _ t W t ε + q t , Q + t must hold for every t , where ε q t , Q t and ε + q t , Q + t denotes the ellipsoid evolution (Definition 3.3.1 on Page 191 and Definition 3.3.2 on Page 193) carried out by internal or external calculation, respectively, and W t denotes the Solvable Domain (Definition 1.4.2 on Page 19 in [22]), which is
W t = ε q , Q t θ ε 0 , P t d t , t t 0 , θ ,
where the denotes Ellipsoid difference as defined on Page (128) in [22].
From Theorem A1, we could check whether Equations (1) and (2) is solvable by checking whether
ε q 0 , Q 0 ε q , Q t θ ε 0 , P t d t
holds for the given dynamic system with corresponding coefficients.

References

  1. Roman, R.-C.; Precup, R.-E.; Petriu, E.M. Hybrid Data-Driven Fuzzy Active Disturbance Rejection Control for Tower Crane Systems. Eur. J. Control 2021, 58, 373–387. [Google Scholar] [CrossRef]
  2. Chi, R.; Li, H.; Shen, D.; Hou, Z.; Huang, B. Enhanced P-type Control: Indirect Adaptive Learning from Set-point Updates. IEEE Trans. Automat. Control 2022, 665–669. [Google Scholar] [CrossRef]
  3. Gustavi, T.; Hu, X.M. Observer-based leader-following formation control using onboard sensor information. IEEE Trans. Robot. 2008, 24, 1457–1462. [Google Scholar] [CrossRef]
  4. Mariottini, G.L.; Morbidi, F.; Prattichizzo, D.; Valk, N.V.; Michael, N.; Pappas, G.; Daniilidis, K. Vision-based localization for leader-follower formation control. IEEE Trans. Robot. 2009, 25, 1431–1438. [Google Scholar] [CrossRef]
  5. Chen, X.; Yan, P.; Serrani, A. On input-to-state stability-based design for leader/follower formation control with measurement delays. Int. J. Robust Nonlin. 2013, 23, 1433–1455. [Google Scholar] [CrossRef]
  6. Panagou, D.; Kumar, V. Cooperative visibility maintenance for leader-follower formations in obstacle environments. IEEE Trans. Robot. 2014, 30, 831–844. [Google Scholar] [CrossRef]
  7. Kownacki, C. Multi-UAV flight using virtual structure combined with behavioral approach. Acta Mech. Autom. 2016, 10, 92–99. [Google Scholar] [CrossRef]
  8. Antonelli, G.; Arrichiello, F.; Chiaverini, S. Experiments of formation control with multirobot systems using the null-space-based behavioral control. IEEE. Trans. Control Syst. Technol. 2009, 17, 1173–1182. [Google Scholar] [CrossRef]
  9. Ren, W.; Beard, R.W. Decentralized scheme for spacecraft formation flying via the virtual structure approach. J. Guid. Control. Dynam. 2004, 27, 73–82. [Google Scholar] [CrossRef]
  10. Sadowska, A.; den Broek, T.V.; Huijberts, H.; van deWouw, N.; Kostić, D.; Nijmeijer, H. A virtual structure approach to formation control of unicycle mobile robots using mutual coupling. Int. J. Control 2011, 84, 1886–1902. [Google Scholar] [CrossRef]
  11. Gazi, V. Swarm aggregations using artificial potentials and sliding-mode control. IEEE Trans. Robot. 2005, 21, 1208–1214. [Google Scholar] [CrossRef]
  12. Mabrouk, M.; McInnes, C. Solving the potential field local minimum problem using internal agent states. Robot. Auton. Syst. 2008, 56, 1050–1060. [Google Scholar] [CrossRef]
  13. Zhang, J.L.; Yan, J.G.; Zhang, P. Fixed-wing UAV formation control design with collision avoidance based on an improved artificial potential field. IEEE. Access 2018, 6, 78342–78351. [Google Scholar] [CrossRef]
  14. Xie, L.; Liu, S. Dynamic obstacle-avoiding motion planning for manipulator based on improved artificial potential filed. IET Control Theory A 2018, 35, 1239–1249. [Google Scholar]
  15. Jiang, A.; Yao, X.; Zhou, J. Research on path planning of real-time obstacle avoidance of mechanical arm based on genetic algorithm. J. Eng. 2018, 16, 1579–1586. [Google Scholar] [CrossRef]
  16. Li, S.L.; He, J.H.; Ao, H.Y.; Liu, Y.B. Intelligent visibility graph algorithm of Mars aircraft based on pigeon-inspired optimization. Flight. Dyn. 2020, 38, 90–94. [Google Scholar]
  17. Chen, G.R.; Guo, S.; Wang, J.Z.; Qu, H.B.; Chen, Y.Q.; Hou, B.W. Convex optimization and A-star algorithm combined path planning and obstacle avoidance algorithm. J. Control Decis. 2020, 35, 2907–2914. [Google Scholar]
  18. Kurzhanski, A.B. Problem of collision avoidance for a team motion with obstacles. Proc. Steklov Inst. Math. 2016, 293, 120–136. [Google Scholar] [CrossRef]
  19. Kurzhanski, A.B.; Mitchell, I.M.; Araiya, P.V. Optimization techniques for state-constrained control and obstacle problems. J. Optim. Theory Appl. 2006, 128, 499–521. [Google Scholar] [CrossRef]
  20. Kurzhanski, A.B. The problem of measurement feedback control. J. Appl. Math. Mech. 2004, 68, 487–501. [Google Scholar] [CrossRef]
  21. Komarov, Y.; Kurzhanski, A.B. Minimax-maximin relations for the problem of vector-valued criteria optimization. Dokl. Math. 2020, 101, 259–261. [Google Scholar] [CrossRef]
  22. Kurzhanski, A.B.; Valyi, I. Ellipsoidal Calculus for Estimation and Control; Springer: Berlin/Heidelberg, Germany, 1997. [Google Scholar]
  23. Kurzhanski, A.B. On the problem of control for ellipsoidal motions. Proc. Steklov Inst. Math. 2012, 277, 160–169. [Google Scholar] [CrossRef]
  24. Kurzhanski, A.B. On a team control problem under obstacles. Proc. Steklov Inst. Math. 2015, 51, 128–142. [Google Scholar] [CrossRef]
  25. Kurzhanski, A.B.; Mesyats, A.I. Optimal control of ellipsoidal motions. Differ. Equ. 2012, 48, 1502–1509. [Google Scholar] [CrossRef]
  26. Kurzhanski, A.B.; Mesyats, A.I. Control of ellipsoidal trajectories: Theory and numerical results. Comput. Math. Math. Phys. 2014, 54, 418–428. [Google Scholar] [CrossRef]
  27. Dockner, E.J.; Jorgensen, S.; Long, N.V.; Sorger, G. Differential Games in Economics and Management Science, 1st ed.; Cambridge University Press: Cambridge, UK, 2000. [Google Scholar]
Figure 1. Elliptical tubes of trajectories for the E c [ t ] with no regard for the obstacle constraint.
Figure 1. Elliptical tubes of trajectories for the E c [ t ] with no regard for the obstacle constraint.
Mathematics 10 03478 g001
Figure 2. Elliptical tubes of trajectories for the E c [ t ] with regard for the obstacle constraint.
Figure 2. Elliptical tubes of trajectories for the E c [ t ] with regard for the obstacle constraint.
Mathematics 10 03478 g002
Figure 3. The comparison diagram of the center of the terminal ellipse and the target ellipse without obstacle constraint.
Figure 3. The comparison diagram of the center of the terminal ellipse and the target ellipse without obstacle constraint.
Mathematics 10 03478 g003
Figure 4. The comparison diagram of the center of the terminal ellipse and the target ellipse with obstacle constraint.
Figure 4. The comparison diagram of the center of the terminal ellipse and the target ellipse with obstacle constraint.
Mathematics 10 03478 g004
Figure 5. The comparison diagram of the configuration matrix of the terminal ellipse and the target ellipse.
Figure 5. The comparison diagram of the configuration matrix of the terminal ellipse and the target ellipse.
Mathematics 10 03478 g005
Table 1. Details of all the parameters.
Table 1. Details of all the parameters.
α 1 Constraining controls of the ellipsoid
α 2 Constraining the volume of the ellipsoid by constraining an upper bound on the norm of the ellipsoid configuration matrix
α 3 Constraining the volume of the ellipsoid by constraining a lower bound on the norm of the ellipsoid configuration matrix
β Constraining the distance between the center of ellipsoid and the obstacle
γ Constraining the distance between the ellipsoid and the target
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Guo, S.; Jing, L.; Dai, Z.; Yu, Y.; Dang, Z.; You, Z.; Su, A.; Gao, H.; Guan, J.; Song, Y. Collision Avoidance Problem of Ellipsoid Motion. Mathematics 2022, 10, 3478. https://doi.org/10.3390/math10193478

AMA Style

Guo S, Jing L, Dai Z, Yu Y, Dang Z, You Z, Su A, Gao H, Guan J, Song Y. Collision Avoidance Problem of Ellipsoid Motion. Mathematics. 2022; 10(19):3478. https://doi.org/10.3390/math10193478

Chicago/Turabian Style

Guo, Shujun, Lujing Jing, Zhaopeng Dai, Yang Yu, Zhiqing Dang, Zhihang You, Ang Su, Hongwei Gao, Jinqiu Guan, and Yujun Song. 2022. "Collision Avoidance Problem of Ellipsoid Motion" Mathematics 10, no. 19: 3478. https://doi.org/10.3390/math10193478

APA Style

Guo, S., Jing, L., Dai, Z., Yu, Y., Dang, Z., You, Z., Su, A., Gao, H., Guan, J., & Song, Y. (2022). Collision Avoidance Problem of Ellipsoid Motion. Mathematics, 10(19), 3478. https://doi.org/10.3390/math10193478

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