Next Article in Journal
Innovative Load Forecasting Models and Intelligent Control Strategy for Enhancing Distributed Load Levelling Techniques in Resilient Smart Grids
Previous Article in Journal
Using Mixed Reality for Control and Monitoring of Robot Model Based on Robot Operating System 2
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dynamic Attitude Inertial Measurement Method for Typical Regions of Deck Deformation

1
School of Instrumentation Science and Engineering, Harbin Institute of Technology, Harbin 150001, China
2
Anshan Industrial Technology Research Institute of Harbin Institute of Technology, Anshan 114000, China
*
Author to whom correspondence should be addressed.
Electronics 2024, 13(17), 3555; https://doi.org/10.3390/electronics13173555
Submission received: 21 July 2024 / Revised: 30 August 2024 / Accepted: 5 September 2024 / Published: 6 September 2024

Abstract

:
Due to the deformation of ships, it becomes difficult to ensure the accuracy of attitude measurement in typical areas on the deck, which seriously impacts the safety and operational efficiency of shipborne equipment. To address this issue, this paper presents a parameter identification method for dynamic deformation models based on angle increment differences and introduces the related vector machine (RVM) algorithm for online estimation of dynamic deformation model parameters. In view of the truncation error and non-Gaussian noise of the model, this article proposes a dynamic attitude measurement method based on model predictive filtering (MPF), constructs a dynamic measurement model using Rodrigues parameters in an inertial frame, and designs a maximum correlation entropy (MCE) robust filter to achieve robust estimation of deck dynamic deformation. The performance of the method is verified through simulation analysis and shipborne experiments. The comparative results indicate that, compared with existing methods, the proposed improved deck dynamic attitude measurement algorithm based on model prediction (IDAM) can substantially enhance the accuracy of attitude measurement in the presence of deck dynamic deformations.

1. Introduction

To measure the attitude in the area equipped with shipborne equipment on the deck, a sub-inertial measurement unit (IMU) can be installed, and the attitude information can be derived from the inertial measurement information [1,2]. But the sub-optimal performance of the sub-IMU fails to meet the accuracy requirement for attitude, and it is necessary to use the ship’s main inertial navigation information for auxiliary measurement [3,4,5]. However, due to the fact that the ship is not an ideal rigid body, it may undergo deformation under external influences such as sunlight and wave impact [6,7]. The deformation of the ship manifests as static and dynamic types [8]. The static deformation of the ship can be estimated and regarded as the installation error angle between the main and sub-inertial navigation systems. The dynamic deformation of the ship will directly affect the gyroscope output of the inertial measurement unit and result in additional errors caused by changes of the dynamic lever arm, thereby affecting the attitude calculation accuracy. In order to achieve accurate measurement of the dynamic attitude in typical deck areas, it is necessary to consider the influence of ship deck deformation.
The method based on inertia matching is easy to implement and can calculate and measure the deformation angle of ships with high accuracy. However, to obtain high-precision measurements, it is necessary to identify the parameters of the dynamic deformation model. The original method of parameter identification for dynamic flexural deformation models involves processing historical data of deformation measurements. However, due to the influence of environmental factors, the parameters of dynamic models are constantly changing, so the parameters obtained from prior data analysis may not meet the actual needs [9]. Zheng et al. [10] proposed an online receding horizon identification method for nonlinear ship motion systems. The Tufts–Kumaresan method was introduced to compute unknown parameters of the dynamic flexure model using the angular increment correlation function [11]. Yu et al. [12] employed the particle swarm optimization (PSO) algorithm to form a complete online identification for ship deformation angle measurement model. The support vector regression (SVR) algorithm-based identification method was presented for the estimation of parameters in the ship dynamic model [13]. Wang et al. [14] designed an improved whale optimization algorithm (WOA) based on logic mapping to perform parameter identification.
The dynamic deformation of the deck can be described by a second-order Markov process. However, the model inevitably has a truncation error, and the dynamic deformation estimated by the filter will also be subject to estimation errors due to the truncation error. In addition, the model is not driven by ideal Gaussian white noise due to the impact interference of the wave environment. Therefore, based on the above analysis, in the process of measuring the dynamic attitude of the deck, it is necessary to consider the effects of model truncation error and non-Gaussian noise. A multi-model adaptive hull deformation estimation algorithm was proposed to address the uncertainty of system parameters and the statistical characteristics of measurement noise [15]. Ren et al. [16] proposed a dynamic filter algorithm combining wavelet and a Kalman filter for dynamic inertial measurement, in which the wavelet method is used to remove the outliers in the acquired data, and the Kalman filter effectively reduces the influence of white noise. An adaptive dynamic particle swarm optimization (ADPSO) algorithm and a bidirectional long short-term memory (BiLSTM) neural network-based prediction model were proposed and applied to accurately predicting ship motion attitude [17]. The neural network Kalman filter (NNKF) was introduced for real-time hull deformation calculation to further reduce the nonlinear error of the system. This method can accurately measure the hull deformation in real time and effectively suppress the nonlinear error caused by the large misalignment angle [18]. Traditional Kalman filtering regards model errors as process noise, which has limitations in its application. To circumvent this problem, the filtering estimator based on MPF was derived and treated as a suitable estimator for random measurement processes [19]. MPF realizes the error estimation of the system model with good robust stability [20,21]. In [22], model predictive filter technology was applied in the transfer alignment, which can estimate deformation and compensate alignment errors. A novel model predictive-based unscented Kalman filter (MP-UKF) was proposed to predict the dynamic model error persistently and correct the filtering procedure of UKF online [23].
However, traditional inertial measurement methods from our previous work cannot be directly applied to the dynamic attitude measurement of a deck with deformation because of the following reasons: (1) The method based on inertia matching can accurately calculate and measure the deformation angle of ships, and the most crucial aspect is how to accurately identify the parameters of the dynamic deformation model. The existing methods cannot meet practical needs. (2) The traditional Kalman filter considers the model error as process noise, and models that approximate dynamic deformation using second-order Markov processes inevitably have truncation errors. In addition, due to the impact interference of the wave environment, the model is not driven by ideal Gaussian white noise, and the criterion of minimum mean square error based on Gaussian noise cannot provide robust estimation performance. (3) The traditional linear error model for small misalignment angles cannot be applied to situations with large misalignment angles. However, nonlinear models based on large misalignment angles also involve excessive redundant calculations in the nonlinear filtering process. Therefore, it is necessary to study the system error model that is insensitive to the magnitude of misalignment angles.
The remaining structure of this paper is organized as follows: The dynamic deformation model of the deck is established in Section 2. The online parameter identification for the dynamic deformation model based on angle increment difference is described in Section 3. The dynamic attitude measurement method for the deck based on model predictive filtering is presented in Section 4. The simulation analysis and shipborne experiment verification of the proposed method are carried out in Section 5 and Section 6, respectively. Finally, Section 7 summarizes the main conclusions.
The contributions of this paper can be summarized as follows: (1) The related vector machine algorithm was applied to the parameter identification of the dynamic deformation model based on angle increment difference. (2) Considering the truncation error of the model and the influence of non-Gaussian noise, this paper designed a maximum correlation entropy robust filter based on MPF. In order to estimate and compensate for filtering model errors, the maximum correlation entropy was introduced as an optimization criterion to achieve robust estimation of deck dynamic deformation angle. (3) This method constructed the filter model based on Rodrigues parameters in an inertial frame, which is suitable for any misalignment angle situation, eliminating the calculation error of the sub-inertial navigation system. (4) All experimental data were obtained from the navigation experimental platform. Both simulation analysis and online verification were conducted to validate the effectiveness of the proposed method.

