Research on the Shearer Positioning Method Based on the MEMS Inertial Sensors/Odometer Integrated Navigation System and RTS Smoother

The shearer positioning method with an inertial measurement unit and the odometer is feasible in the longwall coal-mining process. However, the positioning accuracy will continue to decrease, especially for the micro-electromechanical inertial measurement unit (MIMU). In order to further improve the positioning accuracy of the shearer without adding other external sensors, the positioning method of the Rauch-Tung-Striebel (RTS) smoother-aided MIMU and odometer is proposed. A Kalman filter (KF) with the velocity and position measurements, which are provided by the odometer and closing path optimal estimation model (CPOEM), respectively, is established. The observability analysis is discussed to study the possible conditions under which the error states of KF can be estimated. A RTS smoother with the above-mentioned KF as the forward filter is built. Finally, the experiments of simulating the movement of the shearer through a mobile carrier were carried out, with a longitudinal movement distance of 44.6 m and a lateral advance distance of 1.2 m. The results show that the proposed method can effectively improve the positioning accuracy. In addition, the odometer scale factor and mounting angles can be estimated in real time.


Introduction
Automated mining based on a longwall face has shown significant potential to improve mining productivity, increase personnel safety, and secure environmental sustainability [1]. As shown in Figure 1, the longwall face equipment include a shearer, some hydraulic supports, and an armored face conveyor (AFC). The shearer rides on the AFC to cut coal back and forth. The AFC provides the running track for the shearer while transporting the coal. The hydraulic supports not only support the roof but also push the AFC towards the coal seam. The position of the shearer is directly related to the control of the AFC and hydraulic supports [2]. Hence, the positioning of the shearer is the key technology to realize automated mining. The inertial measurement unit, which contains a 3D inertial sensor, is widely used to estimate the position of the shearer due to its high reliability and autonomy [3]. The micro-electromechanical inertial measurement unit (MIMU) is especially favored by the mine engineers with the advantages of low cost and small size. However, the free inertial position error can grow quickly over time due to the drifting of the inertial devices [1], which includes repeatability biases, the slow-varying drifts, and the fast-varying drifts [4,5], and thus the integrated navigation mode with MIMU as the core component becomes a better choice. The MIMU/Global Positioning System (GPS) integrated navigation system is widely used as a conventional and low-cost positioning method. Unfortunately, GPS cannot be used in underground environments. The positioning system composed by ultra-wideband range measurements can be used in GPS-denied Zero-velocity update (ZUPT)-aided MIMU is a simple and robust integrated st [9,10]. However, it requires short and frequent stops. The motion constraint-aided ZUPT method can reduce the number of stops for simple ZUPT correction [11,12], accuracy cannot meet the actual demand of shearer positioning. The odometer is reg as one of the most potentially useful autonomous speed sensors for land vehicles, a system with the MIMU and odometer has been proven to be autonomous, robu feasible for shearer positioning [13][14][15]. The dead reckoning (DR) with Euler angl vided by MIMU and velocity provided by the odometer is one of the methods to es the position of the shearer. The inertial navigation attitude error, especially the h angle error, is the main error source precluding further improvement of DR [1]. T tion constraint [16] and closing path optimal estimation model (CPOEM) [1] are u improve the accuracy of DR. However, these two methods can only slow down the gence of the position error. The reason is that DR, as an open-loop structure, cann vent the divergence of inertial navigation attitude errors. Establishing a Kalman filt based on MIMU and the odometer is another information fusion method, which is used in the field of conventional land navigation [17][18][19][20]. The advantages of this m lie in the closed-loop correction of the MIMU attitude and the real-time estimation inertial device parameters. The effect of correction and estimation has a lot to do w steering maneuver of the vehicle. In the case of only MIMU and the odometer, the and slow steering maneuver of the longwall shearer is difficult to achieve the idea tioning effect [20,21]. Therefore, we plan to build a KF based on MIMU, odometer, m constraint, and CPOEM to achieve better positioning.
The position error of the MIMU/odometer integrated system repeats "increa crease" changes along with the reciprocating operation of the shearer [22], which typical case of the Rauch-Tung-Striebel (RTS) smoother. RTS smoothing is a tech Zero-velocity update (ZUPT)-aided MIMU is a simple and robust integrated strategy [9,10]. However, it requires short and frequent stops. The motion constraint-aided MIMU ZUPT method can reduce the number of stops for simple ZUPT correction [11,12], but its accuracy cannot meet the actual demand of shearer positioning. The odometer is regarded as one of the most potentially useful autonomous speed sensors for land vehicles, and the system with the MIMU and odometer has been proven to be autonomous, robust, and feasible for shearer positioning [13][14][15]. The dead reckoning (DR) with Euler angles provided by MIMU and velocity provided by the odometer is one of the methods to estimate the position of the shearer. The inertial navigation attitude error, especially the heading angle error, is the main error source precluding further improvement of DR [1]. The motion constraint [16] and closing path optimal estimation model (CPOEM) [1] are used to improve the accuracy of DR. However, these two methods can only slow down the divergence of the position error. The reason is that DR, as an open-loop structure, cannot prevent the divergence of inertial navigation attitude errors. Establishing a Kalman filter (KF) based on MIMU and the odometer is another information fusion method, which is widely used in the field of conventional land navigation [17][18][19][20]. The advantages of this method lie in the closed-loop correction of the MIMU attitude and the real-time estimation of the inertial device parameters. The effect of correction and estimation has a lot to do with the steering maneuver of the vehicle. In the case of only MIMU and the odometer, the small and slow steering maneuver of the longwall shearer is difficult to achieve the ideal positioning effect [20,21]. Therefore, we plan to build a KF based on MIMU, odometer, motion constraint, and CPOEM to achieve better positioning.
The position error of the MIMU/odometer integrated system repeats "increasedecrease" changes along with the reciprocating operation of the shearer [22], which fits the typical case of the Rauch-Tung-Striebel (RTS) smoother. RTS smoothing is a technology that uses all observation information in a certain time interval to re-estimate the state based on the KF algorithm [23]. The RTS smoothing algorithm is one of the core technologies of the position and orientation system, which can enhance the ability of the position and orientation system to be free from external disturbance [24][25][26][27]. Since the estimation accuracy of the RTS smoothing algorithm is superior to that of filtering, it is often used as a reference for post-analysis of the integrated navigation system [28,29]. In the above applications, RTS smoothing works in off-line mode, which not only requires a large storage space, but also limits its use in scenarios with real-time requirements. An on-line smoothing method can overcome the above limitations and has been successfully used in the MIMU/GPS integrated navigation system [30] and pedestrian navigation system [31]. The core idea of on-line RTS smoothing can be summarized as smoothing is executed immediately after a certain time window to achieve a near real-time application effect. Obviously, on-line RTS smoothing technology is more suitable for the needs of longwall mining. Therefore, we propose a positioning method based on MIMU, odometer, motion constraint, CPOEM, and on-line RTS smoother, whose contributions and benefits can be summarized as follows: (1) The RTS smoother is introduced into the MIMU and odometer integrated system, which can improve the shearer positioning accuracy without additional external sensors.
(2) Observability analysis of position measurement is added to theoretically provide the estimation conditions of the main error states on the basis of the previous work on velocity measurement.
(3) The mounting angles between the MIMU frame and the odometer frame can be estimated in real time, avoiding the tedious pre-calibration process.

