# A Mathematical Algorithm of Locomotive Source Localization Based on Hyperbolic Technique

Authors: Homayun Kabir; Jeevan Kanesan; Ahmed Wasif Reza; Harikrishnan Ramiah
Journal: International Journal of Distributed Sensor Networks
 (2015)

---

## Abstract

Recent trend shows that sensors situated on an axis in two-dimensional scenario measuring the time difference of arrival (TDOA) and frequency difference of arrival (FDOA) of the emitting signal from a moving source can estimate the emitting signal’s position and velocity from the intersection point of hyperbola, which derives from TDOA and FDOA. However, estimating the location of an emitter based on hyperbolic measurements is a highly nonlinear problem with inconsistent data, which are created due to the measurement noise, the deviation between assumption model and actual field of the velocity, and so forth. In addition, the coefficient matrix of TDOA and FDOA equations set is singular in the linear sensor array network (LSAN). In this paper, a noniterative and simpler method is proposed to locate the instantaneous position of the moving source in LSAN by estimating the position and velocity based on TDOA and FDOA which does not have the convergence problem. In addition, the method avoids the singularity problem of LSAN by introducing the nuisance variables. The proposed method achieved the theoretical lower bound for near to far field with same and different velocity and different baseline of sensors in low to moderate noise.


---

## Body

## 1. Introduction

In recent decades, localization awareness has aroused great interest in wireless sensor network (WSN) applications, such as in tracking the mobile, homeland security, defense command and control, ship navigation, and disaster management [1–6]. The instantaneous position of signal emitter is determined by using relative position information as bearing [7], distance, received signal strength (RSS) [8], received signal strength ratio (RSSR) [9], time of arrival (TOA) [10], and time difference of arrival (TDOA) [11] of sensors which capture the signal from the source in WSN. In this paper, the source position was estimated through a passive localization technique, called TDOA, represented in Figure 1 as this technique does not require timestamping. Furthermore, it is easy to implement in the practical localization set-up [12]. For the moving source, velocity estimation applies the vital rule for instantaneous localization of the source [13]. Hence, frequency difference of arrival (FDOA) was combined with the TDOA for improving the localization accuracy of mobile source [14, 15]. Two types of geometrical shapes of sensors in WSN were used, namely, arbitrary sensor array network (ASAN) and linear sensor array network (LSAN). In ASAN geometry, the two-step weighted least square (LS) method was derived in [14] to estimate the position and velocity of moving source based on TDOA and FDOA. Furthermore, an accurate and closed-form solution was proposed based on multidimensional scaling (MDS) analysis that optimized a cost function related to the MDS’s scalar product matrix in [16]. In addition, a closed-form of LS estimation was derived for circular array sensor network in [17]. An analytical solution for cuboid network without premeasured wave velocity was proposed in [18]. Another LS method in [19] and an effective closed-form solution in [20] were derived to locate the mobile source position by considering the source and reference sensor motions which create ASAN. Most of the source localization algorithms are for estimating the position of stationary sources. A stationary source localization algorithm using LS estimation that utilized the multiple-input multiple-output radar system with widely separated antennas was derived in [6]. Also, a noniterative method with comparative performance based on TDOA for Global Positioning System (GPS) was proposed in [21, 22]. In this paper, we have focused our research based on the LSAN to locate the position and velocity of moving source at near to far field.Figure 1
TDOA based localization technique in LSAN.It has been found from the literature that the TDOA and FDOA equations with high nonlinearity form hyperbolas, which may not be intersected at a single point due to the inconsistent data, such as measurement error and the deviation between assumption model and actual field of the velocity, in LSAN [14, 23]. Also, the singularity problem arises in the coefficient matrix of hyperbolic equations set in 2D LSAN [23, 24]. To overcome the challenges in the 2D scenario, a mathematical solution by approximate ML estimation based on TOA was developed in [25], which was only applicable for three linear sensors. Additionally, the ML estimation in [26], a closed-form of LS estimation in [19, 23], and the geometric solution in [23] were proposed based on the TDOA to predict the position of stationary source in LSAN.Moreover, the trigonometric mathematical approach in [5] was derived to estimate the position of a moving source in LSAN for only a special case where the range between the source and receivers was large compared to the spacing between receivers (far field). To the best of the authors’ knowledge, no literature has been found with the linear sensor array in a 2D scenario for estimating the position as well as velocity of near to far field moving emitter. In the 2D scenario for LSAN, each element in the array lies in same axis [23]. The coordinate of the source position in that axis is also absent in the nonlinear equation set. Hence, coefficient matrix of the LSAN is singular [24]. For overcoming these issues, we have proposed a mathematical approach based on TDOA and FDOA to estimate the position and velocity of a moving source in LSAN. In this paper, we simulated the proposed method to estimate the position and velocity of moving source at near to far field with the same and different velocity in LSAN where Gaussian noise is considered.The paper is arranged in the following manner. In the following section, the proposed mathematical model is presented. Next, the derivation of Cramer-Rao lower bound (CRLB) is provided. After that, the results are analyzed and performance evaluation is explained. Then, the limitation and future work are discussed. Finally, the concluding remarks are provided.


