Next Article in Journal
A Barcoded Polymer-Based Cross-Reactive Spectroscopic Sensor Array for Organic Volatiles
Next Article in Special Issue
Real-Time Maritime Traffic Anomaly Detection Based on Sensors and History Data Embedding
Previous Article in Journal
Intelligent Framework for Learning Physics with Aikido (Martial Art) and Registered Sensors
Previous Article in Special Issue
Leveraging Visual Place Recognition to Improve Indoor Positioning with Limited Availability of WiFi Scans
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hybrid IRBM-BPNN Approach for Error Parameter Estimation of SINS on Aircraft

Xi’an Research Institute of High Technology, Xi’an 710025, China
*
Author to whom correspondence should be addressed.
Submission received: 9 July 2019 / Revised: 18 August 2019 / Accepted: 22 August 2019 / Published: 24 August 2019
(This article belongs to the Collection Positioning and Navigation)

Abstract

:
To realize the error parameter estimation of strap-down inertial navigation system (SINS) and improve the navigation accuracy for aircraft, a hybrid improved restricted Boltzmann machine BP neural network (IRBM-BPNN) approach, which combines restricted Boltzmann machine (RBM) and BP neural network (BPNN), is proposed to forecast the inertial measurement unit (IMU) instrument errors and initial alignment errors of SINS. Firstly, the error generation mechanism of SINS is analyzed, and initial alignment error model and IMU instrument error model are established. Secondly, an unsupervised RBM method is introduced to initialize BPNN to improve the forecast performance of the neural network. The RBM-BPNN model is constructed through the information fusion of SINS/GPS/CNS integrated navigation system by using the sum of position deviation, the sum of velocity deviation and the sum of attitude deviation as the inputs and by using the error parameters of SINS as the outputs. The RBM-BPNN structure is improved to enhance its forecast accuracy, and the pulse signal is increased as the input of the neural network. Finally, we conduct simulation experiments to forecast and compensate the error parameters of the proposed IRBM-BPNN method. Simulation results show that the artificial neural network method is feasible and effective in forecasting SINS error parameters, and the forecast accuracy of SINS error parameters can be effectively improved by combining RBM and BPNN methods and improving the neural network structure. The proposed IRBM-BPNN method has the optimal forecast accuracy of SINS error parameters and navigation accuracy of aircraft compared with the radial basis function neural network method and BPNN method.

1. Introduction

Strap-down inertial navigation system (SINS) has the advantages of high reliability, continuous outputs, independence and strong anti-interference capability and has been widely used in the field of aircraft navigation. However, the gyroscopes and accelerometers of inertial measurement unit (IMU) inevitably exhibit drift and bias caused by manufacturing process and other reasons. Measurement errors accumulate with time, thereby affecting navigation accuracy. IMU instrument errors seriously affect the stability and accuracy of aircraft flights, especially for high-speed, long-range and long-endurance ones, such as ballistic missiles, hypersonic vehicles.
An aircraft requires initial alignment before launching to determine the attitude matrix, and the initial alignment also has a huge impact on the navigation accuracy of the aircraft. In practice, initial alignment errors are inevitable and should be considered. Overall, the error parameters of SINS, including IMU instrument errors and initial alignment errors, should be compensated and corrected to improve the navigation accuracy and environmental adaptability of aircraft. At present, the estimation and correction methods of SINS error parameters include filtering algorithms, intelligent optimization algorithms and artificial neural network methods.
Filtering algorithms including Kalman filtering (KF) and its variants are widely used in the field of integrated navigation. Zhang et al. [1] proposed an estimation method for spacecraft attitude and position based on cubature KF (CKF), and the simulation results showed that the proposed method had higher attitude and position estimation accuracy than extended KF (EKF) algorithm. Zhao et al. [2] used an extended H filtering algorithm to realize the information fusion of INS, and the simulation results showed that the method could improve the navigation accuracy of INS for aerospace vehicles. Ding et al. [3] proposed an adaptive KF algorithm, that used an adaptive process noise scaling method based on covariance to improve the filter performance, and the simulation results showed that the method was highly reliable and robust for the error estimation of navigation parameters of GPS/INS integrated navigation system. However, filtering algorithms have complicated derivation of state equation and long filtering time, and KF algorithms are required to determine the precise priori knowledge of system models and noise characteristics.
The error propagation of SINS is a complex nonlinear system problem. Intelligent optimization algorithms have certain advantages in solving the error estimation problem of SINS because they are not constrained by search space limitations; does not require to have assumptions, such as continuous and derivative existence; and possess hidden parallelism, huge search space, simple structure and easy to realize. Zhu et al. [4] proposed a particle swarm optimization (PSO) algorithm to solve the parameter selection problem of compass alignment circuit of SINS, and the experimental results showed that the proposed PSO algorithm was effective and feasible for compass alignment of SINS. Guo et al. [5] used an online identification method of initial alignment errors based on adaptive PSO to realize the high-precision initial alignment of SINS for ballistic missile, and the simulation results showed that the adaptive PSO had a fast convergence velocity and high convergence accuracy, that can effectively improve the navigation accuracy of SINS. He et al. [6] employed a genetic algorithm to optimize the compass alignment parameters for SINS, and the experimental results showed that the method could effectively adapt to the compass alignment problem under the large misalignment angle. The above studies have used intelligent optimization algorithms to the optimize the initial alignment errors of INS, without considering the IMU instrument error estimation of SINS, and the initial alignment accuracy should be further improved. Few studies have been conducted on the parameter optimization of IMU instrument errors using intelligent optimization algorithms. Dai et al. [7] proposed an accelerometer error calibration method based on PSO algorithm to reduce the dependence of IMU calibration on high-precision turntable, and the method could achieve rapid accelerometer error calibration. However, the proposed method belonged to static calibration based on turntable and only calibrated the accelerometer error parameters. Of course, intelligent optimization algorithms have disadvantages of huge amount of calculation and relatively poor real-time performance because of the need for the calculation of fitness of each population.
Artificial neural network methods have good real-time calculation and strong fitting ability, without identifying the internal structure of the system. With the continuous development of artificial intelligence technology, artificial neural network methods have played an increasingly important role in the field of inertial navigation [8]. Hasan and Saadeddin et al. [9,10] adopted the adaptive neuro-fuzzy inference system (ANFIS) based on current and previous INS data to forecast the INS errors when GPS signals are interrupted, and the simulation results showed that ANFIS had a good estimation effect. Ning et al. [11] used radial basis function neural network (RBFNN) and enhanced adaptive robust KF method for GNSS/INS integrated navigation system in complex urban suburbs, and field tests showed that the method could improve the position navigation accuracy when the GNSS signal was interrupted. Yao et al. [12] applied an improved multi-layer perceptual network (MLP) to estimate the pseudo GPS position for improving the performance of GPS/INS integrated navigation system during GPS interruption, and the method had better forecast performance than traditional artificial neural network.
However, the network initial values of the above neural network forecast methods are randomly generated, and the neural network training effect is sensitive to the network initial values. A neural network easily falls into local optimal when the selected initial values are not properly determined. In addition, the forecast effect of a single network model is limited. Scholars at home and abroad have combined and improved different models to further improve the navigation accuracy of SINS and overcome the shortcomings of a single neural network model. Wang et al. [13] proposed a method to optimize the structure of BP neural network (BPNN) based on PSO, and used the optimized weight and threshold for the initialization of the neural network. The INS position error of surface ship was forecasted by using BP neural network model when the star sensor of INS/CNS integrated navigation system did not work properly. Fan et al. [14] constructed a hybrid algorithm based on genetic algorithm and BP neural network to estimate the navigation state parameter errors of GPS/DRS simple integrated navigation system. The disadvantages of local optimal and slow convergence speed of traditional algorithms were overcome, and the method had higher precision and stability compared with conventional KF algorithms. Sharaf et al. [15] used wavelet transformation to process the data noise of INS/GPS integrated navigation system and utilized RBFNN to forecast the INS position errors.
With the improvement of neural network, the above research methods have made the network initial values more reasonable and improved the forecast effect of error parameters and navigation accuracy of INS. However, in the navigation field, few studies have been conducted on the optimization of initial weight and threshold of neural network through unsupervised machine learning. The forecast performance of network models can be improved when the initial values of neural network are initialized through machine learning.
Boltzmann machine (BM) is a type of symmetrically connected and random neural network model with no self-feedback [16]. The output of BM neuron is only 0 or 1, where 1 indicates the active state and 0 indicates the inactive state. The state value can be obtained through probability statistics [17]. The BM model has powerful unsupervised learning ability to learn and extract complex information from data. However, due to the existence of full connections between each neuron within the visible and hidden layer, BM also has obvious disadvantages, such as complex structure, long network training time and low efficiency. The expected sample distribution is difficult to obtain. On this basis, Smolensky [18] proposed a Restricted Boltzmann machine (RBM) model to solve the problem on the basis of the BM model.
RBM is a bidirectional probability graph model based on statistical laws. RBM has unsupervised learning ability, easy reasoning, parallelism and fault tolerance because of its excellent characteristics and special structure. RBM is one of the basic structural units that constitute a deep neural network model, and multiple RBMs stacking can be constructed to form a deep belief neural network (DBN) [19]. RBM has good forecast effect and strong feature information extraction capability, and has been successfully utilized to solve regression, classification, image feature extraction, forecast and other problems [20,21,22,23].
A hybrid improved restricted Boltzmann machine BP neural network (IRBM-BPNN) approach combining RBM and BPNN is proposed to forecast and estimate the IMU instrument errors and initial alignment errors of SINS based on the above research and analyses and relevant research achievements. Firstly, the overall scheme of SINS error parameter estimation is designed, and IMU instrument error model and initial alignment error model are established. Then, the RBM-BPNN forecast model is constructed based on SINS/GPS/CNS integrated navigation system measurement information by using the sum of position deviation, the sum of velocity deviation and the sum of attitude deviation as the network inputs and SINS error parameters as the network outputs. The model firstly uses an unsupervised learning method to train the RBM, assigns the weight and bias of the trained RBM to the BPNN for completing the network initialization, and uses a BP algorithm to conduct supervised training on the entire neural network for improving the network forecast performance. The structure of RBM-BPNN is improved to enhance the forecast accuracy of SINS error parameters. The pulse signal is increased as the input of neural network on the basis of the original navigation deviation, and SINS error parameters are forecasted and estimated by the IRBM-BPNN. Finally, the reliability and effectiveness of IRBM-BPNN are verified through simulation experiments. The main innovative contributions of this paper are summarized as follows:
(1)
Traditional neural network methods for integrated navigation system have focused on the estimation and compensation of navigation state parameter error, where the estimation of IMU instrument errors has been rarely reported. What’s more, few studies have been conducted on parameter estimation and compensation of SINS that consider the IMU instrument errors and initial alignment errors. In this study, an artificial neural network method is used to simultaneously estimate and forecast the IMU instrument errors and initial alignment errors, and the research content is in-depth and challenging;
(2)
For the design of neural network structure, previous studies have mainly regarded the navigation state parameters or navigation deviations of integrated navigation system as the inputs or outputs of the neural network and have not utilized the pulse information from IMU. In this paper, the network structure is improved, the gyroscope pulse signal and accelerometer pulse signal are taken as the inputs of the neural network, and the data are more fully utilized, which are conducive to the estimation of SINS error parameters;
(3)
In this paper, an RBM model is used to initialize the BPNN for obtaining good initial network values to avoid the neural network falling into local optimal. An RBM-BPNN forecast model is constructed by combining RBM and BPNN networks to improve the forecast accuracy of SINS error parameters.
(4)
Unsupervised feature learning mechanism is introduced in the field of integrated navigation and SINS error parameters estimation, which enhances the processing and adaptability of the neural network model to the measurement data of integrated navigation system and provides new insight into the error parameter estimation and navigation accuracy improvement of SINS.
The rest structure of this paper is organized as follows. In Section 2, the overall scheme design of SINS error parameter estimation is introduced. The initial alignment error model and IMU instrument error model of SINS are established in Section 3. Section 4 constructs the RBM-BPNN forecast model. In Section 5, the structure of RBM-BPNN is improved, and simulation and comparison experiments are provided to verify the feasibility of the neural network method. Section 6 summarizes the conclusions, and expounds the future research direction of SINS error estimation.

