Next Article in Journal
NE-Motion: Visual Analysis of Stroke Patients Using Motion Sensor Networks
Previous Article in Journal
Numerical Analysis of Protection Method of Metallic Sub-Wavelength Concentric Arrays for Radially Polarized Light Selection and Its Applications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

A Tri-Satellite Interference Source Localization Method for Eliminating Mirrored Location

1
The 54th Research Institute of CETC, Shijiazhuang 050081, China
2
Hebei Key Laboratory of Electromagnetic Spectrum Cognition and Control (ESCC), Shijiazhuang 050081, China
3
School of Electronic Information and Electrical Engineering, Shanghai Jiao Tong University, Shanghai 200240, China
4
National Laboratory of Radar Signal Processing, Xidian University, Xi’an 710071, China
*
Author to whom correspondence should be addressed.
Sensors 2021, 21(13), 4483; https://doi.org/10.3390/s21134483
Submission received: 26 May 2021 / Revised: 16 June 2021 / Accepted: 24 June 2021 / Published: 30 June 2021
(This article belongs to the Section Navigation and Positioning)

Abstract

:
With the increase in satellite communication interference, the tri-satellite time difference of arrival (TDOA) localization technique, which is an effective method to determine the location of the interference using sensors or antennas, has been developed rapidly. The location of the interference source is determined through the intersection of the TDOA lines of position (LOP). However, when the two TDOA LOP have two mirrored intersection points, it is theoretically difficult to determine the real location. Aiming at this problem, a method for eliminating mirrored location based on multiple moment TDOA is proposed in this paper. First, the TDOA results are measured at multiple moments using the cross-ambiguity function (CAF), and the localization equation set is established based on the World Geodetic System (WGS)-84 earth ellipsoid model. Then, the initial location result can be obtained by solving the equation set through the Newton iteration method. Finally, the high-precision location result after eliminating the mirrored location is obtained after the single moment localization based on the initial location. Simulation experiments and real measured data processing results verify the effectiveness of the proposed method. It still has good robustness under the condition of large measurement errors and deviations from the prior initial values.

1. Introduction

With the rapid growth of satellite communication services, the interference signals received are also increasing. Satellite interference source localization technology [1,2], which can receive electromagnetic signals passively to locate the interference source, is the principal method of satellite anti-interference technologies. Aiming at the complex interference environment of satellite communication, several advanced technologies have been proposed to suppress the satellite interference [3,4,5] in recent years. Meanwhile, several satellite interference systems have been developed [6,7], represented by a “NOSS” series of a marine surveillance satellite developed by the US Navy, the Transmitter Location System (TLS) model 2000 developed by the Interferometric Inc., the SatID developed and operated by Defence Research Angency (DRA) of the UK, and the State Space Satellite Monitoring System (SSSMS) developed by the State Radio Regulation of China. The satellite interference source localization technology has the advantages of wide coverage, long operating distance, concealed reception, and strong survivability. Using the passive location, the satellite localization technology has strong survivability and concealed reception. The localization technology is based on the long baseline interferometry method for the satellites to realize wide area coverage. In addition, the detection of very weak signals is achieved through correlation techniques. It can locate interference source targets quickly and accurately, and has an important role and broad application prospects in both civil and military fields.
The current satellite interference source localization methods mainly use the time difference of arrival (TDOA) and the frequency difference of arrival (FDOA) measured by the signals received from the main satellite and the adjacent satellites to determine the location. The localization methods can be classified into analytical methods [8,9,10,11,12,13,14], optimized solutions [15,16,17,18,19,20], and intelligent algorithms [21,22,23] in recent years. The analytical methods derive the closed-form solution of the nonlinear localization equation set based on mathematical approximation and the intermediate variable. The nonlinear equations of localization are transformed into a set of linear equations, in which the closed-form solution can be derived. In the optimization techniques, the localization problem is reformulated as a convex optimization problem and can be solved based on the convex optimization algorithm. However, the local convergence may be produced in the optimization procedure. Considering the nonlinear characteristics of the location equations, the intelligent methods are proposed in recent years to deal with the grid search problem. Compared with the exhaustive method, the intelligent search algorithm has a higher search efficiency. However, when the need for the location accuracy increases, the complexity of the algorithm increases exponentially.
Compared with TDOA measurement, the FDOA measurement requires more reference stations for satellite ephemeris calibration. Therefore, the tri-satellite TDOA localization system has lower cost and more reliable location results. The tri-satellite TDOA localization method is based on two TDOA lines of position (LOP) for intersection positioning. Affected by the ephemeris of the satellite and the location of the interference source, the TDOA LOP have two intersection points at a certain moment, and it is difficult to distinguish the correct location from the mirrored location of the interference source. The research in the existing literature mainly focuses on the optimization of TDOA measurement [24,25] and the solution of localization equations [26,27]. For the mirrored location that appears in the location process, the correct location can only be obtained based on the prior knowledge of the interference location in the conventional methods. When the prior knowledge is missed, the mirrored location cannot be eliminated. Therefore, it is urgent to propose an effective and robust localization method in which the mirrored location results are eliminated.
Aiming to address the problems above, this paper proposes a method to eliminate mirrored locations in tri-satellite TDOA localization based on the TDOA measured at multiple moments. First, the TDOA measure model based on cross-ambiguity function (CAF) is established. Then, a multi-moment TDOA localization equation set is derived based on the World Geodetic System (WGS)-84 earth ellipsoid model. The location results near the correct location of the interference source is directly obtained through the TDOA measurement data at multiple moments using the Newton iteration method. Finally, the precise iterative localizations are obtained using the multi-moment location result as the initial value in the iterative localization moment by moment. This method can obtain more accurate location results on the basis of eliminating the mirrored location. At the same time, it still has good robustness under the conditions of large deviation of the prior initial value and large TDOA measurement error.
The paper is organized as follows: Section 2 introduces the TDOA/FDOA measure model based on the CAF for satellite localization. Section 3 provides the multiple-moment TDOA localization method. Section 4 provides the proposed tri-satellite localization method. In Section 5, simulations are performed to validate the performance of our methods as well as the real data result analysis. We conclude this paper in Section 6.

