A Preferable Airborne Integrated Navigation Method Based on INS and GPS

An integrated navigation method based on INS and GPS was proposed for airborne navigation. The influence of scale factor error and misalignment error of gyroscope and accelerometer on navigation accuracy was analyzed. Compared with traditional INS/GPS integrated navigation method, scale factor error and misalignment error were added to the state model of the integrated navigation system. The observability of scale factor error and misalignment error was analyzed combined with typical airborne movement. Then the integrated system was optimized, and the new navigation model of the integrated system was obtained. The optimized INS/GPS integrated model was validated by numerical simulation and turntable test. Comparing the proposed model with traditional integrated model (integrated system error states do not include scale factor error and misalignment error), the results showed that the proposed integrated navigation method can improve the accuracy from 8% to 28% of the east, north, and upward positions.


Introduction
Inertial navigation system (INS) is an autonomous navigation system that does not depend on any external information [1,2].However, the characteristics of the location error accumulate with time, making it difficult to work independently for a long time.Global Positioning System (GPS) can measure three-dimensional position and velocity accurately, but the disadvantage is susceptible to interference and control [3][4][5].Therefore, INS and GPS have complementary characteristics.Since the 1990s, INS/GPS integrated navigation system has been a great success at home and abroad, and it has developed into a specialized technology [6,7].
INS/GPS integrated navigation system works as follows: when GPS signal is good, the system selects the integrated navigation mode.The precision of integrated navigation basically depends on GPS precision, and the inertial measurement unit (IMU) errors can be estimated and compensated online.When GPS signal is disturbed or shielded, the system automatically shifts into inertial navigation mode.At this point, navigation accuracy basically depends on the precision of IMU [8].Therefore, the estimation accuracy of IMU errors in integrated navigation can affect the accuracy of the subsequent inertial navigation [9,10].In airborne INS/GPS integrated navigation system, however, the errors of IMU only consider the bias of gyroscope and accelerometer without considering scale factor and misalignment at present.A method of the dynamic parameter identification of the scale factor error and misalignment error was designed based on Kalman filter.The observability of the scale factor error and misalignment error with different maneuvers was analyzed in [11].Zhou et al. described the error dynamic system equation and observation equation of inertial navigation system and the singular value of the system states of online calibration [12].Therefore, a more advanced method can be designed to make the IMU errors (including bias, scale factor, and misalignment) be more accurate in estimation and compensation in the integrated navigation process.When entering the inertial navigation mode, it can get higher navigation accuracy.
For the application of airborne navigation, this paper analyzes the influence of scale factor error and misalignment error on the accuracy of integrated navigation.Based on the analysis, the scale factor error and the misalignment error are added to the error model of the integrated navigation system.Then the observability of the scale factor error and misalignment error is analyzed combined with the typical airborne movement.According to the observability analysis results, the integrated system is optimized and the new error model of the integrated navigation system is obtained.Finally, the optimized INS/GPS integrated model proposed in this paper is validated by numerical simulation and turntable test, and then the proposed model is compared with the traditional integrated model.

Error Analysis of Airborne Inertial Measurement Unit
The errors of inertial measurement unit mainly include bias, scale factor error and misalignment error.The error model of INS and IMU is given in the following passage, and the influence of IMU errors on the navigation accuracy is analyzed combined typical airborne movement.

Error Model of Inertial Navigation System
.
where δλ, δL, and δh represent the longitude error, the latitude error, and the altitude error, respectively, and R M and R N are the curvature radius of the meridian and prime vertical.δv n and φ represent velocity error and attitude error, respectively, and δf n represents the accelerometer measurement error, which contains accelerometer bias, scale factor error, and misalignment error.ε n represents the gyroscope measurement error, which contains gyroscope bias, scale factor error, and misalignment error.

Analysis of the Influence of IMU Errors on Navigation
Accuracy.In order to facilitate quantitative analysis, combined with the actual low-precision inertial navigation systems commonly used in airborne navigation, the error parameters of the IMU are set as follows: When the airframe turns, the angular velocity of the Earth and the platform are small relative to the IMU rotation rate.So we ignore the influence of angular velocity of the Earth and the platform on the attitude error; the attitude error equation is further simplified [14]: Journal of Sensors When the system only changes the pitch angle, the roll angle and heading angle change could be assumed to be zero, then Combined with the gyroscope error parameters given above and (8), the influence of gyroscope scale factor error gSF x on attitude error is greater than that of gyroscope bias gB x on attitude error when ω x > 0 93 °/s; the influence of gyroscope misalignment errors gMA yx and gMA zx on attitude error is greater than that of gyro bias gB x and gB z on attitude error when ω x > 1 43 °/s.Therefore, when the body pitch angle changes, the influence of gyroscope scale factor error and misalignment error on attitude error cannot be ignored.
Similarly, when the system roll angle changes or heading angle changes, then ω b nb is Substituting ( 9) into ( 6), we can get the same conclusion.

