Mathematical Problems in Engineering

Volume 2013, Article ID 524284, 12 pages

http://dx.doi.org/10.1155/2013/524284

## A Fast and High-Accuracy Compass Alignment Method to SINS with Azimuth Axis Rotation

Key Laboratory of Micro-Inertial Instrument and Advanced Navigation Technology, School of Instrument Science & Engineering, Ministry of Education, Southeast University, Nanjing 210096, China

Received 30 September 2012; Accepted 10 January 2013

Academic Editor: Gerhard-Wilhelm Weber

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

#### Abstract

Azimuth axis rotating modulation was introduced to improve the alignment accuracy of strapdown inertial navigation system (SINS) through compass algorithm, in which the limit accuracy was determined by equivalent sensor errors in the eastern and northern direction. In this modulation, horizontal sensor errors were modulated into zero mean periodic variables. Furthermore, two methods were introduced to ensure alignment accuracy and speed: (1) shortened rotating cycle and redesigned compass parameters were selected to eliminate or ease the amplification to low-frequency senor error inputs in compass loop caused by rotation and (2) a data repeated calculation method was designed to shorten prolonged alignment time caused by the above redesigned parameters. Based on a certain SINS, turntable test proves that alignment accuracy and time were significantly improved and slightly shortened in comparison with the classical compass alignment.

#### 1. Introduction

Rotating modulation in the strapdown inertial navigation system (SINS), being successfully applied since 1980s, is an effective way to improve the accuracy of SINS with certain accuracy inertial sensors. In 1980, Levinson and Giovanni introduced a strapdown ship compass, in which uniaxial 4-position stop program was used to compensate for the drifts of ring laser gyros, thus significantly improving the system accuracy [1]. Subsequently, Sperry developed MK-49 dual-axis rotating modulated laser SINS with 16-order rotating program and this system was selected as a standard ship’s INS for the North Atlantic Treaty Organization (NATO) in 1989 [2, 3].

At the same time, theories about rotating modulation were concern. The widely accepted views about rotating modulation are as follows [1–4]: with uniaxial rotation, symmetry sensor errors installed at the vertical direction to the rotating axis can be eliminated and with biaxial rotation, all symmetry sensor errors can be eliminated. In 2007 and 2012, Yuan et al. proposed an improved 4-position program for uniaxial rotation and 16-order rotating program for biaxial rotation [5, 6], to further improve the system accuracy compared with the methods in the literatures of [1–3]. In existed researching files, some scholars aim to improve the navigation accuracy by rotating scheme designing [4–6] and others focus upon the added observability caused by rotation when the optimal estimation method based on Kalman filter is used for initial alignment [7–11]; however, less attention was put to the possibility of improving alignment accuracy of compass alignment method brought by rotating modulation. Therefore, apart from study on azimuth rotating modulation, more research goes to oscillation in horizontal alignment, divergence in azimuth, and time-consuming of alignment process and corresponding solutions were offered.

In SINS, through classical compass alignment method, the horizontal alignment accuracy is subjected to the equivalent horizontal accelerometer errors while the azimuth is subject to the equivalent eastern gyro error [11–14]. Applied in engineering, these equivalent errors are determined by senor errors and the carrier’s attitude together. While with azimuth rotation, the horizontal sensor errors are modulated into zero mean periodic variables [5], which brings the possibility to improve the alignment accuracy of SINS through compass method. But the studies in this paper indicated an amplification effect in compass loop upon the inputs of low-frequency sensor errors that are constant values if no rotation is involved. Such effect brought certain oscillation in horizontal alignment and divergence in azimuth. With the accelerated rotating speed and redesigned compass parameters, the amplification effect was eliminated or eased. Meanwhile, due to the above redesigned parameters, the alignment time was greatly prolonged.

