Next Article in Journal
GOES-R Time Series for Early Detection of Wildfires with Deep GRU-Network
Previous Article in Journal
Water Level Change Monitoring Based on a New Denoising Algorithm Using Data from Landsat and ICESat-2: A Case Study of Miyun Reservoir in Beijing
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Dynamic Self-Tuning Maximum Correntropy Kalman Filter for Wireless Sensors Networks Positioning Systems

School of Automation, Beijing Institute of Technology, Beijing 100081, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(17), 4345; https://doi.org/10.3390/rs14174345
Submission received: 26 June 2022 / Revised: 23 August 2022 / Accepted: 29 August 2022 / Published: 1 September 2022

Abstract

:
To improve the accuracy of the maximum correntropy Kalman filter (MCKF) in wireless sensors networks (WSNs) positioning, a dynamic self-tuning maximum correntropy Kalman filter (DSTMCKF) is proposed, where innovation and the sensors information of the WSNs are used to adjust the noise covariance matrices, and the maximum correntropy criterion is the criterion for the filter’s optimality. By dynamically adjusting the noise covariance matrices, the DSTMCKF ensures that the correntropy distribution is accurate in the presence of non-Gaussian noise (NGN), thus improving its ability to handle the NGN. In simulation and real environment positioning experiments, the DSTMCKF is used to compare with the MCKF, variable kernel width–maximum correntropy Kalman filter (VKW-MCKF) and robust minimum error entropy Kalman filter (R-MEEKF). Among the four filters, the DSTMCKF has the highest accuracy, and the error of the DSTMCKF is reduced by 34.5, 42.9 and 40.0%, respectively, compared with the MCKF, VKW-MCKF and R-MEEKF in the real-world environment positioning experiment. The application of the DSTMCKF in WSNs positioning systems improves the stability of the control systems because of the rising positioning accuracy, which makes WSNs positioning systems more widely used in scenarios requiring high stability, such as automatic parking.

Graphical Abstract

1. Introduction

With the development of sensors technology, wireless sensors networks (WSNs) positioning systems are widely used in factories, docks, warehouses and other places and have broad development prospects [1,2,3]. In recent years, Gulo et al. [4] and Duan et al. [5] have applied WSNs positioning technology to occasions other than industrial applications, such as mobile phone positioning and vehicle positioning, further improving the practicability of WSNs positioning technology. Because of a low cost and high positioning accuracy, the range-based positioning method has become the most shared WSNs positioning method. Authors such as Xiong et al. [6], Na et al. [7] and Nguyen et al. [8] have all adopted range-based positioning methods to effectively improve the positioning accuracy.
However, the range-based positioning algorithm is seriously affected by non-Gaussian noise (NGN) [9,10], which results in its inaccuracy in practical application. Huang et al. [11], Fang et al. [12] and Navon et al. [13] all pointed out that WSNs positioning systems are affected by various NGNs, which seriously reduce the positioning accuracy.
To suppress the influence of NGN, various robust Kalman filter methods are applied to WSNs positioning systems [5,14]. The classical robust Kalman filter includes the H∞ filter [15,16] and the strong tracking Kalman filter [17,18]. However, these filters are not effective in dealing with the noise in WSNs positioning problems, and the positioning accuracy decreases seriously.
In recent years, the maximum correntropy Kalman filter (MCKF) was proposed by Chen et al. [19], according to the information theoretic learning criterion. Compared with the traditional robust Kalman filter, the MCKF has a higher accuracy and has the ability to deal with various NGNs [20]. Owing to its excellent performance, the MCKF is rapidly applied to various systems [21,22]. However, in WSNs positioning systems, most of the NGN comes from the electromagnetic radiation of electronic equipment and non-line-of-sight signals due to occlusion and refraction. These noises are usually much larger than the Gaussian white noise in the channel, and the size and duration are difficult to predict. With the MCKF, it is hard to deal with these noises, so the filtering accuracy of the MCKF will be reduced in the WSNs positioning systems [23,24].
To improve the performance of the MCKF in WSNs positioning systems, a dynamic self-tuning maximum correntropy Kalman filter (DSTMCKF) is proposed. In order to deal with the non-Gaussian noise, a dynamic self-tuning (DST) algorithm is proposed to combine with the MCKF. The DST algorithm uses innovation and the information from sensors to modify the noise covariance matrices, and the MCKF uses a maximum correntropy criterion (MCC) for the filter’s optimality.
By using the DST algorithm to adjust the noise covariance matrices in real time, the response speed and accuracy of the filter when encountering non-Gaussian noise are improved. Compared with the minimum mean squared error (MMSE) criterion, the maximum correntropy criterion has stronger robustness to ensure the convergence of the filter when encountering non-Gaussian noise. By combining the DST algorithm and the MCC, the DSTMCKF ensures that the correntropy distribution is accurate in the presence of NGN, thus improving its ability to handle NGN.
To verify the performance of the DSTMCKF, three filters, including the MCKF, variable kernel width–maximum correntropy Kalman filter (VKW-MCKF) and robust minimum error entropy Kalman filter (R-MEEKF), are used to compare with the DSTMCKF in the simulation experiments and real-world dynamic positioning experiments. The simulation experiments consider the different non-Gaussian noise in the motion state and observation, and the real-world dynamic positioning experiment is selected in the real office scene.
The rest of the paper is organized as follows: In Section 1, the principle of the MCKF is presented in detail, along with the limitations of the MCKF. In Section 2, a DSTMCKF is proposed, and the principle of the DSTMCKF is described. In Section 3, the simulation and real-world dynamic experiments results are presented, both comparing the DSTMCKF to other filters and exploring some of its qualities.

2. Maximum Correntropy Kalman Filter

2.1. Principle of MCKF

To suppress the influence of NGN on positioning accuracy, the MCKF is used to deal with NGN and has achieved promising results. Compared with the MMSE, the MCC has a stronger robustness and lower risk when facing NGN [25]. Therefore, compared with the classical Kalman filter, the MCKF is more suitable for dealing with NGN [19].
The MCC is an optimal estimation criterion, also known as the maximum information criterion. For two variables X and Y whose joint distribution function is F X Y ( x , y ) , the definition of correntropy V ( X , Y ) is
V ( X , Y ) = E κ ( X , Y ) = κ ( X , Y ) d F X Y ( x , y ) ,
where E is the expectation operator and κ ( X , Y ) represents the kernel function of correntropy.
Generally, κ ( X , Y ) is
κ ( X , Y ) = G σ ( e ) = exp ( e 2 2 σ 2 ) ,
where G σ is the Gaussian function, e is equal to X Y and σ is the bandwidth of the Gaussian kernel function.
When applying MCC to KF, a linear system is considered first. The state equation and observation equation of linear system are
x k = F k x k 1 + w k y k = H k x k + v k ,
where x k m is the system state at the k-th sampling, w k is the process noise, y k n is the observation at the k-th sampling and v k is the observation noise. State transition matrix and observation matrix are F k m × m and H k m × n respectively.
Then, the process noise covariance matrix and observation noise covariance matrices Q k and R k are
Q k = E w k w k T R k = E v k v k T .
In order to apply the MCC to Kalman filter, it is necessary to establish the correntropy model. The a priori estimate D k and a posteriori estimate W k x ^ k | k in the Kalman filter are used as correlation variables in the correntropy model. Then, construct matrices D k and W k are
B k = chol P k | k 1 0 0 R k ,
D k = B k 1 [ x ^ k | k 1 y k ] T ,
W k = B k 1 [ I H k ] T ,
where P k | k 1 is the a priori estimate of process noise covariance matrix. The prior estimation of the process noise covariance matrix P k | k 1 in this iteration can be calculated by using the state transition matrix F k and the posterior estimation of the process noise covariance matrix P k 1 obtained in the previous iteration. In addition, chol ( x ) is a Cholesky decomposition operator, which means Cholesky decomposition of x.
In order to calculate the correntropy of D k and W k x ^ k | k , E k needs to be calculated as
E k = D k W k x ^ k | k ,
where x ^ k | k is the a priori estimate of the k-th sampling.
Then, the correntropy V ( D k , W k x ^ k | k ) of D k and W k x ^ k | k is
V ( D k , W k x ^ k | k ) = 1 m + n i = 1 m + n G σ ( e ˜ k ( i ) ) ,
where e k ( i ) is the i-th element of vector E k .
According to MCC, when V ( D k , W k x ^ k | k ) is the maximum, x ^ k | k is the optimal value. To calculate the optimal value of x ^ k | k , fixed-point iteration method is used to calculate x ^ k | k . The steps of MCKF are summarized as follows:
(a)
Prediction:
x ^ k | k 1 = F k x ^ k 1 | k 1 P k | k 1 = F k P k 1 F k T + Q k 1 ,
where x k | k 1 is the prior estimate of the current iteration obtained by the state transition matrix F k and the state posterior estimation x k 1 | k 1 of the previous iteration.
(b)
Update required parameters of MCKF:
B p = chol ( P k | k 1 ) B r = chol ( R k ) .
(c)
Fixed-point iteration:
C ˜ x = diag ( G σ ( e ˜ k ( 1 ) ) , , G σ ( e ˜ k ( n ) ) ) C ˜ y = diag ( G σ ( e ˜ k ( n + 1 ) ) , , G σ ( e ˜ k ( n + m ) ) ) ,
R ˜ k = B r C ˜ y 1 B r T P ˜ k | k 1 = B p C ˜ x 1 B p T ,
K ˜ k = P ˜ k | k 1 H T k [ H k P ˜ k | k 1 H k T + R ˜ k ] 1 ,
x ^ k | k ( t ) = x ^ k | k 1 + K ˜ k [ y k H k x ^ k | k 1 ] ,
where K ˜ k is the Kalman filter gain, and t represents the number of iteration rounds of fixed-point iteration.
(d)
The judgement criterion for iteration stop is
x ^ k | k ( t ) x ^ k | k ( t 1 ) x ^ k | k ( t ) ε ,
where ε is the error threshold. It can be seen from document [19], to ensure the accuracy and numerical stability of the algorithm, the value of ε is 10 4 in this paper. When Equation (16) holds, the fixed-point iteration is stopped; otherwise, continue with step (c).
(e)
Update the filtering results and prepare for the next iteration:
x ^ k | k = x ^ k | k ( t ) ,
P k = [ I K ˜ k H k ] P k | k 1 [ I K ˜ k H k ] T .
Through the above steps, the state with maximum correntropy is estimated. MCKF replaces MMSE criteria with MCC, which has higher robustness when the noise size cannot be determined. However, MCKF has limitations in practical application, which result in underdeveloped filtering effect and even divergence of filtering results.

2.2. Limitations of MCKF

In practical application, when the noise type or target motion state is multiple, the filtering accuracy of the MCKF will decrease to a certain extent. In order to explain this phenomenon, a situation is considered when the MCKF is working properly [26].
Considering that the a priori estimate and observation measurement are only one dimension. Supposing the state prior estimate x k | k 1 is 40, the state observation value y k is 60, process noise covariance matrix P k | k 1 is 6, observation noise covariance matrix R k is 1 and the value of Gaussian kernel σ is 5. Where P k | k 1 and x k | k 1 are the prior estimators of the current iteration obtained by the state transition matrix F k and the posterior estimators of the previous iteration, and the values are not affected by the observation noise of the current iteration.
If the values of the estimated noise covariance matrices P k | k 1 and R k are not much different from the values of the noise covariance in reality, the distribution of correntropy is shown in Figure 1a. The correntropy has only one maximum value; hence, the fixed-point iteration ensures convergence to the point with maximum correntropy. In this case, MCKF has strong robustness and has a better ability to deal with NGN.
On the contrary, when the MCKF is interfered by noise of an unpredictable size, the actual observed noise covariance is inconsistent with R k . At this time, the observation y k will also change due to the influence of observation noise. If the noise further increases, the difference between state a priori estimation and state observation will further increase. Assuming that the state view measurement is 80 and other variables remain unchanged, the distribution of correntropy is shown in Figure 1b.
There are three extreme points in the solution space of the fixed-point iteration algorithm under these conditions, which do not satisfy the iteration convergence condition and provoke a large error in the calculation results. Thus, when the observation noise or the system motion state is multiple, the performance of the MCKF is underdeveloped.
One reason behind this phenomenon is that the MCKF is sensitive to the values of P k | k 1 and R k , but the values cannot be dynamically corrected. When the values of P k | k 1 and R k are inconsistent with the reality, there will be errors in the distribution of correntropy, resulting in the degradation of the MCKF performance.
Another reason is that the fixed-point iteration is difficult to converge to the maximum correntropy point when the observation noise or the system motion state is multiple. Under these conditions, there are multiple extreme points in the solution space, which do not meet the convergence condition of the fixed-point iteration. The fixed-point iteration method may converge to the extreme point far away from the maximum correntropy point, resulting in a large error.
To solve this problem, it is necessary to modify P k | k 1 and R k in each iteration [11,27]. The adaptive algorithms realize the dynamic correction of parameters to a certain extent, and yet when dealing with NGN in WSNs positioning systems, the adaptive algorithm has a slow response speed and insufficient robustness. Therefore, several algorithms have been proposed to improve the filtering accuracy of the MCKF in recent years.

2.3. Improved Algorithm of MCKF

With the aim of improving the performance of the MCKF, several algorithms have been proposed in recent years. There are two common ways to improve the accuracy: one is to adjust the Gaussian kernel bandwidth and the other is to optimize the maximum correntropy criterion.
Much of the research has examined how to improve the filtering precision of the MCKF by using a modified Gaussian kernel bandwidth, and the variable kernel width–maximum correntropy Kalman filter (VKW-MCKF) proposed by Huang et al. [24] has achieved positive results. The adaptive algorithm is used to adjust the Gaussian kernel bandwidth in the VKW-MCKF, so that the Gaussian kernel bandwidth can be adjusted adaptively with the size of the noise, which effectively improves the accuracy and robustness of the MCKF. However, the VKW-MCKF does not have the ability to adjust the noise covariance matrix, so there is still an inaccurate distribution of the correntropy when the noise does not match the estimation.
In recent years, a series of studies have expected to improve the filtering accuracy of the MCKF by optimizing the maximum correntropy criterion, which the robust minimum error entropy Kalman filter (R-MEEKF) proposed by Chen et al. [23] performed favorably. The R-MEEKF optimizes the maximum correntropy principle to the minimum entropy principle, which further improves the ability of the filter to process non-Gaussian noise. Nonetheless, the R-MEEKF algorithm requires a trade-off between the convergence speed and steady-state error, and it is often difficult to guarantee that both performances meet the requirements.
In summary, although several algorithms have improved the filtering accuracy of the MCKF in recent years, none of the above improved algorithms can solve the problem that the MCKF is difficult to handle the variable process and observation noise. The above improved algorithms are difficult to apply in positioning systems with complex and variable noise; hence, a dynamic self-tuning algorithm is proposed to adjust the values of the noise covariance matrix P k | k 1 and R k .

3. Dynamic Self-Tuning Maximum Correntropy Kalman Filter

3.1. Determination of Noise Type for WSNs Positioning

In order to dynamically adjust P k | k 1 or R k , it is necessary to judge whether NGN appears in the observation. A WSNs positioning quality index is proposed, which is calculated according to the distance information between the tag and anchors.
The quality index is obtained by trilateral positioning method, and the principle is shown in Figure 2. In Figure 2a, ( X 1 , Y 1 ) , ( X 2 , Y 2 ) and ( X 3 , Y 3 ) are the coordinates of the three anchors. d 1 , d 2 and d 3 are the distances to tag measured by the three anchors [28,29]. Then, the coordinate ( X T , Y T ) of tag is
X T Y T = 1 2 X B X A Y B Y A X C X A Y C X A 1 · d A 2 d B 2 + X B 2 X A 2 + Y B 2 Y A 2 d A 2 d C 2 + X C 2 X A 2 + Y C 2 Y A 2 .
As shown in Figure 2b, when there is an error in the measured distance, the three circles in the triangulation positioning method cannot intersect at the same point. At this time, the distance d ^ i from the positioning result ( X T , Y T ) to the Anchor i is not equal to the measured distance d i . Utilizing the difference between d ^ i and d i , a WSNs positioning quality evaluation index is proposed.
The positioning quality index L is
L = i = 1 N d i d ^ i N ,
where N is the number of anchor used for positioning.
The L obtained by the k-th sampling is denoted as L k . When the L k is large, it indicates that the difference between d ^ i and d i is large in this sampling. In addition, the greater the difference between d ^ i and d i , the less reliable the positioning result. At this time, NGN appears in the observation, and R k needs to be adjusted.
From another perspective, if ( X T , Y T ) is taken as the unbiased estimation of the positioning result, L k is the mean ranging error of each positioning result. To test whether the ranging error satisfies the hypothesis, the hypothesis testing method is used to determine whether L k is within the rejection domain. When L k is not in the rejection domain, it means that the ranging error conforms to the hypothesis; otherwise, it means that the ranging error does not conform to the hypothesis. The critical point of the rejection domain is the threshold L t h r , so when L k is greater than L t h r , it is considered that the ranging error does not conform to the hypothesis, and there is NGN in the observation.

3.2. Dynamic Self-Tuning Algorithm

After determining whether the NGN exists in the observation, the next step is to correct the noise covariance matrices. A dynamic self-tuning (DST) algorithm is presented, which uses innovation to adjust the value of P k | k 1 or R k . Record innovation z k as:
z k = y k H k x k | k 1 .
When the estimation of the model is accurate and the noise meets the Gaussian distribution, the expectation C k of innovation covariance is shown in Equation (22):
C k = Δ H k P k | k 1 H k T + R k .
In the actual filtering process, the covariance C ˜ k is calculated in Equation (23). In addition, according to [30], after weighing the robustness of the algorithm and the response speed to NGN, the forgetting factor λ is taken as 0.95.
C ˜ k = z k z k T , k = 1 λ C ˜ k 1 + z k z k T 1 + λ , k > 1 .
When the estimation of the model is accurate and the noise meets the Gaussian distribution, the value of C ˜ k is similar to C k . On the contrary, when the motion state changes or there is NGN in the observation, the difference between C ˜ k and C k will increase significantly. At this time, P k | k 1 or R k needs to be tuned to force C ˜ k equal to C k .
If the L k or L k 1 in Section 3.1 has an increased value than that of the threshold L thr , there is NGN in the observation, so it is necessary to adjust R k to R ^ k , forcing Equation (24) to hold.
C ˜ k = H k P k | k 1 H k T + R ^ k .
Note that the intermediate variable N k is:
N k = C ˜ k H k P k | k 1 H k T .
Equation (24) holds when R ^ k is equal to N k . In order to ensure that R ^ k is positive definite, the diagonal element of R k is set by the expansion coefficient β i | k , and then the value of non-diagonal element is corrected by using the correlation coefficient, as shown below:
β i | k = max [ 1 , n i i r i i ] ,
R ^ k = β 1 | k r 11 β 1 | k β 2 | k r 12 β 1 | k β n | k r 1 n β 1 | k β 2 | k r 21 β 2 | k r 22 β 2 | k β n | k r 2 n β 1 | k β n | k r n 1 β 2 | k β n | k r n 2 β n | k r n n .
where n i j and r i j are the elements of row i and column j of matrices N k and R k , respectively.
When the L k or L k 1 in Section 3.1 is smaller than the threshold L thr , the possibility of NGN appearing in the observation is little or nothing. If the difference between x k | k 1 and y k is still large under these condition, the motion state has changed, and P k | k 1 needs to be corrected to P ^ k | k 1 .
Theoretically, P ^ k | k 1 needs to make Equation (28) valid. However, because H k is not necessarily reversible, the value of P ^ k | k 1 cannot be directly calculated by Equation (28). Therefore, the expansion coefficient α k is introduced to adjust P ^ k | k 1 , as shown in Equations (29) and (30):
C ˜ k = H k P ^ k | k 1 H k T + R k ,
α k = max [ 1 , tr ( C ˜ k R k ) tr ( H k P k | k 1 H k T ) ] ,
P ^ k | k 1 = α k P k | k 1 .
Based on the above steps, a DST algorithm is proposed, as shown in Figure 3. Because there is a possibility of misjudgment when only L k is used, both L k and L k 1 are required to be smaller than L thr to ensure the robustness of the algorithm.
The DST algorithm uses sensors information as a priori information to effectively adjust the noise covariance matrices. When the noise is non-stationary, the DST algorithm can correct the noise covariance matrix to be closer to the current noise in real time. This means that the DST algorithm can still ensure the accurate estimation of the noise covariance matrix in the presence of NGN. Combined with the DST algorithm, the robustness and accuracy of the MCKF will be improved significantly.

3.3. Dynamic Self-Tuning Maximum Correntropy Kalman Filter

The DST algorithm is combined with the MCKF, and then the dynamic self-tuning maximum correntropy Kalman filter (DSTMCKF) is proposed. The DST algorithm is used for correcting the values of P k | k 1 or R k firstly. Then, the modified noise covariance matrices are used for the Kalman filter, and the filter result x ^ k | k ( 0 ) is used as the starting point of the fixed-point iteration.
When the estimates of P k | k 1 and R k are close to the true value, the value of x ^ k | k ( 0 ) is very close to the maximum correntropy point. Taking x ^ k | k ( 0 ) as the starting point of the fixed-point iteration ensures that there is only one maximum value in the solution space of the fixed-point iteration. Finally, the state quantity with the maximum correntropy is calculated by using the fixed-point iterative algorithm, which is the filtering result of DSTMCKF. The above process is expressed as Equations (31) to (48):
(a)
Prediction:
x ^ k | k 1 = F k x ^ k 1 | k 1 P k | k 1 = F k P k 1 F k T + Q k 1 .
(b)
Adjust noise covariance matrix by DST algorithm:
z k = y k H k x k | k 1 ,
C ˜ k = z k z k T , k = 1 λ C ˜ k 1 + z k z k T 1 + λ , k > 1 .
IF L k L thr | L k 1 L thr :
N k = C ˜ k H k P k | k 1 H k T ,
β i | k = max [ 1 , n i i r i i ] ,
P ^ k | k 1 = P k | k 1 ,
R ^ k = β 1 | k r 11 β 1 | k β 2 | k r 12 β 1 | k β n | k r 1 n β 1 | k β 2 | k r 21 β 2 | k r 22 β 2 | k β n | k r 2 n β 1 | k β n | k r n 1 β 2 | k β n | k r n 2 β n | k r n n .
If L k < L thr   &   L k 1 < L thr :
α k = max [ 1 , tr ( C ˜ k R k ) tr ( H k P k | k 1 H k T ) ] ,
P ^ k | k 1 = α k P k | k 1 ,
R ^ k = R k .
(c)
Calculation of fixed-point iteration starting point by KF:
K k = P ^ k | k 1 H k T [ H k P ^ k | k 1 H k T + R ^ k ] 1 ,
x ^ k | k ( 0 ) = x ^ k | k 1 + K k [ y k H k x ^ k | k 1 ] .
(d)
Parameters required to calculate MCKF:
B p = chol ( P ^ k | k 1 ) ,
B r = chol ( R ^ k ) ,
B k = B p 0 0 B r ,
D k = B k 1 [ x ^ k | k 1 y k ] T ,
W k = B k 1 [ I H k ] T ,
E k = D k W k x ^ k | k .
(e)
Fixed-point iteration:
C ˜ x = diag ( G σ ( e ˜ k ( 1 ) ) , , G σ ( e ˜ k ( n ) ) ) C ˜ y = diag ( G σ ( e ˜ k ( n + 1 ) ) , , G σ ( e ˜ k ( n + m ) ) ) ,
R ˜ k = B r C ˜ y 1 B r T P ˜ k | k 1 = B p C ˜ x 1 B p T ,
K ˜ k = P ˜ k | k 1 H k T [ H k P ˜ k | k 1 H k T + R ˜ k ] 1 ,
x ^ k | k ( t ) = x ^ k | k 1 + K ˜ k [ y k H k x ^ k | k 1 ] .
(f)
When Equation (53) holds, the fixed-point iteration is stopped; otherwise, continue with step (e):
x ^ k | k ( t ) x ^ k | k ( t 1 ) x ^ k | k ( t ) ε .
According to the above steps, the filtering result of DSTMCKF is obtained. The process of DSTMCKF is also summarized as shown in Figure 4, and DST algorithm is in Figure 3.
When combined with the DST algorithm, the ability of the MCKF to process NGN is theoretically improved. The noise covariance matrices, which can be dynamically corrected with the noise size, enable the DSTMCKF to have higher accuracy and robustness than the MCKF. Compared with the VKW-MCKF, the DSTMCKF still works well when the noise covariance matrices estimation is inaccurate. The VKW-MCKF, which can only adjust the size of the Gaussian kernel bandwidth, is theoretically less efficient than the DSTMCKF when estimating the noise covariance matrices inaccurately. In addition, compared with the R-MEEKF, the DSTMCKF does not need to balance accuracy with the response speed and has stronger robustness in theory.

4. Simulation and Application in Office Indoor Environment

4.1. Simulation of Indoor Positioning Scene and Analysis of Filtering Results

