Cooperative Midcourse Guidance Law with Communication Delay

The initial conditions of the cooperative terminal guidance law, which are the terminal conditions of the cooperative midcourse guidance, have a greater impact on its cooperation and guidance precision, so it is worthy of investigating the cooperative midcourse guidance. In addition, the problem of communication delay between the network nodes is inevitable and has a greater impact on the cooperative guidance law. To solve the above problems, a novel distributed cooperative midcourse guidance (DCMG) law with communication delay is proposed by combining the cooperative term with a distributed consensus protocol including communication delay under the directed communication topology. Firstly, a DCMG law with communication delay is designed by combining the trajectory shaping guidance with the distributed protocol including communication delay under the directed communication topology; secondly, the consensus of the proposed DCMG law with communication delay under the directed graph is proved; finally, the effectiveness and superiority of the proposed DCMG law are verified by numerical simulations.


Introduction
As more and more high-value targets are equipped with antimissile systems, the traditional single missile's killing ability of attacking a single target has seriously deteriorated. Therefore, a multimissile cooperative attacking strategy is effective for increasing the killing possibility. For mediumand long-range missiles, the cooperative midcourse guidance phase needs to provide good initial conditions for the cooperative terminal guidance phase, ensuring that the cooperative terminal guidance law can complete the task of hitting the target at the same time within a very limited time. During the last few decades, cooperative guidance has attracted much more attention.
In the early studying stages of cooperative guidance, some scholars designed the impact time control guidance (ITCG) law to realize the salvo strike through an appropriate predesigned impact time for all missiles, which can ensure all missiles arriving at the target simultaneously. For a stationary target, Jeon et al. introduced an ITCG law which can make each missile of the group hit the target individually at the same time [1]. Lee et al. presented an ITCG guidance considering both the time and angle constrain [2]. The authors in [3] designed a cooperative guidance law based on the missile's acceleration. However, the predesigned impact time has a greater impact on the ITCG law. If the predesigned impact time is inappropriate, the missiles cannot reach the target simultaneously due to the limitation of the missile speed and maneuverability. Because no actual cooperative operation is conducted in real time, so the ITCG guidance cannot be viewed as a real-time multimissile cooperative guidance.
Subsequently, with the development of communication networks, the cooperative homing guidance based on communication has attracted much more attention, in which multimissiles attack a target simultaneously by interacting with the cooperative information in the formation [4]. Based on the estimated time-to-go, the authors in [5,6] introduced cooperative guidance laws to achieve consensus. Under time-varying communication topology, Wang et al. [7] presented a guidance with angle constraint in a leader-following coordination strategy. Zhou and Yang [8] designed a cooperative guidance law for moving targets. Furthermore, they used the missile's range-to-go as a cooperative parameter to design a cooperative guidance [9]. Based on finite time theory, the authors in [10][11][12] designed distributed cooperative guidance laws, which can realize finite time consensus. Zhao et al. in [13] combined time-to-go estimation with a cooperative term to design three-dimensional cooperative guidance based on traditional PNG, which is easily implemented. Later, Zhao and Zhou used the receding horizon control technique to design a guidance against a stationary target [14]. The authors in [15] introduced cooperative guidance laws based on the direction and normal direction of LOS considering impact-angle constraint to achieve consensus in finite time. Based on the direction and normal direction of LOS, Zhang and Yang [16] provide a fully distributed, adaptive, and optimal approach to address the problem of simultaneous attack against a maneuvering target with multiple missiles. Furthermore, they designed a cooperative guidance based on a time-varying terminal sliding mode for maneuvering target [17]. Lyu et al. [18] proposed a three-dimensional finite-time cooperative guidance for multiple missiles without radial velocity measurements. Based on the ITCG guidance theory, the designed cooperative guidance law in [19] requires each missile to communicate with other missiles; that is, communication topology must to be strongly connected. The authors in [20,21] designed distributed cooperative guidance considering the impact-angle constraint, which can ensure all missiles attack the target at a desired angle simultaneously. Zhang et al. used a time-varying navigation parameter with traditional PNG to design the cooperative guidance [22]. Based on fixed time theory, Chen et al. provided a fixed-time robust cooperative guidance [23]. Chen et al. [24] designed a cooperative guidance for multiple powered missiles. In view of the characteristics of the communication topology, Zhou et al. proposed cooperative guidance laws in discrete-time communication and sampled data communication, respectively [25,26].
The above results are applicable to the cooperative terminal guidance phase, which can ensure that all missiles hit targets at the same time under the directed graph, and the differences of relative ranges and times-to-go can be zero. However, if the initial conditions of the cooperative terminal guidance for medium-range and long-range missiles do not meet the requirements, the final cooperative strike task cannot be completed. For example, the relative range from the target and the time-to-go of one missile are 10 km and 20 s, respectively, and the relative range from target and the time-to-go of another missile are 20 km and 35 s, respectively. Under this circumstance, it is difficult for cooperative terminal guidance law to guarantee the completion of the task of hitting the target at the same time within the limited fight time. Therefore, cooperation in the midcourse guidance phase is very important and provides good initial conditions for the cooperative terminal guidance phase. The cooperative midcourse guidance (CMG) law is an effective way to ensure the initial conditions of the terminal guidance phase to be satisfied. Based on the Dubins path planning method, Zeng et al. [27] designed a cooperative strategy to synchronize the arrival time of the missiles in the midcourse guidance phase with a distinct leaderfollowing framework, but the missile needs to fly around a circular trajectory, which consumes a lot of energy. In [28], Wu et al. presented a CMG law to realize a cooperative task, but suitable for centralized network. However, none of the above research results consider the influence of communication delay on the cooperative guidance law.
In the real battlefield environment, communication delay is inevitable and has a greater impact on the distributed cooperative guidance law. He et al. introduced a three-dimensional cooperative guidance law with communication delay [29]. Liu et al. proposed a cooperative guidance law with communication delay for a moving target [30]. Zhaohui et al. introduced a fixed-time cooperative guidance law with input delay for a moving target [31]. Zhang et al. presented a cooperative guidance law with communication delay based on the time-to-go [32]. The above studies all aim at undirected communication topology structure which requires two nodes to communicate with each other. Compared with the undirected graph, the directed graph requires a directed path between two nodes, which is more suitable for practical engineering. The authors in [30,33] proposed cooperative guidance laws with communication delay based on the directed graph and obtained the upper bound of communication delay through the linear matrix inequalities (LMI) approach, but the LMI approach is more complicated and may not have a solution. It should be noted that the above research results are all based on the cooperative terminal guidance phase with communication delay and did not pay attention to the cooperative midcourse guidance phase. Therefore, this paper, in the case of cooperative midcourse guidance phase with communication delay, proposes a cooperative midcourse guidance law with communication delay to overcome the influence of communication delay. Furthermore, the LMI method is not convenient to apply in the engineering practice, so it is necessary to study a method that can easily obtain the upper bound of the communication delay.
Based on the above analysis, the distributed cooperative midcourse cooperative guidance law based on the directed graph with communication delay is worth studying, and the upper bound of the communication delay can be easily obtained to facilitate practical applications. Therefore, a novel DCMG law with communication delay is investigated under the directed graph in this paper. The main contributions of this paper are summarized as follows: (A) Compared with Refs. [29][30][31][32], this paper proposes a DCMG law with communication delay against a maneuvering target under the directed graph in the three-dimensional planes. (B) Compared with Refs. [29][30][31], this paper solves the communication delay problem under the directed communication topology and has more practical engineering significance than the undirected graph. (C) Refs. [30,33] need to solve the LMI in real time to obtain the upper bound of the communication delay, which places higher requirements on the computing power of the system. The proposed guidance law can directly give the maximum communication delay through the communication topology, which is more conducive to engineering realization.