2. Overall Error Estimation Scheme Design

At present, multiple source information fusion is an important means to improve navigation accuracy. The overall error estimation scheme of the RBM-BPNN method is designed based on SINS/GPS/CNS integrated navigation system information to effectively estimate and identify SINS error parameters of aircraft, as shown in Figure 1. Where N a and N b represent the accelerometer pulse output and the gyroscope pulse output, respectively [24]; f b and ω b represent the specific force and the angle velocity after error compensation, respectively; Q 0 represents the initial value of quaternion; C i b represents the attitude matrix between the inertial coordinate and body coordinate.
As shown in Figure 1, the movement of an aircraft generates signal excitation to the IMU, in which the gyroscope and accelerometer output pulse signal N b and N a , respectively. The pulse signals are then input into the SINS error model. The SINS error model includes the initial alignment error model, gyroscope error model and accelerometer error model. The specific force f b and angular velocity ω b are calculated and output by the IMU error compensation model. The initial value of quaternion Q 0 is output by the initial alignment error model. SINS navigation calculation module is conducted to calculate the SINS navigation state parameters, including attitude, position and velocity, based on the above output information. In addition, auxiliary navigation system GPS can measure the position and velocity values of the aircraft, and CNS can measure the attitude values of the aircraft.
Then, multiple source information is fused based on the SINS/GPS/CNS integrated navigation system measurements. The specific fusion strategy in this paper can be described as follows. The sum of position deviation and the sum of velocity deviation between SINS and GPS are calculated based on the position output and velocity output of SINS/GPS integrated navigation system. The sum of attitude deviation between SINS and CNS is calculated based on the attitude output of the SINS/CNS integrated navigation system.
Finally, the sum of position deviation, the sum of velocity deviation and the sum of attitude deviation after information fusion are input into the RBM-BPNN forecast model, and then the initial alignment errors and IMU instrument errors of SINS are forecasted and estimated by the RBM-BPNN method. Once the estimation of error parameters is complete, the switch in Figure 1 is disconnected. Then the initial alignment error parameter forecast values are fed back to the initial alignment error model to update the initial value of quaternion Q 0 ; the IMU instrument error parameter forecast values are fed back to the IMU error compensation model to realize the error compensation and correction. The estimation of SINS error parameters is finally completed to meet the navigation accuracy requirements of aircraft.

3. SINS Error Model Establishment

According to the principle of SINS measurement, there are two main sources of navigation error.
The first source is the initial measurement reference error of SINS. The calculation of navigation state parameters of SINS is conducted under the given reference conditions, where the attitude matrix established by SINS will have an error when the given reference has an error, which will cause navigation error in measurement and navigation calculation of SINS. Therefore, the effect of initial alignment errors must be considered.
The second source is the measurement model error of SINS. Theoretically, the measurement error of SINS should be small and the navigation accuracy of aircraft should be high when the IMU instrument error coefficient calibrated at the ground stage is accurate after the real-time error compensation by the IMU error compensation model. However, because the bias and drift of IMU are variable in different start each time, the calibration values of the IMU instrument error coefficient are not applicable. The measurement error of SINS accumulates over time, which affects the navigation accuracy of aircraft. Therefore, the effect of IMU instrument errors must also be considered.

3.1. Initial Alignment Error Model of SINS

