Abstract

The firing accuracy of the projectile has a positive relation with aerodynamic parameters. Due to the complex dynamic characteristics of projectiles, there is an overfitting risk when a single extreme learning machine (ELM) is used to identify the aerodynamic parameters of the projectile, and the identification results oscillate transonic region. To obtain the aerodynamic parameters of the projectile accurately, an aerodynamic parameter identification model based on ensemble learning theory and ELM optimized by improved particle swarm optimization is proposed. The improved particle swarm optimization algorithm (IPSO) with an adaptive update strategy is used to optimize the weight and threshold of ELM. Combined with the ensemble learning theory, the improved ELM neural network is regarded as a weak learner to generate a strong learner. The structural parameters of the strong learner were continuously optimized through training, and an aerodynamic parameter identification model of projectile based on ensemble learning theory is obtained. The simulation results show that the introduction of the IPSO and ensemble learning theory enables the model to exhibit excellent generalization ability. The proposed identification model can accurately describe the variation of aerodynamic parameters with the Mach number.

1. Introduction

Aerodynamics is the decisive factor affecting the ballistic trajectory and flight stability of the projectile [1]. The projectile’s firing accuracy positively relates to aerodynamic parameters [2]. Currently, numerical computation, wind tunnel test, and shooting test are common technical means to obtain aerodynamic parameters [3]. The result of the numerical calculation method depends on the accuracy of the ballistic model. Still, due to the strong coupling of the projectile flight motion and external disturbance, it is difficult to obtain a completely accurate ballistic model. The wind tunnel experiment simulates the genuine flight attitude of the projectile by changing the attitude and velocity of the model, which is generally used to test and correct the shape parameters and aerodynamic characteristics of the projectile. The parameter identification method is used to process the measured data (provided by shooting test) and indirectly extract the aerodynamic parameters of projectile. Research on the aerodynamic parameter identification method has significant practical engineering application value [4]. The least squares method (LSM), maximum likelihood method (ML), Kalman filtering method (KF), and intelligent algorithm are mature algorithms in the field of parameter identification [57].

LSM [8] is a classical estimation method in aircraft parameter identification. Dunkel [9] realized the identification of aerodynamic derivative and stability derivative by LSM. However, LSM is susceptible to extreme outliers. To mitigate the effects of outliers, Su and Song [10] used recursive LSM with fading memory to improve the identification performance. Kamali et al. [11] proposed improved LSM and successfully identified Dutch roll movement parameters and segment period parameters of aircraft. Mu et al. [12] combined the model reduction technology with LSM for parameter identification. Due to the complex motion characteristics of the projectile, the traditional identification equation (based on LSM) needs to be better posed. Yang et al. [13] combined LSM with an engineering test to obtain the aerodynamic parameters of projectile.

ML has been widely used due to its unbiasedness [1417]. Carnduff and Cooke [18] applied ML to reconstruct the aerodynamic model of an unmanned aerial vehicle (UAV) with an unconventional fuselage structure. However, ML highly relies on model accuracy, and the variance is significant at high latitudes, so it is often used in combination with other methods. To improve the performance of ML, Kumar and Rao [19] combined the output error method with ML, which can accelerate the speed of convergence. Zou and Li [20] introduced the interior point optimization method into ML to reduce the severe error caused by second-order numerical differentiation.

By extending the parameters (to be identified) to system states, the problem of parameter identification can be transformed into the issue of optimal state estimation so that KF can be applied to parameter identification [2123]. In nonlinear systems, KF cannot accurately estimate the state matrix. It was later developed into the extended Kalman filter (EKF) and unscented Kalman filter (UKF). Zheng et al. [24] proved that UKF is better than EKF in projectile parameter identification because of unscented transformation (U-transformation). To reduce the computational complexity and improve the accuracy, many other improved KF algorithms have been proposed. Menon et al. [25] combined the differential vortex lattice algorithm with EKF to realize flight path reconstruction. Majeed and Kar [26] proposed adaptive UKF to improve identification accuracy. Shen et al. [27] combined EKF with the aerodynamic semiempirical method to identify derivative residuals.

Due to the complex environment and unknown external interference, getting a completely accurate aircraft model is complex. To solve the modeling error of the traditional identification method, intelligent algorithms and their variants have been widely used in aerodynamic parameter identification. Du et al. [28] combined particle swarm optimization algorithm (PSO) with real-coded genetic algorithm (GA) to obtain the resistance coefficient of the projectile. Based on the maximum likelihood criterion, Li et al. [29] applied the neural network Newton method to extract the zero-lift drag coefficient of the projectile. To improve accuracy and accelerate convergence, Wang et al. [30] introduced an elite crossover strategy into GA. Guan et al. [31] combined GA with ML to identify the zero-lift resistance coefficient of a high-speed rotating projectile based on the speed data of the projectile. Aiming at the problem that gradient descent optimization algorithm is easy to fall into local optima in traditional aerodynamic parameter identification, Han et al. [32] put forward a double backpropagation (BP) neural network, Pu et al. [33] put forward a method of sample expansion and neural network parameter online fast correction based on support vector machine (SVM), and Li et al. [34] proposed an improved teaching-learning-based optimization (ITLBO) for aerodynamic parameter identification. To avoid the initial value estimation, Yan et al. [35] proposed a derivative method for identifying the aerodynamic parameters of aircraft by the three-layer neural network. Hou et al. [36] applied a differential evolution algorithm to weaken initial value sensitivity. Ji-gang et al. [37] combined the advantage of PSO in the initial value section and the advantages of the Newton iteration method in precise iteration and successfully identified the drag coefficient of the projectile. Mohamad et al. [38] put forward the concept of dynamic parameter estimation (DAPE).

Extreme learning machine (ELM) [39] is an algorithm for training single hidden layer feedforward neural networks (SLFNs). The structural parameters (input weights and hidden thresholds) of ELM are generated randomly and require no iterative adjustment. Owing to it, ELM has low computational complexity and good real-time performance and has been widely used in cloud computing, data visualization, and random projection [4042]. Akusok et al. [41] applied ELM to identify the drag coefficient of the projectile for the first time. Affected by uncertain factors such as the actual combat environment and external meteorological conditions, the ballistic trajectory data is characterized by solid nonlinearity, time-dependent nature, and susceptibility to random noise. Randomly generated structural parameters lead to the oscillation of ELM identification results [42]. In addition, when a single ELM is used to identify the aerodynamic parameters of the projectile, all the given training samples are often used to model the global situation. In other words, it is easy to make insufficient use of the sample data and cause overfitting. In summary, using a single ELM to identify aerodynamic parameters has the following limitations: (a)The identification result oscillates (especially in the transonic region) due to randomly generated structural parameters(b)When a single ELM is used to identify projectile aerodynamic parameters, it is hard to make sufficient use of the local information of the sample data and then causes overfitting

To overcome the above problems and then accurately obtain the aerodynamic parameters of the projectile, a large number of documents are referenced. For problem (a), the classical idea is to apply PSO, GA, and other optimization algorithms to optimize the structural parameters of ELM [4349]. However, iterative optimization increases the time complexity of the algorithm. Then, the adaptive update strategy is introduced to improve the performance of PSO. For problem (b), Schapire [50] proved that multiple weak learners could generate a strong learner with good generalization performance by ensemble theory. Considering that the projectile parameter identification problem is a regression problem, AdaBoost. RT algorithm [51] is used as the integration framework. Above all, this paper puts forward an aerodynamic parameter identification model of projectile based on improved ELM and ensemble learning theory (we named IPSO-ELM-AdaBoost). The proposed IPSO-ELM-AdaBoost is a comprehensive application of multiple algorithms. In short, the functions of the hybrid algorithm can be generalized as follows: (a)ELM (function as a weak learner) establishes the mapping relationship between ballistic data and aerodynamic parameters(b)Improved PSO (IPSO) provides ELM with optimized structural parameters(c)AdaBoost. RT is responsible for integrating multiple weak learners into strong learners

The rest of the paper is arranged as follows: in Section 2, the concrete expression of the ballistic equation is given. ELM-AdaBoost. RT aerodynamic parameter identification model based on improved particle swarm optimization (IPSO-ELM-AdaBoost) is described in detail in Section 3. The simulation results under standard meteorological conditions are analyzed in Section 4. Ultimately, conclusions are summarized in Section 5.

2. Ballistic Trajectory Model

Before parameter identification, the ballistic trajectory model (6DOF) [52] must first be solved to obtain the ballistic data. 6DOF treats the projectile motion as a rigid body motion. It considers the three degrees of freedom of the projectile centroid motion and the three degrees of freedom of the angular motion. Ignoring the dynamic unbalance of the projectile, the aerodynamic eccentricity, and the Coriolis inertial force caused by the Earth’s rotation, 6DOF can be mathematically described as where is the speed; is the relative velocity; , , and are the components of relative velocity along , , and axes; is the distance; is the altitude; is the lateral distance; , , and are the projected components of rotational speed along , , and axes; is the elevation angle in coordinate system of projectile axes; is the characteristic area; is the mass of projectile; is the reference projectile length; is the reference projectile diameter; is the acceleration of gravity; is the air density; , , and are the projected components of wind velocity along , , and ; is the relative angle of attack; is the latitude; is the rotational angular velocity of the earth; is the angle of departure; is the relative angle of attack; is the elevation angle in velocity coordinate system; is the direction cosine in velocity coordinate system; is the direction cosine in coordinate system of projectile axes; is the aerodynamic malalignment of additional moment; is the aerodynamic malalignment of additional force; is the horizontal component of angle of attack; is the longitudinal component of angle of attack; is the drag coefficient; is the lift coefficient; is the lift coefficient derivative; is the Magnus force coefficient; is the rolling damping moment coefficient derivative; is the static moment coefficient; is the oscillating damping moment coefficient derivative; and is the Magnus moment coefficient derivative.

3. IPSO-ELM-AdaBoost

The proposed IPSO-ELM-AdaBoost is a comprehensive application of multiple algorithms: ELM (function as a weak learner) is used to extract aerodynamic parameters from ballistic data. PSO variants provide ELM with optimized structural parameters. AdaBoost. RT algorithm (function as an integration framework) is responsible for integrating multiple weak learners into strong learners and outputting the projectile’s final aerodynamic parameter identification results. In this section, ELM and AdaBoost. RT will be briefly introduced as prior knowledge, and the idea of adaptive update strategy in IPSO will be presented in detail.

3.1. ELM

ELM [39] is a special SLFN without iterative adjustment of structural parameters. The working process of ELM can be divided into learning and prediction. For ELM with input layer, hidden neurons, output layer, and activation function (the activation function can be any nonzero function), the structure of ELM is shown in Figure 1.

3.1.1. Training Process

Given input vector ,

The structural parameters of ELM are randomly generated as follows:

Then, the output of the hidden layer is as follows:

The output weight matrix is mathematically described as

The output can be calculated as

Given training samples , the output of ELM is as follows:

The training goal is mathematically described as

The output weight matrix can be calculated as where is the Moore–Penrose generalized inverse of .

In order to improve the generalization ability of ELM and avoid overfitting, based on the ridge regression principle [53], a regularization term [54] is introduced in (21): where is the identity matrix and is the regularization factor.

3.1.2. Prediction Process

When ELM completes the training process, the output matrix can be calculated by equation (22). Given predicting samples , the output of ELM is

3.2. ISPO

PSO is a metaheuristic algorithm proposed by Kennedy and Eberhart [55]. Compared with other metaheuristic algorithms, such as GA [56] and ant colony algorithm [57], the PSO algorithm has a simple structure, easy implementation, global solid searchability, and fast convergence speed [58]. Based on the global best position and the individual best position , the particles are iteratively updated until convergence. The steps of PSO can be described as follows: (1)Set hyperparameter

Hyperparameters that need to be set include the following: population size , particle dimension , velocity inertia weight , learning factors and , maximum iteration , minimum error , maximum particle velocity , and maximum position . (2)Initialization

Initialize position and velocity , . (3)Iterative update

The fitness function is selected to calculate the fitness value of particles and find out the individual and global optimum of particles. For optimization problems, the update rule is where specifies the individual optima of the th particle in the th iteration and specifies global optima in the th iteration.

Particle is iteratively updated by where . and are random numbers subject to uniform distribution. (4)Iteration stop

When the algorithm converges, the optimization result is output. If not, step (3) is transferred to continue the iteration.

The structural parameters of ELM that are optimized by PSO can contain more training sample information than randomly one and can effectively improve the identification accuracy. However, the iterative optimization of PSO increases the time complexity of the algorithm. The in equation (26) can keep the motion inertia of particles and make them tend to expand the search space, which has an important influence on the optimization performance of the algorithm. Dynamic can obtain better optimization results than fixed one [59]. A larger can improve the global search ability of the algorithm, while a smaller one can improve the local search ability of the algorithm. In order to improve the convergence speed of the algorithm, the adaptive update strategy is introduced in IPSO. Formula (28) calculates the average distance from the global optimal particle to other particles in the th iteration and then maps to the interval [0, 1] by formula (29) to obtain the adaptive factor . describes the state of the population. In other words, the larger the , the farther the particle is from the global optimal particle, and the particle needs a larger to quickly approximate the global optimal solution. The smaller the , the closer the particle is to the global optimal particle, and a smaller is required to limit the particle to the vicinity of the global optimal solution and improve the search accuracy. The adaptive update strategy (based on ) of can be described as