2
International Journal of Aerospace Engineering The remainder of this paper is organized as follows: a few preliminaries and some conclusions are introduced in Section 2. In Section 3, the proposed DCMG law with communication delay is proposed, and its time-space consensus under the directed communication topology is proved, and the upper bound of the communication delay is given explicitly. Simulation results are given in Section 4. Finally, the conclusions are drawn in Section 5.

Preliminaries
In this section, a few preliminaries are given. Some definitions and important conclusions are presented, and the investigation goal is introduced.
2.1. Graph Theories. Suppose the graph G with n nodes of missile is the communication topology of the multimissile system, which consists of a pair ðV, EÞ, where V = fv 1 , v 2 , ⋯, v n g represents the set of nodes, and E ∈ V × V indicates edges between nodes. An edge ðv i , v j Þ presents that the ith node is connected with the jth node. The adjacency matrix A = ða ij Þ ∈ R n×n is defined as a ij = 1 if ðv i , v j Þ ∈ E and 0 otherwise. If there exists at least one node having a directed path to all other nodes, then the directed graph is said to include a directed spanning tree. The Laplacian matrix of G is L = ðl ij Þ ∈ R n×n , where Obviously, the matrix L = ðl ij Þ ∈ R n×n has one zero eigenvalue λ 1 ðLÞ = 0, and all the other eigenvalues λ i ðLÞ ði = 2, ⋯, nÞ are positive. In this paper, the following assumption is made about the communication topology. Assumption 1. The directed graph G has a directed spanning tree.