2. TDOA/FDOA Measure Model

The model of the geodetic coordinate system is established as shown in Figure 1. The X-axis is pointing to the 0° longitude direction, the Z-axis is pointing to the center of the earth, and the Y-axis can be obtained according to the right-hand rule. It can be found that the main lobe of the interference transmitter is pointed towards the main satellite. At the same time, the two adjacent satellites can receive the sidelobe signal from the interference source.
Assume that the s 0 ( t ) , s 1 ( t ) and s 2 ( t ) are the signal received by the main satellite, the first adjacent satellite and the second adjacent satellite from the interference source in time t . The received signals from the interference source in three satellites from can be expressed as:
s 0 ( t ) = s ( t ) + n 0 ( t )
s 1 ( t ) = A 1 s ( t D 1 ) e j 2 π f d 1 t + n 1 ( t )
s 2 ( t ) = A 2 s ( t D 2 ) e j 2 π f d 2 t + n 2 ( t )
where s ( t ) is the complex envelope of the signal received by the interference source in the second adjacent satellite, A 1 and A 2 are the amplitude ratio coefficients between the signals, D 1 and D 2 are the delay time of the signal transmission for the first adjacent satellite and the second adjacent satellite, f d 1 and f d 2 are the frequency differences, n 0 ( t ) , n 1 ( t ) , and n 2 ( t ) are observation noise indicated.
Since the received interference signal relayed by the adjacent satellite is the signal from the side lobe of the interference transmitter that points to the main satellite, its signal-to-noise ratio is generally low, even submerged in the noise. The communication terminal downlink signal relayed by the main satellite has a high signal-to-noise ratio generally. Aiming at this problem, the cross-correlation technique is used to realize the detection of weak sidelobe signals.
The second-order complex CAF of the received signal from the main and adjacent satellites is expressed as:
A s 0 s 1 ( s ) ( f , τ ) = A e j 2 π f d 1 t A s ( s ) ( f f d 1 , τ D 1 )
A s 0 s 2 ( s ) ( f , τ ) = A e j 2 π f d 2 t A s ( s ) ( f f d 2 , τ D 2 )
A s ( s ) ( f , τ ) = 1 τ 0 T s ( t ) s * ( t + τ ) e j 2 π f t d t
According to the triangle inequality and the Schwartz inequality, we have:
A s ( s ) ( f , τ ) = A s ( s ) ( 0 , 0 )
The estimation of the time delay D ^ 1 and D ^ 2 and the frequency difference F ^ d 1 and F ^ d 2 can be obtained corresponding to the maximum value of CAF:
( F ^ d 1 , D ^ 1 ) = arg max { | A s 0 s 1 ( s ) ( f , τ ) | }
( F ^ d 2 , D ^ 2 ) = arg max { | A s 0 s 2 ( s ) ( f , τ ) | }