To meet the requirement of SINS’s rapid response, a fast alignment method based on data repeated calculation was designed. In this method, a characteristic of mathematical platform, cutting off the direct relationship between sensor data and feedback signals for alignment loop, in SINS is used, which is different from the relationship in PINS. With this characteristic and coupled with compass alignment control loop, the same set of data can be solved repeatedly to produce different direct cosine matrix (DCM) and different projected sensor data in mathematical platform from measured sensor data in body frame, thus producing different feedback signals. Moreover, a turntable test with one certain type of SINS was used to certify the effectiveness of these above methods.

The article is organized as follows: In Section 2, classical compass alignment to SINS will be introduced and simulation results will be given. In Section 3, azimuth rotating modulation will be introduced, simulation about compass alignment with rotation will be done, and the failure reasons will be analyzed; in the meantime, improved methods proved by simulation will be illustrated. In Section 4, a fast compass alignment based on data repeated calculation will be put forward, and simulation about this method will be done. In Section 5, a turntable test with a certain type of SINS will be taken. Lastly, conclusions will be drawn in Section 6.

#### 2. Classical Compass Alignment Algorithm to SINS

##### 2.1. Classical Compass Alignment Algorithm

With the development of SINS, compass alignment method has been widely used for its self-alignment. The main difference between SINS and platform INS (PINS) lies in the application of mathematical platform in SINS while that of physical platform in PINS, that is to say, the sensor data in PINS are used as feedback signals while the projected sensor data in mathematical platform are used in SINS. But in terms of alignment process, the mechanisms are the same in both INSs.

Figures 1(a) and 1(b) are the northern and azimuth channels of compass fine alignment loop in SINS, respectively. And in Figure 1, , , and represent the inertial frame, the body frame, and the Earth frame, and the subscript represents the control signal from alignment loop, and the “^{?}” and “~” represent the calculation and measurement value of a vector respectively.

In Figure 1(a), the calculated value of DCM contains the eastern misalignment . With this DCM , Earth gravity will be projected into the calculated value of northern acceleration . When is integrated, north velocity error that can be set as a feedback signal in the northern channel to eliminate will be led. In this channel, the damping term (dot lines) is applied to realize weak oscillation from equal amplitude oscillation, the forward-feed term (dash-dot lines) is applied to shorten the oscillation period to , and the integration term (double dash-dot lines) is introduced to eliminate the effect of .

In Figure 1(a), when the rapid and stable requirements of alignment process are comprehensively considered, control parameters can be designed by the following steps: selecting the damping coefficient and the system adjusting time , calculating the damped oscillation angular rate and the attenuation coefficient , and calculating the control parameters as where is angular rate of Schuler.

Similarly, in Figure 1(b), the parameters for azimuth channel can be calculated as

With the above parameters in (1) and (2), the limit accuracy of , and can be achieved as follows: where is the Earth angular rate, is the latitude, and , , and are equivalent sensor errors in northern and eastern direction in the navigation frame.

In SINS equivalent errors in (3) are jointly determined by sensor errors and carrier’s attitude, while those in PINS are determined by sensor errors. Such difference indicates that different carrier’s attitude leads to different alignment accuracy by same sensor in SINS.

##### 2.2. Simulation

*(1) Simulation Setting. *To set the system damping parameter as 0.707 and the adjusting time of horizontal loop as 150?s, , , can be obtained. To set the adjusting time of azimuth loop as 300?s, , , can be gotten. In alignment, the horizontal alignment process lasts for the first 150?s and the azimuth alignment covers the rest.

The horizontal attitudes of the carrier are set as 0°, the azimuth as 0°, 45°, 90°, 135°, 180°, respectively, all initial misalignment as 0°, and initial latitude and longitude as 32° and 118°. The gyroscope constant and random errors are set as 0.05°/h and 0, respectively, and the accelerometer constant and random errors as 500?ug and 0, respectively, and the height in navigation as 0 constantly. With above assumption, the equivalent sensor errors in navigation frame can be calculated as shown in Table 1.