## 2. Proposed Method

TheN linear stationary sensors of LSAN are considered in a 2D space to estimate the locomotive emitter with unknown position p=xyT and the velocity v=vxvyT utilizing the TDOA and the FDOA localization techniques, where matrix transpose operation is denoted by T, represented in Figure 1. The signals emitted from the source is captured by N stationary sensors of LSAN, which are situated at Ri=xiyit where yi=0 and i=1,2,3,…,N. The range between the ith sensor and emitter is(1)di0=p-Ri=p-RiT(p-Ri).The true path difference of sensor pair ith and reference that assumed sensor 1 is (2)di10=cti1=di0-d10.With regard to (2), c is the velocity of signal propagation and ti1 is the time difference of sensor pair ith and reference. Equation (2) is rearranged as di10+d10=di0 and then both sides are squared. Substituting the d10 and di0 from (1), the TDOA equation set obtained is(3)di102+2d10di10=RiTRi-R1TR1-2Ri-R1Tpwherei=2,3,…,N.Equation (3) is a nonlinear set with unknown p and d10, which created N-1 hyperbolic curves with focus Ri=xiyit where i=1,2,3,…,N. Those hyperbolic curves intersect at a point, which gives the estimated position of the source. In the 2D scenario, at least two hyperbolic curves are required to solve the localization problem by utilizing the TDOA, represented in Figure 2.Figure 2
Source position estimation by the hyperbolic technique.Not only the position but also the velocity estimation is essential for determining the instantaneous position of moving emitter. Conversely, TDOA equations set may not be adequate to provide the needed localization accuracy of moving emitter as TDOA estimates only the position of the source. The FDOA measurement which is obtained from the relative velocity between the source [13] and sensors is applied to improve instantaneous localization accuracy of the source. The relation between the path rate and source position parameters is obtained by the time derivative of (1) as follows:(4)vi0=vT(p-Ri)di0.The FDOA is calculated by the time derivative of (3) as follows:(5)2(vi10di10+vi10d10+di10v10)=-2Ri-R1Tvwherei=2,3,…,N.In terms of (5), the path difference rate is denoted by vi10, which is obtained from FDOA. The unknown position p and velocity v of the source are estimated by solving the obtained TDOA and FDOA equations set. The ni1 and nvi1 are additive white Gaussian noise summed with the true range difference and the rate of range difference of sensor pair, respectively, as shown in (6) and (7): (6)di1=di10+ni1,(7)vi1=vi10+nvi1.The vector of noisy TDOA and FDOA, D=d21d31⋯dN1, and V=v21v31⋯vN1 have a covariance matrix [14]. Consider(8)Q=EDTVTTDTVT=σ2Q1O0Q1.Here, σ2 is the variance of zero mean Gaussian noise and 0 is a zero square matrix, and (9)Q1=10.50.5⋯0.50.510.5⋯0.50.50.51⋯0.5⋮⋮⋮⋮⋮0.50.50.5⋯1,whichisa(N-1)×(N-1)matrix.Substituting the di1 and vi1 in (3) and (5), the two nonlinear equation sets obtained are(10)2Ri-R1Tp+2di1d10=-di12+RiTRi-R1TR1,2vi1d10+2Ri-R1Tv+2di1v10=-2di1vi1.In the LSAN, y-axis coordinates of sensors are equal to zero. Hence, the singularity problem arises in (10). For mitigating the singularity problem, these equations can be rewritten as (11)-2xi-x1x-2di1d10=di12+x12-xi2,(12)-2vi1d10-2xi-x1vx-2di1v10=2di1vi1.With regard to (12) and (13), x, d10, vx, and v10 are unknown. The solutions of the above equations are obtained by initially denoting an auxiliary vector θ=xd10vxv10T, which contains the x-axis unknown position parameters and two nuisance variables d10 and v10 (nuisance variables are associated with the variation of the dependent variable as an outcome, which is extraneous to the effects of the independent variables). The reference sensor position is considered at the origin and the error vector of the above equations (11) and (12) is (13)ε=εtεfT=h-gθ,where εt and εf are the error vector of TDOA and FDOA equation sets, respectively, whereas(14)h=d212+x12-x22d312+x12-x32⋮dN12+x12-xN22d21v212d21v31⋮2dN1vN1,g=-2x2-x1d2100x3-x1d3100⋮⋮⋮⋮xN-x1dN1000v21x2-x1d210v31x3-x1d31⋮⋮0vN1xN-x1dN1.To minimize the error vector, (13) can be written as (15)θ=argmin⁡θ(h-gθ).Here, the weight least square method is used for minimizing the error vector. We get the unknown vector as (16)θ=gTWg-1gTWh,where W=Q-1. From the above equation, only x-axis unknown position x and velocity vx of the locomotive source and two nuisance variables d10 and 〈v10〉 are estimated. To determine the y-axis position parameter, 〈x〉, 〈vx〉, 〈d10〉, and v10 values are substituted into (1) and then the equation can be rewritten as (17)y=d102-x2.Finding the appropriate sign of y-axis source position coordinate, the equation cost error function J=ε1TQε1 of (11) has to be minimized. Regarding equation cost error function, ε1=h1-g1θ1, and (18)h1=d212+x12-x22d312+x12-x32⋮dN12+x12-xN2g1=-2x2-x1d21x3-x1d31⋮⋮xN-x1dN1,θ1=xd10T.This solution is identical to the one given in [26]. Finally, the obtained y-axis position parameter value is put into (4) to estimate the y-axis velocity parameter as indicated below: (19)vy=d10v10-xvxy.