Error Analysis of Velocity and Position.
The most common movements of the aircraft are uniform motion and accelerated motion, this paper mainly analyzes the influence of the accelerometer errors on position error and velocity error in these two kinds of motion.To simplify the analysis, this paper also does not consider other error terms, so velocity error equation can be simplified as When the system only changes the pitch angle, the roll angle and heading angle change could be assumed to be zero; the attitude matrix C n b can be simplified as Substituting ( 11) into (10), In uniform motion (f x = f y ≈ 0), according to the accelerometer error parameters and (12), we can conclude that aM A ⋅ g ≈ 194 μg < aB = 500 μg.The influence of accelerometer misalignment errors aMA xz and aMA yz on the east and north velocity errors is about 40% of the bias aB x and aB y .In addition, aSF z ⋅ g ≈ 300 μg < aB z = 500 μg.The influence of the accelerometer scale factor error aSF z on upward velocity error is about 60% of accelerometer bias aB z .Therefore, the influence of the accelerometer scale factor error and misalignment error on velocity error cannot be ignored in uniform motion.
In accelerated motion (f x ≠ 0, f y ≠ 0), according to the accelerometer error parameters and (12), we can conclude that when f y > 1 67g and f z > 1 67g, the influence of the accelerometer scale factor errors aSF y and aSF z on the north velocity error and upward velocity error is greater than that of accelerometer bias aB y and aB z .When f y > 2 57g and f z > 2 57g, the influence of the accelerometer misalignment errors aMA zy and aMA yz on the north velocity error and upward velocity error is greater than that of accelerometer bias aB z and aB y .In the actual motion, the acceleration of aircraft is about 0.3g.According to the analysis above, the influence of the accelerometer scale 3 Journal of Sensors factor error on velocity error is about 18% of the accelerometer bias and the influence of accelerometer misalignment error on velocity error is 12% of the accelerometer bias in accelerated movement.
Similarly, when the system's roll angle changes or heading angle changes, attitude matrix Substituting ( 13) into (10), we can get the same conclusion.
In summary, the influence of scale factor error and misalignment error of gyroscope and accelerometer on navigation cannot be ignored.

INS/GPS Integrated Navigation Model.
From what has been analyzed above, we can conclude that the influence of scale factor error and misalignment error on navigation cannot be ignored in airborne navigation.Therefore, the scale factor error and misalignment error need to be considered in the INS/GPS integrated error model.In this paper, INS/GPS integrated navigation is realized by using Kalman filter.The position and velocity differences between GPS and INS are taken as measurement errors.IMU errors and navigation errors can be online estimated and compensated, so the high navigation accuracy can be acquired.INS/GPS integrated navigation schematic diagram is shown in Figure 1.

Error State Model of Integrated Navigation.
The integrated system error model can be expressed as X t is the 33-dimensional error vector of integrated system which can be expressed as The velocity and position differences between GPS and INS are taken as measurement errors in the process of integrated navigation.So the measurement equation can be expressed as ,  Considering the complexity of the procedure, this paper does not simulate the whole flight process but chooses nine typical stages [18,19].The flight stages of aircraft and simulation trajectory are shown in Table 1 and Figure 2.
With the combined typical flight trajectory above, according to the eigenvalues of the error covariance matrix, the observability of the scale factor error and misalignment error of the gyroscope and accelerometer is analyzed in different maneuvering conditions.And the influence of different maneuvering conditions on observability of each error term is analyzed.The eigenvalues of the covariance matrix can reflect whether the system state estimation is good or bad.The smaller the eigenvalue is, the smaller the estimated variance of the corresponding state is and the better the degree of observability is.Otherwise, the estimation accuracy is low and the degree of observability is poor.
Figures 3 and 4 are the normalized eigenvalues of the corresponding errors in the covariance matrix during the whole motion.Tables 2 and 3 are the normalized covariance of scale error and misalignment error of gyroscope and accelerometer.
As can be seen from Figure 3 and Table 2, the scale factor error and misalignment error of gyroscope are not observable during gliding.When climbing, the degree of observability of gSF x is increased.In a uniform flight, the flight time is long and the observability of gSF x and gMA yx is significantly increased.The degree of observability of gSF y and gMA yz is increased when circling to the left and the degree of observability of gMA xz is increased during circle exit.It can be seen that rotational maneuvering has a good incentive   As can be seen from Figure 4 and Table 3, the degree of observability of scale factor error aSF z is increased during accelerated gliding.The change of the observability degree of scale factor error and misalignment error is not obvious during climbing and transformation level.In a uniform flight, the degree of observability of  7 Journal of Sensors and aMA zy are unobservable.After the unobservable errors are deleted, the dimensions of the error states is reduced from 33 to 24 and the optimized error state vector of the integrated system is So the system noise vector and the noise driving matrix of the integrated system are W t = 0 1×3 aW x aW y aW z gW x gW y gW z 0 1×15 , The measurement matrix is expressed as In this way, the optimized integrated model is obtained, and then the simulation and turntable test are used to verify the optimized model, respectively.