*(2) Simulation Results*. The simulation lasts 1200?s, and simulation results are shown in Figure 2. The curves in Figure 2 indicate that about 200?s is needed to fulfill horizontal alignment and about 600–800?s for yaw. Statistics from 1000?s to 1200?s are shown in Table 2.

?The curves in Figure 2 and statistics in Table 2 indicate that the alignment accuracy was jointly determined by sensor errors and carrier’s attitude. In order to improve alignment accuracy with certain accuracy sensors and to eliminate the effect from the carrier’s attitude, more methods are required.

#### 3. Compass Alignment Algorithm with Azimuth Axis Rotation

##### 3.1. Mechanism of Azimuth Rotational Modulation

In azimuth axis rotational modulation, the inertial measurement unit (IMU) is installed on the rotational platform that can rotate around the -axis. Under this circumstance, the sensors measure the carrier’s movement in the rotational frame . Theses measured angular rate and acceleration are supposed to be initially projected into the carrier body frame to participate in navigation solution:

Neglecting all other interferences and concerning sensor errors only, the above equation can be expanded as

Equation (5) indicates that the sensor errors have been modulated with rotation. When the rotation angular rate of rotational platform is set as , the attitude matrix can be constructed as Then and can be expanded as

The above equations in frame show that the sensor errors along -axis and -axis are modulated into zero mean periodic variables with azimuth axis rotation.

##### 3.2. Simulation for Classical Compass Alignment with Azimuth Rotation

In SINS, the above modulated error in frame must be projected into frame with DCM to get involved in the navigation solution, which will bring some negative effect to the above modulation. For ships or cars, the negative effects caused by the above projection are limited because of their generally small horizontal attitude. Therefore, compared with the nonrotated modulation, the navigation accuracy of SINS can still be significantly improved.

Meanwhile, the above modulation brings the possibility to improve compass alignment accuracy. In the following simulation, azimuth rotation is brought into classical compass alignment method.

*(1) Simulation Setting*. The simulation is conducted under the following conditions: the control loop parameters and the sensor errors set as same as those in Section 2.2, the initial attitude and initial misalignments set as zero, the rotating period of IMU set as 240?s, and no horizontal angle and linear movement.

*(2) Simulation Results*. The simulation lasts 2100?s and the simulation results are shown in Figure 3 as dot-dash lines. In Figure 3, the alignment results in Figure 2 with the carrier’s attitude 0° are redrawn in solid lines as a reference data. From Figure 3, conclusion can be drawn that the horizontal attitude error curves with azimuth rotation are oscillating with the period of 240?s and the oscillation amplitudes are approximate to the reference accuracy with nonrotation. The azimuth curve is oscillating with a period of 240?s, but its oscillation amplitude is larger than that of nonrotation. Such a result indicates the failure of the simple combination of classical alignment and rotating modulation.

From Figure 3, the conclusion can be reached that negatively influential elements that result in the oscillation of horizontal alignment and failure of azimuth are sure to be found in the process of classical compass alignment with azimuth rotation.

##### 3.3. Further Analysis and Improved Methods

###### 3.3.1. The Reasons for Alignment Failure

Take the azimuth alignment loop as an example to analyze the reasons. In Figure 1(b), set as ; then where and is the equivalent upward gyro error.

With the azimuth axis rotation, the horizontal sensor errors are modulated into the periodically changing values; in other words, to the azimuth loop, the inputs are changed from constant values to periodically changing values.

According to (8), TF1 and TF2 are defined as transfer functions from the input and to the output . And the Bode charts of TF1 and TF2 with the parameters in Section 2.2 are shown in Figures 4(a) and 4(b), respectively.

In Figures 4(a) and 4(b), the magnification of the equivalent eastern gyro error is 93.6?dB and that of the equivalent northern accelerometer is 42.2?dB with the errors modulation period 240?s (the angular velocity is 0.0262?rad/sec). Therefore, alignment failure tends to result mostly from large magnification.

###### 3.3.2. Improved Methods

