Method for locating a GNSS jamming source, and associated computer program product and locating device

Abstract

A method for locating a GNSS signal-jamming source, including setting two antennas in rotation about a common axis of rotation so as to form N different respective positions corresponding to various angles of rotation, in each of the N positions, using each antenna to acquire a GNSS signal including a payload signal and a jamming signal, computing a phase offset between the acquired jamming signals, and determining a direction of the jamming source using a maximum value of the N computed phase shifts.

Claims

1. A locating method of a GNSS signal jamming source, comprising: setting in rotation two antennas about a common axis of rotation to form N different respective positions corresponding to different angles of rotation; in each of the N positions, acquiring by each antenna a GNSS signal comprising a useful signal and a jamming signal; and computing a phase shift between the jamming signals acquired; and determining a direction of the jamming source using a maximum value of the N calculated phase shifts.

2. The method according to claim 1, wherein said setting in rotation is performed by a rotating carrier, the antennas being stationary with respect to the carrier.

3. The method according to claim 1, wherein said setting in rotation comprises setting the two antennas in a full turn.

4. The method according to claim 1, wherein the GNSS signal acquired at each location comprises K samples of the signal.

5. The method according to claim 4, wherein said computing comprises calculating a complex coefficient of cross-correlation between the samples of the acquired GNSS signals at the corresponding position.

6. The method according to claim 5, wherein said computing comprises determining each phase shift between the acquired jamming signals by the argument of the complex cross-correlation coefficient.

7. The method according to claim 1, wherein said determining comprises determining an azimuth angle of the jamming source in a local coordinate frame associated with the two antennas, the azimuth angle being determined in a plane of rotation of the two antennas.

8. The method according to claim 7, wherein the azimuth angle is determined as the angle of rotation of the antennas in the respective position of the antennas corresponding to the maximum value of the calculated N phase shifts.

9. The method according to claim 7, further comprising determining a direction of the jamming source in a geographic coordinate frame, from the azimuth angle of the jamming source and inertial data characterizing an angular position of the antennas in the geographic coordinate frame.

10. The method according to claim 1, wherein the direction of the jamming source is specified by reiterating the method from a different geographical position of the antennas.

11. The method according to claim 10, wherein the different geographic position is determined along the direction of the jamming source determined in a preceding iteration of the method.

12. A computer program product including software instructions which, when executed by a computer, implement a a locating method of a GNSS signal jamming source, comprising: setting in rotation two antennas about a common axis of rotation to form N different respective positions corresponding to different angles of rotation; in each of the N positions, acquiring by each antenna a GNSS signal comprising a useful signal and a jamming signal; and computing a phase shift between the jamming signals acquired; and determining a direction of the jamming source using a maximum value of the N calculated phase shifts.

13. A device for locating a jamming source for GNSS signals, comprising apparatus suitable for implementing a locating method of a GNSS signal jamming source, comprising: setting in rotation two antennas about a common axis of rotation to form N different respective positions corresponding to different angles of rotation; in each of the N positions, acquiring by each antenna a GNSS signal comprising a useful signal and a jamming signal; and computing a phase shift between the jamming signals acquired; and determining a direction of the jamming source using a maximum value of the N calculated phase shifts.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0028] The features and advantages of the invention will appear upon reading the following description, given only as an example, but not limited to, and making reference to the enclosed drawings, wherein:

[0029] FIG. 1 is a schematic view of a device for locating a jamming source according to the invention;

[0030] FIG. 2 is a flowchart of a detection method according to the invention, the method being implemented by the device shown in FIG. 1; and

[0031] FIGS. 3 and 4 are views illustrating the implementation of at least certain steps of the locating method shown in FIG. 2.

DETAILED DESCRIPTION OF EMBODIMENTS

[0032] FIG. 1 illustrates a device 10 for locating a GNSS signal jamming source 12.

[0033] The jamming source 12 has e.g. any electronic device serving to transmit radio signals, called jamming signals, preventing normal reception, by a GNSS receiver, of GNSS signals coming from a GNSS system 14. More particularly, as is known per se, the GNSS system 14 is formed by a plurality of satellites configured to transmit GNSS signals to the ground. The GNSS receiver receives the signals from at least certain of the satellites of the GNSS 14 system in order to determine the geographical position thereof. The GNSS 14 system is e.g. the GPS (Global Positioning System) system or the GALILEO system, known per se.