2. Dynamic Deformation Analysis and Modeling of the Deck

The deck undergoes deformation under the influence of the external environment. The deformation angle of deck can generally be divided into static deformation angle and dynamic deformation angle as follows:
φ = Φ + θ ( t ) ,
where Φ is the static deformation angle of the deck and θ ( t ) is the dynamic deformation angle of the deck over time. The static deformation angle is considered to be constant or to slowly change, and the model construction of such deformation is relatively simple. In contrast, dynamic deformation shows higher frequency and greater complexity. The following models were used to describe the static and dynamic deformation angles.
(1)
The static deformation of the deck
Under the combined action of multiple external forces, the deck undergoes bending deformation, resulting in static deformation angles. In fact, the deck deformation measurement system measures the real-time relative misalignment angle between the IMU in the deck area and the ship’s main inertial navigation IMU, which includes the static and dynamic deformation of the ship and the system installation error angle. The installation error angle between two sets of IMUs can be considered a constant value, and its statistical characteristics are similar to the static deformation angle of the deck. Consequently, the static deformation angle and the system installation error angle were combined into the static deformation angle. Generally, the static deformation angle changes slowly and can be approximated as a constant value as described in Equation (2).
Φ ˙ = 0 .
(2)
The dynamic deformation of the deck
Compared with the static deformation, the dynamic deformation changes more rapidly. The dynamic deformation is primarily caused by the wind and waves, which can be regarded as a random variable under random external disturbances, akin to the Markov process for white noise. Therefore, using second-order Markov process models as the ship’s dynamic deformation model is usually suitable for most scenarios. The deformation models in three directions are considered independent of each other, and the dynamic deformation angle can be given as follows [24,25]:
θ ¨ + 2 μ θ ˙ + b 2 θ = w f ,
where θ = [ θ x θ y θ z ] T is the deflection angle; μ = [ μ x μ y μ z ] T is the rate of deflection deformation angle; w f and Var ( w f ) = Q f = 4 β 3 σ 2 ( β = 2.146 / τ ) are the white noise and its variance; τ = [ τ x τ y τ z ] T is the relevant time; and σ = [ σ x σ y σ z ] T is the variance of the deflection deformation angle.
However, the lever arm between the main and sub-IMUs will have an impact on the inertial measurement information. Specifically, the fixed arm length between the two IMUs remains unchanged. Due to the influence of deck dynamic deflection deformation, it can lead to dynamic changes in lever-arm length, resulting in dynamic lever-arm error. Figure 1 shows the schematic diagram of the length variation of the dynamic lever arm.
Assuming the fixed lever-arm vector is r 0 = r x r y r z T and the deflection angle vector is θ = θ x θ y θ z T in the main inertial vehicle coordinate system (m system), we can obtain
δ r z 1 = r y tan θ x r y θ x δ r y 1 = r z tan θ x r z θ x .
Similarly, the change in the component of the lever arm caused by the dynamic lever -arm θ y and θ z is presented in the following Equation
δ r x 1 = r z θ y δ r z 1 = r x θ y δ r x 2 = r y θ z δ r y 2 = r x θ z .
The dynamic lever r m can be described as follows:
r m = r x + δ r x 1 + δ r x 2 r y + δ r y 1 + δ r y 2 r y + δ r z 1 + δ r z 2 T = r x + r z θ y r y θ z r y + r x θ z r z θ x r z + r y θ x r x θ y T .
Equation (7) gives the first and second derivatives of the dynamic lever-arm vector.
r ˙ m = r z θ ˙ y r y θ ˙ z   r x θ ˙ z r z θ ˙ x   r y θ ˙ x r x θ ˙ y T r ¨ m = r z θ ¨ y r y θ ¨ z   r x θ ¨ z r z θ ¨ x   r y θ ¨ x r x θ ¨ y T .

3. Online Parameter Identification of the Dynamic Deformation Model Based on Angle Increment Difference

To achieve a high-precision dynamic deformation model, a dynamic deflection deformation model parameter identification method based on the difference in angular increment is proposed, in which an optimization model was constructed based on the autocorrelation function of the angle increment difference between two sets of IMUs, and the related vector machine algorithm was applied for online parameter estimation. Finally, the estimated dynamic deflection parameters were applied to the filtering model based on inertia matching to estimate the deformation angle. The online identification process of dynamic deformation model parameters is shown in Figure 2.

3.1. Derivation and Analysis of Correlation Functions Based on Angular Increment Difference