Problem Formulation.
In this part, the threedimensional geometry between several missiles and a single maneuvering target is considered, and some assumptions are presented as follows. The engagement geometry in three-dimensional planes is depicted in Figure 1, where ðX I , Y I , Z I Þ denotes the inertial system; V i and V T indicate the speeds of the ith missile and target, respectively; and M i and T indicate the ith missile and target, respectively.
In order to ensure the attack effect, the DCMG law is required to provide better initial conditions for cooperative terminal guidance law. The seeker of the missile cannot work in the midcourse guidance phase, so the missile cannot detect and track the target directly, but through external detector, the missile still receive the location and velocity information of itself and the target, as well as the estimated acceleration information of the target obtained through observation and estimation. Therefore, based on this information, which can be directly obtained in the midguidance phase, we will design the DCMG law in the inertial system. Suppose the dynamics of the missile are given as where ðx mi , y mi , z mi Þ, ðV xi , V yi , V zi Þ, and ða xi , a yi , a zi Þ indicate the position, velocity, and acceleration of the ith missile in the inertial system, respectively. In order to ensure the consistency of the DCMG law, the design of the communication topology plays a crucial role. Compared with the undirected graph, the directed graph can effectively save communication resources, so it is more suitable for practical application. In addition, the communication delay cannot be ignored. Taking into account the above factors, this paper designs a DCMG law with communication delay to ensure that the requirements of cooperative terminal guidance are met. In the cooperative terminal guidance phase, the goal of the cooperative terminal guidance law is to ensure that all missiles reach the target at the same time, that is, where R i ðtÞ presents the relative range between the target and the ith missile, t goi indicates the time-to-go of the ith missile, and t f represents the end of the terminal guidance Figure 1: Engagement geometry in a 3-D plane for a single missile.
3 International Journal of Aerospace Engineering phase. If the differences of R i ðtÞ and t goi ðtÞ of missiles are too large at the beginning of the terminal guidance phase, the cooperative terminal guidance law cannot ensure that all missiles reach the target at the same time in a limited time. Therefore, in order to ensure the effectiveness of the cooperative terminal guidance law, the differences of the R i ðtÞ and t goi ðtÞ of missiles should not be too large at the end of cooperative midcourse guidance phase. Based on the above analysis, the cooperative guidance law must meet the following conditions, that is, where t mf represents the end of the cooperative mid-course guidance phase and σ 1 and σ 2 are positive constants.