Model Verification
4.1.Simulation Verification.In this paper, numerical simulation is used to verify the validity and correctness of the proposed INS/GPS integrated model and then compare it with traditional integrated model.
The error parameters of IMU and GPS are as follows: (1) Gyroscope (i) Bias: 1 °/h; noise: 0 05 °/ h The simulation time is 700 seconds, the first 500 seconds for INS/GPS integrated navigation.GPS failure after 500 seconds, the receiver cannot output velocity and position information.At this time, integrated system cannot update measurement equation and turns into inertial navigation mode.The simulation results are as follows, Figure 5 shows the estimation of gyroscope bias, scale factor error, and misalignment error in integrated navigation.Figure 6 shows the estimation of accelerometer bias, scale factor error, and misalignment error in integrated navigation.Figure 7 shows the true value, estimation of position and position error of the traditional integrated navigation model, and the proposed navigation model in the whole simulation process, including the east position, the north position, and the upward position.
From the simulation results above, the set errors of gyroscope bias, scale factor, and misalignment are 1 °/h, 300 ppm, and 40 ″ , respectively.Figure 3 shows that the error estimation of the gyroscope are 0.96 °/h, 305 ppm, and 42.8 ″ , respectively.The set errors of accelerometer bias, scale factor, and misalignment are 500 μg, 300 ppm, and 40 ″ , respectively.Figure 4 shows that the error estimation of the accelerometer are 482 μg, 297 ppm, and 43 ″ , respectively.Therefore, the observable errors of IMU are effectively estimated in integrated navigation.
GPS signal is lost after 500 seconds; system conducted 200 seconds inertial navigation at 700 seconds.Figure 7 shows that the east, north, and upward position errors of the traditional model are 380 m, 231.4 m, and 126.1 m, respectively.The east, north, and upward position errors of the proposed model are 271.2m, 175.1 m, and 111.5 m, respectively.The navigation accuracy of the east, north, and upward positions of the proposed model improved by 28.6%, 24.3%, and 8.3%, respectively.4.2.Turntable Test.Numerical simulation has proved the validity of the proposed integrated navigation method.In order to further verify the effectiveness of the method in actual integrated navigation system, turntable test is designed on the basis of numerical simulation.Test equipment include INS/GPS integrated navigation system, turntable, two-way DC power supply (0~30 V, 0~3 A), and data acquisition 9 Journal of Sensors computer.Turntable is one of the important equipment in test for simulating the change of attitude and attitude rate of aircraft in space.INS/GPS integrated navigation system consists of inertial navigation system and GPS receiver.The installation of the main test equipment is shown in Figures 8 and 9.
The error values of IMU are shown in Table 4.
Based on the analysis of the typical airborne trajectory, turntable is used to simulate the typical airborne movement.The rotation parameter of the turntable during the test is shown in Table 5.
The simulation time is 700 seconds, the first 500 seconds for INS/GPS integrated navigation.Integrated system turns into inertial navigation mode after 500 seconds.The test results are as follows, Figure 10 shows the estimation of gyroscope bias, scale factor error, and misalignment error in integrated navigation.Figure 11 shows the estimation of accelerometer bias, scale factor error, and misalignment error in integrated navigation.Figure 12 shows the estimation of position error of the traditional integrated navigation model and the proposed navigation model in the whole process, including the east position error, the north position error, and the upward position error.
From Figures 10-12, observable error estimation status of IMU in the integrated navigation process is shown in Table 6.
As can be seen from Table 6, there are 7 observable errors of IMU of which the estimation accuracy is more than 80%.There are 4 observable errors of which the estimation accuracy is between 50% and 80%, and there are 4 observable errors between 20% and 50%.Consequently, most of the observable errors of IMU are effectively estimated in the integrated navigation process.
From the test results above, we can conclude that after 500 seconds integrated navigation and 200 seconds inertial navigation, the position errors of the east, north, and upward of the traditional model are 798.3m, 550.6 m, and 207 m, respectively.The east, north, and upward position errors of the proposed integrated model are 592.2m, 444.8 m, and 176.6 m, respectively.The navigation accuracy of the east, north, and upward positions of the proposed model improved by 25.8%, 19.2%, and 14.7%, respectively.Consequently, the integrated navigation method proposed in this paper can achieve higher navigation accuracy.So the integrated navigation model proposed in this paper is superior and effective.

