Next Article in Journal
V2G Strategies to Flatten the Daily Load Curve in Seoul, South Korea
Next Article in Special Issue
Application and Research on Improved Adaptive Monte Carlo Localization Algorithm for Automatic Guided Vehicle Fusion with QR Code Navigation
Previous Article in Journal
Investigating Polarization-Sensitive Transmission and Reflection Metasurfaces for Advanced Wavefront Manipulation
Previous Article in Special Issue
CERRT: A Mobile Robot Path Planning Algorithm Based on RRT in Complex Environments
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Research on Multi-Sensor Data Fusion Positioning Method of Unmanned Ships Based on Threshold- and Hierarchical-Capacity Particle Filter

1
School of Mechanical Engineering, Jiangsu University of Science and Technology, Zhenjiang 212100, China
2
Department of Intelligent Equipment Research, Zhangjiagang JUST Industrial Technology Research Institute, Zhangjiagang 215600, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2023, 13(18), 10390; https://doi.org/10.3390/app131810390
Submission received: 20 August 2023 / Revised: 10 September 2023 / Accepted: 14 September 2023 / Published: 17 September 2023
(This article belongs to the Special Issue Advances in Robot Path Planning, Volume II)

Abstract

:
To improve the positioning accuracy of unmanned ships, a multi-sensor system including ZigBee, a Global Positioning System (GPS), and BeiDou Navigation Satellite System (BDS) is constructed, and an adaptive multi-sensor data fusion positioning method based on the threshold and hierarchical capacity particle filter (TCPF) is designed. First, the ZigBee-GPS/BDS multi-sensor measurement data is preprocessed to achieve a consistent space–time reference and transformed into the same coordinate system by projection. Then, the fault data is weighted and corrected through the consistency inspection of ZigBee-GPS/BDS multi-sensor positioning data, and the corresponding confidence factor is given according to the confidence distance of the positioning data; furthermore, the confidence factor is associated with stratified sampling. After that, the multi-sensor positioning data is filtered and denoised using a basic particle filter. Finally, a TCPF data fusion algorithm is designed, and the navigation positioning data of the unmanned ship is fused and filtered to obtain its positioning information. Numerical tests show that compared with other filtering algorithms, the mean square root error and standard deviation of the proposed TCPF algorithm decrease by an average of 25.0% and 28.0%, respectively, which verifies its high filtering accuracy and its advantages in suppressing particle degradation and avoiding sample scarcity. The experimental tests show that compared with other fusion algorithms, the proposed TCPF algorithm can not only realize the precise positioning during unmanned ship navigation, but also in the positioning and fault tolerance test, the average positioning error, root-mean-square error, and standard deviation of the former decrease by 36.0%, 38.0%, and 37.0%, respectively, and the corresponding performance indicators of the latter decrease by an average of 20.0%, 19.5%, and 17.5%, which verifies that it has the advantages of high data reliability and good filtering fault tolerance, and helps to improve the positioning accuracy of unmanned ships.

1. Introduction

The water quality of rivers directly affects people’s lives. In the early days, traditional water quality testing mainly relied on manual testing and reporting methods. Currently, water quality testing still relies primarily on manual testing, but data is uploaded through the Internet. Both traditional and current water quality testing face the following problems: (1) high labor input and cost; (2) high manual collection intensity and low efficiency; (3) data easily influenced by subjective factors; and (4) untimely information transmission. In recent years, unmanned ships have been widely used in the field of water quality monitoring due to their strong maneuverability and good controllability [1]. Based on this background, this paper conducts research on modern water quality testing based on unmanned ships. Unmanned ships carry various water quality sensors such as pH value, conductivity, turbidity, dissolved oxygen, etc., and can autonomously conduct regular water quality tests throughout the river area and transmit real-time data to a monitoring center via the Internet. This method has the following advantages: it (1) reduces labor costs; (2) improves detection efficiency; (3) unaffected by subjective factors, improves accuracy; and (4) ensures timely data availability. From the above analysis, it can be seen that unmanned ships play a key role in modern water quality monitoring, with their autonomous navigation positioning accuracy directly affecting safety and efficiency.
How to achieve accurate positioning of unmanned ships based on multiple combined sensors has been the focus of research [2]. Shen et al. [3] used GNSS/GPS combination sensors to obtain coordinates of unmanned ship positions while reducing noise interference caused by long-distance nodes, but they suffer from signal loss due to external environmental interference leading to reduced positioning accuracy. Wu et al. [4] realized the positioning of the unmanned ship by solving its position information collected by the GPS/INS combined sensors, but the INS has the problem of mechanical fault jumping. As time goes by, the GPS/INS combined sensor is likely to cause a large positioning deviation of the unmanned ship. Deng et al. [5] obtained the position information of the unmanned ship by fusing the GPS/IMU combined sensor data, but in the fusion process, if the sensor data is lost at a certain moment, it will cause a large fusion deviation and affect the positioning accuracy. Data fusion based on two kinds of sensors is the current mainstream technology for unmanned ship positioning, but there is a deficiency in that the positioning accuracy and reliability are affected due to the weakening or loss of a certain sensor signal. In recent years, how to improve positioning accuracy by further fusing data from three or more types of sensors has become a research hotspot. Wang et al. [6] achieved the reliable positioning of trains by using the combined positioning method of PPP-GPS/IMU. However, the built multi-sensor positioning system is prone to environmental interference, which reduces the reliability of data samples. Tang et al. [7] realized the elimination of abnormal signals by setting outlier identification in the process of IMU/ODM/UWB multi-sensor data fusion, which improved the positioning accuracy of multi-sensors, but reduced the fault-tolerant performance of the algorithm because the fault data was not weighted. Sofia et al. [8] divided the observation information of different sensors into blocks according to the size of the scale to obtain a multi-scale system model, which effectively solved the problem of measurement delay, but the efficiency and accuracy of algorithm fusion were reduced because no data confidence assignment was performed during data processing. As can be seen from the above, adding sensors helps to improve the reliability and accuracy of positioning. However, there is still a lack of research on how to further enhance the accuracy and efficiency of algorithm fusion through data confidence interval testing, as well as how to achieve effective positioning under low- or weak-signal conditions through data correction. Therefore, a ZigBee-GPS/BDS multi-sensor positioning system is built in this paper, which realizes the accurate and reliable positioning of unmanned ships by complementing and weighted-correcting signals from three types of sensors, as well as assigning confidence values. In multi-sensor data fusion, filtering algorithms are crucial, and particle filtering has been widely used due to its ability to handle various probability distribution models, but it has also become the focus of research due to its heavy dependence on initial state estimation and particle degradation. Ha et al. [9] expanded the search space of particle states by designing crossover operators and mutation operators in the particle filter calculation module, which improved the filtering accuracy of multi-sensor data but also reduced the computational efficiency. Alam et al. [10] processed the prefetched value in the weight storage of the particle-filter algorithm in parallel with the value in the random function generator, which reduced the time required for resampling and improved the fusion efficiency of sensor data, but it did not suppress the degradation of particle samples very well. Wu et al. [11] used the unscented particle-filter algorithm with constrained residuals to fuse the two sets of sensor positioning data, which effectively overcomes external environmental disturbances. However, the filtering accuracy is prone to divergence after multiple iterations. A large number of studies have shown that although the existing multi-sensor data fusion algorithms have improved the positioning accuracy of moving targets to a certain extent, their fault-tolerant performance is not good. In addition, they did not carry out pre-judgment processing on the collected data, and the problem of low filtering accuracy has not been solved when the data collected by the multi-sensor combined positioning system is fused and filtered. Therefore, a multi-sensor data fusion positioning method based on adaptive threshold and hierarchical capacity particle filter is designed in this paper. Firstly, the basic particle filter is introduced to denoise multi-sensor data. Then, the latest multi-sensor positioning data is incorporated into the proposal distribution using unscented transformation. After that, the Gaussian mixture model is constructed and the adaptive threshold is set, and the data confidence factor is associated with hierarchical sampling to further improve the fusion filtering accuracy of the algorithm. Finally, the validity of the method is verified by the numerical test of the model and the experimental test of unmanned ship navigation.

2. Unmanned Ship Positioning System and Its Multi-Sensor Data Fusion Framework

2.1. Construction of Unmanned Ship Positioning System

In order to realize the information complementarity among the multi-sensor measurement data, reduce the interference of the surrounding environment on the working state of the positioning system, improve the reliability of the multi-sensor data samples, and then improve the positioning accuracy of an unmanned ship for water quality detection, a ZigBee-GPS/BDS multi-sensor combined positioning system, as shown in Figure 1, was constructed. The entire positioning system includes a ZigBee multi-node positioning system, GPS positioning system, BDS positioning system, PC communication system, etc.
A GPS/BDS positioning system is mainly composed of space, ground monitoring, and user receiver. First, the ground monitoring part monitors and controls the operation of each satellite through the main control station and monitoring station, traverses the navigation message, and maintains the system time. Then, the GPS and BD satellites in space continuously send their own ephemeris and time information to the GPS/BDS dual module installed in the unmanned ship. Finally, the GPS/BDS dual module calculates the latitude and longitude coordinates of the unmanned ship in real time by analyzing the signal messages sent by the respective satellites [12].
A ZigBee multi-node positioning system is mainly composed of blind nodes, reference nodes, and wireless gateway nodes. First, the communication environment of the full grid network is constructed according to the reference nodes with known position coordinates. Then, a blind node installed on the unmanned ship is placed in a communication network composed of reference nodes [13]. After that, the wireless gateway node is connected to the PC through the serial port. On the one hand, it receives the configuration data of each reference node and mobile node from the monitoring software and forwards them to the corresponding nodes. On the other hand, it receives the valid data fed back by each node and transmits them to the monitoring software. Finally, the RSSI technology is used to determine the distance or direction from the unknown node to the beacon node, and the maximum likelihood estimation method is used to calculate the position of the blind node on the unmanned ship during navigation in the river.
Figure 2 shows the self-developed unmanned ship for water quality testing. The unmanned ship is equipped with a detector with multiple water quality detection sensors connected to the launch and recovery device, as well as a positioning system including Zigbee, GPS, and BDS modules. The controller of the unmanned ship is an industrial computer with a Linux system, which is responsible for receiving and processing the data of each sensor, issuing control instructions, and transmitting signals with the PC.
The lidar of Silan A3 is responsible for collecting the environmental information around the unmanned ship. The ATK1218-GPS-BDS dual-positioning module is selected in the GPS/BDS positioning system to obtain the latitude and longitude information of the unmanned ship. The CC2530 wireless module is selected in the ZigBee positioning system to obtain the plane coordinate information of the unmanned ship. The Hikvision DS-2 camera is selected in the vision system to be responsible for the collection of environmental information of the unmanned-ship channel.

2.2. Framework of Multi-Sensor Data Fusion

In order to improve the reliability, fault tolerance, and accuracy of multi-sensor data fusion for the unmanned ship for water quality detection, a new multi-sensor-data-fusion positioning framework based on particle filtering as, shown in Figure 3, was designed. From the figure, it can be seen that the longitude and latitude information of the unmanned ship was first unified in time and space, and projected into the local coordinate system constructed by the ZigBee multi-node module through coordinate transformation, so as to realize mutual complementarity among sensor data. Then, through the consistency inspection of the multi-sensor data, the fault data are weighted and corrected to improve the fault-tolerant performance of the data fusion algorithm. At the same time, the corresponding confidence factor was given according to the confidence distance of the positioning data, and the confidence factor was associated with hierarchical sampling to improve the fusion efficiency and accuracy of the data fusion algorithm. After that, the basic particle filter was used to de-noise the multi-sensor positioning data to improve the reliability of the data samples. Finally, an improved particle-filter algorithm was designed to fuse and filter the multi-sensor positioning data to achieve accurate positioning of the unmanned ship.

3. Pre-Processing of Multi-Sensor Data from the Unmanned Ship

3.1. System Modeling

3.1.1. Measurement Equation of GPS/BDS

The measurement equation of GPS/BDS consists of nonlinear pseudorange measurement equations [14]. The measurement equations of the BDS pseudorange pB and the GPS pseudorange pG are expressed as follows:
p i B ( k ) = r ( k ) r i B ( k ) + c δ t B + ε i ( k ) + v i ( k ) = ( X j ( k ) X i B ( k ) ) 2 + ( Y j ( k ) Y i B ( k ) ) 2 + ( Z j ( k ) Z i B ( k ) ) 2 + δ t B + ε i ( k ) + v i ( k )
p j G ( k ) = r ( k ) r j G ( k ) + c δ t G + ε j ( k ) + v j ( k ) = ( X j ( k ) X j G ( k ) ) 2 + ( Y j ( k ) Y j G ( k ) ) 2 + ( Z j ( k ) Z j G ( k ) ) 2 + δ t G + ε i ( k ) + v i ( k )
where Xj(k), Yj(k) and Zj(k) are the position coordinates of the receiver at time k; X i B (k), Y i B (k) and Z i B (k) are the position coordinates of the ith BDS satellite at time k; X j G (k), Y j G (k) and Z j G (k) are the position coordinates of the jth GPS satellite at time k; r(k) is the position vector of the receiver; r i B (k) and r j G (k) are the position vectors of the ith BDS and the jth GPS satellite at time k, respectively; p i B (k) and p j G (k) are the pseudoranges of the ith BDS and the jth GPS satellite at time k, respectively; δtB, δtG is the error of the BDS- and GPS-receiving clock, respectively; εi(k) is the non-white noise error of channel i at time k; ν i (k) is the measurement noise of channel i at time k, i = 1, 2, …, nl, j = 1, 2, …, n2, k = 1, 2, …, T. n1 and n2 are the number of BDS and GPS satellites received, respectively. T is the maximum duration.
When the signals of n satellites (n = n1 + n2 ≥ 5) are received, the observation variables of the system are
z B , G ( k ) = p 1 B ( k ) , p 2 B ( k ) , , p n 1 B ( k ) ,   p 1 G ( k ) , p 2 G ( k ) , , p n 2 G ( k ) T
Assuming that the number of GPS and BDS satellites observed at time k are equal (namely n1 = n2), the observation equation of the GPS/BDS satellite positioning system is established as
z B , G k = p B ( k ) p G ( k ) + R ( k )
where p B = { p 1 B , p 2 B , , p n 1 B } is the BDS pseudorange, p G = { p 1 G , p 2 G , , p n 2 G } is the GPS pseudorange, and R is the covariance matrix of observation noise.

3.1.2. Measurement Equation of ZigBee

According to the principle of ZigBee positioning networking in this paper, and taking the sampling time Tk of basic particle filter (BPF) as the time reference, the measurement equation of RSSI signal collected by the ZigBee positioning system was defined [15]. The specific steps were as follows:
Step 1 Establish the signal strength power model P ˜ i R of the ith reference node according to the path loss principle.
P ˜ i R = P 0 10 ϑ lg ( ( x z x i R ) 2 + ( y z y i R ) 2 / d 0 )
where ∀ i ∈ [1, n], ( x i R , y i R ) is the coordinate of the ith reference node; P0 is the received power of the reference node at d0 from the unmanned ship; ϑ is the path loss coefficient, and (xz, yz) is the blind node coordinate.
Step 2 Establish the covariance matrix R R S S I of the observation equation.
R R S S I = diag ( σ d B 1 R 2 , σ d B 2 R 2 , , σ d B n R 2 )
where ∀ i ∈ [1, n], σ d B i R 2 is the initial variance of the signal strength of the ith ZigBee reference node.
Step 3 Define the observation equation for RSSI measurements collected by ZigBee.
z R S S I ( k ) = [ P ˜ R ( k ) ,   R R S S I ( k ) ] T + v ( k )
where P ˜ R = { P ˜ 1 R , P ˜ 2 R , , P ˜ n R } is the signal strength power model; v is the covariance matrix of the observed noise.

3.1.3. Multi-Sensor Data Definition for the Unmanned Ship

In order to facilitate the description of the subsequent algorithm, the following data definition was carried out in this paper: The data collected by the sensor or the known data was defined as the measurement data, which involves the reference node coordinates (xR, yR) in the ZigBee system and the geodetic coordinates (B, L) collected by the GPS/BD system. The data obtained by calculation and conversion was defined as positioning data, which involves the coordinates of blind nodes (xz, yz) in ZigBee system, and the plane positioning coordinates (xG, yG) and (xB, yB) in GPS/BDS system. The GPS positioning data set was defined as Hg = { h 1 g , h 2 g , …, h n g }, the BDS positioning data set as Hb = { h 1 b , h 2 b , …, h n b }, and the ZigBee blind node positioning data set as Hz = { h 1 z , h 2 z , …, h n z }, where h i z = ( x i z , y i z ), h i g = ( x i G , y i G ), h i b = ( x i B , y i B ) are the ith positioning data in the ZigBee, GPS, and BDS data sets, respectively, i = 1, 2, …, n.

3.2. Unification of GPS/BDS Space-Time Reference

In order to obtain the trajectory information of the unmanned ship collected by the GPS/BDS dual-module sensor, it is necessary to unify the time of the GPS and BDS satellite systems. Due to the jump error between GPS time (GPST) and UTC time (UTCT) [16], as time goes by, the deviation between the two gradually increases, and at present the deviation has reached 16 s, that is
GPST UTCT + 16
BDS time (BDST) takes the national standard unit system s as the basic unit for time accumulation. Because there is no jump second in BDST, it is continuous time. The start time of BDST is converted into the number of weeks of GPST and the count of seconds in a week is 1356 weeks—14.000 s. Therefore, the following relationship exists between BDST and GPST
BDST week   number = GPST week   number + 1356
BDST Seconds   of   the   week   count = GPST Seconds   of   the   week   count + 14
Each country has established different spatial benchmarks according to its own satellite navigation system, so it is necessary to unify the spatial benchmarks when performing combined positioning of different satellite systems. WGS-84 was used for GPS, and the 2000 National Geodetic Coordinate System (NGCS2000) was used for the BDS. The origin, scale, orientation, and the scale evolved from the orientation defined by the two coordinate systems were the same, and both are closely related to the Earth Reference System (ERS). After multiple optimizations, the difference between the two is only at the centimeter level [17]. For most non-precision positioning applications such as navigation, the WGS-84 and NGCS2000 coordinate systems do not need to be converted, so they are considered as unified coordinate systems in this paper.

3.3. Network Positioning of Zigbee System

In order to meet the multi-node effective communication of the unmanned-ship Zigbee positioning system, n reference nodes were set in the positioning network coordinate system in this paper, so that the blind node on the unmanned ship is in a network composed of n reference nodes with a known signal strength and position coordinates [18]. The specific steps are as follows:
Step 1 Define the RSSI theoretical value of the signal strength indication of the received power, and calculate the distance between the blind node and the reference node.
R S S I = ( 10 ϑ lg d + A )
From Equation (11), we can obtain
d = 10 A + R S S I 10 ϑ
where ϑ is the path loss coefficient, d is the distance between the blind node and the reference node, RSSI is the signal strength value, and A is the strength of the initial signal at the transmitter.
Step 2 Establish the following distance equations among the blind node and the reference nodes in the plane:
( x 1 R x z ) 2 + ( y 1 R y z ) 2 ( x 2 R x z ) 2 + ( y 2 R y z ) 2 ( x n R x z ) 2 + ( y n R y z ) 2 = d 1 2 d 2 2   d n 2
where ∀ i ∈ [1, n], di is the distance between the blind node and the ith reference node.
Step 3 Based on the reference node’s distance equation set established by Equation (13), the linear transformation equation 2 A x z y z = B of the blind node coordinates of the unmanned ship was constructed, where A and B are expressed as follows:
A = x 1 R x n R y 1 R y n R x 2 R x n R   y 2 R y n R     x n 1 R x n R   y n 1 R y n R   B = d n 2 d 1 2 ( x n R ) 2 + ( x 1 R ) 2 ( y n R ) 2 + ( y 1 R ) 2 d n 2 d 2 2 ( x n R ) 2 + ( x 2 R ) 2 ( y n R ) 2 + ( y 2 R ) 2 d n 2 d n 1 2 ( x n R ) 2 + ( x n 1 R ) 2 ( y n R ) 2 + ( y n 1 R ) 2
Step 4 The coordinates (xz, yz) of the blind node in the ZigBee system network are obtained by the least-squares method.
x z y z = 1 2 ( A T A ) 1 A T B

3.4. Coordinate Transformation

In order to perform data fusion filtering on GPS/BDS positioning data and ZigBee blind node location information, it is necessary to transform the WGS-84/NGCS2000 coordinate system [19] and unify it to the local plane coordinates of the unmanned-ship navigation test with the ZigBee multi-node networking coordinate system. For this reason, the existing geodetic coordinates (B, L) were used in this paper to transform them into plane coordinates through Gauss–Kruger projection, and translate and rotate them to make them unified with the local plane coordinates of the unmanned ship navigation. The specific steps are as follows:
Step 1 convert WGS-84/NGCS2000 geodetic coordinates (B, L) to WGS-84/NGCS2000 Gaussian plane coordinates (xs, ys):
x s = X 0 B + 1 2 N p t m 0 2 + 1 24 ( 5 t 2 + 9 η 2 + 4 η 4 ) N p t m 0 4 + 1 720 ( 61 58 t 2 + t 4 ) N p t m 0 6
y s = N p m 0 + 1 6 ( 1 t 2 + η 2 ) N p m 0 3 + 1 120 ( 5 18 t 2 + t 4 + 14 η 2 58 η 2 t 2 ) N p m 0 5
where m 0 = Δ L γ cos ( B ) , t = tan ( B ) , η 2 = ( e ) 2 cos 2 B and Δ L are the difference between the longitude (L) of the desired point and the longitude (L0) of the central meridian, namely ΔL = LL0. γ is the positioning coefficient error, B is the latitude of the desired point, and its unit is radians, and e is the second eccentricity of the ellipsoid.
Np is the curvature radius of the primary vertical circle passing through this point, which can be described as
N p = a 1 e 2 sin 2 ϕ
where a is the long radius of the ellipsoid, e is the first eccentricity of the ellipsoid, and φ is the geodetic dimension.
X 0 B is the arc length from the equator to the central meridian cut by a parallel circle passing through this point and can be described as
X 0 B = C 0 cos B ( C 1 sin B + C 2 sin 3 B + C 3 sin 5 B + C 4 sin 7 B )
where C0, C1, C2, C3 and C4 are ellipsoid parameters.
Step 2 transform the Gaussian coordinates of the unmanned ship navigation collected by GPS/BDS into the local coordinate system of the unmanned ship test through translation and rotation.
x G ,   B y G ,   B = Δ x Δ y + ( 1 + m ) cos θ sin θ sin θ   cos θ x s y s
where (xG, B, yG, B) are the GPS/BDS plane positioning coordinates of the unmanned ship, Δx and Δy are the translation amount of the coordinate system, m is the scale factor, and θ is the rotation angle of the two coordinate systems.

4. Confidence Determination of Multi-Sensor Data

In order to improve the fusion efficiency and filtering accuracy of multi-sensor data of the unmanned ship, the confidence distance judgment and credibility assignment of ZigBee-GPS/BDS multi-sensor data was carried out in this paper. Firstly, the credibility of the positioning data was determined by performing a confidence interval test on the sensor positioning data. Then, the credibility was divided into the corresponding confidence intervals according to its size, and the corresponding confidence factors were assigned. Finally, the confidence factor was associated with the TCPF algorithm to improve the filtering accuracy of the algorithm.

4.1. Judgement of Confidence Distance

When the location information of the unmanned ship is positioned based on the ZigBee-GPS/BDS multi-sensor positioning system, its sensor characteristic function can be described by the Gaussian distribution curve p (xi) [20], and the positioning data xi follows the normal distribution N(μ, σ), namely
p ( x i ) = 1 2 π σ e ( x i μ ) 2 2 θ 2
where xi is the ith positioning data of the sensor, μ is the true value of the positioning data, θ is the measurement accuracy of the sensor information, and σ is the measurement error of the sensor information, i = 1, 2, …, n.
Let X = {Hz, Hg, Hb}, x i ZGB = { h i z , h i g , h i b }, ∀ i ∈ [1, n], x i ZGB X: the confidence distance of the sensor positioning data is determined using the defined p(xi). The specific steps are as follows:
Step 1 calculate the confidence distance among the positioning data of the ZigBee-GPS/BDS unmanned-ship positioning system at different times.
d i ( k ) = p r x i ZGB ( k ) x i ZGB ( k + 1 ) τ min ( τ i ZGB ( k ) , τ i ZGB ( k + 1 ) )
where pr{·} is a probability function, x i ZGB ( k ) and x i ZGB ( k + 1 ) are the observation values of the ith positioning data at time k and k + 1, respectively, τ i ZGB ( k ) and τ i ZGB ( k + 1 ) are the measurement variances of the ith positioning data at time k and k + 1, and τ is the mean of the measurement variance, k = 1, 2, …, T.
Step 2 determine the credibility of the sensor information at time k according to the confidence level ε of the preset sensor support.
High   data   credibility ,   i f   d i ( k ) 1 ε Low   data   credibility ,   i f   d i ( k ) > 1 ε
Step 3 Divide the positioning data into corresponding confidence intervals according to the confidence distance of the positioning data [21].
d i ( k ) > 1 ε ,   x i ZGB ( k ) [ ε τ ι , ε τ ] d i ( k ) 1 ε ,   x i ZGB ( k ) [ ε τ , ε τ + ι ]
where ι is the variable coefficient of the sampling sample in the probability interval.
In this paper, the credibility between sensor information at different times was determined by calculating the ZigBee-GPS/BDS confidence distance, and the determined data was divided into corresponding confidence intervals according to the probabilistic measurement value between the positioning data, which caused the credible data to be close to the region with higher confidence. Giving priority to the confidence factors of the sensor positioning data is helpful in improving the efficiency of data fusion. The set confidence level ε corresponded to the mean value τ of the sensor measurement variance, which was used to represent the change in the confidence degree of the multi-sensor positioning data at different times. The classification of the corresponding confidence interval also changed when the mean value τ of the measurement variance was different, which improved the assignment accuracy of the data fusion algorithm to the sensor positioning data.

4.2. Credibility Assignment

After the confidence distance was determined for the ZigBee-GPS/BDS positioning data, the multi-sensor positioning data was divided into corresponding confidence intervals. In order to improve the filtering accuracy of the data fusion algorithm, the confidence factor was given according to the confidence distance of the positioning data, and the confidence factor was associated with the hierarchical sampling of the proposed TCPF algorithm. The higher the confidence of the positioning data, the higher the corresponding fusion bias, which helps to improve the filtering accuracy of the TCPF algorithm for the location information of the unmanned ship. The specific steps are as follows:
Step 1 define βi as the comprehensive support degree of the ith positioning data, which is composed of several confidence weight coefficients [22].
β i = y 1 ω 1 + y 2 ω 2 + + y n ω n
ω i = 1 n 1 2 a + 1 n a i = 1 n d i
i = 1 n β i = 1
where ω i is the confidence weight coefficient of the ith positioning data, n is the numericized amount of information collected by multiple sensors, and a is the adjustment parameter and can be expressed as a = (n − 1)/2. y1, y2, …, yn are a set of non-negative numbers.
Step 2 according to the comprehensive support degree of the defined positioning data, the norm equation of the dynamic support factor β i o (k) for multi-sensor information and the probability model pr( x i ZGB ( k ) ) for the positioning data is constructed.
β i o ( k ) = p r ( x i ZGB ( k ) ) F i = 1 n p r ( x i ZGB ( k ) ) F
where ||.||F is the Frobenius norm, k = 1, 2, …, T, i = 1, 2, …, n.
Step 3 Calculate the multi-sensor measurement error w i ZGB according to the dynamic support factor β i o (k):
w i ZGB ( k ) = x i ZGB ( k ) A β i o ( k )
where A is the state transition matrix.
Step 4 The credibility of the positioning data of the unmanned ship positioning system is assigned using w i ZGB :
ς i ZGB ( k ) = ( w i ZGB ( k ) ) 1 ω i i = 1 n w i ZGB ( k )
where ς i ZGB ( k ) is the confidence factor of the ZigBee-GPS/BDS multi-sensor positioning data.

5. Inspection and Weighted Compensation of Multi-Sensor Data

In the data set collected by the ZigBee-GPS/BDS multi-sensor system of the unmanned ship, some data were valid, but some data may have caused measurement deviation due to environmental or noise interference. In order to improve the fault-tolerant performance of the data fusion algorithm, the consistency inspection of the positioning data set collected by the multi-sensor positioning system was carried out in this paper, and the inconsistency fault data was corrected by weighting.

5.1. Consistency Inspection of Sensor Data

For the multi-sensor positioning system composed of ZigBee-GPS/BDS, the ith positioning data x i ZGB can be expressed as
x i ZGB = x i t + ξ i
where x i t is the true value, and ξi is the measurement noise, i = 1, 2, …, n.
The obtained multi-sensor positioning data were arithmetically averaged and inspected for consistency [23]. The specific steps were as follows:
Step 1 the arithmetic mean of the positioning data, x ¯ , k = 1, 2, …, T, was calculated:
x ¯ i = ( k = 1 T x i ZGB ( k ) ) / T
Step 2 according to the arithmetic mean value of the positioning data, the consistency inspection of the multi-sensor positioning data was carried out.
x ZGB ( k ) x ¯ i τ
where xZGB(k) = { x 1 ZGB (k), x 2 ZGB (k), …, x n ZGB (k)} is the multi-sensor positioning data, and τ is the system requirement error.
Through the set system error requirements, the consistency inspection of the multi-sensor positioning data was carried out. If the difference between the sensor positioning data and the arithmetic mean value was less than the error required by the system, the positioning data was determined to be consistent, that is, credible data, and could be denoised using basic particle filtering. On the contrary, the variance of the sampled data needed weighted correction to meet the sampling requirements of the basic particle filter for data samples.

5.2. Weighted Correction for Variance

When the ZigBee-GPS/BDS multi-sensor positioning system was used to locate the navigation path of the unmanned ship at different positions in the same space, the measurement noise of each sensor followed Gaussian distribution [24]. The variance of the ith positioning data can be expressed as
σ i 2 = E ( x i ZGB x i t ) 2
In the real navigation situation of the multi-sensor positioning system of the unmanned ship, since the true positioning value x i t in Equation (34) is unknown, the arithmetic mean value x ¯ i of the positioning data in Equation (32) was taken as the unbiased estimate of the true value x i T , namely
σ i 2 = D ( x i ZGB x ¯ i )
In order to improve the credibility of data samples, the multi-sensor positioning system in this paper measured the position of the unmanned ship m times, and integrated the m times of positioning data into a data set. The jth positioning value of the ith sensor information was xij, and the sensor positioning data x i ZGB was replaced in Equation (35) with xij to obtain the information variance of the data set, which can be expressed as
σ ¯ i 2 = 1 m j = 1 m ( x i j x ¯ i ) 2
where j = 1, 2, …, m.
According to the variance of the data collected by the sensor, the fusion weight κi of the fault data was defined as follows:
κ i = ( σ ¯ i 2 ) 1 i = 1 n 1 σ ¯ i 2
i = 1 n κ i = 1
According to the obtained fusion weight κi, the multi-sensor fault data was weighted and fused.
x ˜ ( k ) = x ^ o ( k ) κ
where x ˜ ( k ) = { x ˜ 1 ( k ) , x ˜ 2 ( k ) ,…, x ˜ n ( k ) } is the multi-sensor positioning data after variance weighting, x o ( k ) = { x 1 o ( k ) , x 2 o ( k ) ,…, x n o ( k ) } is the multi-sensor fault data, and κ = { κ 1 , κ 2 ,…, κ n } is the fusion weight, i = 1, 2, …, n.
By weighting the variance of the fault data, a certain amount of information correction was realized, and the fault-tolerant performance of the data fusion algorithm was improved. Even if the ZigBee signal had measurement deviation due to environmental interference, or the GPS/BDS signal was weakened or lost due to environmental occlusion, data samples with high credibility based on the proposed method could still be obtained.

6. Denoising Processing of Positioning Data Based on Particle Filter

In order to improve the reliability of the information collected by the ZigBee-GPS/BDS multi-sensor positioning system and reduce external noise interference, the multi-sensor data after the consistency inspection and variance weighting were filtered and denoised using the basic particle-filter algorithm [25].

6.1. Principle of Particle-Filter Algorithm

The design idea of the standard particle-filter algorithm is to approximate the probability density function of the system with some discrete random sampling points and replace the integration operation with the sample mean to obtain the minimum variance estimate of the state. First, a set of random particle samples was generated according to the prior condition of the system state vector, and the weight of each sampled particle was calculated. Then, according to the system observation information, the particle weights and positions were continuously corrected, and the corrected particles were used to approximate the posterior probability density function of the target state, so that the approximate posterior probability density function can be used to estimate the target state. Finally, the particle filter state was output, namely
x ^ i ( k ) = i = 1 n ω i ( k ) x i ( k )
P i ( k ) = ω i ( k ) [ x i ( k ) x ^ i ( k ) ] [ x i ( k ) x ^ i ( k ) ] T
where x i ( k ) is the state of ith particle iterated at time k, x ^ i ( k ) is the weighted estimation state of the particle, P i ( k ) is the estimated variance, and ω i ( k ) is the particle weight, k = 1, 2, …, T, i = 1, 2, …, n.

6.2. Denoising Processing of Sampling Data Based on Particle Filter

Particle propagation was achieved by sampling the ZigBee-GPS/BDS multi-sensor state-transition model q ( x ( k )   x ( k 1 ) ,   z ( k ) ) and generating new particle states x i ( k ) . In order to reduce the interference of environmental noise on the positioning data, BPF was used to filter and denoise the multi-sensor positioning data [26], and the latest data collected by sensors was substituted into the observation equation to establish a new particle filter observation equation, which improves the reliability of multi-sensor data samples. The specific steps were as follows:
Step 1 The filtering model of the multi-sensor of the unmanned ship was established, and the state and measurement model of the multi-sensor system was defined as follows:
x ( k ) = f ( x ( k 1 ) ) + λ ( k 1 ) z ( k ) = h ( x ZGB ( k ) , x ˜ ( k ) ) + ν ( k )
where x(k) is the system state, z ( k ) is the sensor recursive positioning data, f(·) is the state transition function, h(·) is the measurement function, x(k − 1) is the system state at the last moment, x ZGB ( k ) is the sensor positioning data at time k, x ˜ ( k ) is the sensor positioning data weighted by the variance at time k, λ(k − 1) is the estimation noise, and ν(k) is the measurement noise.
Step 2 initialization: the prior density p(x0) was randomly sampled and the initialization particle set x(0) was generated.
Step 3 importance sampling.
(a) Randomly select n particle samples that satisfy the following distribution from the importance density function:
x i ( k ) ~ q ( x ( k ) | x ( k 1 ) ,   z ( k ) ) = p ( x ( k ) | x ( k 1 ) )
(b) Calculate the weights of the sampled particles and update them as follows:
ω i ( k ) = ω i ( k 1 ) ( p z ( k ) | x i ( k ) p x i ( k ) | x i ( k 1 ) / q x i ( k ) | x i ( k 1 ) , z ( k )
where ω i ( k ) is the ith particle weight at time k.
(c) Normalized the importance weights.
ω ˜ i ( k ) = ω i ( k ) [ i = 1 n ω i ( k ) ] 1
Step 4 Resampling.
Calculate the effective particle number Neff = ( i = 1 n ( ω i ( k ) ) 2 ) 1 . If Neff is less than the threshold value Nth, the particle set { x i ( k ) , ω ˜ i ( k ) } is resampled; otherwise, no resampling is required.
Step 5 Output the local estimation and covariance matrix of the filtered multi-sensor positioning data:
x ^ i ( k ) = i = 1 n ω ˜ i ( k ) x i ( k )   P i ( k ) = ω ˜ i ( k ) [ x i ( k ) x ^ i ( k ) ] [ x i ( k ) x ^ i ( k ) ] T
where x ^ i ( k ) is the local estimation of particle information, and P i ( k ) is the covariance matrix, k = 1, 2, …, T, i = 1, 2, …, n.

7. Sensor Data Fusion Based on TCPF

When the basic particle-filter algorithm is used to filter the ZigBee-GPS/BDS sensor state-transition function q ( x ( k )   x ( k 1 ) ,   z ( k ) ) , BPF can obtain a posterior probability that is close to the real state estimation because it is based on the Bayesian recursion of the sensor state prediction information according to the sampling idea. However, when filtering and fusing three or more input data, BPF will suffer from a low operating efficiency and lack of particle samples in the resampling process [25]. In order to meet the requirements of fusion filtering processing of the multi-input data model based on ZigBee-GPS/BDS, a particle-filter algorithm based on adaptive threshold and hierarchical capacity was proposed to perform fusion filtering processing of the multi-input data, so as to realize the precise positioning of the unmanned ship.

7.1. Principle of TCPF Algorithm

The proposed TCPF algorithm first used unscented transformation [27] to integrate the latest observation information into the proposal distribution, so that the proposal distribution was close to the real distribution of the probability density function: x i ( k )  ← TO( x i ( a ) ( k 1 ) ). Then, a Gaussian mixture model was constructed and an adaptive threshold was set to reduce the operation steps of clustering similar components in the Gaussian mixture and improve the clustering efficiency. Finally, the stratified sampling proportion capacity was set, the continuous probability density function was layered, and the combination of particle weights from the inferior layer was optimized to improve particle diversity. The specific steps are as follows:
Step 1 integrate the latest observation information into the proposal distribution.
(a) Extract the particle state and covariance matrix processed by the basic particle-filter algorithm, and calculate the sigma point set:
x i ( a ) ( k 1 ) = x ¯ i ( a ) ( k 1 ) ,   x ¯ i ( a ) ( k 1 ) ± ( n a + λ ) P i ( a ) ( k 1 )
where x i ( a ) ( k 1 ) is the sigma point set, x ¯ i ( a ) ( k 1 ) is the sigma point set after the unscented transformation, P i ( a ) ( k 1 ) is the covariance matrix of the sigma point set, na is the dimension of the sigma point set, and λ is the scale parameter.
(b) Integrate ZigBee-GPS/BDS multi-sensor positioning data into the obtained Sigma sampling point set, and update the system status and covariance.
x i ( k ) = x i ( a ) ( k 1 ) + k k z ( k 1 )
p i ( k ) = p i ( a ) ( k 1 ) k k p i ( z ( k 1 ) ) ( k 1 ) k k T
where x i ( k ) is the update state of the system, p i ( k ) is the update covariance of the system, kk is the Kalman gain, and p i ( z ( k 1 ) ) ( k 1 ) is the covariance matrix of the positioning data.
(c) Construct a proposal distribution that is closer to the target probability function using the system state and covariance, and sample from it.
x i ( k ) ~ q ( x ( k )   x ( k 1 ) ,   z ( k ) ) = N ( x i ( k ) , p i ( k ) )
where N(·) is the Gaussian function.
Step 2 construct a Gaussian mixture model.
(a) Generate a posterior probability density function p h x ( k ) | z ( k ) with time step k according to the Gaussian mixture components.
p h x ( k ) | z ( k ) = i = 1 C ( k ) ξ i N ( x ( k ) | m i , v i )
where N(x(k) | mi, vi) is the ith component in the mixed Gaussian model, C(k) is the number of component units of discrete samples, and ξ is the number of component units for discrete samples, k = 1, 2, …, T, i = 1, 2, …, n.
(b) Integrate the discrete sampling points sampled by Equation (50) and their corresponding weight { x i ( k ) , ω i ( k ) } into the Gaussian mixture component unit of Equation (51), and use the reconstructed continuous posterior probability density function p g x ( k ) | z ( k ) to resample the discrete particles:
p g x ( k ) | z ( k ) = i = 1 n ω i ( k ) N ( x i ( k ) | x i ( k 1 ) , h p ( k ) )
where
h = 0.5 N 2 / n x
p ¨ ( k ) = i = 1 n ω i ( k ) ( x i ( k ) x ¨ i ( k ) ) ( x i ( k ) x ¨ i ( k ) ) T
x ¨ i ( k ) = ( i = 1 n ω i ( k ) x i ( k ) ) / n
where p ¨ ( k ) is the covariance of the discrete particle filter distribution, x ¨ i ( k ) is the mean value of the discrete particle filter distribution, h is the normalized constant, and nx is the particle distribution dimension.
(c) Merge similar units of the Gaussian mixture in p g x ( k ) | z ( k ) using cluster analysis.
Step 3 set the adaptive threshold Tc for merging similar units in cluster analysis.
(a) Take the particle xm with the largest weight in the discrete particle sample set as the cluster center, and calculate the Mahalanobis distance Di between the other particles i and xm after selecting the importance sampling process.
D i = ( x m β c i ) T S 1 ( x m β c i )
where β c i is the probability density of particle i, and S is the covariance matrix.
(b) Calculate the number of effective particle samples Ne in the cluster unit.
N e = n 1 + σ β 2 c
where n is the number of particle samples and σ β c 2 is the covariance of particle probability density.
(c) Construct threshold T:
T = T 0 + k e N e R
where T0 is the initial value of the threshold, ke is the proportion coefficient, and R is the classification times.
(d) Substitute Equation (57) into Equation (58) to obtain the adaptive threshold Tc.
T c = T 0 + k e ( 1 + σ β 2 c ) n R
(e) Compare Di with the adaptive threshold Tc. If Di is less than or equal to Tc, the particle is classified into the component unit related to its probabilistic mass; on the contrary, skip the particle and cluster other particles.
(f) Select the particle with the largest weight from the remaining particle samples as the cluster center, and repeat step (e) until the clustering ends.
(g) Substitute the clustered component units into the continuous probability density function p ^ ( x ( k ) | z ( k ) ) of the constructed particle set. The p ^ ( x ( k ) | z ( k ) ) is expressed as follows:
p ^ ( x ( k ) | z ( k ) ) = i = 1 C ( k ) β i N ( x i ( k ) | γ i , p i ) ,   i = 1 C ( k ) β i = 1
where βi is the probability mass of the similar component i, γi is the mean of component i, and pi is the covariance of component i, i = 1, 2, …, n.
Step 4 set the hierarchical sampling proportional capacity.
(a) According to the layering theory [28], the continuous probability density function p ^ ( x ( k ) | z ( k ) ) is divided into l layers, and the probability density function of each layer is defined as pi(x). According to its probability quality, the component layers are divided into a group of weight advantage layers and two groups of disadvantage layers, and are defined as la, lb and lc, respectively.
(b) Set the proportional capacity of the particle number in layers la, lb and lc as n/4, n/3 and n/3, respectively.
(c) The particles whose weights are less than the average value ϖ ( k ) in lb and lc layers are optimized and combined to obtain the optimized particle weights ψ i ( k ) , and the sample data is sampled hierarchically. The ϖ ( k ) and ψ i ( k ) are calculated as follows:
ϖ ( k ) = 1 n i = 1 n ω i ( k )
ψ i ( k ) = η 1 η ω i ( k ) + 1 η ϖ ( k ) , i f   ω i ( k ) < ϖ ( k ) ,   Divide   particles   into   l b   layer ψ i ( k ) = ω i ( k ) , i f   ω i ( k ) ϖ ( k ) ,   Divide   particles   into   l a   layer
where η is the proportional coefficient.
(d) Output the sampling results of multi-sensor data fusion.
In the process of importance sampling, the current positioning data was integrated into the design of the proposal distribution of particle sets, which made the proposal distribution closer to the real posterior probability density and improved the estimation performance of the algorithm. By constructing Gaussian mixture probability density function instead of resampling, and constructing adaptive threshold in the cluster analysis of Gaussian mixture, the discrete particle samples were merged into similar component units, which reduces the complexity of clustering operation, and improves the real-time performance and operational efficiency of system signal processing. In addition, the continuous probability density function was stratified and the proportional capacity was set to ensure that there are enough particles in the inferior layer for weight optimization combination, which improved the diversity of particles, effectively suppressed the lack of particle samples, and improved the fusion processing accuracy of the TCPF algorithm for multi-input data.

7.2. Association of Stratified Sampling with Sensor Confidence

In order to further improve the fusion filtering accuracy of the algorithm, a new design was made for the hierarchical sampling step of the TCPF algorithm in this paper. The confidence factor ς i ZGB ( k ) of Equation (30) was associated with the hierarchical sampling of the TCPF algorithm, and a new weight optimization equation was obtained. The specific steps are as follows:
Step 1 substitute the confidence factor into the calculation of the weight optimization.
ψ i ( k ) = ς i ZGB ( k ) 1 ς i ZGB ( k ) ω i ( k ) + ς i ZGB ( k ) ϖ ( k ) , ω i ( k ) < ϖ ( k )
Step 2 the data samples are stratified according to the weights of multi-sensor sampling particles after optimized combination.
ψ i ( k ) = ς i ZGB ( k ) 1 ς i ZGB ( k ) ω i ( k ) + ς i ZGB ( k ) ϖ ( k ) , i f   ω i ( k ) < ϖ ( k ) ,   Divide   particles   into   l b   layer ψ i ( k ) = ω i ( k ) , i f   ω i ( k ) ϖ ( k ) ,   Divide   particles   into   l a   layer
The confidence factor of the positioning data was associated with the stratified sampling in the TCPF algorithm, so that the sensor positioning data with a large confidence factor would be sampled first when the TCPF algorithm was used to sample and fuse the positioning data, which would improve the sampling efficiency of the TCPF algorithm. The larger the confidence factor data, the higher the credibility, which improved the reference value of samples from sensor information fusion, further improved the fusion filtering accuracy of the TCPF algorithm, and finally realized the precise positioning of the unmanned ship.

7.3. Steps of Sensor Data Fusion Based on TCPF

Step 1 generate an initial particle sample set {xi(0), i = 1, 2, …, n} by sampling from the state transition function q ( x ( k ) |   x ( k 1 ) ,   z ( k ) ) of the ZigBee-GPS/BDS multi-sensor.
Step 2 perform basic UPF operations: x i ( k )  ← TO( x i ( a ) ( k 1 ) ), k ∈ [1,T], i ∈ [1,n].
Step 3 sample particles from the proposal distribution N ( x i ( k ) , p i ( k ) ) :
x i ( k ) ~ q ( x ( k ) |   x ( k 1 ) ,   z ( k ) ) = N ( x i ( k ) , p i ( k ) )
Step 4 construct an adaptive threshold Tc, and a Gaussian mixture continuous probability density function p ^ ( x ( k ) | z ( k ) ) using cluster analysis.
Step 5 the continuous probability density function is sampled hierarchically, and the sampling space is divided into la, lb, and lc layers.
Step 6 set the proportional capacity, and divide the particles into la, lb, and lc layer groups according to the weight of the confidence factor ς i ZGB of the sampled particles. At the same time, the weights of particles in the inferior layer lb and lc layer groups are combined and optimized, and they are added to the la layer to participate in sampling after the optimized particle weight ψ i ( k ) is obtained.
Step 7 output the fusion sampling results.
Step 8 k = k + 1, and return to step 2.

8. Numerical Simulation and Experimental Testing

8.1. Simulation of Data Fusion Algorithm Based on TCPF

In order to verify the performance of the proposed fusion positioning algorithm, 30 independent tests were carried out on the simulation model shown in Equations (63)–(65), and the test results were compared with the results of an Extended Kalman Filter (EKF) [29], Unscented Kalman Filter (UKF) [30], Unscented Particle Filter (UPF) [31] and Basic Particle Filter (BPF) [32] in terms of root-mean-square error (RMSE) and standard deviation (Std).
f 1 = x ( k ) = sin ( 0.6 π ( k 1 ) ) + 2 x ( k 1 ) + 6 + Q × r a n d n ( ) z ( k ) = 0.5 x 2 ( k ) + R × r a n d n ( )
f 2 = x ( k ) = 2.2 x ( k 1 ) + 3 x ( k 1 ) 2 + x 2 ( k 1 ) + cos ( 5 ( k 1 ) ) + Q × r a n d ( ) z ( k ) = 0.35 x 2 ( k ) 8 + R × r a n d ( )
f 3 = x ( k ) = x ( k 1 ) + 6 x ( k 1 ) 1 + x 2 ( k 1 ) + 4 cos ( 5 ( k 1 ) ) + Q × r a n d n ( ) z ( k ) = 0.25 x 2 ( k ) 3 + R × r a n d n ( )
The remaining parameters of TCPF were set as follows: N was 100, σ was 0.75, τ was 2; the process noise variance Q and measurement noise variance R were separately set according to different models, that is, the Q of the f1, f2 and f3 test models were Q1~N (0,0.5), Q2~N (0,1) and Q3~N(0,5), respectively; and the R of the f1, f2 and f3 test models were R1~N(0,1), R2~N(0,5) and R3~N(0,10), respectively. The parameters of the algorithm to be compared were taken from the corresponding references.
Table 1 shows the model test results of five algorithms. From the table, it can be seen that the mean and maximum values of RSME and Std from the proposed TCPF are smaller than the corresponding performances from the other four algorithms, which indicates that the filter fusion accuracy of TCPF was the highest. The superior performances were mainly due to the fact that in the process of importance sampling in this paper, the latest positioning data were integrated into the importance function through unscented transformation, which improved the credibility of particle samples. In addition, the algorithm in this paper clustered discrete particles by constructing an adaptive threshold in the Gaussian mixture and optimized the weights of particles in the inferior layer, which improved the diversity of data samples. The performance comparison shows that compared with the other four algorithms, the mean values of RMSE and Std of data fusion based on the proposed TCPF decreased by 25.0% and 28.0%, respectively, which fully shows that the TCPF data fusion algorithm can effectively suppress the particle shortage problem and improve the filtering accuracy.
Figure 4 shows the results of 30 RSME and Std tests of five algorithms against three types of models (f1, f2, f3). From the figure, it can be seen that when the noise variance was small, the RMSE and Std values of the TCPF data fusion algorithm were smaller than the values of the other four algorithms in each independent simulation experiment. When the noise variance increased, the RMSE and Std values of the five algorithms increased, which indicates that the filtering performance of the algorithms decreased, but the estimation accuracy of TCPF was still significantly higher than those of the other four algorithms. The test results in Figure 4 also verify that the filtering optimization performance of the TCPF data fusion algorithm was the best, and the filtering accuracy was the highest.

8.2. Experimental Testing and Analysis

8.2.1. Establishment of Unmanned-Ship Experimental Environment

In order to further verify the effectiveness of the multi-sensor data fusion positioning method of the unmanned ship for water quality testing based on the proposed TCPF, the ZigBee-GPS/BDS multi-sensor combined positioning system test platform as shown in Figure 1 and Figure 2 was established with the three-body unmanned ship as the carrier, and the unmanned ship navigation positioning experiment was carried out in the river channel of Zhangjiagang Campus of Jiangsu University of Science and Technology. In order to facilitate the experimental analysis and calculation, the starting point coordinate was subtracted from the preprocessed unmanned ship coordinate data, and the relative coordinate was used for testing. After calculation, the starting point coordinate of the unmanned ship was set to (3.8941, 0).
In order to evaluate the filtering accuracy of the unmanned-ship position information based on the proposed TCPF data fusion algorithm, first, the industrial computer equipped with the Linux system collected the position information of the unmanned ship using multiple sensors and outputted it from the serial port to the PC in the form of a log file through the wireless communication module. Then, the PC integrated the received position information of the unmanned ship into a data set [33]. After that, the designed TCPF data fusion algorithm was used to fuse and filter the position information data set of the unmanned ship, and finally realize the precise positioning of the unmanned ship. Figure 5 shows the experimental environment for the unmanned-ship positioning.
In order to verify the validity and superiority of the ZigBee-GPS/BDS multi-sensor positioning system for unmanned-ship navigation in different environments, the test river was divided into two sections, A and B, and two groups of ZigBee reference nodes were set along the A/B river according to the ZigBee networking principle. Among them, section A was located on the west side of the school library, which has an open water surface and was easy for the unmanned ship to navigate. However, the B section of the river was located on the west side of the school gymnasium, with dense vegetation and narrow channels, which had a certain environmental disturbance effect on the signal transmission of the ZigBee-GPS/BDS sensors. Figure 6 shows the schematic diagram of the node layout and navigation trajectory on the A/B river.

8.2.2. Unmanned-Ship Positioning Test Based on Multi-Sensor Data Fusion

In order to verify the superior performance of the proposed TCPF algorithm, some unmanned-ship navigation and positioning experiments based on the multi-sensor positioning system were carried out. According to the unmanned-ship positioning data collected by the sensors in the A/B experimental river, the test results of the six algorithms, namely TCPF, EKF, UKF, UPF, KF, and BPF, were compared against the three performance indicators of average positioning error, RMSE and Std. The filtering data of the proposed TCPF algorithm was ZigBee-GPS/BDS positioning data, and the sensor positioning data of the other five algorithms were all ZigBee-GPS. Figure 7 shows the filtering results of the unmanned-ship navigation trajectory based on six filtering algorithms.
From Figure 7b, it can be seen that when the unmanned ship sailed on the experimental river in section B, the sensor signal transmission was affected to a certain extent by environmental disturbances, such as dense trees and a narrow river. Except for the TCPF algorithm, the filtering and positioning results of other algorithms seriously deviated from the expected trajectory of the unmanned ship in some places; in particular, at the turning of the river, the positioning results were prone to large jumps, which caused the unmanned ship to gradually deviate from the expected trajectory and reduced the positioning accuracy and reliability. From Figure 7a, it can be seen that because the water surface of the experimental river in section A was more open than that in section B, and there were fewer shelters along the river, the positioning results of the other five algorithms showed that the unmanned ship could basically navigate according to the expected trajectory, but there was still a certain jump at the turning.
Combining Figure 7a,b, it can be seen that compared with the other five algorithms, the positioning accuracy of the unmanned ship based on the proposed TCPF algorithm was significantly higher. No matter whether at the turning of the river or in a straight line, the unmanned ship could basically navigate along the expected trajectory. This was mainly because the TCPF data fusion algorithm achieves the weighted correction of multi-sensor fault data through consistency inspection, which improves the fault-tolerant performance of the algorithm, and makes the TCPF algorithm have a high-reliability data sample for filtering. If the data are not inspected and corrected, when a type of sensor is subjected to a large environmental disturbance, the positioning data will be deviated, which will reduce the credibility of the data sample, and then reduce the positioning accuracy of the unmanned ship. Thus, this shows that the proposed TCPF algorithm of ZigBee-GPS/BDS not only ensures the accuracy of fault data inspection but also improves the fault-tolerant performance of the algorithm.
Table 2 and Table 3 show the comparison results of trajectory positioning for six algorithms. From the tables, it can be seen that the test results of the proposed TCPF algorithm were smaller than those of the other five algorithms, regardless of the average positioning error or the maximum positioning error of the A/B section, which indicates that the proposed method has the highest filtering and fusion accuracy. This is mainly because the posterior probability density function representing the state was approximated as a Gaussian mixture distribution in this paper, and particle samples were extracted on it instead of resampling, which ensured the diversity of ZigBee-GPS/BDS multi-sensor sampling particles, and avoided a lack of samples.
From the performances of Std and RMSE in the A/B section, compared with the other five algorithms, the value of the TCPF algorithm was also the smallest, which shows that the estimation performance of the TCPF algorithm is also the best. This is mainly due to the fact that the proposed TCPF algorithm integrated the latest multi-sensor positioning data collected by the unmanned-ship positioning system into the proposal distribution of the particle set in the importance sampling process, which made the proposal distribution closer to the real posterior probability density and improved the estimation performance of the algorithm. In addition, the integration operation in Bayesian estimation was solved according to the Monte Carlo method adopted by the basic particle-filter algorithm, which reduced the interference of environmental noise on the information collected by multi-sensors, and improved the credibility of multi-sensor positioning data samples.
The calculation shows that, compared with the other five algorithms, when the unmanned ship navigated based on TCPF in section A, the average positioning error, RMSE and, Std performances of its trajectory decreased by 35.0%, 36.0%, and 35.0%, respectively. In addition, the average positioning error, RMSE and, Std performance of its trajectory in section B decreased by 37.0%, 40.0%, and 39.0%, respectively. In the experiment involving the whole river, the three performances of TCPF decreased by an average of 36.0%, 38.0%, and 37.0%, respectively, which fully shows that the proposed TCPF algorithm had a high estimation performance and fault-tolerant performance, and the positioning accuracy of the unmanned ship’s navigation trajectory was improved.

8.2.3. Fault-Tolerance Test of the Unmanned-Ship Positioning Algorithm

In order to further verify the fault tolerance and filtering performance of the TCPF data fusion algorithm, two sets of fault-tolerance tests were carried out. One set of test data was from ZigBee-GPS, and the other set of selected test data was from ZigBee-BDS. In addition, the BPF algorithm was selected as the comparison algorithm for the performance test, and the positioning error, RMSE, and Std were used as performance indicators. The test data set of BPF was from ZigBee-GPS/BDS. Figure 8 shows the trajectory fusion results of the unmanned ship based on different algorithms in the A/B river.
From Figure 8, it can be seen that although the positioning results of the unmanned ship, based on the TCPF algorithm for different data sets, basically fitted the expected trajectory, the positioning accuracy of the TCPF algorithm for the ZigBee-GPS/BDS data set was significantly higher, and the positioning jump error at the river turn was also smaller, because the TCPF algorithm realized the information complementarity between the GPS and BDS data set. When a signal in the positioning system was lost, another sensor could supplement the collected position information of the unmanned ship to the sampling data set; thus, the sensor signal-loss problem caused by the large deviation in the unmanned ship position has been effectively overcome, which verifies that the constructed ZigBee-GPS/BDS multi-sensor positioning system performed the function of information complementarity between the data.
Combining the performance indicators of the average positioning error in Figure 8 and Table 4 and Table 5, it can be seen that the BPF algorithm using the ZigBee-GPS/BDS data set improved the positioning accuracy of the unmanned ship to a certain extent, but there was still a certain gap compared with the positioning accuracy of the TCPF algorithm, which was mainly due to the optimization performance of the TCPF algorithm itself. Firstly, the adaptive threshold was constructed in the cluster analysis of Gaussian mixture units, and the discrete particle samples were merged into similar component units, which reduces the complexity of clustering operations and improves the real-time performance of data fusion algorithms for signal processing. Secondly, the samples of the continuous probability density function of the constructed weighted point set were stratified, and the proportional capacity of each sampling layer was set, which ensured the reasonable distribution of the number of sampling particles in the stratification. Next, the layer group la was sampled, and the particle weights in the layer group lb and lc were optimally combined, which increased the probability quality and prevented the lack of samples.
From Table 4 and Table 5, it can be seen that the RMSE and Std performance indicators of the TCPF algorithm were better than those of the BPF algorithm. This was mainly because the confidence factor of multi-sensor positioning data is associated with hierarchical sampling in the TCPF algorithm, so when using the TCPF algorithm to filter and fuse ZigBee-GPS/BDS multi-sensor positioning data, the sensor positioning data with a large confidence factor will be fused and sampled first, which improves the sampling efficiency of the TCPF algorithm. Because the confidence factor reflects the credibility of the data, the larger the confidence factor was, the higher the reliability of the data sample was, which further improved the reference value of the sensor-information fusion sample, and further improved the fusion filtering accuracy of the TCPF algorithm for multi-input data, and finally achieved the accurate positioning of the unmanned ship.
The calculations show that, compared with the other three algorithms, the average positioning error, RMSE and Std performances of the unmanned ship, based on the TCPF algorithm during navigation in river A, decreased by 22.0%, 20.0%, and 17.0%, respectively, and the average positioning error, RMSE, and Std in river B decreased by 18.0%, 19.0%, and 18.0%, respectively. The three performances of the unmanned ship in the whole experimental river decreased by an average of 20.0%, 19.5%, and 17.5%, respectively, which fully demonstrates that the TCPF algorithm has a better optimization performance and fault-tolerant performance compared with the basic particle-filter algorithm, and improves the ship’s positioning accuracy under the interference of the unmanned-ship environment.

9. Discussion

Accurate positioning is the key to achieving the high-precision autonomous navigation of unmanned ships, and the positioning of unmanned ships based on multiple combined sensors has always been the focus of research by scholars. Data fusion based on two types of sensors is the current mainstream technology for unmanned-ship positioning, but there are shortcomings in affecting positioning accuracy and reliability due to the weakening or loss of certain sensor signals. In recent years, how to improve positioning accuracy by further fusing the data of three or more types of sensors has become a research hotspot. Adding sensors can help improve the reliability and accuracy of unmanned-ship positioning, but there is still a lack of research on how to further improve the accuracy and efficiency of data fusion through testing the confidence distance of data, as well as how to achieve effective localization under low or weak signals through data correction. In view of this, for the positioning of unmanned ships, in this paper, a multi-sensor positioning system incorporating ZigBee, GPS, and BDS was first built. Then, the particle-filter algorithm was introduced for denoising multi-sensor data. Next, the latest positioning data from multiple sensors were integrated into the proposed distribution using unscented transform to cause the proposed distribution to be close to the real distribution of the probability density function. Finally, the Gaussian mixture model was constructed, the adaptive threshold was set, and the data-confidence factor was associated with the hierarchical sampling to further improve the fusion filtering accuracy of the positioning data. The numerical test results of the three groups of models and the experimental test results of the navigation path all show that the proposed method not only realizes the mutual complementation, weighted correction, and confidence assignment of multi-sensor signals, but also significantly reduces their root-mean-square error and standard deviation compared with other filtering algorithms, which verifies that the proposed method has the advantages of high data reliability and a good filtering-fault-tolerance performance, and achieves the accuracy and reliable positioning of unmanned ships.

10. Conclusions and Future Work

In this paper, a multi-sensor data fusion positioning method for unmanned ships based on a threshold- and hierarchical-capacity particle filter was proposed to address the issues of significant deviation in data fusion accuracy caused by signal weakening and loss in conventional integrated positioning systems in unmanned ships. Through algorithm simulation and experimental testing, the following conclusions could be obtained:
(1)
The positioning data collected by the ZigBee-GPS/BDS multi-sensor is used to complement the information, which effectively overcomes the problem of sensor signal weakening or loss caused by environmental masking, and improves the accuracy and effectiveness of multi-sensor data fusion.
(2)
By conducting consistency checks on the multi-sensor positioning data and weighted correction of the faulty data, not only does it enhance the fault-tolerance performance of the data fusion algorithm but also strengthens the credibility of the data set samples.
(3)
The latest positioning data of multiple sensors are integrated into the proposal distribution of the particle set, which causes the suggested distribution to be closer to the true posterior probability density and improves the estimation performance of the algorithm. At the same time, the adaptive threshold is constructed in the cluster analysis of Gaussian mixing units, and the discrete particle samples are merged into similar component units, which improves the real-time performance and computing efficiency of the system’s signal processing.
(4)
Through stratified sampling and the setting of proportional capacity, a sufficient number of particles in the disadvantage layer are ensured for weight optimization and combination, and the diversity of particles is improved. Simultaneously, associating confidence factors with the TCPF algorithm’s layered sampling prioritizes the selection of positioning data with larger confidence factors for sample fusion during the fusion filtering process of the unmanned ships position information to enhance the efficiency and accuracy of the data fusion algorithm.
The multi-sensor data fusion method in this paper was mainly oriented to unmanned ships that conduct water quality detection in inland rivers. In recent years, unmanned ships have also been applied in various fields such as maritime rescue, maritime monitoring, and maritime cargo transportation. The marine environment is much more complex than inland rivers. Whether the multi-sensor fusion algorithm proposed in this paper is suitable for the marine environment will be the next research topic.

Author Contributions

Methodology, Y.S. and Z.Z.; validation, Z.Z. and M.Y.; original draft preparation, Y.S. and Z.Z.; writing, Z.Z., M.Y. and S.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the 2022 Research Plan for Innovative Products in the Industrial Chain of Zhangjiagang City (ZKC2206) and the 2022 Industry-University Research Pre-research Fund Project of Zhangjiagang City (ZKYY2253).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank Wei Liu of the T-SEA Marine Technology Co., Ltd. for his help in the sampling and laboratory analysis.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Nie, L.; Zhang, L.L. Key technologies and future development trends of unmanned ships. Mar. Equip. Mat. Mar. 2022, 30, 12–14. [Google Scholar] [CrossRef]
  2. Liu, S.D.; Liang, T.R.; Wang, Y.; Zhang, Y. Underwater moving target localization based on multi-modal information fusion. J. Nav. Pos. 2022, 10, 14–24. [Google Scholar] [CrossRef]
  3. Sheng, J.L. Application research of GNSS combined with unmanned ship technology in water system comprehensive management survey. Sci. Tech. Inf. 2021, 19, 52–54. [Google Scholar] [CrossRef]
  4. Wu, S.Y.; Zhu, H.S.; Fei, X.X.; Ma, Z.J.; Tang, M. Design and implement of unmanned boat system for water quality monitoring. Comput. Era 2022, 10, 93–97. [Google Scholar] [CrossRef]
  5. Deng, Z.L.; Wang, Z.G.; Sheng, W.; Liu, Y.P. Formation control of unmanned surface vehicles with position estimation algorithm. Mech. Sci. Technol. Aero. Eng. 2022, 41, 626–633. [Google Scholar] [CrossRef]
  6. Wang, X.Y.; Yang, C.; Liang, S.; Wang, A. Application of PPP and DGPS in IMU/GPS aided 1:10000 aerial survey project. Beijing Surv. Mapp. 2016, 6, 45–48. [Google Scholar] [CrossRef]
  7. Tang, C.; He, C.Y.; Dou, L.H. An IMU/ODM/UWB joint localization system based on modified cubature kalman filtering. Sensors 2021, 21, 4823. [Google Scholar] [CrossRef]
  8. Yousuf, S.; Kadri, M.B. Information fusion of GPS, INS and odometer sensors for improving localization accuracy of mobile robots in indoor and outdoor applications. Robotica 2020, 39, 48–52. [Google Scholar] [CrossRef]
  9. Hachemi, L.; Guiatni, M.; Nemra, A. Fault diagnosis and reconfiguration for mobile robot localization based on multi-sensors data fusion. Unmanned Syst. 2021, 41, 413–417. [Google Scholar] [CrossRef]
  10. Alam, S.A.; Gustafsson, O. Improved particle filter resampling architectures. J. Signal Process. Syst. 2020, 92, 555–568. [Google Scholar] [CrossRef]
  11. Wu, B.; Tian, Q. Improved indoor moving target localization algorithm based on UPF. Transducer Microsy. Technol. 2021, 40, 153–156+160. [Google Scholar] [CrossRef]
  12. Sun, H.R.; Ge, J.L.; Niu, J.D.; Wu, W.P. Positioning precision analysis of BDS-3, GPS and combined system. Pt. Eng. Technol. 2022, 59, 117–120. [Google Scholar] [CrossRef]
  13. Xu, H.Q.; Huang, K.; Yang, Y.; Zhu, S.N. Research and design of downhole positioning system based on ZigBee. Comput. Knowl. Technol. 2022, 18, 123–124+126. [Google Scholar] [CrossRef]
  14. Guo, N.; Han, H.Z.; Zhang, J. Research on single frequency PPP positioning performance of BDS-3/BDS-2/GPS with high sampling rate. Beijing Surv. Mapp. 2022, 36, 935–939. [Google Scholar] [CrossRef]
  15. Li, N.; Yu, M. Method of GPS coordinate transformation to Beijing 54 coordinate. Northen Commun. 2012, 8, 36–37. [Google Scholar] [CrossRef]
  16. Li, H.Y. Design of general ZigBee end device based on JN5169. Electron. Des. Eng. 2022, 30, 155–158. [Google Scholar] [CrossRef]
  17. Yakovlev, V.V.; Korzhenevskaya, I.V. The accounting possibilities for the earth rotation speed in satellite coordinate definitions. Mat. Sci. Eng. 2019, 698, 044007. [Google Scholar] [CrossRef]
  18. Chen, H.; Xu, C.H.; Gao, J.X.; Song, X.F.; Yuan, L.S. Precision analysis of pseudorange single point positioning by BDS, GPS and combined BDS/GPS. J. Shandong Univ. Sci. Technol. 2015, 34, 72–78. [Google Scholar] [CrossRef]
  19. Kimoto, R.; Ishida, S.; Yamamoto, T.; Tagashira, S.; Fukuda, A. MuCHLoc: Indoor ZigBee localization system utilizing inter-channel characteristics. Sensors 2019, 19, 1645. [Google Scholar] [CrossRef]
  20. Stig, E.; Marie, L. The impact of redundancy on reliability in machinery systems on unmanned ships. WMU J. Marit. Aff. 2022, 14, 312–317. [Google Scholar] [CrossRef]
  21. Zhang, M.; Lv, S.L. A method for determining confidence interval in train positioning. Railw. Signal Commun. 2018, 54, 1–4. [Google Scholar] [CrossRef]
  22. Peng, X.S.; Huang, X.H.; Qiu, W.N.; Wei, K.; Liu, S.W. A weighted fusion WiFi positioning algorithm based on distance rendezvous. Eng. Surv. Mapp. 2017, 26, 72–75. [Google Scholar] [CrossRef]
  23. Xu, D.; Guan, C.; Zhou, C.; Li, C. Research on consistency test of equipment maintainability pre-test data. Aerosp. Control. 2020, 38, 61–66. [Google Scholar] [CrossRef]
  24. Xu, L.; Chen, X.; Ding, Z.X.; Hu, J.; Song, T.C. Outdoor positioning method combining RSSI and TOA based on lora technology. Meas. Control Eng. 2020, 39, 41–46. [Google Scholar] [CrossRef]
  25. Wu, T.Z.; Zhao, T.; Xu, S.Y. Prediction of remaining useful life of the lithium-ion battery based on improved particle filtering. Front. Energy Res. 2022, 10, 95–101. [Google Scholar] [CrossRef]
  26. Zeng, R.; Tian, J.; Jiang, H.; Wang, Z.N. UAV track prediction method based on improved particle filtering. Transducer Microsyst. Technol. 2022, 41, 148–151. [Google Scholar] [CrossRef]
  27. Li, X.J.; Li, J.F.; Zhu, M.; Peng, N.L.; Zuo, S. Research on positioning fusion and verification algorithm based on UKF. Automot. Eng. 2021, 43, 825–832. [Google Scholar] [CrossRef]
  28. Wei, J.X.; Cheng, L.; Han, P.; Zhu, Y.X.; Huang, W.D. Decision tree-based data stratification method for the minimization of the masking effect in adverse drug reaction signal detection. Appl. Sci. 2021, 11, 122–128. [Google Scholar] [CrossRef]
  29. Ding, N.; Prasad, K.; Lie, T.T. State of charge estimation of a composite lithium-based Battery model based on an improved extended kalman filter algorithm. Inventions 2019, 4, 66. [Google Scholar] [CrossRef]
  30. Gong, Y.J.; Xie, L.Y.; Xue, S.T.; Tang, H.S. Adaptive unscented Kalman filter for nonlinear structural identification. J. Bldg. Str. 2023, 44, 182–191+224. [Google Scholar] [CrossRef]
  31. Yu, Z.; Liu, Y.; Wei, F.; Liu, L.J. Prediction of aeroengine exhaust gas temperature based on unscented particle filter algorithm. Aeroengine 2021, 47, 1–6. [Google Scholar] [CrossRef]
  32. Gao, Y.; Mao, Y.H.; Yang, Y. Research on artificial fish school algorithm-particle filter. Mod. Electron. Tech. 2021, 44, 170–174. [Google Scholar] [CrossRef]
  33. Rim, C.M.; Sin, Y.C.; Paek, K.H. A mobile robot localization method based on polar scan matching and adaptive niching chaos optimization algorithm. J. Intell. Robot. Syst. 2022, 106, 19. [Google Scholar] [CrossRef]
Figure 1. ZigBee-GPS/BDS multi-sensor combined positioning system for an unmanned ship.
Figure 1. ZigBee-GPS/BDS multi-sensor combined positioning system for an unmanned ship.
Applsci 13 10390 g001
Figure 2. Test platform for unmanned-ship positioning system.
Figure 2. Test platform for unmanned-ship positioning system.
Applsci 13 10390 g002
Figure 3. Particle filter-based multi-sensor data fusion positioning framework for unmanned ships.
Figure 3. Particle filter-based multi-sensor data fusion positioning framework for unmanned ships.
Applsci 13 10390 g003
Figure 4. RMSE and Std tests of five algorithms against different models: (a) f1 model, Q1~N(0,0.5), R1~N(0,1); (b) f1 model, Q1~N(0,0.5), R1~N(0,1); (c) f2 model, Q2~N(0,1), R3~N(0,5); (d) f2 model, Q2~N(0,1), R3~N(0,5); (e) f3 model, Q3~N(0,5), R5~N(0,10); and (f) f3 model, Q3~N(0,5), R5~N(0,10).
Figure 4. RMSE and Std tests of five algorithms against different models: (a) f1 model, Q1~N(0,0.5), R1~N(0,1); (b) f1 model, Q1~N(0,0.5), R1~N(0,1); (c) f2 model, Q2~N(0,1), R3~N(0,5); (d) f2 model, Q2~N(0,1), R3~N(0,5); (e) f3 model, Q3~N(0,5), R5~N(0,10); and (f) f3 model, Q3~N(0,5), R5~N(0,10).
Applsci 13 10390 g004aApplsci 13 10390 g004b
Figure 5. Experimental environment for unmanned ship positioning: (a) Communication node; (b) PC terminal.
Figure 5. Experimental environment for unmanned ship positioning: (a) Communication node; (b) PC terminal.
Applsci 13 10390 g005
Figure 6. Schematic diagram of node layout and navigation trajectory on A/B river: (a) ZigBee node layout; (b) navigation trajectory of unmanned ship.
Figure 6. Schematic diagram of node layout and navigation trajectory on A/B river: (a) ZigBee node layout; (b) navigation trajectory of unmanned ship.
Applsci 13 10390 g006
Figure 7. Filtering results of the unmanned-ship navigation trajectory based on five filtering algorithms: (a) positioning trajectories on section A; (b) positioning trajectories on section B.
Figure 7. Filtering results of the unmanned-ship navigation trajectory based on five filtering algorithms: (a) positioning trajectories on section A; (b) positioning trajectories on section B.
Applsci 13 10390 g007
Figure 8. Trajectory fusion results of the unmanned ship based on different algorithms in the A/B river: (a) section A positioning trajectory; (b) section B positioning trajectory.
Figure 8. Trajectory fusion results of the unmanned ship based on different algorithms in the A/B river: (a) section A positioning trajectory; (b) section B positioning trajectory.
Applsci 13 10390 g008
Table 1. Model test results of five algorithms.
Table 1. Model test results of five algorithms.
Test ModelMean Value of RMSEMean Value of Std
EKFUKFBPFUPFTCPFEKFUKFBPFUPFTCPF
f10.4930.4210.6460.3900.3450.4730.3930.6220.3810.235
f24.4332.5543.9723.3832.4854.2102.5742.7922.8932.367
f35.4563.5984.8124.5633.4875.4983.7023.7104.2363.396
Test ModelMaximum Value of RMSEMaximum Value of Std
EKFUKFBPFUPFTCPFEKFUKFBPFUPFTCPF
f10.9100.5680.9400.6130.4090.8770.5660.8910.5990.412
f26.6813.4596.1424.6243.3176.3913.4313.4524.2833.283
f37.8564.1926.1216.1423.9987.7314.2144.1725.9313.969
Table 2. Trajectory positioning test results in section A.
Table 2. Trajectory positioning test results in section A.
AlgorithmEKFUKFBPFUPFKFTCPF
Performance Index
Average positioning error3.2654.1103.5382.8213.8562.352
Maximum positioning error5.2216.0105.3194.1145.5892.941
RMSE3.0573.9563.4322.7193.8542.351
Std0.4350.5890.5420.3950.4980.345
Table 3. Trajectory positioning test results in section B.
Table 3. Trajectory positioning test results in section B.
AlgorithmEKFUKFBPFUPFKFTCPF
Performance Index
Average positioning error3.7864.5394.4643.6584.3652.875
Maximum positioning error6.5198.0187.6996.3937.8564.463
RMSE3.6584.3894.3583.5724.2152.723
Std0.9521.8531.2581.1521.5620.876
Table 4. Data positioning test results in section A.
Table 4. Data positioning test results in section A.
AlgorithmTCPF-Z + GTCPF-Z + BBPF-Z + G/BTCPF-Z + G/B
Performance Index
Average positioning error2.558 2.606 2.532 1.828
Maximum positioning error3.4384.0133.6632.676
RMSE2.541 2.648 2.588 1.860
Std0.368 0.605 0.381 0.332
Table 5. Data positioning test results in section B.
Table 5. Data positioning test results in section B.
AlgorithmTCPF-Z + GTCPF-Z + BBPF-Z + G/BTCPF-Z + G/B
Performance Index
Average positioning error2.754 2.864 2.727 2.092
Maximum positioning error4.5614.8754.9453.631
Positioning RMSE2.834 2.973 2.847 2.168
Positioning Std0.572 0.717 0.751 0.506
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

Shen, Y.; Zhao, Z.; Yuan, M.; Wang, S. Research on Multi-Sensor Data Fusion Positioning Method of Unmanned Ships Based on Threshold- and Hierarchical-Capacity Particle Filter. Appl. Sci. 2023, 13, 10390. https://doi.org/10.3390/app131810390

AMA Style

Shen Y, Zhao Z, Yuan M, Wang S. Research on Multi-Sensor Data Fusion Positioning Method of Unmanned Ships Based on Threshold- and Hierarchical-Capacity Particle Filter. Applied Sciences. 2023; 13(18):10390. https://doi.org/10.3390/app131810390

Chicago/Turabian Style

Shen, Yi, Zeyu Zhao, Mingxin Yuan, and Sun Wang. 2023. "Research on Multi-Sensor Data Fusion Positioning Method of Unmanned Ships Based on Threshold- and Hierarchical-Capacity Particle Filter" Applied Sciences 13, no. 18: 10390. https://doi.org/10.3390/app131810390

APA Style

Shen, Y., Zhao, Z., Yuan, M., & Wang, S. (2023). Research on Multi-Sensor Data Fusion Positioning Method of Unmanned Ships Based on Threshold- and Hierarchical-Capacity Particle Filter. Applied Sciences, 13(18), 10390. https://doi.org/10.3390/app131810390

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