To reduce the above magnification effect, two options can be selected: to speed up the rotating angular rate, the input frequency will move from left to right in the Bode chart and to adjust compass parameters will make the high-magnification range move from right to left part in Bode chart.

*(1) Method One*. 120?s is deemed as rotating cycle with respect to the increase of burden of mechanical indexing mechanism and other sensor errors followed by too fast rotating rate. Then the magnification of the equivalent eastern gyro error is 64.3?dB and that of the equivalent northern accelerometer is 18.8?dB as shown in Figures 4(a) and 4(b). The simulation results are presented in Figure 5 as with solid lines. Compared with the results in Figure 3, conclusion can be drawn that the horizontal alignment errors are suppressed and the azimuth error curve is oscillating and its oscillation amplitude is approximate to the reference value. Simulation results through this method indicate that the alignment accuracy is increased slightly.

*(2) Method Two*. ?In order to restrain the oscillation amplitude of horizontal alignment error and realize the convergence of azimuth alignment error, the adjusting time of horizontal loop is set as 600?s and that of azimuth is set as 800?s. Then the parameters for horizontal loop are calculated as , , and , and the those for azimuth are as , , . Besides, the rotating period is still 240?s. The horizontal alignment lasts for the first 600s and the azimuth alignment lasts for the rest in alignment process. With these parameters, the magnification of the equivalent eastern gyro error is 60.2?dB and that of the equivalent northern accelerometer is 8.94?dB as shown in Figures 4(c) and 4(d). The simulation results are presented as dot-dash lines in Figure 5. Compared with those solid lines, the horizontal oscillation is obviously suppressed and its amplitude is decreased. The azimuth error curve converges slowly and tends to 0, but the alignment time is extended to around 1900?s.

*(3) Method Three*. In order to combine the advantages of the above two methods, 120?s is set as the rotating period and the alignment parameters are set as those in method 2. With the above parameters, the magnification of the equivalent eastern gyro error is 30.2?dB and that of the equivalent northern accelerometer is -15.2?dB as shown in Figures 4(c) and 4(d). The simulation results are presented as dotted line in Figure 5. Compared with those solid lines and dot-dash lines, the oscillation amplitude is further decreased and the convergence period of azimuth is shortened to about 1500?s. It is confirmed that the combination of the above two methods is effective.

In Table 3, the mean data from 2000?s to 2100?s in method 3 are listed and the limit accuracy of compass algorithm without rotation when the carrier’s attitudes are 0° is also listed.

Compared results in Table 3 indicate that alignment accuracy can be significantly improved with azimuth rotation. But the curves in Figure 5(c) also show that the alignment time is too long.

In engineering, the alignment time may be longer due to the various interference factors and initial misalignments. This, for the sake of shortening alignment time, requires more available methods.

#### 4. A Fast Compass Alignment Method Based on Data Repeated Calculation

##### 4.1. A Possible Way to Accelerate Compass Alignment Speed of SINS

In Figure 6(a), inertial sensors are installed on the physical navigation platform in PINS and the difference between the sensor measured data and the ideal data from ideal navigation frame is decided by the misalignments. In other words, there are direct corresponding relationships between the sensor data and feedback correction signals.

While in SINS, inertial sensors are installed on the carrier and a mathematical platform is used to simulate the physical navigation platform as shown in Figure 6(b), in which is used to denote mathematical calculated frame from the calculated DCM . Before alignment, there are misalignments between and the ideal navigation frame . Only after the sensor data obtained from frame is projected into frame can it get involved in compass alignment operation. That is, no direct corresponding relationship between the sensor data and feedback signals is found. On such condition different attitude matrixes determine different projection data in frame and further determine the feedback signals with the same sensor data.

This difference between the compass methods in SINS and PINS indicates that a same set of data collected in frame can be used to fulfill alignment by the repeated calculation for , that is to say, sensor data collecting time will be reduced and alignment time will be shortened with use of computer’s high-speed performance at most.

##### 4.2. Repeated Calculation Algorithm in SINS