The dynamic deflection angle information is contained in the angle increment output difference of the IMU and presented in the form of exponentially decaying sine signals (EDS). In order to obtain accurate model parameters for dynamic deformation angle, the relationship between the output of the gyroscope and the dynamic deformation angle is crucial. To estimate the model parameters using online estimation methods, the relationship between the output of the gyroscope and the dynamic deformation angle needs to be derived. The angular velocity is output by the gyroscope, and the relationship between angular increment and dynamic deformation angle parameters can be calculated using the output data of the main and sub-inertial navigation system gyroscopes. Specifically, it is assumed that both the main IMU and sub-IMU of a ship can be sensitive to the projection of the ship’s inertial angular velocity on the o x y z frame and o x y z frame, represented by ω ib b and ω ib x . According to [26], the angular velocity matching relationship is as follows:
ω ib s = C b s ( φ ) ω ib b + θ ˙ .
In the case of small deformation angles, we can obtain
C s b 1 φ z φ y φ z 1 φ x φ y φ x 1 = I + ( φ × ) .
The static deformation angle is treated as a constant value. When the deformation is regarded as a linear model, the difference in angular velocity between the main and sub-IMU is expressed as follows:
Δ ω = ω i b b ω i b s = ( φ × ) ω i b b = ω ib b × φ θ ˙ ,
where
ω i b b × = 0 ω i h z b ω i b y b ω i h z b 0 ω i b x b ω i b y b ω i h x b 0 .
Integrating both sides of Equation (10) in the k-th sampling simultaneously yields
t k t 0 t k ω ib b   d t t k t 0 t k ω i b s d t = t k t 0 t k ω ib b × φ d t t k t 0 t k θ ˙ d t ,
where t 0 is the sampling period.
By sorting out Equation (12), two sets of angular increment matching equations for inertial navigation within k sampling cycles can be obtained.
δ Θ 1 k δ Θ 2 k = δ Θ 1 k × φ k + φ k 1 2 θ k θ k 1 .
For Equation (13), δ Θ 1 k and δ Θ 2 k represent the angular increments of IMU1 and IMU2 within a single sampling period, respectively. By φ k = Φ + θ k , Equation (13) can be rewritten as follows:
Δ Θ k = δ Θ 1 k δ Θ 2 k = δ Θ 1 k × Φ + θ k + θ k 1 2 θ k θ k 1 = δ Θ 1 k × Φ + δ Θ 1 k × 2 I 3 θ k + δ Θ 1 k × 2 + I 3 θ k 1 .
By estimating the static deformation angle Φ , the term δ Θ 1 k × Φ in Equation (14) is compensated online and we obtain
Δ Θ k δ Θ 1 k × 2 I 3 θ k + δ Θ 1 k × 2 + I 3 θ k 1 .
Further, Equation (15) can be approximated as follows:
Δ Θ k θ k + θ k 1 .
When the correlation time τ ≥ 0, the correlation function of the dynamic deformation angle can be expressed as follows:
R θ i ( τ ) = λ i σ i 2 μ i 2 + λ i 2 exp μ i τ sin λ i τ + ε i , i = x , y , z ,
where ε i = arctan λ i / μ i .
The following definition is made in terms of the above Equation.
R ˜ θ i ( τ ) = λ i σ i 2 μ i 2 + λ i 2 exp μ i τ + j λ i τ + ε i , i = x , y , z .
Combining Equation (17) and Equation (18) yields
R θ i ( τ ) = 1 2 j R ˜ θ i ( τ ) R ˜ θ i * ( τ ) , i = x , y , z .
According to the operation rule of the correlation function, the correlation function of discrete angular increments Δ Θ k can be expressed as
R Δ Θ ( n ) = Δ Θ k , Δ Θ k + n = θ k , θ k + n + θ k 1 , θ k + n 1 θ k 1 , θ k + n θ k , θ k + n 1 = 2 R θ ( n ) R θ ( n 1 ) R θ ( n + 1 ) .
Meanwhile, similar to the dynamic deformation angle model, the correlation function of Δ Θ k can be expressed as follows:
R Δ Θ ( n ) = 1 2 j R Δ Θ ( n ) R Δ Θ * ( n ) ,
R ˜ Δ Θ ( n ) = 2 R ˜ θ ( n ) R ˜ θ ( n 1 ) R ˜ θ ( n + 1 ) = 2 λ i σ i 2 μ i 2 + λ i 2 exp j ε i 1 cosh μ i + j λ i exp μ i + j λ i n .
Equation (22) can be further organized to obtain
R ˜ Δ Θ ( n ) = C i exp μ i + j λ i n ,
where
C i = 2 λ i σ i 2 μ i 2 + λ i 2 exp j arctan λ i / μ i 1 cos h μ i + j λ i .
It can be seen that Equation (23) has a similar form to the exponential sinusoidal attenuation signal. If the parameters μ i , λ i , and C i are solved by Equation (21), the first two terms present the attenuation factor and dominant frequency of the dynamic deflection model. By inputting the above parameters into Equation (24), σ i 2 can be calculated.

3.2. Online Identification of Dynamic Deformation Model Parameters Based on the RVM Algorithm

RVM is based on the structure of prior parameters and the automatic relevance determination theory to establish a sparse model under the Bayesian framework. Compared to SVM, the computational complexity of the kernel function is significantly reduced, and there is no need for the selected kernel function to satisfy the Mercer condition [27].
It is given that N samples are x 1 , y 1 , , x N , y N R n × R m , and the regression function is given by
y ( x ; ω ) = i = 1 N ω i K ( x , x i ) + ω 0 ,
where ω = ω 1 , ω 2 , , ω N T is the weight vector; ω 0 is the threshold value; and K ( x , x i ) is the kernel function.
Because the target sample y i is independent, the likelihood function can be denoted as follows:
P ( y | ω , σ 2 ) = ( 2 π σ 2 ) N / 2 exp y Φ ω 2 2 σ 2 ,
where Φ is a N × ( N + 1 ) order structural matrix composed of the kernel function.
According to the Bayesian framework, the prior distribution can be described as
P ( ω | α ) = ( 2 π σ 2 ) N / 2 i = 1 N α i 1 / 2 exp α i ω i 2 2 ,
where α = α 1 , α 2 , , α N T is the super parameter vector.
By Bayesian inference, the posterior distribution satisfies the following Equation
P ( y | ω , α , σ 2 ) ~ N ( μ , Σ ) ,
where the posterior mean μ and variance can be calculated as follows:
= ( σ 2 Φ T Φ + A ) 1 μ = σ 2 Φ T y ,
where A = d i a g ( α 1 , α 2 , , α N ) .
By the edge integration of weights, the likelihood distribution of output is obtained by
P ( y | α , σ 2 ) = P ( y | ω , σ 2 ) P ( ω | α ) d ω .
We can obtain the edge likelihood distribution of hyperparameter, that is
P ( y | α , σ 2 ) ~ N ( 0 , C ) ,
where C = σ 2 I + Φ A 1 Φ T . I is an identity matrix.
Updated α i n e w and ( σ 2 ) n e w are determined by Equation (32).
α i n e w = γ i μ i 2 ( σ 2 ) n e w = y Φ μ 2 N i = 0 N γ i ,
where μ i is the i-th element value of μ ; and γ i = 1 α i N i i and N i i are the i-th diagonal value of .
μ and are updated by Equation (32) until the convergence condition presented as Equation (33) is met.
L ( ω ) = 1 2 σ 2 i = 1 N y i ω T Φ ( x i ) 2 + i = 0 N log ω i ,
where the first term is the sum of measurement errors, and the second term is the regularization term. Figure 3 shows the flowchart of the RVM algorithm.
Using the RVM algorithm for dynamic deformation model parameter identification involves deriving the autocorrelation function form of the difference in the output angle increment of two sets of IMUs. In the process of simplifying the angle increment difference Δ Θ k , online compensation is implemented for δ Θ 1 k × Φ . The autocorrelation function of the angle increment difference is obtained to identify the parameters of the dynamic deformation model by comparing Equations (21) and (23). The related vector machine algorithm is then selected to assist in establishing the algorithm model using the general form of exponentially decaying sine signals, and the parameters are reasonably set to complete the online parameter identification. As shown above, the specific process of online identification of dynamic deformation model parameters is as follows:
(1)
Calculate the angular increment output difference Δ Θ k of two sets of IMUs in each sampling cycle and compensate for δ Θ 1 k × Φ . The first static deformation angle Φ can be roughly determined by the coarse alignment, and then estimated by the deformation algorithm.
(2)
Solve the numerical value R Δ Θ ( n ) of the autocorrelation function Δ Θ k .
(3)
Substitute R Δ Θ ( n ) into the RVM algorithm and use Equation (33) for parameter identification.
(4)
Apply the identified parameters to the angular velocity matching algorithm for deformation angle estimation, substitute the estimated deformation angle back in step (1), and proceed to the next step of calculation.

4. Dynamic Attitude Measurement Method for Deck Based on Model Predictive Filtering

To solve the problems of truncation error and non-Gaussian noise in the dynamic attitude deformation angle model of the deck, this paper proposes a dynamic attitude measurement method for the deck based on model prediction filtering. By designing a robust maximum correlation entropy filter based on model prediction, estimating and compensating for the filtering model error, a robust estimation of the deck deflection deformation is achieved.

4.1. Model Error Analysis and Prediction Estimation

Model predictive filtering is an advanced filtering method proposed in recent years that overcomes the limitation of traditional Kalman filtering that assumes model errors as process noise processing. The biggest advantage of model predictive filtering is its fast calculation speed. It estimates observable low-state variables as model errors and uses observations to correct uncertain errors. In addition, since the model predicts the filtering equation through analytical expressions and does not have an iterative inversion process, it can reduce system dimensionality, accelerate calculation speed, and improve real-time performance while ensuring high system accuracy [22,23].
Specifically, the principle of predictive filtering involves using the predicted output of the filter to track the actual observed output in real time, thereby estimating the model error of the system. Firstly, the system state space differential model is derived from differential equations and is as follows:
x ˙ ( t ) = F ( t ) x ( t ) + G w ( t ) z ( t ) = H ( t ) x ( t ) + v ( t ) ,
where x ( t ) is the system status; F ( t ) is the state transition matrix; w ( t ) is the system noise vector; H ( t ) is the measure matrix; z ( t ) is the external measurement vector; and v ( t ) is the measurement noise. Representing the state estimation value and corresponding measurement prediction value with x ^ ( t ) and z ^ ( t ) , respectively, and the measurement model is obtained by the Taylor expansion
z ^ ( t + Δ t ) z ^ ( t ) + Z ( x ^ ( t ) , Δ t ) + Λ ( Δ t ) S ( x ^ ( t ) ) D ( t ) ,
where
z ^ ( x ^ ( t ) , Δ t ) = k = 1 p 1 Δ t k k ! L f k H 1 k = 1 p n Δ t k k ! L f k H m T ,
Λ ( Δ t ) = diag Δ t p 1 p 1 ! Δ t p m p m ! ,
S ( x ^ ( t ) ) = L g 1 L f p 1 1 H 1 L g q L f p 1 1 H 1 L g 1 L f p m 1 H m L g q L f p m 1 H m m × q ,
where p i ( i = 1 , , m ) is the lowest order that appears in the differentiation of the i-th component; H i ( t ) ; L f k H i is the k-th Lie derivative; and g i ( i = 1 , , q ) is the i-th column in G ( t ) . The definition of the predictive filtering performance indicator function is
J = 1 2 z k + 1 z ^ k + 1 T R 1 z k + 1 z ^ k + 1 + 1 2 D k T A D k = min ,
where A is the positive semi-definite weighting matrix for model error.
The estimated model error D k that satisfies the performance metric of Equation (39) within time t k , t k + 1 can be designed as
D ^ k = Λ ( Δ t ) S x ^ k T R 1 Λ ( Δ t ) S x ^ + A 1             × Λ ( Δ t ) S x ^ T R 1 z ^ k + 1 + Z x ^ , Δ t z k + 1 .
According to Equation (40), D ^ k can be calculated and compensated in the state transition model.

4.2. Design of Maximum Correlation Entropy Robust Filter Based on Model Prediction

In the model of the dynamic deformation angle of the deck mentioned above, it is assumed that the noise follows a Gaussian probability density function (PDF). However, due to the impact of wave environments and other external interferences from the ocean environment, the model is not driven by ideal Gaussian white noise and exhibits a non-Gaussian distribution. Compared to the minimum mean square error (MMSE) criterion based on Gaussian noise, the maximum correlation entropy criterion offers better robust estimation performance in non-Gaussian situations by capturing high-order moments of information. Therefore, this paper combines model prediction filtering and employs MCE as the optimal estimation criterion instead of MMSE to achieve robust and accurate estimation of the dynamic deformation angle of the deck.
The model error vector estimated by Equation (40) is used to compensate for the state prior mean as follows:
x k | k 1 = F k x k 1 + G k D ^ k ,
where G k is the transfer matrix for model error.
Furthermore, the covariance matrix of the state prior mean is calculated as follows:
P k | k 1 = F k P k 1 F k T + G k Q k G k T .
The measurement update process is based on maximizing the correlation entropy between state prior and measurement information, rather than minimizing the mean square error of state estimation. Specifically, the measurement prediction information is calculated as follows:
z k | k 1 = H k x k | k 1 .
Using Equation (44), the cost function is built as
J MCC x k = G σ x k x k | k 1 P k k 1 1 2 + G σ z k z k | k 1 H k x k x k | k 1 R k 1 2 ,
where A 2 = ( ) T A ( ) is the square of the weighted Mahalanobis distance with A as the weight matrix.
The Gaussian kernel function is defined as follows:
G σ ( A 2 ) = exp A 2 2 σ 2 ,
where σ > 0 is the kernel width and the optimal estimation based on maximum correlation entropy is obtained through the following objective function:
x ^ k = arg max x k J MCC x k .
The optimization problem is solved by the following Equation:
J MCC x k x k = 2 P k k 1 1 x k x ^ k k 1 2 L k H k T R k 1 z k z ^ k k 1 H k x k x ^ k k 1 = 0 ,
where
L k = G σ z k z ^ k k 1 H k x k x ^ k k 1 R k 1 2 G σ x k x ^ k k 1 P k k 1 1 2 .
By reorganizing Equation (47), we obtain
x k = P k k 1 1 + L k H k T R k 1 H k 1 L k z H k T R k 1 z k z ^ k k 1 H k x ^ k k 1 + P k k 1 1 x k .
Because the state–transition model in inertial measurement systems is constructed accurately, the prior information about the state is accurate. In addition, a correlation entropy matrix was constructed to replace the correlation entropy scalars of each measurement information. Therefore, through further sorting, the mean and covariance matrix of the posterior estimation of the state can be obtained as follows:
x ^ k = x ^ k | k 1 + K k z k z ^ k | k 1 ,
P k = P k k 1 K k P z z , k k 1 ( K k ) T ,
where
K k = P x z , k k 1 P z z , k k 1 1 = P k k 1 H k T C k H k P k k 1 H k T C k + R k 1 ,
C k = diag G σ k z ˜ 1 , k R 1 , k 1 2 , , G σ k z ˜ m , k R m , k 1 2 .

4.3. Dynamic Measurement Model Based on Rodrigues Parameters in Inertial Frames

