Research Article  Open Access
Hao Liang, Yumin Tao, Meijiao Wang, Yu Guo, Xingfa Zhao, "SystemLevel Temperature Compensation Method for the RLGIMU Based on HHORVR", Journal of Sensors, vol. 2021, Article ID 6613574, 16 pages, 2021. https://doi.org/10.1155/2021/6613574
SystemLevel Temperature Compensation Method for the RLGIMU Based on HHORVR
Abstract
The ring laser gyro inertial measurement unit has many systematic error terms and influences each other. These error terms show a complex nonlinear drift that cannot be ignored when the temperature changes, which seriously affects the stability time and output accuracy of the system. In this paper, a systemlevel temperature modeling and compensation method is proposed based on the relevance vector regression method. First, all temperaturerelated parameters are modeled; meanwhile, the Harris hawks optimization algorithm is used to optimize each model parameter. Then, the system compensation is modeled to stabilize the system output to the desired temperature. Compared with the least square method, the fitting performance comparison and the system dynamic compensation experiment prove this method’s superiority. The root mean square error, the mean absolute error, the squared, and the variance of residual increased by an average of 35.27%, 39.29%, 2.29%, and 30.34%, respectively.
1. Introduction
In the inertial sensors, the ring laser gyroscope (RLG) [1–3] has many advantages such as small random drift, wide dynamic range, fast startup speed, and high reliability. It has a high value in the field of inertial navigation. As a representative of strapdown inertial navigation, the ring laser gyroscope inertial measurement unit (RLGIMU) [4–6], which uses RLG as its core measurement component, is widely used in the field of inertial navigation for its outstanding advantages of small size, high precision, strong environmental adaptability, and low cost. How to improve the performance of the RLG and the RLGIMU has always been the focus of researchers. After years of research, researchers have done much work on various physical fields that affect the system performance, such as temperature, magnetic field, vibration environment, and how to improve the performance of the RLGIMU under the influence of various physical fields.
Temperature is an important factor affecting the output precision and startup stability time of the RLGIMU. On the one hand, the operating temperature range of the RLGIMU is extensive, which requires strong adaptability of the system. On the other hand, many RLGIMU devices have a temperature control system composed of the heating part, temperature feedback part, and temperature control part. At the same time, the RLGIMU has many circuit components, which also generate heat. Other components inside the IMU, such as magnetic coils, also generate heat. After the RLGIMU is started, a large amount of heat will be generated inside the system. The internal temperature changes, and IMU’s output is unstable due to the internal heating part. Such changes in the system’s internal and external temperature affect the RLGs and quartz flexible accelerometers’ performance and the RLGIMU internal structure’s characteristics. The change of structural characteristics leads to the thermal deformation of inertial devices, which causes a lot of fluctuation and error in output data of the RLGIMU and limits the startup time and performance of the RLGIMU. Therefore, it is necessary to study the temperature effect of the IMU system to improve its temperature adaptability and reduce the system stability time.
At present, there are two methods to solve the influence of temperature on the IMU system: temperature control and temperature compensation. The method of temperature control makes the system quickly heat up to reach the temperature control point, to make the system stable as soon as possible. After the RLGIMU is started, the internal temperature control system starts to work. The system temperature generally reaches a relatively stable state after a relatively long time (generally dozens of minutes). At this time, the output of the RLGIMU is stable, and the output accuracy can meet the requirements. Before that, the output of the IMU was unstable due to the drastic temperature change. However, this method’s problem is that temperature control requires a “long setting time” or called a “long startup time.” The accuracy of the temperature control system will also affect the IMU’s performance, and additional hardware needs to be added, and the system’s power consumption increased. The method of temperature compensation only needs to obtain the error model under the temperature, and the compensation can be started at the startup of the IMU. Theoretically, the system’s output can be stable as soon as it is powered on, but the accuracy of the model establishment will affect the accuracy of the compensation. How to reduce the startup stable time by using the system temperature compensation method is the purpose of this paper. By modeling the system measurement equation’s parameters under different temperatures, the output of the RLGIMU at different temperatures can be converted to the output at a stable temperature. In this way, the operation of the RLGIMU will not be affected by the change of temperature so that IMU can meet the requirements without “long setting time.”
Several methods have been used for the RLG’s or quartz flexible accelerometer’s zero position temperature compensation, such as stepwise regression [7, 8], artificial neural network [9–12], support vector machine [13–15], and mean [16]. The IMU system’s various error terms are very complicated, including the zero position and the influence of the scale factor. The system structure changes will also affect the zero position and scale factor of the sensors. These factors not only are affected by temperature but also have a coupling relationship with each other. Compensating for the RLG or quartz flexible accelerometer zero position alone, without considering the effects of other error terms, cannot cover the system’s error, which will affect the system’s compensation accuracy. Therefore, all drifts must be compensated at the system level. At present, the least square method (LSM) is usually used for systemlevel compensation in engineering [17–20]. The IMU has more than a dozen parameter items related to temperature. When these parameters change with temperature in a sophisticated nonlinear manner, this traditional modeling method’s fitting accuracy is limited.
Machine learning has powerful capabilities in the regression and prediction of complex functions. Based on the Bayesian framework, Tipping proposed the relevance vector machine [21–23]. Compared with the support vector machine, the relevance vector machine uses fewer vectors and has stronger sparsity. Although the training time is longer than that of the support vector machine, the prediction time is much less than that of the support vector machine. The kernel function does not need to meet Mercer’s condition: in a finite input space, the function is a map. If the kernel matrix is positive semidefinite, then the function can be a kernel function. In SVM theory, the kernel function must satisfy Mercer’s condition. In relevance vector machine theory, because of the difference between the relevance vector machine and the SVM architecture, the kernel function does not have to satisfy Mercer’s condition, so more kernel functions can be selected. In the case of fewer training data samples, it can ensure excellent generalization ability. The method of using the relevance vector machine for regression is called relevance vector regression (RVR) [24, 25].
In order to solve the problem of highprecision compensation of the RLGIMU systemlevel temperature, this paper proposes a systemlevel temperature error model and compensation method for the RLGIMU based on the RVR. According to the IMU system’s inputoutput model, all the parameters that affect the output are modeled at the system level and compensated so that the systemlevel temperature error compensation is more comprehensively achieved. Since the setting of the relevance vector machine kernel function’s width parameter has a severe impact on the regression accuracy, it is necessary to optimize this parameter. The Harris hawks optimization (HHO) [26] is a novel metaheuristic [27, 28] optimization algorithm. Compared with the genetic algorithm [29] and particle swarm optimization algorithm [30], it has fast optimization speed and high precision. In this paper, the HHO algorithm is used to optimize the kernel width parameter in the relevance vector machine to improve the regression accuracy of the model, so the method is called HHORVR.
For the RLGIMU system, the influence of temperature on the RLGIMU system is multifaceted. Although the influence of the temperature on the RLGIMU is multifaceted, all the effects will be reflected in the measured output pulse at different temperatures. According to the pulseangular velocity equation and the pulseapparent acceleration equation, if the objective angular velocity and apparent acceleration are considered real and do not change with temperature, the pulse change caused by temperature is caused by other equation parameters with the change of temperature. This paper is aimed at obtaining the fitting model of the parameters in the equation varying with temperature. When the system works at all temperatures, the system’s output pulse can be converted into a pulse output at a selected temperature. Thus, the output of the system is more stable when the temperature changes. This compensation method does not need to consider every factor affected by temperature change in the system. Therefore, we call this compensation method to be systemlevel temperature compensation.
This paper’s structure is as follows: Section 1 introduces the background, current problems, and this paper’s work. Section 2 analyzes the influence of temperature on the IMU system. Section 3 introduces the relevance vector machine regression theory and the HHO optimization process and establishes the system parameters’ temperature and compensation models. Section 4 introduces the experimental methods, experimental results, and analyses to verify the method’s effectiveness and superiority. Section 5 summarizes the whole paper.
2. Analysis of Temperature Effects
As shown in Figure 1, a typical RLGIMU system consists of three orthogonal RLGs and three orthogonal quartz flexible accelerometers mounted on the base. The IMU measures angular velocity in three directions with three RLGs (, , and ) and apparent acceleration in three directions with three quartz flexible accelerometers (, , and ).
The influence of temperature on IMU mainly includes three aspects: (1)Impact on the RLG
The laser gyroscope is based on the Sagnac effect principle and uses the optical path difference to measure the rotational angular velocity. The ring laser gyro is essentially an active ring laser. It is a laser source filled with a heliumneon mixture. The effect of temperature on the gyro is comprehensive [31], such as the change of material characteristics of the laser gyroscope, the influence on the length of the resonant cavity and its coplanarity, the change of the gas flow rate, the change of the characteristics of the mirror, the deformation of the capillary, changes of lockin characteristics, Langmuir flow, and discharge symmetry. These changes are reflected in changes in the scale factor and zero bias [32]. (2)Impact on the quartz flexible accelerometer
When the quartz flexible accelerometer is started, the fluctuation of the internal temperature field will cause the differential capacitance sensor to generate a larger current, further increasing the moment coil’s temperature. The increase in temperature affects the magnetic materials and torque coils in the accelerometer, causing changes in the magnetic flux and causing the accelerometer scale factor’s temperature drift. Simultaneously, due to the thermal imbalance in the quartz flexible accelerometer, the arm will produce slight distortion, which will affect the stability of the accelerometer zero bias. Besides, changes in the arm length can affect the scale factor of the accelerometer. (3)Impact on IMU system structure
For the IMU system, the three RLGs and three quartz flexible accelerometers are fixed on the base and cannot be strictly orthogonal. They have relative installation errors. Due to the generally strong rigidity of the base, the installation errors are relatively stable. However, the geometry of the base is also a small change that will occur with the temperature change, and this change will also be reflected in the parameter drift of the single sensor.
From the above analysis, it can be concluded that the influence of temperature on the RLGIMU system is multifaceted and complex. Temperature affects not only the individual sensor but also the system structure. Structural changes are also reflected in the output of the sensor. It is not straightforward to analyze and compensate for these factors separately. Therefore, this paper uses the systemlevel compensation method to integrate all the influencing factors to model and compensate directly.
3. Relevance Vector Regression and Temperature Compensation Modeling
3.1. Function Regression Theory of the Relevance Vector Machine
The relevance vector machine is based on the Bayesian theory. In linear regression problems, the input and output can be described as , where is the weight vector, is the input vector, and is the offset. When the relationship is nonlinear, it can be expressed as or . Given a set of input target data sets: , is the number of data samples. Consider that the target set is data samples superimposed with noise, expressed as , where is the noise of the Gaussian distribution satisfying the mean of 0 and the variance of , that is, ; the conditional probability of the target is [33–35]
is independent of each other; then, the likelihood functions of all corresponding data sets are expressed as where , , and is the matrix with , wherein . is the kernel function.
In order to establish the relationship model between input and output, the weight in the model is given prior probability.
Hyperparameter is used to describe the inverse variance of each .
The posterior probability of the weight is where the covariance and mean are , , and .
Use the maximum likelihood method to obtain the optimal solution of .
For new inputs , the predicted distribution can be calculated.
Since both of the integrands are Gaussian, this can be calculated.
In equation (2), the kernel function is used to map the feature vector to the highdimensional space, which can reduce the calculation complexity. Commonly used are the linear kernel function, polynomial kernel function, sigmoid kernel function, and radial basis kernel function. The linear kernel function is mainly used in the case of linear separability. The application effect of the polynomial kernel function on nonlinear data is not ideal, and there are many parameters to be adjusted, which increases the complexity of the model. The sigmoid kernel function has good performance in dealing with nonlinear data, and two parameters need to be adjusted. The radial basis function is a kind of kernel function with robust localization, which can realize nonlinear mapping. The radial basis function is better than the linear kernel function in nonlinear data processing. The shape of the function is a bellshaped curve, and only one parameter controls the width of the curve. Compared with the polynomial kernel function and the sigmoid kernel function, the radial basis function has fewer parameters, which greatly reduces the model’s complexity. The radial basis function is the most widely used kernel function. It has good performance in dealing with large and small samples. This paper uses the radial basis kernel function, whose expression is , where denotes the width parameter of the kernel function. Different affects the performance of the relevance vector machine, so this parameter must be optimized, and the optimization algorithm based on metaheuristic can significantly improve the optimization efficiency.
3.2. Parameter Optimization Based on HHO
HHO is a novel metaheuristic swarm intelligence algorithm. Compared with the genetic algorithm and the particle swarm optimization algorithm, it has better exploration capabilities and convergence speed [26]. When the hawks catch their prey, they adopt different strategies in different stages. The HHO algorithm can find the optimal value by modeling hawks’ behavior in each stage of the predation process.
The HHO algorithm includes the exploration phase and the exploitation phase. In the process of predation, with the increase in the moving distance of prey, the physical strength of prey will gradually decrease, which is called escape energy. The algorithm can transfer from the exploration phase to the exploitation phase through the value range of prey’s escape energy. The escape energy is modeled as [26] where denotes the escaping energy of the prey, denotes the maximum number of iterations, denotes the current number of iterations, and denotes the initial energy. randomly changes inside the interval at each iteration.
When , enter the exploration phase. The hawks inhabit different positions randomly and use two strategies to detect prey. An equal probability is used to distinguish the two strategies. They perch based on the positions of other hawks and the prey when and perch on random tall trees when . The exploration phase can be expressed as where denotes the position vector of the nextgeneration hawks, denotes the position vector of the currentgeneration hawks, denotes the position of the prey, ,, , , and are random values between and are updated every generation, and denote the boundary of the variable, denotes a random position selected from the currentgeneration hawks, and denotes the average position of the hawks.
When , enter the exploitation phase. At this time, use four siege strategies. is a random number between and indicates the chance for the prey to escape. indicates that the prey can escape, and indicates that the prey cannot escape successfully [26]. The use of soft and hard besiege strategies depends on the prey’s energy . When , it means that the prey still has high energy. At this time, the Hawks adopt the soft besiege strategy, constantly approaching and chasing the prey to consume the energy of the prey. When , it indicates that the prey’s energy is weak. At this time, the Hawks can easily catch the prey by using the hard besiege strategy.
(1) Soft Besiege. When and , it means that the prey has enough escape energy and jumps randomly. At this time, the hawks encircle it softly to exhaust the prey, expressed as where is a random number between , denoting the random jumping energy of prey, which is used to simulate the prey’s activity.
(2) Hard Besiege. When and , it means that the prey is weak and has low escape energy. At this time, the following equation is used to update the position.
(3) Soft Besiege with Progressive Rapid Dives. When and , the prey has enough energy to escape. At this time, equations (11) and (12) are used to update the position, and the Lévy flight is used to simulate the sudden movement of the hawks and the prey. The Lévy flight is expressed as equation (13), which can further enhance the ability of the algorithm to jump out of the local optimum. where and are random values inside , is a default constant set to 1.5, and is the Gramma function: .
Compare the two update methods, and select the most suitable position update method according to the fitness function value.
(4) Hard Besiege with Progressive Rapid Dives. When and , it means that the prey does not have enough energy. At this time, equations (15) and (16) are used to update the position and equation (14) is used to select the appropriate update method.
The relevance vector machine is used for data fitting and regression. HHO is used to optimize the kernel width parameter of the relevance vector machine.
The flow of each HHORVR is as follows:
Step 1. Initialize the number of HHO populations and the number of iterations. The hawk’s position in the HHO algorithm represents the kernel width parameter of the relevance vector machine. The position of each hawk in the firstgeneration group is randomly distributed.
Step 2. Determine the fitness function. The sample data is randomly divided into five parts, 4 of which are used for the relevance vector machine training, and the remaining 1 part is used for verification, and the test accuracy of the relevance vector machine is used as the fitness function of HHO.
Step 3. Calculate the fitness value of each hawk, and use the hawk corresponding to the optimal fitness as the prey.
Step 4. Update the positions of the hawks, according to equations (7)–(16).
Step 5. Perform iterative operations. Repeat steps 3 and 4 until the number of iterations is satisfied and the calculation is stopped. At this time, the prey position is the optimal kernel width parameter.
3.3. RLGIMU Temperature Modeling
The physical input of the RLG is the angular velocity, and the measured output is the number of pulses. The computer counts the number of pulses in the sampling time. According to the measurement equation of the system, the angular velocity of the RLG is obtained. The pulseangular velocity measurement model of the RLGs can be expressed as
is the vector of the number of pulses output by the RLGs on three orthogonal axes within the sampling time . is the scale factor matrix. is the bias matrix. is the angular velocity component vector of the three axes in the IMU measurement coordinate system. is the installation error matrix of the RLG relative inertial coordinate system for each axis.
In equation (17), the output vector of the RLGs changes with temperature and is very little affected by temperature. Because the influence of temperature on the RLG is multifaceted, the change of the sensor’s output to the measured angular velocity caused by the temperature change is directly reflected in the scale factor matrix and bias matrix , and the structural changes will also be reflected on and , so equation (17) can be expressed as
The physical input of the quartz flexible accelerometer is the apparent acceleration, and the measured output is the number of pulses. The computer counts the number of pulses in the sampling time. According to the measurement equation of the system, the apparent acceleration of the quartz flexible accelerometer is obtained. The pulseapparent acceleration measurement model of the quartz flexible accelerometers is expressed as
is the vector of the number of pulses output by the quartz flexible accelerometers on three orthogonal axes within the sampling time . is the scale factor matrix of the quartz flexible accelerometers. Since the quartz flexible accelerometer is very sensitive to positive and negative scale factors, it is divided into positive and negative: the positive scale factor matrix is expressed as and the negative scale factor matrix is expressed as . is the bias matrix of the quartz flexible accelerometers. is the apparent acceleration component vector of the three axes in the IMU measurement coordinate system. is the installation error matrix of the quartz flexible accelerometer relative inertial coordinate system for each axis. Similarly, equation (20) is expressed as
In equations (19) and (22), the input parameter affected by temperature is called the temperaturerelated parameter, and there are 15 temperaturerelated parameters in total.
Define a set to denote 15 temperaturerelated parameters. Make denote each element in and denote each element in .
Since the temperature measurement position of each RLG and quartz flexible accelerometer is different, the temperature measurement curve of each temperaturerelated parameter will be different. Define a temperature set . denotes the temperature change vector of each temperaturerelated parameter and corresponds onetoone with ; then, the corresponding relationship between set and temperature can be expressed as . From this, 15 models can be built as follows: where is the model corresponding to each temperaturerelated parameter of the IMU and is the optimal width parameter of the kernel function of each model.
3.4. System Compensation Model
Ideally, in equations (19) and (22), angular velocity and apparent acceleration are independent of temperature. First, the output models of angular velocity and apparent acceleration can be obtained according to IMU measurement equations (19) and (22) and the models of temperaturerelated parameters. The values of the temperaturerelated parameters at the desired temperature are then selected as constants in the angular/apparent acceleration measurement equation to obtain the compensation pulse. The desired temperature is a hypothetical stable working state. The purpose is that when the system works at any temperature, the measured value of physical quantity can be converted to the value at the desired temperature so that the system will not be affected by the temperature change. The desired temperature is not unique. It can be set according to the actual needs and requirements. No matter how much the desired temperature is set, the measured value of different temperatures can be converted to the measured value at the desired temperature by compensating for the measurement equation’s parameters to realize the temperature compensation.
For the RLGs, according to equation (19), the output model of is
corresponding to each RLG is , , and .
The compensation equation of RLG output pulse at the desired temperature is
is the output pulse of the three RLGs after compensation. and are, respectively, the scale factor matrix and bias matrix of the three RLGs at the desired temperature .
By bringing equation (24) into (25), there is
Corresponding to different RLGs, includes , , and . The data set can be obtained by temperature experiment. The data set was used to train the relevance vector machine. At the same time, the optimization process introduced in Section 3.2 is used to optimize the kernel width parameter of the relevance vector machine, so as to obtain the optimized model. Similarly, the model of can be obtained. Since the temperature models of and have been obtained, it is possible to compensate for the output at all temperatures, thereby compensating for to at different temperatures.
For the quartz flexible accelerometers, according to equation (22), the output model of is
corresponding to each quartz flexible accelerometer is , , and .
The compensation equation of quartz flexible accelerometer output pulse at the desired temperature is
is the output pulse of the three quartz flexible accelerometers after compensation. and are, respectively, the scale factor matrix and bias matrix of the three quartz flexible accelerometers at the desired temperature .
By bringing equation (27) into (28), there is
Corresponding to different quartz flexible accelerometers, includes , , and . Similarly, the models of and can be obtained. Since the temperature models of and have been obtained, it is possible to compensate for the output at all temperatures, thereby compensating for to at different temperatures.
Equations (26) and (29) are the compensation model of the IMU.
4. Experiment and Analysis
4.1. Data Acquisition and Modeling
The IMU used in the experiment uses three RLGs and three quartz flexible accelerometers as measurement sensors. The scale factor of the selected RLG on the IMU system is approximately 2.14 pulses/arcsec. The scale factor of the selected accelerometer on the IMU system is approximately . The RLGIMU is fixed on a twoaxis rotating platform in the temperature chamber. Set the temperature chamber to 5°C, and keep it for 4 hours to ensure the temperature inside the IMU and the chamber temperature are consistent. The temperature chamber is set to 7000 minutes from 5°C to 65°C. In the process of temperature rise in the temperature chamber, the IMU continuously carries out rotation calibration work in the temperature chamber. The IMU uses the systematic calibration method, and the rotation speed is 20 deg/s. Various errors will be excited by the continuous rotation of multiple positions. Using the Kalman filter, various error parameters of the RLGs and accelerometers can be calculated, that is, , , , , , and in Equations (17) and (20). Each calibration takes approximately 40 minutes. A total of 175 calibrations have been completed within the set temperature range and time. The IMU temperature change is caused not only by the temperature change in the temperature chamber but also by the heating part inside the IMU. The temperature sensor measurement shows that the temperature change between every two calibrations is about 0.4°C. The temperature change is tiny. Therefore, the temperature is considered to be stable in one calibration process. Through the experimental process described above, 175 sets of calibration data of the IMU can be obtained over the entire experimental temperature range. After offline calibration data processing, 175 sets of temperaturerelated parameters at 175 temperature points can be obtained.
Through experiments in the temperature chamber, the data set of each temperaturerelated parameter and the corresponding temperature set are obtained. There are 15 such sets in total. It can be expressed as . Each is used as the relevance vector machine training data. Use the HHORVR process in Section 3.2 to train each relevance vector machine. The kernel width parameter of the relevance vector machine is used as the input variable of the HHO. The total population size is set to 30, the max number of iteration is set to 500, and the value range is set to . The kernel width parameter with the highest prediction accuracy in the range of values is found using the 5fold crossvalidation method. Each temperaturerelated parameter is modeled according to this process. The set of models in equation (23) can be obtained.
For comparison, the least square method is used to model the 15 temperaturerelated parameters. The model equation is expressed as where is the output vector of the temperaturerelated parameters, is the coefficient matrix of the equation, and is the temperature input vector.
Fifteen temperaturerelated parameters were modeled by the RVR method and the least square method, respectively. In the RVR model, because the radial basis function is used as the kernel function, there is only one kernel width parameter. According to equation (30), the least square method of a cubic fitting polynomial is used. Each fitting polynomial needs to determine four coefficients: the row vector of matrix in equation (30). The kernel width parameters, the number of relevance vectors, the proportion of relevance vectors of the 15 RVR models, and the coefficients of 15 least square fitting polynomials are listed in Table 1.