###### 4.2.1. Navigation Resolution Algorithm in SINS

The recurrence equations for SINS navigation resolution can be expressed as where is recursive number, is the sensor data sampling and navigation resolution update cycle, , , “” denotes the projection of a motion vector that means the relative motion from to frame frame in frame, is the DCM, is the velocity vector, , , and are the latitude and longitude and height, respectively and denotes the antisymmetric matrix of the vector “”. For ships and cars, the height can be set as zero constantly.

###### 4.2.2. A Repeated Calculation Algorithm for Compass Alignment

The above equations (9) constitute the basic update equations of navigation calculation. The sensor data are of the “real-time” meaning when involved in the recursion operation or losing their meaning in the next period after the completion of the current navigation calculation. But as the analysis in Section 4.1, initial alignment of SINS can be completed by the repeated calculation for with a set of data.

When the carrier is in the angular but nonline motion, which is common state for compass alignment, a repeated calculation algorithm is designed as shown in Figure 7.

The normal collection for sensor data and normal navigation calculation is done in the period from to and the sensor data is stored at the same time. In addition, the forward repeated calculation is done in the subsequent period from to or the period from to according to Figure 7.

In the period from to , is the starting time and the forward calculation has the same calculation process as in the normal navigation resolution. If the carrier is in nonline motion, the velocity and the position of it at can be set as the initial velocity and the position at . However, the attitude at differs from that at by reason of its angular movement. What’s more, an instantaneous solidification body frame is introduced to take the carrier body frame at as the solidification body frame , which means . Then, the can be updated by gyro data:

With (10), initial DCM at can be obtained by the matrixes and :

In the period from to , similar processes can be conducted as the period from to . With the above process, in the period from to , the vehicle can move from the origin to the end, just as that in the period from to .

###### 4.2.3. Real-Time Attitude Calculation at the End of Repeated Calculation

When the above repeated calculation ends, calculated DCM is for the time , not for the time . In order to ensure the real-time property of SINS navigation, DCM for the time must be obtained. To meet this requirement, solidification body frame in Section 4.2.2 is also used. A new instantaneous solidification body frame is defined at time and can be updated with gyro data. Therefore, when repeated calculation ends, the real-time DCM can be obtained as

##### 4.3. Simulation for Compass Alignment Based on Repeated Calculation

*(1) Simulation Setting*. The simulation is carried out with the following conditions: the gyro constant and random errors set as 0.05°/h and the accelerometer constant error and random errors as 500?ug, the amplitude of pitch, rolling, and azimuth set as 2°, 5°, and 2.5° and the swing period of pitch rolling, and azimuth as 4?s, 6?s, and 4?s, respectively, the initial azimuth of carrier set as 0 and its initial horizontal attitude as 0, the initial misalignment of pitch, rolling and azimuth set as -1°, -1°, and 3°, respectively, the initial velocity and acceleration set as 0 with its initial latitude and longitude as 32° and 118°, respectively, the initial height as 0 and the height set as 0 constantly in the calculation, the first 120?s data stored to do the repeated calculation (120?s cycle period), and the compass loop parameters set as those in Section 2.2 and the update cycle as 10?ms.

*(2) Simulation Results*. The simulation lasts 1200 s, during which two different simulation scheme results are compared. Moreover, the second scheme is based on classical compass alignment with data repeated calculation while the first one is not. In Figure 8, (a1), (b1), and (c1) are the alignment error of pitch, roll, and yaw of Scheme 1, respectively, and those of Scheme 2 are likewise.

In Figure 8, the curves in (a1) and (a2), (b1) and (b2), and (c1) and (c2) present the same trends. In other words, the effect of the compass alignment based on the date repeated calculation matches that of the classical compass alignment. Thus, alignment time can be shortened by conducting the repeated calculation to a set of sensor data with high-speed computer calculation.

#### 5. Turntable Test