Due to the presence of the lever arm between the main IMU and sub-IMU, the dynamic deformation angle will introduce an error in the dynamic lever arm. Therefore, it was necessary to consider the influence of the dynamic lever arm. The measurement relationship between the main IMU and sub-IMU is provided in the inertial frame. The projection relationship of the specific force sensitive to the main IMU and sub-IMU under the influence of the lever-arm effect in the geocentric inertial frame is presented as follows:
f s i = f m i + C m i [ ω i m m × ( ω i m m × r m ) + ω ˙ i m m × r m + 2 ω i m m × r ˙ m + r ¨ m ] ,
where C m i is the coordinate transformation matrix from the main inertial coordinate system to the geocentric inertial coordinate system; ω i m m is the output of the main IMU’s gyroscope; ω i m m × ω i m m × r m is the acceleration caused by the rigid lever arm; and 2 [ ω i m m × ] r ˙ m + r ¨ m is the acceleration caused by the dynamic lever arm. According to Equation (54), after compensating for the acceleration of the rigid lever arm, we obtain
f s i = f m c i + C m i f r d f r d = { [ ( ω i m m × ) ( ω i m m × ) + ( ω ˙ i m m × ) ] L 0 + L 1 } θ + [ 2 ( ω i m m × ) L 0 + L 2 ] μ + L 0 w f ,
where f m c i is the specific force information of main IMU after compensating for the acceleration of the rigid lever arm; f r d is the dynamic arm acceleration caused by deflection deformation; and L 0 , L 1 and L 2 are represented as follows:
L 0 = 0 r z r y r z 0 r x r y r x 0 L 1 = 0 β y 2 r z β z 2 r y β x 2 r z 0 β z 2 r x β x 2 r y β y 2 r x 0 L 2 = 0 2 β y r z 2 β z r y 2 β x r z 0 2 β z r x 2 β x r y 2 β y r x 0 .
Considering the influence of the lever-arm effect, the angular velocity relationship measured by the main IMU’s and sub-IMU’s gyroscope in the inertial frame is defined as follows:
ω i s i = ω i m i + C m i ω m s m = ω i m i + C m i μ ,
where ω i s i is the projection of the sub-IMU’s gyroscope in the inertial frame; ω i m i is the value of the main IMU’s gyroscope; ω m s m is the rotational angular velocity of the sub-IMU relative to the main IMU; and C m i ω m s m can be regarded as the deflection deformation angular velocity.
The main inertial vehicle coordinate system (m) and sub-inertial vehicle coordinate system (s) were fixed at the initial moment of transfer alignment, and the inertial solidification coordinate system, i m and i s , could be obtained, respectively. The coordinate transformation matrix C m i m between them is a time-varying matrix that can be calculated from the output of the main IMU’s gyroscope ω i m m . According to the directional cosine matrix differential equation, we can obtain
C ˙ m i m = C m i m ( ω i m m × ) C m i m ( 0 ) = I 3 .
Similarly, the transformation matrix C s i ^ s can be represented as follows:
C ˙ s i ^ s = C s i ^ s ( ω ^ i s s × ) C s i ^ s ( 0 ) = I 3 ,
where i ^ s is the computational inertial solidification coordinate system and ω ^ i s s is the output of sub-IMU’s gyroscope.
For the short-term alignment process, this tracking error ϕ s i s can be considered a small amount, and there are
C s i ^ s = C i s i ^ s C s i s = [ I 3 ( ϕ s i s × ) ] C s i s ϕ ˙ s i s = C s i s [ ( ε s + w g s ) ] ,
where ε s is the constant zero bias of the sub-IMU’s gyroscope and w g s is the random noise of the sub-IMU’s gyroscope.
After accurately fixing the length of the arm of the main and sub-inertial navigation systems, the rigid arm acceleration compensation was applied to the main inertial navigation system. The specific force after compensation is denoted as f m c m . f m c m can be projected to the i m system by C m i m as follows:
f m c i m = C m i m f m c m .
Similarly, the output of sub-IMU’s accelerometer can be transformed and described in i ^ s system as
f ^ s i ^ s = C s i ^ s f ^ s s = [ I 3 ( ϕ s i s × ) ] C s i s ( f s s + s + w a s ) ,
where s is the constant zero bias of the sub-IMU’s accelerometer and w a s is the random noise of the sub-IMU’s accelerometer.
If Equation (62) is expanded and high-order small quantities are ignored, we can obtain
f ^ s i ^ s = f s i s + f ^ s i ^ s × ϕ s i s + C s i ^ s s + C s i ^ s w a s .
Combining Equations (55), (61), and (63) yields
f ^ s i ^ s f ^ s i ^ s × ϕ s i s C s i ^ s s C s i ^ s w a s = C i m i s ( f m c i m + C m i m f r d ) ,
where C i m i s is a constant matrix and represents the directional cosine matrix corresponding to the installation error angle between the main and sub-inertial vehicle coordinate systems at the initial moment of alignment.
Similarly, for the output of the main and sub-IMU’s gyroscope, we obtain
ω ^ s i ^ s ω ^ s i ^ s × ϕ s i s C s i ^ s ε s C s i ^ s w g s = C i m i s ( ω i m i m + C m i m μ ) ,
where
ω ^ s i ^ s = C s i ^ s ω ^ i s s ω i m i m = C m i m ω i m m .
C i m i s can be represented by Rodrigues vectors l i m i s . According to the Cayley transformation, the relationship between the two satisfies the following Equation
C i m i s = [ I 3 + ( l i m i s × ) ] 1 [ I 3 ( l i m i s × ) ] .
By substituting Equation (67) into Equations (64) and (65), it can be concluded that
f ^ s i ^ s f m c i m = ( f ^ s i ^ s + f m c i m ) × l i m i s + [ I 3 + ( l i m i s × ) ] δ f s i ^ s + [ I 3 ( l i m i s × ) ] C m i m f r d ,
ω ^ s i ^ s ω i m i m = ( ω ^ s i ^ s + ω i m i m ) × l i m i s + [ I 3 + ( l i m i s × ) ] δ ω s i ^ s + [ I 3 ( l i m i s × ) ] C m i m μ ,
where δ f s i ^ s = f ^ s i ^ s × ϕ s i s + C s i ^ s s + C s i ^ s w a s and δ ω s i ^ s = ω ^ s i ^ s × ϕ s i s + C s i ^ s ε s + C s i ^ s w g s .
The right-hand side of Equations (68) and (69) contain the Rodrigues vectors corresponding to the installation error angles between the main IMU and sub-IMU. Therefore, the specific force ( f ^ s i ^ s , f m c i m ) and angular velocity ( ω ^ s i ^ s , ω i m i m ) of the main IMU and sub-IMU can be used as matching variables to construct the corresponding matching equations. Since there is no small angle assumption for the installation error angle during the derivation process, this matching equation is suitable for the transfer alignment process with a large misalignment angle.
The specific force and angular velocity outputs of the sub-IMU contain random noise disturbances, such as dynamic deflection deformation and measurement noise, which affect the matching accuracy in Equations (68) and (69). The smoothing characteristics of the integration process can attenuate random noise, improve the signal-to-noise ratio of measurement information, and improve the performance of misalignment angle estimation. However, although integral operation can enhance the anti-interference ability against random noise to a certain extent, it also complicates the transfer process from measurement information to state variables. Therefore, only the specific forces of the main IMU and sub-IMU in the inertial solidification coordinate system were selected for integration.
Specifically, integrating Equation (68) yields
V s i s V m i m = ( V s i s + V m i m ) × l i m i s + [ I 3 + ( l i m i s × ) ] δ V s i s + [ I 3 ( l i m i s × ) ] V r i m ,
where
V s i ^ s = 0 t f ^ s i ^ s ( τ ) d τ = 0 t C s i ^ s ( τ ) f ^ s i s ( τ ) d τ V m i m = 0 t f m c i m ( τ ) d τ = 0 t C m i m ( τ ) f m c m ( τ ) d τ δ V s i ^ s = 0 t ( f ^ s i ^ s × ϕ s i s + C s i ^ s s + C s i ^ s w a s ) d τ V r i m = 0 t C m i m f r d d τ .
Because C i m i s is a constant matrix, the Rodrigues parameters corresponding to C i m i s satisfy the following differential Equation:
l ˙ i m i s = 0 .
Except for the constant drift of sub-inertial gyroscopes ε s and the constant bias s of accelerometers, the integration of dynamic arm acceleration V r i m , deflection angle θ , and angular velocity μ in system i m were selected as state variables. Equation (73) gives the definition of state vector.
X = [ ( l i m i s ) T ( δ V s i s ) T ( ϕ s i s ) T ( ε s ) T ( s ) T ( V r i m ) T ( θ ) T ( μ ) T ] T .
Therefore, based on the above derivation, the system state differential equation was constructed as follows:
X ˙ = F X + G w ,
where the process noise w = [ ( w g s ) T ( w a s ) T ( w f ) T ] T , and the state transition matrix F is represented as
F = 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 ( f ^ s i ^ s × ) 0 3 × 3 C s i ^ s 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 C s i ^ s 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 C m i m A r 1 C m i m A r 2 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 I 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 A w 1 A w 2 ,
where
A r 1 = [ ( ω i m m × ) ( ω i m m × ) + ( ω ˙ i m m × ) ] L 0 + L 1 A r 2 = 2 ( ω i m m × ) L 0 + L 2 A w 1 = diag ( β 2 ) A w 2 = diag ( 2 β ) .
The process noise driving matrix G is represented as follows:
G = 0 3 × 3 0 3 × 3 ( C s i ^ s ) T 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 ( C s i ^ s ) T 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 ( C m i m L 0 ) T 0 3 × 3 I 3 × 3 T .
According to Equations (69) and (70), the measurement equation can be constructed as follows:
Z = V s i s V m i m ω ^ s i ^ s ω i m i m = h ( X ) + w k = V s i s + V m i m ω ^ s i ^ s + ω i m i m × l i m i s + I 3 + ( l i m i s × ) 0 3 × 3 0 3 × 3 I 3 + ( l i m i s × ) δ V s i s ω ^ s i ^ s × ϕ s i s + C s i ^ s ε s + I 3 ( l i m i s × ) 0 3 × 3 0 3 × 3 I 3 ( l i m i s × ) V r i m C m i m μ + w k ,
where w k is the measurement noise caused by the integration of residual errors and inertial device errors. Compared with the ship’s deformation angle measurement model based on angular velocity matching, this model does not require the initial parameters binding and strapdown calculation of the main IMU. It reduces the computational workload while suppressing the influence of compensating deflection deformation and dynamic arm errors, thereby enhancing the real-time performance of the calculation.