4.2. Analysis of the Regression Performance of the Models
To compare the fitting and regression performance of HHORVR and LSM, HHORVR and LSM were used to model the temperature value of 15 temperaturerelated parameters, with a total of 30 models. Three indicators are used to measure the fitting and regression performance of the two modeling methods:
(1) Root Mean Square Error [36]. It is used to measure the error and dispersion between the regression value and the real value. The smaller the value, the smaller the estimation error and the smaller the error dispersion. The equation is as follows:
(2) Mean Absolute Error [37]. It is used to measure the average absolute errors. Reflect the actual situation of the predicted value error. The smaller the value, the smaller the estimated error. The equation is as follows:
(3) RSquared. It is used to measure the fitting accuracy of the model; the closer to 1, the higher the fitting accuracy. The equation is as follows:
In equations (32), (33), and (34), is the sample input, is the regression model, is the regression output, is the sample output value, is the mean of the sample output, and is the number of samples.
The results of the three indicators are shown in Table 2. Compared with the least square method, all the regression test performances of the 15 models were improved by using the HHORVR method. The root mean square error value increased by a maximum of 80.12% and an average of 35.27%. The mean absolute error value increased by a maximum of 83.51% and an average of 39.29%. The squared value increased by a maximum of 9.15% and an average of 2.29%. The results show that the RVR has better performance than the least square method in terms of error, error dispersion, and fitting accuracy.
