A Distributed Cooperative Guidance Law with Prescribed-Time Consensus Performance


3.1. Design Strategy

Inspired by the mathematical expression of the total flight time in Equation (15), this paper proposes a two-step control strategy to achieve the distributed cooperative guidance law.

First, it is assumed that the aircraft design is based on Equation (7), with c i = 1 . Once the initial conditions are given, the total flight time can be calculated using Equation (15). Next, each time step is treated as the initial time, and the instantaneous state is considered as the initial state to recalculate the flight time using Equation (15). The equation is updated as follows:

t i ′ = σ i R i V i sin σ i ,   i = 1 , 2 , … , n .

Combining with Equation (16), the relative impact time error of the i-th aircraft is defined as follows:

ϵ i = ∑ j = 1 n a i j t i − t j ,   i = 1 , 2 , … , n .

where

t i = t i ′ + t = σ i R i V i sin σ i + t

Let ϵ = [ ϵ 1 , … , ϵ n ] T and t i = [ t 1 , … , t n ] T ; then, the compact vector form can be expressed as:

Second, the control parameters are adjusted to achieve consensus on the total flight time. From the previous discussion, the consensus on the total flight time allows for simultaneous arrival. Once the consensus error calculated using Equation (17) becomes zero, all aircraft will set their control parameters to c i = 1 and maintain that value.

3.2. Undirected Topology Cooperative Guidance Law

To accomplish the task of cooperative arrival, each aircraft must estimate its own remaining flight time while also obtaining the remaining flight times of other aircraft. By analyzing this information, they can continuously adjust their flight speeds to achieve cooperative guidance. During the flight, if aircraft can exchange information with one another, they can maximize the use of flight data, enhancing cooperative efficiency, which represents the optimal communication state. Based on these conditions, this subsection considers the ideal scenario of designing cooperative guidance under an undirected topology communication state.

Lemma 3.

A first-order multi-agent system with n agents is described by the following dynamic equation:

x ˙ i ( t ) = u i ( t ) ,   i = 1 , 2 , … , n .

where u i ( t ) ∈ R denotes the designed protocol, and x ˙ i ( t ) ∈ R represents the state of the i-th agent. If the graph Ω connecting the agents is connected and undirected, each agent follows the protocol defined as:

u i ( t ) = − ξ ˙ ( t ) 2 λ 2 ( L ) ( 1 − ξ ( t ) + δ ) + α ∑ j = 1 n a i j x i ( t ) − x j ( t )

where λ 2 ( L ) is defined by Lemma 1, 0 < δ ≪ 1 , α > 0 are constant parameters, and ξ ( t ) is a time-varying function with the following properties:

1.

ξ ( t ) is twice continuously differentiable on the interval 0 , ∞ .

2.

ξ ( t ) is monotonically increasing from ξ ( 0 ) = 0 to ξ ( t p r e ) = 1 , where t p r e < ∞ represents the prescribed time.

3.

ξ ′ ( 0 ) = ξ ′ ( t p r e ) = 0 .

4.

ξ ′ ( t ) = 0 and ξ ( t ) = 1 when t > t p r e .

Thus, the system will achieve consensus within the prescribed time, expressed as:

lim t → t pre | x i ( t ) − x j ( t ) | ≤ ( n − 1 ) 2 δ 1 + δ V 1 ( 0 ) | x i ( t ) − x j ( t ) | ≤ ( n − 1 ) 2 δ 1 + δ V 1 ( 0 ) , ∀ t > t pre lim t → + ∞ | x i ( t ) − x j ( t ) | = 0   i = 1 , 2 , … , n .

where V 1 ( 0 ) is the initial value of a non-negative function V 1 ( t ) = 1 2 x T ( t ) L x ( t ) .

For the cooperative arrival problem of multiple aircraft under an undirected topology, this paper makes the following assumptions:

Assumption 1.

The communication topology between the aircraft is undirected and connected.

Assumption 2.

During the flight, the changes are small-angle variations; thus,  
sin σ i ≈ σ i , cos σ i ≈ 1 − σ i 2 / 2 .

In this section, the control variable is selected as σ ˙ i , and a prescribed-time guidance law for stationary targets under an undirected graph is proposed, taking the following form:

σ ˙ i = − ( 1 + 2 ϵ i ξ ˙ ( t ) 2 λ 2 ( L ) ( 1 − ξ ( t ) + δ ) + α σ i 2 ) V i R i sin σ i ,   i = 1 , 2 , … , n .

where λ 2 ( L ) satisfies the definition in Lemma 1, δ , ξ ( t ) , and α meet the conditions outlined in Lemma 3, and ϵ i follows the definition in Equation (17).

Theorem 1.

If Assumption 1 and Assumption 2 are satisfied, the guidance law defined by Equation (23) guarantees that the actual flight time reaches prescribed-time consensus. Moreover, when  t > t p r e , due to  ϵ i ≈ 0 , it can be assumed that  σ ˙ i = V t p r e , i R t p r e , i sin σ t p r e , i = C i , where  C i  is a constant.

Proof 1.

Taking the derivative with respect to time yields:

t ˙ i = R ˙ i σ i sin σ i + σ ˙ i ( R i sin σ i − R i σ i cos σ i ) V i sin 2 σ i + 1

Combining Assumption 2 and substituting into Equations (2), (4), and (23) yields:

t ˙ i = R ˙ i σ i V i sin σ i − ( 1 + 2 ϵ i ξ ˙ ( t ) 2 λ 2 ( L ) ( 1 − ξ ( t ) + δ ) + α σ i 2 ) ( sin σ i − σ i cos σ i sin σ i ) + 1 = − σ i cos σ i sin σ i + σ i cos σ i sin σ i − 2 ϵ i ξ ˙ ( t ) 2 λ 2 ( L ) ( 1 − ξ ( t ) + δ ) + α σ i 2 ( 1 − σ i cos σ i sin σ i ) = − 2 ϵ i ξ ˙ ( t ) 2 λ 2 ( L ) ( 1 − ξ ( t ) + δ ) + α σ i 2 ( 1 − 1 + σ i 2 / 2 ) = − ∑ j = 1 n a i j ( t i − t j ) ξ ˙ ( t ) 2 λ 2 ( L ) ( 1 − ξ ( t ) + δ ) + α

Equation (25) satisfies the form of Equation (21), which confirms that the total flight time can reach consensus within the prescribed time.

Consider the Lyapunov function of the form:

V 1 ( t ) = 1 4 ∑ i = 1 n ∑ j = 1 n a i j t i − t j 2 = 1 2 t T L t

It is clear that V 1 ( t ) is positive and semi-definite, and V 1 ( t ) = 0 if and only if t 1 = t 2 = ⋯ = t n . This indicates that once V 1 ( t ) = 0 , consensus on the remaining flight time is achieved. Combining Property 3 and Lemma 3 allows for the derivation of the following:

lim t → t pre | t i − t j | ≤ ( n − 1 ) 2 δ ( 1 + δ ) V 1 ( 0 )   i , j = 1 , 2 , … , n .

where V 1 ( 0 ) is the initial value of V 1 ( t ) in Equation (26). Now, consider the case when t > t p r e . By applying Property 3 from Lemma 3 and Equation (25), it can be concluded that:

t ˙ i = − α ϵ i ,   t > t pre

Using Equations (19) and (28), the first-order time derivative of V 1 ( t ) can be expressed as:

V 1 . ( t ) = t T L t = − α ∑ i = 1 n ϵ i 2 ≤ 0 ,   t > t pre

This indicates that V 1 ( t ) monotonically decreases for t > t p r e . Combining with Equation (27), it follows that:

| t i − t f | ≤ ( n − 1 ) 2 δ ( 1 + δ ) V 1 ( 0 ) ,   t > t pre   i , j = 1 , 2 , … , n .

This means that when δ is chosen to be sufficiently small, | t i − t f | ensures that for all i , j ∈ V , the values will approach to a sufficiently small region within the prescribed time and remain in that neighborhood until interception, thereby ensuring the accuracy of the impact timing. □

3.3. Directed Topology Cooperative Guidance Law

In practical applications, the communication channels between aircraft are often unreliable due to the varying performance of communication devices. A more general scenario is examined to further lessen the communication burden, with the communication topology between the aircraft being directed. It is assumed that in addition to the n aircraft, there is one leader aircraft (indexed as 0). The communication topology among the n + 1 aircraft is characterized by a directed graph Ω ˜ ( ν ˜ , ε ˜ ) , where ε ˜ denotes the communication links and ν ˜ denotes the set of aircraft.

Lemma 4.

For a multi-agent system with a stationary leader (indexed as 0), let the dynamics of the leader be represented by  x ˙ 0 ( t ) = u 0 ( t ) = 0 , where  u 0 ( t ) ∈ R and  x 0 ( t ) ∈ R  represent the control input and state of the leader, respectively. It is assumed that the communication topology  Ω ˜  contains a directed spanning tree with the leader as the root. The rule established for each follower is:

u i ( t ) = − h max ξ ˙ ( t ) λ min ( U ) ( 1 − ξ ( t ) + δ ) + α ∑ j = 0 n a i j x i ( t ) − x j ( t )

For Equation (31), if δ and α satisfy the conditions of Lemma 3, h m a x and λ m i n ( U ) meet the definitions of Lemma 2; then:

lim t → t pre | x i ( t ) − x 0 ( t ) | ≤ 2 δ h min ( 1 + δ ) V 2 ( 0 ) | x i ( t ) − x 0 ( t ) | ≤ 2 δ h min ( 1 + δ ) V 2 ( 0 ) lim t → + ∞ | x i ( t ) − x 0 ( t ) | = 0 t i > t pre   i = 1 , 2 , … , n .  

where h min satisfies the definition of Lemma 2, and V 2 ( 0 ) denotes the initial value of the non-negative function V 2 ( t ) , defined as follows:

V 2 ( t ) = 1 2 e T ( t ) H e ( t )

For the functions V 2 ( t ) , H satisfies the definitions in Lemma 2 and has the following properties: e ( t ) = [ e 1 ( t ) , … , e n ( t ) ] T = [ x 1 ( t ) − x 0 ( t ) , … , x n ( t ) − x 0 ( t ) ] T .