In order to verify the performance of the DSTMCKF in WSNs positioning systems, the shared situations in the WSNs positioning systems are considered for simulation firstly. Two motion states and three kinds of noise between the anchors and tag are considered.
The simulation sampling frequency is set to 100 Hz. The two motion states are the uniform linear motion process and the acceleration–deceleration process. During the uniform linear motion process, the speed of the tag in the x-direction is 1 m/s, the speed in the y-direction is 0 m/s and the process noise on the speed in each sampling is ω k N ( 0 , 0.01 ) . During the acceleration–deceleration process, the acceleration of the tag in the x-direction is zero, the acceleration in the y-direction is 5   m / s 2 in the first 0.5 s and 5   m / s 2 in the last 0.5 s. The process noise on the speed in each sampling is ω k N ( 0 , 0.025 ) .
The simulation environment is shown in Figure 5. The simulation uses three anchors with coordinates of (0,0), (20,0) and (10,17.32) to form an equilateral triangle with a side length of 20 m. There are three kinds of shared noises between the anchors and tag. Note that the three kinds of noise are S 1 , S 2 and S 3 , and the distribution of the noise is shown in Table 1, where A i represents the channel between the i-th anchor and the tag.
Under the above simulation conditions, the MCKF, VKW-MCKF, R-MEEKF and DSTMCKF are used for the comparative experiments. Assuming that there is only noise S 1 , the significance level of the hypothesis testing is taken as 0.05, and then the threshold value L thr is calculated as 0.13. The state transition matrix F k is obtained according to the kinematic characteristics of the simulation object, as shown in Equation (54). Because the WSNs can only measure the position of the simulation object, the observation matrix H k is shown in Equation (55).
F k = 1 T 0 0 0 1 0 0 0 0 1 T 0 0 0 1 ,
H k = 1 0 0 0 0 0 1 0 .
Four cases are discussed in the simulation experiment, as shown below:
Case (1): In Case (1), there is an acceleration–deceleration phase and noise S 1 in Table 1. The Gaussian kernel σ of the MCKF, R-MEEKF and DSTMCKF is 1.2, and because the VKW-MCKF will diverge when the Gaussian kernel σ is 1.2, the Gaussian kernel σ of the VKW-MCKF in the figure is 5.0. The simulation sampling results and experimental results are shown in Figure 6, where Figure 6e is the cumulative distribution function (CDF) of each filter.
The experimental results indicate that when the motion states are multiple, the DSTMCKF has the highest accuracy, while the filtering accuracy of the MCKF has decreased to a certain extent. In order to discuss the influence of bandwidth σ on the filtering results, the value of σ is modified, and the results are shown in Table 2.
When the value of σ is small, the dynamic tracking ability of the MCKF is underdeveloped. In addition, the error is large or even divergent when the motion state is multiple. With the increase of σ , the MCKF gradually converges to KF, and the dynamic performance is improved.
The response speed of the VKW-MCKF and R-MEEKF is improved to some extent compared with the MCKF, so the positioning accuracy is higher in the presence of multiple motion states. However, because neither of the two algorithms have the ability to adjust the noise covariance values adaptively, the accuracy decreases when the observed noise does not match the estimation.
For the DSTMCKF, a small σ also leads to a decrease in the dynamic performance, and an increase of σ will invalidate the MCC. Thus, the DSTMCKF will perform best only when σ is moderate. Nonetheless, the dynamic tuning of parameters reduces the sensitivity of the DSTMCKF to σ , and the RMSE of the DSTMCKF changes little when σ changes.
Case (2): In Case (2), the noise S 1 in Table 1 occurs in the whole process and the noise S 2 in Table 1 occurs during iterations 500 to 600 s. The Gaussian kernel σ of the MCKF, R-MEEKF and DSTMCKF is 1.2. In addition, because the VKW-MCKF will diverge when the Gaussian kernel σ is 1.2, the Gaussian kernel σ of the VKW-MCKF in the figure is 5.0. The simulation sampling results and experimental results are shown in Figure 7.
The DSTMCKF performs better than the MCKF when σ is 1.2. However, because the motion state is single at this time, σ takes a smaller value to improve the performance of the MCKF. The experimental results are shown in Table 3. When σ is small enough, the accuracy of the MCKF is improved to a certain extent. Nonetheless, to ensure the dynamic performance of the filter, it is difficult for σ to obtain such a small value.
Due to the variable bandwidth of the Gaussian kernel, the VKW-MCKF has a better ability to deal with NGN in observation, and the positioning accuracy is almost the same as that of the DST-MCKF. The steady-state error of the R-MEEKF is lower than that of the MCKF, but the corresponding convergence speed is slower, and the anti-jamming ability is underdeveloped; therefore, the positioning accuracy of the R-MEEKF is not better than the MCKF. The DSTMCKF has the ability to adjust the noise covariance matrices during operation; hence, it has a better ability to handle the NGN in the observations and has the highest accuracy.
Case (3): In Case (3), the noise S 1 in Table 1 occurs in the whole process and the noise S 3 in Table 1 occurs during iterations 700 to 800. The Gaussian kernel σ of the MCKF, R-MEEKF and DSTMCKF is 1.2. However, because the VKW-MCKF will diverge when the Gaussian kernel σ is 1.2, the Gaussian kernel σ of the VKW-MCKF in the figure is 5.0. The simulation sampling results and experimental results are shown in Figure 8.
It can be seen from the figure that in case (3), the accuracy of the MCKF is also lower than that of the DSTMCKF. And not similar to case (2), the performance of the MCKF is not improved by reducing the value of σ , as shown in Table 4.
In Table 4, even if σ is small enough, the MCKF will be severely affected by the noise. Conversely, the DSTMCKF still ensures a high accuracy no matter how the value of σ changes. This means that the DSTMCKF will still converge when the motion state changes and has higher robustness and numerical stability than the MCKF. This means that the DSTMCKF has a stronger performance in processing NGN with offset.
The VKW-MCKF has a high accuracy when the Gauss kernel bandwidth σ is small, yet the numerical stability of the VKW-MCKF is underdeveloped. In practice, the Gauss kernel bandwidth σ is usually larger. The R-MEEKF is relatively insensitive to the change of the Gaussian kernel σ and has strong numerical stability, but the accuracy of the R-MEEKF is the lowest.
Case (4): In Case (4), there is an acceleration–deceleration phase in the motion state. The noise S 1 in Table 1 occurs in the whole process, the noise S 2 in Table 1 occurs during iterations 500 to 600 and the noise S 3 in Table 1 occurs during iterations 700 to 800. The Gaussian kernel σ of the MCKF, R-MEEKF and DSTMCKF is 1.2. In addition, because the VKW-MCKF will diverge when the Gaussian kernel σ is 1.2, the Gaussian kernel σ of the VKW-MCKF in the figure is 5.0. The simulation sampling results and experimental results are shown in Figure 9.
When the motion state and noise type are multiple, it is difficult for σ to obtain an appropriate value to make the MCKF have better accuracy. In order to discuss the influence of bandwidth σ on the filtering results, the value of σ is modified, and the results are shown in Table 5.
When the value of σ is small, the dynamic tracking ability of the MCKF is underdeveloped. On the contrary, when the value of σ increases, the ability of the MCKF to deal with NGN decreases. Therefore, when σ is equal to 2.0, the MCKF has the best accuracy. Nevertheless, the DSTMCKF has the ability to self-tune parameters, and it has the highest accuracy no matter how σ changes. The DSTMCKF also performs better in a simulation than the VKW-MCKF and R-MEEKF.
In addition, in each iteration, the averaged simulation time for the DSTMCK, MCKF and Kalman filter is 0.144, 0.096 and 0.025 ms, respectively. For the control systems, in order to meet the control requirements, it is required that the frequency of the uploading positioning results can reach more than 100 Hz [31]. Although the DSTMCKF is more computationally burdensome, it is still within the acceptable range.
In summary, the DSTMCKF outperforms the MCKF and its improvement method in all four cases and effectively handles situations where there is no single motion state or noise type. Meanwhile, the influence of σ on the DSTMCKF is limited, which indicates that the DSTMCKF has better numerical stability and robustness than the MCKF. The simulation results show that the DSTMCKF is more suitable for WSNs positioning systems than the MCKF and its improvement method.

4.2. Application in Office Indoor Environment

In order to prove that the DSTMCKF is applied to the actual WSNs positioning systems, experiments are carried out with an ultra-wideband (UWB) positioning network in the real environment, and the experimental environment in the real world is shown in Figure 10. The indoor positioning environment selected in this experiment is a common office environment, with an area of about 10 × 7.5 m2. In the experiment, three UWB anchors and one UWB tag are selected to form a UWB positioning network. The positions of the three anchors are shown in Figure 11a, respectively, (0, 0), (8.9, 0) and (4.45, 7.4). The tag is placed on the robot and travels around the indoor aisle. There are several mark points in the room, and the robot counts the time when passing through the mark points. The recorded time and the coordinates of the mark points are used as features for Gaussian smoothing to obtain an ideal path. The ideal path is shown by the orange line in the figure, and the sampled data are shown by the blue line in the figure.
The positioning result of the UWB is affected by various NGNs. Most of these NGNs come from the electromagnetic radiation of electronic equipment and non-line-of-sight signals due to occlusion and refraction. To deal with these NGNs, the VKW-MCKF, R-MEEKF, MCKF and DSTMCKF are used in the experiment, and the results are shown in Figure 11b–d.
It is seen from the figure that when NGN does not exist, the performance of the three filters is almost the same. When NGN exists, the accuracy and robustness of the DSTMCKF are better than the MCKF, VKW-MCKF and R-MEEKF. The RMSE and maximum error (ME) of the three filtering results are shown in the Table 6. The experiment results indicate that the DSTMCKF also has better performance than the MCKF and its improvement method, and the error is reduced by 34.5% compared to the MCKF.
To sum up, the DSTMCKF solves the problem of the underdeveloped performance of the MCKF in the noise type or motion state to a certain extent. It has a higher accuracy and robustness and significantly improves the positioning accuracy of the systems in the actual WSNs positioning systems.