In order to verify the effectiveness of the IPSO, the standard test functions and are selected to conduct the test independently for 100 times and compared with the PSO algorithm. The expressions of test functions and are shown in equations (31) and (32), and the related parameter settings of the improved particle swarm are shown in Table 1. The results after 100 independent tests are given in Table 2. Test results show that the introduction of adaptive updated in IPSO can effectively improve the accuracy and convergence speed of the algorithm. In the case of the same number of particles, for the same test function, the results of IPSO are closer to the actual value, and the average convergence speed is faster.

The structure of ELM optimized by IPSO (IPSO) is shown in Figure 2.

3.3. AdaBoost. RT

AdaBoost. RT [60] is proposed for regression problems. The idea is to filter out examples with a relative estimation error higher than the preset threshold value and then follow the AdaBoost procedure. Its basic idea can be as follows: preset threshold , train weak learners (IPSO-ELM) based on training samples, and update the weight of training samples according to the current prediction error of IPSO-ELM. In the next round of training, the training samples with significant prediction error will have larger weights and continue to train weak learners based on the new weight. After training rounds, weak learners are obtained, and finally, the output of all weak learners is weighted to get the final prediction result. The implementation process of IPSO-ELM-AdaBoost is described as follows: (1)Initialization

Randomly select training samples from the dataset , initialize the weight of the training samples, and set the threshold of the algorithm, the initial prediction error rate , and the iteration round number . (2)Train weak learners

The th weak learner is trained by training data, and calculate the error of each training sample () and weak learner (). where is the prediction result of the th weak learner (th round) in the th training data and is the actual value. (3)Update the weight of training samples

The updated formula is as follows: where is the normalization factor. (4)Repeat the training

Repeat rounds of the step to obtain weak learners and weighting to obtain the output of the strong learner:

The flow chart of IPSO-ELM-AdaBoost is shown in Figure 3.

4. Simulation Verification

In this section, a series of simulation tests are carried out under standard weather conditions to validate the feasibility and excellence of IPSO-ELM-AdaBoost in aerodynamic parameter identification. In simulation test 1, IPSO-ELM-AdaBoost is compared with three other machine learning algorithms (ELM, IPSO-ELM, and ELM-AdaBoost), and the performance of mentioned machine learning algorithms is stated in the quantitative ways. In simulation test 2, IPSO-ELM-AdaBoost, LSM, ML, and UKF are used to reconstruct the trajectory, respectively. The performance of the four algorithms is evaluated from the point of fall and side deflection.

4.1. Data Preprocessing

Under standard meteorological conditions, 6DOF in Section 2 is solved, and 10000 ballistic data are obtained. Table 3 shows the initial launch parameters of 6DOF, and the variation laws of the , , , and of the projectile with time are shown in Figures 47.

The ballistic trajectory data contains the flight velocity, position, and attitude of the projectile, and different information has different dimensions. To eliminate the influence of different dimensions on data analysis, the min-max normalization is used to preprocess 10000 original datasets. The formula is as follows: where represents the original data, is the minimum value of data, represents the maximum value of data, and represents the normalized data.

4.2. Simulation Test 1

Figures 812 are graphs of aerodynamic parameter identification results. In the figure, the abscissa represents the Mach number (), and the ordinate represents aerodynamic parameters to be identified. The identification result of IPSO-ELM-AdaBoost can accurately fit the real aerodynamic parameter curves in the whole trajectory range. For the same parameter to be identified, the introduction of the IPSO and ensemble learning theory enables the model to exhibit excellent generalization ability. Take the identification results of (Figure 12) as an example, the identification result curve of ELM oscillates, especially in the transonic region . The identification results of IPSO-ELM and ELM show that the structural parameters optimized by IPSO contain more sample information, which leads to higher curve fitting degree. Comparing the identification results of ELM-AdaBoost and ELM, it is found that the identification result curve considering multiple weak learners can better describe the variation law of aerodynamic parameters with . It is worth mentioning that (by comparing the identification results of IPSO-ELM and ELM-AdaBoost), although IPSO can improve the performance of the model, respectively, ISPO-ELM still has the risk of overfitting because it only takes a single learner into consideration. For different aerodynamic parameters (to be identified), the linearity between aerodynamic parameters and affects the fitting degree. The stronger the linearity between aerodynamic parameters to be identified and , the higher the curve fitting degree of identification results.

