Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
On the Stability of Rubble Mound Structures under Oblique Wave Attack
Previous Article in Journal
A Deep Learning Method for Ship Detection and Traffic Monitoring in an Offshore Wind Farm Area
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Underwater Transmitter Localization Based on TDOA and FDOA Considering the Unknown Time-Varying Emission Frequency

System Engineering Department, Sejong University, Seoul 05006, Republic of Korea
J. Mar. Sci. Eng. 2023, 11(7), 1260; https://doi.org/10.3390/jmse11071260
Submission received: 24 May 2023 / Revised: 14 June 2023 / Accepted: 20 June 2023 / Published: 21 June 2023
(This article belongs to the Section Ocean Engineering)

Abstract

:
This article considers locating a noncooperative underwater transmitter utilizing multiple receivers, such that each receiver can measure the frequency difference of arrival (FDOA) as well as the time difference of arrival (TDOA) of the transmitter’s sound. This article considers the case where the unknown emission frequency of the transmitter changes as time goes. This article addresses hybrid TDOA-FDOA localization, under the assumption that the transmitter’s maximum speed is known in advance. To the best of our knowledge, this article is unique in tackling hybrid TDOA-FDOA localization, considering the case where the unknown emission frequency changes as time goes on. Under MATLAB simulations, this article shows that the proposed hybrid localization method is comparable to the ideal case, where the time-varying emission frequency is known in advance. Furthermore, we show that the proposed localization approach outperforms the case where the emission frequency is estimated as a wrong value.

1. Introduction

In underwater environments, electromagnetic signals are easily dissipated; hence, sound is mainly analyzed for underwater transmitter localization. Underwater transmitter localization utilizing sonar measurements is important in passive transmitter tracking.
This article considers locating a noncooperative underwater transmitter utilizing multiple receivers. This article assumes that each receiver can measure the frequency difference of arrival (FDOA) as well as the time difference of arrival (TDOA) of the transmitter’s signal.
To locate a stationary transmitter, one common technique is to measure the TDOAs of the transmitter signal at multiple receivers. Each TDOA defines a hyperbola in which the transmitter can exist. The intersection of the hyperbolas gives the transmitter location estimate. When the transmitter is moving, FDOA measurements can be used in addition to TDOAs to accurately estimate the transmitter’s position and velocity [1].
This article handles the case where the emission frequency of the transmitter is not known in advance. Moreover, the emission frequency can change as time goes on [2]. This is feasible, since one considers tracking a noncooperative transmitter, which can generate a time-varying emission frequency. In practice, a transmitter cannot move with infinite speed. A transmitter’s information can be used to estimate the maximum speed of the transmitter. For instance, in the case where one tracks an underwater transmitter, one has a priori information on the maximum speed of an underwater vehicle. Thus, this article addresses hybrid TDOA-FDOA localization under the assumption that the transmitter’s maximum speed is known in advance.
Various papers [3,4,5] have tackled tracking a moving transmitter based on TDOA-only measurements of multiple sensors. Many papers [6,7,8,9,10] have handled transmitter localization utilizing hybrid TDOA-AOA measurements. Reference [6] considered tracking a transmitter, in the case where a receiver can measure the TDOA, as well as the angle of arrival (AOA) of the transmitter’s sound. Hybrid TDOA-AOA requires only two sensors to locate a transmitter [7,8]. Reference [10] addressed a TDOA-AOA hybrid location estimation utilizing only two sensors, whose bearing noise is sufficiently small. Assuming that the AOA measurement noise is sufficiently small, weighted least-squares estimation (WLSE) can be applied for TDOA-AOA localization [9].
There are many papers on transmitter localization utilizing hybrid TDOA-FDOA measurements [1,11,12,13,14,15,16]. The authors of [11] proposed an pseudolinear method for locating a moving transmitter utilizing TDOA, FDOA, and differential Doppler rate measurements. Their method was based on the pseudolinear equations and the two-step weighted least-squares estimator. Reference [1] addressed an algebraic solution for moving transmitter localization utilizing TDOA and FDOA measurements. Reference [13] addressed the transmitter localization by utilizing TDOA-FDOA measurements, by reformulating the localization problem as a weighted least-squares (WLS) problem and performing semidefinite relaxation (SDR) to obtain a convex semidefinite programming (SDP) problem. Reference [14] developed a robust SDR method for transmitter localization utilizing TDOA-FDOA measurements. Considering TDOA-FDOA measurements, [17] developed a localization method based on an improved invasive weed optimization algorithm for locating a transmitter. Reference [15] addressed a transmitter localization approach based on hybrid TDOA-FDOA, which copes with unknown propagation speed and sensor parameter errors.
As far as we know, no paper on hybrid TDOA-FDOA localization considered the case where the emission frequency is not known in advance. Moreover, no paper on hybrid TDOA-FDOA localization has considered the case where the emission frequency changes as time goes on. However, considering a noncooperative transmitter, its emission frequency may not be available in practice. Moreover, emission frequency can vary as time goes on. This inspired us to handle hybrid TDOA-FDOA localization, considering the case where the unknown emission frequency changes.
In practice, a transmitter vehicle cannot move with infinite speed. This article thus addresses hybrid TDOA-FDOA localization under the assumption that the transmitter’s maximum speed is known in advance. Considering the transmitter’s maximum speed, one derives the emission frequency range where the emission frequency can exist. By discretizing the emission frequency range, one obtains a set of feasible emission frequencies. Then, each feasible emission frequency leads to a TDOA-FDOA least-squares solution. Among all solutions, one searches for a solution that leads to the minimum residual. In this way, one can calculate the most feasible emission frequency as well as the transmitter solution.
To the best of our knowledge, this article is novel in tackling hybrid TDOA-FDOA localization, considering the case where the unknown emission frequency changes. Under MATLAB simulations, this manuscript verifies that the proposed hybrid localization method is comparable to the ideal case, where the varying emission frequency is known in advance. We show that the proposed localization approach outperforms the case where the emission frequency is estimated as a wrong value.
Section 2 states the problem handled in this work. Section 3 provides the TDOA-FDOA estimation in the case where the emission frequency is not known in advance. Section 4 discusses MATLAB simulations. Section 5 provides Conclusions.