[0034] In one embodiment, the jamming source 12 is intended to deliberately impair the proper operation of the GNSS receiver. In another embodiment, the jamming source 12 unintentionally impairs the proper operation of the GNSS receiver.

[0035] The locating device 10 according to the invention serves to localize the jamming source 12. Once localized, the jamming source 12 can be deactivated to re-establish the proper operation of the GNSS receiver.

[0036] With reference to FIG. 1, the locating device 10 comprises an input module 21, a processing module 22 and an output module 23. In certain cases, the locating device 10 further comprises a GNSS receiver serving to determine the position thereof in the absence of jamming signals.

[0037] The input module 21 serves to receive radio signals, in particular GNSS signals, which comprise payload signals coming from the GNSS system 14 and jamming signals coming from the jamming source 12. The input module 21 also serves to transmit the received signals to the processing module 22.

[0038] To receive GNSS signals, the input module 21 comprises an antenna array comprising at least two antennas spaced apart from each other. In the example shown in FIG. 1, two antennas 31 and 32 are illustrated. In a generic case, the antenna array may comprise a number of antennas strictly greater than 2.

[0039] As shown in FIG. 1, the antennas 31, 32 are arranged on a carrier 35 in the same plane and are spaced apart from each other in the plane by a distance d. The carrier 35 is advantageously an aircraft, in particular a drone.

[0040] According to the preferred embodiment of the invention, the antennas 31, 32 are stationary with respect to the carrier 35. In such a case, the carrier 35 has a rotating carrier apt to rotate the plane comprising the antennas 31, 32 about an axis of rotation perpendicular to the plane.

[0041] According to another embodiment, the antennas 31, 32 are mounted on a rotating platform which is apt to rotate with respect to the carrier 35. In such a case, the carrier 35 is configured to move in space along e.g. a substantially rectilinear trajectory or has a stationary carrier.

[0042] The processing module 22 is configured to process the GNSS signals received by the input module 21 in order to determine the direction of the jamming source 12, as will be explained in greater detail thereafter.

[0043] The processing module 22 takes the form e.g. of one or a plurality of software stored in a memory and executable by one or a plurality of processors. In a variant or in addition, the processing module 22 is at least partially in the form of a programmable logic circuit, such as an FPGA (Field-Programmable Gate Array) circuit.

[0044] In certain embodiments, the processing module 22 is further configured to control the operation of the antennas 31, 32 and, if appropriate, of the carrier 35. For example, the processing module 22 is configured to control the rotation of the antennas 31, 32 as explained hereinabove. According to other embodiments, the control of the carrier 35 and in particular the setting in rotation of the antennas 31, 32 are carried out from a dedicated control module embedded in the carrier 35 or remote therefrom. Such a control module may also be a part of the locating device 10.

[0045] In the example shown in FIG. 1, the processing module 22 is embedded in the carrier 35, just like the input module 21. According to another embodiment, the processing module 22 is offset from the carrier 35. In such a case, same is apt to receive the signals received by the input module 21 by any appropriate means.

[0046] The output module 23 is configured for delivering the processing performed by the processing module 22. More particularly, the output module 23 is configured for delivering the direction of by the jamming source 12 determined by the processing module 22. For example, the direction of the jamming source 12 is delivered in the form of a heading angle of the jamming source 12 in a geographical coordinate frame the axes of which are e.g. formed by the North, East and vertical directions. According to another embodiment, the direction of the jamming source 12 is delivered in the form of an angle between the direction of movement of the carrier 35 and the direction of the jamming source 12. In the first case, thereof is thus an absolute direction of the jamming source 12 and in the second case, a relative direction.

[0047] The output module 23 is e.g. suitable for supplying the absolute and/or relative direction of the jamming source 12 to an operator and/or to any other system that can be used e.g. to control the carrier 35, such as the control module mentioned hereinabove.

[0048] Finally, just like the processing module 22, the output module 23 can be carried on board the carrier 35 or then offset therefrom.

[0049] The locating device 10 serves to implement the method of locating 100 according to the invention, which will henceforth be explained with reference to FIG. 2 which shows a flowchart of the steps thereof.

[0050] During an initial step 110, the antennas 31, 32 are rotated about the axis of rotation to form N different respective positions corresponding to different angles of rotation.

