#### Abstract

In order to conduct saturation attacks on a static target, the cooperative guidance problem of multimissile system is researched. A three-dimensional guidance model is built using vector calculation and the classic proportional navigation guidance (PNG) law is extended to three dimensions. Based on this guidance law, a distributed cooperative guidance strategy is proposed and a consensus protocol is designed to coordinate the time-to-go commands of all missiles. Then an expert system, which contains two extreme learning machines (ELM), is developed to regulate the local proportional coefficient of each missile according to the command. All missiles can arrive at the target simultaneously under the assumption that the multimissile network is connected. A simulation scenario is given to demonstrate the validity of the proposed method.

#### 1. Introduction

Saturation attack, which involves simultaneous attack from different missiles in a communication network, is an important combat manner to penetrate the missile defence system. In fact, a group of well-organized missiles of low cost and poor performance may yield better results than a single excellent one. The key to cooperative guidance of multimissile system is that all missiles reach the target at the same time in the case of saturation attacking.

Cooperative control theories have been researched broadly with respect to different agents, such as unmanned aerial vehicle [1–3], satellite [4], and some abstract objects [5, 6]. A cooperative control strategy for achieving cooperative timing among teams of vehicles which is based on coordination variables and functions was developed [7]. However, there are only a few existing literatures considering the cooperative timing problem of missiles. In [8], with a combination of the PNG law and the feedback of the impact time error, an impact time control guidance (ITCG) law for salvo attack of antiship missiles was presented and could be used to guide multiple missiles to hit a stationary target simultaneously at a desirable impact time. Based on this law, a cooperative PNG law was proposed by introducing a new concept of the variance of time-to-go of multiple missiles in [9]. With the weighting average consensus algorithm [10], a cooperative guidance scheme [11] was developed based on the ITCG law.

Different from the unmanned aerial vehicle or satellite, the maneuver flight of missile is mainly based on aerodynamic force and its guidance does not include task assignment, cooperative path planning, and path tracking. This makes it particularity hard to design a cooperative strategy for the multimissile system. In the case of a group of missiles intercepting a single maneuver target, an optimal cooperative guidance law [12] based on comprehensive cost function of missiles was derived with the constraint of a relative intercept angle. A time-cooperative guidance architecture which is a centralized coordination control form was proposed based on leader-follower strategy [13]. Using the dynamic surface control theory and disturbance observation technology, [14] developed a novel integrated guidance and control law, which can realize cooperation of impact time and flight position for multiple missiles during their cooperative attack. All these methods adopt a centralized control manner in which one missile must exchange information with all the other missiles.

Many classical literatures [15–18] described and analysed the PNG law in two-dimensional plane, and most of the three-dimensional trajectory simulations [19] generally adopted the dimensionality reduction method to decompose the movement of missile into longitudinal and lateral plane. Then PNG law was utilized to obtain the trajectory in two planes, respectively, and thereby compounded the desired three-dimensional trajectory. However, this method needs to solve some antitrigonometric functions [20] on account of the direction of coordinate axis; meanwhile the coupled problems have adverse effects on the precision. Furthermore, geometric methods mentioned above are on a necessary condition that the velocities of all missiles are constant. In order to settle these problems, this paper introduces the space vector method to generate the PNG law in three dimensions.

In order to achieve the simultaneous attack of missiles, the consensus of time-to-go (remaining flight time, i.e., arriving time) is considered in this paper. From the previous literature, it is known that the time-to-go of missiles can be changed by adjusting proportional coefficient of the PNG law. However, the relationship between the command and the local proportional guidance law is nonlinear and hard to obtain directly, especially in three dimensions, and then the fitting methods are considered. An artificial neural network (ANN), with a fixed number of numeric inputs and outputs, can be regarded as a complex nonlinear function. By learning iteratively large numbers of samples, the ANN stores the mapping relation between inputs and outputs. Due to this ability of self-learning and self-adaption, ANN is widely used in expert systems. Although ANN is of high precision, the learning speed of traditional ANN is so slow that it cannot be used in online cooperative guidance. ANN with extreme learning machine (ELM) is introduced in the paper to overcome this problem. It is notable that ELM is of strong generalization ability [21], which can learn the rules more accurately from training data; meanwhile it can obtain the appropriate output even if the input data is out of training set. Therefore, the problem of the inaccurate results caused by some special inputs of ANN is solved.