2. Problem Formulation

This section discusses the problem handled in this article. Beforehand, one introduces several concepts. Considering a vector A , A [ n ] indicates the n-th element in A . Let A T define the transpose of a vector A . Let A [ n : m ] denote the vector ( A [ n ] , A [ n + 1 ] , , A [ m ] ) T . For example, A [ 1 : 3 ] defines the vector composed of first 3 elements in A . Considering a matrix M , let M [ i , i ] define the i-th diagonal element of M . Let diag ( a , b , c , ) denote a diagonal matrix, whose diagonal elements are a , b , c , in this order. Let ( A ; B ) define a matrix generated by stacking a matrix A over another matrix B .
Let N denote the total number of receivers. Let s i denote the i-th receiver ( i { 1 , 2 , , N } ). Let s i = ( x i , y i , z i ) T present the 3D position of the i-th receiver ( i { 1 , 2 , , N } ). Let s ˙ i = ( x ˙ i , y ˙ i , z ˙ i ) T denote the 3D velocity of the i-th receiver ( i { 1 , 2 , , N } ). In addition, let e = ( x e , y e , z e ) T present the transmitter’s 3D position. Let e ˙ = ( x ˙ e , y ˙ e , z ˙ e ) T present the transmitter’s 3D velocity.
The range between the i-th receiver and the transmitter is:
r i = e s i .
In addition, the range rate is:
r ˙ i = ( e ˙ s ˙ i ) T ( e s i ) r i .
Assume that V m a x , the maximum speed of the transmitter, is known a priori. This implies that:
V m a x x ˙ e 2 + y ˙ e 2 + z ˙ e 2 .
In the case where one tracks an underwater transmitter, one can set V m a x as the maximum speed of an underwater vehicle.

2.1. TDOA Formulation

This subsection introduces the TDOA formulation used in this article. Let t i denote the time of arrival (TOA), when a signal is measured to arrive at s i . In practice, the measurement noise exists in time arrival measurements. In other words, one utilizes:
t i = e s i C + n i t o .
Here, n i t o indicates time measurement noise, and C indicates the signal speed. It is assumed that C is known a priori. This assumption is widely applied in TOA-based target tracking [18,19,20].
In Equation (4), n i t o has a Gaussian distribution with mean 0 and standard deviation σ t o . This Gaussian noise has been widely applied in TOA-based localization [21,22,23,24].
This article considers a noncooperative transmitter whose signal emission time is not known in advance. Therefore, t i in Equation (4) cannot be measured by s i .
Let t i , 1 = t i t 1 define the TDOA between s i and s 1 . Utilizing Equation (4), one can measure:
t i , 1 = e s i C e s 1 C + n i , 1 t o ,
where t i , 1 = t i t 1 and n i , 1 t o = n i t o n 1 t o . Since n i t o has a Gaussian distribution with mean 0 and standard deviation σ t o , n i , 1 t o has a Gaussian distribution with mean 0 and standard deviation 2 σ t o .
In Equation (5), t i , 1 defines the TDOA formulation measured by both s 1 and s i . This time difference of arrival measurement, t i , 1 , is applied in TDOA localization problems [25,26]. The TDOA measurement constructs a hyperbolic curve, indicating feasible transmitter locations [6].
Utilizing Equation (5), the transmitter range difference between s 1 and s i is:
r i , 1 = C × ( t i , 1 n i , 1 t o ) = r i r 1 .

2.2. FDOA Formulation

This subsection introduces the FDOA formulation used in this article. The sound of the transmitter is measured by s i , and s i can analyze the frequency of the sound utilizing the fast Fourier transform (FFT).
To the best of our knowledge, no paper on hybrid TDOA-FDOA localization [1,11,12,13,14,15,16] has considered the case where the emission frequency is not known in advance. As presented in [27], the frequency measurement of s i is derived as follows. Let f e denote the unknown emission frequency, which can change as time goes on. Let f i denote the frequency measurement of s i , given as:
f i = f e × ( 1 r ˙ i C ) + n i f ,
where n i f defines the frequency measurement noise, having a Gaussian distribution with mean 0 and standard deviation σ f .
Equation (7) leads to:
f i , 1 = f i f 1 = f e ( r ˙ 1 r ˙ i C ) + n i , 1 f ,
where f i , 1 defines the FDOA formulation measured by both s 1 and s i . In addition, n i , 1 f = n i f n 1 f has a Gaussian distribution with mean 0 and standard deviation 2 σ f .

2.3. Least-Squares Estimation