5. Simulation Experiment and Analysis

To verify the proposed method, we simulated and generated the ship’s motion attitude, IMU data, dynamic deformation angle, and other data through the following simulation conditions.
(1)
The initial latitude was set to 45° and the initial longitude was set to 126°. The motion state of the ship was to move forward with a constant speed of 20 knots. The initial roll angle was 0°, the initial pitch angle was also 0°, and the initial heading angle was 45°. The ship’s sway motion satisfies Equation (79)
θ = 10 ° sin ( π 4 t ) γ = 10 ° sin ( π 4 t ) γ = 45 ° + 10 ° sin ( π 5 t ) .
(2)
The constant the gyroscope drift was set to 0.01°/h, the angle random walk to 0.001°/ h , the constant zero deviation of the accelerometer to 50 ug, and the acceleration random walk to 10 ug/ h . The installation error angle of IMU is 5″, and the scale factor error was 5 ppm. The calculation frequency was set to 100 Hz, and the initial horizontal misalignment angle and azimuth misalignment angle were 1′, and 5′, respectively.
(3)
The static deformation angle was set to a constant value [10′ 10′ 30′], and the dynamic deformation angle was generated using a second-order Markov model as follows:
θ i ( k + 2 ) = 2 1 μ i T θ i ( k + 1 ) + 2 μ i T b 2 T 2 θ i ( k ) + 2 b D i μ i T 3 w ( k ) ,
where w ( k ) is zero-mean Gaussian white noise. The variance of the deflection angle was set to σ = [ 6 6 10 ] 2 , and the time related to the deflection angle was set to τ 1 = [ 10 s 10 s 30 s ] T . The fixed arm length was set to r b = [ 10 m 10 m 2 m ] T .
(4)
Compare the proposed improved deck dynamic attitude measurement algorithm (IDAM) based on model prediction with the following common existing methods to verify the effectiveness of the proposed algorithm regarding the support vector regression-based algorithm for deck dynamic attitude measurement (SVRAM) [12] and the whale optimization-based algorithm for deck dynamic attitude measurement (WOAM) [13]. The root mean square error (RMSE) was used as a quantitative indicator to evaluate the accuracy of the different measurement methods and was defined as follows:
RMSE t = 1 T k = 1 T x k x ^ k 2 ,
where x ^ k and x k represent the estimated state value and the true value of the state.
The CPU configuration of simulation computer was an Intel Core i5-1135G7 (2.4 GHz) and the simulations were run using Microsoft Visual Studio 10.0 software.
The estimation of deflection angle using different methods is shown in Figure 4. Figure 5 and Figure 6 give the attitude estimation errors and its RMSE values with the different methods, respectively.
It can be concluded from Figure 4 that the proposed algorithm can track and estimate the dynamic deformation angle more effectively compared to the SVRAM and WOAM methods, especially for the estimation of the z-axis dynamic deformation angle. Furthermore, as shown in Figure 5 and Figure 6, the proposed IDAM method significantly improves the accuracy of deck attitude measurement with accurate estimation of dynamic deformation angles and effective implementation of model prediction filtering. For the z-axis attitude, the estimated RMSE of the IDAM was 4.260′. Compared with the 6.560′ and 5.032′ of SVRAM and WOAM methods, respectively. Thus, the estimation accuracy decreased by 35.06% and 15.34%, respectively, indicating that the proposed algorithm can effectively improve the attitude measurement accuracy under deck dynamic deformation.

6. Shipborne Experiment and Analysis

