A receiver of signals likely to be disturbed by an interfering signal, the receiver includes an interference rejection filter comprising a multiplying gain applied respectively to an in-phase channel and an in-quadrature channel of a signal received by the receiver, the value of the multiplying gain being determined by a nonlinear function ki(ρ) of the modulus ρ of the received signal, the function ki(ρ) being determined so as to maximize the signal-to-noise ratio after application of the multiplying gain and demodulation of the received signal.
|
1. A receiver of signals likely to be disturbed by an interfering signal, the receiver comprising an interference rejection filter comprising a multiplying gain applied respectively to an in-phase channel and to an in-quadrature channel of a signal received and digitized by the receiver, the value of the multiplying gain being determined by a nonlinear function ki(ρ) of the modulus ρ of the received signal, the function ki(ρ) being determined so as to maximize the signal-to-noise ratio after application of the multiplying gain and demodulation of the received signal.
16. A method for receiving signals likely to be disturbed by an interfering signal, comprising the steps of:
receiving respectively an in-phase channel and an in-quadrature channel of a signal,
computing a multiplying gain from a nonlinear function ki(ρ) of the modulus ρ of the received signal, the nonlinear function ki(ρ) being determined so as to maximize the signal-to-noise ratio after application of the multiplying gain and demodulation of the received signal,
applying the multiplying gain respectively to the in-phase channel and to the in-quadrature channel of the digitized signal.
2. The receiver of signals according to
3. The receiver of signals according to
4. The receiver of signals according to
5. The receiver of signals according to
6. The receiver of signals according to
7. The receiver of signals according to
8. The receiver of signals according to
9. The receiver of signals according to
10. The receiver of signals according to
11. The receiver of signals according to
12. The receiver of signals according to
13. The receiver of signals according to
14. The receiver of signals according to
15. The receiver of signals according to
17. The method for receiving signals according to
18. The method for receiving signals according to
19. The method for receiving signals according to
|
This application claims priority to foreign French patent application No. FR 1701348, filed on Dec. 21, 2017, the disclosure of which is incorporated by reference in its entirety.
The invention relates to the field of receivers of signals which are likely to be disturbed by an interfering signal. The invention relates in particular to the field of the receivers of satellite radionavigation signals, or GNSS receivers. It relates to a receiver comprising an interference rejection filter.
A general objective of the invention is to allow a correct demodulation of the signals received by a receiver even when the latter are disturbed by an interfering signal with constant envelope, that is to say with constant amplitude. More specifically, an objective of the invention is to enhance the efficiency, that is to say the signal-to-noise ratio, after demodulation, when the received signal superimposed on an interfering signal with constant envelope.
One possible solution for demodulating a signal disturbed by an interfering signal with constant envelope is described in the European patent application EP 2336800. This solution is based on a quantization of the signal on three levels: +1, 0 and −1. It entails quantizing the analogue signal specifically on the extrema and the 0 crossings of the signal to maximize the sensitivity relative to the useful signal without increasing the power of the noise. This solution presents the drawback of requiring a very precise locking of the automatic gain control in order to accurately position the amplitude of the received analogue signal relative to the quantization thresholds of the analogue-digital converter. Moreover, this solution offers sub-optimal efficiency.
Other receivers use a technique consisting in removing, from baseband complex digitized received signal, a component whose modulus is equal to the estimated mean modulus of the received signal and direction is equal to that of the received signal. This technique presents the drawback of requiring an architecture comprising algorithms of CORDIC (coordinate rotation digital computer) type in cascaded arrangement, which are complex to implement.
The invention proposes implementing an interference rejection filter which essentially comprises the application of a nonlinear function in the form of a multiplying gain, to the in-phase and in-quadrature channels of the received and digitized signal. The proposed solution offers reduced implementation complexity and optimal efficiency.
The subject of the invention is a receiver of signals likely to be disturbed by an interfering signal, the receiver comprising an interference rejection filter comprising a multiplying gain applied respectively to an in-phase channel I and an in-quadrature channel Q of a signal received and digitized by the receiver, the value of the multiplying gain being determined by a nonlinear function Ki(ρ) of the modulus ρ of the received signal, the function Ki(ρ) being determined so as to maximize the signal-to-noise ratio after application of the multiplying gain and demodulation of the received signal.
According to a particular aspect of the invention, the nonlinear function Ki(ρ) depends on the probability density f(ρ) of the received signal, a function of the modulus ρ of the signal.
According to a particular aspect of the invention, the nonlinear function Ki(ρ) depends on the ratio between the derivative f′(ρ) of the probability density of the received signal as a function of the modulus and the probability density f(ρ) of the received signal, divided by the modulus ρ of the received signal.
According to a particular aspect of the invention, the interference rejection filter comprises a unit for computing the modulus of the received signal from the in-phase channel and from the in-quadrature channel of the received signal, the multiplying gain being implemented by a memory in which are saved a plurality of values of the nonlinear function Ki(ρ), the memory having a first input linked to the output of the modulus computation unit.
According to a particular aspect of the invention, the multiplying gain is implemented by a memory in which are saved a plurality of values of the nonlinear function Ki(ρ), the memory having a first input intended to receive the absolute value of the in-phase channel of the signal, a second input intended to receive the absolute value of the in-quadrature channel of the signal.
According to a particular aspect of the invention, the nonlinear function Ki(ρ) depends on an a priori estimate of a ratio J/N between the power of an interfering signal and the power of a thermal noise.
According to a particular embodiment, the receiver of signals according to the invention further comprises an estimator of a ratio J/N between the power of an interfering signal and the power of a thermal noise, the nonlinear function Ki(ρ) depending on said ratio J/N.
According to a particular embodiment of the invention, the estimator of a ratio between the power of an interfering signal and the power of a thermal noise comprises a unit for computing the fourth order moment of the modulus of the signal.
According to a particular embodiment of the invention, the unit for computing the fourth order moment comprises an adder for summing the in-phase channel of the signal, raised to the power of 2, and the in-quadrature channel of the signal, raised to the power of 2, a multiplier for raising to the power of 2 the output of the adder and an averaging filter.
According to a particular embodiment of the invention, the estimator of a ratio between the power of an interfering signal and the power of a thermal noise further comprises a member for computing the ratio between the fourth order moment of the signal and the second order moment of the signal raised to the power of 2.
According to a particular embodiment, the receiver of signals according to the invention comprises at least one automatic gain control device configured to control the amplitude of the received signal in accordance with a power setpoint, the second order moment of the signal being taken to be equal to the power setpoint of the automatic gain control.
According to a particular aspect of the invention, the probability density f(ρ) of the modulus of the received signal is determined, for a plurality of values of the ratio J/N between the power of an interfering signal and the power of a thermal noise, from the value of the amplitude of the interfering signal and from the standard deviation of the thermal noise.
According to a particular aspect of the invention, the memory has an additional input linked to the output of the estimator of the ratio J/N.
According to a particular aspect of the invention, the interfering signal has a constant envelope.
According to a particular aspect of the invention, the received signals are satellite radionavigation signals.
Also a subject of the invention is a method for receiving signals likely to be disturbed by an interfering signal, comprising the steps of:
According to a particular embodiment, the method for receiving signals according to the invention further comprises a step of estimation of a ratio J/N between the power of an interfering signal and the power of a thermal noise, the nonlinear function Ki(ρ) depending on said ratio J/N.
Other features and advantages of the present invention will become more apparent on reading the following description in relation to the attached drawings which represent:
The invention consists in applying an anti-interference filter to the signal received by a receiver, on the in-phase channel and the in-quadrature channel of the signal, after digitization of the signal.
The received analogue signal is first of all transposed into baseband on an in-phase channel I and an in-quadrature channel Q then digitized by means of an analogue-digital converter.
According to one embodiment of the invention, the receiver also comprises an automatic gain control device, on each channel I, Q, which makes it possible to ensure a constant average power of the signal on each channel. In another embodiment of the invention, the automatic gain control is not present.
The anti-interference filter consists in the application of a nonlinear function to the complex signal represented by the channels I, Q.
The theory on which the invention is based is described first of all.
In the presence of non-Gaussian noise, it is possible to enhance the signal-to-noise ratio of the signal after demodulation by applying a nonlinear function f to the signal upstream of the demodulation.
In the presence of an interfering signal with constant envelope, the noise affecting the received signal is non-Gaussian, so this theory can therefore be applied.
Ioptimal and Qoptimal denote the in-phase (respectively in-quadrature) signals after application of the anti-interference filter. The signal-to-noise ratio measured on these signals is optimal when a nonlinear function is applied to the signal as defined by the following relationships:
Ioptimal(I,Q)=−∂P/∂I/P(I,Q) (1)
Qoptimal(I,Q)=−∂P/∂Q/P(I,Q) (2)
P(I,Q) is the probability density of the signal according to the two dimensions I and Q.
∂P/∂I is the partial derivative of P relative to I.
∂P/∂Q is the partial derivative of P relative to Q.
The relationships (1) and (2) are demonstrated in appendix 1 at the end of the description.
It is assumed that the probability density follows a symmetry of revolution in the plane (I, Q) relative to the point (0, 0). This hypothesis is legitimate insofar as the channels I and Q play an equivalent, interchangeable role, linked to the unknowns of the clock of the receiver at the moment of the transition to baseband.
It is assumed that P(I,Q)=f(ρ), with ρ the modulus of the complex signal whose real part corresponds to the in-phase channel I and whose imaginary part corresponds to the in-quadrature channel Q.
Moreover, the following apply: ∂P/∂I=∂P/∂ρ.∂ρ/∂I+∂P/∂θ.∂θ/∂I.
In polar notation, the following applies: I+jQ=ρ exp(jθ)
Knowing that ∂P/∂θ=0, ∂P/∂ρ=f′(ρ) and ∂ρ/∂I=I/ρ, the relationship (1) becomes:
Ioptimal(I,Q)=−∂P/∂I/P(I,Q)=−∂ρ/∂I.f′(ρ)/f(ρ)=−I/ρ.f′(ρ)/f(ρ)=I.K(ρ).
Likewise, the relationship (2) becomes:
Qoptimal(I,Q)=−∂P/∂Q/P(I,Q)=−∂ρ/∂Q.f′(ρ)/f(ρ)=−Q/ρ.f′(ρ)/f(ρ)=Q.K(ρ)
Thus, the nonlinear function to be applied to the signal is:
K(ρ)=−1/ρ.f′(ρ)/f(ρ) (3)
To compute K(ρ), it is first necessary to determine the probability density f(ρ). The probability density of a signal composed of an interfering signal with constant amplitude A and a Gaussian noise of variance σ2 is given by the following formula:
The probability density f(ρ) therefore depends only on the ratio J/N=A2/σ2 between the power of the interfering signal and the power of the Gaussian noise.
From the principles described above, the invention proposes a receiver 200 comprising an interference rejection filter which consists in applying, to each of the channels I and Q of the digital signal, a gain K(ρ) dependent on the modulus of the signal and on the ratio J/N.
The receiver 200 thus comprises two channels for receiving the signal in phase land in quadrature Q resulting from the transposition of the baseband analogue signal (not represented in
The receiver 200 comprises in particular, on each of the channels I, Q, an amplifier 201, 202, an analogue-digital converter 210, 220 and an automatic gain control device 211, 221 for guaranteeing a power of the signal that is constant on each channel before digitization of the signal.
The receiver 200 also comprises an anti-interference filter, comprising a multiplying gain Ki(ρ) 250 which receives as input, on the one hand, the modulus ρ of the signal computed by a computation unit 240 and, on the other hand, an estimate i of the ratio J/N determined by an estimator 230.
According to a particular embodiment of the invention, the ratio J/N can be estimated by calculating the ratio between the fourth order moment of the signal and the second order moment raised to the power of 2.
To this end, the estimator 230 comprises an adder 231 for summing the signals raised to the power of 2 at the output of the respective automatic gain control devices 211, 221 of the channels I and Q and a multiplier 232 for raising the output of the adder 231 to the power of 2.
The estimator 230 also comprises a filter 233 for producing an average of the output (I2+Q2)2 of the multiplier 232 in order to calculate the fourth order moment. The filter 233 is also configured to produce an average of the output I2+Q2 of the adder 231, then raise the average to the power of 2 in order to calculate the second order moment raised to the power of 2. The estimate of the ratio J/N is then obtained by establishing the ratio between the fourth order moment of the signal and the second order moment raised to the power of 2.
The calculation of the second order moment is in particular necessary in the absence of an automatic gain control device 211, 221 or, more generally, when the power of the received signal cannot be considered constant.
Alternatively, when an automatic gain control device 211, 221 is used, as illustrated in
In this case, the calculation of the second order moment of the signal is not necessary.
The computation unit 240 essentially comprises a unit 243 for calculating the modulus ρ of the signal which can be implemented, for example, by means of a table which encodes the values of the square root of I2+Q2, or by a direct calculation, depending on the receiver complexity constraints.
Optionally, in order to limit the number of bits and therefore the capacity of the table used to implement the calculation of the modulus ρ, the computation unit 240 comprises a member 244 for truncating the output of the adder 231 to limit the number of bits at the input of the unit 243.
In the case where the receiver does not comprise an automatic gain control device 211, 221 or more generally the power of the signal cannot be assumed to be constant, the unit 243 for computing the modulus ρ of the signal comprises a computation of the average power of the signal and a normalization of the modulus ρ by the square root of the computed average power.
The multiplying gain 250 Ki(ρ) can be implemented in the form of a table which stores different values of the nonlinear function K(ρ) as a function of the possible values of the modulus ρ and of the ratio J/N.
The receiver 200 finally comprises two multipliers 251, 252 on each of the channels I, Q of the digitized signal, for multiplying the signal by the gain supplied by the table 250.
Without departing from the context of the invention, any other implementation of the anti-interference filter making it possible to produce the same correction function for the channels I and Q of the signal by a corrective gain computed from the relationship (3), is possible.
In this variant, the computation unit 240 comprises a first unit 241 for computing the absolute value of the in-phase channel I of the signal and a second unit 242 for computing the absolute value of the in-quadrature channel Q of the signal which feed the table 243 for computing the modulus ρ of the signal whose values are directly pre-computed as a function of the values of |I| and |Q|.
In this variant, the table 250 for computing the gain Ki(ρ) comprises three inputs instead of two. The values of the gain Ki(ρ) are directly precomputed as a function of the ratio J/N and of the absolute values of I and Q supplied by the computation units 241,242.
In a particular embodiment of the invention not represented in the figures, the estimator 230 of the ratio J/N is eliminated and an a priori value is set for the ratio J/N, for example a maximum value corresponding to a worst case. This a priori value is used to compute the nonlinear function K(ρ). In this particular embodiment, the table 250 storing the values of the multiplying gain comprises only a single input linked to the output of the unit 243 for computing the modulus ρ. This particular embodiment is applicable to the architectures described in
From the curves of
The table 250 thus contains different precomputed values of the gain Ki(ρ) using the curves of
A gain in performance levels is noted when using the invention from a ratio of J/N greater than 0 dB. It is also seen that the invention makes it possible to greatly improve the efficiency of the receiver for ratios between the power of the interfering signal and the power of the useful signal greater than 0 dB. By using the invention, the signal-to-noise ratio losses are limited to 3 dB whatever the level of the interfering signal.
The invention applies for any type of receiver of signals, in particular for receivers of satellite radionavigation signals.
The receiver according to the invention can be implemented by software or hardware means or by a combination of software and hardware techniques.
For example, the receiver according to the invention or each of the elements of which it is composed, can be implemented by means of a processor which can be a generic processor, a specific processor, an application-specific integrated circuit (also known by the acronym ASIC) or a field-programmable gate array (also known by the acronym FPGA). The invention can use one or more dedicated electronic circuits or a general-purpose circuit. The technique of the invention can be produced on a reprogrammable computation machine (a processor or a microcontroller for example) executing a programme comprising a sequence of instructions, or on a dedicated computation machine (for example a set of logic gates such as an FPGA or an ASIC, or any other hardware module).
According to the implementation of this method, in a step 901, a signal is received respectively on an in-phase channel I and an in-quadrature channel Q, for example after transposition of a received signal into baseband. In a step 902, a ratio J/N is estimated between the power of an interfering signal and the power of a thermal noise as described previously. In a step 903, the multiplying gain Ki(ρ) is computed in the manner described above, then, in a step 904, the multiplying gain is applied respectively to the in-phase channel I and to the in-quadrature channel Q of the signal.
In a particular embodiment of the invention, the step 902 of estimation of the ratio J/N is omitted and a predetermined value of the ratio J/N is a priori set for computing a single curve of the multiplying gain K(ρ) as a function of the modulus ρ of the signal.
Appendix 1
This appendix aims to demonstrate the principle given by the relationships (1) and (2) whereby the signal-to-noise ratio measured on the signals Ioptimal and Qoptimal is optimal when a nonlinear function is applied to the signal I, Q as defined by the following relationships, in which P(I,Q) is the probability density of the signal in two dimensions:
Ioptimal(I,Q)=−∂P/∂I/P(I,Q) (1)
Qoptimal(I,Q)=−∂P/∂Q/P(I,Q) (2)
To demonstrate this principle, X is considered to be a non-centred random variable, of average m which corresponds to a real signal, that is to say to a single channel I or Q. The odd (or anti-symmetrical) nonlinear function f is sought that is independent of m and which maximizes the signal-to-noise ratio of the variable Y=f(X).
PX,m(x) denotes the probability density of the variable X as a function of x and PY,m(y) denotes the probability density of Y as a function of y.
The following hypotheses are also set:
∫[−∞,+∞] PX,0(x) dx=1 because PX,0 is a probability density.
∫[−∞,+∞]x PX,0(x) dx=0 which reflects the fact that the average of PX,0 is nil.
∫[−∞,+∞]x2 PX,0(x) dx=σ2 with σ2 the variance of X and σ the standard deviation.
The following properties can be formulated:
PX,m(x)=f′(x).PY,m(y)=dy/dx.PY,m(y)
PX,m(x).dx=PY,m(y).dy
The average of the variable X respects the following properties:
MX(m)=∫[−∞,+∞]x.PX,m(x) dx
MX(m)=∫[−∞,+∞]x.PX,0(x−m) dx
MX(m)=∫[−∞,+∞](x−m).PX,0(x−m) dx+m.∫[−∞,+∞]PX,0(x−m) dx
MX(m)=∫[−∞,+∞]u.PX,0(u) du+m.∫[−∞,+∞]PX,0(u) du (change of variable: u=x−m)
MX(m)=m
The variance of the variable X respects the following properties:
VX(m)=∫[−∞,+∞](x−m)2.PX,m(x) dx
VX(m)=∫[−∞,+∞](x−m)2.PX,0(x−m) dx
VX(m)=∫[−∞,+∞]u2.PX,0(u) du (change of variable: u=x−m)
VX(m)=σ2
The average m characterizes the share of useful signal in the variable X.
The variance σ2 characterizes the share of noise in the variable X.
The average of the variable Y respects the following properties:
MY(m)=∫[−∞,+∞]y.PY,m(x) dy
MY(m)=∫[−∞,+∞]f(x).PX,m(x) dx (change of integration variable)
The variance of the variable Y respects the following properties:
VY(m)=∫[−∞,+∞]y2.PY,m(x) dy
VY(m)=∫[−∞,+∞]f(x)2.PX,m(x) dx (change of integration variable)
From these properties, it is possible to calculate the signal-to-noise ratio of the variables X and Y.
The signal-to-noise ratio for the variable X is: S/NX,m=(MX,m)2/VX,m=m2/σ2
The signal-to-noise ratio for the variable Y is: S/NY,m=(MY,m)2/VY,m
The average of the variable Y is:
MY(m)=∂MY/∂m (m).m=∂MY/∂m (0).m for m«1 and MY(0)=0 (because the function f is odd)
There are also:
∂MY/∂m (m)=∫[−∞,+∞]f(x).∂PX,m/∂m (x) dx
∂MY/∂m (m)=−∫[−∞,+∞]f(x).PX,0′(x−m) dx
∂PX,m/∂m (x)=−PX,0′(x−m) because PX,m(x)=PX,0(x−m)
∂MY/∂m (0)=−∫[−∞,+∞]f(x).PX,0′(x) dx
The variance of the variable Y is:
VY(m)=∫[−∞,+∞]f(x)2.PX,0(x) dx for m«1
For m«1, the signal-to-noise ratio gain R between the variable Y and the variable X is given by the following relationship:
R=(S/NY,m)/(S/NX,m)=[∫[−∞,+∞]f(x).PX,0′(x) dx]2/[∫[−∞,+∞]f(x)2.PX,0(x) dx].σ2.
To optimize the signal-to-noise ratio R, the law PX,0(x) being given, the function f which maximizes the efficiency R is sought.
If it is taken that R=C2/ψ.σ2.
The optimization of R amounts to determining the optimal function f which minimizes the quadratic criterion: ψ=∫[−∞,+∞]f(x)2.PX,0(x) dx subject to the linear constraint: C=∫[−∞,+∞]f(x).PX,0′(x) dx=1
The solution to this optimization problem is: foptimal(x)=λ.PX,0′(x)/PX,0(x)
And, the maximum signal-to-noise ratio is:
Rmax(PX,0)=C2/ψ.σ2=∫[−∞,+∞]PX,0′(x)2/PX,0(x) dx.σ2.
Martin, Nicolas, Perre, Jean-Michel, Bastien, Nicolas
Patent | Priority | Assignee | Title |
Patent | Priority | Assignee | Title |
5291525, | Apr 06 1992 | Freescale Semiconductor, Inc | Symmetrically balanced phase and amplitude base band processor for a quadrature receiver |
6321073, | Jan 31 2000 | Google Technology Holdings LLC | Radiotelephone receiver and method with improved dynamic range and DC offset correction |
6816539, | Jul 20 1998 | Samsung Electronics Company, Limited | Multiple-channel digital receiver for global positioning system |
20020032570, | |||
20070052478, | |||
20070140391, | |||
20070170984, | |||
20080254347, | |||
20090004990, | |||
20090154599, | |||
20110188614, | |||
20120314821, | |||
20130101300, | |||
20130287410, | |||
20140098908, | |||
EP2336800, |
Executed on | Assignor | Assignee | Conveyance | Frame | Reel | Doc |
Dec 17 2018 | Thales | (assignment on the face of the patent) | / | |||
Dec 17 2018 | MARTIN, NICOLAS | Thales | ASSIGNMENT OF ASSIGNORS INTEREST SEE DOCUMENT FOR DETAILS | 047822 | /0081 | |
Dec 17 2018 | PERRE, JEAN-MICHEL | Thales | ASSIGNMENT OF ASSIGNORS INTEREST SEE DOCUMENT FOR DETAILS | 047822 | /0081 | |
Dec 17 2018 | BASTIEN, NICOLAS | Thales | ASSIGNMENT OF ASSIGNORS INTEREST SEE DOCUMENT FOR DETAILS | 047822 | /0081 |
Date | Maintenance Fee Events |
Dec 17 2018 | BIG: Entity status set to Undiscounted (note the period is included in the code). |
May 03 2023 | M1551: Payment of Maintenance Fee, 4th Year, Large Entity. |
Date | Maintenance Schedule |
Nov 19 2022 | 4 years fee payment window open |
May 19 2023 | 6 months grace period start (w surcharge) |
Nov 19 2023 | patent expiry (for year 4) |
Nov 19 2025 | 2 years to revive unintentionally abandoned end. (for year 4) |
Nov 19 2026 | 8 years fee payment window open |
May 19 2027 | 6 months grace period start (w surcharge) |
Nov 19 2027 | patent expiry (for year 8) |
Nov 19 2029 | 2 years to revive unintentionally abandoned end. (for year 8) |
Nov 19 2030 | 12 years fee payment window open |
May 19 2031 | 6 months grace period start (w surcharge) |
Nov 19 2031 | patent expiry (for year 12) |
Nov 19 2033 | 2 years to revive unintentionally abandoned end. (for year 12) |