In this subsection, one assumes that f e in Equation (7) is known in advance. One addresses a least-squares estimation (LSE) to calculate the transmitter solution based on TDOA and FDOA measurements.
In practical underwater environments, the noise covariance of measurements may not be known accurately. Thus, the LSE solution ignores the noise terms in the measurements. Moreover, the LSE solution runs fast, since it does not consider the noise covariance.
We acknowledge that one can utilize the noise covariance for deriving the transmitter solution. Considering the ideal case, where the emission frequency is known in advance, the authors of [1] addressed an algebraic solution for a transmitter’s location utilizing hybrid TDOA-FDOA measurements. The algebraic solution in [1] used the noise covariance explicitly. For comparison with the proposed hybrid tracking method, we ran the algebraic solution in [1]. In MATLAB simulation section (Section 4), a l g e b indicates a simulation result of [1]. We show that the location accuracy of the proposed estimation approach is comparable to a l g e b .
Assuming that f e in Equation (7) is known a priori, Equation (8) results in:
r ˙ i , 1 = r ˙ i r ˙ 1 = C × f i , 1 f e .
Since f i , 1 is accessible, r ˙ i , 1 is available, as long as f e is accessible.
Equation (6) leads to:
( r i , 1 + r 1 ) 2 = r i 2 .
In addition, Equation (1) leads to:
r i 2 = ( e s i ) T × ( e s i ) .
Equations (10) and (11) result in:
r i , 1 2 + s 1 T × s 1 s i T × s i + 2 ( s i s 1 ) T × e + 2 × r i , 1 × r 1 = 0 .
This leads to:
0.5 × ( r i , 1 2 + s 1 T × s 1 s i T × s i ) = ( s i T + s 1 T ) × e r i , 1 × r 1 .
Let S = ( e , r 1 , e ˙ , r ˙ 1 ) T denote the transmitter state to be estimated. Let h t be defined as:
h t = 0.5 × ( r 2 , 1 2 + s 1 T × s 1 s 2 T × s 2 ) 0.5 × ( r 3 , 1 2 + s 1 T × s 1 s 3 T × s 3 ) 0.5 × ( r N , 1 2 + s 1 T × s 1 s N T × s N ) .
In addition, let G t be defined as:
G t = s 2 T + s 1 T r 2 , 1 0 0 s 3 T + s 1 T r 3 , 1 0 0 s N T + s 1 T r N , 1 0 0 .
Then, Equation (13) leads to:
h t = G t × S .
As one performs the time derivative of Equation (12), one obtains:
r i , 1 × r ˙ i , 1 + s ˙ 1 T × s 1 s ˙ i T × s i = ( s ˙ 1 T s ˙ i T ) × e r ˙ i , 1 × r 1 + ( s 1 s i ) T × e ˙ r i , 1 × r ˙ 1 .
Let h f be defined as:
h f = r 2 , 1 × r ˙ 2 , 1 + s ˙ 1 T × s 1 s ˙ 2 T × s 2 r 3 , 1 × r ˙ 3 , 1 + s ˙ 1 T × s 1 s ˙ 3 T × s 3 r N , 1 × r ˙ N , 1 + s ˙ 1 T × s 1 s ˙ N T × s N .
In addition, let G f be defined as:
G f = ( s ˙ 1 T s ˙ 2 T ) r ˙ 2 , 1 ( s 1 s 2 ) T r 2 , 1 ( s ˙ 1 T s ˙ 3 T ) r ˙ 3 , 1 ( s 1 s 3 ) T r 3 , 1 ( s ˙ 1 T s ˙ N T ) r ˙ N , 1 ( s 1 s N ) T r N , 1 .
Then, Equation (17) leads to:
h f = G f × S .
Utilizing Equations (20) and (16), one obtains:
h = G × S .
where h = ( h t ; h f ) , and G = ( G t ; G f ) .
Then, one obtains the LSE solution as:
S = ( G T G ) 1 G T h .
Note that in the solution S = ( e , r 1 , e ˙ , r ˙ 1 ) T , r 1 and e are correlated utilizing Equation (1). In addition, r ˙ 1 and e ˙ are correlated utilizing Equation (2).
In order to avoid the case where S is ill-defined, we require that G in Equation (22) has full ranks. G has size ( 2 N 2 ) × 8 . For having full ranks, we require that 2 N 2 8 , i.e., N 5 . This implies that we require at least five receivers.

3. Hybrid TDOA-FDOA Estimation as f e Is Not Known in Advance

This section presents how to tackle hybrid TDOA-FDOA localization, considering the case where the emission frequency f e in Equation (7) is not known in advance. Considering the transmitter’s maximum speed V m a x , one derives the emission frequency range where the emission frequency can exist. By discretizing the emission frequency range, one calculates a set of feasible emission frequencies. Then, each feasible emission frequency leads to a TDOA-FDOA least-squares solution in Equation (22). Among all solutions, one searches for a solution that leads to the minimum residual. In this way, one can calculate the most feasible emission frequency as well as the transmitter solution.
One derives the range where f e in Equation (7) can exist. Utilizing Equation (7), the lower bound for f e is:
f m i n e = m i n i { 1 , 2 , , N } f i × C C + V m a x .
Utilizing Equation (7), the upper bound for f e is:
f m a x e = m a x i { 1 , 2 , , N } f i × C C V m a x .
The emission frequency range is [ f m i n e , f m a x e ] . By discretizing the emission frequency range, one calculates a set of feasible emission frequencies. One calculates the set as [ f m i n e , f m i n e + Δ Q , f m i n e + 2 Δ Q , , f m a x e ] . Here, Q + 1 denotes the number of elements in this set. In addition, Δ = f m a x e f m i n e . Let f ^ e ( j ) denote the j-th element in the set. One utilizes:
f ^ e ( j ) = f m i n e + ( j 1 ) × Δ Q .
Utilizing the LSE approach in Section 2.3, f ^ e ( j ) leads to a TDOA-FDOA solution in Equation (22). Instead of Equation (9), one utilizes:
r ˙ ^ i , 1 ( j ) = C × f i , 1 f ^ e ( j ) .
Here, the a ^ operator is used to indicate that a ^ is an estimation of an arbitrary variable a. In Equation (26), f ^ e ( j ) is used, since the true emission frequency is not known in advance.
Utilizing r ˙ ^ i , 1 ( j ) in Equation (26), Equation (22) leads to a transmitter solution:
S ^ ( j ) = ( G T G ) 1 G T h .
Note that in S , r 1 = S [ 4 ] and e = S [ 1 : 3 ] are correlated utilizing Equation (1). In addition, r ˙ 1 = S [ 8 ] and e ˙ = S [ 5 : 7 ] are correlated utilizing Equation (2). The residual associated with f ^ e ( j ) is:
r e s ( j ) = r e s 1 ( j ) + r e s 2 ( j ) ,
where
r e s 1 ( j ) = ( S ^ ( j ) [ 1 : 3 ] s 1 S ^ ( j ) [ 4 ] ) 2
and
r e s 2 ( j ) = ( S ^ ( j ) [ 5 : 7 ] s ˙ 1 S ^ ( j ) [ 8 ] ) 2 .
Among all solutions S ^ ( j ) where j { 1 , 2 , , N } , one searches for a solution that leads to the minimum residual. In other words, one finds S ^ ( j * ) , where
j * = a r g m i n j { 1 , 2 , , N } r e s ( j ) .
Then, S ^ ( j * ) is set as the transmitter solution, in the case where the emission frequency is not available.

