sensors
Letter Method of Calculating Desynchronization of DVB-T Transmitters Working in SFN for PCL Applications
Karol Klincewicz * and Piotr Samczy ´nski Institute of Electronic Systems, Warsaw University of Technology, 00-665 Warsaw, Poland; [email protected] * Correspondence: [email protected]
Received: 8 September 2020; Accepted: 7 October 2020; Published: 12 October 2020
Abstract: This paper presents a novel method of calculating desynchronization between transmitters working in a single frequency digital video broadcasting-terrestrial (DVB-T) network. The described method can be a useful tool for enhancing passive radar operations and improving passive coherent location (PCL) sensors to correct their measurements of target localization. The paper presents the problem of localizing DVB-T transmitters utilized by passive radars, and proposes a novel method based on Time Difference of Arrival (TDoA) techniques to solve the problem. The proposed technique has been validated using real signals collected by a PCL sensor receiver. The details of the experiment and extensive result analysis are also contained in this article.
Keywords: passive radar; TDoA; SFN; DVB-T
1. Introduction The idea of using non-cooperative transmitters as an illuminator of opportunity (IO) for a passive radar receiver is not a new one. The history of passive radar goes back to the 1930s, when Robert Watson-Watt made the first demonstration of this technology utilizing a BBC radio transmitter as an illuminator of opportunity [1]. The technology was further developed during World War II, when the Germans built a passive radar named Klein Heidelberg which used Chain Home British radars as illuminators of opportunity [1]. After the Second World War, passive radar had no major research interest until the 1990s due to the limited computing power available. Nowadays, after nearly three decades, passive radar technology is in a stage of maturity. It provides many applications, from air surveillance [2–4] to target imaging [5,6], where different kinds of IOs are used, from commercial transmitters (e.g., FM radio, DAB, DVB-T, GSM, among others) [7] to utilize other radars as IOs [8]. Recently, DVB-T transmitters of opportunity became one of the most frequently used as IOs for passive coherent localization (PCL) systems [7,9–11], as they are characterized by relatively wide bandwidth (8 MHz) and relatively high power transmitters, which allows for medium range operation up to 100km with a bistatic range resolution of ca. 36 m [7]. DVB-T can operate in multi-frequency network (MFN) or single frequency network (SFN) configurations. One of the problems of an SFN in passive radars is ghost targets. This phenomenon is a result of the number of signal copies acquired by a single receiver. There are several scientific works in the literature that focus on problems of passive localization systems utilizing SFNs [9,12,13]. A common assumption in most of these works is that SFN transmitters are time-synchronized. In reality, this condition is almost never fulfilled. Most commercial SFN transmitters experience a time shift between each other by design [14]. This fact seems unexplored by passive radar studies, therefore no attempts to cope with it were made. This is unfortunate because SFN desynchronization may introduce major errors in passive radar target localization. The authors of this article recognize a gap in passive radar signal processing studies, and propose a novel method of
Sensors 2020, 20, 5776; doi:10.3390/s20205776 www.mdpi.com/journal/sensors SensorsSensors 20202020,, 2020,, 5776x FOR PEER REVIEW 22 ofof 1212
and propose a novel method of calculating the time shift between transmitters to compensate for this calculating the time shift between transmitters to compensate for this desynchronization, which might desynchronization, which might be used in passive radar signal processing to correct target be used in passive radar signal processing to correct target localization error. localization error. In previous publications [15], the authors described a low-cost system designed and developed to In previous publications [15], the authors described a low-cost system designed and developed localize any source of a signal (an emitter) which might be used for the localization of transmitters of to localize any source of a signal (an emitter) which might be used for the localization of transmitters opportunity within the area of passive radar operations. Based on the results of this work, further of opportunity within the area of passive radar operations. Based on the results of this work, further analysis was made. While positioning a chosen DVB-T transmitter, additional detections occurred analysis was made. While positioning a chosen DVB-T transmitter, additional detections occurred through using Time Difference of Arrival (TDoA) techniques. It was discovered that there were through using Time Difference of Arrival (TDoA) techniques. It was discovered that there were two two transmitters broadcasting in the SFN, and additional detections were caused both by another transmitters broadcasting in the SFN, and additional detections were caused both by another transmittertransmitter workingworking inin thethe SFN,SFN, andand thethe correlationcorrelation byproductbyproduct ofof twotwo transmitters.transmitters. InIn thisthis article,article, anan analysisanalysis ofof thethe originorigin ofof allall thethe detectionsdetections isis made.made. Furthermore,Furthermore, thethe impactimpact ofof transmittertransmitter desynchronizationdesynchronization onon thesethese detectionsdetections isis explored.explored. Finally, aa novelnovel methodmethod ofof calculatingcalculating thethe desynchronizationdesynchronization betweenbetween transmitterstransmitters working in an SFNSFN basedbased onon detecteddetected TDoATDoA peakpeak positionspositions isis presented.presented. The authors believe that this method is a solid addition to anyany passivepassive radarradar operatingoperating inin an an SFN, SFN, and and will will be a be triggering a triggering factor forfactor further for studiesfurther ofstudies the impact of the of SFNimpact transmitter of SFN desynchronizationtransmitter desynchronization on PCL applications. on PCL Inapplications. the second chapterIn the ofsecond this article, chapter passive of this radar article, geometry passive is described.radar geometry The third is described. chapter provides The third a brief chapter description provides of ana brief SFN. description In the next chapter,of an SFN. the In problem the next of localizingchapter, the transmitters problem of used localizing in passive transmitters radars is formulated.used in passive The nextradars two is chaptersformulated. provide The anext concise two explanationchapters provide of TDoA. a concise In the explanation seventh chapter, of TDoA. the experimentIn the seventh scenario chapter, and the setup experiment are presented. scenario In and the eighthsetup are chapter, presented. the results In the of eighth the experiment chapter, arethe shown.results of The the next expe tworiment chapters are shown. present The the proposednext two methodchaptersand present impact the ofproposed measured method desynchronization and impact of on measured the passive desynchronization radar localization on of the the passive target. Theradar last localization chapter concludes of the target. this article.The last chapter concludes this article. 2. Passive Radar Geometry 2. Passive Radar Geometry Passive radar operates on a different basis from active radar. In traditional active radar, the signal Passive radar operates on a different basis from active radar. In traditional active radar, the used for illumination is generated by the transmitting part of the radar system. Most commonly, signal used for illumination is generated by the transmitting part of the radar system. Most a single antenna is used for both transmission and reception, switching between those two states. commonly, a single antenna is used for both transmission and reception, switching between those This configuration is described as monostatic. In a passive solution, the transmitter, or in this case the two states. This configuration is described as monostatic. In a passive solution, the transmitter, or in IO, works independently of the radar system. This creates a receiver–transmitter (Rx–Tx) pair working this case the IO, works independently of the radar system. This creates a receiver–transmitter in a bistatic configuration. (Rx–Tx) pair working in a bistatic configuration. The independence of the exploited transmitter forces each receiver to gather two signals The independence of the exploited transmitter forces each receiver to gather two signals simultaneously: the reference signal gathered from the directional antenna facing the IO; and the simultaneously: the reference signal gathered from the directional antenna facing the IO; and the measurement signal composed of all echo signals (see Figure1). Both signals are then cross-correlated, measurement signal composed of all echo signals (see Figure 1). Both signals are then providing the bistatic range and velocity of targets. In terms of localization (with 2D geometry cross-correlated, providing the bistatic range and velocity of targets. In terms of localization (with assumption), instead of forming a circle (as in a monostatic configuration), the bistatic pair (Rx–Tx) 2D geometry assumption), instead of forming a circle (as in a monostatic configuration), the bistatic forms an ellipse on which the target may be found (see Figure2a). By adding more receivers and /or pair (Rx–Tx) forms an ellipse on which the target may be found (see Figure 2a). By adding more transmitters, one can form more ellipses. The target is located at the point of their intersection receivers and/or transmitters, one can form more ellipses. The target is located at the point of their (Figure2b). intersection (Figure 2b).
Figure 1. Passive radar principle visualization. Figure 1. Passive radar principle visualization.
SensorsSensors 20202020,, 2020,, 5776x FOR PEER REVIEW 33 ofof 1212
Figure 2. Passive radar geometry: (a) for a single ellipse (bistatic configuration), Rx (receiver) and Figure 2. Passive radar geometry: (a) for a single ellipse (bistatic configuration), Rx (receiver) and Tx Tx (transmitter) are ellipse foci, and the sum of RR and RT distances is constant; (b) by adding more ellipses(transmitter) (multistatic are ellipse configuration), foci, and the one sum can findof R targetR and localization.RT distances is constant; (b) by adding more ellipses (multistatic configuration), one can find target localization. Target localization precision is highly dependent on the accuracy of the transmitter and receiver positionTarget data. localization precision is highly dependent on the accuracy of the transmitter and receiver position data. 3. Single Frequency Network Description 3. SingleAn SFN Frequency is a network Network in which Description all transmitters broadcast the same signal simultaneously. In this article,An the SFN DVB-T is a network SFN is in at which the center all transmitters of interest. br Whenoadcast designing the same a signal DVB-T simultaneously. SFN, there are In many this variablesarticle, the to DVB-T consider. SFN One is at of the which center is the of ineliminableinterest. When problem designing of interference a DVB-T SFN, caused there by are multiple many transmittersvariables to consider. broadcasting One at of the which same is time the inelimin in the sameable frequency. problem of It interference can be experienced caused byby amultiple DVB-T networktransmitters user broadcasting at nearly every at the position. same time To prevent in the same any signal frequency. loss or It intersymbol can be experienced interference, by a DVB-T usesnetwork Orthogonal user at nearly Frequency-Division every position. Multiplexing To prevent any (OFDM) signal with loss Guard or intersymbol Interval (GI), interference, preventing DVB-T data lossuses in Orthogonal transmission. Frequency-Divisi The values of GIon vary Multiplexing between 7 µ(OFDM)s and 224 withµs for Guard an 8 MHz Interval channel. (GI), This preventing method solvesdata loss the in problem transmission. of the multipathThe values for of the GI user,vary thereforebetween no7μs strict and synchronization224μs for an 8 MHz of transmitters channel. This in anmethod SFN, norsolves precise the knowledgeproblem of ofthe its multipath value is needed for th [e16 user,]. However, therefore for no passive strict radarssynchronization using DVB-T of transmitterstransmitters asin illuminatorsan SFN, nor ofprecise opportunity, knowledge the exactof its valuevalue ofis theneeded time [16]. shift betweenHowever, transmitters’ for passive broadcastsradars using is valuableDVB-T transmitters information. Aas slightilluminators time diff oference, opportunity, assuming the the exact PCL radarvalue system of the processingtime shift thebetween broadcast transmitters’ is using broadcasts a single transmitter is valuable as information. a reference, willA slight change time target difference, localization assuming ellipses. the InPCL some radar cases, system this shiftprocessing is significant the broadcast enough is for using PCL processinga single transmitter algorithms as nota reference, to consider will the change given ellipsestarget localization as intersecting ellipses. within In an some acceptable cases, errorthis shift margin, is significant and the measurement enough for will PCL be processing discarded asalgorithms non-conclusive. not to consider the given ellipses as intersecting within an acceptable error margin, and the measurement will be discarded as non-conclusive. 4. Transmitter Localization in Passive Radars 4. TransmitterAs mentioned Localizati in previouson in Passive sections, Radars the precise localization of illuminators of opportunity in passiveAs radarsmentioned is of utmostin previous importance. sections, Any the inaccuracy precise localization is propagated of illuminators onto target localization.of opportunity If the in positionspassive radars of the is ellipses’ of utmost focal importance. points (Rx Any and inaccuracy Tx) are uncertain, is propagated the whole onto ellipse target positionlocalization. may If vary. the Inpositions extreme of cases, the ellipses’ when twofocal or points more (Rx ellipses and crossTx) are on uncertain, sections nearly the whole parallel ellipse to eachposition other, may a slight vary. ellipseIn extreme shift cases, may when result two in a or disproportionately more ellipses cross large on sections target localization nearly parallel error, to each or the other, rejection a slight of theellipse measurement. shift may result in a disproportionately large target localization error, or the rejection of the measurement.There are many databases provided by broadcast infrastructure operators, government agencies, and noncommercialThere are many websites, databases although provided transmitters’ by broadcast positions infrastructure are frequently operators, inaccurate, government out of date, oragencies, simply and missing. noncommercial In order to websites, validate dataalthough acquired transmitters’ from external positions sources, are frequently it is possible inaccurate, to use a TDoA-basedout of date, or algorithm simply missing. [15]. Since In order passive to validate radars commonlydata acquir consisted from of external more than sources, one receiverit is possible [17], into mostuse a cases, TDoA-based no additional algorithm hardware [15]. is Since needed passive to gather radars data commonly for TDoA processing.consist of more than one receiver [17], in most cases, no additional hardware is needed to gather data for TDoA processing.
Sensors 2020, 20, 5776 4 of 12
5.Sensors TDoA 2020 Geometry, 20, x FOR PEER REVIEW 4 of 12 TDoA is a technique used in multilateration localization. Multiple receiver stations gather signals 5. TDoA Geometry from a single emitter. Considering two-dimensional space, the distance between the nth station position Rxn (xTDoARxn, yRxn is )a and technique the emitter used position in multilateration (xTx, yTx) can localization. be described Multiple as receiver stations gather signals from a single emitter. Consideringq two-dimensional space, the distance between the nth 2 2 station position Rxn (xRxn, yRxn) andDn =the emitter(xRxn positionxTx) + ((xyRxnTx, yTx)y Txcan) be described as (1) − − 𝐷 𝑥 𝑥 𝑦 𝑦 (1) Measuring the TDoA between two stations gives information about the difference of distances betweenMeasuring the emitter the andTDoA corresponding between two stations: stations gives information about the difference of distances between the emitter and corresponding stations: Dn Dm t = 𝐷− 𝐷, (2) 𝑡 c , (2) 𝑐 where c is the speed of light. where c is the speed of light. The TDoA from single-pair receivers forms a hyperbola on which the source of emission is located The TDoA from single-pair receivers forms a hyperbola on which the source of emission is (see Figure3a). By adding more receivers, a set of non-linear equations is created. The solution of located (see Figure 3a). By adding more receivers, a set of non-linear equations is created. The this set is all the possible locations of the emitter. The problem of solving this equation is well known, solution of this set is all the possible locations of the emitter. The problem of solving this equation is and multiple approaches have been proposed [18]. well known, and multiple approaches have been proposed [18].
Figure 3. TDoA (Time Difference of Arrival) geometry: (a) for single hyperbola pair of receivers (Rx1Figure and 3. Rx2), TDoA Tx (Time are hyperbola Difference foci of and Arrival) the di geometry:fference of R(a1) andfor single R2 distances hyperbola is constant; pair of (receiversb) by adding (Rx1
moreand Rx2), hyperbolae Tx are onehyperbola can find fo theci and localization the difference of the sourceof R1 and of emission.R2 distances is constant; (b) by adding more hyperbolae one can find the localization of the source of emission. For any number of receivers, N, greater than one, it is possible to create a set of N-1 equations and form N-1For any hyperbolae number (Figure of receivers,3b). N, greater than one, it is possible to create a set of N-1 equations and Theform exact N-1 valuehyperbolae of the TDoA(Figure between 3b). two stations (Rxn and Rxm) is calculated by cross-correlating the acquiredThe exact continuous value of signals the TDoA (sn and betweensm, respectively). two stations Cross-correlation (Rxn and Rx ism a) mathematicalis calculated tool by usedcross-correlating in signal processing the acquired to analyze continuous similarities signals between (sn and two s signalsm, respectively). in the function Cross-correlation of time delay is (τ ).a Formathematical continuous tool signals, used cross-correlation in signal processing can be writtento analyze as similarities between two signals in the function of time delay (τ). For continuous signals, cross-correlation can be written as Z∞ Cnm(τ) = sn∗ (t∗) sm(t + τ) dt, (3) 𝐶 𝜏 𝑠 𝑡 𝑠 𝑡 𝜏 dt, (3) −∞ wherewhere * denotes* denotes complex complex conjugation. conjugation. The expected The expect resulted of cross-correlationresult of cross-correlation in a bistatic configurationin a bistatic shouldconfiguration consist should of a distinct consist peak of a corresponding distinct peak corresponding to the TDoA value, to thet. TDoA value, t.
6.6. TDoATDoA inin SFNSFN InIn thethe casecase ofof multiplemultiple emittersemitters workingworking inin anan SFN,SFN, thethe resultsresults ofof cross-correlationcross-correlation areare didifferent.fferent. ConsideringConsidering two two SFN SFN fully-synchronized fully-synchronized transmitters transmitters as sourcesas sources of emission,of emission, more more than than one peakone peak will will occur. Since each receiver acquires the same signal from two transmitters simultaneously (see Figure 4), there are four possible peaks on cross-correlation:
Sensors 2020, 20, 5776 5 of 12 occur. Since each receiver acquires the same signal from two transmitters simultaneously (see Figure4), Sensorsthere are2020 four, 20, x possibleFOR PEER peaks REVIEW on cross-correlation: 5 of 12
1.1. CorrelationCorrelation between between Signal Signal 1-1 1-1 and and Signal Signal 1-2. 1-2. This This peak peak is is the the result result of of the the TDoA TDoA for for the the signal signal emittedemitted by by Tx1. Tx1. 2.2. CorrelationCorrelation between between Signal Signal 2-1 2-1 and and Signal Signal 2-2. 2-2. This This peak peak is is the the result result of of the the TDoA TDoA for for the the signal signal emittedemitted by by Tx2. Tx2. 3.3. CorrelationCorrelation between between Signal Signal 1-1 1-1 and and Signal Signal 2-2. 2-2. This This peak peak is is a a ghost ghost detection. detection. No No physical physical interpretation.interpretation. 4.4. CorrelationCorrelation between between Signal Signal 1-2 1-2 and and Signal Signal 2-1. 2-1. This This peak peak is is a a ghost ghost detection. detection. No No physical physical interpretation.interpretation.
Figure 4. SFN (single frequency network) signal propagation. Assuming the geometry of two receivers and two transmitters, there are four different signals to consider in cross-correlation analysis. Figure 4. SFN (single frequency network) signal propagation. Assuming the geometry of two receiversAssuming and that two the transmitters, precise position there ofare transmitters four different (x Tisignals, yTi) is to well consider known, in thecross-correlation calculation of the exactanalysis. time of occurrence for all four peaks is possible. Writing distance between point A(xA, yA) and point B(x , y ) as B B q Assuming that the precise position of transmitters2 (xTi, yTi) is 2well known, the calculation of the RA B = (xA xB) + (yA yB) , (4) exact time of occurrence for all four −peaks is possible.− Writing− distance between point A(xA, yA) and pointand the B(x timesB, yB) ofas appearance for all described peaks are