Mathematical Models of Velocity and Position
The coordinate systems involved in this paper are shown in Figure 2, which are defined as follows: b-frame, the MIMU frame, which originates at the sensitive center of the MIMU, with the axes pointing to the AFC advance direction (right), shearer moving direction (forward), and upward; m-frame, the odometer frame, whose axes point right, forward, and upward; n-frame, the local-level east-north-up coordinate.
based on the KF algorithm [23]. The RTS smoothing algorithm is one of the core technologies of the position and orientation system, which can enhance the ability of the position and orientation system to be free from external disturbance [24][25][26][27]. Since the estimation accuracy of the RTS smoothing algorithm is superior to that of filtering, it is often used as a reference for post-analysis of the integrated navigation system [28,29]. In the above applications, RTS smoothing works in off-line mode, which not only requires a large storage space, but also limits its use in scenarios with real-time requirements. An on-line smoothing method can overcome the above limitations and has been successfully used in the MIMU/GPS integrated navigation system [30] and pedestrian navigation system [31]. The core idea of on-line RTS smoothing can be summarized as smoothing is executed immediately after a certain time window to achieve a near real-time application effect. Obviously, on-line RTS smoothing technology is more suitable for the needs of longwall mining. Therefore, we propose a positioning method based on MIMU, odometer, motion constraint, CPOEM, and on-line RTS smoother, whose contributions and benefits can be summarized as follows: (1) The RTS smoother is introduced into the MIMU and odometer integrated system, which can improve the shearer positioning accuracy without additional external sensors.
(2) Observability analysis of position measurement is added to theoretically provide the estimation conditions of the main error states on the basis of the previous work on velocity measurement.
(3) The mounting angles between the MIMU frame and the odometer frame can be estimated in real time, avoiding the tedious pre-calibration process.

Mathematical Models of Velocity and Position
The coordinate systems involved in this paper are shown in Figure 2, which are defined as follows: b-frame, the MIMU frame, which originates at the sensitive center of the MIMU, with the axes pointing to the AFC advance direction (right), shearer moving direction (forward), and upward; m-frame, the odometer frame, whose axes point right, forward, and upward; n-frame, the local-level east-north-up coordinate. The aim of the strapdown inertial navigation system (SINS) alignment process is to determine the transition matrix from b-frame to n-frame, denoted by C n b . One of the purposes of MIMU and odometer joint calibration is to determine the mounting angles between m-frame and b-frame, which can be expressed by a vector, α . The direction cosine matrix from b-frame to m-frame is denoted by C m b . The aim of the strapdown inertial navigation system (SINS) alignment process is to determine the transition matrix from b-frame to n-frame, denoted by C n b . One of the purposes of MIMU and odometer joint calibration is to determine the mounting angles between m-frame and b-frame, which can be expressed by a vector, α. The direction cosine matrix from b-frame to m-frame is denoted by C m b .

Measured Velocity Model
The moving speed of the shearer is denoted by v m y , which can be measured by the odometer. Taking into account the scale factor error, δk D , and measurement noise, w OD , of the odometer, the actual output, v m y , of the odometer can be expressed as [20,21,32]: According to the motion constraint, there is no sideslip along the AFC advance direction and no motion normal to the AFC under ideal conditions, so the velocities along the x m axis and z m axis are regarded as zero [11]. The ideal velocity of the shearer in m-frame can be expressed as: Integrating (1) and (2), the measured velocity model in m-frame is: where w OD/MC is the noise vector, defined as: where w MC,x and w MC,z are the motion constraint noise along the x m axis and the z m axis, respectively.

Measured Position Model
A typical shearer operation process is to repeat "straight cutting-oblique cuttingreverse straight cutting" to form a closed path, which is shown in Figure 3 [1]. The shearer runs in the order of A-B-C-D-E-F-G-H-I-J. Assuming that the shearer moves from the right to the left of Figure 3 during the first cutting cycle, we define the end corresponding to A, F, and I as the near end of the longwall face, and the corresponding end of the B, E, and J as the far end of the longwall face. The lengths of j − 1, j, and j + 1 cutting cycles, corresponding to the lengths of AB, EF, and IJ, respectively, are the same. The advance distance between two adjacent cutting cycles is also the same, denoted by d, which can be measured by the displacement sensor fixed in the push arm of the hydraulic support. Some points with ∆ interval in each cutting cycle are selected as optimal points for information fusion of KF. The ideal advance displacement of an optimal point can be expressed as: where D m i,j−1|j represents the advance displacement of point i from the j − 1 cutting cycle to the j cutting cycle.
According to the CPOEM principle and the longwall mining process [1], the positions of the optimal points during the next cutting cycle can be predicted through the positions of the current cutting cycle and the advance displacement. The position of point i in the j + 1 cutting cycle is expressed as: where p i,j and p i,j+1 denote the positions of point i in the j and j + 1 cutting cycles, respectively, which are expressed in the form of longitude, λ, latitude, L, and height, h, and matrix C, whose function is to convert the position increment in n-frame into the form of longitude, latitude, and height, is: where R N and R M denote the transverse and meridian radius of curvature, respectively, and they are the parameters used to describe the earth ellipsoid model, which are regarded as constant values in this paper.
where N R and M R denote the transverse and meridian radius of curvature, resp and they are the parameters used to describe the earth ellipsoid model, which garded as constant values in this paper. The initial positions of the optimal points are provided in the first cutting cy specific process is as follows:


The two optimal points corresponding to both ends of the longwall face can rately measured in advance, as mentioned in [16,33].


To avoid additional surveying and mapping, the initial positions of the re optimal points can be obtained using the position estimates of the integrated in the first cutting cycle. This initial value assignment method is mentioned Taking into account the existence of measurement noise, the measured p  p / CPO EM i , of point i predicted by CPOEM can be approximated as:

Error State Equation of Integrated Navigation System
The mounting angles between m-frame and b-frame are inevitable, even if calibration process is carried out. The residual installation errors after calibratio regarded as random constants, denoted by δ δα δα δα The scale factor the odometer can also be regarded as a random constant. Thus, the following e can be obtained: The direction cosine matrix,  C m b , with the residual error δα satisfies: where I is a third-order unit matrix and δ × ( ) α denotes the skew symmetric m δα with 3 rows and 3 columns.
When the mounting angles, α , are controlled within a small range by a pre chanical installation, the equation can be obtained directly without the perform the pre-calibration process. Some related error state vectors of the SINS satisfy the following equation: The initial positions of the optimal points are provided in the first cutting cycle. The specific process is as follows:

•
The two optimal points corresponding to both ends of the longwall face can be accurately measured in advance, as mentioned in [16,33].

•
To avoid additional surveying and mapping, the initial positions of the remaining optimal points can be obtained using the position estimates of the integrated system in the first cutting cycle. This initial value assignment method is mentioned in [1].
Taking into account the existence of measurement noise, the measured position, p i/CPOEM , of point i predicted by CPOEM can be approximated as: where w i/CPOEM is the CPOEM noise vector of point i.

Error State Equation of Integrated Navigation System
The mounting angles between m-frame and b-frame are inevitable, even if the precalibration process is carried out. The residual installation errors after calibration can be regarded as random constants, denoted by δα = [δα x δα y δα z ] T . The scale factor error of the odometer can also be regarded as a random constant. Thus, the following equation can be obtained: δ where 0 a×b is a a × b zero matrix. The direction cosine matrix, C m b , with the residual error δα satisfies: where I is a third-order unit matrix and (δα×) denotes the skew symmetric matrix of δα with 3 rows and 3 columns. When the mounting angles, α, are controlled within a small range by a precise mechanical installation, the equation C m b = I can be obtained directly without the need to perform the pre-calibration process.
Some related error state vectors of the SINS satisfy the following equation: where v n , C b n , and p are the true velocity, attitude matrix, and position, respectively, and v n SINS , C b n , and p SINS are the error-contaminated velocity, attitude matrix, and position, respectively, calculated by the SINS. δv n = [δv E δv N δv U ] T and δp are the velocity errors and position errors of the SINS, respectively, ϕ n = [ϕ E ϕ N ϕ U ] T is the misalignment angles of C b n in n-frame, and the subscripts E, N, and U are the east, north, and up directions in n-frame.
A 19-dimensional error state vector is defined as: where ε b and ∇ b are the gyro and accelerometer biases, respectively. Taking into account the low-speed motion characteristics of the shearer, the velocityrelated terms in the SINS error equation can be ignored. Then, the SINS error equation can be simplified to: where ω n ie denotes the rotation rate vector of the earth and f n = C n b f b , in which f b is the specific force measured by the accelerometers.
The error state equation of the integrated navigation system can be expressed as: .
where w(t) is the noise vector of the integrated system and the 15 × 15 transition matrix F SINS is denoted by: The discretized error state equation corresponding to (14) is: where X k and X k−1 are the discretized error state vectors at t k and t k−1 , respectively, W k−1 denotes the system noise matrix at t k−1 , and Φ k,k−1 is the discretized state transition matrix that satisfies the following equation: where T denotes the filtering period.

Measurement Equations of Velocity and Position
The velocity, v m SINS , in m-frame calculated by the SINS can be expressed as: where We find the difference between (3) and (18) to obtain the velocity measurement equation as: (19) where the measurement matrix, H v , of velocity is expressed as: Similarly, the position measurement equation can be expressed as: where the measurement matrix, H p , of the position is given by: The update of the position measurement depends on the mileage calculated by the odometer, which means that the update period of the position measurement is an integer multiple of that of the velocity measurement. Therefore, the overall measurement model includes two forms: simultaneous velocity and position measurement update and separate velocity measurement update. The specific update conditions and model equations are as follows.
If the shearer is not in the first cutting cycle and the position measurement is judged to be valid by the mileage calculated by the odometer, then: where z and H, respectively, represent the overall measurement vector and measurement matrix. Otherwise: The discretized measurement equation is: where Z k , H k , and V k are the discretized measurement vector, the discretized measurement matrix, and the measurement noise sequence at t k , respectively.

Observability Analysis of Integrated System
The main purpose of observability analysis is to study the observability and estimatable conditions of the error state. An observability analysis method directly relies on the state and measurement equations to investigate the observability.

Observability Analysis Based on the Velocity Measurement
The observability analysis based on the velocity measurement has been discussed in detail in our previous work [34], so here, we directly summarize the observability conclusions about the velocity measurement with the motion characteristics of the longwall shearer.

•
The velocity errors, δv n , are observable, and the estimation accuracy is related to the estimation degree of other error states.

•
The position errors, δp, are unobservable, but the estimation accuracy will still be improved with the effective estimation of the velocity errors.

•
The acceleration and deceleration process of the shearer is the premise of exciting the error states δα x , δα z , and δk D , which contribute to the positioning errors. The error δα y is unobservable.

•
The error ∇ b z is observable. The separation of ∇ b x and ∇ b y , and the distinction of ε b x and ε b y , improving the estimation accuracy of ϕ n , depend on the turning motion of the shearer. The azimuth error, ϕ U , is directly related to the lateral positioning error, thus, restricting the estimation accuracy of error δα z .
It can be concluded that frequent turning of the shearer is necessary to improve the estimation accuracy of the error states. However, limited by working conditions, the longwall shearer has very few steering maneuvers. Therefore, it is difficult to achieve high shearer positioning accuracy using only the integrated system of MIMU and the odometer.

Observability Analysis Based on the Position Measurement
In (21), the measurement values are constructed by the position errors. Therefore, the position errors, δp, are observable. The estimation accuracy of δp depends on the position accuracy of the optimal points predicted by CPOEM. The vectors that make up the measurement values have the same initial position errors, so the initial position errors cannot be estimated.
Taking the time derivative of (21) obtains: Equation (26) indicates that the velocity errors, δv n , are observable, and the estimation accuracy is determined by the position errors.
The terms related to v n in the mathematical models can be ignored due to the low-speed characteristics of the shearer. Therefore, the time derivative of (26) can be simplified as: and ω ie is the rotation rate of the earth. The acceleration and deceleration process of the shearer makes at least one of f E and f N non-zero, which means that the coefficient of ϕ U will not be 0 in (27). That is to say, the acceleration and deceleration process can improve the observability of the azimuth error, ϕ U .
Let us assume a high-precision position predicted by CPOEM is obtained. Then, the errors δp and δv n will be accurately estimated. The estimation accuracy of the azimuth error, ϕ U , depends on the separating degree from the errors ϕ E , ϕ N , ∇ E , and ∇ N . Although the state errors δk D , δα x , and δα z are not directly observable in (21), (26), and (27), their accuracy will be improved with the estimation of error states δp and ϕ n . Similarly, the estimation accuracy of ε b will also be improved with the estimation of other state errors.

RTS Smoothing
RTS smoothing uses all the measurement values obtained in a time interval to estimate the error states at every epoch in this interval. A typical RTS smoothing process is divided into two steps: forward filtering and backward smoothing.
The forward filtering is obtained through standard KF, whose basic equation for discretization can be expressed as: where X k,k−1 and P k,k−1 are the one-step predicted states and covariance at t k , calculated from the information at t k−1 , P k−1 and P k denote the state estimate covariance at t k−1 and t k , respectively, and Q k−1 is the variance matrix of the system noise sequence W k−1 at t k−1 . K k represents the filter gain matrix at t k , and R k is the variance matrix of the measurement noise sequence V k at t k . The forward filtering in this paper refers to the KF based on the velocity and position measurements mentioned above. In the forward filtering process, it is required to save X k,k−1 , X k , Φ k,k−1 , P k,k−1 , and P k at every epoch in the time interval. After completing the forward filtering in the time interval, the backward smoothing is performed. The procedure of the backward smoothing is broken down into the following steps:

Backward Smoothing Initialization
Define the time interval as [t j ,t j+N ], then the initialization equation is given as: where X j+N and P j+N are the error state vector and the mean square error matrix of the forward filtering process at t j+N , and X rts,j+N and P rts,j+N , which are the initial values required for backward smoothing, denote the error state vector and the mean square error matrix of the backward smoothing process at t j+N .

Backward Smoothing Update
The recursive equations of the backward smoothing update are given by: where A k is the smoother gain, determined as: Sorting out the process of KF and the RTS smoother, the flowchart is shown in Figure 4. If mod(S, ∆) = = 0, where S denotes the mileage of the shearer calculated by the odometer, the optimal point from the CPOEM data is valid and the position measurement update of KF is performed. Backward smoothing is only executed when the shearer is at both ends of the longwall face and not in the first cutting cycle. The positions of the optimal points of the next cycle will be updated immediately after the end of backward smoothing.
Sorting out the process of KF and the RTS smoother, the flowchart is shown in Figure  4. If mod(S, Δ) = = 0, where S denotes the mileage of the shearer calculated by the odometer, the optimal point from the CPOEM data is valid and the position measurement update of KF is performed. Backward smoothing is only executed when the shearer is at both ends of the longwall face and not in the first cutting cycle. The positions of the optimal points of the next cycle will be updated immediately after the end of backward smoothing.

Experiments
To evaluate the performance of the proposed positioning method, experiments were carried out, as shown in Figure 5. A mobile carrier equipped with the MIMU (Xsens MTi-G-700) and the odometer simulated the movement of the shearer. The MIMU was installed on the mobile carrier through an adapter plate and the odometer was connected to its wheel. A GPS receiver with an antenna was also installed on the mobile carrier. The GPS receiver can output centimeter-level positioning results through the network differential technology. The positioning result of the GPS receiver only provides an evaluation basis for tests. The specifications of the MIMU and the initial errors are listed in Table 1, which are related to the initial parameter configuration of the Kalman filter. The specifications and the initial attitude error refers to the MTI user manual [35]. The initial position error refers to the network differential positioning accuracy [36]. The initial parameters of the filter are set as Appendix A.
The mobile carrier simulated four cutting processes of the shearer with a reciprocating travel distance of 44.6 m and an advance distance of 1.2 m, as shown in Figure 6. Figure  6 is drawn by the network differential results provided by the GPS receiver. The symbols "*" and "o" denote the start and end of the trajectory, respectively.