## 3. Cramer-Rao Lower Bound

It is important to know the optimum achievable localization accuracy that can be attained with the available measurement set. The CRLB provides a lower bound on the covariance that is asymptotically achievable by any unbiased estimation algorithm [14]. Therefore, the CRLB sets a benchmark of an unbiased estimation, which has been compared with the proposed method. However, the CRLB is equal to the sum of the diagonal elements of covariance matrix estimation. To estimate the covariance matrix of source position and velocity in LSAN based on the TDOA and FDOA, we perturbed the random quantities in θ mentioned as auxiliary vector and proceed as before to obtain the following [23]:(20)EΔθΔθT=g0Tφ-1g0-1.With regard to (20), φ=B1QB1, B1=BOB0B, B=2diagd20d30d40⋯dN0, B0=2diagv20v30v40⋯vN0, and O is a zero square matrix. Following (1), θ1=ΔxΔyT and Δd10 are related to (21)Δd10=Δθ1Tp0-R1d10.According to (4), θ2=ΔvxΔvyT and Δv10 are related to (22)Δv10=Δθ2Tp0-R1d10.Equation (22) is formulated due to ∂d10/∂x=∂v10/∂vx and ∂d10/∂y=∂v10/∂vy. From (21) and (22), we obtain(23)Δθ=KΔθ1Δθ2=1000x0-x1d10x0-x1d1000001000x0-x1d10x0-x1d10ΔxΔyΔvxΔvy.Therefore, the required CRLB is(24)CRLBθ=K-1g0Tφ-1g0-TK-T=KTg0TB1-1Q-1g0TK-1.


## 4. Results and Discussion

