Table of Contents Author Guidelines Submit a Manuscript
Discrete Dynamics in Nature and Society
Volume 2016 (2016), Article ID 5098784, 9 pages
http://dx.doi.org/10.1155/2016/5098784
Research Article

A Novel Optimal Control Method for Impulsive-Correction Projectile Based on Particle Swarm Optimization

1Nanjing University of Science and Technology, Nanjing 210094, China
2China Academy of Engineering Physics, Mianyang, Sichuan 621900, China

Received 25 May 2016; Revised 9 September 2016; Accepted 5 October 2016

Academic Editor: Seenith Sivasundaram

Copyright © 2016 Ruisheng Sun et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

This paper presents a new parametric optimization approach based on a modified particle swarm optimization (PSO) to design a class of impulsive-correction projectiles with discrete, flexible-time interval, and finite-energy control. In terms of optimal control theory, the task is described as the formulation of minimum working number of impulses and minimum control error, which involves reference model linearization, boundary conditions, and discontinuous objective function. These result in difficulties in finding the global optimum solution by directly utilizing any other optimization approaches, for example, Hp-adaptive pseudospectral method. Consequently, PSO mechanism is employed for optimal setting of impulsive control by considering the time intervals between two neighboring lateral impulses as design variables, which makes the briefness of the optimization process. A modification on basic PSO algorithm is developed to improve the convergence speed of this optimization through linearly decreasing the inertial weight. In addition, a suboptimal control and guidance law based on PSO technique are put forward for the real-time consideration of the online design in practice. Finally, a simulation case coupled with a nonlinear flight dynamic model is applied to validate the modified PSO control algorithm. The results of comparative study illustrate that the proposed optimal control algorithm has a good performance in obtaining the optimal control efficiently and accurately and provides a reference approach to handling such impulsive-correction problem.

1. Introduction

Impulsive-correction projectile, as a type of precision-guided munitions, produces directive force for trajectory control, which results in quick-maneuvering and precise-correcting via small impulsive rockets equipped around the cross section of the airframe. This has also attracted significant interests of research due to its short-time response, high ratio of effectiveness-consumption, and simpleness of guidance [15]. Considerable work on the design of magnitude of impulse thrust and its position arrangement [4], real-time computational algorithm [6], flight stability of impulsive-correction trajectory [7], and so on has been developed. For instance, Liu and Willms [8] provided a novel approach to obtain necessary and sufficient conditions for impulsive controllability of continuous linear dynamics exercising discrete-time actions for applying into maneuvers of spacecraft. Rempala and Zabczyk [9] developed a simple and direct proof of a version of Blaquiere’s maximum principle for deterministic fixed-time impulsive control problems. Prado [10] considered the problem of optimal maneuvers to insert a satellite in a constellation in the application of two impulses to the satellite and the objective to perform this maneuver with minimum fuel consumption. Nevertheless, due to its nondifferentiable and discontinuous characteristics of impulsive control and its limit of finite working number of impulses, that is, discrete and finite-times correction, it is impossible to implement continuous control like conventional aerodynamic-fin control for the airframe.

Nowadays, numerical optimization controls are categorized into two different classes with their own advantages and characteristics, that is, direct method based on mathematical programming and parameterization of state and control histories and indirect method grounded on solution of two-point boundary value problem (TPBVP) using optimal control principle [1113]. They are of probability to be applied for achieving an optimal performance of such problem. In general, direct method is more popular in application than indirect method due to the difficulty to obtain analytical solutions of indirect approach for nonlinear complex system. Hp-adaptive pseudospectral method, as one kind of the most popular and efficient direct methods, is combining Legendre pseudospectral method [14, 15] and Hp-adaptive method [16], which discretizes state variables and control variables on a series of Legendre-Gauss-Lobatto (LGL) points. What limits its application in impulsive system is that impulsive control variable is nondifferentiable with flexible-time intervals, which does not satisfy the Karush-Kuhn-Tucker (KKT) conditions.

To tackle such problems effectively, swarm intelligence- (SI-) based methods among those evolutionary algorithms [1720] such as genetic algorithms (GAs), simulated annealing (SA), and ant-colony optimization are becoming more popular due to their speed and accuracy qualities. They are inspired by natural phenomena, for example, the behavior of groups of birds, ant colonies, herds of animals, and even social connections between human beings [19]. As a type of SI-based methods, the particle swarm optimization (PSO), primarily introduced in 1995 by Eberhart and Kennedy [21] and then extended by other researchers [22, 23], has been showing brilliant effect in optimizing discontinuous problems for its briefness in concept, easiness to implement, and high computational efficiency. Some mathematical approaches also make contributions to the PSO. For example, Couceiro and Sivasundaram provided a modified PSO algorithm to overcome traditional PSO algorithm’s drawbacks by considering a fraction calculus approach [24]. Pires et al. proposed a novel method for controlling the convergence rate of the PSO algorithm using fractional calculus concepts and observed the relationship between the fractional order velocity and the convergence of the algorithm [25]. Reference [26] interpreted PSO as a finite difference scheme for solving a system of stochastic ordinary differential equations (SODE) and proposed a class of modified PSO iteration methods based on local attractors of the SODE which behaved differently for different problems. Unlike the traditional optimization techniques, the PSO does not rely on the rigid mathematical characteristics (continuity, derivability) of the optimization problem itself and constraints in the optimization process. Reference [27] developed a PSO approach with a punish function to design the static parameters such as the working number of impulses, the magnitude and axial eccentricity of each single impulse thrust, and the fin oblique angle for impulsive-correction projectile. Yang et al. [28] presented a new approach to a fuel-optimal impulsive control problem of the guided projectile by using an improved PSO technique; this method did not consider the optimal working modes and flexible-time intervals in detail. Reference [29] proposed a new method for solving an optimal control problem applied to spacecraft reentry trajectory by using a PSO method and avoiding the calculations needed in the common analytical approaches.

Therefore, the aim of this paper is to present a new method for solving an optimal impulsive control problem with discrete, flexible-time interval, and finite-times correction using a modified PSO method, where the Hp-adaptive pseudospectral method has not been fully solved. The remaining of this paper is organized as follows. Section 2 deduces the mathematic model of impulsive-correction projectile system and states the optimal impulsive control problem. In Section 3, a modification based on the basic PSO is presented in detail for the impulsive optimization control design and the structure and parameter design of the controller, which solve the optimal setting of impulsive control for impulsive-correction projectile. Moreover, a suboptimal control and guidance law based on PSO technique are developed for the real-time consideration of the online design in practice. In Section 4, a simulation case is shown to demonstrate impulsive-correction projectile by implementing the modified PSO algorithm. In order to validate the performance of the impulsive control design, a specific nonlinear flight dynamic model coupled with different conditions such as optimal and suboptimal algorithm is also carried out. Finally, some conclusions are presented in Section 5.

2. Problem Statement

For the convenience of discussion, the motion of impulsive-correction projectile in longitude plane is chosen in this paper under assumption that the Earth is flat and motionless. The flight dynamic equations of impulsive-correction projectile arewhere the state variables include true airspeed of the projectile , trajectory inclination angle , pitch angle , pitch angular ratio , horizontal position , height position , and mass of the projectile . and are the magnitude and axial eccentricity of the impulsive thrust, respectively; means the pitching moment; is the moment of inertia about -axis; denotes the mass flow rate; represents the attack of angle; and the drag force and lift force are both the functions of dynamic pressure , reference area , and the drag and lift coefficients , as shown in (2).

Considering the wind disturbance during the flying period of the projectile, additional equations can be governed aswhere the subscript means the relative movement between the projectile and the air stream. (+ for downwind and − for upwind) is the horizontal velocity of wind during the flight, and Mach number is the function of and sonic speed .

