1. Introduction
The initial alignment (IA) accuracy of a launch vehicle (LV) directly affects the navigation accuracy, and the initial alignment time directly affects the launch vehicle test and launch process. Due to the influence of disturbance factors, such as wind and waves at sea [
1], offshore launches of launch vehicles are subject to more severe swaying than land-based launches. It is necessary to solve the self-alignment problem of the launch vehicle under motion and swaying conditions and achieve considerable accuracy in a very short launch period. The alignment is one of the most important stages of a navigation system [
2]. The optical sighting or self-sighting algorithm used for initial alignment of land-based launch vehicles cannot be directly applied to the initial alignment at sea [
3]. In case of land-based launch, the optical sighting equipment is fixed and the launch vehicle is slightly swayed. If the launch is on an offshore launching platform, both the optical aiming equipment and the launch vehicle are swaying at the same time, it will be impossible to calibrate. How to propose a high-precision initial alignment method suitable for sea launch in response to the characteristics of the marine motion environment is one of the technical difficulties for the launch of launch vehicles [
3].
At present, optical sighting is usually used at home and abroad for azimuth reference transfer, and filter compensation or traditional transfer alignment method is used to complete the horizontal reference transfer [
3]. The fundamental principle of the common method for initial alignment of sea launches of launch vehicles is shown in
Figure 1. The principle of this method is simple and reliable, but the application cost is high, and the accuracy is slightly worse.
Conventional self-alignment methods are performed sequentially in chronological order, and data processing is often one-off; the data of the previous stage are not used in the later stage, and there is no storage requirement, resulting in inadequate information mining [
4]. In contrast, the computing power of current navigation computers is strong enough and the capacity is large enough to save the SINS data, which can then be repeatedly processed and analyzed, providing another method of improving the accuracy and speed of an LV’s self-alignment at sea.
The coarse alignment (CA) process mainly accomplishes the coarse estimation of the unknown attitude angle [
5]. Qin et al. [
6] proposed a coarse alignment algorithm based on gravitational acceleration, which effectively isolated the angular swaying interference by introducing the assumption of a solidification coordinate system (CS) and tracking the attitude change in real time according to the gyro output. Bai et al. [
7] demonstrated the effectiveness of the gravitational acceleration-based alignment algorithm through verification experiments. The measured angular velocity of the earth’s rotation and gravitational acceleration can be used in the coarse alignment stage. Swaying is a perturbation situation that mainly affects orientation changes [
8]. With the help of the idea of decomposing the inertial coordinate system, the gravity vectors in the inertial space at different moments are not co-linear, which can be solved by the double-vector attitude determination [
9,
10] or multi-vector attitude determination algorithm, which isolates the swaying interference [
11,
12,
13]. To improve the accuracy and efficiency of latitude and attitude estimation by fully utilizing inertial measurements, an anti-swaying self-alignment method based on the inertia tensor of gravity observation vector was presented in Reference [
2].
Coarse alignment results can be improved by applying fine alignment (FA) at the completion of the coarse alignment phase. Fine alignment improves the accuracy of the fusion process by using external sensors or information [
14,
15]. Tsukerman et al. [
16] derived an analytic assessment of the steady-state characteristics of the fine alignment procedure.
To improve the alignment performance by taking full advantage of all sensor data in a limited amount of time [
17,
18,
19], a real-time information-reusing alignment strategy was proposed for rotating an inertial navigation system in Reference [
20]. Li et al. [
21] designed a backtracking navigation scheme for INS (Inertial Navigation System)/DVL (Doppler Velocity Log) integration. Yan et al. [
22] proposed an inter-travel alignment scheme based on reverse data processing, which organically combined the initial alignment of the moving base with the position navigation through the reverse navigation algorithm, which first carried out the initial alignment of the moving base attitude and the sensor sampling data storage of the SINS (Strapdown Inertial Navigation System) gyrocompass and then used the reverse and forward dead reckoning algorithms to obtain the position information, which solved the previous problem that the position needed to be corrected using landmark points after inter-travel alignment. Liang et al. [
23] applied the reverse navigation algorithm to the static base alignment of the SINS, which effectively shortened the alignment time and improved the alignment accuracy. Wang et al. [
24] introduced the reverse Kalman filtering technique, which processed the collected historical data several times using the forward and reverse Kalman filtering algorithms and completed the reverse fine alignment and forward combined navigation on the basis of the coarse alignment results. In Reference [
25], aiming at the problem of long alignment time for traditional methods, a forward-forward backtracking initial alignment method under swaying base conditions was proposed to achieve the purpose of rapid initial alignment, which did not consider linear vibration.
The incremental output SINS is taken as the research object. Enlightened by [
19], a self-alignment strategy based on the reverse data processing technology is proposed, which uses pure mathematical analysis to complete the alignment of azimuth and horizontal directions simultaneously without the need for complex optical transmission equipment, thus achieving the purpose of shortening the alignment time and improving the alignment accuracy. The main innovations of this paper are as follows:
Reverse Kalman filtering algorithm is designed, and the backtracking navigation is completed by the combination of the reverse navigation solution and reverse Kalman filtering;
A cycle-index control function is designed for approximate strict backtracking navigation;
In addition to the attitude correction at the initial moment (IM), attitude error compensation is completed twice during the subsequent alignment.
Compared to the existing alignment methods, firstly, using of the gravity projection along the inertial system, with reverse Kalman filter technique, could avoid the effects of disturbed accelerates and angular velocities on the SINS alignment. Secondly, attitude error compensation, to a certain extent, can achieve real-time tracking, which is the most obvious difference from the existing methods. Then, the designed cycle-index control function enables the control of the number of forward and reverse navigation solutions, thus controlling the amount of computation. Simulations and experiments show that the strategy proposed is effective. In addition, the strategy can be applied to engineering practice, as it effectively improves alignment accuracy and reduces alignment time compared with the traditional alignment methods.
2. Problem Description
Coordinate systems and Attitude Angles are described as follows:
Geocentric Inertial CS ( system): . The origin of the CS is located at the center () of the earth. The -axis points to the vernal equinox, and the -axis is along the earth’s axis of rotation. The -axis and the -axis and -axis form a right-handed CS.
Earth CS ( system): . The CS is fixed to the earth, and the origin is at the center () of the earth. The -axis passes through the intersection of the prime meridian and the equator, the -axis passes through the north pole and the -axis passes through the intersection of the meridian 90°E and the equator.
Geographic CS ( system): . The origin of the CS is located at the center of mass of the carrier, and the -axis, -axis and -axis point to the east, north and up, respectively.
Navigation CS ( system): . The CS is the reference CS for the navigation solution, and the geographic CS is used as the navigation CS in this paper.
Calculation of Navigation CS ( system): . The CS is the navigation CS established according to the navigation parameters obtained from the inertial navigation solution.
Inertial Measurement Unit (IMU) CS ( system): . The origin of the CS is located at the center of mass of the IMU, and the three axes point to the direction of the IMU sensitive axis. It is assumed that the body CS is rigidly connected to the IMU CS.
IMU CS at the IM ( system) coincides with the IMU CS ( system) at the instance of the IA, and there is no rotation relative to inertial space.
Navigation Inertial CS at the IM ( system) coincides with the navigation CS ( system) at the instance of the IA, and there is no rotation relative to inertial space.
Launch CS ( system): . The origin of the CS is located at the center of the launch point, the -axis is horizontally forward along the ballistic direction, the -axis is vertically upward, the -axis is perpendicular to the ballistic plane to the right and the -axis and -axis form a right-handed CS. Obviously the -plane is the ballistic plane. The angle between the ballistic plane (or -axis) and the local geographic north direction is usually called the launch angle, which is recorded as .
The Inertial Navigation System (IMU) is installed on the LV, and its axial definition is shown in
Figure 2. When the LV is “lying down” horizontally, the three axes are along the lateral, longitudinal and vertical directions respectively, i.e., “right-front-up”. Euler angles are defined as “132”:
Pitch angle (): After a coordinate translation so that the origin of the navigation CS coincides with the origin of the IMU CS, the pitch angle is the angle between the projection line of the longitudinal axis (-axis) of the LV in the ballistic plane and the -axis. The angle is positive if the longitudinal axis of the offshore platform points above the horizontal plane; otherwise, it is negative, and the value range is .
Yaw angle (): After a coordinate translation so that the origin of the navigation CS coincides with the origin of the IMU CS, the yaw angle is the angle between the longitudinal axis of the LV and the ballistic plane. Looking at the angle plane, if the rotation from the -axis to the -axis is counterclockwise, the angle is positive; otherwise, it is negative, and the value range is .
Roll angle (): After a coordinate translation so that the origin of the navigation CS coincides with the origin of the IMU CS, the roll angle is the angle between the vertical axis of the LV and the plumb plane containing the longitudinal axis of the LV. Looking forward along the longitudinal axis from the tail of the LV, if the -axis is on the right side of the vertical plane, the angle is positive; otherwise, it is negative. The value range is .
The conversion relationship applied throughout the article is described as follows: is the representation of the attitude of B system with respect to N system. It is equivalent to say that it is the transformation matrix from N system to B system.
The differential equation for the specific force of SINS in a launch vehicle can be described as:
where
is the ideal output real-time specific force of the accelerometer;
represents the real-time angular rate of rotation of the navigation system due to Earth’s rotation;
represents the
n system rotation of the SINS moving near the Earth’s surface due to the curvature of the Earth’s surface;
represents the gravity vector in the
n system [
27].
Assuming the offshore platform is moored. Influenced by the sea state, the motion state of the launching platform can be decomposed into swaying motion and linear motion. Under such conditions, the initialization of the SINS in a launch vehicle will be obtained by the self-alignment strategy designed later.
3. Design of Self-Alignment Method
To realize the rapid self-alignment process of the LV launched on the offshore platform, this research established a backtracking system error model, and a reverse Kalman filtering algorithm was designed. The self-alignment process is divided into three parts: ASCA, backtracking navigation FA and forward FA.
3.1. Self-Alignment Strategy of LV
The designed self-alignment strategy of the LV launched on the offshore platform is shown in
Figure 3. In
Figure 3,
represents the initial attitude direction cosine matrix calculated at the end of the CA.
represents the initial attitude direction cosine matrix calculated at the end of the backtracking FA.
represents the initial attitude direction cosine matrix obtained from each backtracking FA estimation.
Firstly, at the IM of the alignment, the CA process is started by the ASCA method while storing IMU data in the mooring condition. The CA process is completed at the moment . The computed attitude matrix , which is obtained by the CA, is used as the initial value of the first backtracking FA, and the backtracking FA uses the collected IMU data information to return the state of the LV to the IM . At the same time, the attitude information estimated by the BN is combined with the known velocity information at the starting point to correct the inertial navigation. The first backtracking FA process is completed, and the attitude information with a certain accuracy is obtained. Immediately after the fast forward navigation to return the system state to the moments, the attitude errors are compensated using the data stored from to moment. Then, the backtracking-forward cyclic FA starts from moment . When the last forward FA returns to moment , the data stored after moment are used to catch up with the normal real-time computation until the end of the alignment. The navigation computer maintains high speed computing during the whole process. The entire FA process is completed, and higher accuracy attitude information can be obtained.
The BN algorithm is formally identical to the forward navigation algorithm. As long as the forward saved gyro data and the angular rate of Earth rotation, etc., are inverted and the forward final value is set to the backtracking initial value, BN can be realized. During the backtracking solution, all sensor data must be stored. Therefore, in the IA, if the data processed by the backtracking continues to be used as a new set of data, or if the forward and backtracking attitude solutions are combined and the data are used repeatedly, this is equivalent to increasing the quantity of data.
3.2. Anti-Swaying Coarse Alignment Algorithm
Through the idea of decomposition of inertial solidification CS, the initial attitude matrix
is decomposed into
represents the transformation matrix from the navigation inertial CS (-system) at the IM to the navigation CS (-system); represents the transformation matrix from the IMU CS (-system) at the IM to the navigation inertial CS (-system) at the IM; represents the transformation matrix from the IMU CS (-system) to the IMU CS ( system) at the IM.
According to the angular velocity of the earth’s rotation and the time interval between the two moments, the relationship between the two moments of the navigation CS can be obtained; according to the projection of the gravitational acceleration vector at the two moments in the navigation CS at the IM of alignment, and the projection of the specific force measurement value in the initial IMU CS, the conversion relationship between the IMU CS and the navigation CS at the IM can be solved by the double-vector attitude determination algorithm. From the projection and conversion relationship represented by the above three matrices, the conversion relationship between the IMU CS and the navigation CS at the end of the alignment can be determined, and then the CA is completed.
3.3. Backtracking Navigation Algorithm
In order to realize the sampling and utilization of the output signals of the gyroscope and accelerometer without an omission in time, the system collects the angular incremental output of the gyroscope and the velocity incremental output of the accelerometer within a defined time interval, rather than the angular velocity output and the specific force output at discrete time points. The attitude, velocity and position update equations of the incremental SINS inertial guide can be generally expressed as follows.
3.3.1. Forward Navigation [27]
where
and
represent the SINS attitude matrix at
and
(
), respectively. For the navigation update period
,
.
denotes the midpoint moment.
represents the rotation of the
n system with respect to the
i system at the moment of the midpoint in the forward navigation. The subscript
indicates that they are the variables at the midpoint moment.
denotes the angular increment in the time interval
, and
denotes the angular increment in the time interval
. The equivalent rotation vector can be obtained using the two-subsample conical error compensation algorithm:
can be obtained from the equivalent rotation vector
(
):
where
and
are the specific force velocity increment of the navigation system and the velocity increment of the detrimental acceleration in the time period
, respectively.
is the rotational error compensation amount of the speed, which is caused by the rotation change of the specific force direction in space during the solution time period;
is the paddling error compensation amount.
where
,
,
denotes the velocity increment in the time interval
, and
denotes the velocity increment in the time interval
.
where
,
can be obtained by linear extrapolation [
28].
3.3.2. Backtracking Navigation
To obtain the inertial navigation attitude, velocity and position update equations for the backtracking solution, Equations (3)–(7) are transformed by shifting the terms, by taking the reverse of the angular velocity of ground rotation, the angular velocity of rotation of the navigation CS relative to the earth CS, the acceleration of gravity, the angular increment of the gyro output and the velocity increment of the angular velocity meter output to obtain:
- (1)
Attitude update during backtracking
where
and
represent the SINS attitude matrix at
and
(
), respectively. For the backtracking navigation update period
,
.
denotes the midpoint moment.
represents the rotation of the
n system with respect to the
i system at the moment of the midpoint in the backtracking navigation. The subscript
indicates that they are the variables at the midpoint moment.
denotes the angular increment in the time interval
, and
denotes the angular increment in the time interval
.
- (2)
Velocity update during backtracking
where
and
are the specific force velocity increment of the backtracking navigation system and the velocity increment of the detrimental acceleration in the time period
, respectively.
is the rotational error compensation amount of the speed, which is caused by the rotation change of the specific force direction in space during the solution time period;
is the paddling error compensation amount in the backtracking navigation system.
where
,
,
denotes the velocity increment in the time interval
, and
denotes the velocity increment in the time interval
.
- (3)
Position update during backtracking
Through the above derivation, it can be seen that the attitude, velocity and position can at the same time be basically consistent during the forward and backward navigation process, the same piece of navigation data can be solved repeatedly and the known information at the starting point can be used for correction to improve the alignment accuracy.
3.4. Reverse Kalman Filtering
The error equation of the BN system is deduced according to the inertial navigation error equation [
25].
The actual calculated attitude matrix
has a small angle of attitude error angle with the ideal attitude matrix
, then
can be expressed as
The calculated attitude matrix
is generally solved by the following differential equation:
From the above, it is clear that the BN process requires the inversion of the gyro output angular velocity
and the rotation angular velocity
of the navigation system, i.e.,
Substituting the differentiation of Equations (15), (16) and (13) into Equation (14), we get
Multiplying both sides of Equation (17) simultaneously right by
and omitting the second-order small quantity,
The inertial navigation attitude error equation for BN is obtained by simplifying Equation (18),
According to the velocity update equation for BN, the differential equation of inertial navigation velocity for BN can be expressed as
The actual calculated velocity
is obtained by Equation (21).
where
,
,
,
.
The velocity error
can be obtained by subtracting Equations (21) and (20)
Substituting Equation (13) into Equation (22) and expanding it gives:
The Equation (23) is simplified by omitting the second-order small quantity error and the gravity error
, and the velocity error equation for the BN is obtained as follows:
The position update differential equation for BN can be written as:
The position error equation for the BN is obtained from the position update differential equation Equations (25)–(27) as follows:
Integrating the error equation of the BN derived above, the corresponding systematic error equation of the BN system is as follows:
Define:
where
,
and
represent the attitude error, velocity error and position error of the SINS in the three directions of East, North and Up, respectively;
and
are the gyroscopic drift and the accelerometer zero error, respectively.
According to Equation (31), the reverse Kalman filtering model of the navigation system can be deduced as:
where
is a 14-dimensional state variable;
is the observation vector;
and
represent the system matrix and observation matrix, respectively;
is the noise allocation matrix;
and
represent the system noise and observation noise, respectively.
System matrix in reverse Kalman filtering:
Observation vector:
where
denotes the velocity component obtained from the SINS solution.
According to the above equations, more accurate attitude information can be obtained by retrospectively solving the system state through the backtracking Kalman filtering algorithm.
Remark 1. The continuous-time reverse Kalman filter model derived above is discretized, and then the corresponding terms are replaced byand.
For a launch vehicle, the changes of
and
are small in a short period of time
, so it is generally considered that the quadratic functions of
and
are slow quantities of time. The calculation error caused is generally small, and the relatively simple gradient integration method can be used to obtain
where
denotes the output angular velocity of the gyroscope at the moment of the midpoint
;
denotes the output specific force of the accelerometer at the moment of the midpoint
.
3.5. Cycle-Index Control Function
For incremental SINS, repeated forward and backtracking analysis to calculate the stored data is beneficial to improving the navigation accuracy. In this process, in order to simplify the calculation, some second-order small quantities are often omitted, or small angle approximation is adopted, so that the error of each iteration is accumulated. When the number of forward and backtracking is small, the effect of error on the results is not obvious. When the number of forward and backtracking is large and reaches a certain number of times, the error is amplified to the same multiple as the number of times due to repeated calculations, resulting in oscillation convergence or even the drift of the result not converging.
For the case where the swaying amplitude is small, the backtracking process approximates the strict reverse derivation process, and the solution is more accurate. Multiple iterations of forward and backtracking solutions can guarantee the accuracy and avoid the accumulation of calculation errors. At this point, the cycle-index control function is established as follows:
where the whole total sampling time is recorded as
, the first backtracking time is recorded as
, and
is the frequency control coefficient.
3.6. Attitude Error Compensation
During the self-alignment process of the LV launched on the offshore platform, when returning to the starting point of the BN, the velocity and attitude angle change due to the influence of waves, etc., so attitude error compensation is required.
The first attitude error compensation: after obtaining the relatively attitude matrix at moment , the attitude update is solved with its initial value using the SINS data stored after moment . It is updated until the latest sampled data moment , thereby obtaining the attitude matrix estimated at moment .
The second attitude error compensation: the cyclic solution starts at the moment , while the sampled data after are stored. After the last forward navigation to obtain the exact matrix at time , the SINS data stored after moment are used as the initial value for attitude update solving. Until the IA ends, the latest attitude matrix is obtained.
4. Experiments and Results Analysis
Remark 2. The self-alignment of the LV launched on OP is obviously affected by the line vibration. The oscillation in the convergence section of the estimated curve reflects the characteristics of the motion and swaying of the LV. For the influence of line vibration disturbance, considering that the period of the gravitational acceleration vector is very small with inertia, low-pass filtering and other methods can be used, but this is not within the scope of this study.
Remark 3. The orientation attitude of the LV’s self-alignment is the roll angle rather than the yaw angle.
Remark 4. The form of attitude matrix in the table is: [pitch, yaw, roll].
The simulation process simulates the self-alignment process of the LV in the mooring condition. In the Xisha and towards the equatorial sea area, the swaying of the LV and the launch platform disturbed by the waves is approximately sinusoidal with equal amplitude oscillation in sea state level 4, the swaying amplitude does not exceed and the swing period is .
The hardware environment for simulation was an Intel(R) Core(TM) i5-8400 CPU with 2.80 GHz, 8.00 GB RAM and Windows 10 operating system.
The inertial sensor error parameters [
29] were set as shown in
Table 1. The data output rate of the inertial sensor was set to 200 Hz.
The swaying parameters [
30] are shown in
Table 2.
In the mooring condition, in addition to swaying, there are also lateral, longitudinal and vertical linear velocities. The linear velocity parameters [
30] are shown in
Table 3.
The installation height of the SINS in the LV does not exceed 60 m. Combining the swaying of the offshore launch platform with the linear velocity, the trajectory of the LV in the relative coordinate system is shown in
Figure 4 and it is bounded by the scope. It has a displacement of approximately one meter, which is characteristic of the engineering application environment.
The simulation scheme is as follows:
Scheme A: The conventional method uses the first 300 s of data for ASCA, and then uses the next 360 s of data for conventional Kalman filtering FA.
Scheme B: ASCA uses the first 300 s of data, the first backtracking FA uses the first 300 s of data, the first forward FA uses the first 330 s of data, repeated backtracking and forward FA uses the first 330 s of data; and the final forward FA uses the first 360 s of data.
The time taken for scheme A and B is shown in
Figure 5.
The time taken for scheme A is 660 s, while Scheme B takes 360 s. Therefore, the proposed scheme shortens the time by 45.5% compared with Scheme A.
4.1. Estimation of the Misalignment Angle at the First Cycle
Using the strategy proposed in this paper, a stochastic simulation was performed.
Figure 6 shows the estimated curves of the misalignment angles during the first backtracking navigation FA and the first forward navigation FA. The arrow in
Figure 6 indicates the estimated direction.
After using the reverse Kalman filtering back to the IM and correcting the SINS with the estimated results, most of the attitude misalignment angles are corrected, and the filter curve converges during the forward alignment process, which can guarantee the alignment accuracy.
4.2. Analysis of the Last Forward Navigation Alignment Results
The CA with an alignment time of 300 s and the FA with an alignment time of 360 s were studied to simulate the conventional self-alignment based on scheme A and the self-alignment strategy based on scheme B, and the results of attitude error obtained for scheme A and scheme B are shown in
Figure 7.
A stochastic simulation was performed. In the conventional approach, the roll angle error enters the error band after completing the alignment, while the roll angle error enters the error band in the proposed method. Since the proposed method utilizes the SINS data from the CA stage in the FA, more information is obtained at the same time. The alignment time of the proposed strategy depends on the convergence time of the FA filter, and the CA no longer occupies the alignment time; the SINS data of the CA stage are used twice with less storage and calculation, which not only increases the quantity of data but also effectively shortens the alignment time.
The convergence segment of the simulation results from 300 s to 360 s was selected for analysis. The statistics in
Table 4 are for the convergence segment.
As can be seen from
Table 4, the convergence quality of the proposed strategy is significantly better than that of the conventional method. The proposed strategy makes use of the SINS data in the CA stage from the beginning and obtains more information in a limited time with better results; with the increase in the frequency of backtracking, there is no obvious adjustment in the alignment process of the method. Through analysis, the algorithm can achieve accuracy after one or two backtracking processes.
4.3. Simulation Analysis of Multiple Reverse Navigation
Under the above simulation conditions, four of the reverse processes were selected and the alignment error curves are shown in
Figure 8. Clearly, the orientation attitude angle is continuously corrected as the number of reversals increases.
4.4. Monte Carlo Simulation Results and Discussion
The statistics of the results of 50 simulations under the alignment conditions of scheme A and scheme B are shown in
Figure 9.
The statistical characterization of the converged three-axis attitude angle errors was carried out, and the corresponding characterization is shown in
Table 5.
As can be seen from
Table 5, the alignment accuracy of scheme B is obviously better than that of scheme A, its attitude angle error is less than 3.9′ and the average value is within 0.82′. The standard deviations of pitch error, roll error and yaw error are 0.0853′, 1.4095′ and 0.0123′, respectively. It can be seen that the stability of the convergence section of the attitude angle error in the horizontal direction is improved, and the stability of the convergence section of the roll angle error decreases.
4.5. The Effect of CA Time on Alignment Accuracy of the Proposed Method
At present, there are few studies on the BN that focus on the effect of CA time in alignment accuracy. Based on the proposed method in this paper, 50 Monte Carlo simulations were conducted using the CA time of 120, 240, 300 and 330 s. The statistical analysis of the alignment accuracy results is shown in
Table 6.
In general, the ASCA algorithm has a relatively weak anti-interference ability within a short period of time, and if the CA time is extended appropriately, good CA results can often be achieved. However, it can be seen from
Table 6 that appropriately extending the CA time and keeping the FA time unchanged does not necessarily improve the final alignment accuracy, which is related to the accumulation of errors during the backtracking.
4.6. Validation Test
To compare the proposed strategy with conventional optical aiming method, the laser IMU is mounted on a six-degree-of-freedom turntable, and the relationship between the IMU orientation and the frame angle of the turntable is determined using the prism on the IMU to establish the reference orientation. The update frequency of IMU’s raw data is set to 100 Hz. The bias stability of gyroscopes and accelerometers are better than 0.01°/h and 20 µg, respectively. Following the laboratory confidentiality requirements, we will not show the physical picture.
Different attitude commands were executed on the turntable to simulate the shaking conditions during the launch of an LV, and the alignment results obtained are shown in
Table 7.
As can be seen from
Table 7, the experimental results under various conditions are basically consistent. The proposed alignment strategy can track the sway of the SINS well, and the alignment accuracy is within 3′.