4. MATLAB Simulations

MATLAB simulations are utilized to show the effectiveness of the proposed localization scheme. The scenario runs for 500 s.
In MATLAB simulations, one considers tracking an underwater transmitter utilizing sound measurements. One utilizes the sound speed as C = 1400 m/s. In Equation (4), n i t o has a Gaussian distribution with mean 0 and standard deviation σ i t o = 0.1 C . This implies that the standard deviation for range error is 0.1 m. In Equation (7), n i f has a Gaussian distribution, with mean 0 and standard deviation σ f = 10 C .
The sampling interval is T = 5 seconds. As the maximum speed of a transmitter, one utilizes V m a x = 50 meters per second. Moreover, Q = 100 samples are used in the feasible frequency range.
One runs M c = 100 Monte Carlo simulations. L k c ( c { 1 , 2 , , M c } ) is the transmitter’s 3D location estimate at time-step k under the c-th Monte Carlo simulation. Let L k represent the true transmitter’s location at time-step k.
Considering the transmitter’s position, the following RMSE (in meters) is used:
R M S E k = c = 1 M c L k c L k 2 M c .
Let V k c ( c { 1 , 2 , , M c } ) denote the transmitter’s 3D velocity estimate at time-step k under the c-th Monte Carlo simulation. Let V k represent the true transmitter’s velocity at time-step k. Considering the transmitter’s velocity, the following RMSE (in m/s) is used:
R M S E k V = c = 1 M c V k c V k 2 M c .
As time goes on, the emission frequency f e changes, as presented in Figure 1. Note that f e is not known in advance; thus, it needs to be estimated under the proposed tracking method.
Assuming that the emission frequency is known in advance, the authors of [1] addressed an algebraic solution for a transmitter’s location, utilizing hybrid TDOA-FDOA measurements. For comparison with the proposed hybrid tracking method, one runs the algebraic solution in [1]. For convenience, let a l g e b denote the algebraic solution in [1]. Note that a l g e b considers the ideal case, where the emission frequency is known a priori. MATLAB simulations will show that the proposed hybrid localization method in Section 3 is comparable to the ideal case, where the emission frequency is known in advance.

4.1. Posterior Cramer–Rao Lower Bound (PCRLB)

This article utilizes the posterior Cramer–Rao lower bound (PCRLB) as the lowest estimation error of any unbiased estimator [28]. Let e = ( x e , y e , z e ) T denote the transmitter’s 3D position. Let e ˙ = ( x ˙ e , y ˙ e , z ˙ e ) T define the transmitter’s 3D velocity. Let X k denote the estimate of ( e ; e ˙ ; f e ) at time-step k. Here, recall that f e denotes the emission frequency that needs to be estimated. Let B k | k define the error covariance (uncertainty) of the state vector X k .
At every time-step k, one utilizes the extended Kalman filter (EKF) in [29] to update B k | k . According to [29], one obtains the recursive form of the PCRLB by merging the prediction and the update step for the error covariance of the EKF.
The transmitter moves utilizing the following motion model:
X k + 1 = X k + n k ,
where n k indicates the process noise.
One ignores the effect of process noise in the PCRLB. Then, the error covariance (uncertainty) of the state vector X k + 1 is predicted as
B k + 1 | k = B k | k .
This implies that the PCRLB considers the case where there is no uncertainty in the motion of the transmitter.
Utilizing Equations (5) and (8), the hybrid TDOA-FDOA measurement equation is:
M ( X k ) = ( t 2 , 1 ; t 3 , 1 ; t N , 1 ; f 2 , 1 ; f 3 , 1 ; f N , 1 ) .
Utilizing Equations (5) and (8),
R = diag ( 2 σ t o 2 , 2 σ t o 2 , , 2 σ f 2 , 2 σ f 2 )
indicates the error covariance in hybrid TDOA-FDOA measurements. 2 σ t o 2 denotes the error covariance associated with t i , 1 , where i N . In addition, 2 σ f 2 denotes the error covariance associated with f i , 1 , where i N .
Utilizing Equations (35)–(37), the recursive form of the PCRLB is:
( B k + 1 | k + 1 ) 1 = ( B k | k ) 1 + M J T ( R ) 1 M J .
Here, M J is the Jacobian matrix of M ( X k ) in Equation (36) given by M J = M ( X ) X | X = X k .
One utilizes the PCRLB in Equation (38) as the lowest estimation error of any unbiased estimator [28]. In Equation (38), the initial error covariance utilizes:
B 0 | 0 = diag ( δ 2 , δ 2 , δ 2 , V m a x 2 , V m a x 2 , V m a x 2 , ( δ f ) 2 ) .
Here, δ > 0 and δ f > 0 denote the initial position uncertainty and the initial frequency uncertainty respectively. In simulations, we utilize δ = δ f = 1 .
B k | k [ j , j ] defines the j-th diagonal element in B k | k . B k | k [ j , j ] indicates the lower bound for the estimation error of X k [ j ] . Since X k denotes the estimate of [ e ; e ˙ ; f e ] at time-step k, the PCRLB for transmitter’s position at time-step k is set as:
P C R L B k = j = 1 3 B k | k [ j , j ] .
Since X k denotes the estimate of [ e ; e ˙ ; f e ] at time-step k, the PCRLB for transmitter’s velocity at time-step k is set as:
P C R L B k V = j = 4 6 B k | k [ j , j ] .

4.2. Scenario 1

