1. INTRODUCTION
Developed by the US Department of Defense for military purpose, Global Positioning System (GPS) is a radio navigation system which provides location, speed and time information to users. Recently due to the increasing demand for realtime location information, GPS is utilized in various fields not only in the military but also in private sectors such as the medical, port, and logistics industries. However, devices such as Personal Privacy Devices (PPD) are employed to avoid the detection of one’s location by GPS. Available on the internet and much affordable than GPS receiver, PPD can be easily purchased. As the electric power of GPS signals received by a receiver is about -158 dBW, PPD may interfere GPS signals even though it transmits signals with low electric power of about 1 dBW. For instance, a truck driver used PPD to avoid detection of his route near Newark International Airport in the US in 2010, and the signals transmitted from the PPD caused a failure in the high-tech instrument landing system of the airport (Pullen et al. 2012). TDOASuch a failure in a system using GPS may lead to a large-scale accident. Thus, recently, some institutions analyzed the type of signals transmitted by PPD and their impact on Global Navigation Satellite System (GNSS) receiver (Kraus, Bauernfeind, & Eissfeller 2011, Mitch et al. 2011, Grabowski 2012). University of Federal Armed Forces in Germany investigated on interferences such as PPD. According to the study results, it was found that PPD transmitted swept Continuous Wave (CW) signals of the frequency adjacent to GPS L1 signal band, and had a bandwidth smaller than 1 kHz. As the signals transmitted from PPD may cause jamming in a precise localization system based on GPS, many institutions are engaged in vigorous research activities to reduce interferences by PPD (Gromov et al. 2000, Balaei, Dempster & Barnes 2006, Brown & Reynolds 2010, Cetin, Thompson & Dempster 2011, Chang & Huang 2011, Xu & Trinkle 2011). Most studies, however, were mainly about analyzing signals transmitted by interferers and categorizing them based on their characteristics, so there was no study on localization of interferences. Therefore a study on a technique to locate interference is needed to directly reduce the demage of interferences.
Most typical techniques to locate interference are Array of Arrival (AOA) method which uses the arrival angle of signal and Time Difference of Arrival (TDOA) method which uses arrival time differences of signal. AOA method has restrictions that it needs to use array antenna and synchronize phases between RF/IF channels (Niculescu & Nath 2003), and TDOA method has the restriction of requiring time synchronization between sensors (Smith & Abel 1987, Chan & Ho 1994a,b). Moreover, the AOA method has the drawback of causing a huge error if signals are scattered by surrounding environment or sensor is located near a source, this paper using the TDOA method estimates the location of interferences which transmit swept CW signals. However, if the integration time for cross-correlation to obtain a TDOA measurement is longer than sweep time, problem of ambiguity arises as many TDOA measurements are obtained. Because of such ambiguity, general TDOA method cannot be used for localization. Thus, this paper design a localization algorithm based on TDOA method that can accurately locate an interferer transmitting swept CW signals by resolving problem of ambiguity. In addition, feasibility of the designed algorithm has been verified by simulation results.
2. CHARACTERISTIC OF TDOA MEASUREMENT OF SWEPT CW SIGNALS
2.1 Swept CW Signal Model
Kraus analyzed that PPD transmit swept CW signals, and generally swept CW signals are sinusoidal waves in which frequencies vary by time. Accordingly, the model of swept CW signals is calculated according to Eq. (1) as follows (Kraus, Bauernfeind, & Eissfeller 2011).
\(x(t)=a\sin{\left[2\pi\left(f_0+kt/2\right)t\right]},(0\le t<T_{SW}) \) (1)
Here, \(a\) is amplitude and has a constant value, \(f_0\) is initial frequency, \(k\) is chirp rate, and \(T_{SW}\) is sweep time. The Eq. (1) can be changed into Eq. (2) using discrete time.
\(x(n)=a\sin{\left[2\pi\left(f_0+knT_{\mathrm{SP}}/2\right)nT_{\mathrm{SP}}\right]}, \left (\mathrm{when} \ n \ \mathrm{is \ an \ integer,}\ 0\le n<\frac{T_{\mathrm{SW}}}{T_{\mathrm{SP}}} \right )\) (2)
Here, \(T_{SP}\) is sampling period, and B stands for bandwidth. The bandwidth can be defined according to Eq. (3).
\(Β≜{kT_{sw} \over 2}\) (3)
Eq. (2) can be represented as Fig. 1.
Fig. 1. Swept CW interference signal model.
2.2 Characteristics of TDOA Measurement of Swept CW Signals
In order to locate interferences which transmit swept CW signals by using the TDOA method, TDOA measurement need to be obtained. The TDOA measurement for two arbitrary sensors can be calculated with cross-correlation function, and the function is defined in Eq. (4).
\(R_{ri}(\tau)\overset{\triangle}= E\left[s_r(t)s_i(t-\tau)\right]=\frac{1}{T}\int_{0}^{T}{s_r(t)s_i(t-\tau})d\tau\) (4)
Here, \(S_r (t)\) is a signal received at reference sensor, \(S_i (t)\) is a signal received at \(i\) th sensor, which is transmitted signal from an interferer, \(T\) is the integration time, and \(\tau\) is time delay. As for interference signal, the paper sets to 1 ms as the period of interference signal cannot be known.
As mentioned in the introduction, when cross-correlating swept CW signals, an issue of ambiguity may arise as the peak of cross-correlation appears repeatedly at every sweep time. Ambiguity arises since sweep time is shorter than integration time. Fig. 2 shows the result of cross-correlating four sensors, while fixing one as a reference sensor with the sweep time of 200 us and the sampling period of 20 ns in Fig. 1.
Fig. 2. Cross-correlation value of swept CW signals.
It can be seen from Fig. 2 that the peak of cross-correlation appears repeatedly at every sweep time. We cannot know true peak of them. Therefore, TDOA measurement cannot be obtained.
3. DESIGN OF AN ENHANCED TDOA METHOD FOR LOCATION OF SWEPT CW INTERFERENCES
3.1 TDOA Measurement Model
As identified in Section 2, in case of swept CW interference signals, the characteristic of the TDOA measurement has an ambiguity that the cross-correlation value appear repeatedly at every sweep time. Thus, a TDOA measurement model is necessary to track the location of swept CW interferer.
Fig. 3. Delay on initial peak.
Fig. 3 shows delays of the first peak for samples from 0 to 100. \(\Delta \tau_{12}\) is TDOA measurement obtained by sensor 1 (reference sensor) and sensor 2, and a delay on cross-correlation value which happens first among repeated cross-correlation values. Therefore TDOA measurement of swept CW interference signals can be calculated as Eq. (5).
\(\tau_{ri}=\Delta\tau_{ri}+N_{ri}T_{SW}\) (5)
Here, \(\tau_{ri}\) is TDOA measurement obtained by reference sensor and \(i\) th sensor, \(\Delta \tau_{ri}\) is the first TDOA measurement among those repeatedly obtained from reference sensor and \(i\) th sensor, and \(N_{ri}\) is an integer ambiguity which needs to estimate as ambiguity.
3.2 Localization Algorithm
As a localization method using the TDOA method, there are direct solutions and indirect solutions. Spherical Intersection (SX), Spherical Interpolation (SI), Quadratic Correction Least Square (QCLS) are included in direct solutions and Least Squares (LS) is included in indirect solutions. This paper use LS method which is known to find more accurate value compared to direct solutions (Foy 1976, Smith & Abel 1987, Chan & Ho 1994a).
Eq. (5) that is equation with TDOA measurement can be expressed as Eq. (6) that is equation with range measurement.
\(c\tau_{ri}=c\Delta\tau_{ri}+cN_{ri}T_{SW}\) (6)
Here, \(c\) is propagation velocity. A measurement used in the designed algorithm \(c\Delta\tau_{ri}\) is that is calculated by using the delay obtained from the first peak happened during the integration time. When the number of sensors is N, N-1 number of equations are gained by substituting \(r\) with 1 and \(i\) with 2~N in Eq. (6). Thus those equations can be arranged by matrix and can be expressed as Eq. (7).
\(c\left[\begin{matrix}\Delta\tau_{12}\\\Delta\tau_{13}\\\vdots\\\Delta\tau_{1N}\\\end{matrix}\right]=\left[\begin{matrix}G_{x2}&G_{y2}\\G_{x3}&G_{y3}\\\vdots&\vdots\\G_{xN}&G_{yN}\\\end{matrix}\right]\left[\begin{matrix}\delta x\\\delta y\\\end{matrix}\right]-cT_{sw}\mathbf{I}_{\left((N-1)\times(N-1)\right)}\left[\begin{matrix}N_{12}\\N_{13}\\\vdots\\N_{1N}\\\end{matrix}\right]\) \(=\left[\begin{matrix}G_{x2}&G_{y2}&-cT_{sw}&0&\cdots&0\\G_{x3}&G_{y3}&0&-cT_{sw}&\cdots&0\\\vdots&\vdots&0&0&\ddots&0\\G_{xN}&G_{yN}&0&0&\cdots&-cT_{sw}\\\end{matrix}\right]\left[\begin{matrix}\delta x\\\delta y\\N_{12}\\N_{13}\\\vdots\\N_{1N}\\\end{matrix}\right]\)
\(\mathbf{Z}_{SCW}=\mathbf{G}_{SCW}\boldsymbol{\delta}_{SCW}\) (7)
Here, \(G_{xN}\) is \( \frac{\hat{x}-x_N}{\sqrt{\left(\hat{x}-x_N\right)^2+\left(\hat{y}-y_N\right)^2}}-\frac{\hat{x}-x_1}{\sqrt{\left(\hat{x}-x_1\right)^2+\left(\hat{y}-y_1\right)^2}}\), \(G_{yN}\) is \(\frac{\hat{y}-y_N}{\sqrt{\left(\hat{x}-x_N\right)^2+\left(\hat{y}-y_N\right)^2}}-\frac{\hat{y}-y_1}{\sqrt{\left(\hat{x}-x_1\right)^2+\left(\hat{y}-y_1\right)^2}}\), \((\delta x, \delta y)\) is delta position, and \(\mathbf{I}\) is unit matrix. \(\boldsymbol{\delta}_{SCW}\) consisted of delta position and integer ambiguities in Eq. (7) can be calculated by LS method shown in Eq. (8).
\(\boldsymbol{\delta}_{SCW}=\left(\mathbf{G}_{SCW}^T\mathbf{Q}^{-1}\mathbf{G}_{SCW}\right)^{-1}\mathbf{G}_{SCW}^T\mathbf{Q}^{-1}\mathbf{Z}_{SCW}\) (8)
Here, \(\mathbf{Q}\) is covariance matrix of measurement error. Thus, to estimate the solution, \(\boldsymbol{\delta}_{SCW}\) computed in Eq. (8) add to initial values as Eq. (9).
\( \left [ \begin{matrix} x \\ y \\ N_{12} \\ N_{13} \\ ⋮ \\ N_{1N}\end{matrix} \right ] = \left [ \begin{matrix} x_0 \\ y_0 \\ N_{(12)_0} \\ N_{(13)_0} \\ ⋮ \\ N_{1N_0} \end{matrix} \right ] + \left [ \begin{matrix} \delta_x \\ \delta_y \\ \hat{N}_{12} \\ \hat{N}_{13} \\ ⋮ \\ \hat{N}_{1N} \end{matrix} \right ] \) (9)
Finally, the position of an interferer and integer ambiguity can be calculated by repeating from Eq. (3) to Eq. (9) and updating until the magnitude of satisfies specific stop condition.
Fig. 4 shows the flowchart of the designed algorithm in this paper to estimate the position of an interferer which transmits the swept CW signal. This paper presets the stop condition with the magnitude of \(\boldsymbol{\delta}_{SCW}\) below 0.001 or repetition over 10 times.
Fig. 4. Flow chart.
4. SIMULATION AND RESULTS
In order to estimate position of interferer by designed algorithm, MATLAB was utilized for computer simulation. Table 1 shows the simulation environment.
Table 1. Simulation environment.
Parameter | Set value |
Sensor arrangement Location of interferer Type of noise signal Power of noise signal Sampling frequency IF frequency Sweep time Integration time Initial position on algorithm Number of repeating algorithm |
Arrange four sensors as square form (2000, 2000), (-2000, 2000), (-2000, -2000), (2000, -2000) m (868.2, 4924) m AWGN (Additive White Gaussian Noise) 1 mW 5.714 MHz 1.134 MHz 10 µs 1 ms Place on (0,0) centered on sensor arrangement within 10 times |
Fig. 5 shows cross-correlation value of swept CW signals created from set sampling frequency and sweep time for the simulation. The TDOA measurement to locate an interferer from several cross-correlation peaks generated at sweep times were found by calculating the integer ambiguities from the designed algorithm, and Table 2 shows actual TDOA, TDOA measurement, and the difference between the two.
Fig. 5. Cross-correlation value for swept CW signals created for simulation.
Table 2. Actual TDOA, TDOA measurement, error.
Actual TDOA (Actual distance/ minute difference) |
TDOA measurement (Estimated measure of distance/minute difference) |
TDOA error (Error of distance/ minute difference) |
|
\(\tau_{12}\) \(\tau_{13}\) \(\tau_{14}\) |
3.202 ms (960.5 m) 14.531 ms (4359.2 m) 12.935 ms (3880.5 m) |
3.481 ms (1044.2 m) 14.777 ms (4433.1 m) 13.018 ms (3905.4 m) |
0.279 ms (83.7 m) 0.246 ms (73.9 m) 0.083 ms (24.9 m) |
It can be seen from Fig. 5 and Table 2 that the designed algorithm can find TDOA measurement by solving the integer ambiguities for happened problem of ambiguity with the cross-correlation peaks repeatedly at every sweep time.
Fig. 6 shows the localization of interferer transmitting swept CW signal by using the designed TDOA method.
Fig. 6. Results of interference localization using the designed algorithm.
As shown in Fig. 6 and Table 3, the designed algorithm can be estimated for location of an interferer transmitting swept CW signal. Thus, It can be confirmed that the designed algorithm can estimate location of an interferer by solving ambiguity of TDOA measurements.
Table 3. True position, estimated position, position error.
True position | Estimated position | Position error | |
x m y m |
868.2 4924 |
972.6 5034.4 |
104.4 110.4 |
5. CONCLUSION AND FURTHER STUDY
This paper modeled the TDOA measurement for swept CW interference signals, and designed the algorithm to locate position of an interferer transmitting swept CW signal using TDOA method. In addition, we confirmed that the designed algorithm can determine the integer ambiguity to resolve problem of ambiguity and can locate position of an interferer transmitting swept CW interference signal.
A further study will be conducted to analyze changes of performance according to sampling frequency, and design an advanced method with higher accuracy to locate swept CW interferers.
ACKNOWLEDGMENTS
This work was supported by a grant from Korea Research Council of Fundamental Science & Technology funded by Ministry of Education, Science and Technology in 2012 [Project: A Study on Satellite based Position Tracking Technology for Calamity Prevention and Public Safety Improvement].
참고문헌
- Balaei, A. T., Dempster, A. G., & Barnes, J. 2006, A novel approach in detection and characterization of CW interference of GPS signal using receiver estimation of C/N0, Position, Location, and Navigation symposium, 2006 IEEE/ION, pp.1120-1126.
- Brown, A. & Reynolds, D. 2010, Civil Applications of the GPS Jamming Detection and Location (JLOC) System, Presentations in National Physical Laboratory.
- Cetin, E., Thompson, J. R., & Dempster, G. 2011, Interference Localisation within the GNSS Environmental Monitoring System (GEMS), Proceedings of IGNSS Symposium 2011.
- Chan, Y. T. & Ho, K. C. 1994a, A Simple and Efficient Estimator for Hyperbolic Location, IEEE Transactions On Signal Processing, 42, 1905-1915. https://doi.org/10.1109/78.301830
- Chan, Y. T. & Ho, K. C. 1994b, An efficient closed-form Localization solution from time difference of arrival measurements, IEEE Transactions on Acoustics, Speech, and Signal Processing, 2, 393-396
- Chang, C.-L. & Huang, G.-S. 2011, A Compressive Sampling Approach to Narrowband Interference Elimination for GNSS, Proc. IGNSS Symposium 2011.
- Foy, W. H. 1976, Position-location Solution by Taylor-series Estimation, IEEE Transactions on Aerospace and Electronic Systems, AES-12(2), pp. 187-194. https://doi.org/10.1109/TAES.1976.308294
- Grabowski, J. 2012, Field Observations of Personal Privacy Devices, Proceedings of ION ITM 2012.
- Gromov, K., Akos, D., Pullenm, S., Enge, P., & Parkinson, B. 2000, GIDL: Generalized Interference Detection and Localization System, Proceedings of ION GPS 2000, pp. 19-22.
- Kraus, T., Bauernfeind, R., & Eissfeller, B. 2011, Survey of In-Car Jammers - Analysis and Modeling of the RF signals and IF samples (suitable for active signal cancellation), Proceedings of ION GNSS 2011, pp. 430-435.
- Mitch, R. H., Dougherty, R. C., Psiaki, M. L., Powell, S. P., & O'Hanlon, B. W. 2011, Signal Characteristics of Civil GPS Jammers, Proceedings of ION GNSS 2011, pp. 1907-1919.
- Niculescu, D. & Nath, B. 2003, Ad Hoc Positioning System (APS) Using AOA, INFOCOM 2003, Twenty-Second Annual Joint Conference of the IEEE Computer and Communications, IEEE Societies, 3, 1734-1743.
- Pullen, S., Gao, G. Tedeschi, Carmen, T., & Warburton, J. 2012, The Impact of Uninformed RF Interference on GBAS and Potential Mitigations, Proceedings of ION International Technical Meeting.
- Smith, J. O. & Abel, J. S. 1987, Closed-Form Least-Squares Source Location Estimation from Range-Difference Measurements, IEEE Transactions on Acoustics, Speech and Signal Processing, 35, 1661-1669. https://doi.org/10.1109/TASSP.1987.1165089
- Xu, Z. & Trinkle, M. 2011, Interference Angle of Arrival Estimation within the GNSS Environmental Monitoring System (GEMS) using Antenna Arrays, proceedings of IGNSS Symposium 2011.