Experiments
To evaluate the performance of the proposed positioning method, experiments were carried out, as shown in Figure 5. A mobile carrier equipped with the MIMU (Xsens MTi-G-700) and the odometer simulated the movement of the shearer. The MIMU was installed on the mobile carrier through an adapter plate and the odometer was connected to its wheel. A GPS receiver with an antenna was also installed on the mobile carrier. The GPS receiver can output centimeter-level positioning results through the network differential technology. The positioning result of the GPS receiver only provides an evaluation basis for tests. The specifications of the MIMU and the initial errors are listed in Table 1, which are related to the initial parameter configuration of the Kalman filter. The specifications and the initial attitude error refers to the MTI user manual [35]. The initial position error refers to the network differential positioning accuracy [36]. The initial parameters of the filter are set as Appendix A.
The mobile carrier simulated four cutting processes of the shearer with a reciprocating travel distance of 44.6 m and an advance distance of 1.2 m, as shown in Figure 6. Figure 6 is drawn by the network differential results provided by the GPS receiver. The symbols "*" and "o" denote the start and end of the trajectory, respectively.
In order to study the influence of the selection of optimal points on the estimation of error states, we set the optimal points interval, ∆, as 2, 6, and 10 m, in turn. Figure 7 shows the positioning errors of the proposed integrated system without performing RTS smoothing. A1-A4 corresponded to the time periods of the first to fourth cutting cycles, respectively. The remaining time periods corresponded to the process of the shearer processing the end coal seam. It can be seen that during the first cutting cycle, the curves of different values of ∆ overlapped in the east, north, and height directions, respectively, since position measurement filtering was not performed. During the second to fourth cutting cycles, the curves of the east and north errors were smoother as ∆ decreased. Figure 8 shows the positioning errors of the proposed integrated system with performing RTS smoothing. It was straightforward to see that the positioning accuracy of all three axes with performing RTS smoothing was higher than that of not performing, regardless of the value of ∆. In addition, it can be seen that there were larger burrs in the east and north directions as ∆ increased. Therefore, we can choose a smaller value of ∆ to obtain smoother position estimations.   In order to study the influence of the selection of optimal points on the estimation of error states, we set the optimal points interval, Δ, as 2, 6, and 10 m, in turn. Figure 7 shows the positioning errors of the proposed integrated system without performing RTS smoothing. A1-A4 corresponded to the time periods of the first to fourth cutting cycles, respec-      positioning errors of the proposed integrated system with performing RTS smoothing was straightforward to see that the positioning accuracy of all three axes with performi RTS smoothing was higher than that of not performing, regardless of the value of Δ. addition, it can be seen that there were larger burrs in the east and north directions as increased. Therefore, we can choose a smaller value of Δ to obtain smoother position es mations.   Figures 9 and 10, respectively, show the estimations of the accelerometer and gy biases with and without performing RTS smoothing, while we set Δ as 2 m. It can be se that the estimations of MIMU biases were equivalent with and without performing R smoothing when the estimators tended to be stable. Although the true values of MIM biases cannot be accurately obtained, the results of the pure navigation before and af the bias compensation can reflect their estimation accuracy. We used the first 300 s of da to perform the pure navigation calculations before and after the bias compensation, a the results are shown in Figure 11. Compared with the results before compensation of t was straightforward to see that the positioning accuracy of all three axes with performi RTS smoothing was higher than that of not performing, regardless of the value of Δ. addition, it can be seen that there were larger burrs in the east and north directions as increased. Therefore, we can choose a smaller value of Δ to obtain smoother position e mations.   Figures 9 and 10, respectively, show the estimations of the accelerometer and gy biases with and without performing RTS smoothing, while we set Δ as 2 m. It can be se that the estimations of MIMU biases were equivalent with and without performing R smoothing when the estimators tended to be stable. Although the true values of MIM biases cannot be accurately obtained, the results of the pure navigation before and af the bias compensation can reflect their estimation accuracy. We used the first 300 s of d to perform the pure navigation calculations before and after the bias compensation, a the results are shown in Figure 11. Compared with the results before compensation of  Figures 9 and 10, respectively, show the estimations of the accelerometer and gyro biases with and without performing RTS smoothing, while we set ∆ as 2 m. It can be seen that the estimations of MIMU biases were equivalent with and without performing RTS smoothing when the estimators tended to be stable. Although the true values of MIMU biases cannot be accurately obtained, the results of the pure navigation before and after the bias compensation can reflect their estimation accuracy. We used the first 300 s of data to perform the pure navigation calculations before and after the bias compensation, and the results are shown in Figure 11. Compared with the results before compensation of the MIMU biases, the positioning accuracy was greatly improved after compensation. Therefore, the MIMU biases can be effectively estimated. The RTS smoothing technology only reduced the error fluctuations in the estimation process, and did not affect the final estimation accuracy.
MIMU biases, the positioning accuracy was greatly improved after compensation. Therefore, the MIMU biases can be effectively estimated. The RTS smoothing technology only reduced the error fluctuations in the estimation process, and did not affect the final estimation accuracy.   MIMU biases, the positioning accuracy was greatly improved after compensation. Therefore, the MIMU biases can be effectively estimated. The RTS smoothing technology only reduced the error fluctuations in the estimation process, and did not affect the final estimation accuracy.   The estimation results of the odometer scale factor and mounting angles are shown in Figure 12. We know that δα y is always unobservable from the observability analysis.
Therefore, its estimation curve is not drawn in this paper. The estimated scale factor and The estimation results of the odometer scale factor and mounting angles are shown in Figure 12. We know that δα y is always unobservable from the observability analysis. Therefore, its estimation curve is not drawn in this paper. The estimated scale factor and mounting angles, whose record order was δkα = [δk D δα X δα z ] T , with and without performing RTS smoothing, were δkα 1 = [0.017 0.333 • − 0.307 • ] T and δkα 2 = [0.029 0.352 • 0.093 • ] T , respectively. The positioning accuracy of the DR algorithm is restricted by δkα. In other words, the DR results can reflect the estimation effect of these three error states. The DR navigation calculations were performed after the compensation of δkα 1 and δkα 2 using the data that has been compensated for the MIMU biases, and the results are shown in Figure 13. The horizontal positioning accuracy after compensation of δkα 1 was better than that after compensation of δkα 2 , which means that the estimation accuracy of the error states δk D and δα z with performing RTS smoothing was improved. According to the theory of observability analysis, the estimation accuracy of error states (such as δk D and δα z ) can be improved by feedback when the position errors of the optimal points were reduced. Since the positioning accuracy with performing the RTS smoothing was better than that without performing (shown in Figures 7 and 8), the position errors of the corresponding optimal points were also smaller. Therefore, the results of Figure 13 are consistent with the observability analysis.  The shearer is a long and narrow machine with limited space for the installation of external sensors. The small size of MEMS inertial sensors makes them a goal pursued by researchers. However, small-sized and low-cost inertial sensors often have low measurement accuracy. Therefore, we studied the impact of MIMU, which has a lower accuracy  The shearer is a long and narrow machine with limited space for the installation of external sensors. The small size of MEMS inertial sensors makes them a goal pursued by researchers. However, small-sized and low-cost inertial sensors often have low measurement accuracy. Therefore, we studied the impact of MIMU, which has a lower accuracy The shearer is a long and narrow machine with limited space for the installation of external sensors. The small size of MEMS inertial sensors makes them a goal pursued by researchers. However, small-sized and low-cost inertial sensors often have low measurement accuracy. Therefore, we studied the impact of MIMU, which has a lower accuracy than MTi-G-700, on positioning accuracy of the mobile carrier. The idea of MIMU data generation is to superimpose errors on the original gyroscope and accelerometer data corresponding to Figure 6 to simulate lower-precision MEMS inertial sensor data. The error components of the gyroscope and accelerometer include repeatability biases, slow-varying drifts, and fast-varying drifts [4,5]. The repeatability biases can be regarded as a random constant. The slow-varying drift can be approximated as white noise due to the short correlation time of MIMU. The fast-varying drift is often abstracted as a white noise process. The white noise process is usually evaluated by random walk. In summary, the error model of the gyroscope and accelerometer can be expressed as: where ε b s and ∇ b s represent the total gyroscope and accelerometer errors, respectively, ε b and ∇ b denote the constant biases of the gyroscope and accelerometer, respectively, and ε b w and ∇ b w are the random walk noises of the gyroscope and accelerometer, respectively. A series of MIMU errors were designed, as listed in Table 2. The errors in Table 2 were added to the original data corresponding to Figure 6. The position errors before and after RTS smoothing are shown in Figures 14 and 15, respectively. Similar to the phenomenon in Figures 7 and 8, even if the accuracy of the inertial sensors is reduced, the positioning accuracy can still be improved after performing RTS smoothing. It can be seen from Figure 15 that the east errors of Par 3 and Par 2 are much larger than those of other parameters, and the east error of Par 3 is greater than that of Par 2. This is because the constant biases of the gyroscope are not fully estimated, and the residual constant component is larger as the constant biases increase. Correspondingly, the attitude and position errors are bound to be greater. In summary, it can be concluded that the proposed method can still improve the positioning accuracy even if the low-precision inertial sensors are used, and the positioning accuracy is closely related to the sensor parameters, especially the gyroscope.   In order to further verify the performance of the method proposed in this paper, it is compared with the traditional method. The method mentioned in [1] is state-of-the-art based on the SINS and the odometer. However, this method does not have the ability to autonomously estimate the biases of inertial sensors, mounting angles between MIMU and the odometer, and the scale factor error of the odometer. Considering that the above parameters have a great influence on the traditional method, it was compensated with the parameter values obtained by the proposed method. The positioning results of the traditional method before and after the estimated parameter compensation and the method proposed in this paper are shown in Figure 16. It can be seen from Figure 16 that after the parameters are compensated, the positioning errors of the traditional method are significantly reduced. It can also be seen that the accuracy of the method proposed in this paper is better than the traditional method. In order to visually describe the accuracy improvement range of the proposed method, the spherical error probable (SEP) is calculated. The SEP is a universal evaluation method of 3D positioning accuracy [37], which is listed in Table 3. A1-A4 correspond to the first to fourth cutting cycles, respectively. It can be seen from Table 3 that compared with the traditional method without compensation, the positioning accuracy of the proposed method increases by 71.43%, 83.92%, 94.30% and 92.64% in turn from A1 to A4. Since the position measurement is not performed in A1, and the optimal point position of A2 is related to that of A1, the SEP of the proposed method is In order to further verify the performance of the method proposed in this paper, it is compared with the traditional method. The method mentioned in [1] is state-ofthe-art based on the SINS and the odometer. However, this method does not have the ability to autonomously estimate the biases of inertial sensors, mounting angles between MIMU and the odometer, and the scale factor error of the odometer. Considering that the above parameters have a great influence on the traditional method, it was compensated with the parameter values obtained by the proposed method. The positioning results of the traditional method before and after the estimated parameter compensation and the method proposed in this paper are shown in Figure 16. It can be seen from Figure 16 that after the parameters are compensated, the positioning errors of the traditional method are significantly reduced. It can also be seen that the accuracy of the method proposed in this paper is better than the traditional method. In order to visually describe the accuracy improvement range of the proposed method, the spherical error probable (SEP) is calculated. The SEP is a universal evaluation method of 3D positioning accuracy [37], which is listed in Table 3. A1-A4 correspond to the first to fourth cutting cycles, respectively. It can be seen from Table 3 that compared with the traditional method without compensation, the positioning accuracy of the proposed method increases by 71.43%, 83.92%, 94.30% and 92.64% in turn from A1 to A4. Since the position measurement is not performed in A1, and the optimal point position of A2 is related to that of A1, the SEP of the proposed method is larger in A1 and A2. It is significantly reduced after A3. Compared with the traditional method with compensation, the positioning accuracy of the proposed method increases by 60.61% in A3 and is equivalent in A4. The above phenomenon not only proves the superiority of the proposed method in positioning accuracy, but also further shows that the method proposed in this paper can effectively estimate some parameters, such as biases of the MIMU, odometer scale factor, etc.