Scenario 1 is as follows. One utilizes N = 6 receivers in total. The receivers are located at [100, 0, 100], [−100, 1500, 1], [200, 1000, 10], [100, 3000, 1], [−5100, −1000, 10], and [5100, 5000, 1] in meters, respectively. The initial transmitter’s location is [ 1000 , 1000 , 200 ] in meters. Initially, the transmitter’s velocity is [10, 0, 0] in m/s. The transmitter performs a coordinates turn in the xy-plane with turn rate ω = 0.5 degree per second. The transmitter moves utilizing the following coordinate turn:
X k + 1 = F c X k ,
where
F c = 1 0 0 d t × S w d t × C w 0 0 1 0 d t × C w d t × S w 0 0 0 1 0 0 d t 0 0 0 cos ( ω × d t ) sin ( ω × d t ) 0 0 0 0 sin ( ω × d t ) cos ( ω × d t ) 0 0 0 0 0 0 1 .
Here, S w = sin ( ω × d t ) ω × d t , and C w = 1 cos ( ω × d t ) ω × d t .
Figure 2 illustrates Scenario 1. In Figure 2, the transmitter’s path is depicted with red circles. Deployed receivers are marked with black diamonds.
Considering Scenario 1 in Figure 2, Figure 3 illustrates the RMSE for transmitter’s position ( R M S E k ) as k varies. The emission frequency f e changes, as presented in Figure 1. In Figure 3, u n k n o w n E m i t F r e q represents the case where the emission frequency is not known in advance. Under u n k n o w n E m i t F r e q , the proposed location scheme in Section 3 is used for deriving the transmitter solution. In Figure 3, k n o w n E m i t F r e q indicates the case where the emission frequency is known a priori. Under k n o w n E m i t F r e q , the LSE solution in Section 2.3 is used for deriving the solution. In addition, a l g e b defines the algebraic solution in [1]. In Figure 3, w r o n g E m i t F r e q indicates the case where the emission frequency is set as f 2 , 1 in Equation (8). In this case, the LSE solution in Section 2.3 is used to derive the solution. Note that the proposed hybrid localization method in Section 3 is comparable to the ideal case, where the emission frequency is known in advance.
Regarding Scenario 1 in Figure 2, Figure 4 illustrates the RMSE for the transmitter’s velocity ( R M S E k V in Equation (33)) as k varies. The emission frequency f e changes, as presented in Figure 1. w r o n g E m i t F r e q does not appear in this plot, since its RMSE is too large. Note that the proposed hybrid localization method in Section 3 is comparable to the ideal case, where the emission frequency is known in advance.
Considering Scenario 1, Table 1 presents the computational time for all MC simulations, considering every algorithm in Figure 4. Note that the proposed location scheme u n k n o w n E m i t F r e q in Section 3 runs as fast as k n o w n E m i t F r e q .

4.3. Scenario 2

Scenario 2 is as follows. One utilizes N = 7 receivers in total. Every receiver is located at [100, 0, 100], [−100, 1500, 1], [200, 1000, 10], [100, 3000, 1], [−5100, −1000, 10], [5100, 5000, 1], and [510, 0, 10], respectively.
The initial transmitter’s location is [ 1000 , 1000 , 200 ] in meters. The transmitter moves with a constant velocity [10, 0, 0] in m/s.
Figure 2 illustrates Scenario 2. In Figure 2, the transmitter’s path is marked with red circles. Deployed receivers are marked with black diamonds.
Regarding Scenario 2 in Figure 5, Figure 6 illustrates the RMSE for the transmitter’s position ( R M S E k in Equation (32)) as k changes. The emission frequency f e changes, as presented in Figure 1. In Figure 6, u n k n o w n E m i t F r e q indicates the case where the emission frequency is not known in advance. In Figure 6, k n o w n E m i t F r e q represents the ideal case, where the emission frequency is known a priori. a l g e b denotes the algebraic solution in [1]. w r o n g E m i t F r e q indicates the case where the emission frequency is set as f 2 , 1 in Equation (8). w r o n g E m i t F r e q shows worse performance than other solutions, since it utilizes a wrong emitting frequency. The worst performance of w r o n g E m i t F r e q indicates that the unknown emission frequency can significantly degrade the localization performance. Figure 6 shows that the proposed hybrid localization method in Section 3 is comparable to the ideal case, where the emission frequency is known in advance.
Considering Scenario 2, Figure 7 illustrates the RMSE for the transmitter’s velocity ( R M S E k V in Equation (33)) as k varies. The emission frequency f e changes, as presented in Figure 1. w r o n g E m i t F r e q does not appear in this plot, since its RMSE is too large.
Considering Scenario 2, Table 2 presents the computational time for all MC simulations, considering every algorithm in Figure 7. Acknowledge that the proposed location scheme u n k n o w n E m i t F r e q in Section 3 runs slightly slower than other algorithms. In the Conclusions, we address a parallel processing approach to improve the time efficiency of the proposed location scheme.

4.4. The Effect of Changing the Emission Frequency

We check the effect of changing the emission frequency. As time goes on, the emission frequency f e changes, as presented in Figure 8. See that one has A large variation in f e as time goes on.
Regarding Scenario 2 in Figure 5, Figure 9 illustrates the RMSE for the transmitter’s position ( R M S E k in Equation (32)) as k changes. As plotted in Figure 8, one has a large variation in f e as time goes on. The worst performance of w r o n g E m i t F r e q indicates that the unknown emission frequency can significantly degrade the localization performance. Note that the proposed hybrid localization method in Section 3 is comparable to the ideal case, where the emission frequency is known in advance.
Regarding Scenario 2, Figure 10 illustrates the RMSE for the transmitter’s velocity ( R M S E k V in Equation (33)) as k varies. As plotted in Figure 8, one has a large variation in f e as time goes on. w r o n g E m i t F r e q does not appear in this plot, since its RMSE is too large. Note that the proposed hybrid localization method in Section 3 is comparable to the ideal case, where the emission frequency is known in advance.
One next considers the case where the emission frequency does not change as time goes on. One has f e = 1000 Hz at all times. Regarding Scenario 2 in Figure 5, Figure 11 illustrates the RMSE for the transmitter’s position ( R M S E k in Equation (32)) as k changes. In Figure 11, one has f e = 1000 Hz at all times. The worst performance of w r o n g E m i t F r e q indicates that the unknown emission frequency can significantly degrade the localization performance. Note that the proposed hybrid localization method in Section 3 is comparable to the ideal case, where the emission frequency is known in advance.
Regarding Scenario 2, Figure 12 illustrates the RMSE for the transmitter’s velocity ( R M S E k V in Equation (33)) as k varies. In Figure 12, one has f e = 1000 Hz at all times. w r o n g E m i t F r e q does not appear in this plot, since its RMSE is too large.