5. Discussion

The range-based positioning method is the most shared positioning method in WSNs positioning systems, and the positioning results are affected by NGN [32,33,34,35,36]. To improve the accuracy of WSNs positioning systems, a DSTMCKF is proposed. Four simulation scenarios are considered to test the performance of the DSTMCKF, and the DSTMCKF has a significantly improved performance in various scenarios compared to other common filtering algorithms. To further prove the reliability of the DSTMCKF, the algorithm is used in real-world indoor scenes. The results show that the DSTMCKF has a higher accuracy than other filtering algorithms, can adapt to a complex noise environment in real scenes and effectively improves the accuracy of a WSNs positioning system.
The main idea of this paper is to improve the accuracy of a WSNs positioning system with a DSTMCKF. We mainly discuss the effects of the non-Gaussian process and observation noise on the positioning results and use the DSTMCKF to mitigate these effects. Because the NGN in the real-world scene is not as extreme as the NGN in the simulation scene, the advantages of the DSTMCKF are not as obvious as those in the simulation experiments. If the positioning experiments can be conducted in more complex indoor scenes, the advantages of the DSTMCKF will be further highlighted. In future work, we plan to investigate the improvements and extensions of the method for determining the type of noise. This helps reduce the risk of filter divergence due to misidentified noise types, improves filter accuracy in complex environments and enhances the robustness of the filter.

6. Conclusions

A dynamic self-tuning maximum correntropy Kalman filter is proposed to improve the performance of an MCKF. The DSTMCKF realizes the dynamic self-tuning of noise covariance matrices on the basis of innovation and a priori information from sensors. Applying the adjusted noise covariance matrices to the MCC effectively improves the performance of the MCKF when there is NGN or a motion state that does not conform to the estimate. In this way, the DSTMCKF solves the problem that the MCKF performs underdeveloped when there is NGN or a motion state that does not conform to the estimate. The simulation and experiments in the real-world environment indicate that the DSTMCKF also has a better performance than the MCKF, and the error of the DSTMCKF is reduced by 34.5, 42.9 and 40.0%, respectively, compared with the MCKF, VKW-MCKF and R-MEEKF in the real-world environment.
The experimental results prove that the DSTMCKF improved the response speed to NGN and the accuracy by using the DST algorithm to adjust the noise covariance matrices in real time. It also proves that the DSTMCKF improved the ability to process non-Gaussian noise by combining the DST algorithm with the MCKF.
The proposed DSTMCKF improves the filter precision in a WSNs positioning system, which improves the positioning accuracy of the WSNs positioning system. A higher positioning accuracy brings more scenarios for WSNs positioning systems. This means that WSNs positioning systems using a DSTMCKF are qualified for use in scenarios with higher accuracy requirements, such as automatic parking.

Author Contributions

This paper is a collaborative work by all the authors. T.L. proposed the main idea, designed the experiments, performed the experiments, analyzed the data and wrote the manuscript. K.H. and Y.D. added some ideas, gave suggestions and revised the rough draft. X.W. and S.S. assisted with certain experiments, and all authors proofread the paper. All authors have read and agreed to the published version of the manuscript.

Funding

This research was sponsored by the Beijing Municipal Natural Science Foundation Committee (L191020) and the China National Railway Corporation (P2021T002).

Data Availability Statement

Not applicable.

Acknowledgments

We would like to acknowledge the School of Automation, Beijing Institute of Technology, for the support in our research.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
WSNsWireless Sensors Networks
MCKFMaximum Correntropy Kalman Filter
DSTMCKFDynamic Self-Tuning Maximum Correntropy Kalman Filter
VKW-MCKFVariable Kernel Width–Maximum Correntropy Kalman Filter
R-MEEKFRobust Minimum Error Entropy Kalman Filter
DSTDynamic Self Tuning
NGNNon-Gaussian noise
MMSEMinimum Mean Squared Error
CDFCumulative Distribution Function
MCCMaximum Correntropy Criterion
RMSERoot Mean Squared Error
MEMaximum Error