In this paper, we design a cooperative guidance strategy to achieve simultaneous attack based on expert system using ELM, which just requires that the communication network is connected. Considering that the communication between missiles might be incomplete on account of the disturbances from defense system and other environment factors, the centralized cooperative guidance strategy is liable to fail and the distributed strategy is more effective. Via a distributed protocol through the connected network, which aims at asymptotical consensus of time-to-go commands, the commands of all missiles are coordinated. Then the local proportional guidance law of each missile is regulated by the expert system according to the commands. Through the consensus algorithm of distributed cooperative guidance strategy, the proportional coefficients of all missiles are adjusted to achieve saturation attack.

The remainder of this paper is organized as follows. Section 2 describes the cooperative guidance problem of multiple missiles on three dimensions. Section 3 proposes a new distributed cooperative guidance strategy based on the expert system with ANNs using ELM. Section 4 shows the simulation results to illustrate the validity of the proposed algorithm. Finally, some concluding remarks to this paper are presented in Section 5.

#### 2. Cooperative Guidance Problem Formulation

Consider the cooperative guidance of similar missiles which are denoted by . The objective is to make all missiles arrive at a static target at the same time. First, the coordinate system which is used to describe the relative movement of missile and target in space is built, and the three-dimensional guidance model is obtained. Then, the PNG law is extended to three dimensions. At last, the cooperative guidance problem of multiple missiles is described.

##### 2.1. Movement Coordinate System

Set as the relative velocity vector of missile to the target and as the relative position vector (sight vector). Assume that the plane composed of and is instantaneous sight transfer plane at any time. and are the projection elements of on direction and vertical direction of , respectively. The relative movement coordinate system is depicted in Figure 1. The origin is located on the target. The plane is instantaneous sight transfer plane and axis is on the horizontal plane of target which points to the missile. Axis is on the plumb line of target which points up and axis is perpendicular to the plane. is the line of sight (LOS) angle and is the included angle of and . is the LOS angular velocity vector and is the command acceleration vector of missiles.

Let and denote the position vector of missiles and target and and denote the velocity vector, respectively. Then, we obtain relative position vector and relative velocity vector . Assume that the target is static in the following design of cooperative guidance law; that is, . It only needs to consider the motion of missile as follows:

##### 2.2. Three-Dimensional PNG Law

In the relative movement coordinate system, the command acceleration vector of PN is perpendicular to LOS, and the magnitude is proportional to the approach velocity and the LOS angular rate; that is, where denotes the proportional coefficient. In general, large results in sensitivity to noises, while small leads to slow the response of the navigation against a quick target. Besides, the overload imposed on the missile is limited; the choice of is usually [22].

From the three-dimensional model, it is easy to know that the LOS angular velocity vector is perpendicular to the sight transfer plane and along the direction. Analyzing the two components of , it can be found that the rotation of LOS is caused by , while only makes the missile approach the target. Therefore, the magnitude of has the direct relationship with . The two input parameters of the guidance law are both determined by relative velocity vector , so the problem is transformed into the decomposition of . It can be represented by the dot product and cross product of the relative velocity vector and position vector:

Suppose that the mass of the missile is , and the momentum is . So the angular momentum to the origin is described by velocity as

Additionally, the moment of inertia to the origin can be expressed as , and angular momentum is represented by angular velocity as

From expression (4) and (5), we can obtain the angular velocity vector as follows:

Since the magnitude of is obtained, then it is needed to discuss its direction. The aim of command acceleration is to eliminate the LOS angular speed and make the relative velocity along the LOS direction. Moreover, from the above analysis, it is known that the LOS angular rate is caused by . Thus, the command acceleration can be used to produce a velocity of opposite direction to , in order to make the relative speed of this direction equal to zero. The best choice of is parallel to and in the opposite direction. It is noted that is perpendicular to the plane and is perpendicular to LOS. So, the direction of is the same as the cross product of and . Additionally, is perpendicular to , and then it has

Compared to expression (2), the PNG law can be extended to the three dimensions as follows:

On account of the fact that the direction of is opposite to relative position vector , combined with expression (3), vector can be obtained as follows:

Considering that the cooperative guidance problem is to make a series of missiles attack the target simultaneously, the distance between missile and target is decreasing and the symbol of is negative all the time. Thus, the sign of absolute value can be removed.

Substitute (6) and (9) into guidance law (8); it has

Next, the cooperative guidance problem of multiple missiles is introduced. Almost all of the variables including position and velocity vectors are time-variant and the subscript of a variable denotes that it belongs to the th missile for convenience.

##### 2.3. Cooperative Guidance Problem of Multiple Missiles

How multimissiles cooperatively attack the ground static target is shown in Figure 2. missiles fly towards the target from different distances and directions, respectively; besides, their communication topology is undirected and strongly connected. Namely, every missile can only exchange information with its neighbors.

Then the objective of cooperative guidance is to find the appropriate acceleration command and make all missiles reach the target simultaneously; that is, . The time-to-go of missile is a monotonic decreasing function of , so we can regulate to adjust the arrival time.

In order to achieve simultaneous attack, a distributed cooperative guidance strategy is adopted as shown in Figure 3. The protocol makes the expected of all missiles achieve consensus asymptotically. Then, consisting of ANNs trained offline, these local expert systems transform the command to corresponding proportional coefficient . The details will be introduced in the next section.

#### 3. Algorithm Description

Denote of missile by for convenience. Since the closed loop system of (1) is autonomous, if and are known and when , then is determined by the function

Similarly, when the required time-to-go is given, there exists corresponding that can be obtained by

Figure 4 shows some curves of with respect to . Although it is hard to derive the analytic form of the function and its inverse , fortunately it can resort to the ANNs.

##### 3.1. Extreme Learning Machine

ELM is a simple learning algorithm for single-hidden layer feedforward neural networks (SLFNs) which achieves fast learning through increasing the number of hidden nodes and obtains good generalization performance [21]. In the training process with ELM, the input weights (linking the input layer to the hidden layer) and hidden layer biases of SLFNs can be assigned arbitrarily and do not need to adjust. After being chosen randomly, SLFNs can be simply considered as a linear system and the output weights (linking the hidden layer to the output layer) can be analytically determined through simple generalized inverse operation of the hidden layer output matrices. Only by setting the number of hidden nodes is it easy to get the single optimal solution of SLFNs with ELM which has a learning speed much faster than traditional feedforward neural network learning algorithms like back-propagation (BP) algorithm while obtaining better generalization performance. ELM not only tends to reach the smallest training error but also obtains the smallest norm of weights.

Compared with BP, ELM only needs to adjust the number of hidden nodes and avoids the multiple iterations; hence, it avoids the problem of local minimum or infinite training iteration and can reach the minimum training error. Owing to the strong generalization ability of ELM, the precision of training results can be guaranteed as long as the range of sample data is appropriate.

##### 3.2. Distributed Cooperative Guidance Strategy

In order to approximate and , we build two SLFNs and train them with plentiful simulation data beforehand under ELM. Then introduce as coordination variables; the structure of cooperative guidance strategy is built as shown in Figure 3. These two SLFNs get assembled into an expert system as shown in Figure 5. Net 1 is the approximation of and used to evaluate the as the initial state under the state of missile at the beginning of cooperative guidance; Net 2 is the approximation of and used to calculate according to in real time.

The time-variant undirected graph [23] of the multimissile system is defined as , where is the set of all missiles and is a set of edges, in which the edge means that missiles and can exchange with each other. The Laplacian matrix is defined aswhere is the degree matrix and is the adjacency matrix of the graph. If is connected, there exists .

If there is no cooperation among missiles, the original satisfies

In order to make , achieve consensus, we choose the distributed cooperative protocol below for missile : where is the neighborhood (the set of missiles that can exchange information with ) of missile and is specified as

The purpose of is that all converge to and thatwhich means .

From Figure 4, we can see that decreases with increasing. Through decreasing of the missile with less time-to-go, the practical flight time of it can be added. The range of traditional proportional guidance coefficient is [2, 6]. If is specified as 6, then .

Theorem 1. *System (15) will achieve consensus asymptotically.*

*Proof. *Let and ; then (15) can be rewritten asDefine energy functionthen, the derivate isBy substituting (18) into (20), we can obtainBecause can be regarded as the eigenvector of matrix under the eigenvalue of zero, . Besides, ; then we havewhere the condition of is ; hence,which means that system (15) will achieve consensus asymptotically.

Usually, we should choose big to accelerate the convergence rate of system (15). Furthermore, when max decreases into the allowable tolerance scope at , we can stop cooperative guidance and fix . The reason is that a longer convergence process may require a wider varying range of and of the samples for these two SLFNs. This will lead to longer training time and more complicated topological structure.

#### 4. Simulation and Results

In this section, a scenario is given to illustrate the proposed cooperative guidance strategy [9]. Consider the saturation attack of missiles on a static target. The initial states of these missiles and target are listed as

The topological structure of the communication network among them is shown in Figure 6 and the corresponding Laplacian matrix is

With the algorithm of ELM and Neural Network Toolbox of Matlab, these two SLFNs can be built and trained conveniently. Each SLFN includes 250 neurons in the hidden layer and adopts activation functions “tansig” and “purelin” for the hidden layer and output layer, respectively. In order to generate enough samples, the guidance of a single missile is simulated with the fixed proportional guidance coefficient iteratively, with the initial conditions , , and varying at points

The training results of the SLFN with 250 hidden neurons using ELM are quite good. The error of mean square is 0.0069, which is close to the result of BP. However, the training time is only 4.3 seconds, which is much less than BP of 51.1 seconds. When increasing the number of hidden neurons to 500, the training time with ELM still needs only 9.5 seconds but the training results improve significantly. In conclusion, the training speed and fitting precision of ELM both have a great advantage over other ANN algorithms.

With the two trained SLFNs, setting the parameters and proportional coefficients , the simulation is conducted and results are depicted in Figures 7–9. When s, the time-to-go commands have already converged to . So we stop the cooperative guidance and fix the proportional guidance coefficients. The result is that all missiles arrive at the target almost at the same time s. It is worth noting that always. Compared with results of [9], it is obvious that the convergence rate of time-to-go is much faster, and the variation range of proportional coefficients is smaller (not more than 6), without achieving the same end value. In all, the distributed cooperative guidance strategy provides better flexibility and autonomy for each missile.

#### 5. Conclusion

In this paper, the three-dimensional cooperative guidance for simultaneous attack on a static target of multimissile system is investigated. The classic PNG law is extended to three dimensions using vector operation; besides, an indirect and distributed cooperative proportional guidance strategy, based on expert systems and first-order consensus protocol, is proposed. This strategy can achieve simultaneous attack by adjusting commands cooperatively and guarantee that the overload satisfies limitation of the missile body since the proportional guidance coefficient is always less than 6 in the simulation. Simulation results demonstrate the performance of this proposed cooperative guidance strategy. Future studies will be conducted on the cooperative guidance of simultaneous attack on a maneuvering target within finite time, control saturation, and control deviation.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgments

This work was supported in part by the National Nature Science Foundation of China (nos. 61473124, 61203081 and 61174079), Doctoral Fund of Ministry of Education of China (no. 20120142120091), Fundamental Research Funds for the Central Universities of HUST (no. 2013054), and Precision Manufacturing Technology and Equipment for Metal Parts (no. 2012DFG70640).