Conclusions
This paper proposed a positioning method of the shearer based on an integrated system and RTS smoothing technology. Performing RTS smoothing on the basis of the Kalman filter is a major feature of this paper. An experiment was carried out to verify the performance of the proposed positioning method. The experimental results showed that the positioning accuracy after performing RTS smoothing was significantly improved, which was closely related to the sensor parameters, and the estimatable ability of some error states was improved. In addition, a comparison with traditional methods was also carried out. The result shows that the positioning accuracy of the proposed method can be improved by at least 60.61%.

Conclusions
This paper proposed a positioning method of the shearer based on an integrated system and RTS smoothing technology. Performing RTS smoothing on the basis of the Kalman filter is a major feature of this paper. An experiment was carried out to verify the performance of the proposed positioning method. The experimental results showed that the positioning accuracy after performing RTS smoothing was significantly improved, which was closely related to the sensor parameters, and the estimatable ability of some error states was improved. In addition, a comparison with traditional methods was also carried out. The result shows that the positioning accuracy of the proposed method can be improved by at least 60.61%.
According to the existing theories and the experimental results in this paper, it can be seen that RTS smoothing technology has a significant improvement effect on the jump phenomenon of the error states. Therefore, not only the position measurement, but other excellent measurement information that can be captured may also cause the error state to jump. At this time, RTS smoothing can still play an important role, which can be further studied.

Conflicts of Interest:
The authors declare no conflict of interest.