References

  1. Yuan, C.; Shuangshuang, Z.; Guizhi, Y.; Tao, W. Research on high precision tracking method of guided transport vehicle based on autonomous combination positioning. Chin. J. Electron. 2020, 29, 779–785. [Google Scholar]
  2. Jin, Y.; Zhou, L.; Zhang, L.; Hu, Z.; Han, J. A Novel Range-Free Node Localization Method for Wireless Sensor Networks. IEEE Wirel. Commun. Lett. 2022, 11, 688–692. [Google Scholar] [CrossRef]
  3. Xu, Y.; Tian, G.; Chen, X. Enhancing INS/UWB Integrated Position Estimation Using Federated EFIR Filtering. IEEE Access 2018, 6, 64461–64469. [Google Scholar] [CrossRef]
  4. Gulo, E.; Sohn, G.; Ahmad, A. Indoor positioning using wlan fingerprint matching and path assessment with retroactive adjustment on mobile devices. Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 2018, 42, 321–327. [Google Scholar] [CrossRef]
  5. Duan, J.; Shi, H.; Liu, D.; Yu, H. Square Root Cubature Kalman Filter-Kalman Filter Algorithm for Intelligent Vehicle Position Estimate. Procedia Eng. 2016, 137, 267–276. [Google Scholar] [CrossRef]
  6. Xiong, H.; Peng, M.; Gong, S.; Du, Z. A Novel Hybrid RSS and TOA Positioning Algorithm for Multi-Objective Cooperative Wireless Sensor Networks. IEEE Sensors J. 2018, 18, 9343–9351. [Google Scholar] [CrossRef]
  7. Na, H.J.; Yoo, S.J. PSO-Based Dynamic UAV Positioning Algorithm for Sensing Information Acquisition in Wireless Sensor Networks. IEEE Access 2019, 7, 77499–77513. [Google Scholar] [CrossRef]
  8. Nguyen, C.L.; Georgiou, O.; Gradoni, G.; Di Renzo, M. Wireless Fingerprinting Localization in Smart Environments Using Reconfigurable Intelligent Surfaces. IEEE Access 2021, 9, 135526–135541. [Google Scholar] [CrossRef]
  9. Kang, S.; Kim, T.; Chung, W. Target Localization with Unknown Transmit Power and Path-Loss Exponent Using a Kalman Filter in WSNs. Sensors 2020, 20, 6582. [Google Scholar] [CrossRef]
  10. Xie, J.; Li, C.; Zhang, J.; Lv, A. Research on train positioning based on extended kalman particle filter. Int. J. Electr. Eng. 2019, 26, 111–117. [Google Scholar]
  11. Huang, Y.Y.; Jing, Y.W.; ShI, Y.B. Non-parametric mobile node localization for IOT by variational Bayesian approximations adaptive Kalman filter. Cogn. Syst. Res. 2018, 52, 27–35. [Google Scholar] [CrossRef]
  12. Fang, X.; Nan, L.; Jiang, Z.; Chen, L. Noise-aware fingerprint localization algorithm for wireless sensor network based on adaptive fingerprint Kalman filter. Comput. Netw. 2017, 124, 97–107. [Google Scholar] [CrossRef]
  13. Navon, E.; Bobrovsky, B. An efficient outlier rejection technique for Kalman filters. Signal Process. 2021, 188, 108164. [Google Scholar] [CrossRef]
  14. Liu, J.; Guo, G. Pseudolinear kalman filters for target tracking using hybrid measurements. Signal Process. 2021, 188, 108206. [Google Scholar] [CrossRef]
  15. Chen, B.S.; Yang, K.C.; Lee, M.Y. Robust H-infinity NLOS-Tolerant Localization Filter and NLOS-Tolerant Remote Reference Tracking Control of Mobile Robot in Wireless Sensor Networks. IEEE Access 2021, 9, 164801–164819. [Google Scholar] [CrossRef]
  16. Wang, M.; Dong, X.; Qin, C.; Liu, J. Adaptive H-infinity Kalman filter based random drift modeling and compensation method for ring laser gyroscope. Measurement 2021, 167, 108170. [Google Scholar] [CrossRef]
  17. Lu, D.H.; Wang, J.Q.; Xiong, K.; Hou, B.W.; He, Z.M. Strong tracking Kalman filter for non-Gaussian observation. Kongzhi Lilun Yu Yingyong/Control Theory Appl. 2019, 36, 1997–2004. [Google Scholar]
  18. Xie, S.B.; Zhang, Y.; Wen, K.R.; Zhang, S.; Liu, Z.M.; Qi, N.M. Motion estimation for non-cooperative target based on strong tracking cubature Kalman filter. Jilin Daxue Xuebao (Gongxueban)/J. Jilin Univ. 2021, 51, 1482–1489. [Google Scholar]
  19. Chen, B.; Liu, X.; Zhao, H.; Principe, J.C. Maximum correntropy Kalman filter. Automatica 2017, 76, 70–77. [Google Scholar] [CrossRef]
  20. Deng, Z.; Shi, L.; Yin, L.; Xia, Y.; Huo, B. UKF Based on Maximum Correntropy Criterion in the Presence of Both Intermittent Observations and Non-Gaussian Noise. IEEE Sens. J. 2020, 20, 7766–7773. [Google Scholar] [CrossRef]
  21. Fan, X.; Wang, G.; Han, J.; Wang, Y. Interacting Multiple Model Based on Maximum Correntropy Kalman Filter. IEEE Trans. Circuits Syst. II Express Briefs 2021, 68, 3017–3021. [Google Scholar] [CrossRef]
  22. Wang, G.; Zhang, Y.; Wang, X. Maximum Correntropy Rauch–Tung–Striebel Smoother for Nonlinear and Non-Gaussian Systems. IEEE Trans. Autom. Control 2021, 66, 1270–1277. [Google Scholar] [CrossRef]
  23. Wang, G.; Chen, B.; Yang, X.; Peng, B.; Feng, Z. Numerically stable minimum error entropy Kalman filter. Signal Process. 2021, 181, 107914. [Google Scholar] [CrossRef]
  24. Huang, F.; Zhang, J.; Zhang, S. Adaptive Filtering Under a Variable Kernel Width Maximum Correntropy Criterion. IEEE Trans. Circuits Syst. II Express Briefs 2017, 64, 1247–1251. [Google Scholar] [CrossRef]
  25. Xiong, W.; Schindelhauer, C.; So, H.C.; Wang, Z. Maximum Correntropy Criterion for Robust TOA-Based Localization in NLOS Environments. Circuits Syst. Signal Process. 2021, 40, 6325–6339. [Google Scholar] [CrossRef]
  26. Zhou, N.; Lau, L.; Bai, R.; Moore, T. A genetic optimization resampling based particle filtering algorithm for indoor target tracking. Remote Sens. 2021, 13, 132. [Google Scholar] [CrossRef]
  27. Scardua, L.A.; da Cruz, J.J. Complete offline tuning of the unscented Kalman filter. Automatica 2017, 80, 54–61. [Google Scholar] [CrossRef]
  28. Li, C.; Zhen, J.; Chang, K.; Xu, A.; Zhu, H.; Wu, J. An indoor positioning and tracking algorithm based on angle-of-arrival using a dual-channel array antenna. Remote Sens. 2021, 13, 4301. [Google Scholar] [CrossRef]
  29. Xu, S.; Rice, M.; Rice, F. Optimal TOA-Sensor Placement for Two Target Localization Simultaneously Using Shared Sensors. IEEE Commun. Lett. 2021, 25, 2584–2588. [Google Scholar] [CrossRef]
  30. Liu, X.; Guan, J.; Jiang, R.; Gao, X.; Chen, B.; Ge, S.S. Robust strong tracking unscented Kalman filter for non-linear systems with unknown inputs. IET Signal Process. 2022, 16, 351–365. [Google Scholar] [CrossRef]
  31. Xhafa, A.; Del Peral-Rosado, J.A.; Lopez-Salcedo, J.A.; Seco-Granados, G. Evaluation of 5g positioning performance based on utdoa, aoa and base-station selective exclusion. Sensors 2022, 22, 101. [Google Scholar] [CrossRef] [PubMed]
  32. Cui, Z.; Liu, T.; Tian, S.; Xu, R.; Cheng, J. Non-Line-of-Sight Identification for UWB Positioning Using Capsule Networks. IEEE Commun. Lett. 2020, 24, 2187–2190. [Google Scholar] [CrossRef]
  33. Zhao, Y.; Li, Z.; Hao, B.; Shi, J. Sensor Selection for TDOA-Based Localization in Wireless Sensor Networks with Non-Line-of-Sight Condition. IEEE Trans. Veh. Technol. 2019, 68, 9935–9950. [Google Scholar] [CrossRef]
  34. Yin, Z.; Jiang, X.; Yang, Z.; Zhao, N.; Chen, Y. WUB-IP: A high-precision UWB positioning scheme for indoor multiuser applications. IEEE Syst. J. 2019, 13, 279–288. [Google Scholar] [CrossRef] [Green Version]
  35. Kanhere, O.; Rappaport, T.S. Position location for futuristic cellular communications: 5G and beyond. IEEE Commun. Mag. 2021, 59, 70–75. [Google Scholar] [CrossRef]
  36. Geng, C.; Yuan, X.; Huang, H. Exploiting Channel Correlations for NLOS ToA Localization with Multivariate Gaussian Mixture Models. IEEE Wirel. Commun. Lett. 2020, 9, 70–73. [Google Scholar] [CrossRef]