Based on the azimuth axis rotating modulation and repeated calculation, the fast and high-accuracy compass alignment algorithm of SINS proposed in this paper has been verified by the triaxial swing turntable test. The test was run as semiphysical simulation with data collected from the turntable and SINS.

##### 5.1. Test Condition

In this test, the rate controlling accuracy of the turntable is ±0.0005°/s and angle measuring accuracy is ±0.0001°/s.

Developed by Southeast University, flexible gyros and quartz accelerometers are used in this type of SINS. The accuracy of gyros and accelerometers are provided in Table 4. Besides, the frequency of navigation computer is 333?MHz. The IMU is mounted on the turntable as shown in Figure 9. The update frequency of the turntable data and the SINS sensor data are 100?Hz. The navigation solution frequency of SINS is 100?Hz as well.

Sensors’ zero bias, scale factor, coupling coefficient, installing error, and the variable of the flexible gyro related to can be calculated by the exact calibration [15]. A conclusion can be gotten that no above-mentioned errors exist and the table frame is overlapped with the sensor frame if the test is conducted followed by the calibration. In order to simulate the sensor errors caused by a long time of conservation and sensor unsuccessive turning-on error, the sensors’ errors are artificially added as shown in Table 5 after the calibration.

##### 5.2. Description of Test Process

The rotating period is set as 120s and the parameters for compass loop as those of in Section 3.3.2.

During the test, the inner frame and the intermediate frame of the turntable are kept in level states constantly. At , the outer frame rotates to 35° (this value can be selected randomly); then take the table frame at this time as the carrier body frame , which means the carrier theoretical attitude value is [0 0 35°]. At 14.6?s, the turntable begins to rotate by 10°/s^{2} angular rate acceleration to get the angular rate 3°/s. After 1.4?s, the angular rate is stable enough. Moreover, at the moment of 16 s is defined as and the table as frame .

During the test, the incremental angle information from turntable is adopted and the angular rate is considered as a constant value. Under such condition, analytical coarse alignment [7] is executed to get the DCM before the turntable starts to rotate afterwards; incremental angle of turntable from to is used to calculate DCM , and thus the initial DCM at can be calculated as . When compass alignment based on repeated calculation ends, fine DCM can be gotten and fine DCM can be calculated as??.

##### 5.3. Test Results and Discussion

Two compass alignment schemes are compared with data repeated calculation applied in Scheme 1 instead of Schema 2. Data from 0?s to 14.4?s are used for analytical coarse alignment in both schemes. In Scheme 1, real-time navigation solution and compass alignment are run from 16?s to 136?s, with the data in this period used for repeated calculation. After 136?s, the 19 cycles of repeated calculations and compass alignment are completed. In Scheme 2, real-time navigation solution and compass alignment are run from 16?s to 2416?s. However, the update times of both navigation solution and alignment operation are 240000. The alignment results come to the DCM and the true attitude [0 0 37.5287°] can be obtained from turntable output.

###### 5.3.1. Analysis of Sensor Data

The process of recording the turntable data and sensor data takes about 3700?s. Solid line and dotted line represent sensor data and projected data from frame to the frame both in the time of 16?s to 3516?s, respectively, in Figure 10. Therefore, with azimuth rotation, sensor constant errors along -axis and -axis are modulated into zeros mean periodic variables while sensor error along -axis is not modulated. Apart from this, the gyro output caused by turntable rotation is removed after projection from frame to , but gyro constant error still exists as presented in Figure 10(c).

###### 5.3.2. Alignment Accuracy

In accordance with alignment results in Figure 11, the enlarged portions are the attitude curves of the last repeated calculation cycle for Scheme 1 and of the last rotation cycle for Scheme 2, respectively.

In Figure 11, the true attitude and the limit attitude are represented by dotted lines and dot-dash lines. With the true value of DCM and sensor error in Table 5, equivalent northern and eastern accelerometer errors in frame are 0.0140?m/s^{2} and -0.0020?m/s^{2} respectively and equivalent eastern gyro is -0.1221°/s; thus the limit accuracy can that be gotten by classical compass alignment without rotation is ; in other words, the limit attitude is .