To further verify the effectiveness and reliability of the proposed method, a shipborne experiment was conducted. Figure 7a shows the overall view of the shipborne experiment platform. All experimental equipment is annotated in Figure 7b including the high-precision fiber optic inertial navigation system (PHINS), the main IMU, the sub-IMU and so on. PHINS, the main IMU, and sub-IMU were connected with the synchronous recording device and data record computer to synchronously record all navigation information.
The detailed information of data collected from the experiment is listed in Table 1. The sampling time was 10 ms.
In order to complete the precise alignment of the PHINS as quickly as possible during the experiment, approximately 20 min of maneuvering were carried out until the alignment of the PHINS was completed. After the alignment process was finished, static data were collected for 2 h. The fixed installation deviations between the testing IMUs and PHINS were determined by integrated alignment using the static data. The calculation results are shown in Table 2.
As shown in Table 2, installation deviation1 represents the installation deviation between the main IMU and PHINS-1 on the left side in Figure 7b, and installation deviation2 denotes the installation deviation between the sub-IMU and PHINS-2 on the right side in Figure 7b.
Figure 8 gives the navigation trajectory of the experimental platform.
Figure 9 gives the attitude error of the deck under different methods, and its RMSE is shown in Figure 10.
From Figure 9, the following conclusions can be drawn: (1) Compared with the SVRAM method, the attitude errors of the WOAM and proposed IDAM methods were significantly reduced; (2) in terms of horizontal attitude error, the difference between the WOAM methods and the proposed IDAM methods was minimal; (3) in terms of heading error, the WOAM method had a large error, around 150 s, and the proposed IDAM method had better performance. Finally, the obtained conclusions are consistent with the simulation results.
It can be seen from Figure 10 that different methods have different attitude measurement performance. For the method based on SVRAM, the values of RMSE for roll, pitch, and heading are 5.155′, 5.098′, and 25.143′, respectively. For the method based on WOAM, the values of RMSE for roll, pitch, and heading are reduced to 0.828′, 1.914′, and 10.735′, respectively. Compared with SVRAM, the accuracy is improved by 83.94%, 62.46%, and 57.30%, respectively. Furthermore, the proposed IDAM achieves high-accuracy dynamic angle parameters, and the model prediction method further improves the estimation effect. Therefore, the estimated RMSE of roll, pitch, and heading angles are reduced to 0.797′, 1.852′, and 9.589′, which are 3.74%, 3.24%, and 10.68% better than WOAM, respectively, showing better attitude estimation performance and verifying the effectiveness of the proposed algorithm.

7. Conclusions

In order to address the difficulty of achieving accurate attitude information of typical deck areas, firstly, this paper carried out an online estimation of dynamic deformation model parameters by using the RVM algorithm to adapt to the complex changes in dynamic deformation of the deck. Secondly, an MPF-based dynamic attitude measurement method using maximum correlation entropy as the optimization criterion is proposed to deal with the influence of model truncation error and non-Gaussian noise. Finally, this paper designs a dynamic measurement model based on Rodrigues parameters in an inertial frame by integrating the application of the dynamic deformation model with adaptive parameters. Meanwhile, there is no restriction on the range of the misalignment angle and no need for sub-inertial navigation calculation. The experiment results and comparison analysis demonstrate that the proposed dynamic attitude inertial measurement method has higher estimation accuracy than the previous methods. Moreover, from an algorithm adaptability point of view, the proposed method significantly improves the navigation performance of inertial navigation systems in extreme cases.
In the future, we plan to test more types of algorithms and their combinations applied to this field and conduct the verification through shipborne experiments.

Author Contributions

Conceptualization: B.Z.; methodology: B.Z. and X.X.; validation: T.W.; software: T.W.; formal analysis: B.Z.; investigation: X.X.; resources: W.G.; writing—review and editing: B.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by the National Natural Science Foundation of China under Grant 52271315 and Grant 51909048, the China Postdoctoral Science Foundation under Grant 2019T120260, the Postdoctoral Foundation of Heilongjiang Province Government Grant LBH-TZ1015, and the Postdoctoral Foundation of Heilongjiang Province Government Grant LBH-Z22188.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Acknowledgments

The authors of this paper would like to express their sincere gratitude to the editors and reviewers for their hard work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wang, J.; Zhang, T.; Jin, B.; Zhu, Y.; Tong, J. Student’s t-Based Robust Kalman Filter for a SINS/USBL Integration Navigation Strategy. IEEE Sens. J. 2020, 20, 5540–5553. [Google Scholar] [CrossRef]
  2. Dai, H.; Lu, J.; Guo, W.; Wu, G.; Wu, X. IMU based deformation estimation about the deck of large ship. Optik 2016, 127, 3535–3540. [Google Scholar] [CrossRef]
  3. Wang, Y.; Sun, F.; Zhang, Y.; Liu, H.; Min, H. Central difference particle filter applied to transfer alignment for SINS on missiles. IEEE Trans. Aerosp. Electron. Syst. 2012, 48, 375–387. [Google Scholar] [CrossRef]
  4. Liu, X.; Xu, X.; Liu, Y.; Wang, L.H. A fast and high-accuracy transfer alignment method between M/S INS for ship based on iterative calculation. Measurement 2014, 51, 297–309. [Google Scholar] [CrossRef]
  5. Yang, G.; Wang, Y.; Yang, S. Assessment approach for calculating transfer alignment accuracy of SINS on moving base. Measurement 2014, 52, 55–63. [Google Scholar] [CrossRef]
  6. Peng, Y.; Zhang, A.; Ming, F.; Wang, S. A meshfree framework for the numerical simulation of elasto-plasticity deformation of ship structure. Ocean Eng. 2019, 192, 106507. [Google Scholar] [CrossRef]
  7. Seo, C.; Jeong, B.; Kim, J.; Song, M.; Noh, J.; Lee, J. Determining the influence of ship hull deformations caused by draught change on shaft alignment application using FE analysis. Ocean Eng. 2020, 210, 107488. [Google Scholar] [CrossRef]
  8. Xu, B.; Duan, T.; Wang, Y. An inertial measurement method of ship deformation basesd on IMM filtering. Optik 2017, 140, 601–609. [Google Scholar] [CrossRef]
  9. Wu, W.; Qin, S.; Chen, S. Coupling influence of ship dynamic flexure on high accuracy transfer alignment. Int. J. Model. Identif. Control 2013, 19, 224–234. [Google Scholar] [CrossRef]
  10. Zheng, J.; Yan, D.; Yan, M.; Yi, Y.; Zhao, Y. An unscented kalman filter online identification approach for a nonlinear ship motion model using a self-navigation test. Machines 2022, 10, 312. [Google Scholar] [CrossRef]
  11. Wu, W.; Chen, S.; Qin, S. Online estimation of ship dynamic flexure model parameters for transfer alignment. Ocean Eng. 2012, 21, 1666–1678. [Google Scholar] [CrossRef]
  12. Yu, D.; Yang, G.; Xie, Z. Online Estimation of Ship Dynamic Flexure Model Parameters Based on Particle Swarm Optimization. Nav. Control. 2019, 18, 71–79. [Google Scholar]
  13. Zhu, M.; Hahn, A.; Wen, Y.; Sun, W. Optimized support vector regression algorithm-based modeling of ship dynamics. Appl. Ocean. Res. 2019, 90, 101842. [Google Scholar] [CrossRef]
  14. Wang, Y.; Zhang, Y.; Xu, D.; Miao, W. Improved whale optimization-based parameter identification algorithm for dynamic deformation of large ships. Ocean Eng. 2022, 245, 110392. [Google Scholar] [CrossRef]
  15. Wang, Y.; Zhang, Y.; Wang, K.; Wang, Z.; Chang, J.; Xu, D. Research on multi-model adaptive hull deformation measurement algorithm. In Proceedings of the 2020 IEEE/ION Position, Location and Navigation Symposium (PLANS), Portland, OR, USA, 20–23 April 2020; pp. 718–722. [Google Scholar]
  16. Ren, B.; Li, T.; Li, X. Research on Dynamic Inertial Estimation Technology for Deck Deformation of Large Ships. Sensors 2019, 19, 4167. [Google Scholar] [CrossRef] [PubMed]
  17. Zhang, G.; Tan, F.; Wu, Y. Ship Motion Attitude Prediction Based on an Adaptive Dynamic Particle Swarm Optimization Algorithm and Bidirectional LSTM Neural Network. IEEE Access 2020, 8, 90087–90098. [Google Scholar] [CrossRef]
  18. Xu, D.; Peng, X.; Zhang, G.; Hu, X. An Improved Inertial Matching Algorithm for Hull Deformation with Large Misalignment Angle. IEEE Access 2021, 9, 36634–36644. [Google Scholar] [CrossRef]
  19. Crassidis, J.; Markley, F. Predictive filtering for nonlinear systems. J. Guid. Control Dyn. 1997, 20, 566–572. [Google Scholar] [CrossRef]
  20. Chatzis, S.; Kosmopoulos, D.; Varvarigou, T. Robust Sequential Data Modeling Using an Outlier Tolerant Hidden Markov Model. IEEE Trans. Pattern Anal. Mach. Intell. 2009, 31, 1657–1669. [Google Scholar] [CrossRef]
  21. Winters-Hilt, S.; Jiang, Z. A Hidden Markov Model with Binned Duration Algorithm. IEEE Trans. Signal. Proces. 2010, 58, 948–952. [Google Scholar] [CrossRef]
  22. Mu, H.; Han, L.; Li, A. Application of Model Prediction Filter in Shipborne Weapons Initial Alignment. Ship Electron. Eng. 2021, 41, 47–50. [Google Scholar]
  23. Hu, G.; Ni, L.; Gao, B.; Zhu, X.; Wang, W.; Zhong, Y. Model Predictive Based Unscented Kalman Filter for Hypersonic Vehicle Navigation with INS/GNSS Integration. IEEE Access 2019, 8, 4814–4823. [Google Scholar] [CrossRef]
  24. Lu, J.; Xie, L.; Li, B. Applied quaternion optimization method in transfer alignment for airborne AHRS under large misalignment angle. IEEE Trans. Instrum. Meas. 2015, 65, 346–354. [Google Scholar] [CrossRef]
  25. Gong, X.; Fan, W.; Fang, J. An innovational transfer alignment method based on parameter identification UKF for airborne distributed POS. Measurement 2014, 58, 103–114. [Google Scholar] [CrossRef]
  26. Mochalov, A.; Kazantsev, A. Use of ring laser units for measurement of moving object deformations. Symmetry 2002, 4680, 85–92. [Google Scholar]
  27. Tipping, M. Sparse Bayesian learning and the relevance vector machine. J. Mach. Learn. Res. 2001, 1, 211–244. [Google Scholar]