[0051] More particularly, in the example of FIG. 3 illustrating a coordinate frame (X.sub.Ant, Y.sub.Ant, Z.sub.Ant) linked to the array of antennas, the antenna 31 is placed at the center of the coordinate frame and the antenna 32 is initially placed at a distance d from the antenna 32 along the axis OY.sub.Ant. The plane (X.sub.Ant, Y.sub.Ant) thereby corresponds to the plane of rotation of the antennas 31, 32 and the axis OZ.sub.Ant corresponds to the axis of rotation of the antennas 31, 32.

[0052] Advantageously, during the step 110, a complete revolution is performed about the axis OZ.sub.Ant.

[0053] In each respective position of the antennas 31, 32 during the rotation thereof, the line connecting the centers of the two antennas forms an angle .sub.Ant with respect to the axis OY.sub.Ant. The angle .sub.Ant thus defines each respective position of the antennas 31, 32 during the rotation thereof, called the angle of rotation. Given the initial position of the antenna 32, the angle .sub.Ant varies from 0 to 360 during a full revolution.

[0054] The following step 120 is implemented in parallel with step 110.

[0055] During the step 120, in each position, each antenna 31, 32 acquires a GNSS signal comprising, as explained hereinabove, a payload signal and a jamming signal.

[0056] More particularly, each GNSS signal is acquired in the form of K samples.

[0057] Thereby, by noting s.sub.1(k) the sample k acquired by the antenna 31 and s.sub.2(k) the sample k acquired by antenna 32 in a given position, the samples can be written in the following form:

[00001]$\begin{array}{cc}\begin{array}{c}{s}_1\left(k\right)=s_{1,\mathrm{GNSS}}\left(k\right)+s_{1,B}\left(k\right)\\ s_2\left(k\right)=s_{2,\mathrm{GNSS}}\left(k\right)+s_{2,B}\left(k\right)\end{array}\},& k=1..K\end{array}$

[0058] where S.sub.p,GNSS(k) refers to the payload signal and S.sub.p,B(k) to the jamming signal of the sample k from the corresponding antenna p (p=1, 2).

[0059] During the same step, the processing module 22 then determines a phase shift .sub.est between the jamming signals acquired in the corresponding position.

[0060] To this end, the processing module 22 first calculates the complex cross-correlation coefficient R.sub.xx of the acquired samples, according to the following expression:

[00002]$R_{\mathrm{xx}}\left(s_1,s_2\right)=\frac{1}{K}\underset{k=1}{\overset{K}{.Math.}}{s}_1\left(k\right).Math.s_2^{*}\left(k\right)$

[0061] where (.) * is the operator of complex conjugation.

[0062] The cross-correlation coefficient R.sub.xx is a complex number, i.e. a number with a real part Re(R.sub.xx) and an imaginary part Im(R.sub.xx).

[0063] The phase shift .sub.est between the two jamming signals received on the two antennas is then given by the angle (or argument) of the complex number R.sub.xx, i.e.:

.sub.est=atan 2(Im(R.sub.xx),Re(R.sub.xx)).

[0064] It is thus clear that during the step 120, a phase shift value .sub.est is calculated for each of the N positions defined by the angle of rotation .sub.Ant.

[0065] During the next step 130, the processing module 22 determines the relative direction of the jamming source 12 by using a maximum value of the N calculated phase shifts .sub.est

[0066] More particularly, during the step 130, the processing module 22 determines an azimuth angle AZ.sub.Ant of the jamming source 12 in the plane (X.sub.Ant, Y.sub.Ant). According to the invention, the azimuth angle corresponds to the maximum value of all the phase shifts .sub.est between the two jamming signals determined during the preceding step.

[0067] More particularly, it is clear that the phase difference .sub.est between the two jamming signals received by the two antennas 31, 32 is related to the azimuth Az.sub.Ant and to the site of Si.sub.Ant the jamming source 12 by the following relation:

[00003]${}_{\mathrm{est}}=\frac{2}{}+b_{}$

[0068] where corresponds to the path difference between the antenna 32 and the antenna 31 as shown in FIG. 3, according to which:

[00004]$=d.Math.\cos \left({\mathrm{Az}}_{\mathrm{Ant}}-{}_{\mathrm{Ant}}\right).Math.\sin \left({\mathrm{Si}}_{\mathrm{Ant}}\right)$