Assumption 3.

Among the n + 1 aircraft, the communication topology is structured as a leader–follower graph  Ω ˜ ( ν ˜ , ε ˜ ) , where the leader aircraft acts as the root of a directed spanning tree. The communication network among the other aircraft, excluding the leader, is represented by  Ω ( ν , ε ) .

Remark 1.

It should be emphasized that the leader aircraft does not receive any information from other aircraft. Consequently, it adheres to the guidance law specified in Equation (14) to navigate towards the target.

The relative impact time error of the i-th aircraft (as a follower) in the leader–follower setup is defined as:

ϵ i = ∑ j = 0 n a i j t i − t j ,   i = 1 , 2 , … , n .

where t i = t i ′ + t = σ i R i V i sin σ i + t . Additionally, e ( t ) = [ e 1 ( t ) , … , e n ( t ) ] T and e i = t i − t 0 are defined.

Then, for a stationary target, the actual prescribed-time guidance law designed under the directed graph is expressed as:

σ ˙ i = − ( 1 + 2 ϵ i h max ξ ˙ ( t ) λ min ( U ) ( 1 − ξ ( t ) + δ ) + α σ i 2 ) V i R i sin σ i ,   i = 1 , 2 , … , n .

where δ , ξ ( t ) , and α satisfy the conditions outlined in Lemma 3, while h max and λ min ( U ) meet the conditions specified in Lemma 2.

Theorem 2.

If Assumption 2 and Assumption 3 hold, the guidance law defined by Equation (35) ensures that the actual flight time achieves prescribed-time consistency. Additionally,  | e i | , i = 1 , 2 , … , n .  is bounded for  t > t p r e .

Proof 2.

Under Remark 1, it is known that the leader’s acceleration variation satisfies a ˙ ≡ 0 . According to Equation (15), the leader’s total flight time can be expressed as:

t 0 = σ 0 R 0 V 0 sin σ 0

Then, the flight time t i is differentiated for the followers, and Equations (2), (4), and (35) are substituted into it to obtain:

t ˙ i = R ˙ i σ i V i sin σ i − ( 1 + 2 ϵ i h max ξ ˙ ( t ) λ min ( U ) ( 1 − ξ ( t ) + δ ) + α σ i 2 ) ( sin σ i − σ i cos σ i sin σ i ) + 1 = − σ i cos σ i sin σ i + σ i cos σ i sin σ i − 2 ϵ i h max ξ ˙ ( t ) λ min ( U ) ( 1 − ξ ( t ) + δ ) + α σ i 2 ( 1 − σ i cos σ i sin σ i ) = − 2 ϵ i h max ξ ˙ ( t ) λ min ( U ) ( 1 − ξ ( t ) + δ ) + α σ i 2 ( 1 − 1 + σ i 2 / 2 ) = − h max ξ ˙ ( t ) λ min ( U ) ( 1 − ξ ( t ) + δ ) + α ϵ i

Equation (37) satisfies the form of Equation (31), which indicates that the flight times t i can achieve consensus within the prescribed time.

Consider the Lyapunov function of the form:

V 2 ( t ) = 1 2 e T ( t ) H e ( t )

where H is a diagonal matrix with positive entries and satisfies the conditions of Lemma 2. Clearly, V 2 ( t ) is positive and definite, and for Equation (38), V 2 ( t ) = 0 if and only if all e i = 0 .

This indicates that once V 2 ( t ) = 0 , the total flight times of all aircraft can be expressed as t i = t 0 = σ 0 R 0 V 0 sin σ 0 , achieving consensus in flight times.

According to Lemma 4 and Equation (37), it can be concluded that when t → t p r e :

lim t → t pre | e i | ≤ 2 δ h min ( 1 + δ ) V 2 ( 0 ) ,   i = 1 , 2 , … , n .

where h min satisfies the definition in Lemma 2, and V 2 ( 0 ) is the initial value of V 2 ( t ) in Equation (38). Now consider the case when t > t p r e . According to Property 3 in Lemma 3 and Equation (37), the time derivative of e can be expressed as:

e ˙ = − α M e , t > t p r e

where M is defined by Lemma 2. Differentiating V 2 ( t ) and substituting Equation (40) yields:

V ˙ 2 ( t ) = − α e T H M e = − 1 2 α e T ( H M + M T H ) e = − 1 2 α e T U e ≤ − 1 2 α λ min ( U ) e T e ≤ 0

where from Lemma 2, U = H M + M T H > 0 . Therefore, V 2 ( t ) is monotonically decreasing for t > t p r e , and for Equation (39), it follows that:

| e i | ≤ 2 δ h min ( 1 + δ ) V 2 ( 0 ) , t > t pre , i , j = 1 , 2 , … , n .

This means that within the prescribed time, the errors of all aircraft can converge to a sufficiently small neighborhood around the origin and remain within that neighborhood before arrival occurs, thereby ensuring the accuracy of cooperative arrival timing. â–¡



Source link

Chao Ou www.mdpi.com