Conclusions
The influence of scale factor error and misalignment error on navigation accuracy is analyzed in this paper.Based on the analysis, scale factor error and misalignment error are considered in the error state vector.Then the observability of scale factor error and misalignment error is analyzed combined with typical airborne movement.The integrated system is optimized according to the observability analysis results.Finally, this method is verified by numerical simulation and turntable test.The results all show that the INS/ GPS integrated navigation model proposed in this paper

4 2. 3 . 1 .gSF x gMA xy gMA xz gMA yx gSF y gMA yz gMA zx gMA zy gSF z
gSF x = gSF y = gSF z = 300 ppm, gMA xy = gMA xz = gMA yx = gMA yz = gMA zx = gMA zy = 40″, aB = 500 μg, aSF x = aSF y = aSF z = 300 ppm, aMA xy = aMA xz = aMA yx = aMA yz = aMA zx = aMA zy = 40 ″ Analysis of Attitude Error.In order to more clearly and easily analyze the influence of gyroscope scale factor error and misalignment error on attitude error, we temporarily do not consider other terms.The attitude error equation can be simplified as φ = −C n b

3. 2 . 1 .
Observability Analysis.In this paper, the movement of aircraft is divided into the following three typical processes: straight flight, climbing flight, and turning flight.Straight flight includes uniform motion and accelerated motion.Climbing flight includes preparation for climbing, climbing, and transformation level.Turning flight includes tilting and spiral.Flight trajectory consists of the three processes above.

Figure 2 :
Figure 2: The flight trajectory of aircraft.

Figure 4 :
Figure 4: Eigenvalues of covariance matrix corresponding to accelerometer errors.
gB x , gB y , and gB z are the gyroscope biases.ω b ibx , ω b iby , and ω b ibz are the gyroscope ideal outputs., aSF y , and aSF z are the accelerometer scale factor errors of the x-axis, y-axis, and z-axis.aMA xy , aMA xz , aMA yx , aMA yz , aMA zx , and aMA zy are the accelerometer misalignment errors.aB x , aB y , and aB z are the ibz , 2 where gSF x , gSF y , and gSF z are the gyroscope scale factor errors of the x-axis, y-axis, and z-axis.gMA xy , gMA xz , gMA yx , gMA yz , gMA zx , and gMA zy are the gyroscope misalignment errors.
x aB y aB z gB x gB y gB z aSF x aSF y aSF z aMA xy aMA xz aMA yx aMA yz aMA zx aMA zy gSF x gSF y gSF z gMA xy gMA xz gMA yx gMA yz gMA zx gMA zy T , 15 W t is the system noise model which can be expressed as W t = 0 1×3 aW x aW y aW z gW x gW y gW z 0 1×24 T , 16 where aW x , aW y , and aW z are accelerometer noises, gW x , gW y , and gW z are gyroscope noises.G t is the noise driving matrix of integrated system, Scale factor error and misalignment error are added in the new integrated model.However, the scale factor error and misalignment error cannot be observed completely.It is difficult to obtain the desired estimation results if the estimated states cannot be observed.So the observability of the scale factor error and misalignment error is analyzed by a combined typical airborne trajectory.Based on the analysis of observability, we delete unobservable error states and get the optimized integrated navigation model.

5
Journal of Sensors for gSF y , gMA yz , and gMA xz .So gSF x , gSF y , gMA xz , gMA yx , and gMA yz are observable.Moreover, at the end of the whole movement, the normalized covariance of gSF z , gMA xy , gMA zx , and gMA zy are decreased not obviously.In summary, gSF z , gMA xy , gMA zx , and gMA zy of gyroscope errors are unobservable.
x aB y aB z gB x gB y gB z aSF y aSF z aMA xy aMA zx gSF x gSF y gMA xz gMA yx gMA yz 199

Table 3 :
Normalized covariance of accelerometer scale factor and misalignment.

Table 2 :
Normalized covariance of gyroscope scale factor and misalignment.

Table 6 :
Observable error estimation status of IMU.