This application is a National Stage of International patent application PCT/EP2015/071932, filed on Sep. 24, 2015, which claims priority to foreign French patent application No. FR 1402548, filed on Nov. 12, 2014, the disclosures of which are incorporated by reference in their entirety.
The invention relates to a method for detecting signals in a frequency-ambiguous digital receiver. It also relates to a receiver implementing such a method. The invention lies in particular in the field of broadband digital receivers.
In certain applications, owing to the technological limitations of analog-digital coding circuits, digital receivers use sampling frequencies that are much lower than the total reception band, for example a sampling frequency of the order of a GHz for a total band to be processed included in a frequency band of the order of about ten Ghz.
Under these conditions the spectral analysis, performed by digital filtering after signal sampling and coding, gives an ambiguous measurement of the frequency of the signal, owing to spectral aliasings.
To process the spectral aliasings resulting from sub-sampling, these receivers use several reception pathways operating at different sampling frequencies. It is then possible to remove the ambiguity in the frequency measurement on the basis of the various frequency measurements obtained on each of the reception pathways, on the condition that the sampling frequencies are chosen judiciously and that the signal-to-noise ratio is sufficient to allow detection of the signal in each reception pathway. Such a device is in particular described in document EP 1 618 407 A1.
In the cases where the signal-to-noise ratio at the output of the filtering is insufficient, the detection and a fortiori the estimation of the frequency of the signal become impossible. Such is the case for example when the signal to be detected results from a continuous or quasi-continuous emission with low peak power.
A problem to be solved is in particular the detection or the measurement of frequency of a signal with long pulse duration and of low peak power on the basis of a digital receiver having broad frequency band and multiple sub-samplings.
An aim of the invention is in particular to allow the resolution of such a problem.
For this purpose, the subject of the invention is a method for detecting signals in a frequency-ambiguous digital receiver by aliasing of the frequency spectra, said receiver comprising at least two reception pathways, with a sampling frequency for said signals being specific to each pathway, said method comprises at least the following steps:
determining a frequency-wise search band bounded above by a frequency Fmax which is less than or equal to half the smallest sampling frequency from among the sampling frequencies of said pathways, said search band being contained in one and the same ambiguity rank in relation to the sampling frequencies of said pathways;
determining the sampling frequencies of said pathways in such a way that the aliased frequencies corresponding to said pathways are a monotonic function of the true frequency of said signals in said search band;
digitally filtering said signals in said reception pathways in banks of filters of like central frequency which is a multiple of a quantity 1/T and of like width equal to said quantity 1/T, Fechl and Fechm being respectively the sampling frequency of any pathway of order l and of any pathway of order m from among said pathways, Fechl/L=Fechm/M=1/T, L and M being integer numbers;
in said search band, the frequency aliasing of the signal in the pathway of order l giving a response in the filter of rank i and in a pathway of order m a response in the filter of rank k+i or i−k, carrying out a computation of inter-correlation between the signals arising from the filters of rank i of central frequency Fechl/L of the pathway of order l with at least the signals arising from the filters of rank i+k or i−k of the pathway of order m, of central frequency (i+k)Fechm/M or (i−k)Fechm/M;
carrying out the detection of the signals by comparing the power of the signal on output from the inter-correlation computation with respect to a given threshold, the frequency of a detected signal being identified by the knowledge of the ranks of said filters.
The sampling frequencies of said pathways are for example determined in such a way that the disparities between the aliased frequencies corresponding to said pathways are constant in said search band.
The frequency-wise search domain is for example displaced by modifying the sampling frequencies of said pathways in the course of time.
Advantageously, said signals can be signals of low peak power and of long pulse duration.
Advantageously, said method can be applied to an amplitude goniometer comprising several antennas, at least two reception pathways being connected to each antenna. For a given antenna, the inter-correlation computation is for example performed between the signals received on the two reception pathways of the antenna, the result of said computation affording access to the amplitude of the signals which is necessary for estimating their direction of arrival. An inter-correlation computation for example is performed between the received signals of two adjacent antennas.
Advantageously, said method can also be applied to a frequency-ambiguous interferometer, a reception pathway being linked to each antenna of said interferometer. Each reception pathway being associated with a different antenna, the relative phase of a signal is extracted on the basis of the various inter-correlation computations performed between the various pathways taken pairwise.
Said method is for example performed in parallel with a detection processing suitable for detecting pulsed signals.
The subject of the invention is also a digital receiver implementing said method.
Other characteristics and advantages of the invention will become apparent with the aid of the description which follows, given in relation to appended drawings, in which:
In the case of an array of amplitude goniometry antennas comprising several antennas as illustrated by
In the case of an interferometry antenna as illustrated by
In the case of a computational beamforming antenna as illustrated by
For all these devices in particular, it is necessary to solve the problem of the detection and of the measurement in terms of frequency of the signals of long duration and of low peak power on the basis of a broadband digital receiver with multiple sub-sampling. Currently, the sensitivity of broadband receivers, be they analog or digital, is insufficient to detect these signals of long duration and low peak power.
According to the invention, as will be described subsequently, specific means allowing the detection of the signals of low peak power and of long duration by carrying out an inter-correlation of the received signals sampled at different frequencies are integrated into a broadband digital receiver.
In such a receiver, the various frequencies Fech1, Fech2, Fech3, Fech4 are chosen so that the ambiguity in the frequency measurement is removed with a sufficient margin of safety in the presence of noise, thus requiring a sufficient spacing between the various frequencies, typically of the order of a few tens of MHz. These various sampling frequencies must also be chosen so that their lowest common multiple is greater than the total analysis band in respect of the received signals. Moreover, because of Shannon's theorem, the instantaneous bandwidth being limited to a value of less than half the lowest sampling frequency, the sampling frequencies must be chosen to be as high as possible. Finally, it is practical to use a constant spacing δF between the sampling frequencies since this makes it possible in particular to use simple algorithms to remove the distance-wise ambiguities.
These constraints and facilities lead for example to choosing sampling frequencies as follows, M being an integer number:
Fech1=MδF;Fech2=(M+1)δF;Fech3=(M+2)δF;Fech4=(M+3)δF (1)
For example, if the sampling frequency of the coders is of the order of 2 GHz, it is possible to choose, by taking M=52:
For each sampling frequency, this ambiguous frequency, forming an aliased frequency, can be written in the following manner:
Fambi=Ftrue−kambi×Fechi if Ftrue>kambi×Fechi (1)
Fambi=kambi×Fechi−Ftrue if Ftrue<kambi×Fechi (2)
where:
Returning to
To remove the ambiguity and preserve the relative phase between the pathways 1, 2, 3, 4, the same analysis resolution is used on the various pathways. For example, if a frequency resolution of 10 MHz is desired, it is possible to choose:
This mode of detection and estimation of the signal frequency as described hereinabove in relation to
Moreover, the frequency of the input signal of the receiver being unknown, it is not possible to know the way in which the aliasings 51, 52, 53, 54 are performed and therefore the filters in which the signal is present on the various pathways 1, 2, 3, 4. This prohibits a priori the possibilities of recombining of the received signals between the various pathways aimed at increasing the signal-to-noise ratio, so as to allow detection. Such is the case in particular in the presence of long signals of low peak power.
According to the invention, on the basis of a determined and limited frequency domain 61 in which it is desired to seek to detect a continuous emission or a long pulse emission of low peak power, the various sampling frequencies are chosen in such a way that the disparities between aliased frequencies corresponding to the various reception pathways 1, 2, 3, 4 are constant in this determined frequency domain 61.
This is obtained, if:
For any true frequency value, Ftrue, contained in the frequency domain 61 defined by bounds Fmin, Fmax, we then have:
Fechmin is the smallest of the sampling frequencies of the various pathways. In the present example, Fechmin is Fech1. Likewise, Fechmax is the largest of the sampling frequencies and Fechmax is Fech4.
For example, by choosing the increasing direction of variation of the ambiguous frequency, it is possible to choose:
Fechmax≥Fmin/kamb.
For a frequency domain 61 limited to ΔF, it follows that:
ΔF≤Fmax−Fmin, i.e.:
ΔF≤kamb×(Fechmin−Fechmax)+Fechmin/2.
And by using relationship (3), we obtain for four sampling frequencies:
For example, if the low bound of the search domain is fixed at
Fmin=9 GHz,
it is possible to choose, having regard to the maximum accessible sampling frequencies of the order of 2 GHz, by taking as sampling frequency step size δF=40 MHz:
Hence:
And for a frequency search domain bounded by ΔF≤580 MHz, we obtain Fmax=9.58 GHz.
After having fixed the sampling frequencies in accordance with the foregoing, the disparity of the aliased frequencies between the various reception pathways 1, 2, 3, 4 in the search domain 61 is determined in a following step.
In this domain 61, and as illustrated by
According to relationship (3) specifying that:
Fech1/N1=Fech2/N2=Fech3/N3=Fech4/N4=1/T
it follows that:
i2−i1=−kamb.δF.T
i3−i1=−2kamb.δF.T
i4−i1=−3kamb.δF.T.
For example, for the above example, if the spectral analysis resolution is fixed at 10 MHz, corresponding to T=100 ns, for kamb=4 and δF=40 MHz:
i2−i1=−16
i3−i1=−32
i4−i1=−48.
Knowing the rank i1 involved in the first bank of filters 701, the ranks i2, i3, i4 of the filters involved in the other banks 702, 703, 704 are deduced therefrom.
Having thus identified the differences of index of the filters involved 71, 72, 73, 74, a following step consists in performing complex inter-correlations 75, 76, 77 of the signals arising from these filters between the various pathways.
This inter-correlation is performed over a long time, typically of the order of 100 microseconds, in accordance with the type of inter-correlation described for example in patent application FR 1400514.
It makes it possible not only to obtain a sufficient signal-to-noise ratio for detection, but also to extract the phase difference between the reception pathways, for example in the case of an interferometer, or else the amplitude difference between two adjacent antennas, in the case of an amplitude goniometer.
In this process, the output of the filter of index i1 of pathway 1 will be correlated with the output of the filter of index i2 of pathway 2 and for example, for 4 reception pathways, the output of the filter of pathway 1 will also be correlated with the output of the filter i3 of pathway 3 and i4 of pathway 4, as illustrated by the example of
This process is performed for all the indices of filters whose central frequencies are contained in the search domain.
In a following step 78, the results of the inter-correlations are thereafter compared in terms of amplitude with a threshold so as to ensure detection of the signals and identify the indices of the filters corresponding to the signal so as to estimate the frequency thereof. The amplitude and the phase of the inter-correlations are also stored at this level, for example so as to extract the direction of arrival of the signal.
After having thus utilized a search domain 61, it is thus possible to choose a new domain by defining a new set of sampling frequencies.
The invention has been described by way of example in respect of a digital receiver with four pathways. It applies more generally in respect of a receiver comprising at least two different pathways corresponding to different sampling frequencies Fechl, Fechm in banks of filters 701, 702, 703, 704 of like central frequency which is a multiple of a quantity 1/T, where 1/T=Fechl/L=Fechm/L, and of like widths 1/T=Fechl/L=Fechm/L, L and M being integer numbers, T being the duration of the observation window.
Advantageously, the invention can in particular be applied to an interferometer or to a receiver with amplitude goniometry.
The digital signals arising from the coders are thereafter processed in accordance with the method according to the invention, as output from spectral analyses 45, 46, 47, 48. The inter-correlation as described hereinabove is performed in particular between the received signal of the first antenna 21 with the received signals of the other antennas 22, 23, 24, the filter of index being associated with the first reception pathway linked to the first antenna 21. Stated otherwise, each reception pathway, of different sampling frequency, being associated with a different antenna, the relative phase of the signals is extracted on the basis of the various inter-correlation computations performed between the various pathways taken pairwise. The direction of arrival of the signals is obtained with the aid of this phase.
The removals of ambiguity of the frequency measurement and of the interferometry phase measurement are thus performed in a single operation resulting from the inter-correlation products, as is described in relation to
On a given antenna 104, the inter-correlation is performed on two pathways 1, 2 of different sampling frequencies. This inter-correlation computation affords access to the measurement of amplitude of the signal necessary for estimating the direction of arrival. The removal of ambiguity in the frequency measurement is for example performed with the aid of the frequency-aliased signals obtained on the basis of the same signal sampled with the aid of two different frequency pairs. This signal may originate from one and the same antenna 104, or two adjacent antennas 103, 104 receiving it simultaneously. Coverage over 360° is obtained with the aid of 12 reception pathways for a goniometric device with 6 antennas, each antenna being linked to two pathways.
Number | Date | Country | Kind |
---|---|---|---|
14 02548 | Nov 2014 | FR | national |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/EP2015/071932 | 9/24/2015 | WO | 00 |
Publishing Document | Publishing Date | Country | Kind |
---|---|---|---|
WO2016/074841 | 5/19/2016 | WO | A |
Number | Name | Date | Kind |
---|---|---|---|
20020093342 | Ivanov | Jul 2002 | A1 |
20060195301 | Mueller | Aug 2006 | A1 |
20070109213 | Takahashi | May 2007 | A1 |
20110035168 | Lelong | Feb 2011 | A1 |
20140225591 | Knierim | Aug 2014 | A1 |
20160226505 | Auston | Aug 2016 | A1 |
Number | Date | Country |
---|---|---|
2 725 707 | Apr 2014 | EP |
2 897 441 | Aug 2007 | FR |
2 939 988 | Jun 2010 | FR |
3 018 136 | Sep 2015 | FR |
2004097450 | Nov 2004 | WO |
Entry |
---|
M. Yaghoobi et al., “A low-complexity sub-Nyquist sampling system for wideband Radar ESM receivers,” 2014 IEEE International Conference on Acoustics, Speech and Signal Processing, May 2014, pp. 1788-1792, XP032618103. |
McCormick, “A high resolution, near real-time frequency estimator for sub-microsecond pulses,” IEEE International Conference on Systems Engineering, Jan. 1, 1989, pp. 33-38, XP032132227. |
Number | Date | Country | |
---|---|---|---|
20170336450 A1 | Nov 2017 | US |