As analyzed in Section 2, the equivalent sensor errors are determined by sensor errors and the carrier’s attitude in SINS. On the above condition, the equivalent error of northern accelerometer is relatively small and that of eastern gyro is relatively large. Due to this, the accuracy of classical compass alignment method cannot objectively reflect the effects of sensor errors to the alignment accuracy. Assume that the carrier’s attitude is 0; then the limit accuracy obtained from classical compass alignment without rotation is , that is to say, limit attitude is . In Figure 11, this limit attitude is shown as solid lines.

The mean, variance, and peak-peak values from no. 229600 to no. 241600 of Scheme 1 and 2296?s to 2416?s of Scheme 2 are listed in Table 6. The statistical results indicate that the alignment accuracy is significantly improved in both schemes.

In Table 6, compared with that acquired from classical compass alignment method, which is , the horizontal alignment accuracy is improved by 6 times and the azimuth alignment accuracy is improved by 8 times, if the mean attitude of Scheme 1 is chosen as alignment result, which is .

In Table 6, compared with the limit alignment error, the horizontal alignment accuracy is improved by 1 time and the azimuth alignment accuracy is improved by 4 times, if the maximum/minimum values of the alignment curves of Scheme 1 is chosen as alignment result, which is .

With the fine alignment result of DCM and DCM , DCM can be calculated.

###### 5.3.3. Alignment Time

In this compass alignment, which is based on azimuth rotation and data repeated calculation, the alignment time is determined by the performance of the navigation computer after coarse alignment, the accelerated rotation, and data collection in one modulation cycle. As for a computer with 333?MHz frequency about 1?ms is required for a single navigation solution and compass alignment operation. To fulfill 22800 update operations (19 cycles of repeated calculations) 228?s is needed. The whole process for alignment takes 364?s (16?s + 120?s + 228?s) which is superior to the classical alignment time (600–800?s) as shown in Section 2.2.

In engineering, the conflict among the real-time sensor data collecting, real-time update for DCM in Section 4.2.3,and the repeated calculation can be resolved by multi-tasking real-time operating system (RTOS). In this type of SINS, VxWorks, an RTOS, is applied. When RTOS’ spending time for tasking scheduling and the time consumed for are concerning, the total consumed time for alignment is a little longer than 364 s and the measured time in the above SINS is 402?s.

#### 6. Conclusions

In this paper, in order to improve the alignment accuracy with certain accuracy sensors and eliminate the effect of carrier’s attitude to alignment accuracy, azimuth rotating modulation method is introduced to classical compass alignment for SINS.

Two methods are adopted to cope with the amplification effort in classical compass loop to low-frequency input, which does not exist without rotating modulation. The first method is to increase the rotating frequency and the second is to redesign the parameters of compass loop. The simulation results indicate that alignment accuracy is significantly increased through the above two methods; however, results turn out that alignment time is too long to be accepted.

To meet the requirement of rapid response, a fast compass alignment method based on data repeated calculation is designed. In this method, the same set of data is solved repeatedly to produce different DCM and different projected sensor data in mathematical platform from measured sensor data in body frame, thus producing different feedback signals. With the higher performance of navigation computer and less time for sensor data collecting, alignment time is significantly decreased.

Compared with the classical compass alignment, the results of turntable test can be summarized as follows: the horizontal alignment accuracy is improved by 6 times and the azimuth alignment accuracy is improved by 8 times when the mean of the alignment is taken as the evaluation criteria and the horizontal alignment accuracy is improved by 1 time and the azimuth alignment accuracy is improved by 4 times when the peak to peak value of the alignment curve is taken as the evaluation criteria. The alignment time measured by the computer with 333?MHz frequency is 402?s, which is slightly shortened in comparison with that in classical compass alignment.

#### Acknowledgments

This work was supported in part by the National Natural Science Foundation (61004125, 61273056) and Chinese university research and operation expenses (104.205.2.5).