[0069] and where corresponds to the wavelength of the jamming signal, b.sub. to a phase shift due to the defect of the antennas and analog channels of the electronic components, and Si.sub.Ant is an elevation angle calculated with respect to the axis OZ.sub.Ant.

[0070] The phase shift .sub.est is thus written as:

[00005]${}_{\mathrm{est}}=\frac{2}{}d.Math.\cos \left({\mathrm{Az}}_{\mathrm{Ant}}-{}_{\mathrm{Ant}}\right).Math.\sin \left({\mathrm{Si}}_{\mathrm{Ant}}\right)+b_{}$ [0071] Since the elevation angle Si.sub.Ant remains constant during the rotation of the antennas, the last relation can be written as:

[00006]${}_{\mathrm{est}}C.Math.\cos \left({\mathrm{Az}}_{\mathrm{Ant}}-{}_{\mathrm{Ant}}\right)+b_{},$

[0072] where C is a constant value.

[0073] In other words, the phase shift .sub.est has a sinusoidal curve. An example of curves .sub.est is shown in FIG. 4. In particular, the FIG. 4 shows in the left-hand part thereof a sinusoidal curve of .sub.est for the elevation value Si.sub.Ant=90 and in the right-hand part thereof a sinusoidal curve of .sub.est for the elevation value Si.sub.Ant=30. In both cases, it is considered that d= and b.sub.=30.

[0074] As the two examples show, it is clear that the phase shift .sub.est reaches the maximum value thereof when Az.sub.Ant=.sub.Ant.

[0075] Thereby, during the step 130, the processing module 22 analyzes all the pairs {.sub.Ant(k), .sub.est(k)}.sub.k=1 . . . N, (.sub.Ant.sub.(N).sub.Ant.sub.(1))2 obtained during the rotation during the previous step and obtains the arrival azimuth estimated Az.sub.Ant.sup.est by:

[00007]$\{\begin{array}{c}{\mathrm{Az}}_{\mathrm{Ant}}^{\mathrm{est}}=\left(m\right)\\ m=\underset{k}{\arg \max }\left({\left\{{}_{\mathrm{est}}\left(k\right)\right\}}_{k=1..N,\left({}_{\mathrm{Ant}}\left(N\right)-{}_{\mathrm{Ant}}\left(1\right)\right)2}\right)\end{array}.$

[0076] In order to improve the accuracy of the determination of the index m of the position of the maximum of the function, in certain embodiments, the processing module 22 determines the intersections of the function .sub.est with a median value (straight line Lm in FIG. 4) which is halfway between the maximum (straight line Lmax in FIG. 4) and the minimum (straight line Lmin in FIG. 4). The processing module 22 then determines the maximum of the function .sub.est which is located at the center of the obtained two intersections framing the first estimated position of the maximum.

[0077] During the next step 140, if need be, the processing module 22 determines the direction of the jamming source 12 in the geographical coordinate frame.

[0078] To this end, the processing module 140 uses e.g. inertial data characterizing the angular position of the carrier 35 with respect to the geographical coordinate frame.

[0079] For example, the processing module 140 can associate an angular position of the carrier 35 with each phase shift value .sub.est measured during the step 120 and can then determine the angular position of the carrier 35 corresponding to the maximum value of the phase shifts .sub.est.

[0080] The direction of the jamming source 12 in the geographical coordinate frame (such as a heading) can then be obtained by transforming, in the geographical coordinate frame, the azimuth and elevation angles Az.sub.Ant Si.sub.Ant determined during the preceding step.

[0081] In the preferred embodiment of the invention, at least steps 110 to 130 and advantageously step 140 are repeated to specify the direction of the jamming source 12.

[0082] For example, steps 110 to 140 can be reiterated several times from different positions of the carrier 35 and then the direction of jamming source 12 is specified by cross-checking the results obtained during the different iterations.

[0083] According to another embodiment, only steps 110 to 130 are reiterated several times. In such case, for each subsequent iteration, the carrier 35 is directed along the direction of the jamming source 12 obtained during the preceding iteration. It is thus clear that in such case, only the relative direction of the jamming source 12 with respect to the carrier 35 is needed. The advantage of such solution is that even if the results of the first iterations are approximate, the carrier 35 will always end up converging in the right direction and the closer the carrier is, the more accurate the results will be.