5. Conclusions

This article tackles locating an underwater transmitter utilizing hybrid TDOA-FDOA measurements, in the case where the unknown emission frequency varies as time goes on. Utilizing MATLAB simulations, this article demonstrates the superiority of the proposed location scheme in Section 3, by comparing it with the ideal case, where the emission frequency is known a priori. Furthermore, MATLAB simulations show that the proposed localization approach outperforms the case where the emission frequency is estimated as a wrong value. The proposed location scheme can be used to locate a ground or aerial target that emits RF signals. In the future, experiments will be performed utilizing real sonar sensors, in order to demonstrate the effectiveness of the proposed hybrid localization more rigorously.
We can further improve the computational efficiency of the proposed localization methods using parallel processing. Here, parallel processing is a computing technique where multiple streams of data processing tasks co-occur through multiple central processing units (CPUs) working concurrently. Recall that in Section 3, the emission frequency range is [ f m i n e , f m a x e ] . By discretizing the emission frequency range, one calculates a set of feasible emission frequencies, given as [ f m i n e , f m i n e + Δ Q , f m i n e + 2 Δ Q , , f m a x e ] . Here, Q + 1 denotes the number of feasible emission frequencies. In the case where we can use multiple CPUs, we can apply parallel processing to improve the time efficiency. In the case where we have M CPUs, each CPU can process Q + 1 M frequencies in a parallel manner. In this way, we can improve the time efficiency of the proposed localization methods.
In practice, receiver position and velocity errors can occur [30]. For instance, in underwater environments, sea current leads to receiver position errors. The uncertain receiver location can seriously deteriorate the localization performance. The authors of [31] addressed a Taylor series method using hybrid TDOA-FDOA to estimate both the transmitter state and receiver position simultaneously. Reference [12] considered locating multiple transmitters under hybrid TDOA-FDOA measurements, in the presence of sensor position and velocity errors. In the future, we will handle hybrid TDOA-FDOA target tracking considering both receiver position and velocity errors.

Funding

