Therefore, goal of this work is to propose techniques to overcome the existing limitations in antenna array processing and snapshot processing for modern GNSS receivers. The proposed techniques not only reduce the implementation cost but also leverage the distributed data processing ability.
MINISTRY OF EDUCATION AND TRAINING HANOI UNIVERSITY OF SCIENCE AND TECHNOLOGY NGUYEN DINH THUAN ROBUST SIGNAL PROCESSING TECHNIQUES FOR MODERN GNSS RECEIVERS Major: Computer Engineering Code No.: 9480106 COMPUTER ENGINEERING DISSERTATION SUPERVISORS: Assoc Prof Ta Hai Tung Prof Letizia Lo Presti Hanoi - 2019 STATEMENT OF ORIGINALITY AND AUTHENTICITY I confirm that my dissertation is an original and authentic piece of work written by myself The data, results in the thesis is reliable and has never been published by others I further confirm that I have fully referenced and acknowledged all material incorporated as secondary resources in accordance with the regulations Hanoi, SUPERVISORS PHD STUDENT PGS.TS Tạ Hải Tùng Nguyễn Đình Thuận Prof Letizia Lo Presti ACKNOWLEDGEMENTS I would like to express my gratitude to Hanoi University of Technology, Graduate School, School of Information and Communication Technology, Department of Computer Engineering and Politecnico di Torino, NavSaS group for creating favorable conditions for me to work and study I would like to express my special thanks to my supervisors, Assoc Ta Hai Tung and Prof Letizia Lo Presti The supervisors have always been helpful, giving great advice, scientific orientations so that I can develop and complete my research Sincerely thank the lecturers, colleagues in the Department of Computer Engineering, School of Information and Communication Technology, Hanoi University of Science and Technology where I work, study and carry out research projects for the enthusiastic in helping and encouraging me during the research With gratitude to teachers, scientists, colleagues and close friends for encouraging and supporting me in the research process Finally, I would like to express my deep gratitude to my family for encouraging me to overcome all obstacles to complete this thesis Nguyen Dinh Thuan TABLE OF CONTENTS STATEMENT OF ORIGINALITY AND AUTHENTICITY ACKNOWLEDGEMENTS TABLE OF CONTENTS LIST OF ACRONYMS LIST OF TABLES LIST OF FIGURES INTRODUCTION 13 FUNDAMENTAL BACKGROUND 18 1.1 GNSS positioning principle 18 1.2 History and development of GNSS 19 1.3 GNSS Threats 20 1.3.1 Multipath 21 1.3.2 Atmosphere 21 1.3.3 Interference 21 1.3.4 Spoofing 21 1.3.5 GNSS Segment errors 21 1.3.6 Cyber Attacks 22 1.4 1.4.1 Signal Conditioning and Sampling 22 1.4.2 Acquisition 23 1.4.3 Tracking and Data Demodulation 23 1.4.4 Positioning Computation 24 1.5 Countermeasures to GNSS Threats 25 1.5.1 Antenna array processing techniques 25 1.5.2 Frontend and Digital Signal Conditioning based techniques 28 1.5.3 Correlator/Tracking and PVT based techniques 29 1.6 GNSS Receiver Architecture 22 GNSS Simulator and effect of sampling frequency 30 GNSS SIGNAL SIMULATOR DESIGN AND IMPLEMENTATION 32 2.1 Modeling methodology 32 2.2 Overview of the modeling of antenna array signals in GNSS receivers 32 2.2.1 General model of the received signal in GNSS receivers 33 2.2.2 Interference 37 2.2.3 Multipath 38 2.2.4 Noise 39 2.3 Effect of sampling frequency on the positioning performance 39 2.3.1 Residual code phase estimation 40 2.3.2 Correlation output calculation 40 2.3.3 Effect of sampling frequency on correlation shape and DLL discriminator function 42 2.3.4 Effect of the sampling frequency and the integration period selection 42 2.3.5 Effect on the presence of Doppler and local oscillator (LO) clock drift 45 2.3.6 Theoretical code tracking loop error estimate 46 2.3.7 Theoretical results evaluation by simulated, and numerical models 49 2.3.8 Effect of Doppler and coherent integration period 50 2.4 Sampling Frequency Effect Mitigation Technique 53 2.4.1 2.5 Performance verification 57 2.5.1 Verification of the simulated antenna array signals 58 2.5.2 Antenna distortion simulation 64 2.5.3 Verification of multipath simulation 66 2.6 Receiver implementation 55 Conclusion 67 ANTENNA ARRAY PROCESSINGS FOR GNSS RECEIVERS 69 3.1 The proposed solution for synchronizing separated antenna array element 69 3.1.1 Determining the samples difference 70 3.1.2 Determining the clock phase shift 71 3.2 Implementation a low-cost antenna array 75 3.3 Antenna array frontend verification 76 3.3.1 Phase difference between frontends 76 3.3.2 Carrier to noise ration improvement 77 3.4 Conclusion 78 GNSS SNAPSHOT PROCESSING TECHNIQUE FOR GNSS RECEIVERS 80 4.1 Proposed Design of GNSS Snapshot Receiver 80 4.1.1 GNSS Grabber 80 Implementation of GNSS Grabber 80 Firmware Architecture 81 4.2 Server Software 81 4.2.1 GNSS signal acquisition 81 4.2.2 Combined Doppler and Snapshot Algorithm 84 4.3 Loosely coupled Snapshot GNSS/INS 89 4.4 Tightly coupled Snapshot GNSS/INS 96 4.5 Results 97 4.5.1 Standalone Snapshot GNSS Receiver 97 4.5.2 Snapshot GNSS/INS Integration 102 4.6 Conclusion 104 CONCLUSIONS AND FUTURE WORKS 105 PUBLICATIONS 107 REFERENCES 109 APPENDIX 116 A Correlation output calculation 116 B Error analysis for coherent early minus late DLL 117 LIST OF ACRONYMS Acronym Meaning ADC Analog to Digital Converter AGC Automatic Gain Control AWGN Additive White Gaussian Noise BB BaseBand BOC Binary Offset Carrier BPSK Binary Phase Shift Keying C/A Coarse/Acquisition C/N0 Carrier-to-Noise-Density Ratio CDC Conventional Differential Combination CDMA Code Division Multiple Access CRC Cyclic Redundancy Check CS Commercial Service DLL Delay Lock Loop DFT Discrete Fourier Transform DSP Digital Signal Processor EGNOS European Geostationary Navigation Overlay Service EU European Union FEC Forward Error Correction FFT Fast Fourier Transform FPGA Field Programmable Gate Array FOC Full Operational Capability GLONASS Global Orbiting Navigation Satellite System I Inphase IF Intermediate Frequency Q Quadrature PVT Position Velocity Time SDR Software Defined Radio LIST OF TABLES Table 2.1: GNSS Simulator Features 57 Table 2.2: The coordinate of elements 58 Table 2.3: The direction of visible satellites 59 Table 2.4: The carrier phase relative to the first element of each satellite at the four elements of the array 59 Table 2.5: The simulation scenario 60 Table 2.6: Estimated carrier phase using the post-correlator beamforming tracking loop 62 Table 4.1: Configuration of the GPS grabber 97 Table 4.2: Information of acquired satellites 99 LIST OF FIGURES Figure 1.1: Satellite navigation principle 18 Figure 1.2: Typical GNSS Threats 20 Figure 1.3: Signal conditioning and sampling stage 22 Figure 1.4: Acquisition Architecture 23 Figure 1.5: Tracking Architecture 23 Figure 1.6: Transmission time estimation in GNSS receivers 24 Figure 1.7: Interference mitigation techniques in GNSS receivers 25 Figure 1.8: The traditional low-cost architecture of antenna array for GNSS applications 27 Figure 1.9: The correlation between GPS signal grabbed by antenna array 28 Figure 1.10: Spectrum and histogram of GNSS signal in the absence of interference 28 Figure 1.11: Snapshot positioning architecture 29 Figure 2.1: Geometry of antenna array 33 Figure 2.2: The model of the received signal for a single antenna 33 Figure 2.3: GPS multi-antenna frontend 34 Figure 2.4: Flowchart of the simulator 35 Figure 2.5: Bandlimited Gaussian interference model 38 Figure 2.6: Multipath model 38 Figure 2.7: Effect of sampling frequency on the positioning performance 39 Figure 2.8: Residual code phases versus the number of samples per code chip with 4fc < fs < 5fc 40 Figure 2.9: Normalised correlator and EML discriminator functions for different sampling frequencies Results are obtained by correlating the incoming signal with various local generated replica signals that have the time delay from−Tc to Tc with step = 10-2Tc 42 Figure 2.10: Correlation shapes for ms integration with various sampling frequencies 43 Figure 2.11: Ambiguous synchronization between a local PRN code and two different incoming analog signals of the same PRN sequence, but with slightly differing code phase offset 43 Figure 2.12: Correlation shapes and their errors with respect to the ideal correlation at a sampling frequency fs =16.3676 MHz using various coherent integration periods 44 Figure 2.13: Representation of code tracking loop [54] 46 Frequency on GNSS Code Tracking." The Journal of Navigation 71, no (2018): 257280 (SCIE) 108 REFERENCES [1] ICAO, "Report of the Special Communications/Operations Divisional Meeting," 1995 [2] K Borre, D M Akos, N Bertelsen, P Rinder and S H Jensen, A software-defined GPS and Galileo receiver: a single-frequency approach, Springer Science & Business Media, 2007 [3] G ICD, "Global Positioning Systems Directorate System Engineering \& Integration Interface Specification IS-GPS-200H," Navstar GPS Space Segment/Navigation User Interfaces, 2013 [4] I Glonass, "Glonass interface control document," Russian Institute of Space Device Engineering: Moscow, Russia, 2008 [5] I BeiDou, "BeiDou navigation satellite system signal in space interface control document open service signal B1I (Version 1.0)," BeiDou, ICD, 2012 [6] Ignacio Fernández Hernández, "Resilient Position, Navigation and Timing," ITSNT, 2017 [7] M Cuntz, A Konovaltsev, L Kurz, C Hättich, G Kappen,, "Interference and Countermeasures for GNSS Receivers," TUM Navigation Colloquium, 2011 [8] S Alban, Design and performance of a robust GPS/INS attitude system for automobile applications, Citeseer, 2004 [9] D W a A N a H C C Allan, The science of timekeeping, Hewlett-Packard, 1997 [10] D N a A M a A D M Aloi, "A methodology for the evaluation of a GPS receiver performance in telematics applications," IEEE Transactions on Instrumentation and Measurement, vol 1, pp 11-24, 56 [11] S a A D Backen, "Antenna array calibration using live GNSS signals," in ESA Workshop on Satellite Navigation User Equipment Technologies, 2006 [12] J a C Y.-H a D L D S a L S a E P a A D a L J Seo, "A real-time capable software-defined receiver using GPU for adaptive anti-jam GPS Sensors," Sensors, vol 11, pp 8966 8991 , 2011 109 [13] Chen, Yu-Hsuan, "A study of geometry and commercial off-the-shelf (COTS) antennas for controlled reception pattern antenna (CRPA) arrays," in Proceedings of ION GNSS , 2012 [14] Chen, Yu-Hsuan and Juang, Jyh-Ching and Seo, Jiwon and Lo, Sherman and Akos, Dennis M and De Lorenzo, David S and Enge, Per, "Design and implementation of real-time software radio for anti-interference GPS/WAAS sensors," Sensors, vol 12, pp 13417 13440, 2012 [15] Akos, Dennis M, "Who's afraid of the spoofer? GPS/GNSS spoofing detection via automatic gain control (AGC)," Navigation, pp 281 290, 2012 [16] I J a B J R a E S W a P H.-G a O H S a K M G Gupta, "An experimental study of antenna array calibration," IEEE Transactions on antennas and propagation, vol 51, pp 664-667, 2003 [17] B C a S C M S Ng, "Sensor-array calibration using a maximum-likelihood approach," IEEE Transactions on Antennas and Propagation, vol 44, pp 827-835, 1996 [18] C M a G I J Church, "Calibration of GNSS adaptive antennas," in Proceedings of the 22nd International Technical Meeting of The Satellite Division of The Institute of Navigation (ION GNSS 2009), 2009 [19] T M a R T a M D King, "Detection and reduction of periodic jamming signals in GPS receivers and methods therefor" Patent US Patent 8,253,624, 28 2012 [20] Borio, Daniele, "A multi-state notch filter for GNSS jamming mitigation," in Localization and GNSS (ICL-GNSS), 2014 International Conference on, 2014 [21] Borio, Daniele, "Swept GNSS jamming mitigation through pulse blanking," in Navigation Conference (ENC), 2016 European, 2016 [22] Balaei, Asghar Tabatabaei and Motella, Beatrice and Dempster, Andrew, "A preventative approach to mitigating CW interference in GPS receivers," GPS solutions, vol 12, pp 199 209, 2008 [23] Angrisano, Antonio and Gaglione, Salvatore and Gioia, Ciro, "RAIM algorithms for aided GNSS in urban scenario," in Ubiquitous Positioning, Indoor Navigation, and Location Based Service (UPINLBS), 2012, 2012 110 [24] W a C K H a L J a K L a S H a L H K Yoo, "Coarse-time Positioning without Continuous GPS Signal Tracking," in International Global Navigation Satellite Systems Society, 2016 [25] I F Hernandez, SNAPSHOT AND AUTHENTICATION TECHNIQUES FOR SATELLITE NAVIGATION, 2015: Aalborg University [26] F S T Van Diggelen, A-GPS: Assisted GPS, GNSS, and SBAS, Artech House, 2009 [27] S C Wu, W I Bertiger, D Kuang, S M Lichten, S Nandi, L J Romans and J M Srinivasan, "MicroGPS for low-cost orbit determination," Space, vol 5, pp 2-3, 1997 [28] N Othieno and S Gleason, "Combined Doppler and time free positioning technique for low dynamics receivers," in Position Location and Navigation Symposium (PLANS), 2012 IEEE/ION, 2012 [29] S a A D a N M L Backen, "Post-processing dynamic GNSS antenna array calibration and deterministic beamforming," in International Technical Meeting of the Satellite Division of the Institute of Navigation: 16/09/2008-19/09/2008, 2008 [30] D S a o De Lorenzo, Navigation accuracy and interference rejection for GPS adaptive antenna arrays, Stanford University, 2007 [31] E a H C Kaplan, Understanding GPS: principles and applications, Artech house, 2005 [32] I Galileo, "Galileo open service, signal in space interface control document (OS SIS ICD)," European space agency/European GNSS supervisory authority, 2008 [33] Spirent, "Fundamentals of GPS Threats," White Paper, 2015 [34] B a S J Parkinson, Global Positioning System: Theory and Applications, vol 1; Progress in Astronautics and Aeronatics, vol 163, American Institute of Aeronautics and Astronautics, Inc, 1996 [35] O'Brien, Andrew J and Gupta, Inder J, "Comparison of output SINR and receiver C/N\_0 for GNSS adaptive antennas," IEEE Transactions on Aerospace and Electronic Systems, vol 45, 2009 [36] Rougerie, S{\'e}bastien and Carrie, Guillaume and Vincent, Fran{\c{c}}ois and Ries, Lionel and Monnerat, Michel, "A new multipath mitigation method for GNSS 111 receivers based on an antenna array," International Journal of Navigation and Observation , 2012 [37] Arribas, Javier and Closas, Pau and Fern{\'a}ndez-Prades, Carles and Cuntz, Manuel and Meurer, Michael and Konovaltsev, Andriy, "Advances in the theory and implementation of GNSS antenna array receivers," Microwave and Millimeter Wave Circuits and Systems: Emerging Design, Technologies, and Applications, pp 227-273, 2013 [38] Sahmoudi, Mohamed and Amin, Moeness G, "Optimal robust beamforming for interference and multipath mitigation in GNSS arrays," Acoustics, Speech and Signal Processing, 2007 ICASSP 2007 IEEE International Conference on, vol 3, pp III-693, 2007 [39] B a M P a G L a G B Widrow, "Adaptive antenna systems," in Proceedings of the IEEE, 1967 [40] R Compton, "The power-inversion adaptive array: Concept and performance," IEEE Transactions on Aerospace and Electronic Systems, pp 803 814, 1979 [41] Dabak, {\"O}mer Can and Erdem, Fatih and S{\"o}nmez, Tolga and Alatan, Lale and Ko{\c{c}}, S Sencer, "Interference suppression in a GPS receiver with element array design and implementation of beamforming algorithms," in Position, Location and Navigation Symposium (PLANS), 2016 IEEE/ION, 2016 [42] Blas, Jose M and Blazquez, Raul and Alonso, Jose I, "A low cost GPS adaptive antenna array," in Institute of Navigation, National Technical Meeting'Vision 2010: Present and Future', San Diego, CA, 1999 [43] S a A D Backen, "GNSS antenna arrays: hardware requirements for algorithm implementation," Lule{\aa} tekniska universitet, 2006 [44] S W a o Kr{\"u}ger, An inexpensive hyperbolic positioning system for tracking wildlife using off-the-shelf hardware, North-West University (South Africa), Potchefstroom Campus, 2017 [45] Schreiber, Rudolf and Bajer, Josef, "Time difference measurement algorithm for TDOA positioning system using RTL-SDR," in Military Technologies (ICMT), 2017 International Conference on, 2017 112 [46] J Liu, B Priyantha, T Hart, Y Jin, W Lee, V Raghunathan, H S Ramos and Q Wang, "CO-GPS: Energy Efficient GPS Sensing with Cloud Offloading," IEEE Transactions on Mobile Computing, vol 15, pp 1348-1361, 2016 [47] H S Ramos, T Zhang, J Liu, N B Priyantha and A Kansal, "LEAP: a low energy assisted GPS for trajectory-based services," in Proceedings of the 13th international conference on Ubiquitous computing, 2011 [48] Brown, Alison and Gerein, Neil and Taylor, Keith, "Modeling and simulation of GPS using software signal generation and digital signal reconstruction," in Proceedings of ION Technical Meeting, 2000 [49] Corbell, Phillip M, "esign and validation of an accurate GPS signal and receiver truth model for comparing advanced receiver processing techniques," in AIR FORCE INST OF TECH WRIGHT-PATTERSONAFB OH SCHOOL OF ENGINEERING, 2000 [50] Corbell, Phillip M and Miller, Mikel M, "A configurable GPS accumulated I and Q signal component simulator in matlab," in Proceedings of ION NTM, Long Beach, CA, 2001 [51] Akos, Dennis M and Pini, Marco, "Effect of sampling frequency on GNSS receiver performance," Navigation , pp 85 95, 2006 [52] M a L P L Berardo, "On the use of a signal quality index applying at tracking stage level to assist the raim system of a gnss receiver," Sensors, vol 16, p 1029, 2016 [53] Trees Van HL, Detection, Estimation, and Modulation Theory Optimum Array Processing, John Wiley \& Sons, 2002 [54] J W a K K R Betz, "Generalized theory of code tracking with an early-late discriminator part II: Noncoherent processing and numerical results," IEEE Transactions on Aerospace and Electronic Systems, vol 45, 2009 [55] Jang, Jaegyu and Paonni, Matteo and Eissfeller, Bernd, "CW interference effects on tracking performance of GNSS receivers," IEEE Transactions on Aerospace and Electronic Systems, vol 48, pp 243 258, 2012 [56] L a M C a L G Dong, "Implementation and verification of a software-based IF GPS signal simulator," in Proceedings of the US Institute of Navigation NTM (San Diego, Jan 26-28), 2004 113 [57] A a G N Brown, "Advanced gps hybrid simulator architecture," in Proceedings of ION 57th Annual Meeting, 2001 [58] A D a B K Arnold, "Distributed Open Source Software-Defined GPS," in Proceedings of the ASWEC 2015 24th Australasian Software Engineering Conference, 2015 [59] D a G C a D F a B M a F J a B G a B M Borio, "An experimental evaluation of the GNSS jamming threat," in 24th Electrotechnical and Computer Conference ERK, 2015 [60] J B.-Y Tsui, Fundamentals of global positioning system receivers: a software approach, vol 173, John Wiley & Sons, 2005 [61] S Godha, "Performance evaluation of low cost MEMS-based IMU integrated with GPS for land vehicle navigation application," UCGEreport, 2006 [62] E.-H Shin, Accuarcy improvement of low cost INS/GPS for land applications, University of Calgary, 2001 [63] G a E G A a M J T a D F Falco, "Performance analysis of constrained loosely coupled GPS/INS integration solutions," Sensors, vol 12, pp 15983 16007, 2012 [64] P a E P Misra, Global positioning system: Signals, measurements and performance second edition, Massachusetts: Ganga-Jamuna Press, 2006 [65] E a E R Ziglioli, "Extending goGPS for Snapshot Positioning," in FOSS4G-Europe, 2015 [66] S Lannelongue and P Pablos, "Fast acquisition techniques for GPS receivers," in Proceedings of the 54th Annual Meeting of the Institute of Navigation, 1998 [67] D M Truong and T H Ta, "Development of real multi-GNSS positioning solutions and performance analyses," in Advanced Technologies for Communications (ATC), 2013 International Conference on, 2013 [68] W a B M Stempfhuber, "A precise, low-cost RTK GNSS system for UAV applications," International Archives of Photogrammetry, Remote Sensing and Spatial Information Science, vol 38, pp 1-22, 2011 114 [69] D a E A a T J Akos, "A prototyping platform for multi-frequency GNSS receivers," Institute of Navigation, The, 2003 [70] petovello2015can, "Can you list all the properties of the carrier-smoothing filter?," INSIDE GNSS, vol 10, no 4, pp 32-37, 2015 [71] Keshvadi, Hatef and Broumandan, Ali and Lachapelle, G{\'e}rard, "Spatial characterization of GNSS multipath channels," International Journal of Antennas and Propagation, 2012 115 APPENDIX A Correlation output calculation The local replica PRN can be illustrated as: 𝑇0 𝑁−1 𝑀𝑘 −1 𝐶𝑙𝑜𝑐𝑎𝑙 (𝑛) = ∑ 𝑘=0 ∑ 𝐶{𝑘}𝑁 𝛿[𝑛 − 𝑚 − 𝑁𝑀𝑘 ] (A.1) 𝑚=0 𝑗=𝑘 where 𝑁𝑀𝑘 = ∑ 𝑗=0 𝑀𝑗 Applying discrete time Fourier transform on the local replica PRN with ω= 2πfTs: ∞ 𝐶(𝜔) = ∑ 𝑛=−∞ 𝑇0 𝑁−1 𝑇0 𝑁−1 𝑀𝑘 −1 ∑ ∑ 𝐶{𝑘}𝑁 𝛿[𝑛 − 𝑚 − 𝑁𝑀𝑘 ]𝑒 −𝑖𝜔𝑛 𝑘=0 𝑚=0 = ∑ 𝐶{𝑘}𝑁 𝑒 𝑀𝑘 −1 −𝑖𝜔𝑁𝑀𝑘 𝑘=0 𝑇0 𝑁−1 ∞ ∑ 𝑒 −𝑖𝜔𝑚 𝑚=0 = ∑ 𝐶{𝑘}𝑁 𝑒 −𝑖𝜔𝑘𝑁𝑀𝑘 𝑘=0 ∑ 𝛿[𝑛] 𝑛=−∞ 𝑠𝑖𝑛(𝑀𝑘 𝜔/2) −𝑖(𝑀𝑘 +1)𝜔 𝑒 𝑠𝑖𝑛(𝜔/2) (A.2) so it can be illustrated in frequency domain as: 𝑇0 𝑁−1 𝐶(𝑓) = ∑ 𝐶{𝑘}𝑁 𝑒 −𝑖2𝜋𝑓𝑇𝑠 𝑁𝑀𝑘 𝑘=0 𝑠𝑖𝑛(𝑀𝑘 𝜋𝑓𝑇𝑠 ) −𝑖(𝑀 +1)𝜋𝑓𝑇 𝑘 𝑠 𝑒 𝑠𝑖𝑛(𝜋𝑓𝑇𝑠 ) (A.3) and its reverse Fourier transform is: 𝛽𝑟 ⁄2 𝑇0 𝑁−1 𝑀𝑘 −1 ∑ ∑ 𝐶𝑘 𝛿[𝑛 − 𝑚 − 𝑁𝑀𝑘 ] = 𝑇𝑠 𝑘=0 𝑚=0 ∫ 𝐶(𝑓)𝑒 𝑖2𝜋𝑓𝑇𝑠 𝑛 𝑑𝑓 (A.4) −𝛽𝑟 /2 where βr is the receiver front-end bandwidth Replacing (A.4) to (6), it yields to: (A.5) 116 𝑁𝐿 −1 𝑇0 𝑎𝑁−1 𝑀𝑘 −1 𝑅(𝜏) 𝐴𝑠 𝑇𝑠 = ∑ 𝑁𝐿 ∑ 𝐶{𝑘}𝑁 𝛿[𝑛 − 𝑚 − 𝑁𝑀𝑘 − ⌊𝜏 + 𝜃𝑁𝐶𝑂 (𝑘)𝑇𝑐 ⌋ 𝑇𝑠 ] ∫ ∑ 𝑛=0 𝑘=0 𝛽𝑟 ⁄2 𝑁𝐿 −1 = 𝐴𝑠 𝑇𝑠 ∫ 𝑁𝐿 𝐴𝑠 𝑇𝑠 ≅ 𝑁𝐿 𝐴𝑠 𝑇𝑠 = 𝑁𝐿 ∑ 𝐶{𝑘}𝑁 𝛿[𝑛 − 𝑚 − 𝑁𝑀𝑘 − ⌊𝜏 + 𝜃𝑁𝐶𝑂 (𝑘)𝑇𝑐 ⌋ 𝑇𝑠 ]𝑒 𝑖2𝜋𝑓𝑇𝑠 𝑛 𝑑𝑓 ∑ 𝐶(𝑓) ∑ ∫ −𝛽𝑟 /2 𝑘=0 𝑚=0 𝑇0 𝑁−1 𝑠𝑖𝑛(𝜋𝑓𝑀𝑘 𝑇𝑠 ) 𝑖2𝜋𝑓𝑇𝑠 (𝑁𝑀 +⌊𝜏+𝜃𝑁𝐶𝑂 (𝑘)𝑇𝑐 ⌋𝑇𝑠 ) 𝑘 ∑ ( ) 𝐶{𝑘}𝑁 𝑒 ∑ 𝐶{𝑙}𝑁 𝑒 −𝑖2𝜋𝑓𝑇𝑠 𝑙𝑀𝑙 𝑑𝑓 𝑠𝑖𝑛(𝜋𝑓𝑇𝑠 ) 𝑘=0 𝑙=0 𝑖2𝜋𝑓𝑇𝑠 (⌊𝜏+𝜃𝑁𝐶𝑂 (𝑘)𝑇𝑐 ⌋𝑇𝑠 ) |𝐶{𝑘}𝑁 | 𝑒 𝑇0 𝑁−1 𝑠𝑖𝑛(𝜋𝑓𝑀𝑘 𝑇𝑠 ) ∑ ( ) 𝑖2𝜋𝑓𝑇𝑠 ((𝑁𝑀 −𝑁𝑀 )+⌊𝜏+𝜃𝑁𝐶𝑂(𝑘)𝑇𝑐 ⌋𝑇𝑠 ) 𝑑𝑓 𝑠𝑖𝑛(𝜋𝑓𝑇𝑠 ) 𝑘 𝑙 + ∑ 𝐶{𝑘}𝑁 𝐶{𝑙}𝑁 𝑒 𝑘=0 𝛽𝑟 ⁄2 𝑇0 𝑁−1 ∫ −𝛽𝑟 /2 ( 𝑇0 𝑁−1 = 𝐴𝑠 ∑ 𝑁𝐿 𝑘=0 𝐶(𝑓)𝑒 𝑖2𝜋𝑓𝑇𝑠𝑛 𝑑𝑓 −𝛽𝑟 /2 𝑚=0 𝑇0 𝑁−1 𝑀𝑘 −1 𝑛=0 𝛽 − 𝑟 𝛽𝑟 ⁄2 𝑇0 𝑁−1 𝛽𝑟 ⁄2 𝛽𝑟 ⁄2 ∫ −𝛽𝑟 /2 𝑙=0,𝑙≠𝑘 ) sin(𝜋𝑓𝑀𝑘 𝑇𝑠 ) 𝑖2𝜋𝑓𝑇 (⌊𝜏+𝜃 (𝑘)𝑇 ⌋ ) 𝑠 𝑁𝐶𝑂 𝑐 𝑇𝑠 𝑑𝑓 𝑇𝑠 ( ) 𝑒 sin(𝜋𝑓𝑇𝑠 ) for a pseudorandom noise code sequence Ck, k=0, 1, 2, , N, all misaligned chips with k = l have equal probability of having +1 or −1 values Their product sum approaches zero B Error analysis for coherent early minus late DLL Define an error signal at the output of the real-part operator in the coherent early minus late TOA estimator as 𝑒(𝜖) = ℜ{ 𝑅(𝜖 − Δ/2) − 𝑅(𝜖 + Δ/2) 𝑇0 𝑁−1 𝑀𝑘 −1 𝑁𝐿 −1 √𝐶𝑠 + ∑ 𝑤(𝑛) ( ∑ 𝑁𝐿 𝑛=0 ∑ 𝐶{𝑘}𝑁 𝛿[𝑛 − 𝑚 − 𝑘𝑀𝑘 − ⌊𝜖 − Δ/2 + 𝜃𝑁𝐶𝑂 (𝑘)𝑇𝑐 ⌋ 𝑇𝑠 ] 𝑘=0 𝑚=0 𝑇0 𝑁−1 𝑀𝑘 −1 − ∑ 𝑘=0 (B.1 ) ∑ 𝐶{𝑘}𝑁 𝛿[𝑛 − 𝑚 − 𝑘𝑀𝑘 − ⌊𝜖 + Δ/2 + 𝜃𝑁𝐶𝑂 (𝑘)𝑇𝑐 ⌋ 𝑇𝑠 ])} 𝑚=0 where ε = t0 −ts is the code phase misalignment, to is the mean value of the un-smoothed TOA estimate (tu), ts is the smoothed TOA estimate, and ∆ is the Early - Late code space (0