#### References

- E. Levinson and C. S. Giovanni, “Laser gyro potential for long endurance marine navigation,” in
*Proceedings of the IEEE Position Location and Navigation Symposium*, pp. 115–129, Piscataway, NJ, USA, December 1980. View at Scopus - E. Levinson, J. ter Horst, and M. Willcocks, “The next generation marine inertial navigator is here now,” in
*Proceedings of the IEEE Position Location and Navigation Symposium*, pp. 121–127, Las Vegas, Nev, USA, April 1994. View at Scopus - J. I. Lahham and J. R. Brazell, “Acoustic noise reduction in the MK 49 ship's inertial navigation system (SINS),” in
*Proceedings of the IEEE Position Location and Navigation Symposium (PLANS '92)*, pp. 32–39, Charlottesville, Va, USA, March 1992. View at Scopus - H. N. Weng, Q. C. Lu, K. Huang et al., “Rotation scheme design for rotary optical gyro SINS,”
*Journal of Chinese Inertial Technology*, vol. 17, no. 3, pp. 8–14, 2009 (Chinese). View at Google Scholar - B. L. Yuan,
*Research on Rotating Inertial Navigation System with Four-Frequency Differential Laser Gyroscope*, Graduate School of National University of Defense Technology, 2007. - B. L. Yuan, D. Liao, and S. L. Han, “Error compensation of an optical gyro INS by multi-axis rotation,”
*Measurement Science & Technology*, vol. 23, no. 2, pp. 1–9, 2012. View at Google Scholar - Y. Yang and L. J. Miao, “Fiber-optic strapdown inertial system with sensing cluster continuous rotation,”
*IEEE Transactions on Aerospace and Electronic Systems*, vol. 40, no. 4, pp. 1173–1178, 2004. View at Publisher · View at Google Scholar · View at Scopus - Y. N. Gao, J. B. Chen, and T. P. Yang, “Error analysis of strapdown optic fiber gyro compass,”
*Transaction of Beijing Institute of Technology*, vol. 25, no. 5, pp. 423–426, 2005. View at Google Scholar · View at Scopus - J. C. Yu, J. B. Chen, and J. H. Han, “Multiposition observability analysis of strapdown inertial navigation system,”
*Transaction of Beijing Institute of Technology*, vol. 24, no. 2, pp. 150–153, 2004 (Chinese). View at Google Scholar - X. L. Wang and G. X. Sheng, “Fast and precision multiposition initial alignment method of inertial navigation system,”
*Journal of Astronautics*, vol. 23, no. 4, pp. 81–84, 2002 (Chinese). View at Google Scholar · View at Scopus - J. C. Fang and D. J. Wan, “A fast initial alignment method for strapdown inertial navigation system on stationary base,”
*IEEE Transactions on Aerospace and Electronic Systems*, vol. 32, no. 4, pp. 1501–1505, 1996. View at Publisher · View at Google Scholar · View at Scopus - H. T. David and L. W. John,
*Strapdown Inertial Navigation Technology*, Institute of Electrical Engineers, Hampshire, UK, 2nd edition, 2004. - A. Jamshaid and J. C. Fang, “Alignment of strapdown inertial navigation system: a literature survey spanned over the last 14 years,” ftp://labattmot.ele.ita.br/ele/alessandro/Leitura/04%20-%20PAPERS/ARTIGOS/KALMAN/kalmanlib/StrapdownAlignment.pdf.
- Y. Li, X. X. Xu, and B. X. Wu, “Gyrocompass self-alignment of SINS,”
*Journal of Chinese Inertial Technology*, vol. 16, no. 4, pp. 386–389, 2008 (Chinese). View at Google Scholar - X. Liu and X. Xu, “System calibration techniques for inertial measurement units,”
*Journal of Chinese Inertial Technology*, vol. 17, no. 5, pp. 568–571, 2009 (Chinese). View at Google Scholar