This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIT) (Grant Number: 2022R1A2C1091682). This research was supported by the faculty research fund of Sejong University in 2023.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ho, K.C.; Xu, W. An accurate algebraic solution for moving source location using TDOA and FDOA measurements. IEEE Trans. Signal Process. 2004, 52, 2453–2463. [Google Scholar] [CrossRef]
  2. Witzgall, H. Ground vehicle doppler geolocation. In Proceedings of the 2014 IEEE Aerospace Conference, Big Sky, MT, USA, 1–8 March 2014; pp. 1–8. [Google Scholar]
  3. Kim, J. Fast non-line-of-sight receivers conjecturing method in TDOA localisation using obstacle information. IET Radar Sonar Navig. 2019, 13, 347–351. [Google Scholar] [CrossRef]
  4. Wu, P.; Guo, Q.; Zhang, X.; Bo, Y. Maneuvering target tracking using passive TDOA measurements. In Proceedings of the 33rd Chinese Control Conference, Nanjing, China, 28–30 July 2014; pp. 758–762. [Google Scholar] [CrossRef]
  5. Chen, L.; Li, S. IMM tracking of a 3D maneuvering target with passive TDOA system. In Proceedings of the 2003 International Conference on Neural Networks and Signal Processing, Nanjing, China, 14–17 December 2003; Volume 2, pp. 1611–1614. [Google Scholar] [CrossRef]
  6. Beerens, S.P.; van Ijsselmuide, S.; van Koersel, A.C. Passive ranging with flank and towed array sensors. In Proceedings of the 2003 Undersea Defence Technology (UDT) Conference, Malmo, Sweden, 24–26 June 2003. [Google Scholar]
  7. Li, W.; Chen, R.; Guo, Y.; Fu, C. Closed form algorithm of double-satellite TDOA+AOA localization based on WGS-84 model. Chin. J. Aeronaut. 2019, 32, 2354–2367. [Google Scholar] [CrossRef]
  8. Yin, J.; Wan, Q.; Yang, S.; Ho, K.C. A Simple and Accurate TDOA-AOA Localization Method Using Two Stations. IEEE Signal Process. Lett. 2016, 23, 144–148. [Google Scholar] [CrossRef]
  9. Li, W.; Tang, Q.; Huang, C.E.A. A New Close Form Location Algorithm with AOA and TDOA for Mobile User. Wirel. Pers. Commun. 2017, 97, 3061–3080. [Google Scholar] [CrossRef]
  10. Cong, L.; Zhuang, W. Hybrid TDOA/AOA Mobile User Location for Wideband CDMA Cellular Systems. IEEE Trans. Wirel. Commun. 2002, 1, 439–447. [Google Scholar] [CrossRef] [Green Version]
  11. Liu, Z.; Hu, D.; Zhao, Y.; Zhao, Y. An algebraic method for moving source localization using TDOA, FDOA, and differential Doppler rate measurements with receiver location errors. EURASIP J. Adv. Signal Process. 2019, 2019, 25. [Google Scholar] [CrossRef] [Green Version]
  12. Sun, M.; Ho, K.C. An Asymptotically Efficient Estimator for TDOA and FDOA Positioning of Multiple Disjoint Sources in the Presence of Sensor Location Uncertainties. IEEE Trans. Signal Process. 2011, 59, 3434–3440. [Google Scholar]
  13. Wang, G.; Li, Y.; Ansari, N. A Semidefinite Relaxation Method for Source Localization Using TDOA and FDOA Measurements. IEEE Trans. Veh. Technol. 2013, 62, 853–862. [Google Scholar] [CrossRef]
  14. Wang, Y.; Wu, Y. An Efficient Semidefinite Relaxation Algorithm for Moving Source Localization Using TDOA and FDOA Measurements. IEEE Commun. Lett. 2017, 21, 80–83. [Google Scholar] [CrossRef]
  15. Zhang, B.; Hu, Y.; Wang, H.; Zhuang, Z. Underwater Source Localization Using TDOA and FDOA Measurements With Unknown Propagation Speed and Sensor Parameter Errors. IEEE Access 2018, 6, 36645–36661. [Google Scholar] [CrossRef]
  16. Musicki, D.; Kaune, R.; Koch, W. Mobile Emitter Geolocation and Tracking Using TDOA and FDOA Measurements. IEEE Trans. Signal Process. 2009, 58, 1863–1874. [Google Scholar] [CrossRef]
  17. Lin, Y.; Zhang, Z.; Najafabadi, H.E. Underwater source localization using time difference of arrival and frequency difference of arrival measurements based on an improved invasive weed optimization algorithm. IET Signal Process. 2022, 16, 299–309. [Google Scholar] [CrossRef]
  18. Su, J.; Li, Y.; Ali, W.; Li, X.; Yu, J. Underwater 3D Doppler-Angle Target Tracking with Signal Time Delay. Sensors 2020, 20, 3869. [Google Scholar] [CrossRef] [PubMed]
  19. Wang, X.; Musicki, D.; Ellem, R.; Fletcher, F. Enhanced Multi-Target Tracking with Doppler Measurements. In Proceedings of the 2007 Information, Decision and Control Conference, Adelaide, SA, Australia, 12–14 February 2007; pp. 53–58. [Google Scholar]
  20. Kim, J. Three dimensional tracking of a maneuvering emitter utilizing doppler-bearing measurements of a constant velocity observer. Signal Process. 2021, 189, 108246. [Google Scholar] [CrossRef]
  21. Park, C.H.; Chang, J.H. Robust TOA source localisation algorithm using priori location. IET Radar Sonar Navig. 2019, 13, 384–391. [Google Scholar] [CrossRef]
  22. Laaraiedh, M.; Yu, L.; Avrillon, S.; Uguen, B. Comparison of Hybrid Localization Schemes using RSSI, TOA, and TDOA. In Proceedings of the 17th European Wireless 2011—Sustainable Wireless Technologies, Vienna, Austria, 27–29 April 2011; pp. 1–5. [Google Scholar]
  23. Wann, C.; Chin, H. Hybrid TOA/RSSI Wireless Location with Unconstrained Nonlinear Optimization for Indoor UWB Channels. In Proceedings of the 2007 IEEE Wireless Communications and Networking Conference, Hong Kong, China, 11–15 March 2007; pp. 3940–3945. [Google Scholar]
  24. Guvenc, I.; Chong, C.C. A Survey on TOA Based Wireless Localization and NLOS Mitigation Techniques. IEEE Commun. Surv. Tutor. 2009, 11, 107–124. [Google Scholar] [CrossRef]
  25. Torrieri, D.J. Statistical Theory of Passive Location Systems. IEEE Trans. Aerosp. Electron. Syst. 1984, 20, 183–198. [Google Scholar] [CrossRef] [Green Version]
  26. Bucher, R.; Misra, D. A Synthesizable VHDL Model of the Exact Solution for Three-dimensional Hyperbolic Positioning System. VLSI Des. 2002, 15, 507–520. [Google Scholar] [CrossRef] [Green Version]
  27. Kim, J. 3D path planner of an autonomous underwater vehicle to track an emitter using frequency and azimuth–elevation angle measurements. IET Radar Sonar Navig. 2020, 14, 1236–1243. [Google Scholar] [CrossRef]
  28. Logothetis, A.; Isaksson, A.; Evans, R.J. An information theoretic approach to observer path design for bearings-only tracking. In Proceedings of the 36th IEEE Conference on Decision and Control, San Diego, CA, USA, 12 December 1997; Volume 4, pp. 3132–3137. [Google Scholar] [CrossRef]
  29. Ristic, B.; Arulampalam, S.; Gordon, N. Beyond the Kalman Filter: Particle Filters for Tracking Applications; Artech House Radar Library Series; Artech: London, UK, 2004. [Google Scholar]
  30. Deng, B.; Sun, Z.B.; Peng, H.F.; Xiong, J.Y. Source localization using TDOA/FDOA/DFS measurements with erroneous sensor positions. In Proceedings of the 2016 CIE International Conference on Radar (RADAR), Guangzhou, China, 10–13 October 2016; pp. 1–4. [Google Scholar]
  31. Hao, W.; Wei-Min, S.; Hong, G. A novel Taylor series method for source and receiver localization using TDOA and FDOA measurements with uncertain receiver positions. In Proceedings of the 2011 IEEE CIE International Conference on Radar, Chengdu, China, 24–27 October 2011; Volume 2, pp. 1037–1040. [Google Scholar]