Compared with the time-of-flight of the entire trajectory, the working time of the impulse is so transient that it can be treated as instantaneous mutation. Meanwhile, in consideration of low-cost design, impulsive rocket is of open loop, which cannot change the magnitude of its thrust force; that is, impulsive control only stands two statuses:where denotes the constant magnitude of the impulsive force (see Figure 1) and represents the flag of the impulse actuated with order number , to the contrary . Here is the total number of the impulses equipped.

Figure 1: Schematic of impulsive control variables.

Figure 2 gives the control and guidance process of the impulsive-correction projectile. The task of this paper can be described as follows: (1) In terms of the reference model, find an optimal solution of the impulsive design variables () for minimum consumption of impulsive control and minimum control error under the condition of the flight dynamic equations described by (1)-(2), the relative terminal states and control constraints. (2) Based on the optimal numeric solution, find an online control structure and algorithm which will be illustrated in Section 3, to obtain the suboptimal control for the practical realization.

Figure 2: Control and guidance structure of impulsive-correction projectile.

3. Optimal Control Design

3.1. Reference Model

For the optimal flight control problem under consideration of (1) and Figure 2, a mathematical model based on the linear disturbance motion can be taken as reference model controlledwhereHere, , , and represent transfer coefficient, damping, and time constant of the projectile, respectively.

3.2. Objective Function

Considering the impulsive-correction projectile belongs to a kind of low-cost munitions, the purpose of optimizing impulsive control is to satisfy the requirements of the cost as well as the precision of the impulsive projectile system, that is, the minimum amount of the working number of impulses and the control error. Thus the performance index used for minimization of control energy and control error iswhere denotes the working number of actuated impulses; is the coefficient of weight; and the control error . Here and represent the initial and end time, respectively.

3.3. Flight Constraints

For the given control model, the constraints corresponding to the flight optimization are as follows. The constraint on the state control iswhere the subscripts lower and upper represent lower and upper limit, respectively.

In terms of the limit of total number of impulses, the inequality constraint is given as

Similarly, due to the discontinuity of impulsive control, the time interval (see Figure 1) should be subject to

In accordance with the above formulation for the impulsive-correction projectile, the optimization issue in question is to achieve an effective algorithm to set the schematic of the impulses of application for a minimum impulsive control number and minimum control error. Due to the discrete characteristics of impulsive control exhibited in (3) and Figure 1, which is different with that of conventional continuous control, it is well to select , time intervals between the two neighboring impulse forces and as design variables.

3.4. Optimization Design Based on Modified PSO

From the control and guidance schematic of the impulsive-correction projectile system outlined in Figure 2, the control structure can be described in detail as in Figure 3, where is the impulsive controller to be designed in terms of the reference model. In order to complete the design, two steps can be implemented as follows. (1) Obtain the numeric solution of the optimal impulsive control, which will be deducted in this section by using a numeric optimization approach. (2) Determine the suboptimal online controller via establishing the relationship between the error and the control , as expressed in Section 3.6.

Figure 3: Structure of impulse control system.

PSO belonging to a category of SI methods mimics the unpredictable motion of a group of birds while searching for food, trying to take advantage of the information sharing mechanism that affects the overall behavior of the swarm. The main strength of PSO is its fast convergence due to cooperation of all individuals in finding the best solution. The initial population that composes the swarm is generated randomly at the first iteration of the process. Considering an -dimensional search space, the swarm maintains a population of particles. Each particle in the swarm is associated with a position vector as representing a possible solution to the problem and corresponding to a specific value of the objective or fitness function, a memory of its previous best position vector as remembering its best value so far () and the corresponding position in the search space, and velocity vector as denoting the varying speed of the current solution in the search space.

Moreover, each particle knows where the best value for the fitness function has occurred so far in the group (). At the end of the process the best particle in the swarm (corresponding to the best solution with reference to the fitness function) is selected. The velocity and position update equations are given in the form of a single iteration as follows [13]: where ; is the iteration counter; and , are random numbers, generated uniformly in the range responsible for imparting randomness to the flight of the swarm. The inertial weight and the acceleration coefficients and are all nonnegative constant real parameters. Meanwhile, for improving the local ability of convergence of the particles during the range of the error, the inertial weight is proposed to be reduced linearly from the start value to the end.