Guidance Law Design
In this section, the DCMG law with communication delay is presented and proof of consensus is given.
3.1. Guidance Law Scheme. According to the Figure 2, the DCMG law with communication delay mainly consists of the basic guidance law and the cooperative term based on the distributed consensus protocol with communication delay, which can be expressed as where a bi ðtÞ and a ci ðt − τÞ represent the basic guidance law and the cooperative term with communication delay, respectively. Considering the limitation of the type of information, the accelerations are designed based on the X, Y, and Z axes in the inertial coordinate system. In order to provide the better initial conditions for the terminal guidance phase, we take the location of the virtual target collision points as cooperative parameters, which are gradually becoming consistent through the multiagent consistency protocol with communication delay. In the midcourse guidance phase, due to the impact of the cooperative terms, the proposed DCMG law with communication delay can ensure that the R i ðtÞ and t goi ðtÞ are consistent when the terminal guidance phase begins. Due to the limited space, this paper only discusses the acceleration design of the X axis, and the same applies to the Y axis and the Z axis.

Guidance Law Design.
In the midcourse guidance phase, we use the virtual target collision point as the final target point. The x d mi ðtÞ represents the position of virtual target collision point corresponding to the ith missile and can be expressed as where x T ðtÞ represents the position of the target and V Tx ðtÞ denotes the speed of the target. Inspired by the trajectory shaping guidance law in [28,34], the acceleration of the ith missile can be given as where V Tx ðt f Þ and V xi ðt f Þ present the terminal velocities of the target and the missile, respectively; x mi ðtÞ indicates the position of the ith missile; a Tx ðtÞ indicates the acceleration of the target; P i ðtÞ is a cooperative coefficient for the cooperative term; p i is a coefficient to be determined; and t goi ðtÞ can be expressed as   International Journal of Aerospace Engineering x d mi ðtÞ is the cooperative parameter expressed as where N i ðtÞ presents the number of the ith node and its neighbor.
Remark 3. When the missile is in the mid-course guidance phase, it is considered that the angle between the speed of the missile and the line of sight of the missile is very small, so the time-to-go estimation as described by Equation (9) is feasible for midguidance guidance phase. In addition, the cooperative guidance law based on Equation (9) in Refs. [12,15,18] can also effectively complete the task of striking maneuvering targets at the same time.
From the Equation (7) where a ij represents the ði, jÞ element in the adjacency matrix and τ indicates communication delay. Combining Equations (7) and (8) with (11) ðtÞ − x mi d ðtÞÞ converges to 0 with t ⟶ ∞, Equation (7) can ensure the ith missile and its neighbor arriving at the handover area simultaneously, that is, Based on the Taylor series, we can get equations as follows: x mi Neglecting the high-order terms of Δt and substituting Equations (15) and (16) into Equation (14), we can obtain Based on Equation (7), x mi d ðtÞ and _ x mi d ðtÞ can be presented as 5 International Journal of Aerospace Engineering where t goi ðtÞ can be expressed as From Equations (17)- (19), we have where R i ðtÞ can be expressed as Based on the Equations (7) and (23) and where From Equation (24), we can obtain where (21) and (26), we have Obviously, Equation (7) can also be rewritten as where From Equations (27)- (29), we can obtain When ðtÞ ⟶ x mi d ðtÞ, so we can have where As shown in Figure 1, x T ðtÞ > x mi ðtÞ can be guaranteed easily in the midcourse guidance phase, and x mi ðtÞ − x T 6 International Journal of Aerospace Engineering ðtÞ < −R t < 0, where R t represents terminal guidance range. V Tx ðtÞ t goi ðtÞ is smaller than jx mi ðtÞ − x T ðtÞj, so we can get Since V xi − V T − a Tx t go > 0 and V Tx ðtÞ > 0, we can obtain Obviously, during the midcourse guidance phase, _ t goi < 0. From Equation (35), we can get According to Equation (24) and _ t goi < 0, we have From Equation (28), when a ci ⟶ 0, a i and can be approximated as a i ≈ a bi , so we can get From Equations (35) and (38), we can obtain Obviously, according to the Equations (34), (36), and (39), we can get Naturally, when p i > 0, it satisfies Equation (31). So when p i is a positive constant, then η i ðt + ΔtÞ < η i ðtÞ, that means η i ðtÞ is positive monotone decreasing function. Furthermore, η i ðtÞ ⟶ 0 when t ⟶ ∞, according to Equation (12) (7) and (8) can ensure that x d mi ðtÞ ⟶ x mi d ðtÞ. However, Equation (29) is more complicated and not convenient for engineering application. For medium-and long-range missiles, the speed is faster than the target (e.g., airborne warning and control system) in the midcourse guidance phase, and the acceleration of the target is smaller, so V xi − V T − a Tx t go > 0 can be guaranteed. When t ⟶ t mf , the influence of the cooperative term a ci ðtÞ gradually decreases, then a ci ðtÞ ⟶ 0. Therefore, Equation (29) can be simplified to Equa-tion (40), which is convenient for practical engineering applications.
Based on Theorem 4, we can see that selecting the appropriate p i can ensure that x d mi ðtÞ ⟶ x mi d ðtÞ, but in order to ensure multiple missiles to arrive at the handover area simultaneously, it still needs to ensure that x mi d ðtÞ ⟶ x mj d ðtÞ. When there is a directed spanning tree in the communication topology, the first-order multiagent consensus protocol can ensure that x mi d ðtÞ ⟶ x mj d ðtÞ. However, in the actual environment, communication delay must exist and its impact cannot be ignored. Therefore, considering the communication delay, we have the following theorem and lemma.
Lemma 6 (see [35]). Suppose that Q is any n × n positive definite matrix, X and Y are n-dimension vectors, then the following matrix inequality holds: Theorem 7. Suppose that there exists a directed spanning tree in graph G, the adjacency matrix A = ða ij Þ ∈ R n×n is a constant matrix, the diagonalized matrix of its Laplace matrix L responding to A is J ∈ R n×n , and the communication delay τ between any two nodes i and j. If τ < 1/λ max ðLÞ, then the consensus protocol with communication delay as Equation (11) can reach agreement, that is, x mi d ðtÞ ⟶ x mj d ðtÞ. Further, the controller as Equation (11) can make x d mi ðtÞ converge into the kernel space fðτJ 2 − JÞTg, where T is a conversion matrix for the diagonalization.
Proof. Based on the graph theory, Equation (11) can also be expressed as where Define the Lyapunov-Krasovskii function as where P is a n × n positive matrix to be determined and YðtÞ can be expressed as According to the Equation (43), _ VðtÞ can be expressed as From the Equation (45), we have Similarly, Y T ðtÞP _ YðtÞ can also be obtained. Furthermore, the third term of the right side of Equation (46) can be calculated as follows: Let θ′ = t + θ, Equation (48) can further be transformed as Substituting Equations (47) and (49) into Equations (46), we can obtain According to the Lemma 6, we have where Q is any appropriate dimensional positive-definite matrix. Hence, we can let Q = P, then Equation (51) can be expressed as where Since the Laplacian matrix L is a diagonalizable matrix, so L can be expressed as L = T −1 JT, where J ∈ R n×n is a diagonal matrix composed of the eigenvalues of the Laplacian matrix L. Because P is a positive matrix to be determined, we can let P = T T T, and substituting it into Equation (53), we can further calculate _ VðXÞ as Since the Laplacian matrix L is a positive semidefinite matrix and has a eigenvalue is 0, then diagonal matrix J ∈ R n×n can be expressed as Since λ i ðLÞ > 0ði = 2, ⋯, nÞ, if τ < 1/λ max ðLÞ, we have From Equations (55) and (58), we can obtain ☐ According to the Lasselle invariant set theorem, the state of the system as Equation (42) converges into the set Ω = fX ∈ R n | _ VðXÞ = 0g. In the following, we will divide the inequality as Equality (59) into three cases to prove the system as Equality (11) to reach an agreement.
Case 1. When ðTXðtÞÞ T ½−J + τJ 2 TXðtÞ < 0, then _ VðXÞ < 0. According to the Lyapunov stability theorem, the system as Equation (42) is stable, which means the system as Equation (11) is stable. Hence, the protocol can ensure the multiagent system as Equation (11) to reach an agreement.
where ζ ∈ R. According to Equation (60), we can obtain where span Π f g = span T −1 On the other hand, spanfΠg is a one-dimensional subspace, so we can have This completes the proof.
Remark 8. In Refs. [29,31], the communication topology is required to be an undirected connected graph. However, in this paper, the proposed consensus protocol with communication delay as Equation (11)) requires the existence of a directed spanning tree in the communication topology and has lower requirements for communication resources than the undirected graph. Furthermore, the method proposed in this paper can directly obtain the upper bound of the communication delay through the communication topology, which is simpler and more conducive to practical engineering applications than the LMI method.
As can be seen from Theorems 4 and 7, the proposed DCMG law with communication delay as Equation (12) ensures that multiple missiles reach the same virtual target collision point. In the following, we will prove the time con-sistency of the proposed DCMG law with communication delay. Proof. From Equations (7) and (9), we can obtain According to Equation (64), we can get So, t goi ðtÞ − t goj ðtÞ can be expressed as According to Theorem 7, when t ⟶ ∞, we can obtain x d mi ðtÞ − x d mj ðtÞ ⟶ 0. Therefore, from Equation (66), we can get t goi ðtÞ − t goj ðtÞ ⟶ 0. This completes the proof.☐ Based on the discussion and analysis of the above theorems, we can obtain Theorem 10.
Theorem 10. Suppose that the conditions are same as Theorem 7, the proposed DCMG with communication delay is given as Equation (12) and p i > 0; then, all missiles can arrive at the mid-to-terminal handover area simultaneously.
Proof. According to Theorem 4, when p i > 0, the proposed DCMG law can ensure that x d mi ðtÞ ⟶ x mi d ðtÞ. Due to the directed graph G including a directed spanning tree, we can obtain x mi d ðtÞ ⟶ x mj d ðtÞ with communication delay from Theorem 7. Furthermore, according to Theorem 9, we can obtain t goi ðtÞ − t goj ðtÞ ⟶ 0 when t ⟶ ∞.So from Theorems 4-9, the proposed DCMG law with communication delay can ensure that ðx d mi ðtÞ − x d mj ðtÞÞ ⟶ ð x mi d ðtÞ − x mj d ðtÞÞ ⟶ ð t goi ðtÞ − t goj ðtÞÞ ⟶ 0, which presents multimissile arrive at the mid-to-end handover area simultaneously. This completes the proof.☐

Simulation
In this section, the performance of the proposed DCMG law with communication delay is verified through numerical simulations. Considering the situation of which five missiles cooperate to strike the same target, three cases simulation experiments were carried out. In Case 1, there is no communication delay between missiles. In Case 2 and Case 3, the communication delays are different; one is lower than the     The initial location and speed of the target are 80000, 5000, and 1000 m and 300, 20, and 10 m/s, respectively. Assuming when the distance from the missile to the target reaches 20 km, the midcourse guidance phase will terminate. The initial conditions of missiles are shown in Table 1, and the target is supposed to be a maneuvering target with a T = ð20 sin ð0:2tÞ, 20 sin ð0:2tÞ, 20 sin ð0:2t + 1:57075ÞÞ, respectively. The directed communication topology is shown in Figure 3; the upper bound of the communication delay can be obtained as τ * = 0:33 s from Theorem 7. In the cooperative midcourse guidance phase, the acceleration commands of the missiles are restricted as ja i j < 10g, i = 1, 2, ⋯, n. Since the speed of the missile is fast in the midcourse guidance phase (greater than 1000 m/s), it is assumed that at the end of the midcourse guidance phase, the maximum difference of relative ranges is less than 1 km ðσ 1 = 1000 mÞ and the maximum difference of times-to-go is less than 1 s ðσ 2 = 1sÞ, which meets the requirements of the cooperative midcourse guidance phase. In the following simulation experiment, the delay caused by the actuator is considered, so the acceleration command can be expressed as where a rt ðtÞ is the actual acceleration command and T is the time constant and setting as T = 0:1 s.

Case 1 (without Communication Delay).
The simulation results, including the fight trajectory, the relative ranges between the missiles and the target, the times-to-go of missiles, and the accelerations of missiles, are depicted in Figures 4(a)-4(f). The proposed DCMG without communication delay can ensure that all missiles almost simultaneously arrive at the handover area about 20 km from the target seen in Figures 4(a)-4(c). Furthermore, it can be seen from Figure 4 that due to the influence of the virtual collision point, at the end of the midcourse guidance phase, all the missiles form a surrounding situation on the target, which provides good initial conditions for the terminal guidance phase. From Figures 4(b) and 4(c), it can be seen that the times-to-go and the relative ranges of missiles converge gradually. As can be seen from Table 2, the maximum differences of relative ranges and times-to-go gradually decrease from 5238 m to 470 m and 21.20 s to 0.54 s, respectively, which meets the requirements of the cooperative midguidance phase. According to Figures 4(d)-4(f), in order to achieve consistency, the acceleration commands of all missiles are larger in the initial phase.

Case 2 (with Communication
Delay τ = 0:2 s). The simulation results, including the fight trajectory, the relative range between the missile and the target, the times-to-go of missiles, and the accelerations of missiles, are depicted in Figures 5(a)-5(f). Because the selected communication delay is less than its upper bound, the proposed DCMG with communication delay can ensure that all missiles almost simultaneously arrive at the handover area about 20 km from the target shown in Figures 5(a)-5(c). Furthermore, it can be seen from Figure 5(a) that due to the influence of the virtual collision point, at the end of the midcourse guidance phase, all the missiles form a surrounding situation on the target, which provides good initial conditions for the terminal guidance phase. From Figures 5(b) and 5(c), it can be seen that the times-to-go and the relative ranges of missiles converge gradually. As can be seen from Table 3, the maximum differences of relative ranges and times-to-go gradually decrease from 5238 m to 512 m and 21.20 s to 0.71 s, respectively, which meets the requirements of the cooperative midguidance phase. Through Tables 2 and 3, the proposed DCMG law with communication delay can not only meet the requirements of the mid-course guidance section but also the maximum differences of relative ranges and times-to-go of missiles in Case 2 are almost the same as in the Case 1 without communication delay. Comparing , it can be seen that the differences of times-to-go and the relative ranges of missiles are significantly larger. As can be seen from Table 4, the maximum differences of relative ranges and times-to-go gradually decrease from 5238 m to 1201 m and 21.20 s to 1.32 s, respectively, which cannot meet the requirements of the cooperative midguidance phase. Through Tables 3 and 4, the proposed DCMG law with communication delay cannot meet the requirements of the cooperative midguidance phase when the communication delay exceeds the upper bound. Comparing Figures 4(d)-4(f) with Figures 6(d)-6(f), it can be seen that in order to overcome the influence of communication delay, the acceleration command with communication delay in Case 3 is larger than that in Case 1, but not enough to meet the requirements of the cooperative midcourse guidance phase.

Conclusions
In order to provide better initial conditions for the cooperative terminal guidance phase, a DCMG law with communication delay against a maneuvering target was proposed by combining the trajectory shaping guidance with the cooperative term based on the distributed consensus protocol with communication delay. Then, the time-space consensus of

Data Availability
No data were used to support this study.

Conflicts of Interest
The authors declare that they have no conflicts of interest.