Figure 1. The change of f e as time goes on.
Figure 1. The change of f e as time goes on.
Jmse 11 01260 g001
Figure 2. Scenario 1. The initial transmitter’s location is [ 1000 , 1000 , 200 ] in meters. The transmitter’s path is marked with red circles. Deployed receivers are depicted with black diamonds.
Figure 2. Scenario 1. The initial transmitter’s location is [ 1000 , 1000 , 200 ] in meters. The transmitter’s path is marked with red circles. Deployed receivers are depicted with black diamonds.
Jmse 11 01260 g002
Figure 3. Scenario 1. RMSE (Equation (32)) plot for transmitter’s position. The emission frequency f e changes, as presented in Figure 1. u n k n o w n E m i t F r e q indicates the case where the emission frequency is not known a priori. k n o w n E m i t F r e q represents the ideal case, where the emission frequency is known in advance. In addition, a l g e b defines the algebraic solution in [1]. w r o n g E m i t F r e q indicates the case where the emission frequency is set as f 2 , 1 in Equation (8).
Figure 3. Scenario 1. RMSE (Equation (32)) plot for transmitter’s position. The emission frequency f e changes, as presented in Figure 1. u n k n o w n E m i t F r e q indicates the case where the emission frequency is not known a priori. k n o w n E m i t F r e q represents the ideal case, where the emission frequency is known in advance. In addition, a l g e b defines the algebraic solution in [1]. w r o n g E m i t F r e q indicates the case where the emission frequency is set as f 2 , 1 in Equation (8).
Jmse 11 01260 g003
Figure 4. Scenario 1. RMSE (Equation (33)) plot for transmitter’s velocity. The emission frequency f e changes, as presented in Figure 1. w r o n g E m i t F r e q does not appear in this plot, since its RMSE is too large.
Figure 4. Scenario 1. RMSE (Equation (33)) plot for transmitter’s velocity. The emission frequency f e changes, as presented in Figure 1. w r o n g E m i t F r e q does not appear in this plot, since its RMSE is too large.
Jmse 11 01260 g004
Figure 5. Scenario 2. The initial transmitter’s location is [ 1000 , 1000 , 200 ] in meters. The transmitter moves with a constant velocity [10, 0, 0] in m/s. The transmitter’s path is depicted with red circles. Deployed receivers are marked with black diamonds.
Figure 5. Scenario 2. The initial transmitter’s location is [ 1000 , 1000 , 200 ] in meters. The transmitter moves with a constant velocity [10, 0, 0] in m/s. The transmitter’s path is depicted with red circles. Deployed receivers are marked with black diamonds.
Jmse 11 01260 g005
Figure 6. Scenario 2. RMSE (Equation (32)) plot for transmitter’s position. The emission frequency f e changes, as presented in Figure 1. u n k n o w n E m i t F r e q indicates the case where the emission frequency is not known a priori. w r o n g E m i t F r e q indicates the case where the emission frequency is set as f 2 , 1 in Equation (8).
Figure 6. Scenario 2. RMSE (Equation (32)) plot for transmitter’s position. The emission frequency f e changes, as presented in Figure 1. u n k n o w n E m i t F r e q indicates the case where the emission frequency is not known a priori. w r o n g E m i t F r e q indicates the case where the emission frequency is set as f 2 , 1 in Equation (8).
Jmse 11 01260 g006
Figure 7. Scenario 2. RMSE (Equation (33)) plot for transmitter’s velocity. The emission frequency f e changes, as presented in Figure 1. w r o n g E m i t F r e q does not appear in this plot, since its RMSE is too large.
Figure 7. Scenario 2. RMSE (Equation (33)) plot for transmitter’s velocity. The emission frequency f e changes, as presented in Figure 1. w r o n g E m i t F r e q does not appear in this plot, since its RMSE is too large.
Jmse 11 01260 g007
Figure 8. Large change of f e as time goes on.
Figure 8. Large change of f e as time goes on.
Jmse 11 01260 g008
Figure 9. Scenario 2. As plotted in Figure 8, one has large variation in f e as time goes on. RMSE (Equation (32)) plot for transmitter’s position.
Figure 9. Scenario 2. As plotted in Figure 8, one has large variation in f e as time goes on. RMSE (Equation (32)) plot for transmitter’s position.
Jmse 11 01260 g009
Figure 10. Scenario 2. As plotted in Figure 8, one has large variation in f e as time goes on. w r o n g E m i t F r e q does not appear in this plot, since its RMSE is too large. RMSE (Equation (33)) plot for transmitter’s velocity.
Figure 10. Scenario 2. As plotted in Figure 8, one has large variation in f e as time goes on. w r o n g E m i t F r e q does not appear in this plot, since its RMSE is too large. RMSE (Equation (33)) plot for transmitter’s velocity.
Jmse 11 01260 g010
Figure 11. Scenario 2. One has f e = 1000 Hz at all times. RMSE (32) plot for transmitter’s position.
Figure 11. Scenario 2. One has f e = 1000 Hz at all times. RMSE (32) plot for transmitter’s position.
Jmse 11 01260 g011
Figure 12. Scenario 2. One has f e = 1000 Hz at all times. w r o n g E m i t F r e q does not appear in this plot, since its RMSE is too large. RMSE (Equation (33)) plot for transmitter’s velocity.
Figure 12. Scenario 2. One has f e = 1000 Hz at all times. w r o n g E m i t F r e q does not appear in this plot, since its RMSE is too large. RMSE (Equation (33)) plot for transmitter’s velocity.
Jmse 11 01260 g012
Table 1. Computation time for all MC simulations (Scenario 1).
Table 1. Computation time for all MC simulations (Scenario 1).
Alg.Comp. Time (s)
w r o n g E m i t F r e q 134
k n o w n E m i t F r e q 133
a l g e b 130
u n k n o w n E m i t F r e q 124
Table 2. Computation time for all MC simulations (Scenario 2).
Table 2. Computation time for all MC simulations (Scenario 2).
Alg.Comp. Time (s)
w r o n g E m i t F r e q 130
k n o w n E m i t F r e q 138
a l g e b 141
u n k n o w n E m i t F r e q 161
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

Kim, J. Underwater Transmitter Localization Based on TDOA and FDOA Considering the Unknown Time-Varying Emission Frequency. J. Mar. Sci. Eng. 2023, 11, 1260. https://doi.org/10.3390/jmse11071260

AMA Style

Kim J. Underwater Transmitter Localization Based on TDOA and FDOA Considering the Unknown Time-Varying Emission Frequency. Journal of Marine Science and Engineering. 2023; 11(7):1260. https://doi.org/10.3390/jmse11071260

Chicago/Turabian Style

Kim, Jonghoek. 2023. "Underwater Transmitter Localization Based on TDOA and FDOA Considering the Unknown Time-Varying Emission Frequency" Journal of Marine Science and Engineering 11, no. 7: 1260. https://doi.org/10.3390/jmse11071260

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