Table 4 shows the model structure and average identification time (after 100 independent tests) of the four algorithms. From the perspective of model structure, the introduction of IPSO and ensemble learning theory can effectively reduce the number of neurons in the hidden layer and simplify the model. From the point of identification time, although the IPSO algorithm and the introduction of ensemble learning theory simplify the complexity of the model, the iterative optimization process of IPSO and the serial training process of the weak learner by AdaBoost. RT both lead to the increase in identification time.

To further assess the performance of the model, three well-known measures are taken [6164]. The description of these measures is as follows: (a)Correlation Coefficient (). is a statistical measure of the strength of a linear relationship between two variables. Its values can range . shows a perfect positive correlation, or a direct relationship. describes a perfect negative, or inverse, correlation, with values in one series rising as those in the other decline and vice versa. Moreover, when , there is no possible agreement between experimental results and numerical ones(b)Mean Square Error (MSE). As a common model performance evaluation function, MSE can reflect the difference between the observed and the predicted ones. The smaller the order of magnitude of MSE, the higher the identification accuracy(c)Mean Absolute Percentage Error (MAPE). MAPE is a relative error measure that uses absolute values to avoid positive and negative errors canceling out. For each data series, MAPE value is a positive ratio of error value (difference between predicted output and observed one) to observed value

The specific expression of the above-mentioned measures can be described as where represents the number of predicted samples, represents the th actual observation value, represents the th model prediction value, is the mean of , and is the mean of .

The statistical results are collected in Table 5. Under the combined action of ensemble theory and structural parameters (optimized by IPSO), IPSO-ELM-AdaBoost has the best generalization performance among the four machine learning algorithms. The structural parameters of ELM optimized by IPSO can make the input weight and threshold contain more information of input samples, which can effectively improve the identification accuracy. AdaBoost can also improve the identification accuracy by comprehensively considering the learning results of multiple weak learners. The results of IPSO-ELM and ELM-AdaBoost show that ensemble theory improves more on the performance of the algorithm. This is because structural parameters (optimized by IPSO) contain more sample information, but on the other hand, they also increase the risk of overfitting.

4.3. Simulation Test 2

In order to further verify the excellent performance of the proposed algorithm in projectile aerodynamic parameter identification, IPSO-ELM-AdaBoost is compared with the traditional projectile parameter identification algorithm (LSM, ML, and UKF). Figure 13 shows the trajectory reconstruction results of the four algorithms. The trajectory reconstructed by IPSO-ELM-AdaBoost has the highest fitting degree with the actual trajectory. Among the three traditional algorithms, the reconstructed result of UKF is the closest to the actual trajectory. The fitting degree of the four algorithms from low to high is as follows: LSM, ML, UKF, and IPSO-ELM-AdaBoost. In practical engineering, more attention is paid to the landing point and side deflection of the projectile, so the landing point () and side deflection () of the four algorithms are also presented in Table 6. Compared with the traditional algorithm, the proposed algorithm has the highest accuracy. The drop point error of IPSO-ELM-AdaBoost is 6 m, and the sideway error is 2 m. In terms of the drop point error, the drop point error of LSM and ML is more than 50 meters, which can no longer meet the need of actual engineering. From the perspective of side deviation error, UKF and IPSO-ELM-AdaBoost have good performance, and their side deviation error is less than 10 m.

5. Conclusion

Accurate aerodynamic parameters can effectively improve the firing accuracy of the projectile. Affected by uncertain factors such as actual combat environment and external meteorological conditions, the actual ballistic trajectory data is characterized by strong nonlinearity, time-dependent nature, and susceptibility to random noise. Take a single ELM for aerodynamic parameters, and the identification result curve oscillates in transonic region. To accurately obtain aerodynamic parameters of projectile, a new aerodynamic parameter identification model based on improved ELM and ensemble theory is constructed. Under standard weather conditions, the model is trained by trajectory data, and the mapping relationship between trajectory data and aerodynamic parameters is established. The simulation results show that the generalization ability of the model can be effectively improved by optimizing the structural parameters of ELM with IPSO and integrating several weak learners into one strong learner. Although the identification time is increased because of the iterative optimization of IPSO and serial training of AdaBoost. RT, the identification time of IPSO-ELM-AdaBoost is still the same order as that of ELM. The proposed IPSO-ELM-AdaBoost has excellent performance in the aerodynamic parameter identification of projectile and can be extended to the other aircraft.

Data Availability

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.