Simulation results to estimate the position and velocity based on TDOA and FDOA in LSAN for 2D scenario are presented in this section. The sensors positions are (xi=i and yi=0) and (xi=-(i-1) and yi=0) when i (i=1,2,…,10) is even and odd number, respectively. The near to far field sources are situated at A (8 m, 22 m), B (0 m, 50 m), C (−30 m, 25 m), and D (−50 m, 250 m) [23]. The source velocities are V1 (−2 m/s, 1.5 m/s), V2 (1 m/s, 2 m/s), and V3 (−4 m/s, 2 m/s). In all cases, the proposed method and Taylor’s series [27, 28] are compared against CRLB [23] calculated from (24) and zero mean Gaussian noise is considered as dB=10log⁡⁡(σ2). The MSE of the proposed method is calculated via MSEp=∑1Mp-p02/M and MSEv=∑1Mv-v02/M for position and velocity of the moving source where M=1E5 is the quantity of random generation to maintain the covariance of Gaussian noise, p0 is A or B or C or D, and v0 is V1 or V2 or V3. The deviation percentages of the proposed method and Taylor’s series from the CRLB are measured via ((PMSEp-CMSEp)/CMSEp)×100, ((TMSEp-CMSEp)/CMSEp)×100, ((PMSEv-CMSEv)/CMSEv)×100, and ((TMSEv-CMSEv)/CMSEv)×100 for position and velocity where the symbols P, C, and T represent the proposed, CRLB, and Taylor’s methods, respectively.The MSE of position and velocity estimation of the proposed method against Taylor’s series method and CRLB is represented in Figure3 at the noise levels −100 dB to −5 dB. The position and velocity estimation through the proposed method reaches the CRLB, when the noise is below −30 dB. The proposed method diverges from CRLB with an increase in the noise level starting from −30 dB. The MSEs in both position and velocity of the proposed method are 1.005, 1.21, and 2.85 times higher than the CRLB at noise level −80 dB, −25 dB, and −5 dB, respectively. Below −35 dB noise, the ratio of the proposed method varies between 1.005 and 1.05 against CRLB. It should be noted here that Taylor’s series method is also simulated for comparison with the proposed method in the same simulation environment. This method uses truncated Taylor’s series expansion (avoiding the higher order terms) to linearize the TDOA and FDOA nonlinear equations with iterative solution. The MSEs of the emitter by Taylor’s series depended on the initial guess. Also, they only converges to the local minimum solution. Hence, it gives good position and velocity accuracy when the initial guess is also approximate with actual position and velocity of the emitter at very low level noise. On the other hand, the MSEs are large compared to proposed method. Now, to avoid the initialization limitation of Taylor’s series, 1.2 times of actual values of position and velocity are assumed as the initial values of source position and velocity. In this case, a minimum of 4 iterations is needed for each solution when the noise is less than −40 dB. The deviation of Taylor’s series results from CRLB starts from −40 dB noise. The iteration number and the deviation increased with the increment of the noise level. The MSE of Taylor’s series is higher than the proposed method and CRLB, which is shown in Figure 3 for position and velocity estimation of moving source.Figure 3
Comparison of (a) position and (b) velocity MSE of the proposed method with Taylor’s series and CRLB for near field.
(a)
(b)The MSE comparison of source position and velocity estimation at near field for the proposed method, Taylor’s series, and CRLB at noise −40 dB is represented in Tables1 and 2, respectively, where the number of sensors in LSAN is varied. At 10 sensors in LSAN, the MSE of Taylor’s series is slightly less than the CRLB and proposed method due to the low effective noise level. Tables 1 and 2 show that the MSE at the near field source position and velocity reduced with the increment of a number of sensors. In addition, the deviation rate of Taylor’s series from CRLB is higher than the proposed method at 7 to 3 sensors and 8 to 3 sensors for position and velocity estimation, respectively.Table 1
MSE comparison of source position at near field.
Number of sensors
Proposed method
CRLB
Taylor’s series
(m2)
(m2)
(m2)
10
0.009532
0.0095923
0.0094996
9
0.017638
0.0176473
0.0176473
8
0.028882
0.0288823
0.0289751
7
0.062159
0.0619726
0.062521
6
0.1219
0.1208001
0.1228523
5
0.356023
0.3557193
0.3567235
4
1.121039
1.0944217
1.182211
3
8.29721
7.5660109
8.4193549Table 2
MSE comparison of source velocity at near field.
Number of sensors
Proposed method
CRLB
Taylor’s series
(m2/s2)
(m2/s2)
(m2/s2)
10
0.001085
0.001082
0.001072
9
0.001954
0.001933
0.001975
8
0.00318
0.003146
0.003204
7
0.006872
0.00681
0.00694
6
0.013721
0.013588
0.013842
5
0.039103
0.038908
0.039025
4
0.128801
0.122769
0.133223
3
1.017037
0.923471
1.076407Figure4 shows the obtained MSE estimation of position and velocity at far field using the proposed method, Taylor’s series, and CRLB at noise ranges from −100 dB to −20 dB. The MSEs of the proposed method, Taylor’s series, and CRLB are almost same, that is, below −50 dB noise. After −50 dB, the MSE of the proposed method and Taylor’s series becomes higher with an increment in noise. However, the rising slope of the proposed method (position and velocity estimation of distant sources) is less than Taylor’s series, which is clearly depicted in Figure 4.Figure 4
Comparison of (a) position and (b) velocity MSE of the proposed method with Taylor’s series and CRLB for far field.
(a)
(b)The MSEs of far field source position and velocity estimation differ through variation of the baseline of LSAN, as illustrated in Tables3 and 4. In this simulation, −50 dB noises are considered. The true baseline of LSAN is decreased by reducing the number of sensors from the network. Also, the comparative baseline (ratio of the true baseline and the range between source and sensor network) of LSAN decreases. In addition, direction control becomes weaker and weaker due to the reduction of true baseline. Here, it describes the accuracy of the position and the velocity for a pair of sensors and is mostly limited to one direction, which is parallel to these two sensors and perpendicular to the LSAN [29, 30]. Hence, the large variation of position and velocity is observed. Most interestingly, it can be observed from Tables 3 and 4 that the MSEs of position and velocity are too large when the number of sensors is three due to the too weak direction control. In addition, the MSEs of Taylor’ series are significantly larger than the proposed method when the number of sensors is less (4 or 5) due to its linearization errors [23, 31]. The baseline of the network increases with increasing the number of sensors; as a result, the linearization error reduces. In addition, favorable initial guess is also needed for Taylor’s series. In practice, this is not possible and solution divergence may occur.Table 3
Comparison of MSE source position for the proposed method, CRLB, and Taylor’s series at far field.
Number of sensors
Proposed method
CRLB
Taylor’s series
(m2)
(m2)
(m2)
10
10.856
10.896
10.846
9
18.907
18.909
18.908
8
35.365
35.361
35.371
7
70.232
70.093
70.493
6
157.73
155.62
160.21
5
430.68
409.71
449.71
4
1570.9
1449.89
1813.8Table 4
Comparison of MSE source velocity for the proposed method, CRLB, and Taylor’s series at far field.
Number of sensors
Proposed method
CRLB
Taylor’s series
(m2/s2)
(m2/s2)
(m2/s2)
10
1.1007
1.1006
1.1005
9
1.8332
1.8328
1.8343
8
3.4095
3.4071
3.4195
7
6.7421
6.7679
6.8792
6
15.728
15.359
16.283
5
43.499
40.464
44.937
4
182.59
163.37
203.41The comparative position and velocity MSE of the proposed method and Taylor’s series with respect to theoretical position and velocity MSE are presented in Table5 where different source positions A, B, C, and D and same velocity V1 for all source positions are considered. First of all, the deviation percentages of positions are 2.15% and 6.28% (when noise is −50 dB and the number of sensors is 5 in LSAN) and the deviation percentages of velocity are 2.47% and 7.56% for the proposed and Taylor’s series, respectively, at source position C. In addition, the results of position and velocity obtained for the proposed method are 0.5 and 0.6 times Taylor’s series at position C when noise is −30 dB. Secondly, six and seven numbers of sensors in LSAN at noise −40 dB are also considered to measure the comparative position and velocity MSE in Table 5. Here, the deviation of MSE percentage of Taylor’s series is 2 to 5 times compared to the proposed method. Moreover, the deviation percentage goes higher with the reduction of the number of sensors in LSAN.Table 5
The comparative MSE of the proposed method and Taylor’s series with respect to CRLB at near to far source position with same velocity.
Position
Method (%)
5 numbers of sensors
Noise −40 dB
Noise −50 dB
Noise −30 dB
7 numbers of sensors
6 numbers of sensors
Position
Velocity
Position
Velocity
Position
Velocity
Position
Velocity
A (8 m, 22 m)
Proposed methodCRLB (MSE)Taylor’s series
0.2146
0.2475
1.2784
0.9871
0.4251
0.5101
0.9123
0.9588
0.0356
0.00389
3.5636
0.3885
0.0619
0.0068
0.1208
0.0136
1.1024
1.0125
4.1358
3.5784
1.2146
1.1572
1.7265
1.857
B (0 m, 50 m)
Proposed methodCRLB (MSE)Taylor’s series
1.0001
0.9578
5.1279
6.8795
1.0127
1.9245
1.5723
2.4782
0.5645
0.0589
57.557
5.8953
0.96703
0.0999
2.152
0.2227
3.5789
3.4978
11.2345
13.5789
2.6987
1.8912
2.6987
2.7682
C (−30 m, 25 m)
Proposed methodCRLB (MSE)Taylor’s series
2.1579
2.4798
8.1256
9.5879
2.1987
1.0124
3.2783
3.4324
1.2233
0.1452
139.561
13.906
2.0938
0.2368
5.5342
0.6121
6.2879
7.5679
16.248
15.871
4.5789
6.1234
7.2453
8.7152
D (−50 m, 250 m)
Proposed methodCRLB (MSE)Taylor’s series
5.1183
7.5013
Large
Large
9.9821
11.278
17.563
18.201
409.71
40.464
841.84
61.733
1856.51
135.4
9.763
11.005
16.782
15.721
31.298
29.458The position and velocity MSE’s results of the proposed method and Taylor’s series with respect to the theoretical MSE are observed at −40 dB noise and 6 numbers of sensors in LSAN where combination of different positions and velocity of sources are considered in Table6. For all cases, the proposed method’s results are found to be 2 to 3 times better than Taylor’s series. In addition, the deviation percentage of position MSE is marginally increased through the increment of the source velocity. It is to be noted that the deviation percentage of velocity MSE is higher than the deviation percentage of position MSE.Table 6
The comparative MSE of the proposed method and Taylor’s series with respect to CRLB at near to far source position with different velocity.
Position
Method (%)
Nose −40 dB and  6 numbers  of sensors
V
1 (2 m/s, −1.5 m/s)
V
2 (1 m/s, 2 m/s)
V
3 (−4 m/s, 2 m/s)
Position
Velocity
Position
Velocity
Position
Velocity
A (8 m, 22 m)
Proposed methodCRLB (MSE)Taylor’s series
0.9123
0.9588
0.8826
0.9009
0.932
1.1029
0.1208
0.0136
0.1189
0.0129
0.12247
0.01889
1.7265
1.857
1.638
1.4289
1.987019
2.1748
B (0 m, 50 m)
Proposed methodCRLB (MSE)Taylor’s series
1.5723
2.4782
1.3978
1.4123
1.71423
1.8567
2.152
0.2227
2.1222
0.2191
2.21976
0.2719
2.6987
2.7682
2.6127
2.6212
2.97845
3.0784
C (−30 m, 25 m)
Proposed methodCRLB (MSE)Taylor’s series
3.2783
3.4324
3.2315
3.3012
3.29723
3.9256
5.5342
0.6121
5.4994
0.6023
5.5685
0.9373
7.2453
8.7152
7.2021
8.5278
7.37258
9.0157
D (−50 m, 250 m)
Proposed methodCRLB (MSE)Taylor’s series
17.563
18.201
16.928
13.859
19.28794
17.758
1856.51
135.4
1.85604
135.3
1856.69
136.51
31.298
29.458
30.661
27.972
33.56872
31.257The most interesting observation from Tables5 and 6 is that the MSE at source position C is larger than at position B despite the long distance between the source position B and network. This happens because (1) the source C is situated behind the outside to LSAN where the control of x-direction is extremely weak and (2) the source B is on the perpendicular line of the LSAN. Therefore, the source position and velocity will be undefined when it is on the axis, but outside (either left or right) of the LSAN. On the other hand, the source position and velocity will be effectively estimated when its position is on the perpendicular line of the LSAN [29, 30].In conclusion, it is apparent that the MSE of position and velocity of the far field source is higher than the near field because of the geometric spreading that is the ability to estimate the position and the velocity of emitter becomes weaker and weaker as the position moves away from the sensor network [29, 30]. In our simulation results, the proposed method yields better results than Taylor’s series due to the initialization problem, local minimum solution, linearization errors, and so forth, of Taylor’s series. Therefore, the proposed method in close proximity with the CRLB from near to far field source with same and various velocities and different baseline of network at varying noise levels.


## 5. Conclusion

The nonlinear localization equations set, measurement noise, and singularity problem in LSAN pose the challenges to locate the position and velocity of the locomotive source in the 2D scenario based on TDOA and FDOA measurements. To overcome these challenges, nuisance variables are introduced in this study. These variables have contributed to avoidance of the singularity problem of LSAN in nonlinear localization equations set and to improvement of the instantaneous source location estimation. The proposed method is found to be noniterative, of low complexity, and attractive and does not have convergence problem and initialization problems, as in Taylor’s series. The proposed method accomplished the CRLB for low to moderate noise level in case of moving source, which is positioned at near to far field with same and different velocity under the Gaussian noise.



---

*Source: 384180-2015-10-19.xml*