3.5. Modified PSO Algorithm Process

From the above modifications, Figure 4 describes the whole computational process of the algorithm in the following steps.(1)Initialize the population of particles in the solution space of PSO such as swarm size, neighborhood size, iteration number, inertia weight, acceleration coefficient, random swarm locations, and velocities.(2)Calculate the fitness value for each particle in terms of the fitness function for updating individual optimal fitness and the swarm optimal fitness .(3)Evaluate whether the current fitness is individual optimal via comparing the fitness value of its current location to that of the local best . If the current location has a higher fitness value, should be replaced with its current location. Otherwise, we keep the previous .(4)Similarly, compare the fitness value of its current location to that of the global best at each particle along its path. If the current location has a higher fitness value, should be replaced with its current location.(5)Update the velocity of the variables in particles in terms of (10) and the position by using (11).(6)Evaluate whether the target or maximum epochs are reached; that is, judge whether the new particle generated by iteration arrives at the minimum fitness value or whether the time of iteration approaches the predefined maximum value. If any of the two conditions is satisfied, it is indicated that the optimum solution of the problem is corresponding to the global best of the particle, and the criteria should be stopped. Otherwise go to step (2) for next iterations.

Figure 4: Process of proposed PSO algorithm.
3.6. Suboptimal Online Controller Design

Regarding the above optimization process of the impulsive control, the computational time is long so that it cannot be applied to the real-time control. Consequently, for seeking an online control, the controller is considered as the form of inertial element with gains and ; that is,subject to the constraints as follows.(1)Due to the finite amount of number of the impulses, the time interval is substituted to (9).(2)In accordance with the discrete characteristic of the impulsive control, some treatments are governed aswhere represents the threshold value deciding whether a single impulse works.

It is easy to understand the goal of selecting controller parameters , , and is to make the control effect approach to that of the PSO method as far as possible. Similarly, let the error between the control outputs of the two methods as the objective function, and controller parameters can be optimized via utilizing the previously described PSO method.

3.7. Guidance Law Design

In order to make the projectile hit the target precisely, it is also necessary to design an appropriate guidance law; that is, formulate the control command with the relative motion between the projectile and the target. Traditionally, the guidance law is determined to approach the convergence of the distance and the rate of the line of sight (LOS) between the projectile and the target. For instance, in terms of classic proportional navigation (PN) law it is concluded that the projectile guidance command is in proportion to the rate of LOS.

In consideration of the limit of finite energy of impulsive projectile, it is predicted that guidance trajectory should be as similar to the nature trajectory as possible under the standard condition. Thus it is helpful to save the amount of number of the impulses which are not imperative. Here, a compensation of the influence of the gravity is considered based on classic PN method for the guidance law which is given bywhere denotes guidance command of trajectory inclination angle; is the rate of LOS; and , are PN and gravity compensation gains, respectively. Usually is chosen between 2 and 6.

According to (14), the guidance law can be written in the form of integration as follows:

4. Simulation Example and Analysis

4.1. Performance of the PSO Online Controller

In this section, a simulation example using the above method is demonstrated for the impulsive control of a certain impulsive-correction projectile system. Here, in order to validate the performance of the optimal controller, the initial conditions for the design of control optimization are set as follows. The values of , , and are set to be 0.00048, 0.156, and 0.0455. The initial and end time are and , respectively. The coefficient of weight , and the total numbers . Additionally, is chosen as unit step response with the constraints of , .

The parameters of the adopted PSO control are , , . The inertial weight is set to be reduced linearly, from 1.0 to 0.4, and the search will be terminated if the fitness value is less than 0.4, or the number of iterations reaches 150. Controller gains are , , and the threshold value is = 10000 N.

