# Ship Tracking Based on Underwater Electric Potential

Authors: Baoquan Sun; Bing Yan; Jiawei Zhang; Shouwei Hu
Journal: Mathematical Problems in Engineering
 (2018)

---

## Abstract

Underwater Electric Potential is an important signal characteristic of a ship. The signal contains location information which can be used to track the ship. This research tries to study the possibility of ship tracking using Underwater Electric Potential. Aiming at the problems existing in the traditional Kalman filters under large initial errors, a new nonlinear filter is proposed. State space model of ship tracking is established; the problem existing in the ordinary Kalman filters is analyzed from the perspective of Kullbeck-Leibler Divergence; the new algorithm is proposed based on progressive Bayesian; simulations are designed. Simulation results show: it is feasible to use underwater electric potential to track the ship; the new method can effectively improve filter performance under large initial error and can effectively track the ship with preferable precision and convergence, which has great practical value.


---

## Body

## 1. Introduction

Underwater Electric Potential (UEP) resulting from corrosion or anticorrosion current of a ship is an important signal that is not negligible. There is obvious distribution characteristic with UEP [1] which can be used to track the ship. Currently, studies on UEP mainly focus on the model construction [2, 3], and some literatures begin to explore the feasibility of using Static Electric Field (SE) to track ship [4]. With the development of vibration absorption, noise reduction, and electric propulsion technology, ship radiation noise level has been greatly reduced. Coupled with sound line bending, sound reverberation, and other natural conditions, passive acoustic tracking is becoming more difficult. However, UEP signal is less affected by the environment. Ship tracking using UEP can be an effective supplement to acoustic tracking. The present research aims to study the possibility of tracking ship using UEP.The most common method of target tracking is filtering, which has been studied extensively. Filtering method is an on-line recursive method based on the dynamic model and the observation model that can estimate the target state from the noisy observation data in real time. The main difficulty in the filtering method for ship tracking using UEP is that it is not easy to set filter initial value. However, the initial value is a key to affect the filter performance. In the ship tracking problem, the azimuth-distance information is hidden in the UEP signal, resulting in the initial value being difficult to set, and usually with large initial errors.Extended Kalman filter (EKF) usually provides good estimation performance, and it has been widely used. However, EKF fails in some applications when the observation model is strong nonlinear [5]. There are several other filters designed to deal with nonlinear problems. Unscented Kalman filter (UKF) [6] and cubature Kalman filter (CKF) [7], typical of the sigma-point class, are able to retain higher-order terms of the Taylor series by choosing various sigma points spread around the distribution function and constructing the filter gain from the sigma points. Another important set of filters are the high-order filters. The second-order filters (SOF) [8] are the best known among them. These filters approximate the models by a high-order Taylor approximation. However, the complexity and computational load of these filters make their implementation unfeasible, particularly in the case where the state variable is high dimensional [9, 10].Particle filters [11] approximate the full probability distribution, by using a finite set of samples. These filters use sequential sampling and resampling according to an “importance function” to provide much better performance than Kalman filters. The improved performance is especially pronounced when using nonlinear models or non-Gaussian distributions. However, the high computational cost associated with using the large number of samples required makes them unsuitable for underwater tracking applications.One way to solve the inherent problem of traditional filters is the use of progressive Gaussian filtering (PGF) [12], which gradually includes the measurement information instead of using all the information in one step. According to this idea, a new GA filter is derived through approximating intermediate progressive joint PDFs of state and measurement as Gaussian in reference [13]. Another way is the use of particle flow filters which propose an alternate convenient way of performing Bayesian updates [14] through probabilistic inference. Particle flow filters rely on particles to carry out, which leads to the same problem as particle filters. A progressive Bayesian method without particles is derived under the Gaussian assumption in [15]. Comparing [13] and [15], we can find that, under the Gauss hypothesis, two methods are consistent in form. In fact, both class filters are just different manifestations of progressive Bayesian.According to the previous analysis, Kalman filters are more suitable for underwater applications than particle filters. Therefore, consider applying the progressive Bayesian idea to the Kalman filters to improve their performance. Ship motion is of low maneuverability, and its dynamic model can be assumed as a uniform linear motion model. Therefore, the essence of applying progressive Bayesian to Kalman filters lies in redesigning measurement update. Kalman filters perform one measurement update when a measurement is received. The idea behind progressive Bayesian is to introduce the measurement information progressively in pseudo time which equals performing the measurement update of Kalman filters multiple times. So the new methods are named as progressive update Kalman filters (PU-KFs).The remainder of this paper is organized as followed. In Section2, we present the ship tracking problem with UEP and analyze the problem existing in the ordinary Kalman filters from the perspective of Kullbeck-Leibler Divergence (KLD). In Section 3, we derive the proposed algorithm based on progressive Bayesian. In Section 4, we design simulations to verify related theory. In Section 5, we draw a simple conclusion for the paper.


## 2. Problem Statement

### 2.1. State Space Model

The general tracking model consists of the dynamic model and observation model(1)xk=axk-1+wk-1yk=hxk+vkwhere xk∈Rn is the n-dimensional state vector; yk∈Rm is the m-dimensional observation vector; a and h are some known functions; wk-1 and vk are white noises and are independent, wk-1~N(0,Qk-1), vk~N(0,Rk).Ship UEP is mainly modeled by boundary element method [16], finite element method [17] and equivalent source method [18]. Boundary element method and finite element method cannot reverse the source information from electric field information. By comparison, the equivalent source method has an exact mathematical model. Among the various equivalent source models, point current (a point where current flows in or out) model is closest to the actual ship model and has higher precision [19]. Therefore, the ship UEP signal is modeled by point current method.According to point current method the ship can be described as a uniform linear array (referred to as point current array) withN(N≥2) point currents at equal distance. The current density of each point current is Ipi and the distance is ld. In general, the larger N, the higher accuracy of the target’s signal fitting, but at the same time it will lead to an increase in computational complexity. To this end, the choice of N is related to the scale of the target and our requirements for computational accuracy. The UEP at a filed point is the sum of the UEP produced by the N point currents at this field point.(2)U=∑i=1NIpiKIpi,Pfwhere K(Ipi,Pf) is the distance function between the point current and the field point. In the air-seawater-seabed three-layer uniform media (as shown in Figure 1), K(Ipi,Pf) is (3)KIpi,Pf=14πσ1r2+z-zi2+1r2+z+zi-h2+∑m=1∞km1r2+z+zi-h+2mH2+1r2+z+zi-h-2mH2+1r2+z-zi+2mH2+1r2+z-zi-2mH2where H is the depth of seawater; h is the depth of sensor; r2=(x-xi)2+(y-yi)2; σ is the seawater conductivity; σ1 is the seabed conductivity; k=(σ-σ1)/(σ+σ1) is the seabed reflection coefficient; m is the number of reflections; and the limit in calculation is usually 10~20 [20].Figure 1
Air-seawater-seabed three-layer uniform media.UEP measured at the N0.j electric field sensor can be modeled as(4)yjxk=hjxk+vkj=Ukj+vkjwhere vk(j) is the measured noise.And then, the observation equation is(5)yk=hxk+vk=Uk1Uk2⋮UkJ+vk1vk2⋮vkJThe state vector contains the target position, velocity, and point current array, while the motion of the ship’sZ direction is neglected.(6)xk=rk,Vk,Ip1:NTwhere rk=[x,y,z]T, Vk=[Vx,Vy]T, and Ip1:N=[Ip1,Ip2,…,IpN]T.Ship target has low dynamic characteristics and its motion state can be modeled as Discrete White Noise Acceleration (DWNA) model [21]. The motion equation is(7)xkr,V=Φmxk-1r,V+Γmwk-1r,Vwhere Φm=I3×3TsI2×201×202×3I2×2 is motion state transition matrix, Γm=Ts2/2I2×201×2TsI2×2 is noise gain matrix, and Ts is sampling interval.The DWNA model assumes that the acceleration is Gaussian white noise acceleration, i.e.,(8)Σ=Ewk-1r,Vwk-1r,VT=diagσx2,σy2where σx and σy are the process acceleration noise intensities.Point current array is a constant matrix with a small noise. The whole state equation is(9)xk=Fxk-1+wk-1where F=Φm05×N0N×5IN×N; noise covariance matrix is Q=ΓmΣΓmαIN×N; and α is a small constant.


### 2.2. Moment Approximation Analysis

For nonlinear problems, the key that affects the Kalman-type filter algorithm is the approximation of the following observation characteristic moments.(10)h^-=∫RmhxNx;x0,P0dxPhh-=∫Rmhx-h^-hx-h^-TNx;x0,P0dxwhere h^- is the prior observation; Phh- is the corresponding covariance. Formula (10) is the mean and covariance of the approximate Gauss distribution after the nonlinear transformation. Two typical moment approximation methods are considered:(  1) First Order (T1) Taylor Expansion:(11)h~-=hx0P~hh-=H0P0H0Twhere H0=∂h/∂xx=x0 is the Jacobian matrix.(  2) Unscented Transformation (UT) or Cubature Transformation (CT) (12) h ~  - = ∑ j = - n n ω j Y j P ~ h h - = ∑ j = - n n ω j Y j - h ~  - Y j - h ~  - Twhere Y±j=h(X±j),j=0,1,⋯,n; X±j is the deterministic sampling point; and (13)Xo=x0,X±i=x0±n+λP0ω0=λn+λ,ω±i=12n+λwhere i=1,⋯,n; for UT, λ=3-n; for CT λ=0.We can get EKF, UKF, and CKF through these ways of approximation. For a general nonlinear model, suppose the Gaussian moment in the formula (10) could be solved accurately and the KLD between the joint Gaussian hypothesis distribution and the true distribution is [22](14)Ip,N=log⁡I+R-1Pk∣k-1yy-Pk∣k-1xyTPk∣k-1-1Pk∣k-1xyWhen the priori estimation errorPk∣k-1 increases, the KLD between the approximate distribution and true distribution will increase, meaning that the Gaussian distribution is not a valid assumption any more, which will cause the filter to degrade or even divergence under large prior error conditions.


## 2.1. State Space Model

The general tracking model consists of the dynamic model and observation model(1)xk=axk-1+wk-1yk=hxk+vkwhere xk∈Rn is the n-dimensional state vector; yk∈Rm is the m-dimensional observation vector; a and h are some known functions; wk-1 and vk are white noises and are independent, wk-1~N(0,Qk-1), vk~N(0,Rk).Ship UEP is mainly modeled by boundary element method [16], finite element method [17] and equivalent source method [18]. Boundary element method and finite element method cannot reverse the source information from electric field information. By comparison, the equivalent source method has an exact mathematical model. Among the various equivalent source models, point current (a point where current flows in or out) model is closest to the actual ship model and has higher precision [19]. Therefore, the ship UEP signal is modeled by point current method.According to point current method the ship can be described as a uniform linear array (referred to as point current array) withN(N≥2) point currents at equal distance. The current density of each point current is Ipi and the distance is ld. In general, the larger N, the higher accuracy of the target’s signal fitting, but at the same time it will lead to an increase in computational complexity. To this end, the choice of N is related to the scale of the target and our requirements for computational accuracy. The UEP at a filed point is the sum of the UEP produced by the N point currents at this field point.(2)U=∑i=1NIpiKIpi,Pfwhere K(Ipi,Pf) is the distance function between the point current and the field point. In the air-seawater-seabed three-layer uniform media (as shown in Figure 1), K(Ipi,Pf) is (3)KIpi,Pf=14πσ1r2+z-zi2+1r2+z+zi-h2+∑m=1∞km1r2+z+zi-h+2mH2+1r2+z+zi-h-2mH2+1r2+z-zi+2mH2+1r2+z-zi-2mH2where H is the depth of seawater; h is the depth of sensor; r2=(x-xi)2+(y-yi)2; σ is the seawater conductivity; σ1 is the seabed conductivity; k=(σ-σ1)/(σ+σ1) is the seabed reflection coefficient; m is the number of reflections; and the limit in calculation is usually 10~20 [20].Figure 1
Air-seawater-seabed three-layer uniform media.UEP measured at the N0.j electric field sensor can be modeled as(4)yjxk=hjxk+vkj=Ukj+vkjwhere vk(j) is the measured noise.And then, the observation equation is(5)yk=hxk+vk=Uk1Uk2⋮UkJ+vk1vk2⋮vkJThe state vector contains the target position, velocity, and point current array, while the motion of the ship’sZ direction is neglected.(6)xk=rk,Vk,Ip1:NTwhere rk=[x,y,z]T, Vk=[Vx,Vy]T, and Ip1:N=[Ip1,Ip2,…,IpN]T.Ship target has low dynamic characteristics and its motion state can be modeled as Discrete White Noise Acceleration (DWNA) model [21]. The motion equation is(7)xkr,V=Φmxk-1r,V+Γmwk-1r,Vwhere Φm=I3×3TsI2×201×202×3I2×2 is motion state transition matrix, Γm=Ts2/2I2×201×2TsI2×2 is noise gain matrix, and Ts is sampling interval.The DWNA model assumes that the acceleration is Gaussian white noise acceleration, i.e.,(8)Σ=Ewk-1r,Vwk-1r,VT=diagσx2,σy2where σx and σy are the process acceleration noise intensities.Point current array is a constant matrix with a small noise. The whole state equation is(9)xk=Fxk-1+wk-1where F=Φm05×N0N×5IN×N; noise covariance matrix is Q=ΓmΣΓmαIN×N; and α is a small constant.


## 2.2. Moment Approximation Analysis

For nonlinear problems, the key that affects the Kalman-type filter algorithm is the approximation of the following observation characteristic moments.(10)h^-=∫RmhxNx;x0,P0dxPhh-=∫Rmhx-h^-hx-h^-TNx;x0,P0dxwhere h^- is the prior observation; Phh- is the corresponding covariance. Formula (10) is the mean and covariance of the approximate Gauss distribution after the nonlinear transformation. Two typical moment approximation methods are considered:(  1) First Order (T1) Taylor Expansion:(11)h~-=hx0P~hh-=H0P0H0Twhere H0=∂h/∂xx=x0 is the Jacobian matrix.(  2) Unscented Transformation (UT) or Cubature Transformation (CT) (12) h ~  - = ∑ j = - n n ω j Y j P ~ h h - = ∑ j = - n n ω j Y j - h ~  - Y j - h ~  - Twhere Y±j=h(X±j),j=0,1,⋯,n; X±j is the deterministic sampling point; and (13)Xo=x0,X±i=x0±n+λP0ω0=λn+λ,ω±i=12n+λwhere i=1,⋯,n; for UT, λ=3-n; for CT λ=0.We can get EKF, UKF, and CKF through these ways of approximation. For a general nonlinear model, suppose the Gaussian moment in the formula (10) could be solved accurately and the KLD between the joint Gaussian hypothesis distribution and the true distribution is [22](14)Ip,N=log⁡I+R-1Pk∣k-1yy-Pk∣k-1xyTPk∣k-1-1Pk∣k-1xyWhen the priori estimation errorPk∣k-1 increases, the KLD between the approximate distribution and true distribution will increase, meaning that the Gaussian distribution is not a valid assumption any more, which will cause the filter to degrade or even divergence under large prior error conditions.


## 3. Progressive Update Kalman Filters

### 3.1. New Design for Measurement Update of Kalman Filters

The new algorithm focuses on redesigning measurement update. Time interval[tk-1,tk] is divided into Npu subintervals, and the interval of each subinterval is δλ=(tk-tk-1)/Npu, which amount to obtaining Npu observations sequentially in [tk-1,tk]. A measurement update process is performed in each of the subinterval. For time interval changes, observation noise covariance will also change.It is easy to shown that the variance of the mean ofNpu random independent zero-mean samples is equal to the variance of the individual samples divided by Npu, that is [23],(15)E1Npu∑i=1Npuri2=1Npu2Er1+r2+…rNpu2=1Npu2Er12+r22+…rNpu2=1Npuσr2where Eri=0, Erirj=σr2δij. Then the new observation noise covariance and the original noise covariance have the following relation: (16)Rioriginal=1NpuRinew⇔Rinew=NpuRioriginalIn the traditional Kalman filters, the EKF is the simplest in form. Therefore take PUEKF as a representative for PU-KF. Measurement update for PUEKF in[(i-1)δλ,iδλ] is as follows (δλ is omitted in the subscript):(17)Si-1=H^i-1Pi-1H^i-1T+NpuRkKi=Pi-1H^i-1TSi-1-1x^i=x^i-1+Kiyk-hx^i-1Pi=I-KiH^i-1Pi-1In[tk-1,tk] the measurement update will be executed Npu times. Then the PUEKF measurement update is as follows.Steps 1. Initialization, x^0=x^k∣k-1, P0=Pk∣k-1.Steps 2. i=1,2,⋯,NpuFollow formula (17) to carry out measurement update.Steps 3. x^k∣k=x^Npu, Pk∣k=PNpu.WhenNpu=1, the PUEKF algorithm degrades to the EKF algorithm.From formula (17), we can see (18)Pk∣k-1Npu-1=∏I-KHPk∣k-1≤Pk∣k-1Through the analysis of formula (14) we know a large Pk∣k-1 is not conducive to filter stability, and the new algorithm reduces KLD by reducing Pk∣k-1, increasing the stability of the algorithm.


### 3.2. Linear Constraint for PUEKF

To match the point current array model to the actual situation, it is necessary to add linear constraint its state variablesx: (19)Dx=dwhere D=diag01×5,I1×N, d≈0.The state is estimated as [24](20)x~=argminEx-x~2∣Y,Dx~=dFormula (20) ignores the time subscripts for random variables. The condition expectation is (21)Ex-x~2∣Y=∫x-x~Tx-x~pdfx∣Ydx=∫xTxpdfx∣Ydx-2x~∫xpdfx∣Ydx+x~Tx~The Lagrange expression is(22)L=Ex-x~2∣Y+2λLTDx~-d=∫xTxpdfx∣Ydx-2x~∫xpdfx∣Ydx+x~Tx~+2λLTDx~-dTo get the extreme value, the Lagrange partial derivative is set to be 0, i.e.,(23)∂L∂x~=-2x^+2x~+2DTλL=0∂L∂λ=Dx~-d=0where x^=E(x∣Y)=∫xpdf(x∣Y)dx.λ L and x~ are available(24)λL=DDT-1Dx^-dx~=x^-DTDDT-1Dx^-dAdd Steps 4 to PU-EKF measurement update.Steps 4. x^k∣k=x^k∣k-DTDDT-1(Dx^k∣k-d)We get the complete measurement update steps for the linear constraint PUEKF (LPUEKF).


## 3.1. New Design for Measurement Update of Kalman Filters

The new algorithm focuses on redesigning measurement update. Time interval[tk-1,tk] is divided into Npu subintervals, and the interval of each subinterval is δλ=(tk-tk-1)/Npu, which amount to obtaining Npu observations sequentially in [tk-1,tk]. A measurement update process is performed in each of the subinterval. For time interval changes, observation noise covariance will also change.It is easy to shown that the variance of the mean ofNpu random independent zero-mean samples is equal to the variance of the individual samples divided by Npu, that is [23],(15)E1Npu∑i=1Npuri2=1Npu2Er1+r2+…rNpu2=1Npu2Er12+r22+…rNpu2=1Npuσr2where Eri=0, Erirj=σr2δij. Then the new observation noise covariance and the original noise covariance have the following relation: (16)Rioriginal=1NpuRinew⇔Rinew=NpuRioriginalIn the traditional Kalman filters, the EKF is the simplest in form. Therefore take PUEKF as a representative for PU-KF. Measurement update for PUEKF in[(i-1)δλ,iδλ] is as follows (δλ is omitted in the subscript):(17)Si-1=H^i-1Pi-1H^i-1T+NpuRkKi=Pi-1H^i-1TSi-1-1x^i=x^i-1+Kiyk-hx^i-1Pi=I-KiH^i-1Pi-1In[tk-1,tk] the measurement update will be executed Npu times. Then the PUEKF measurement update is as follows.Steps 1. Initialization, x^0=x^k∣k-1, P0=Pk∣k-1.Steps 2. i=1,2,⋯,NpuFollow formula (17) to carry out measurement update.Steps 3. x^k∣k=x^Npu, Pk∣k=PNpu.WhenNpu=1, the PUEKF algorithm degrades to the EKF algorithm.From formula (17), we can see (18)Pk∣k-1Npu-1=∏I-KHPk∣k-1≤Pk∣k-1Through the analysis of formula (14) we know a large Pk∣k-1 is not conducive to filter stability, and the new algorithm reduces KLD by reducing Pk∣k-1, increasing the stability of the algorithm.


## 3.2. Linear Constraint for PUEKF

To match the point current array model to the actual situation, it is necessary to add linear constraint its state variablesx: (19)Dx=dwhere D=diag01×5,I1×N, d≈0.The state is estimated as [24](20)x~=argminEx-x~2∣Y,Dx~=dFormula (20) ignores the time subscripts for random variables. The condition expectation is (21)Ex-x~2∣Y=∫x-x~Tx-x~pdfx∣Ydx=∫xTxpdfx∣Ydx-2x~∫xpdfx∣Ydx+x~Tx~The Lagrange expression is(22)L=Ex-x~2∣Y+2λLTDx~-d=∫xTxpdfx∣Ydx-2x~∫xpdfx∣Ydx+x~Tx~+2λLTDx~-dTo get the extreme value, the Lagrange partial derivative is set to be 0, i.e.,(23)∂L∂x~=-2x^+2x~+2DTλL=0∂L∂λ=Dx~-d=0where x^=E(x∣Y)=∫xpdf(x∣Y)dx.λ L and x~ are available(24)λL=DDT-1Dx^-dx~=x^-DTDDT-1Dx^-dAdd Steps 4 to PU-EKF measurement update.Steps 4. x^k∣k=x^k∣k-DTDDT-1(Dx^k∣k-d)We get the complete measurement update steps for the linear constraint PUEKF (LPUEKF).


## 4. Simulation Results and Analysis

### 4.1. A Demonstration Example

First, consider a demonstration example of tracking a point-current dipole consisting of two point currents. Consider the simulation scene shown in Table1. Two electric field sensors are located at o1,2. The initial value and initial mean square error of the filter are shown in Table 2. In addition, the process acceleration noise intensities σx=0.1, σy=0.1, and α=0.0001; observation noise covariance Rk=σ2I2×2, σ=1e-5; progressive step of LPUEKF δλ=0.05; seawater depth H=60m; simulation time TN=41. Choose the new algorithm LPUEKF and traditional Kalman filters (EKF, UKF, and CKF) as the tracking algorithms.Table 1
Simulation scene parameters.
Parameter/(unit) Magnitude r 0 / ( m ) - 30 , - 30,12.5 V 0 / ( m / s ) 10,10 I 0 / ( A ) - 5.05,5.05 l d / ( m ) 10 o 1,2 [ 0 , ± 5,0 ]Table 2
Filter initial condition. The initial error is set to maximum to test the performance of the algorithm under large initial error conditions.
Parameter/(unit) Initial value Initial mean square error r ^ 0 / ( m ) [ x 0 , y 0 , 12.51 ] d i a g { p 0 2 , p 0 2 , 0.0 1 2 } V ^ 0 / ( m / s ) 10.2,10.2 d i a g { 0 . 2 2 , 0 . 2 2 } I ^ p / ( A ) 0 1 × 2 d i a g { 6 2 , 6 2 } ∗ p 0 = 5,10,20,40   ( x 0 , y 0 ) = ( - 35 , - 25 ) , ( - 40 , - 20 ) , ( - 50 , - 10 ) , ( - 70,10 ).The root mean square error ofr position component at k moment (RMSEkr) is chosen as performance metric, and it is formulated as(25)RMSEkr=1MC∑i=1MCr^ki-ri22;k=1,⋯,TNFigure2 shows RMSEkr under different initial error conditions.Figure 2
R M S E k r under different initial error conditions.
(a)
p 0 = 5 (b)
p 0 = 10 (c)
p 0 = 20 (d)
p 0 = 40As shown in Figure2,(1)
LPUEKF obtains better results under different initial error conditions. In contrast, performance of traditional algorithms become poor as the initial error increases.(2)
when the target is far away from electric sensors,RMSEkr increases with the decrease of SNR (signal to noise ratio), which is due to the fact that the dynamic estimation depends on the current observation information.


### 4.2. A Simulation with Beasy Data

#### 4.2.1. Beasy Software Simulation of Ship UEP

Beasy software is used to evaluate the ship UEP during the ship design phase [25]. Consider a ship: the length 140m, the width 20m, and the draft depth 5m. The hull material is steel and the hull has a corrosion-resistant coating. The propeller material is copper-nickel-aluminum alloy. Assuming that the hull coating is somewhere damaged and the propeller is exposed in the sea, so propeller and hull damage area is corroded area. Seawater conductivity σsea=4s/m, seabed conductivity σsbed=0.1s/m, seawater depth H=60m, and computed plane depth h=23.2m (relative to the sea level). Ship UEP data obtained from Beasy software is approximately consistent with the real ship UEP. Use these data as sensor observation information.


#### 4.2.2. Simulation Conditions, Results, and Analysis

Consider the simulation scene shown in Figure3. Three sensors array lay in line as Y-axis. O is in the middle sensor. Point current array (the ship) is moving parallel to the X-axis. Z-axis is perpendicular to XOY and points to the sea surface. We also use RMSEkr as the performance metric. At the same time, the length and width of the ship have great diversity, so the estimation error of x and y components should also be concerned. The root mean square errors of x,y component at k moment (RMSEkx, RMSEky) are the same as formula (25) in form.Figure 3
Projected diagram of ship passing through the electric field sensor array.The number of point current is 9, and the distanceld=16m. The center point of the point current array is taken as the reference point of tracking. The number of sensor is 3; distance between sensors is 20m; Rk=σ2I3×3, σ=1×10-4; progressive step of LPUEKF is 0.02; seawater depth H=60m; simulation time TN=232; MC=100. The other parameters are shown in Tables 3 and 4.Table 3
Simulation scene parameters.
Parameter/(unit) Magnitude r 0 / ( m ) - 140,10,23.2 V 0 / ( m / s ) 15,0Table 4
Filter initial condition.
Parameter/(unit) Initial value Initial mean square error r ^ 0 / ( m ) - 420,90,23 d i a g { 28 0 2 , 8 0 2 , 0 . 2 2 } V ^ 0 / ( m / s ) 1 7 , 0.26 d i a g { 3 2 , 1 2 } I ^ p / ( A ) 0 1 × 9 I 9 × 9The results are shown in Figure4. UKF allows integral weights to be negative, and in the high dimensional tracking simulation, P is easy to be nonpositive, resulting in filtering anomalies. Therefore, UKF tracking result is not shown in the figure.Figure 4
Tracking results of each algorithm.As shown in Figure4, LPUEKF has better performance than the other algorithms. To further observe the LPUEKF tracking effect, the simulation results are displayed as shown in Figure 5. Figure 5(a) shows RMSEkpos for r, x, and y; Figure 5(b) shows one tracking trajectory of a simulation result.Figure 5
LPUEKF tracking results.
(a)
R M S E k p o sof r,x,y (b)
One of LPUEKF tracking trajectoryIn Figure5(a), at moment k=141, RMSEkr=30.09m, RMSEkx=30.09m, and RMSEky=0.4668m. Compared with the ship scale in direction x(140m), y(20m), and the initial error(280m,60m), the tracking accuracy of the algorithm is enough to meet the demand. We can see from Figure 5(b) that LPUEKF can converge to the correct target position quickly after limited iterations.The ship is described by 9 point currents whose initial value are01×9 and the center point of the point current array is taken as tracking reference point. As a result, local optimal solution inevitably exists in the filtering process. That is why the estimation of the x direction cannot be as accurate as the estimate in the y direction.


## 4.1. A Demonstration Example

First, consider a demonstration example of tracking a point-current dipole consisting of two point currents. Consider the simulation scene shown in Table1. Two electric field sensors are located at o1,2. The initial value and initial mean square error of the filter are shown in Table 2. In addition, the process acceleration noise intensities σx=0.1, σy=0.1, and α=0.0001; observation noise covariance Rk=σ2I2×2, σ=1e-5; progressive step of LPUEKF δλ=0.05; seawater depth H=60m; simulation time TN=41. Choose the new algorithm LPUEKF and traditional Kalman filters (EKF, UKF, and CKF) as the tracking algorithms.Table 1
Simulation scene parameters.
Parameter/(unit) Magnitude r 0 / ( m ) - 30 , - 30,12.5 V 0 / ( m / s ) 10,10 I 0 / ( A ) - 5.05,5.05 l d / ( m ) 10 o 1,2 [ 0 , ± 5,0 ]Table 2
Filter initial condition. The initial error is set to maximum to test the performance of the algorithm under large initial error conditions.
Parameter/(unit) Initial value Initial mean square error r ^ 0 / ( m ) [ x 0 , y 0 , 12.51 ] d i a g { p 0 2 , p 0 2 , 0.0 1 2 } V ^ 0 / ( m / s ) 10.2,10.2 d i a g { 0 . 2 2 , 0 . 2 2 } I ^ p / ( A ) 0 1 × 2 d i a g { 6 2 , 6 2 } ∗ p 0 = 5,10,20,40   ( x 0 , y 0 ) = ( - 35 , - 25 ) , ( - 40 , - 20 ) , ( - 50 , - 10 ) , ( - 70,10 ).The root mean square error ofr position component at k moment (RMSEkr) is chosen as performance metric, and it is formulated as(25)RMSEkr=1MC∑i=1MCr^ki-ri22;k=1,⋯,TNFigure2 shows RMSEkr under different initial error conditions.Figure 2
R M S E k r under different initial error conditions.
(a)
p 0 = 5 (b)
p 0 = 10 (c)
p 0 = 20 (d)
p 0 = 40As shown in Figure2,(1)
LPUEKF obtains better results under different initial error conditions. In contrast, performance of traditional algorithms become poor as the initial error increases.(2)
when the target is far away from electric sensors,RMSEkr increases with the decrease of SNR (signal to noise ratio), which is due to the fact that the dynamic estimation depends on the current observation information.


## 4.2. A Simulation with Beasy Data

### 4.2.1. Beasy Software Simulation of Ship UEP

Beasy software is used to evaluate the ship UEP during the ship design phase [25]. Consider a ship: the length 140m, the width 20m, and the draft depth 5m. The hull material is steel and the hull has a corrosion-resistant coating. The propeller material is copper-nickel-aluminum alloy. Assuming that the hull coating is somewhere damaged and the propeller is exposed in the sea, so propeller and hull damage area is corroded area. Seawater conductivity σsea=4s/m, seabed conductivity σsbed=0.1s/m, seawater depth H=60m, and computed plane depth h=23.2m (relative to the sea level). Ship UEP data obtained from Beasy software is approximately consistent with the real ship UEP. Use these data as sensor observation information.


### 4.2.2. Simulation Conditions, Results, and Analysis

Consider the simulation scene shown in Figure3. Three sensors array lay in line as Y-axis. O is in the middle sensor. Point current array (the ship) is moving parallel to the X-axis. Z-axis is perpendicular to XOY and points to the sea surface. We also use RMSEkr as the performance metric. At the same time, the length and width of the ship have great diversity, so the estimation error of x and y components should also be concerned. The root mean square errors of x,y component at k moment (RMSEkx, RMSEky) are the same as formula (25) in form.Figure 3
Projected diagram of ship passing through the electric field sensor array.The number of point current is 9, and the distanceld=16m. The center point of the point current array is taken as the reference point of tracking. The number of sensor is 3; distance between sensors is 20m; Rk=σ2I3×3, σ=1×10-4; progressive step of LPUEKF is 0.02; seawater depth H=60m; simulation time TN=232; MC=100. The other parameters are shown in Tables 3 and 4.Table 3
Simulation scene parameters.
Parameter/(unit) Magnitude r 0 / ( m ) - 140,10,23.2 V 0 / ( m / s ) 15,0Table 4
Filter initial condition.
Parameter/(unit) Initial value Initial mean square error r ^ 0 / ( m ) - 420,90,23 d i a g { 28 0 2 , 8 0 2 , 0 . 2 2 } V ^ 0 / ( m / s ) 1 7 , 0.26 d i a g { 3 2 , 1 2 } I ^ p / ( A ) 0 1 × 9 I 9 × 9The results are shown in Figure4. UKF allows integral weights to be negative, and in the high dimensional tracking simulation, P is easy to be nonpositive, resulting in filtering anomalies. Therefore, UKF tracking result is not shown in the figure.Figure 4
Tracking results of each algorithm.As shown in Figure4, LPUEKF has better performance than the other algorithms. To further observe the LPUEKF tracking effect, the simulation results are displayed as shown in Figure 5. Figure 5(a) shows RMSEkpos for r, x, and y; Figure 5(b) shows one tracking trajectory of a simulation result.Figure 5
LPUEKF tracking results.
(a)
R M S E k p o sof r,x,y (b)
One of LPUEKF tracking trajectoryIn Figure5(a), at moment k=141, RMSEkr=30.09m, RMSEkx=30.09m, and RMSEky=0.4668m. Compared with the ship scale in direction x(140m), y(20m), and the initial error(280m,60m), the tracking accuracy of the algorithm is enough to meet the demand. We can see from Figure 5(b) that LPUEKF can converge to the correct target position quickly after limited iterations.The ship is described by 9 point currents whose initial value are01×9 and the center point of the point current array is taken as tracking reference point. As a result, local optimal solution inevitably exists in the filtering process. That is why the estimation of the x direction cannot be as accurate as the estimate in the y direction.


## 4.2.1. Beasy Software Simulation of Ship UEP

Beasy software is used to evaluate the ship UEP during the ship design phase [25]. Consider a ship: the length 140m, the width 20m, and the draft depth 5m. The hull material is steel and the hull has a corrosion-resistant coating. The propeller material is copper-nickel-aluminum alloy. Assuming that the hull coating is somewhere damaged and the propeller is exposed in the sea, so propeller and hull damage area is corroded area. Seawater conductivity σsea=4s/m, seabed conductivity σsbed=0.1s/m, seawater depth H=60m, and computed plane depth h=23.2m (relative to the sea level). Ship UEP data obtained from Beasy software is approximately consistent with the real ship UEP. Use these data as sensor observation information.


## 4.2.2. Simulation Conditions, Results, and Analysis

Consider the simulation scene shown in Figure3. Three sensors array lay in line as Y-axis. O is in the middle sensor. Point current array (the ship) is moving parallel to the X-axis. Z-axis is perpendicular to XOY and points to the sea surface. We also use RMSEkr as the performance metric. At the same time, the length and width of the ship have great diversity, so the estimation error of x and y components should also be concerned. The root mean square errors of x,y component at k moment (RMSEkx, RMSEky) are the same as formula (25) in form.Figure 3
Projected diagram of ship passing through the electric field sensor array.The number of point current is 9, and the distanceld=16m. The center point of the point current array is taken as the reference point of tracking. The number of sensor is 3; distance between sensors is 20m; Rk=σ2I3×3, σ=1×10-4; progressive step of LPUEKF is 0.02; seawater depth H=60m; simulation time TN=232; MC=100. The other parameters are shown in Tables 3 and 4.Table 3
Simulation scene parameters.
Parameter/(unit) Magnitude r 0 / ( m ) - 140,10,23.2 V 0 / ( m / s ) 15,0Table 4
Filter initial condition.
Parameter/(unit) Initial value Initial mean square error r ^ 0 / ( m ) - 420,90,23 d i a g { 28 0 2 , 8 0 2 , 0 . 2 2 } V ^ 0 / ( m / s ) 1 7 , 0.26 d i a g { 3 2 , 1 2 } I ^ p / ( A ) 0 1 × 9 I 9 × 9The results are shown in Figure4. UKF allows integral weights to be negative, and in the high dimensional tracking simulation, P is easy to be nonpositive, resulting in filtering anomalies. Therefore, UKF tracking result is not shown in the figure.Figure 4
Tracking results of each algorithm.As shown in Figure4, LPUEKF has better performance than the other algorithms. To further observe the LPUEKF tracking effect, the simulation results are displayed as shown in Figure 5. Figure 5(a) shows RMSEkpos for r, x, and y; Figure 5(b) shows one tracking trajectory of a simulation result.Figure 5
LPUEKF tracking results.
(a)
R M S E k p o sof r,x,y (b)
One of LPUEKF tracking trajectoryIn Figure5(a), at moment k=141, RMSEkr=30.09m, RMSEkx=30.09m, and RMSEky=0.4668m. Compared with the ship scale in direction x(140m), y(20m), and the initial error(280m,60m), the tracking accuracy of the algorithm is enough to meet the demand. We can see from Figure 5(b) that LPUEKF can converge to the correct target position quickly after limited iterations.The ship is described by 9 point currents whose initial value are01×9 and the center point of the point current array is taken as tracking reference point. As a result, local optimal solution inevitably exists in the filtering process. That is why the estimation of the x direction cannot be as accurate as the estimate in the y direction.


## 5. Conclusions

UEP is an important feature of a ship, which contains the location information of the ship. We study the possibility of using UEP to track the ship in this paper. In the current tracking algorithm, the filtering algorithm is the most suitable for real-time tracking. Among them, the Kalman filters are more suitable for underwater applications because of their easy programming and less computation. Aiming at the limitations of traditional Kalman filters in high-dimensional tracking problems, we propose progressive update Kalman filters based on progressive Bayesian. Taking PUEKF as a representative, the simulation shows that it has better convergence than traditional Kalman filters. The simulation with Beasy data demonstrates the good performance of the new algorithm and the possibility of using UEP to track a ship.



---

*Source: 2797621-2018-12-12.xml*