3. Multi-Moment TDOA Localization Model

As shown in Figure 1, the hyperboloid surfaces can be obtained based on the TDOA. The first hyperboloid surface is formed based on the TDOA measured from the main satellite and the first adjacent satellite at some moment. The second hyperboloid surface is formed based on the TDOA measured from the main satellite and the second adjacent satellite at some moment. At each moment, the two TDOA measurements of the three satellites including the main satellite and two adjacent satellite can be measured to form the TDOA hyperboloid surface 1. Under the assumption of zero elevation of the interference source, the location result of the interference source is the intersection of the two TDOA hyperboloid surfaces and the earth ellipsoid. The intersection is the result of tri-satellite location. In the situation shown in Figure 1, it can be found that the TDOA surface intersects at a straight line, and the straight line intersects with the earth’s ellipsoid surface at two points. One point is the correct location of the interference source target, and the other is the mirrored position. Without any prior information, it is difficult to distinguish the correct location from the mirrored location only by the TDOA of a single moment.
The following Figure 2 shows the intersection of the TDOA LOP when there are mirrored points in tri-satellite localization. It can be found that the LOP at multiple times gradually converge to one point at the real position, with the TDOA LOP at the mirrored point position intersect being more divergent. Therefore, the TDOA measurement results at multiple times can be used for localization to determine the true location of the interference source.
Assuming that the coordinate of the interference source is ( x , y , z ) , the coordinate of the main satellite at the first moment is ( x p 0 , y p 0 , z p 0 ) , the coordinate of the first adjacent satellite is ( x p 1 , y p 1 , z p 1 ) , and the second coordinate of the adjacent satellite is ( x p 2 , y p 2 , z p 2 ) . The measured TDOA between the interference source and the first adjacent satellite and second adjacent satellite at the first moment are τ p 1 and τ p 2 , respectively, which can be given by:
τ p 1 = ( r p 0 r p 1 ) / c
τ p 2 = ( r p 0 r p 2 ) / c
where c is the speed of light, r p 0 is the distance between the interference source of the main satellite, r p 1 is the distance between the interference source of the first adjacent satellite, r p 2 is the distance between the interference source of the second satellite, and p represents the p th moment.
According to the distance calculation method and the earth ellipsoid model, the tri-satellite localization equation set based on multiple moments can be given as follows:
{ ( x x 11 ) 2 + ( y y 11 ) 2 + ( z z 11 ) 2 ( x x 10 ) 2 + ( y y 10 ) 2 + ( z z 10 ) 2 = τ 11 ( x x 12 ) 2 + ( y y 12 ) 2 + ( z z 12 ) 2 ( x x 10 ) 2 + ( y y 10 ) 2 + ( z z 10 ) 2 = τ 12 ( x x P 1 ) 2 + ( y y P 1 ) 2 + ( z z P 1 ) 2 ( x x P 0 ) 2 + ( y y P 0 ) 2 + ( z z P 0 ) 2 = τ P 1 ( x x P 2 ) 2 + ( y y P 2 ) 2 + ( z z P 2 ) 2 ( x x P 0 ) 2 + ( y y P 0 ) 2 + ( z z P 0 ) 2 = τ P 2 x 2 N 2 + y 2 N 2 + z 2 [ N ( 1 e 2 ) ] 2 = 1
where N = a / 1 e 2 sin 2 B ^ k is the radius of the target local circle, a = 6378137   m is the long axis of the earth, and e 2 = 0.00669437999013 is the square of the first eccentricity of the WGS-84 earth ellipsoid [28].

4. Improved Tri-Satellite Localization Method

Consider that it is difficult to directly solve the nonlinear localization equations formed by the measurement results at multiple moments. The Newton iterative method is used to solve it. The above equations can be expressed as follows:
{ f 11 = ( x x 11 ) 2 + ( y y 11 ) 2 + ( z z 11 ) 2 ( x x 10 ) 2 + ( y y 10 ) 2 + ( z z 10 ) 2 τ 11 f 12 = ( x x 12 ) 2 + ( y y 12 ) 2 + ( z z 12 ) 2 ( x x 10 ) 2 + ( y y 10 ) 2 + ( z z 10 ) 2 τ 12   f P 1 = ( x x P 1 ) 2 + ( y y P 1 ) 2 + ( z z P 1 ) 2 ( x x P 0 ) 2 + ( y y P 0 ) 2 + ( z z P 0 ) 2 τ P 1 f P 2 = ( x x P 2 ) 2 + ( y y P 2 ) 2 + ( z z P 2 ) 2 ( x x P 0 ) 2 + ( y y P 0 ) 2 + ( z z P 0 ) 2 τ P 2 f e = x 2 N 2 + y 2 N 2 + z 2 [ N ( 1 e 2 ) ] 2 1 ,   where N = a 1 e 2 sin 2 B ^ k
We assume that F ( x , y , z ) = [ f 11 f 12 f P 1 f P 2 f e ] T , then the localization equation set can be formulated as:
F ( x , y , z ) = 0
The corresponding Hessian matrix is given as follows:
F ( x , y , z ) = [ f 11 x f 11 y f 11 z f 12 x f 12 y f 12 z f P 1 x f P 1 y f P 1 z f P 2 x f P 2 y f P 2 z f e x f e y f e z ] = [ x x 11 r 11 x x 10 r 10 y y 11 r 11 y y 10 r 10 z z 11 r 11 z z 10 r 10 x x 12 r 12 x x 10 r 10 y y 12 r 12 y y 10 r 10 z z 12 r 12 z z 10 r 10 x x P 1 r P 1 x x P 0 r P 0 y y P 1 r P 1 y y P 0 r P 0 z z P 1 r P 1 z z P 0 r P 0 x x P 2 r P 2 x x P 0 r P 0 y y P 2 r P 2 y y P 0 r P 0 z z P 2 r P 2 z z P 0 r P 0 2 x N 2 2 y N 2 2 z [ N ( 1 e 2 ) ] 2 ]
If the k th interference source coordinate in this iteration is P k = ( x k , y k , z k ) , it can be iterated according to the following formula:
P k + 1 = P k [ F ( x k , y k , z k ) ] 1 F ( x k , y k , z k )
The process of the tri-satellite interference source localization method based on multiple moments is summarized as follows:
Step 1
Calculate the distance difference between the corresponding three satellites and the interference source ( r p 1 r p 0 ) and ( r p 2 r p 0 ) according to the TDOA at P moments;
Step 2
Initialization, set the number of iterations k = 0 , and the initial position of the interference source P 0 ;
Step 3
Calculate the Hessian matrix F ( x k , y k , z k ) in the process of this iteration, and obtain the estimation result of the interference source position P k + 1 in the next iteration according to formula (16);
Step 4
Define the position error δ = ( x k + 1 x k ) 2 + ( x k + 1 x k ) 2 + ( x k + 1 x k ) 2 of the two iterations as the cost function in the iteration process. If it is less than a certain minimum threshold, stop the iteration and output the positioning result; otherwise, let k = k + 1 , and return to step 3 to continue the iteration.
Step 5
Precise localization moment by moment. The output location result is used as the initial value, and the single moment iteration method is used to obtain the precise value of the positioning result for the measurement result at each moment.
The flow chart of the proposed tri-satellite interference localization is shown in Figure 3.

5. Results

In the following simulation, the number of localization moments in single-moment positioning is p = 1 and the number of positioning moments in multi-moment positioning satisfies p > 1 . The proposed method is based on the multi-moment location results to remove the mirrored location. Then, the obtained location result near the correct location is used as the initial iteration for the single-moment localization method. The proposed method is the improved tri-satellite localization method based on the multi-moment TDOA localization.

5.1. Experiment on the Influence Factors of the Multi-Moment Localization

The performance of multi-moment interference source localization is mainly related to the number of selected moments and the total length of time. The following will explore its change relationship through simulation experiments. The location of the interference source is 125° E, 30° N, and the TDOA measurement root mean square error is 0.5 μs. First, the total time is simulated for 20 min, and the relationship between multi-moment location method and the location error of the proposed method varies with the number of moments as shown in Figure 4a. It can be found that the performance of the multi-moment localization improves with the increase of the number of times, and continues to approach the proposed method. The location error of the proposed method is basically unchanged. Then, the relationship between multi-moment localization and the location error of the proposed method with the total time is simulated when the number of the moments is 10 as shown in Figure 4b. It can be found that the performance of the multi-time localization improves with the increase of the total time, and it continues to approach the proposed method. The location error of the proposed method also remains stable and can achieve high-precision and rapid localization.

5.2. Experiment on the Influence of Iterative Initial Location

The following simulation experiments are used to compare the effect of the initial location of the iteration on the performance of the algorithm. The location of the interference source is 125° E, 30° N, and the TDOA measurement errors are independent and identically distributed random variables distributed uniformly over [ 0.5   μ s , 0.5   μ s ] . The number of Monte Carlo is 500. The time difference measured at three consecutive times is used to remove mirrored points for target localization. The time interval is five minutes, and the error of the initial iteration value is the root mean square error of the prior longitude and latitude deviation of the first iteration. Figure 5 shows the relationship between the location error of single-moment location, multi-moment moment, and the precise positioning method based on multi-moment measurement with the iterative initial value error when the iterative initial value error is between 0° and 5°. It can be found that, when the initial location error of the iteration increases gradually, the location error of the single-moment location degrades seriously. The multi-moment location method is basically not affected by the initial location error of the iteration, but with a poor performance compared with the single-moment method when the initial location error of the iteration is less than 3°. On the basis of multi-moment localization, the proposed method can obtain better location accuracy than the other two methods.

5.3. Experiment on the TDOA Measure Error

The following simulation experiments are used to compare the performance of the three methods over the TDOA measurement error. The location of the interference source is 125° E, 30° N, and the multi-moment location uses the TDOA measured at five consecutive times to remove mirrored points. The time interval is 4 min, and the root mean square of the initial location error of the iteration is 1°. The number of Monte Carlo is 500. It can be found in Figure 6 that, when the TDOA measurement error increases, the location error of the three methods increases. Among them, the single-moment location is mainly affected by the initial location error of the iteration when the TDOA measurement error is small, and it is not significantly affected by the TDOA measurement error. The multi-moment localization method is greatly affected by the time difference measurement error due to the method of solving the equations. The single-moment positioning method based on multi-moment measurement has good robustness when the time difference measurement error increases.

5.4. Real Data Result Analysis

In this section, the satellite real data processing results are presented to validate the effectiveness of the proposed algorithm. A signal from a cooperative interfering transmitter is illuminating the main satellite in synchronous orbit. The location of the interfering transmitter is [108.79° E, 32.32° N]. The two adjacent satellites can receive and relay the sidelobe signal of the interfering transmitter to the ground. This signal relayed from the three satellites is received by three C-band antennas, respectively. The CAF processing is used to measure the TDOA and FDOA of the main satellite and the adjacent satellites. The CAF processing results at a single moment are shown in Figure 7.
After measuring the TDOA and FDOA of the received data from the satellite, the multi-moment localization and the proposed method are used to solve the location of the interference transmitter. It can be seen from the change of the cost function in the iterative process that the required results can be obtained by iterating 2–3 in Figure 8a. From the intersection of the time difference lines in Figure 8b, it can be seen that, due to the error in the time difference measurement, the time difference lines at the real point position at multiple times cannot completely intersect at one point. Compared with the real point position, the time difference line at the mirrored point is more divergent, so the method based on multi-time measurement can be used to remove the mirrored point.
The location results and location errors of the multiple-moment localization method and the proposed method in four moments are shown in Table 1. It can be found that the location errors vary with the location time, which is caused by the motion of the satellite. Meanwhile, the performance of localization of the proposed method is better than the multiple-moment location method with smaller location errors. The multi-time localization algorithm can locate near the real position of the radiation source target, but the accuracy is affected by the time difference measurement error, and it is difficult to obtain high localization accuracy. The proposed method can greatly improve the localization accuracy at multiple moments on the basis of removing mirrored points through subsequent point-by-point localization.

6. Conclusions

The tri-satellite TDOA localization uses the intersection of two TDOA LOP to determine the location of the interference source. However, if there are two intersection points between the two TDOA LOP, it is difficult to determine the correct point location. This paper proposes a mirrored point elimination method based on multi-moment measurement. This method uses the TDOA measurement data at multiple times to determine the area where the interference source is located in reality based on the convergence of the TDOA LOP at multiple moments. Then, the result is used as the initial value of the Newton iteration method to solve the single-time positioning equation. The location results with a high precision can be obtained after removing the mirrored points. The simulation experiment analyzes the influencing factors on the positioning performance, and verifies the robustness under the larger measurement error and the deviation of the prior initial value.
Meanwhile, the location accuracy of the proposed method is verified in real testing. Compared with the conventional method, the proposed method can effectively and quickly obtain higher accuracy. The ability of the algorithm is confirmed and can provide guidance on the real tri-satellite location projects.

Author Contributions

All authors contributed extensively to this paper. L.H. provided the main idea and wrote the paper; M.J. and J.C. review and edit the paper, R.B., B.C. and G.L. modified the manu-script and provided many valuable suggestions; L.H. and P.H. conceived and designed the experiments. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Central Government Guides Local Technology Development Funds, Grant No. 206Z0701G.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors wish to thank the anonymous reviewers for their valuable comments on improving this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dempster, A.G.; Cetin, E. Interference Localization for Satellite Navigation Systems. Proc. IEEE 2016, 104, 1318–1326. [Google Scholar] [CrossRef]
  2. Ho, K.C.; Lu, X.; Kovavisaruch, L. Source Localization Using TDOA and FDOA Measurements in the Presence of Receiver Location Errors: Analysis and Solution. IEEE Trans. Signal Process. 2007, 55, 684–696. [Google Scholar] [CrossRef]
  3. Lin, Z.; Lin, M.; Cola, T.D.; Wang, J.B.; Cheng, J. Supporting IoT With Rate-Splitting Multiple Access in Satellite and Aerial Integrated Networks. IEEE Internet Things J. 2021. [Google Scholar] [CrossRef]
  4. Lin, Z.; Lin, M.; Champagne, B.; Zhu, W.P.; Al-Dhahir, N. Secure Beamforming for Cognitive Satellite Terrestrial Networks With Unknown Eavesdroppers. IEEE Syst. J. 2020, 15, 1–4. [Google Scholar] [CrossRef]
  5. Lin, Z.; Lin, M.; Champagne, B.; Zhu, W.P.; Al-Dhahir, N. Secure and Energy Efficient Transmission for RSMA-Based Cognitive Satellite-Terrestrial Networks. IEEE Wirel. Commun. Lett. 2020, 10, 251–255. [Google Scholar] [CrossRef]
  6. Haworth, D.P.; Smith, N.G.; Bardelli, R.; Clement, T. Interference localization for EUTELSAT satellites-the first european transmitter location system. Int. J. Satell. Commun. 2015, 15, 155–183. [Google Scholar] [CrossRef]
  7. José, F.; Guillermo, R.; Juan, M.T.; Brian, S. Survey on the Performance of Source Localization Algorithms. Sensors 2017, 17, 2666. [Google Scholar]
  8. Yu, H.; Huang, G.; Gao, J.; Liu, B. An Efficient Constrained Weighted Least Squares Algorithm for Moving Source Location Using TDOA and FDOA Measurements. IEEE Trans. Wirel. Commun. 2012, 11, 44–47. [Google Scholar] [CrossRef]
  9. Yang, K.; An, J.; Bu, X.; Sun, G. Constrained Total Least-Squares Location Algorithm Using Time-Difference-of-Arrival Measurements. IEEE Trans. Veh. Technol. 2010, 59, 1558–1562. [Google Scholar] [CrossRef]
  10. Wu, P.; Su, S.; Zuo, Z.; Guo, X.; Wen, X. Time Difference of Arrival (TDoA) Localization Combining Weighted Least Squares and Firefly Algorithm. Sensors 2019, 19, 2554. [Google Scholar] [CrossRef] [Green Version]
  11. Ming, S.; Le, Y.; Ho, D. Efficient Joint Source and Sensor Localization in Closed-Form. IEEE Signal Process. Lett. 2012, 19, 399–402. [Google Scholar]
  12. Lei, W.; Chen, B. High-precision hyperboloid location method using passive time-difference-of-arrival measurements. IET Radar Sonar Navig. 2013, 7, 710–717. [Google Scholar] [CrossRef]
  13. Huang, J.; Xue, Y.; Le, Y. An efficient closed-form solution for joint synchronization and localization using TOA. Future Gener. Comput. Syst. 2013, 29, 776–781. [Google Scholar] [CrossRef]
  14. Ho, K. Bias reduction for an explicit solution of source localization using TDOA. IEEE Trans. Signal Process. 2012, 60, 2101–2114. [Google Scholar] [CrossRef]
  15. You, K.H.; Lee, K.; Lee, H. An Optimized Solution for Hybrid TDOA/AOA based Geolocation using Nelder-Mead Simplex Method. IET Radar Sonar Navig. 2019, 13, 992–997. [Google Scholar]
  16. Sidorenko, J.; Schatz, V.; Bulatov, D.; Scherer-Negenborn, N.; Hugentobler, U. Self-Calibration for the Time Difference of Arrival Positioning. Sensors 2020, 20, 2079. [Google Scholar] [CrossRef] [Green Version]
  17. Qu, X.; Xie, L. An efficient convex constrained weighted least squares source localization algorithm based on TDOA measurements. Signal Process. 2016, 119, 142–152. [Google Scholar] [CrossRef]
  18. Lee, K.; Kim, S.; You, K. Iterative Regression Based Hybrid Localization for Wireless Sensor Networks. Sensors 2021, 21, 257. [Google Scholar] [CrossRef]
  19. Simonetto, A.; Leus, G. Distributed maximum likelihood sensor network localization. IEEE Trans. Signal Process. 2014, 62, 1424–1437. [Google Scholar] [CrossRef] [Green Version]
  20. Shen, H.; Ding, Z.; Dasgupta, S.; Zhao, C. Multiple source localization in wireless sensor networks based on time of arrival measurement. IEEE Trans. Signal Process. 2014, 62, 1938–1949. [Google Scholar] [CrossRef]
  21. Robles, G.; Fresno, J.; Martínez-Tarifa, J. Separation of Radio-Frequency Sources and Localization of Partial Discharges in Noisy Environments. Sensors 2015, 15, 9882. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Li, Q.; Chen, B.; Yang, M. Time Difference of Arrival Passive Localization Sensor Selection Method Based on Tabu Search. Sensors 2020, 20, 6547. [Google Scholar] [CrossRef] [PubMed]
  23. Cui, H.; Shu, M.; Song, M.; Wang, Y. Parameter Selection and Performance Comparison of Particle Swarm Optimization in Sensor Networks Localization. Sensors 2017, 17, 487. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Masmoudi, A.; Bellili, F.; Affes, S.; Stephenne, A. A Non-Data-Aided Maximum Likelihood Time Delay Estimator Using Importance Sampling. IEEE Trans. Signal Process. 2011, 59, 4505–4515. [Google Scholar] [CrossRef]
  25. Lee, B.; Kalker, T.; Schafer, R.W. Correlogram template matching for time-delay estimation. In Proceedings of the Acoustics, Speech and Signal Processing (ICASSP), 2011 IEEE International Conference, Prague, Czech Republic, 22–27 May 2011. [Google Scholar]
  26. Mikhalev, A.; Ormondroyd, R.F. Comparison of Hough Transform and Particle Filter Methods of Emitter Geolocation using Fusion of TDOA Data. In Proceedings of the Workshop on Positioning, Navigation and Communication, Dresden, Germany, 7–8 April 2011. [Google Scholar]
  27. 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]
  28. NIMA. Department of Defense World Geodetic System 1984, Its Definition and Relationships with Local Geodetic Systems; National Imagery and Mapping Agency: Bethesda, MD, USA, 1997. [Google Scholar]
Figure 1. Configuration of the tri-satellite interference localization.
Figure 1. Configuration of the tri-satellite interference localization.
Sensors 21 04483 g001
Figure 2. The intersection of the TDOA LOP in multiple moments.
Figure 2. The intersection of the TDOA LOP in multiple moments.
Sensors 21 04483 g002
Figure 3. The flow chart of the tri-satellite interference localization.
Figure 3. The flow chart of the tri-satellite interference localization.
Sensors 21 04483 g003
Figure 4. Experiments on the influence factors of location error. (a) Location errors versus the number of the location moments; (b) location error versus the length of location time.
Figure 4. Experiments on the influence factors of location error. (a) Location errors versus the number of the location moments; (b) location error versus the length of location time.
Sensors 21 04483 g004
Figure 5. Localization errors versus the error of the initial location error.
Figure 5. Localization errors versus the error of the initial location error.
Sensors 21 04483 g005
Figure 6. Localization errors versus the TDOA measure errors.
Figure 6. Localization errors versus the TDOA measure errors.
Sensors 21 04483 g006
Figure 7. Results of the CAF processing using the real data in actual tests. (a) The CAF processing result in TDOA-FDOA plane; (b) the mesh of the CAF processing results in TDOA-FDOA; (c) the FDOA slice of maximum CAF amplitude; (d) the TDOA slice of maximum CAF amplitude.
Figure 7. Results of the CAF processing using the real data in actual tests. (a) The CAF processing result in TDOA-FDOA plane; (b) the mesh of the CAF processing results in TDOA-FDOA; (c) the FDOA slice of maximum CAF amplitude; (d) the TDOA slice of maximum CAF amplitude.
Sensors 21 04483 g007aSensors 21 04483 g007b
Figure 8. Results of the localization using the real data in actual test. (a) The location error versus the number of iterations; (b) description of what is contained in the second panel. Figures should be placed in the main text close to the first time they are cited. A caption on a single line should be centered.
Figure 8. Results of the localization using the real data in actual test. (a) The location error versus the number of iterations; (b) description of what is contained in the second panel. Figures should be placed in the main text close to the first time they are cited. A caption on a single line should be centered.
Sensors 21 04483 g008
Table 1. Location results and location errors of the multiple time method and the proposed method.
Table 1. Location results and location errors of the multiple time method and the proposed method.
Localization MethodLocation TimeLocation ResultLocation Error
Multiple moment method20 December 2020 14:29:54108.7986° E 32.4250° N11.66 km
20 December 2020 15:29:54108.7949° E 32.3725° N5.84 km
20 December 2020 16:29:54108.8000° E 31.9776° N37.93 km
20 December 2020 17:29:54108.8281° E 31.7897° N58.86 km
Proposed method20 December 2020 14:29:54108.7899° E 32.2679° N5.80 km
20 December 2020 15:29:54108.7937° E 32.3174° N0.48 km
20 December 2020 16:29:54108.7938° E 32.2880° N3.57 km
20 December 2020 17:29:54108.7938° E 32.3015° N2.08 km
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Huo, L.; Bai, R.; Jiang, M.; Chen, B.; Chen, J.; Huang, P.; Liao, G. A Tri-Satellite Interference Source Localization Method for Eliminating Mirrored Location. Sensors 2021, 21, 4483. https://doi.org/10.3390/s21134483

AMA Style

Huo L, Bai R, Jiang M, Chen B, Chen J, Huang P, Liao G. A Tri-Satellite Interference Source Localization Method for Eliminating Mirrored Location. Sensors. 2021; 21(13):4483. https://doi.org/10.3390/s21134483

Chicago/Turabian Style

Huo, Lihuan, Rulong Bai, Man Jiang, Bing Chen, Jianfeng Chen, Penghui Huang, and Guisheng Liao. 2021. "A Tri-Satellite Interference Source Localization Method for Eliminating Mirrored Location" Sensors 21, no. 13: 4483. https://doi.org/10.3390/s21134483

APA Style

Huo, L., Bai, R., Jiang, M., Chen, B., Chen, J., Huang, P., & Liao, G. (2021). A Tri-Satellite Interference Source Localization Method for Eliminating Mirrored Location. Sensors, 21(13), 4483. https://doi.org/10.3390/s21134483

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