Figures 57 give the optimization results and dynamic response qualities of the reference impulsive system. From Figures 5 and 6, it is shown that the dynamic response of the suboptimal controller is similar to that of PSO approach. Meanwhile, the criterion in Figure 7 that represents the fitness value’s trend over iterations remains at the same value for a long iteration period. This indicates that the particles have reached the global optimum point in the feasible area. Although the speed of computations of PSO is rather slow, the suboptimal controller design is well rounded to make up its weakness and make it feasible in the online design and application.

Figure 5: History between and .
Figure 6: History between impulsive thruster and .
Figure 7: PSO results for minimum criterion.
4.2. Trajectory Simulation for Validation

A trajectory simulation case is also presented to validate the nonlinear impulsive control system in detail, in which the numerical integration of the trajectory model uses the fourth-order Runge-Kuta method. The simulation parameters are listed as the impulsive force = 4000 N and the working time of each impulse per time equals 0.015 s. The initial state values of the projectile system are = 0 s, = 40 kg, = 650 m/s, = 0 m, = 0 m, = 0 rad/s, = = 45 deg, = 0 deg, = 0.015 m, and = −5 m/s, and the impulsive control starts at the altitude = 3000 m during the declining period of the trajectory.

According to the above PSO method for the impulsive-correction projectile, the simulation results are illustrated in Figures 812. In the figures, the standard trajectory means natural trajectory of free-flying without considering disturbance and impulsive control; the disturbed trajectory happens under the consideration of wind disturbance, but without any impulsive control. Simultaneously, the optimal trajectory under the PSO algorithm and the suboptimal trajectory with the online algorithm are exhibited for comparative study.

Figure 8: (a) Curve of trajectory under different conditions. (b) Enlarged curve of terminal trajectory.
Figure 9: History between and .
Figure 10: History between α and .
Figure 11: History between θ and .
Figure 12: History between and .

In terms of the results shown in Figures 812, the following can be concluded:(1)From the disturbed trajectory compared with the standard trajectory depicted in Figures 8(a) and 8(b), the distance of displacement resulting from the disturbance of the wind is 218.9 m, and the miss distances corrected by the impulses are both less than 0.5 m under the optimal and suboptimal approaches, which illustrates impulsive control effectively corrects the error of range resulting from the wind disturbance.(2)Figures 911 suggest that the results of suboptimal method are similarly close to those of the modified PSO algorithm, which demonstrate the online design has a good performance.(3)As shown in Figure 12, only the amount of working number of impulses is 11, which is far less than the limit of 32 for satisfying the requirements of correction distance and miss distance. Simultaneously, it is found that there counts more number of impulses with short-time intervals working in the high altitude than that in the lower altitude with long-time intervals because the impulses working in the high altitude get more sensitive to the correction of the error, which is helpful to sufficiently utilize the maneuver ability of the projectile (i.e., save the consumption of the impulses for the same requirement of the maneuver). Also the impulses working in the lower altitude are benefit to improve the precision for the insensitivity of the maneuver.

5. Conclusions

In this paper, a novel modified PSO method is proposed to implement optimal control for a class of impulsive-correction projectile systems based on the minimum impulsive energy consumption as the optimization objective. In consideration of such a nonlinear global optimization problem including the discontinuous cost function, the modified PSO is applied to search for the optimal settings of the impulse control via selecting time intervals between the two neighboring impulse forces and as control variables. Synchronously, rounding operation is applied to solve the discrete variable of the particle for the integral requirement. At the same time, linear-varying inertial weight is adopted to improve the local ability of convergence of the particles during the range of the error. To show effectiveness of the presented PSO control method, the PSO controller is coupled in the flight dynamic model, and comparative simulations, in which wind disturbance has been taken into consideration, are conducted. Simulation results demonstrate that the proposed algorithm can obtain the optimization solution of the impulsive control of the guided projectile and make the close-loop system have a good performance with satisfying working number of impulses and miss distances.

Meanwhile, the proposed PSO control method can be extended to an impulsive and aerodynamic compound control missile, but the decoupling between the impulse force and the aerodynamic force needs further research. The real-time online programming of the PSO for impulsive control problem is also attracting more attentions, while the convergence rate of the PSO algorithm needs to be improved. These issues will be taken into consideration in our future work of developing new mathematical approaches.

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (NSFC) under Grant no. 11176012 and the Research Innovation Project for Graduate Student of Jiangsu-Provincial Ordinary University under Grant no. KLYX15-0394.

References

  1. D. Corriveau, C. Berner, and V. Fleck, “Trajectory correction using impulse thrusters for conventional artillery projectiles,” in Proceedings of the IBC 23rd International Symposium on Ballistics, pp. 639–646, Tarragona, Spain, 2007.
  2. A. J. Calise and H. A. El-Shirbiny, “An analysis of aerodynamic control for direct fire spinning projectiles,” in Proceedings of the AIAA Guidance, Navigation and Control Conference, AIAA 2001-4217, Montreal, Canada, August 2001. View at Publisher · View at Google Scholar
  3. S. Doraiswamy and G. V. Candler, “Detached eddy simulations and Reynolds-averaged Navier-Stokes calculations of a spinning projectile,” Journal of Spacecraft and Rockets, vol. 45, no. 5, pp. 935–945, 2008. View at Publisher · View at Google Scholar · View at Scopus
  4. D. Lin, J. Xu, and J. Song, “Design of pulse control parameter for projectile with ballistic correction,” Armament Automation, vol. 24, no. 3, pp. 40–43, 2005. View at Google Scholar
  5. C.-W. Li, M. Gao, W.-D. Song, and Q.-W. Guo, “Optimization of trajectory correction method of rocket projectile based on pulse jet control,” Journal of Ballistics, vol. 26, no. 4, pp. 24–29, 2014. View at Google Scholar · View at Scopus
  6. X. Tian, G. Chen, and C. Xin, “External ballistic real-time algorithm of trajectory correction projectile,” Journal of Test and Measurement Technology of NCIT, vol. 14, no. 1, pp. 44–47, 2000. View at Google Scholar
  7. Z. Wang, S. Ding, and L. Wang, “Flight stability condition for a projectile of corrected trajectory with impulse moments,” Journal of Nanjing University of Science and Technology, vol. 24, no. 4, pp. 322–325, 2000. View at Google Scholar
  8. X. Liu and A. R. Willms, “Impulsive controllability of linear dynamical systems with applications to maneuvers of spacecraft,” Mathematical Problems in Engineering, vol. 2, no. 4, pp. 277–299, 1996. View at Publisher · View at Google Scholar
  9. R. Rempala and J. Zabczyk, “On the maximum principle for deterministic impulse control problems,” Journal of Optimization Theory and Applications, vol. 59, no. 2, pp. 281–288, 1988. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  10. A. F. B. A. Prado, “Bi-impulsive control to build a satellite constellation,” Nonlinear Dynamics and Systems Theory, vol. 5, no. 2, pp. 169–175, 2005. View at Google Scholar · View at MathSciNet
  11. P. K. A. Menon, V. H. L. Cheng, C. A. Lin, and M. M. Briggs, “High performance missile synthesis with trajectory and propulsion system optimization,” Journal of Spacecraft and Rockets, vol. 24, no. 6, pp. 552–557, 1987. View at Publisher · View at Google Scholar · View at Scopus
  12. R.-S. Sun and C. Ming, “Velocity control for coning motion missile system using direct discretization method,” Discrete Dynamics in Nature and Society, vol. 2015, Article ID 716547, 10 pages, 2015. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  13. J. T. Betts, Practical Methods for Optimal Control Using Nonlinear Programming, SIAM, 1st edition, 2001.
  14. F. Fahroo and I. M. Ross, “Costate estimation by a Legendre pseudospectral method,” Journal of Guidance, Control, and Dynamics, vol. 24, no. 2, pp. 270–277, 2001. View at Publisher · View at Google Scholar · View at Scopus
  15. M. M. Khader, N. H. Sweilam, and W. Y. Kota, “Cardinal functions for Legendre pseudo-spectral method for solving the integro-differential equations,” Journal of the Egyptian Mathematical Society, vol. 22, no. 3, pp. 511–516, 2014. View at Publisher · View at Google Scholar · View at MathSciNet
  16. C. L. Darby, W. W. Hager, and A. V. Rao, “Direct trajectory optimization using a variable low-order adaptive pseudospectral method,” Journal of Spacecraft and Rockets, vol. 48, no. 3, pp. 433–445, 2011. View at Publisher · View at Google Scholar · View at Scopus
  17. K. Graichen and N. Petit, “Constructive methods for initialization and handling mixed state-input constraints in optimal control,” Journal of Guidance, Control, and Dynamics, vol. 31, no. 5, pp. 1334–1343, 2008. View at Publisher · View at Google Scholar · View at Scopus
  18. N. Yokoyama and S. Suzuki, “Modified genetic algorithm for constrained trajectory optimization,” Journal of Guidance, Control, and Dynamics, vol. 28, no. 1, pp. 139–144, 2005. View at Publisher · View at Google Scholar · View at Scopus
  19. K. Premalatha and A. M. Natarajan, “Hybrid PSO and GA for global maximization,” International Journal of Open Problems in Computer Science and Mathematics, vol. 2, no. 4, pp. 597–608, 2009. View at Google Scholar · View at MathSciNet
  20. M. Pontani and B. A. Conway, “Particle swarm optimization applied to space trajectories,” Journal of Guidance, Control, and Dynamics, vol. 33, no. 5, pp. 1429–1441, 2010. View at Publisher · View at Google Scholar · View at Scopus
  21. R. Eberhart and J. Kennedy, “New optimizer using particle swarm theory,” in Proceedings of the 6th International Symposium on Micro Machine and Human Science, pp. 39–43, October 1995. View at Scopus
  22. Y. Shi and R. C. Eberhart, “A modified particle swarm optimizer,” in Proceedings of the IEEE International Conference on Evolutionary Computation, IEEE World Congress on Computational Intelligence, pp. 69–73, Anchorage, Ala, USA, May 1998. View at Publisher · View at Google Scholar
  23. E. C. Laskari, K. E. Parsopoulos, and M. N. Vrahatis, “Particle swarm optimization for integer programming,” in Proceedings of the Congress on Evolutionary Computation (CEC '02), pp. 1582–1587, Honolulu, Hawaii, USA, May 2002. View at Publisher · View at Google Scholar · View at Scopus
  24. M. Couceiro and S. Sivasundaram, “Novel fractional order particle swarm optimization,” Applied Mathematics & Computation, vol. 283, pp. 36–54, 2016. View at Publisher · View at Google Scholar
  25. E. J. S. Pires, J. A. T. Machado, P. B. de Moura Oliveira, J. B. Cunha, and L. Mendes, “Particle swarm optimization with fractional-order velocity,” Nonlinear Dynamics, vol. 61, no. 1, pp. 295–301, 2010. View at Publisher · View at Google Scholar
  26. W. Yang, W. Wu, Y. Fan, and Z. Li, “Particle swarm optimization based on local attractors of ordinary differential equation system,” Discrete Dynamics in Nature and Society, vol. 2014, Article ID 628357, 10 pages, 2014. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  27. S. Chang, X. Cao, Z. Wang, and C. Yuwen, “Approach of optimal design of parameters for impulse-correction projectile,” Journal of Ballistics, vol. 25, no. 1, pp. 36–54, 2013. View at Google Scholar
  28. H. Yang, L. Dou, and M. Gan, “A particle swarm optimization for fuel-optimal impulsive control problems of guided projectile,” in Proceedings of the 25th Chinese Control and Decision Conference (CCDC '10), pp. 3034–3038, IEEE, Xuzhou, China, May 2010. View at Publisher · View at Google Scholar · View at Scopus
  29. A. Rahimi, K. D. Kumar, and H. Alighanbari, “Particle swarm optimization applied to spacecraft reentry trajectory,” Journal of Guidance, Control and Dynamics, vol. 36, no. 1, pp. 307–310, 2013. View at Google Scholar