Figure 1. Distribution of correntropy in different scenes. (a) Distribution of correntropy when parameters are accurate; (b) distribution of correntropy when NGN increases.
Figure 1. Distribution of correntropy in different scenes. (a) Distribution of correntropy when parameters are accurate; (b) distribution of correntropy when NGN increases.
Remotesensing 14 04345 g001
Figure 2. Different quality positioning result. (a) High-quality positioning result; (b) low-quality positioning result.
Figure 2. Different quality positioning result. (a) High-quality positioning result; (b) low-quality positioning result.
Remotesensing 14 04345 g002
Figure 3. DST Algorithm.
Figure 3. DST Algorithm.
Remotesensing 14 04345 g003
Figure 4. Flowchart of DSTMCKF algorithm.
Figure 4. Flowchart of DSTMCKF algorithm.
Remotesensing 14 04345 g004
Figure 5. Simulated environment.
Figure 5. Simulated environment.
Remotesensing 14 04345 g005
Figure 6. Simulation Results in Case (1). (a) Positioning result of X-label when σ = 1.2 ; (b) error of X-label when σ = 1.2 ; (c) positioning result of Y-label when σ = 1.2 ; (d) error of Y-label when σ = 1.2 ; (e) CDF of each filter when σ = 1.2 .
Figure 6. Simulation Results in Case (1). (a) Positioning result of X-label when σ = 1.2 ; (b) error of X-label when σ = 1.2 ; (c) positioning result of Y-label when σ = 1.2 ; (d) error of Y-label when σ = 1.2 ; (e) CDF of each filter when σ = 1.2 .
Remotesensing 14 04345 g006
Figure 7. Simulation Results in Case (2). (a) Positioning result of X-label when σ = 1.2 ; (b) error of X-label when σ = 1.2 ; (c) positioning result of Y-label when σ = 1.2 ; (d) error of Y-label when σ = 1.2 ; (e) CDF of each filter when σ = 1.2 .
Figure 7. Simulation Results in Case (2). (a) Positioning result of X-label when σ = 1.2 ; (b) error of X-label when σ = 1.2 ; (c) positioning result of Y-label when σ = 1.2 ; (d) error of Y-label when σ = 1.2 ; (e) CDF of each filter when σ = 1.2 .
Remotesensing 14 04345 g007
Figure 8. Simulation Results in Case (3). (a) Positioning result of X-label when σ = 1.2 ; (b) error of X-label when σ = 1.2 ; (c) positioning result of Y-label when σ = 1.2 ; (d) error of Y-label when σ = 1.2 ; (e) CDF of each filter when σ = 1.2 .
Figure 8. Simulation Results in Case (3). (a) Positioning result of X-label when σ = 1.2 ; (b) error of X-label when σ = 1.2 ; (c) positioning result of Y-label when σ = 1.2 ; (d) error of Y-label when σ = 1.2 ; (e) CDF of each filter when σ = 1.2 .
Remotesensing 14 04345 g008
Figure 9. Simulation Results in Case (4). (a) Positioning result of X-label when σ = 1.2 ; (b) error of X-label when σ = 1.2 ; (c) positioning result of Y-label when σ = 1.2 ; (d) error of Y-label when σ = 1.2 ; (e) CDF of each filter when σ = 1.2 .
Figure 9. Simulation Results in Case (4). (a) Positioning result of X-label when σ = 1.2 ; (b) error of X-label when σ = 1.2 ; (c) positioning result of Y-label when σ = 1.2 ; (d) error of Y-label when σ = 1.2 ; (e) CDF of each filter when σ = 1.2 .
Remotesensing 14 04345 g009
Figure 10. Experimental environment in real world.
Figure 10. Experimental environment in real world.
Remotesensing 14 04345 g010
Figure 11. Results in real world. (a) Sampling results in real world; (b) filtering results in real world; (c) filtering error in real world; (d) CDF of each filter.
Figure 11. Results in real world. (a) Sampling results in real world; (b) filtering results in real world; (c) filtering error in real world; (d) CDF of each filter.
Remotesensing 14 04345 g011
Table 1. Three kinds of noise between the tag and the anchors.
Table 1. Three kinds of noise between the tag and the anchors.
NoiseDistribution (m)Channel with Noise
S 1 S k N ( 0 , 0.02 ) A 1 , A 2 , A 3
S 2 l n S k N ( 0 , 0.05 ) A 1 , A 2 , A 3
S 3 S k U ( 0 , 0.3 ) + | N ( 0 , 0.2 ) | A 3
Table 2. RMSE (m) with different σ (Case 1).
Table 2. RMSE (m) with different σ (Case 1).
σ DSTMCKFMCKFVKW-MCKFR-MEEKF
0.40.0210.055Divergence0.029
0.80.0190.028Divergence0.023
1.20.0180.027Divergence0.023
2.00.0160.026Divergence0.023
5.00.0150.0260.0330.023
10.00.0150.0260.0320.024
20.00.0150.0260.0320.025
Table 3. RMSE (m) with different σ (Case 2).
Table 3. RMSE (m) with different σ (Case 2).
σ DSTMCKFMCKFVKW-MCKFR-MEEKF
0.40.0070.019Divergence0.032
0.80.0070.041Divergence0.035
1.20.0070.043Divergence0.042
2.00.0080.044Divergence0.046
5.00.0090.0450.0130.047
10.00.0090.0450.0190.047
20.00.0100.0450.0300.047
Table 4. RMSE (m) with different σ (Case 3).
Table 4. RMSE (m) with different σ (Case 3).
σ DSTMCKFMCKFVKW-MCKFR-MEEKF
0.40.0080.073Divergence0.072
0.80.0230.073Divergence0.072
1.20.0270.074Divergence0.073
2.00.0320.074Divergence0.074
5.00.0340.0740.0630.074
10.00.0340.0740.0690.074
20.00.0340.0740.0720.074
Table 5. RMSE (m) with different σ (Case 4).
Table 5. RMSE (m) with different σ (Case 4).
σ DSTMCKFMCKFVKW-MCKFR-MEEKF
0.40.0420.104Divergence0.094
0.80.0440.102Divergence0.100
1.20.0460.102Divergence0.105
2.00.0480.101Divergence0.103
5.00.0480.1030.0790.104
10.00.0480.1030.0860.105
20.00.0510.1040.1010.105
Table 6. RMSE and maximum error of three filters.
Table 6. RMSE and maximum error of three filters.
RMSE (m)ME (m)
DSTMCKF0.0360.17
MCKF0.0550.18
VKW-MCKF0.0630.33
R-MEEKF0.0600.35
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liao, T.; Hirota, K.; Wu, X.; Shao, S.; Dai, Y. A Dynamic Self-Tuning Maximum Correntropy Kalman Filter for Wireless Sensors Networks Positioning Systems. Remote Sens. 2022, 14, 4345. https://doi.org/10.3390/rs14174345

AMA Style

Liao T, Hirota K, Wu X, Shao S, Dai Y. A Dynamic Self-Tuning Maximum Correntropy Kalman Filter for Wireless Sensors Networks Positioning Systems. Remote Sensing. 2022; 14(17):4345. https://doi.org/10.3390/rs14174345

Chicago/Turabian Style

Liao, Tianrui, Kaoru Hirota, Xiangdong Wu, Shuai Shao, and Yaping Dai. 2022. "A Dynamic Self-Tuning Maximum Correntropy Kalman Filter for Wireless Sensors Networks Positioning Systems" Remote Sensing 14, no. 17: 4345. https://doi.org/10.3390/rs14174345

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