Abstract
To improve the probability of missile coordinated detection, this paper investigates the control method of multimissile cooperative formation. Based on the missile’s nonlinear dynamics modeling, the formation control algorithm based on the position of the leader missile (or virtual leader missile) is studied. Meanwhile, considering the influence of the communication time delay, a robust adaptive cooperative controller satisfying the time delay network is designed. The simulation verification result indicates that cooperative formation can be achieved quickly with small velocity and position tracking errors when the missiles are fully connected. Therefore, the multimissile cooperative formation control algorithm designed in this paper can meet the requirements of multimissile coordinated detection, thus improving the missile’s cooperative attack capability.
1. Introduction
With the widespread application of information and network technology in modern warfare, the cooperation of multiaircraft weapon systems has attracted increasing attention. The research and application of multimissile cooperative attacks in modern warfare become crucial. Various countries are actively studying multiaircraft coordinated detection and target tracking, multisource information fusion, cooperative guidance, and formation control technology.
As for multimissile coordinated detection, to better complete combat tasks, it is necessary to consider the mutual coordination and cooperation between multiple missiles. Therefore, the optimal design of multimissile coordinated detection is essential to multimissile cooperative attacks. In the midcourse guidance phase, the satellite system can only obtain the area where the target may appear instead of the motion information of the target, and due to the limitation of maneuverability, the detectable area of a single missile is limited, which makes it difficult to intercept the target successfully. If multimissile coordinated detection is applied, based on the cooperative design of missiles in space position or detection time, the detectable areas of multiple missiles can cover the area where the target may appear as much as possible, thereby ensuring that the target can be intercepted during the handover phase of the mid-terminal guidance. Overall, the formation flight of air-to-air missiles can overcome the shortcomings of a single missile and complete the tasks that a single missile cannot accomplish, such as expanding the detection range, improving the detection accuracy, and enhancing the coordinated detection capability. This has an important influence on the successful interception of the target during the handover phase of the mid-terminal guidance.
Multimissile formation tasks include several processes, such as initial construction, maintenance, shrinkage, expansion, and reconstruction. Moreover, there exist severe nonlinear couplings in the translational dynamics and rotational dynamics of the missiles with external disturbances. Furthermore, the communication delay is also involved in the multi-missile formation tasks. These difficulties make it challenging to design an adaptive control algorithm to address the problem of multi-missile formation flying. In order to achieve high precision formation flying mission, it is necessary to carry out effective cooperative control for each formation member. In the research field of unmanned aerial vehicle (UAVs) formation control, scholars from various countries have proposed many methods in recent years, such as the leader-follower method [1, 2], behavior-based control method [3, 4], and virtual structure method [5, 6]. At present, these methods have been widely used in robot formation [7], aircraft formation [8–10], spacecraft formation [11, 12], and missile formation [13].
For the problem of missile formation coordination, most literatures mainly centered on the cooperative guidance based on the consistent arrival time [14–16]. The main idea is to use the consistency algorithm to design the guidance law so that all missiles can reach the target at the same time, that is, to realize the coordinated attack on the target. In recent years, the research on cooperative guidance of multiple missiles has attracted wide attention of scholars. In [17], the cooperative attack with multiple missiles on a stationary target was considered, and a two-step procedure was proposed. In [18], Jeon et al. combined the classical proportional guidance law with the additional attack time error correction term and proposed the guidance law of cooperative attack target. In [19], Lee et al. designed a cooperative guidance law with attack time and angle constraints, but the attack time and angle should be given in advance. However, Jeon et al. adjusted the difference of the time-to-go of multiple missiles in the process of attacking the target and took attack time as the coordination variable to ensure that multiple missiles intercept the target at the same time in [15]. The coordination variable is the minimum amount of information required to achieve coordination between missiles. At present, in the cooperative guidance design of multiple missiles, the selection of coordination variable is mainly time [20, 21]. In [22], Wei et al. used linear quadratic optimal control theory to design the formation optimal holding controller for multimissile formations. The designed controller can ensure the rapid and stable adjustment and maintenance of missile formation during missile cooperative operation. Zhang et al. [23] designed a nonlinear optimal controller for missile cooperative attack formation by solving Riccati matrix differential equation. The formation maintenance of multiple unmanned aerial vehicles based on proximity behavior was explored in [24]. Individual decision-making was conducted according to the expected UAV formation structure and the position, velocity, and attitude information of other UAVs in the azimuth area. The problem of formation tracking control for the multiple flight vehicle (MFV) system considering actuator saturation constraints was investigated in [25]. Cooperative guidance problems for multiple missiles with directed interaction topologies are dealt with in [16]. In [26], an adaptive output-feedback neural tracking controller was designed for a class of multiple-input multiple-output system. In [27], the problem of asymptotic tracking control for a type of nonlinear stochastic systems was investigated. Although the research achievements of cooperative guidance are increasing year by year, many problems remain to be solved due to the late start.
In the design of missile formation controller mentioned above, the interference of complex external environment is not considered. However, in the course of missile flight, the external disturbances such as airflow and strong wind are often not negligible; so, the requirement of suppressing the disturbances must be considered in the design of missile formation controller [28, 29]. Considering the influence of such external disturbance is one of the problems to be studied in this paper. Meanwhile, in order to implement formation, missiles often carry out cooperative control through communication, and communication delay may occur in this process. Therefore, the influence of communication delay must be considered when designing the controller [30].
For the missile formation flight system, a robust adaptive missile formation cooperative control strategy is proposed, and the Lyapunov stability of the closed-loop system is analyzed strictly. Firstly, an adaptive law is proposed to estimate the upper bound of the disturbance without communication delay, and a missile formation control algorithm is further designed. On this basis, considering the external disturbance and communication delay, the conditions for asymptotic stability of the system are given by Lyapunov stability analysis. Finally, the effectiveness of the proposed robust adaptive controller is verified by numerical simulation.
2. Cooperative Formation Control Scheme for Multiple Missiles
To ensure that the multimissile cooperative search strategy has the highest interception probability at the beginning of the mid-terminal guidance handover, it is necessary to maintain a certain formation flight before the mid-terminal handover and during the handover process. Thus, based on the specific missile formation with the highest probability of target interception (Figure 1), the missile group formation flight control law can be derived.