Ideally, a vertically launched aircraft has no initial attitude errors before take-off. However, initial alignment errors, including initial pitch angle error Δ φ 0 , initial yaw angle error ψ 0 and initial roll angle error γ 0 , exist during the actual launch because of installation error, aiming error and other effects. The initial alignment errors affect the establishment of the attitude matrix C i b , and C i b can be described by quaternion, which is shown as follows:
C i b = [ q 1 2 + q 0 2 q 2 2 q 3 2 2 ( q 1 q 2 q 0 q 3 ) 2 ( q 0 q 2 + q 1 q 3 ) 2 ( q 1 q 2 + q 0 q 3 ) q 0 2 + q 2 2 q 1 2 q 3 2 2 ( q 2 q 3 q 0 q 1 ) 2 ( q 1 q 3 q 0 q 2 ) 2 ( q 0 q 1 + q 2 q 3 ) q 0 2 + q 3 2 q 1 2 q 2 2 ]
where q 0 , q 1 , q 2 and q 3 represent the quaternion, and the initial value of quaternion Q 0 = ( q 0 , q 1 , q 2 , q 3 ) is calculated as:
{ q 0 = q 00 γ 0 2 q 20 q 1 = q 10 + γ 0 2 q 30 q 2 = q 20 + γ 0 2 q 00 q 3 = q 30 γ 0 2 q 10
where
{ q 00 = 2 2 ( 1 Δ φ 0 2 ) q 10 = ψ 0 2 q 00 q 20 = ψ 0 2 q 30 q 30 = 2 2 ( 1 + Δ φ 0 2 )
As shown in Equations (1)–(3), it can be seen that the initial alignment errors determine the initial value of quaternion Q 0 , and the initial value of quaternion determines the calculation result of the attitude matrix C i b . Therefore, as long as the initial alignment error parameter values are estimated accurately, the attitude matrix can be corrected by the above initial alignment error model, thus the navigation accuracy of aircraft can be improved.

3.2. IMU Instrument Error Model of SINS

As the core component of SINS, IMU is mainly composed of gyroscope and accelerometer. The drift term of gyroscope error coefficient mainly includes the zero-order error coefficient and first-order error coefficient, and the installation error coefficient of gyroscope remains unchanged after the IMU is manufactured and stored in a stable period. The zero-order error coefficient and the first-order error coefficient of accelerometer will change, except that the installation error coefficient of accelerometer basically remains same. Therefore, the error model of the accelerometer and gyroscope can be obtained as follows [25].
The accelerometer error model can be formulated as:
{ δ f b x = Δ K 0 x + Δ K 1 x f b x + ε a x δ f b y = Δ K 0 y + Δ K 1 y f b y + ε a y δ f b z = Δ K 0 z + Δ K 1 z f b z + ε a z
where δ f b x , δ f b y and δ f b z represent the apparent acceleration errors of aircraft in body coordinate; Δ K 0 x , Δ K 0 y and Δ K 0 z represent the zero-order coefficient errors of accelerometer; Δ K 1 x , Δ K 1 y and Δ K 1 z represent the first-order coefficient errors of accelerometer; f b x , f b y and f b z represent the component of the specific forces in body coordinate; ε a x , ε a y and ε a z represent the measurement white noises of accelerometer.
The gyroscope error model can be formulated as:
{ δ ω b x = Δ D 0 x + Δ D 1 x f b x + ε g x δ ω b y = Δ D 0 y + Δ D 1 y f b x + ε g y δ ω b z = Δ D 0 z + Δ D 1 z f b x + ε g z
where δ ω b x , δ ω b y and δ ω b z represent the angle velocity errors of aircraft in body coordinate; Δ D 0 x , Δ D 0 y and Δ D 0 z represent the zero-order drift coefficient errors of gyroscope; Δ D 1 x , Δ D 1 y and Δ D 1 z represent the first-order coefficient errors of gyroscope; ε g x , ε g y and ε g z represent the measurement white noises of gyroscope.
In addition, the establishment and description of the IMU error compensation model and the SINS navigation calculation model can be referred in [5], which will not be repeated here. The navigation state parameters of SINS, such as attitude, position and velocity information, can be calculated based on the IMU error compensation model and SINS navigation calculation model.

4. RBM-BPNN Method Design for SINS Error Estimation

4.1. Introduction of Basic Principles of RBM

As an undirected graph model, RBM is a special topology from BM based on energy functions. An RBM consists of a visible layer and a hidden layer. The visible layer represents the input data, and the hidden layer represents the feature information extracted from the input data. Unlike the BM model, there is no connection between the neurons within the same layers for the RBM model. Consequently, RBM has the advantages of simple structure and easy learning. The RBM network structure is shown in Figure 2.
Where n represents the number of neurons in the visible layer, and v = { v 1 , v 2 , , v n } represents the state of the visible units; m represents the number of neurons in the hidden layer, and h = { h 1 , h 2 , , h m } represents the state of the hidden units; W represents the connection weight between the visible layer and the hidden layer; a = { a 1 , a 2 , , a n } and b = { b 1 , b 2 , , b m } represent the biases of the visible layer and the hidden layer, respectively. For a given state ( v , h ) , the energy function of RBM network model is defined as follows:
E ( v , h ) = a T v b T h h T W v
When RBM network parameter θ = ( a , b , W ) is determined, joint probability distribution P ( v , h ) can be obtained based on the energy function given by Equation (6), that is:
{ P ( v , h ) = 1 Z e E ( v , h ) Z = v , h e E ( v , h )
where Z represents the normalization factor, also known as the partition function. For practical engineering problems, visible layer v represents the actual measurement data, and its edge probability distribution P ( v ) can be calculated based on joint probability distribution P ( v , h ) , that is:
P ( v ) = h P ( v , h ) = 1 Z h e E ( v , h )
The conditional probability distributions of visible layer and hidden layer neurons should be determined to estimate parameter θ of the RBM model. On the basis of related theory of conditional probability distribution and reference [16], when the visible layer is known, condition probability distribution P ( h j = 1 | v ) of hidden layer neuron h j = 1 can be expressed as:
P ( h j = 1 | v ) = 1 1 + e ( b j + i v i W i j ) = f ( b j + i v i W i j )
Similarly, when the hidden layer is known, condition probability distribution P ( v i = 1 | h ) of visible layer neuron v i = 1 can be expressed as:
P ( v i = 1 | h ) = 1 1 + e ( a i + j W i j h j ) = f ( a i + j W i j h j )
As shown in Equation (9) and (10), f ( x ) = 1 1 + e x is the sigmoid activation function.
To determine the RBM model, parameter θ = ( a , b , W ) should be calculated to ensure that the probability distribution represented by RBM is consistent with the sample data distribution. The edge probability distribution P ( v ) is shown in Equation (8), and its log-likelihood function l ( θ ) can be described as:
l ( θ ) = log ( P ( v ) )
By maximizing the log-likelihood function l ( θ ) , we can get the parameter estimation θ * of RBM model, that is:
θ * = arg max θ l ( θ ) = arg max θ log ( P ( v ) ) = arg max θ log ( 1 Z h e E ( v , h ) )
The partial derivative of the log-likelihood function l ( θ ) is taken with respect to parameters ( a , b , W ) , and we obtain:
{ l ( θ ) a i = v i d a t a v i mod e l l ( θ ) b j = h j d a t a h j mod e l l ( θ ) W i j = v i h j d a t a v i h j mod e l
where P represents the mathematical expectation of the distribution P , and ‘data’ and ‘model’ represent the probability distribution P ( h | v ) and P ( v , h ) respectively. From Equation (13) we can see that the calculation of the partial derivative needs to obtain a joint probability distribution P ( v , h ) , and it can be seen from Equation (7) that the acquisition of P ( v , h ) requires to calculated the normalization factor Z , while the calculation complexity of Z value is large. Aimed at the problem, Gibbs sampling is usually used to obtain an approximation of the probability distribution P ( v , h ) , but Gibbs sampling requires a large number of sampling steps, which results in low training efficiency of RBM.
To overcome this difficulty, Hintion proposed a contrast divergence (CD) algorithm to effectively solve the problem of RBM learning efficiency [26]. The expectations in Equation (13) are approximately estimated by using the CD algorithm, and the calculation results of the partial derivative are expressed as follows:
{ Δ a i = v i d a t a v i r e c o n Δ b j = h j d a t a h j r e c o n Δ W i j = v i h j d a t a v i h j r e c o n
where Δ a i , Δ b j and Δ W i j are the approximate values of partial derivatives l ( θ ) a i , l ( θ ) b j and l ( θ ) W i j , respectively, and ‘recon’ represents the reconstructed RBM model distribution. Finally, the estimated values of parameters ( a i , b j , W i j ) of the RBM model can be obtained as follows.
{ a i = a i + η Δ a i b j = b j + η Δ b j W i j = W i j + η Δ W i j
where η represents learning rate of RBM model.

4.2. RBM-BPNN Forecast Model Construction

RBM and BPNN theories state that an output layer is added to the two-layer RBM network, and RBM is expanded into a three-layer NN, which is the RBM-BPNN model.
As shown in Figure 3, the input and output data of the neural network should be determined before determining the RBM-BPNN structure. As previously mentioned in Section 2, the inputs of RBM-BPNN structure can be obtained through the fusion of multiple information sources of the SINS/GPS/CNS integrated navigation system, that is, the sum of position deviation, the sum of velocity deviation and the sum of attitude deviation.
The sum of position deviation and the sum of velocity deviation can be calculated based on the position and velocity output values of the SINS/GPS integrated navigation system as follows:
{ δ x = t = 0 T k ( x S I N S ( t ) x G P S ( t ) ) δ y = t = 0 T k ( y S I N S ( t ) y G P S ( t ) ) δ z = t = 0 T k ( z S I N S ( t ) z G P S ( t ) )
where δ x , δ y and δ z are the sum of position deviations between SINS and GPS from time 0 to time T k ; x G P S ( t ) , y G P S ( t ) and z G P S ( t ) are the position parameters measured by GPS at t moment; x S I N S ( t ) , y S I N S ( t ) and z S I N S ( t ) are the position parameters outputted by SINS at t moment. T k is the flight time of aircraft; time step t is 0.1 s.
The sum of velocity deviations is calculated as follows:
{ δ v x = t = 0 T k ( v x S I N S ( t ) v x G P S ( t ) ) δ v y = t = 0 T k ( v y S I N S ( t ) v y G P S ( t ) ) δ v z = t = 0 T k ( v z S I N S ( t ) v z G P S ( t ) )
where δ v x , δ v y and δ v z are the sum of velocity deviations between SINS and GPS from time 0 to time T k ; v x G P S ( t ) , v y G P S ( t ) and v z G P S ( t ) are the velocity parameters measured by GPS at t moment; v x I N S ( t ) , v y S I N S ( t ) and v z S I N S ( t ) are the velocity parameters outputted by SINS at t moment.
The sum of attitude deviation can be calculated based on the attitude output value of the SINS/CNS integrated navigation system as follows:
{ δ φ = t = 0 T k ( φ S I N S ( t ) φ C N S ( t ) ) δ ψ = t = 0 T k ( ψ S I N S ( t ) ψ C N S ( t ) ) δ γ = t = 0 T k ( γ S I N S ( t ) γ C N S ( t ) )
where δ φ , δ ψ and δ γ are the sum of attitude deviations between SINS and CNS from time 0 to time T k ; φ C N S ( t ) , ψ C N S ( t ) and γ C N S ( t ) are the attitude parameters measured by CNS at t moment; φ S I N S ( t ) , ψ S I N S ( t ) and γ S I N S ( t ) are the attitude parameters outputted by SINS at t moment.
As shown in Equation (1) to (5), the initial alignment errors ( Δ φ 0 , ψ 0 , γ 0 ) and gyroscope instrument errors ( Δ D 0 x , Δ D 0 y , Δ D 0 z , Δ D 1 x , Δ D 1 y , Δ D 1 z ) and accelerometer instrument errors ( Δ K 0 x , Δ K 0 y , Δ K 0 z , Δ K 1 x , Δ K 1 y , Δ K 1 z ) are the error parameters that need to be forecasted and estimated in this paper, and also the outputs of RBM-BPNN forecast model. Therefore, RBM-BPNN model is a nine-input and fifteen-output neural network structure, the mapping relationship between the inputs and outputs can be expressed as follows:
( Δ φ 0 , ψ 0 , γ 0 , Δ D 0 x , Δ D 0 y , Δ D 0 z , Δ D 1 x , Δ D 1 y , Δ D 1 z , Δ K 0 x , Δ K 0 y , Δ K 0 z , Δ K 1 x , Δ K 1 y , Δ K 1 z ) = f ( δ x , δ y , δ z , δ v x , δ v y , δ v z , δ φ , δ ψ , δ γ )
The mapping relationship of Equation (19) through analytic calculation is difficult to obtain because the error propagation of SINS is a complex nonlinear process, and the resolution between IMU instrument errors and navigation state parameter deviation is difficult to derive. The RBM-BPNN method has powerful mapping ability for nonlinear systems, which does not require to master and derive the error propagation within SINS, and can directly fit the physical model based on input and output data. The RBM-BPNN structure can be determined by training the sample data, and the mapping relationship between the inputs and outputs is directly established, which is beneficial to solving the problem in this paper.
The specific flow of the error parameters forecast of SINS based on the RBM-BPNN method is described as follows:
(1)
Data preprocess. The sample data are collected based on SINS/GPS/CNS integrated navigation system measurement and SINS error parameter information. Then, the sample data are normalized and mapped to the interval of [0, 1], which is expressed as:
x ˜ = ( x ˜ x ˜ min ) / ( x ˜ max x ˜ min )
where x ˜ max and x ˜ min are the maximum value and minimum value of the sample data, respectively;
(2)
The training sample is input into the RBM model, the RBM structure is trained by using a greedy learning algorithm without supervision, and the deep feature information of the input data is extracted. Then, the weight and bias of the RBM structure are determined;
(3)
The BPNN is initialized based on the RBM training results. The weight and threshold between input layer and hidden layer of BPNN are replaced by the weight and bias of RBM;
(4)
The weight and threshold value of the entire network of RBM-BPNN are fine-tuned by using a BP algorithm, and the mapping relationship between the input layer and output layer is obtained, that is, the RBM-BPNN structure;
(5)
The trained RBM-BPNN model is evaluated by testing samples, and the generalization ability and error parameter forecast effect of the model are verified.
We adopt two criteria to evaluate the forecast performance of SINS error parameters based on the RBM-BPNN model, including mean absolute percentage error (MAPE) and mean absolute error (MAE).
The calculation equations of MAPE and MAE are shown in Equation (21) and (22), respectively.
M A P E = 1 n s i = 1 n s R i = 1 n s i = 1 n s | x i x ˜ i x ˜ i | × 100 %
M A E = 1 n s i = 1 n s | x i x ˜ i |
where n s represents the total number of samples; x i represents the error parameter forecast value of ith sample; x ˜ i represents the error parameter actual value of ith sample; R i represents the relative error of the error parameter of ith sample. It is noted that the relative error is set to 100% when the relative error of the error parameter is more than 100%. The smaller the MAPE and MAE values are, the higher the estimation accuracy, and the better the performance of the RBM-BPNN model will be.

5. Simulation Experiments and Result Analysis

5.1. Simulation Condition Settings

The simulation conditions of SINS/GPS/CNS integrated navigation system and initial alignment error parameters are set as follows:
(1)
GPS and CNS navigation errors. Position measurement error and altitude measurement error of GPS are 10.0 m, velocity measurement error of GPS is 0.1 m/s, and attitude measurement error of CNS is 10″ [2,27,28];
(2)
The initial launch parameters of aircraft. The coordinate of the launch point is (118°E, 32°N, 0.0 m), and the launch azimuth is 90° [28];
(3)
The initial alignment error parameters of SINS. The standard deviation of the initial pitch angle error is 60″, the standard deviation of the initial yaw angle error is 60″, and the standard deviation of the initial roll angle error is 90″ [29];
(4)
IMU instrument error parameter settings are shown in Table 1 [24].
This study aims to take high-speed and high-dynamic aircrafts as the research objects, such as ballistic missiles, hypersonic vehicles, and realize the error parameters estimation of SINS during the ascent phase flight of aircraft. Suppose that the aircraft employs a three-stage rocket to achieve flight, aircraft flight time of ascent phase is T k = 160 s, and SINS navigation cycle is 0.02 s. To generate signal excitation to SINS and realize the effective estimation of error parameters, the standard flight path of ascent phase of aircraft is designed in geocentric coordinate, as shown in Figure 4. Figure 5 shows the velocity curve of the aircraft in inertial coordinate.

5.2. Sample Generation and RBM-BPNN Structure Training

Training sample data should be collected to the train RBM-BPNN structure. The initial alignment error parameters of SINS and IMU instrument error parameters set in Section 5.1 are taken as 1 standard deviation σ , and 5000 groups of SINS error parameters are randomly generated in the form of normal distribution N ( 0 , σ ) , which are taken as the training sample outputs of the neural network. We can obtain 5000 groups of navigation state parameters with navigation error, including position, velocity and attitude, by injecting 5000 groups of error parameters into the flight simulation software of the aircraft. The IMU error compensation model and SINS navigation calculation model of the software are consistent with the onboard computer of the aircraft to ensure consistency and accuracy of the model. The sum of position deviation ( δ x , δ y , δ z ) , the sum of velocity deviation ( δ v x , δ v y , δ v z ) and the sum of attitude deviation ( δ φ , δ ψ , δ γ ) can be calculated based on SINS/GPS/CNS integrated navigation system measurement information and Equation (16)–(18). Thus, 5000 groups of navigation state parameter errors can be obtained, which are taken as the training sample inputs of the neural network.
According to the sample inputs and sample outputs and Equation (19), it can be seen that the number of input layer nodes is n = 9, and the number of the output layer nodes is q = 15. The number of hidden layer nodes m can be determined by using equation m = n + q + a , where a is a constant in the interval [1,10], and m is chosen as 14 in this paper. The neural network activation function is chosen as the sigmoid function, RBM learning rate is η = 0.0001 , and BPNN learning rate is η = 0.01 . RBM and BPNN error training times are set to 1000. RBM-BPNN structure can be trained based on the above information, and the training flow chart of RBM-BPNN is shown in Figure 6.

5.3. Error Parameter Forecast Results by RBM-BPNN

A total of 500 groups of data are randomly generated similar to Section 5.2 to verify the generalization capability of the trained RBM-BPNN structure, which are used as testing samples to evaluate the trained RBM-BPNN structure. According to simulation and calculation, the simulation results of initial alignment residual error using RBM-BPNN method is shown in Figure 7, the simulation results of gyroscope residual error using RBM-BPNN method is shown in Figure 8, and the simulation results of accelerometer residual error using RBM-BPNN method is shown in Figure 9. The residual error is the deviation between the actual value and the forecast value. Statistical results of error parameter forecast values using RBM-BPNN method are shown in Table 2.
As shown in Figure 7, Figure 8 and Figure 9 and Table 2, the RBM-BPNN method has high forecast accuracy for the initial alignment error parameters. The maximum value of initial pitch angle error and initial yaw angle error is not more than 30″, and the maximum value of initial roll angle error is not more than 50″. The MAE values of the initial alignment error are less than 20″, and the MAPE values are less than 30%. The RBM-BPNN method has poor forecast effect on gyroscope error parameters, and the MAPE values of gyroscope errors are more than 60%. For accelerometer error parameter, the RBM-BPNN method forecast effect is general, the MAPE values of accelerometer errors are generally more than 60%, but the forecast accuracy of accelerometer error parameter Δ K 1 x is relatively high, and its MAPE value is 21.05%. This condition is because the axial acceleration of the body is the largest and has a good incentive effect on the estimation of the parameter Δ K 1 x . In general, the RBM-BPNN method has certain forecast effect for error parameters of SINS, especially for the initial alignment error parameters. However, the forecast accuracy of IMU instrument error parameters is poor, the navigation accuracy requirements of high dynamic and high-speed aircraft are difficult to meet, and the forecast accuracy of SINS error parameters should be improved.

5.4. RBM-BPNN Method Improvement and Simulation

As discussed in Section 5.3, the RBM-BPNN method has poor accuracy in forecasting SINS error parameters. The RBM-BPNN structure is improved to enhance the forecast effect of error parameters.
The inputs of RBM-BPNN structure are the sum of position deviation, the sum of velocity deviation and the sum of attitude deviation, which are obtained through the information fusion of SINS/GPS/CNS integrated navigation system. It should be pointed out that the inputs of neural network structure only utilize the navigation state parameter information. The gyroscope pulse and accelerometer pulse outputs of SINS are increased to the inputs of RBM-BPNN structure by utilizing the integrated navigation system information to realize the effective estimation of error parameters. This strategy is mainly adopted because the variation of SINS error parameters will generate the navigation state parameter deviation and change the output value of pulse signals. Therefore, the navigation state parameter deviation and pulse signal are used as the inputs of RBM-BPNN to obtain sufficient excitation information, which is beneficial in forecasting SINS error parameters.
The pulse signal should be processed before utilizing it as the input of RBM-BPNN, and the processing strategy of gyroscope pulse signal is shown in Equation (23):
{ N B x 1 = t = 0 T k N b x ( t ) N b x ( t ) < 0 N B x 2 = t = 0 T k N b x ( t ) N b x ( t ) 0 N B y 1 = t = 0 T k N b y ( t ) N b y ( t ) < 0 N B y 2 = t = 0 T k N b y ( t ) N b y ( t ) 0 N B z 1 = t = 0 T k N b z ( t ) N b z ( t ) < 0 N B z 2 = t = 0 T k N b z ( t ) N b z ( t ) 0
where N b x ( t ) , N b y ( t ) and N b z ( t ) denote the gyroscope pulse outputs at t moment; N B i 1 denotes the sum of gyroscope pulse N b i ( t ) with its pulse value less than 0; N B i 2 denotes the sum of gyroscope pulse N b i ( t ) with its pulse value greater than or equal to 0; i = x , y , z denote x, y and z axes, respectively.
Similarly, the processing strategy of accelerometer pulse signal is shown in Equation (24):
{ N A x 1 = t = 0 T k N a x ( t ) N a x ( t ) < 0 N A x 2 = t = 0 T k N a x ( t ) N a x ( t ) 0 N A y 1 = t = 0 T k N a y ( t ) N a y ( t ) < 0 N A y 2 = t = 0 T k N a y ( t ) N a y ( t ) 0 N A z 1 = t = 0 T k N a z ( t ) N a z ( t ) < 0 N A z 2 = t = 0 T k N a z ( t ) N a z ( t ) 0
where N a x ( t ) , N a y ( t ) and N a z ( t ) denote the accelerometer pulse outputs at t moment; N a i 1 denotes the sum of accelerometer pulse N a i ( t ) with its pulse value less than 0; N a i 2 denotes the sum of accelerometer pulse N a i ( t ) with its pulse value greater than or equal to 0.
It should be noted that the positive and negative values of pulse are distinguished when the pulse signal is processed by using Equation (23) and (24), this strategy can prevent the pulse summation calculation from offsetting the positive and negative pulses, so as to obtain considerable useful pulse information, which is conducive in the estimation of SINS error parameters.
Therefore, the outputs of improved RBM-BPNN structure remain the same, and the inputs increase the pulse information in addition to the navigation state parameter deviations. The mapping relationship between the inputs and outputs of improved RBM-BPNN is adjusted as follows:
( Δ φ 0 , ψ 0 , γ 0 , Δ D 0 x , Δ D 0 y , Δ D 0 z , Δ D 1 x , Δ D 1 y , Δ D 1 z , Δ K 0 x , Δ K 0 y , Δ K 0 z , Δ K 1 x , Δ K 1 y , Δ K 1 z ) = f ( δ x , δ y , δ z , δ v x , δ v y , δ v z , δ φ , δ ψ , δ γ , N B x 1 , N B x 2 , N B y 1 , N B y 2 , N B z 1 , N B z 2 , N A x 1 , N A x 2 , N A y 1 , N A y 2 , N A z 1 , N A z 2 )
This method is called the IRBM -BPNN method, which uses navigation state parameter deviation and pulse signal as inputs.
Next, the error parameters of the above 500 groups of test samples are forecasted by the IRBM-BPNN method to verify the feasibility of the proposed method. After simulation experiments, the residual errors of the error parameter forecast values using RBM-BPNN method and IRBM-BPNN method are shown in Figure 10, Figure 11 and Figure 12. Table 3 shows the statistical results of error parameter forecast values using IRBM-BPNN method.
As shown in Figure 10, Figure 11 and Figure 12, the IRBM-BPNN method has smaller residual error on the error parameter forecast values of SINS, and the forecast accuracy is obviously improved compared with the RBM-BPNN method. As can be seen from Table 3, the IRBM-BPNN method has high forecast effect on the initial alignment error parameter. The maximum value of the initial alignment errors is not more than 45″, the MAE values are less than 10″, and the MAPE values are less than 20%. The forecast accuracy of gyroscope error of the IRBM-BPNN method is better than that of RBM-BPNN method, and the MAPE values of gyroscope zero-order error parameters are less than 55% by using the IRBM-BPNN method. For the accelerometer error parameters, the IRBM-BPNN method shows high forecast level, and the MAPE values of parameters Δ K 0 y , Δ K 1 y , Δ K 1 z are less than 40%. Especially, the MAPE value of parameter Δ K 1 x is less than 20%. Therefore, it is feasible to improve the forecast accuracy of SINS error parameters by improving the RBM-BPNN structure and increasing the pulse signal as the network inputs. The IRBM-BPNN method has better forecast effect compared with the RBM-BPNN method.

5.5. Design of Simulation Comparison Experiments

Comparative experiments are performed under the same simulation conditions, and the inputs and outputs of the training sample are consistent with Section 5.4 to further verify the effectiveness and feasibility of the IRBM-BPNN method in forecasting SINS error parameters. The RBFNN method and BPNN method are introduced as the comparative reference, and the radial base function of RBFNN is chosen as the Gaussian function. The statistical results of error parameter forecast values of testing samples by using the RBFNN method, BPNN method and IRBM-BPNN method are shown in Table 4. The residual error results using the three methods are shown in Figure 13, Figure 14 and Figure 15.
From Table 4 and Figure 13, Figure 14 and Figure 15, we can see that the three neural network methods have a good estimation effect on the initial alignment error, and the MAE values are less than 30″. This condition indicates that the artificial neural network method has certain advantages in forecasting SINS error parameters. For the forecast effect of IMU instrument error parameters, the BPNN method is better than the RBFNN method, whereas the IRBM-BPNN method is better than the BPNN method. The application of RBM to the initialization of BPNN can prevent the neural network from falling into local optimal, and the combination of RBM and BPNN methods can improve the forecast accuracy of SINS error parameters. At the same time, it is shown from Table 4 that the forecast effect of gyroscope first-order error parameters Δ D 1 x , Δ D 1 y , Δ D 1 z and accelerometer zero-order error parameters Δ K 0 x , Δ K 0 z by the three neural network methods are not ideal. Their MAPE values are more than 60%, and research and improvement should be further conducted.
In summary, the IRBM-BPNN method has the highest forecast accuracy of SINS error parameters compared with the RBFNN method and BPNN method, which is feasible and effective in forecasting the initial alignment errors and IMU instrument errors of SINS for aircraft.
Simulation experiments are performed again to validate the improved effect of SINS error forecast and aircraft navigation accuracy based on the IRBM-BPNN method. The initial alignment error parameters and the IMU instrument error parameters in Section 5.1 are set as the actual error parameters, and the actual error parameters are estimated under the same simulation conditions by the RBFNN method, BPNN method and IRBM-BPNN method. The forecast results of SINS error parameter by the RBFNN, BPNN and IRBM-BPNN methods through simulation are shown in Table 5.
From Table 5, we can find that RBFNN method has a good forecast effect on error parameters Δ φ 0 , γ 0 , Δ D 0 x , Δ D 0 y , Δ D 0 z , Δ D 1 y , Δ D 1 z and Δ K 0 x , in which their relative errors are less than 40%. The BPNN method has a good forecast effect on error parameters Δ φ 0 , ψ 0 , γ 0 , Δ D 0 y , Δ D 0 z , Δ K 0 y , Δ K 0 z , Δ K 1 x , Δ K 1 y and Δ K 1 z , in which their relative errors are less than 30%. The IRBM-BPNN method has a good forecast effect on error parameters Δ φ 0 , ψ 0 , γ 0 , Δ D 0 x , Δ D 0 z , Δ D 1 z , Δ K 0 y , Δ K 0 z , Δ K 1 x , Δ K 1 y and Δ K 1 z , in which their relative errors are less than 30%. Obviously, the IRBM-BPNN method has the optimal forecast accuracy of SINS error parameters.
According to the overall scheme design shown in Figure 1, the IMU instrument error parameter forecast values of the artificial neural network method are fed back to the IMU error compensation model for compensation, and the initial alignment error parameter forecast values are fed back to the initial alignment error model for compensation. After error compensation and correction, the attitude error, velocity error and position error of the aircraft using the artificial neural network method can be obtained, as shown in Figure 16, Figure 17, and Figure 18. Table 6 shows the root mean square (RMS) error statistics of navigation parameters.
As shown in Figure 16, Figure 17 and Figure 18, the attitude error, velocity error and position error are smaller than the pure SINS navigation error after the compensation and correction of the error parameters by the artificial neural network method. Therefore, the artificial neural network method can effectively improve the navigation accuracy of aircraft. As shown in Table 6, the RMS errors of attitude of the IRBM-BPNN method are less than 15″, which are smaller than the attitude errors of pure SINS navigation and the smallest among the attitude errors of the RBFNN method and BPNN method. The RMS errors of position and velocity of the IRBM-BPNN method are less than 15 m and 0.25 m/s, respectively. Obviously, the navigation errors of the IRBM-BPNN method are the smallest compared with the RBFNN method and BPNN method. Therefore, the IRBM-BPNN method has the optimal navigation accuracy.

6. Conclusions

In this paper, a hybrid IRBM-BPNN approach that combines RBM and BPNN is proposed to estimate the IMU instrument error parameters and initial alignment error parameters of SINS for aircraft. From the view of machine learning, the network initialization of BPNN is realized by RBM. The RBM-BPNN forecast model takes the navigation deviations as the inputs and SINS error parameters as the outputs based on the measurement information of SINS/GPS/CNS integrated navigation system. The RBM-BPNN structure is improved to enhance the forecast accuracy of SINS error parameters, and the error parameters are forecasted and compensated by the IRBM-BPNN method. The simulation results show that the forecast performance of the neural network can be effectively improved by combining RBM and BPNN methods and improving the RBM-BPNN network structure. The IRBM-BPNN method has the highest forecast accuracy for SINS error parameters and can effectively improve the navigation accuracy of aircraft compared with the RBFNN method and BPNN method.
This paper provides new insight into the estimation of initial alignment error parameters and IMU instrument error parameters of SINS by introducing the unsupervised machine learning mechanism. This study mainly conducts a detailed analysis and demonstrates for SINS error parameter estimation of aircraft at the theoretical and simulation levels, and the reliability and stability of the proposed method should be further verified through flight testing in the future. On the whole, with the rapid development of artificial intelligence and deep learning technology, the advantages of artificial neural network will increase, and the estimation of SINS error parameters for aircraft based on artificial neural network and deep learning will become a hot research topic in the future.

Author Contributions

W.G. and Y.X. proposed the idea of this work; D.Z. and B.L. conceived and designed the experiments; W.G. performed the experiments; L.R. reviewed and edited the manuscript. All authors read and approved this manuscript.

Funding

This work was supported in part by the National Natural Science Foundation of China (No. 61603398) and the Young Talent fund of University Association for Science and Technology in Shaanxi, China (No. 20170109).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhang, L.J.; Yang, H.B.; Lu, H.P.; Zhang, S.F.; Cai, H.; Qian, S. Cubature Kalman filtering for relative spacecraft attitude and position estimation. Acta Astronaut. 2014, 105, 254–264. [Google Scholar] [CrossRef]
  2. Zhao, H.; Xiong, Z.; Shi, L.J.; Yu, F.; Liu, J.Y. A robust filtering algorithm for integrated navigation system of aerospace vehicle in launch inertial coordinate. Aerosp. Sci. Technol. 2016, 58, 629–640. [Google Scholar] [CrossRef]
  3. Ding, W.D.; Wang, J.L.; Chris, R.; Doug, K. Improving adaptive Kalman estimation in GPS/INS integration. J. Navig. 2007, 60, 517–529. [Google Scholar] [CrossRef]
  4. Zhu, B.; Xu, J.N.; He, H.Y.; Guo, S.L. Initial alignment method of strapdown gyrocompass based on particle swarm optimization algorithm. J. Chin. Inert. Technol. 2017, 25, 47–51. [Google Scholar]
  5. Guo, W.L.; Xian, Y.; Li, B.; Ren, L.L. Initial alignment error on-line identification based on adaptive particle swarm optimization algorithm. Math. Probl. Eng. 2018, 2018, 3486492. [Google Scholar] [CrossRef]
  6. He, H.Y.; Xu, J.N.; Qin, F.J.; Li, F. Genetic algorithm based fast alignment method for strap-down inertial navigation system with large azimuth misalignment. Rev. Sci. Instrum. 2015, 86, 1930–1941. [Google Scholar] [CrossRef] [PubMed]
  7. Dai, S.W.; Wang, K.H.; Dai, H.D. A rapid calibration method for accelerometer based on PSO algorithm. Electron. Opt. Control 2014, 21, 57–60. [Google Scholar]
  8. Chiang, K.W.; Chang, H.W. Intelligent sensor positioning and orientation through constructive neural network-embedded INS/GPS integration algorithms. Sensors 2010, 10, 9252–9285. [Google Scholar] [CrossRef]
  9. Hasan, A.M.; Samsudin, K.; Ramli, A.R. GPS/INS integration based on dynamic ANFIS network. Int. J. Control Autom. 2012, 5, 1–22. [Google Scholar]
  10. Saadeddin, K.; Abdel-Hafez, M.F.; Jaradat, M.A.; Jarrah, M.A. Optimization of intelligent approach for low-cost INS/GPS navigation system. J. Intell. Robot. Syst. 2014, 73, 325–348. [Google Scholar] [CrossRef]
  11. Ning, Y.P.; Wang, J.; Han, H.Z.; Tan, X.L.; Liu, T.J. An optimal radial basis function neural network enhanced adaptive robust Kalman filter for GNSS/INS integrated systems in complex urban areas. Sensors 2018, 18, 3091. [Google Scholar] [CrossRef]
  12. Yao, Y.Q.; Xu, X.S.; Zhu, C.C.; Chan, C.Y. A hybrid fusion algorithm for GPS/INS integration during GPS outages. Measurement 2017, 103, 42–51. [Google Scholar] [CrossRef]
  13. Wang, Q.Y.; Li, Y.B.; Diao, M.; Gao, W.; Qi, Z. Performance enhancement of INS/CNS integration navigation system based on particle swarm optimization back propagation neural network. Ocean Eng. 2015, 108, 33–45. [Google Scholar] [CrossRef]
  14. Fan, Y.Z.; Hang, Y.N.; Ma, H.K. Application of the hybrid BP/GA algorithm in simple integrated navigation system. J. Beijing Univ. Aeronaut. Astronaut. 2005, 31, 535–538. [Google Scholar]
  15. Sharaf, R.; Noureldin, A.; Osman, A.; Ei-Sheimy, N. Online INS/GPS integration with a radial basis function neural network. IEEE Aerosp. Electron. Syst. Mag. 2005, 20, 8–14. [Google Scholar] [CrossRef]
  16. Zorzi, M.; Testolin, A.; Stoianov, I.P. Modeling language and cognition with deep unsupervised learning: A tutorial overview. Front. Psychol. 2013, 4, 515. [Google Scholar] [CrossRef]
  17. Hinton, G.E.; Sejnowski, T.J. Learning and relearning in Boltzmann machines. Cogn. Sci. 1986, 9, 147–169. [Google Scholar]
  18. Smolensky, P. Information processing in dynamical systems: Foundations of harmony theory. In Parallel Distributed Processing: Explorations in the Microstructure of Cognition; MIT Press: Cambridge, MA, USA, 1986; Volume 1, pp. 194–281. [Google Scholar]
  19. Hinton, G.E.; Salakhutdinov, R.R. Reducing the dimensionality of data with neural networks. Science 2006, 313, 504–507. [Google Scholar] [CrossRef]
  20. Chopra, P.; Yadav, S.K. Restricted Boltzmann machines and softmax regression for fault detection and classification. Complex Syst. 2017, 4, 67–77. [Google Scholar] [CrossRef]
  21. Liao, L.X.; Jin, W.J.; Pavel, R. Enhanced restricted Boltzmann machine with prognosability regularization for prognostics and health assessment. IEEE Trans. Ind. Electron. 2016, 63, 7076–7083. [Google Scholar] [CrossRef]
  22. Qin, F.C.; Guo, J.M.; Sun, W.D. Object-oriented ensemble classification for polarimetric SAR imagery using restricted Boltzmann machines. Remote Sens. Lett. 2017, 8, 204–213. [Google Scholar] [CrossRef]
  23. Kuremoto, T.; Kimura, S.; Kobayashi, K.; Obayashi, M. Time series forecasting using a deep belief network with restricted Boltzmann machines. Neurocomputing 2014, 137, 47–56. [Google Scholar] [CrossRef]
  24. Wei, S.H.; Xiao, L.X. Research on ballistic missile laser SIMU error propagation mechanism. J. Syst. Eng. Electron. 2008, 19, 356–362. [Google Scholar]
  25. Xian, Y.; Li, G. Research on SINS error transfer model of ballistic missile. Acta Armamentarii 2009, 30, 338–341. [Google Scholar]
  26. Hinton, G.E.; Geoffrey, E. Training products of experts by minimizing contrastive divergence. Neural Comput. 2002, 14, 1771–1800. [Google Scholar] [CrossRef]
  27. Peng, H.; Xiong, Z.; Wang, R.; Liu, J.Y.; Zhang, C. A new dynamic calibration method for IMU deterministic errors of the INS on the hypersonic cruise vehicles. Aerosp. Sci. Technol. 2014, 32, 121–130. [Google Scholar] [CrossRef]
  28. Pan, J.L.; Xiong, Z.; Wang, L.N.; Yu, F.; Zhao, H.; Liu, A.J. A simplified UKF algorithm for SINS/GPS/CNS integrated navigation system in launch inertial coordinate system. Acta Armamentarii 2015, 36, 484–491. [Google Scholar]
  29. Zhou, B.D. Research on Integrated Navigation Technology for Hypersonic Vehicle; National University of Defense Technology: Changsha, China, 2010; pp. 46–50. [Google Scholar]
Figure 1. Overall error estimation scheme of RBM-BPNN method.
Figure 1. Overall error estimation scheme of RBM-BPNN method.
Sensors 19 03682 g001
Figure 2. RBM network structure.
Figure 2. RBM network structure.
Sensors 19 03682 g002
Figure 3. RBM-BPNN network structure.
Figure 3. RBM-BPNN network structure.
Sensors 19 03682 g003
Figure 4. Standard flight path in geocentric coordinate.
Figure 4. Standard flight path in geocentric coordinate.
Sensors 19 03682 g004
Figure 5. Velocity curve of aircraft in inertial coordinate.
Figure 5. Velocity curve of aircraft in inertial coordinate.
Sensors 19 03682 g005
Figure 6. Training flow chart of RBM-BPNN.
Figure 6. Training flow chart of RBM-BPNN.
Sensors 19 03682 g006
Figure 7. Simulation results of initial alignment residual error using RBM-BPNN method.
Figure 7. Simulation results of initial alignment residual error using RBM-BPNN method.
Sensors 19 03682 g007
Figure 8. Simulation results of gyroscope residual error using RBM-BPNN method.
Figure 8. Simulation results of gyroscope residual error using RBM-BPNN method.
Sensors 19 03682 g008
Figure 9. Simulation results of accelerometer residual error using RBM-BPNN method.
Figure 9. Simulation results of accelerometer residual error using RBM-BPNN method.
Sensors 19 03682 g009
Figure 10. Contrast results of initial alignment residual error using RBM-BPNN and IRBM-BPNN method.
Figure 10. Contrast results of initial alignment residual error using RBM-BPNN and IRBM-BPNN method.
Sensors 19 03682 g010
Figure 11. Contrast results of gyroscope residual error using RBM-BPNN and IRBM-BPNN method.
Figure 11. Contrast results of gyroscope residual error using RBM-BPNN and IRBM-BPNN method.
Sensors 19 03682 g011
Figure 12. Contrast results of accelerometer residual error using RBM-BPNN and IRBM-BPNN method.
Figure 12. Contrast results of accelerometer residual error using RBM-BPNN and IRBM-BPNN method.
Sensors 19 03682 g012
Figure 13. Contrast results of initial alignment residual error using three methods.
Figure 13. Contrast results of initial alignment residual error using three methods.
Sensors 19 03682 g013
Figure 14. Contrast results of gyroscope residual error using three methods.
Figure 14. Contrast results of gyroscope residual error using three methods.
Sensors 19 03682 g014
Figure 15. Contrast results of accelerometer residual error using three methods.
Figure 15. Contrast results of accelerometer residual error using three methods.
Sensors 19 03682 g015
Figure 16. Contrast results of attitude error using three methods.
Figure 16. Contrast results of attitude error using three methods.
Sensors 19 03682 g016
Figure 17. Contrast results of velocity error using three methods.
Figure 17. Contrast results of velocity error using three methods.
Sensors 19 03682 g017
Figure 18. Contrast results of position error using three methods.
Figure 18. Contrast results of position error using three methods.
Sensors 19 03682 g018
Table 1. Instrument error parameters of IMU.
Table 1. Instrument error parameters of IMU.
SenorError ParametersError ValueUnits
GyroscopeZero-order error Δ D 0 x 0.1 ° / h
Zero-order error Δ D 0 y 0.1 ° / h
Zero-order error Δ D 0 z 0.1 ° / h
First-order error Δ D 1 x 0.01 ° / ( h g )
First-order error Δ D 1 y 0.01 ° / ( h g )
First-order error Δ D 1 z 0.01 ° / ( h g )
White noise ε g 0.01 ° / h
AccelerometerZero-order error Δ K 0 x 1.0 × 10−4 g
Zero-order error Δ K 0 y 1.0 × 10−4 g
Zero-order error Δ K 0 z 1.0 × 10−4 g
First-order error Δ K 1 x 1.0 × 10−4-
First-order error Δ K 1 y 1.0 × 103-
First-order error Δ K 1 z 1.0 × 103-
White noise ε a 1.0 × 105 g
Table 2. Statistical results of error parameter forecast values using RBM-BPNN method.
Table 2. Statistical results of error parameter forecast values using RBM-BPNN method.
Error ParameterMaximum ValueMinimum ValueMAEMAPE/(%)
Δ φ 0 /(″)12.8593−14.35434.084717.69
ψ 0 /(″)23.1374−26.54277.106927.72
γ 0 /(″)41.0735−46.372312.074228.41
Δ D 0 x /( ° / h )2.8065 × 101−2.8068 × 1017.5766 × 10285.63
Δ D 0 y /( ° / h )2.7374 × 101−2.3656 × 1017.3499 × 10278.21
Δ D 0 z /( ° / h )2.3289 × 101−2.0062 × 1015.3618 × 10262.39
Δ D 1 x /( ° / ( h g ) )2.7299 × 102−2.9003 × 1027.9876 × 10392.51
Δ D 1 y /( ° / ( h g ) )2.6547 × 102−2.7958 × 1027.8484 × 10387.86
Δ D 1 z /( ° / ( h g ) )3.0782 × 102−2.5271 × 1028.1292 × 10381.94
Δ K 0 x /( g )3.2797 × 104−3.1048 × 1047.5387 × 10592.22
Δ K 0 y /( g )2.4197 × 104−3.0791 × 1048.1093 × 10581.02
Δ K 0 z /( g )2.4728 × 104−3.0946 × 1048.6272 × 10577.54
Δ K 1 x 4.7426 × 105−2.9022 × 1051.0119 × 10521.05
Δ K 1 y 3.0925 × 103−2.6130 × 1038.8878 × 10483.01
Δ K 1 z 3.1416 × 103−2.0996 × 1037.4163 × 10479.66
Table 3. Statistical results of error parameter forecast values using IRBM-BPNN method.
Table 3. Statistical results of error parameter forecast values using IRBM-BPNN method.
Error ParameterMaximum ValueMinimum ValueMAEMAPE/(%)
Δ φ 0 /(″)15.1057−25.29374.275217.32
ψ 0 /(″)23.1410−19.36524.238417.57
γ 0 /(″)42.4832−13.42446.380517.24
Δ D 0 x /( ° / h )1.1449 × 101−1.1188 × 1013.2479 × 10252.06
Δ D 0 y /( ° / h )9.7487 × 102−1.2024 × 1013.4799 × 10253.14
Δ D 0 z /( ° / h )1.1065 × 101−1.2577 × 1013.2075 × 10247.47
Δ D 1 x /( ° / ( h g ) )2.4441 × 102−2.7048 × 1027.2427 × 10379.74
Δ D 1 y /( ° / ( h g ) )2.6639 × 102−2.2403 × 1027.9801 × 10376.89
Δ D 1 z /( ° / ( h g ) )2.4528 × 102−2.5661 × 1026.9553 × 10372.46
Δ K 0 x /( g )2.8906 × 104−3.1930 × 1047.6031 × 10582.18
Δ K 0 y /( g )4.8948 × 105−6.0603 × 1051.4845 × 10532.22
Δ K 0 z /( g )1.9109 × 104−2.0527 × 1045.7869 × 10567.25
Δ K 1 x 3.1575 × 105−4.0429 × 1057.2550 × 10617.17
Δ K 1 y 7.3813 × 104−7.6925 × 1042.0235 × 10438.87
Δ K 1 z 4.6047 × 104−4.4187 × 1041.2891 × 10430.32
Table 4. Statistical results of error parameter forecast values of testing samples.
Table 4. Statistical results of error parameter forecast values of testing samples.
Error parameterRBFNNBPNNIRBM-BPNN
MAEMAPE/(%)MAEMAPE/(%)MAEMAPE/(%)
Δ φ 0 /(″)14.471639.476.716923.554.275217.32
ψ 0 /(″)25.394955.026.985525.244.238417.57
γ 0 /(″)26.109344.4712.099230.686.380517.24
Δ D 0 x /( ° / h )5.1322 × 10−263.573.7878 × 10−256.653.2479 × 10−252.06
Δ D 0 y /( ° / h )4.4252 × 10258.663.2558 × 10−251.093.4799 × 10−253.14
Δ D 0 z /( ° / h )4.8451 × 10−258.593.3190 × 10−247.743.2075 × 10−247.47
Δ D 1 x /( ° / ( h g ) )8.3080 × 10−380.477.3119 × 10−379.457.2427 × 10−379.74
Δ D 1 y /( ° / ( h g ) )8.2976 × 10−379.677.1242 × 10−376.237.9801 × 10−376.89
Δ D 1 z /( ° / ( h g ) )8.0855 × 10−377.008.0056 × 10−376.196.9553 × 10−372.46
Δ K 0 x /( g )8.6713 × 10−582.337.7971 × 10−580.147.6031 × 10−582.18
Δ K 0 y /( g )6.0672 × 10−570.411.5862 × 10−534.221.4845 × 10−532.22
Δ K 0 z /( g )7.6700 × 10−578.396.7325 × 10−570.405.7869 × 10−567.25
Δ K 1 x 3.7834 × 10−552.521.4592 × 10−532.657.2550 × 10617.17
Δ K 1 y 4.4197 × 10−459.162.0157 × 10−438.752.0235 × 10−438.87
Δ K 1 z 4.4192 × 10−460.161.5053 × 10−433.971.2891 × 10−430.32
Table 5. Forecast results of error parameter settings by three methods.
Table 5. Forecast results of error parameter settings by three methods.
Error ParameterActual ValueRBFNNBPNNIRBM-BPNN
Estimated ValueRelative Error
/(%)
Estimated ValueRelative Error
/(%)
Estimated ValueRelative Error
/(%)
Δ φ 0 /(″)6069.25084115.4266.37281210.6264.1528526.92
ψ 0 /(″)6022.49255562.5147.19459721.3450.51992615.80
γ 0 /(″)9071.94186320.0666.50860826.1078.08291613.24
Δ D 0 x /( ° / h )0.10.1049014.900.14541645.420.11593215.93
Δ D 0 y /( ° / h )0.10.07445225.550.12027320.270.13698336.98
Δ D 0 z /( ° / h )0.10.07639123.610.07747122.530.08549714.50
Δ D 1 x /( ° / ( h g ) )0.01−0.001552100.000.00383361.670.00484451.56
Δ D 1 y /( ° / ( h g ) )0.010.01337933.790.00264873.520.00127287.28
Δ D 1 z /( ° / ( h g ) )0.010.00239776.030.01304430.440.01251425.13
Δ K 0 x /( g )1.0 × 1046.4540 × 10535.461.2107 × 10587.891.4009 × 10585.99
Δ K 0 y /( g )1.0 × 1045.1649 × 10548.357.0472 × 10529.531.2214 × 10422.14
Δ K 0 z /( g )1.0 × 104−7.7172 × 105100.001.0760 × 1047.607.1746 × 10528.25
Δ K 1 x 1.0 × 1041.3814 × 10586.191.1293 × 10412.939.5368 × 1054.63
Δ K 1 y 1.0 × 1031.9599 × 10480.407.0701 × 10429.308.4009 × 10415.99
Δ K 1 z 1.0 × 1031.5141 × 10351.419.0074 × 1049.938.2878 × 10417.12
Table 6. RMS error statistics of navigation parameters.
Table 6. RMS error statistics of navigation parameters.
Error TypeNavigation ParameterPure SINS ErrorRBFNN ErrorBPNN ErrorIRBM-BPNN Error
Attitude error/(″) φ 56.83647615.9273498.3088035.122294
ψ 108.66551932.57582222.70249112.773308
γ 53.97281324.86488012.4649855.130052
Position error/(m) x 35.62981011.3016812.9679750.636830
y 37.1279532.3932307.2089553.692268
z 90.43341628.69769718.21980110.700961
Velocity error/(m/s) v x 0.6493700.1950830.1082400.028812
v y 0.7094280.1247320.1464940.078337
v z 1.8243310.5273730.3543470.208763

Share and Cite

MDPI and ACS Style

Guo, W.; Xian, Y.; Zhang, D.; Li, B.; Ren, L. Hybrid IRBM-BPNN Approach for Error Parameter Estimation of SINS on Aircraft. Sensors 2019, 19, 3682. https://0-doi-org.brum.beds.ac.uk/10.3390/s19173682

AMA Style

Guo W, Xian Y, Zhang D, Li B, Ren L. Hybrid IRBM-BPNN Approach for Error Parameter Estimation of SINS on Aircraft. Sensors. 2019; 19(17):3682. https://0-doi-org.brum.beds.ac.uk/10.3390/s19173682

Chicago/Turabian Style

Guo, Weilin, Yong Xian, Daqiao Zhang, Bing Li, and Leliang Ren. 2019. "Hybrid IRBM-BPNN Approach for Error Parameter Estimation of SINS on Aircraft" Sensors 19, no. 17: 3682. https://0-doi-org.brum.beds.ac.uk/10.3390/s19173682

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop