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 , indicates the n-th element in . Let define the transpose of a vector . Let denote the vector . For example, defines the vector composed of first 3 elements in . Considering a matrix , let define the i-th diagonal element of . Let denote a diagonal matrix, whose diagonal elements are in this order. Let define a matrix generated by stacking a matrix over another matrix .
Let N denote the total number of receivers. Let denote the i-th receiver (). Let present the 3D position of the i-th receiver (). Let denote the 3D velocity of the i-th receiver (). In addition, let present the transmitter’s 3D position. Let present the transmitter’s 3D velocity.
The range between the
i-th receiver and the transmitter is:
In addition, the range rate is:
Assume that
, the maximum speed of the transmitter, is known a priori. This implies that:
In the case where one tracks an underwater transmitter, one can set as the maximum speed of an underwater vehicle.
2.1. TDOA Formulation
This subsection introduces the TDOA formulation used in this article. Let
denote the time of arrival (TOA), when a signal is measured to arrive at
. In practice, the measurement noise exists in time arrival measurements. In other words, one utilizes:
Here,
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),
has a Gaussian distribution with mean 0 and standard deviation
. 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,
in Equation (
4) cannot be measured by
.
Let
define the TDOA between
and
. Utilizing Equation (
4), one can measure:
where
and
. Since
has a Gaussian distribution with mean 0 and standard deviation
,
has a Gaussian distribution with mean 0 and standard deviation
.
In Equation (
5),
defines the TDOA formulation measured by both
and
. This time difference of arrival measurement,
, 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
and
is:
2.2. FDOA Formulation
This subsection introduces the FDOA formulation used in this article. The sound of the transmitter is measured by , and 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
is derived as follows. Let
denote the unknown emission frequency, which can change as time goes on. Let
denote the frequency measurement of
, given as:
where
defines the frequency measurement noise, having a Gaussian distribution with mean 0 and standard deviation
.
Equation (
7) leads to:
where
defines the FDOA formulation measured by both
and
. In addition,
has a Gaussian distribution with mean 0 and standard deviation
.
2.3. Least-Squares Estimation
In this subsection, one assumes that
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),
indicates a simulation result of [
1]. We show that the location accuracy of the proposed estimation approach is comparable to
.
Assuming that
in Equation (
7) is known a priori, Equation (
8) results in:
Since is accessible, is available, as long as is accessible.
In addition, Equation (
1) leads to:
Equations (
10) and (
11) result in:
Let
denote the transmitter state to be estimated. Let
be defined as:
In addition, let
be defined as:
Then, Equation (
13) leads to:
As one performs the time derivative of Equation (
12), one obtains:
In addition, let
be defined as:
Then, Equation (
17) leads to:
Utilizing Equations (
20) and (
16), one obtains:
where
, and
.
Then, one obtains the LSE solution as:
Note that in the solution
,
and
are correlated utilizing Equation (
1). In addition,
and
are correlated utilizing Equation (
2).
In order to avoid the case where
is ill-defined, we require that
in Equation (
22) has full ranks.
has size
. For having full ranks, we require that
, i.e.,
. This implies that we require at least five receivers.
3. Hybrid TDOA-FDOA Estimation as Is Not Known in Advance
This section presents how to tackle hybrid TDOA-FDOA localization, considering the case where the emission frequency
in Equation (
7) is not 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 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
in Equation (
7) can exist. Utilizing Equation (
7), the lower bound for
is:
Utilizing Equation (
7), the upper bound for
is:
The emission frequency range is
. By discretizing the emission frequency range, one calculates a set of feasible emission frequencies. One calculates the set as
. Here,
denotes the number of elements in this set. In addition,
. Let
denote the
j-th element in the set. One utilizes:
Utilizing the LSE approach in
Section 2.3,
leads to a TDOA-FDOA solution in Equation (
22). Instead of Equation (
9), one utilizes:
Here, the
operator is used to indicate that
is an estimation of an arbitrary variable
a. In Equation (
26),
is used, since the true emission frequency is not known in advance.
Utilizing
in Equation (
26), Equation (
22) leads to a transmitter solution:
Note that in
,
and
are correlated utilizing Equation (
1). In addition,
and
are correlated utilizing Equation (
2). The residual associated with
is:
where
and
Among all solutions
where
, one searches for a solution that leads to the minimum residual. In other words, one finds
, where
Then, 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
m/s. In Equation (
4),
has a Gaussian distribution with mean 0 and standard deviation
. This implies that the standard deviation for range error is 0.1 m. In Equation (
7),
has a Gaussian distribution, with mean 0 and standard deviation
.
The sampling interval is seconds. As the maximum speed of a transmitter, one utilizes meters per second. Moreover, samples are used in the feasible frequency range.
One runs Monte Carlo simulations. is the transmitter’s 3D location estimate at time-step k under the c-th Monte Carlo simulation. Let represent the true transmitter’s location at time-step k.
Considering the transmitter’s position, the following RMSE (in meters) is used:
Let
denote the transmitter’s 3D velocity estimate at time-step
k under the
c-th Monte Carlo simulation. Let
represent the true transmitter’s velocity at time-step
k. Considering the transmitter’s velocity, the following RMSE (in m/s) is used:
As time goes on, the emission frequency
changes, as presented in
Figure 1. Note that
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
denote the algebraic solution in [
1]. Note that
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
denote the transmitter’s 3D position. Let
define the transmitter’s 3D velocity. Let
denote the estimate of
at time-step
k. Here, recall that
denotes the emission frequency that needs to be estimated. Let
define the error covariance (uncertainty) of the state vector
.
At every time-step
k, one utilizes the extended Kalman filter (EKF) in [
29] to update
. 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:
where
indicates the process noise.
One ignores the effect of process noise in the PCRLB. Then, the error covariance (uncertainty) of the state vector
is predicted as
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:
Utilizing Equations (
5) and (
8),
indicates the error covariance in hybrid TDOA-FDOA measurements.
denotes the error covariance associated with
, where
. In addition,
denotes the error covariance associated with
, where
.
Utilizing Equations (
35)–(
37), the recursive form of the PCRLB is:
Here,
is the Jacobian matrix of
in Equation (
36) given by
.
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:
Here, and denote the initial position uncertainty and the initial frequency uncertainty respectively. In simulations, we utilize .
defines the
j-th diagonal element in
.
indicates the lower bound for the estimation error of
. Since
denotes the estimate of
at time-step
k, the PCRLB for transmitter’s position at time-step
k is set as:
Since
denotes the estimate of
at time-step
k, the PCRLB for transmitter’s velocity at time-step
k is set as:
4.2. Scenario 1
Scenario 1 is as follows. One utilizes
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
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
degree per second. The transmitter moves utilizing the following coordinate turn:
where
Here, , and .
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 (
) as
k varies. The emission frequency
changes, as presented in
Figure 1. In
Figure 3,
represents the case where the emission frequency is not known in advance. Under
, the proposed location scheme in
Section 3 is used for deriving the transmitter solution. In
Figure 3,
indicates the case where the emission frequency is known a priori. Under
, the LSE solution in
Section 2.3 is used for deriving the solution. In addition,
defines the algebraic solution in [
1]. In
Figure 3,
indicates the case where the emission frequency is set as
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 (
in Equation (
33)) as
k varies. The emission frequency
changes, as presented in
Figure 1.
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
in
Section 3 runs as fast as
.
4.3. Scenario 2
Scenario 2 is as follows. One utilizes 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 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 (
in Equation (
32)) as
k changes. The emission frequency
changes, as presented in
Figure 1. In
Figure 6,
indicates the case where the emission frequency is not known in advance. In
Figure 6,
represents the ideal case, where the emission frequency is known a priori.
denotes the algebraic solution in [
1].
indicates the case where the emission frequency is set as
in Equation (
8).
shows worse performance than other solutions, since it utilizes a wrong emitting frequency. The worst performance of
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 (
in Equation (
33)) as
k varies. The emission frequency
changes, as presented in
Figure 1.
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
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
changes, as presented in
Figure 8. See that one has A large variation in
as time goes on.
Regarding Scenario 2 in
Figure 5,
Figure 9 illustrates the RMSE for the transmitter’s position (
in Equation (
32)) as
k changes. As plotted in
Figure 8, one has a large variation in
as time goes on. The worst performance of
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 (
in Equation (
33)) as
k varies. As plotted in
Figure 8, one has a large variation in
as time goes on.
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
Hz at all times. Regarding Scenario 2 in
Figure 5,
Figure 11 illustrates the RMSE for the transmitter’s position (
in Equation (
32)) as
k changes. In
Figure 11, one has
Hz at all times. The worst performance of
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 (
in Equation (
33)) as
k varies. In
Figure 12, one has
Hz at all times.
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
. By discretizing the emission frequency range, one calculates a set of feasible emission frequencies, given as
. Here,
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
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.