Generally, formation control refers to the application of multiple missiles to generate a certain formation configuration during the formation flight and maintain or shrink it in real-time. Formation members communicate with each other and conduct cooperative work, controlling the relative distance and direction between the formation members and maintaining the required geometric configuration. From the aforementioned analysis, the distance and azimuth information of the follower missile based on the leader missile has been obtained from the aforementioned optimization problem before the mid-terminal handover. This also means that the expected formation of the missile group is known. Furthermore, through the leader-follower missile control strategy, the missile group formation control is converted to the design of a controller for each missile so that it can track the expected position and velocity information to achieve the purpose of formation maintenance.
The leader-follower missile control strategy needs to adjust the multimissile cooperative search strategy in time according to the number of missiles in the missile group, and then it changes the selection of leader missile. The specific control process is as follows.
Based on the herringbone formation of five missiles with the highest probability of target interception, the position that always points to the center of the target maneuver domain is selected as the center point for the formation control. The velocity is the expected velocity of the entire formation, and it needs to meet the velocity requirement of the seeker for interception. Other missiles track the relative position and velocity of the formation center following the cooperative search strategy to achieve the purpose of multimissile formation flight. For example, when five missiles are used for cooperative search, is selected as the leader missile, and are selected as the follower missiles. always points to the center of the target maneuvering domain, and the position of is calculated based on the optimization problem of the maximum interception probability. In the multimissile formation flight control, let be the leader missile, and its expected information tracks the center point of the formation. separately tracks the velocity and the expected relative position of . For instance, when the four missiles are used for cooperative search, the maximum probability of intercepting the target may be achieved through the overlap and coverage of the seeker’s field of view. At this time, the virtual leader missile is selected to always point to the center of the target maneuvering domain, and other missiles track the velocity and the expected relative position of the virtual leader missile.
3. Design of the Control Algorithm for Multimissile Cooperative Formation
3.1. Relative Nonlinear Dynamic Model of Formation Movement
In three-dimensional space, a formation flight system is composed of missiles with close velocity. The movement of the -th missile can be described as where represents the position of the missile; , , and represent the missile velocity, ballistic inclination, and deflection angles of the missile, respectively; is the mass of the missile; is the resistance suffered by the missile, and it is related to the resistance coefficient , the current air density , the velocity of the missile, and the characteristic area of the missile. The calculation formula of is
Let represent the missile control force vector, and represent the disturbance suffered by the missile. Assume that is bounded and satisfies , also, let represent the position of the missile, and represent the velocity of the missile. Then, formula (1) can be expressed as where
For the convenience of analysis, in the above formula represents the interference component on each channel.
3.2. Description of Formation Tracking Control Problems
In missile formation flight control, the formation is usually maintained according to actual conditions. Assume that the expected position of the -th formation missile is , is the expected position of the formation center, i.e., the position of the leader missile or virtual leader missile, and it corresponds to the center of the target maneuvering domain. is the expected position of the -th missile relative to the center of the formation. It should be noted that and mean that the tracking requirements are fulfilled. Meanwhile, we have and , indicating that the formation maintaining requirements are fulfilled. Therefore, based on this principle, the controller is designed for the missile formation motion model, such that, for , when , we have and . Meanwhile, in the transient convergence process of the tracking errors of and , it is guaranteed to a certain extent that and .
In the missile formation system, the missile obtains the status of other missiles through information interaction to complete the formation task. This information interaction can be described by graph theory. The graph is composed of a node set , an edge set , and a weighted adjacency matrix . In the graph, the edge refers to the information connection from node to node . As for the element in , if , then ; otherwise, . As for the path in the graph, for any , . If the edges in the graph are not directed, and any two nodes can be connected in turn by paths, this graph is called a connected undirected graph. If the edges in the graph are directed, and any two nodes can be connected in turn by paths, this graph is called a strongly connected directed graph. If the direction of the edges is not considered, and any two nodes can be connected at the same time by paths, this graph is called a weakly connected directed graph.
In the problem of missile formation control, there may be time delays in the communication between multiple missiles under the impact of situation-level information. Therefore, the formation controllers are designed for the case of no communication delay and communication delay.
To facilitate the analysis, the following lemmas about graph theory are provided.
Lemma 1. The Laplace matrix of the graph G is defined as If is a strongly connected directed graph, there is a vector , in which each vector element satisfies , such that .
3.3. Formation Controller Design
To obtain the expected configuration of the missile formation to complete special tasks, each missile needs to track its expected position and expected velocity . Firstly, a robust adaptive cooperative controller with no communication time delay is given. Then, considering the influence of the communication delay, a robust adaptive cooperative controller with time delay is designed.
Denote the tracking error as , then we have . The linear sliding surface is selected as
where is a constant. For the convenience of stability analysis, the following lemma is provided.
Lemma 2. For any real number and nonzero real number y, the following inequality holds. where . The minimum value satisfies , and satisfies the equation .
Assume that the mass of the missile is known and there is no communication delay in the information interaction between the missiles, that is, the missile can obtain the information of adjacent missiles in real-time. According to formulas (3) and (6), we have
To deal with unknown disturbances in the multimissile system, an adaptive parameter is introduced, and a robust adaptive controller is designed as where is a constant, represents the element in row and column of the weighted adjacency matrix , is the estimated value of the disturbance upper bound , and , is a constant, is a time-varying parameter, and , and the constant is as defined in Lemma 2. Substituting formula (9) into formula (8), the closed-loop equation of the system can be obtained as
Theorem 3. For the missile formation system (3), robust controllers (9)–(11) are designed. If the communication topology between missiles is a strongly connected directed graph, then the system states and are globally consistent and bounded. Besides, when , and hold. Also, and hold. In this case, the missile can track the expected trajectory while forming and maintaining the expected formation.
Proof. The inequality (7) in Lemma 2 can be written as
According to formula (10) and , for any , holds. Considering that , , and formula (13), we have
where , , and represent the -th components of vectors and (), respectively.
The selected Lyapunov function is
where is a constant, and it is defined in Lemma 1. Based on the derivation of formula (15) and the use of formulas (10)–(12), we have
Substituting formula (14) into formula (16), equation (17) can be obtained.
where . Since the communication topological graph between missiles is a strongly connected directed graph, the following formula can be obtained from Lemma 1.
It can be seen from formula (18) that is bounded. Thus, the system states and are globally consistent and bounded.
Because is bounded, . Since the disturbance is bounded, can be obtained from formula (12). Besides, since is bounded and , we have and . According to the Barbalat Lemma, when and , it can be obtained from the definition of that and . Since and , holds, and so does . Thus, the robust formation controllers (9)–(11) can track the expected trajectory of the missile, while forming and maintaining the expected formation. It can be known from formula (11) that the parameter term decreases monotonously. Too small makes the hyperbolic tangent function of the controller close to the sign function. However, the real-time switching of the sign function can easily cause high-frequency vibrations in the unmodeled dynamics of the system, thereby resulting in system instability. In this case, a smaller can be selected to make smaller so that the parameter decreases more slowly. Based on this, the designed controller is a continuous function with good robustness, and it can simultaneously suppress disturbances and model parameter uncertainty.
Remark 1. In this section, the robust adaptive control scheme is designed to address the multimissile formation problem subject to nonlinear dynamics, uncertain parameters, and communication delay. An adaptive parameter is introduced to identify unknown configurations of the control system. Moreover, the closed-loop stability of the derived control system is proven based on Lyapunov theory. Unlike the common sliding mode method or other traditional robust control methods, the proposed approach is continuous which is simple to be implemented.
4. Simulation Verification of Formation Controller
To verify the validity of the controller, the simulation verification is conducted under , , , , , and the formation parameter . The formation parameter is converted to the coordinates of the ballistic coordinate system of the leader missile . Here, the formation flight of five missiles is taken as an example. Let , where the subscript represents the -th missile. For the convenience of the simulation, the resistance suffered by the missile is ignored. Besides, considering the influence of gravity, the maximum overload of each missile in the missile group is 10 g.
The weights of the missiles are , , , , and . The missile’s initial ballistic inclination angle and ballistic deflection angle are both 0°. The other initial conditions are listed in Table 1.
In the process of missile formation flight, to quantitatively describe the accuracy of formation tracking and formation maintenance, the formation tracking error and the formation keeping error; are defined as
It can be seen that a smaller and make the formation configuration performance of the missile meet the requirements better.
The initial state and velocity of the formation missile are set to the following:
The expected position of the formation missile relative to the formation center is set to the following:
The expected position and velocity of the formation center is , .
Then, the expected position and velocity to be tracked by the formation missile is , , and .
The parameters of the controller and the corresponding initial values are , , , , , and .
The missiles in the missile group are fully connected, and the weighted adjacency matrix is as follows:
The trajectories of five missiles are shown in Figure 2. The position and velocity tracking error are depicted in Figures 3 and 4, respectively. It can be observed from Figures 3 and 4 that the position and velocity tracking error converge to the origin in 10 s. The adaptive parameters and curves are shown in Figure 5. It can be obtained from Figure 5 that converges to a certain constant ultimately and monotonically decreased. Since declines much slower than the position and velocity tracking errors converge to zero, the robustness is guaranteed, and the oscillation of the control system is avoided. Moreover, the formation tracking error defined in (20) is depicted in Figure 6. It is shown that the formation tracking error converges sharply and is less than 2 m in Figure 6, which suggests that the proposed control scheme is of high accuracy. Finally, the missile control force is depicted in Figure 7.






5. Conclusions
To meet the requirements of multimissile cooperative guidance technology, missiles perform cooperative search under a specific space configuration. Meanwhile, it is also necessary to design a multimissile formation control method that meets the constraints of the specific space configuration. In this paper, a robust adaptive cooperative controller is designed according to the requirement of coordinated detection under a specific space formation configuration. The controller can meet the formation control requirements in the fully connected time-delay network and achieve the cooperative formation fast with small velocity and position tracking errors. The future work will further investigate the design of the formation controller for the local and intermittently connected network in a complex environment.
Data Availability
No data were used to support this study.
Conflicts of Interest
The authors declare that they have no conflicts of interest.