Figure 1. Schematic diagram of changes in the dynamic lever-arm length.
Figure 1. Schematic diagram of changes in the dynamic lever-arm length.
Electronics 13 03555 g001
Figure 2. Online identification process of dynamic deformation model parameters.
Figure 2. Online identification process of dynamic deformation model parameters.
Electronics 13 03555 g002
Figure 3. The flowchart of the RVM algorithm.
Figure 3. The flowchart of the RVM algorithm.
Electronics 13 03555 g003
Figure 4. Estimation results of the deflection angle with different methods.
Figure 4. Estimation results of the deflection angle with different methods.
Electronics 13 03555 g004
Figure 5. The error of attitude error angle estimation with different methods; (a) overall view of attitude error angle estimation; (b) partial enlarged drawing of attitude error angle estimation.
Figure 5. The error of attitude error angle estimation with different methods; (a) overall view of attitude error angle estimation; (b) partial enlarged drawing of attitude error angle estimation.
Electronics 13 03555 g005
Figure 6. Attitude error angle estimation error represented by RMSE.
Figure 6. Attitude error angle estimation error represented by RMSE.
Electronics 13 03555 g006
Figure 7. Installation diagram of shipborne experiment; (a) overall view of shipborne experiment platform; (b) installation and layout of all equipment in the experiment platform.
Figure 7. Installation diagram of shipborne experiment; (a) overall view of shipborne experiment platform; (b) installation and layout of all equipment in the experiment platform.
Electronics 13 03555 g007
Figure 8. The trajectory of the experimental platform.
Figure 8. The trajectory of the experimental platform.
Electronics 13 03555 g008
Figure 9. Attitude measurement error.
Figure 9. Attitude measurement error.
Electronics 13 03555 g009
Figure 10. RMSE of attitude measurement error.
Figure 10. RMSE of attitude measurement error.
Electronics 13 03555 g010
Table 1. Detailed information of data collected from the experiment.
Table 1. Detailed information of data collected from the experiment.
IndexDeviceInformationUnitSampling Frequency
1main IMUgyroscope°/h100 Hz
accelerometerm/s2100 Hz
2Sub-IMUgyroscope°/h100 Hz
accelerometerm/s2100 Hz
3PHINS-1gyroscope°/h100 Hz
accelerometerm/s2100 Hz
4PHINS-2gyroscope°/h100 Hz
accelerometerm/s2100 Hz
5GNSSlatitude°1 Hz
longitude°1 Hz
Table 2. Installation deviation between testing IMUs and PHINS.
Table 2. Installation deviation between testing IMUs and PHINS.
ItemInstallation Deviation1Installation Deviation2
Pitch (°)−0.041−0.029
Roll (°)0.4300.525
Heading (°)−0.2500.330
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhao, B.; Xia, X.; Wang, T.; Gao, W. Dynamic Attitude Inertial Measurement Method for Typical Regions of Deck Deformation. Electronics 2024, 13, 3555. https://doi.org/10.3390/electronics13173555

AMA Style

Zhao B, Xia X, Wang T, Gao W. Dynamic Attitude Inertial Measurement Method for Typical Regions of Deck Deformation. Electronics. 2024; 13(17):3555. https://doi.org/10.3390/electronics13173555

Chicago/Turabian Style

Zhao, Bo, Xiuwei Xia, Tianyu Wang, and Wei Gao. 2024. "Dynamic Attitude Inertial Measurement Method for Typical Regions of Deck Deformation" Electronics 13, no. 17: 3555. https://doi.org/10.3390/electronics13173555

APA Style

Zhao, B., Xia, X., Wang, T., & Gao, W. (2024). Dynamic Attitude Inertial Measurement Method for Typical Regions of Deck Deformation. Electronics, 13(17), 3555. https://doi.org/10.3390/electronics13173555

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop