EP4334742A1 - Methods for computer estimation of a single tone and their application to radar systems - Google Patents
Methods for computer estimation of a single tone and their application to radar systemsInfo
- Publication number
- EP4334742A1 EP4334742A1 EP22726142.7A EP22726142A EP4334742A1 EP 4334742 A1 EP4334742 A1 EP 4334742A1 EP 22726142 A EP22726142 A EP 22726142A EP 4334742 A1 EP4334742 A1 EP 4334742A1
- Authority
- EP
- European Patent Office
- Prior art keywords
- estimate
- dft
- complex
- iteration
- frequency
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 95
- 238000004364 calculation method Methods 0.000 claims abstract description 36
- 238000001228 spectrum Methods 0.000 claims abstract description 13
- 230000004044 response Effects 0.000 claims abstract description 11
- 239000013598 vector Substances 0.000 claims description 54
- 238000012545 processing Methods 0.000 claims description 7
- 238000005070 sampling Methods 0.000 claims description 7
- 238000004891 communication Methods 0.000 claims description 4
- URWAJWIAIPFPJE-YFMIWBNJSA-N sisomycin Chemical compound O1C[C@@](O)(C)[C@H](NC)[C@@H](O)[C@H]1O[C@@H]1[C@@H](O)[C@H](O[C@@H]2[C@@H](CC=C(CN)O2)N)[C@@H](N)C[C@H]1N URWAJWIAIPFPJE-YFMIWBNJSA-N 0.000 claims 2
- 238000003491 array Methods 0.000 abstract description 2
- 238000000194 supercritical-fluid extraction Methods 0.000 description 30
- 238000004422 calculation algorithm Methods 0.000 description 12
- 238000007476 Maximum Likelihood Methods 0.000 description 11
- 238000011156 evaluation Methods 0.000 description 9
- 238000013459 approach Methods 0.000 description 8
- 230000008569 process Effects 0.000 description 8
- 238000010586 diagram Methods 0.000 description 7
- 238000009795 derivation Methods 0.000 description 6
- 238000001514 detection method Methods 0.000 description 4
- 230000003595 spectral effect Effects 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000002592 echocardiography Methods 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 230000000737 periodic effect Effects 0.000 description 3
- 239000000654 additive Substances 0.000 description 2
- 230000000996 additive effect Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 238000004590 computer program Methods 0.000 description 2
- 230000001143 conditioned effect Effects 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 230000002045 lasting effect Effects 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 208000030514 Leukocyte adhesion deficiency type II Diseases 0.000 description 1
- 108010076504 Protein Sorting Signals Proteins 0.000 description 1
- GVBNSPFBYXGREE-CXWAGAITSA-N Visnadin Chemical compound C1=CC(=O)OC2=C1C=CC1=C2[C@@H](OC(C)=O)[C@@H](OC(=O)[C@H](C)CC)C(C)(C)O1 GVBNSPFBYXGREE-CXWAGAITSA-N 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 230000010267 cellular communication Effects 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- 230000010356 wave oscillation Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/35—Details of non-pulse systems
- G01S7/352—Receivers
- G01S7/354—Extracting wanted echo-signals
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/06—Systems determining position data of a target
- G01S13/08—Systems for measuring distance only
- G01S13/32—Systems for measuring distance only using transmission of continuous waves, whether amplitude-, frequency-, or phase-modulated, or unmodulated
- G01S13/34—Systems for measuring distance only using transmission of continuous waves, whether amplitude-, frequency-, or phase-modulated, or unmodulated using transmission of continuous, frequency-modulated waves while heterodyning the received signal, or a signal derived therefrom, with a locally-generated signal related to the contemporaneously transmitted signal
- G01S13/345—Systems for measuring distance only using transmission of continuous waves, whether amplitude-, frequency-, or phase-modulated, or unmodulated using transmission of continuous, frequency-modulated waves while heterodyning the received signal, or a signal derived therefrom, with a locally-generated signal related to the contemporaneously transmitted signal using triangular modulation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/35—Details of non-pulse systems
- G01S7/352—Receivers
- G01S7/356—Receivers involving particularities of FFT processing
Definitions
- the present invention relates to a deterministic method for the detection of a tone and the estimation of its parameters in a colocated multiple-input multiple-output radar system operating in a slowly changing propagation scenario.
- a tone as wellknown, is a wave oscillation carachterized by a certain frenquency, amplitude and phase.
- the estimation of the single tone means astimating its frequency, aplitude and phase.
- Coarse estimation is always based on the maximization of the periodogram of the observed signal, whereas fine estimation can be accomplished in an open loop fashion or through an iterative procedure.
- all the open loop estimators exploit spectral interpolation to infer the frequency residual from the analysis of the fast Fourier transform (FFT) coefficients at the maxima of the associated periodogram and at frequencies adjacent to it as in "B. Quinn, "Estimation of frequency, amplitude, and phase from the DFT of a time series," IEEE Trans. Signal Process., vol. 45, no. 3, pp. 814-817, Mar. 1997.” and “Q. Candan, "Analysis and Further Improvement of Fine Resolution Frequency Estimation Method From Three DFT Samples," IEEE Signal Process.
- FFT fast Fourier transform
- a (
- , xpi A ⁇ Ai)
- w c n is the n-th sample of an additive white Gaussian noise (AWGN) sequence (whose elements have zero mean and variance 2s 2 )
- w r n 3 ⁇ 4 ⁇ w cn )
- N is the overall number of samples
- AWGN additive white Gaussian noise
- the signal models (1) and (2) appear in a number of problems concerning radar systems, wireless communications and biomedical applications.
- MIMO multiple- input multiple output
- TDM time division multiplexing
- TX transmit
- RX receive
- N T elements and N R elements respectively; these allow to radiate orthogonal waveforms from different antennas and to receive distinct replicas of the electromagnetic echoes generated by multiple targets.
- TX and RX arrays are made of distinct antenna elements, placed at different positions. Moreover, TX antennas are close to the RX ones and, in particular, are usually placed on the same shield.
- Orthogonality of the transmitted waveforms is achieved by radiating them through distinct TX antennas over disjoint time intervals.
- Fig. 1 The architecture of the radar system considered in this work is illustrated in Fig. 1.
- VCO voltage controlled oscillator
- TX array radiated by its TX array after power amplification.
- FMCW frequency modulated continuous wave
- the VCO is fed by a periodic ramp generator; this produces a chirp FM signal, whose instantaneous frequency evolves periodically, as illustrated in Fig. 2.
- the parameters T, T R and T 0 represent the chirp interval, the reset time and the pulse period (or pulse repetition interval), respectively, whereas the parameters f 0 and B are the start frequency and the bandwidth, respectively, of the transmitted signal.
- LNA low noise amplifier
- ADC analog-to-digital converter
- the couple of the involved physical TX and RX antennas (namely, the p-th TX antenna and the q-th RX antenna) of the considered bistatic radar is usually replaced by a single virtual antenna of an equivalent monostatic radar.
- the abscissa x v and the ordinate y v of the 17-th virtual antenna element associated with the p-th TX antenna and the q-th RX antenna can be computed as (e.g., see "J. Li and P. Stoica, Eds., MIMO Radar Signal Processing. Hoboken, NJ: J. Wiley & Sons, 2009. ") and
- nL and w rn the n-th sample of the AWGN sequence affecting the received signal, is a Gaussian random variable having zero mean and variance s 2 (assumed to be independent of v).
- the real signal model (10) is identical to that expressed by eq. (2) and can be adopted in all the FMCW radar systems that extract only the in-phase component of the signal captured by each RX antenna; b) some commercial MIMO radar devices provide both the in-phase and quadrature components of the received RF signals (e.g., see "M. A. Richards, Fundamentals of Radar Signal Processing. New York: McGraw-Hill, 2005.”).
- the complex model equivalent to that expressed by eq. (1) must be adopted in place of its real counterpart (10) for any n; here, for any v and Z.
- SFCW stepped frequency continuous wave
- a t are still expressed by eqs. (14) and (16), respectively, is the normalised delay characterizing the Z-th target and observed on the v-th virtual antenna, and the parameters (a ; , have exactly the same meaning as the one illustrated for the received signal model (15).
- the signal observed on each virtual channel can be represented a superposition of L real or complex oscillations; moreover, the value of the parameter L has to be considered unknown.
- Distinct radar echoes can be characterized by substantially different signal-to-noise ratios (SNRs).
- the number N of samples acquired over each virtual channel usually ranges from few hundreds to few thousands.
- Scope of the present invention is to provide for methods for estimation of a single tone, which has a low computation complexity and that can be easily implemented in low computation power such has radars and in the base stations of cellular communication networks.
- the basic idea of the present invention is based on the approximate maximum likelihood estimation of the single tone.
- the present invention proposes a new method that can be declined into different ways according to the reference domain, such as real or complex frequency or time domain. This is the reason for which the same inventive concept is described with the aid of the claims 1 and 10. More in detail, when the single tone represents a frequency, the present invention proposes a method for computer estimation of the frequency, amplitude and phase of the single tone defining a signal acquired by means of a receiver, wherein the tone is real or complex; the method including
- the present invention proposes a method for computer estimation of the delay, amplitude and phase of a single complex tone defining a signal acquired by means of a receiver; the method including
- the receiver can be part of a radar or of any device capable to acquire a single tone signal.
- the present invention can be implemented also in 5G base stations.
- the present invetion aims at providing a total of three approaches about the estimation of a tone and their exploitation in colocated frequency modulated continuous wave (FMCW) and stepped frequency continuous wave (SFCW) multiple-input multiple-output (MIMO) radar systems operating at millimeter waves.
- FMCW colocated frequency modulated continuous wave
- SFCW stepped frequency continuous wave
- MIMO multiple-input multiple-output
- a novel ML-based iterative estimator of a single real and complex tone is presented.
- this estimator is based on: a) expressing the dependence of the ML metric on the frequency residual in an approximate polynomial form through standard approximations of trigonometric functions; b) exploiting the alternating minimization technique for the maximization of this metric. Moreover, its most relevant feature is represented by the fact that it requires the evaluation of spectral coefficients that are not exploited by all the other related estimation methods available in the technical literature.
- Figure 1 shows an example of architecture of a MIMO FMCW radar system
- Figure 2 discloses the time diagram of the instantaneous frequency of the RF signal radiated in a FMCW radar system
- Figure 3 shows an example of architecture of a MIMO SFCW radar system
- Figure 4 discloses the time evolution of the instantaneous frequency of the RF signal radiated in a SFCW radar system
- FIG. 5 shows a block diagram of the SFE#1 algorithm
- Figure 6 shows a block diagram of the SFE#2 algorithm
- Figure 7 shows a block diagram of the CSFE#1 algorithm
- Figure 8 shows a block diagram of the CSFE#2 algorithm
- Figure 9 shows a block diagram of the CSDE#1 algorithm
- Figure 10 shows a block diagram of the CSDE#2 algorithm
- MSE mean square error
- X R ⁇ P) ⁇ M ⁇ X (F) ⁇ , X,(P)A 3 ⁇ X(P) ⁇ , is, up to the scale factor 1/iV, the Fourier transform of the sequence ⁇ x r,n ⁇ ,
- F c F c + dF dft , (40)
- F c is a given coarse estimate of F
- d is a real variable called residual
- residual is the normalized fundamental frequency associated with the N 0 - th order discrete Fourier transform (DFT) of the zero padded version of the vector collecting all the elements of the sequence ⁇ x r,n ⁇
- M is a positive integer (dubbed oversampling factor)
- 0 D is a D —dimensional (column) null vector
- iV 0 M ⁇ N .
- N7 exp(-; ' 4reF)-l’ which is easily derived from its definition (39).
- 7(F) can be computed exactly through its expression (38) or, in an approximate fashion, through a computationally efficient procedure based on the fact that the vector
- X s M X 0 (51) collects iV 0 uniformly spaced samples of the function X(F) (38).
- the k-th element of the vector X 0 (42) can be expressed as where
- F F c + d F DFT (55) and on a novel method for estimating the real residual d, i.e. for accomplishing the fine estimation of F .
- This method is derived as follows. Representing the trial normalized frequency F according to eq. (40) allows us to express the variable f h (35) as
- n (F, C) X 2 in + 2 [c 2 + Cf] +2(C
- SFE #1 single frequency estimator #1
- the new estimate D® of D is computed through eq. (78) (eq. (79)); in the evaluation of the coefficients (a(p), b(p), c(p) ⁇ ( ⁇ b(p) , c(p) ⁇ ) appearing in the RHS of these equations, is assumed. Then, is evaluated.
- the SFE#1 can be employed even if the single tone appearing in the RHS of eq. (27) is replaced by the superposition of L distinct tones (see eq. (3)). In this case, the strongest (i.e., the dominant) tone is detected through the periodogram method (see eq. (82)) and the parameters of this tone are estimated in the presence of both Gaussian noise and the interference due to the remaining tones. Therefore, the estimation accuracy of the SFE#1 is affected by both the amplitudes and the frequencies of the other (L— 1) tones.
- the second step of SFE#2, instead, is the same as that of SFE#1.
- the method for the estimation of the frequency, the amplitude and the phase of a real sequence can be described as follows:
- Steps 2.1.1 - 2.1.2 of Fig. 5 can be replaced by an alternative method including (see Steps 2.2.1 - 2.2.2 of Fig. 6):
- IDFT iV 0 order inverse DFT
- CSDEs complex single delay estimators
- the CSDEs compute an estimate d of the delay residual d by solving eq. (74) in the variable D (57); in this case, the coefficients a(p) and c(p) are still expressed by eqs. (105) and (107), respectively, whereas b(p)A -3 ⁇ 4 ⁇ i * X 2,p ⁇ , (109) where A is the estimate of the complex gain A computed on the basis of eq. (102).
- the quantity X k p(i-q is replaced by where pO -1 - 1 is still defined by eq. (94).
- the method for the estimation of the delay, the amplitude and the phase of a complex sequence (see Steps 2.5.0
- said methods can be employed for angle of arrival (Q) calculation of a point target having range (R) detected through a multiple-input multiple-output stepped frequency continuous wave radar system equipped with a plurality of transmitting (TX) and receiving (RX) antennas, wherein each couple of said TX and RX antennas is replaced with the equivalent virtual antenna; said virtual antennas forming a virtual array composed by one uniform linear array; said angle of arrival being computed as: here, l is the radar wavelength, d is the distance between two adjacent virtual antennas of said uniform linear array, and T is said estimated normalized delay according to said methods.
- This invention can be implemented advantageously in a computer program comprising program code means for performing one or more steps of such method, when such program is run on a computer.
- the patent shall also cover such computer program and the computer-readable medium that comprises a recorded message, such computer-readable medium comprising the program code means for performing one or more steps of such method, when such program is run on a computer, such as a computing device associated to a radar system or to a base station of cellular telecommunication, such as 5G and similar ones.
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Signal Processing (AREA)
- Radar Systems Or Details Thereof (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
Method for estimating the angular coordinates of point targets whose range has been estimated through a MIMO FMCW radar equipped with a plurality of transmitting (TX) and receiving (RX) antennas, wherein each couple of said TX and RX antennas is replaced with the equivalent virtual antenna; said virtual antennas form a virtual array composed by one or multiple horizontal uniform linear arrays (HULAs); said MIMO FMCW radar being arranged to generate real or complex signals in response to a propagation scenario including a plurality of point targets; said method including acquisition of a spectrum of said signal and its first derivatives, acquisition of a list of a predetermined number of discretized ranges, named frequency bins, of said previously detected point targets, and sequential execution for each of said frequency bins; said execution including: estimation of the horizontal spatial frequency and complex amplitude of the most dominant point target, the estimation including: computation of the spectrum its first three derivatives through EFT calculation, checking whether combination of the spectra associated with said Reference HULA, generates an amplitude peak in the resulting spectrum and when the checking is positive execution of the following steps: refinement of said range and complex amplitude of said most dominant point target; cancellation of said most dominant point target, calculation of a residual energy in the considered frequency bin and its comparison with a predetermined threshold, and when a residual energy exceeds is below the predetermined threshold computation of the spatial coordinates, otherwise updating of said acquired list of said frequency bins and return analyze a further frequency bin; finally, generation of an overall image of the propagation scenario.
Description
Methods for Computer Estimation of a Single Tone and Their Application to Radar Systems
Field of the invention
The present invention relates to a deterministic method for the detection of a tone and the estimation of its parameters in a colocated multiple-input multiple-output radar system operating in a slowly changing propagation scenario.
A tone, as wellknown, is a wave oscillation carachterized by a certain frenquency, amplitude and phase. Thus, the estimation of the single tone, means astimating its frequency, aplitude and phase.
State of the art
The problem of estimating the amplitude, phase and frequency of a tone in additive white Gaussian noise (AWGN) has received significant attention for a number of years because of its relevance in various fields, including radar systems and wireless communications. It is well known that the maximum- likelihood (ML) approach to this problem leads to a complicated nonlinear optimization problem. In practice, the most accurate ML-based single tone estimators available in the technical literature achieve approximate maximisation of this metric through a two-step procedure; the first step consists in a coarse search of tone frequency, whereas the second one in a fine estimation generating an estimate of the so called frequency residual (i.e., of the difference between the real frequency and its coarse estimate). Coarse estimation is always based on the maximization of the periodogram of the observed signal, whereas fine estimation can be accomplished in an open loop fashion or through an iterative procedure. On the one hand, all the open loop estimators exploit spectral interpolation to infer the
frequency residual from the analysis of the fast Fourier transform (FFT) coefficients at the maxima of the associated periodogram and at frequencies adjacent to it as in "B. Quinn, "Estimation of frequency, amplitude, and phase from the DFT of a time series," IEEE Trans. Signal Process., vol. 45, no. 3, pp. 814-817, Mar. 1997." and "Q. Candan, "Analysis and Further Improvement of Fine Resolution Frequency Estimation Method From Three DFT Samples," IEEE Signal Process. Lett., vol. 20, no. 9, pp. 913-916, Sep. 2013.". Unlike iterative estimators, the accuracy they achieve is frequency dependent and gets smaller when the signal frequency approaches the center of one of the FFT bins. On the other hand, the iterative estimation techniques available in the technical literature are based on various methods, namely on: a) numerical methods for locating the global maximum of a function (e.g., the secan method "D. Rife and R. Boorstyn, "Single tone parameter estimation from discrete-time observations," IEEE Trans. Inf. Theory, vol. 20, no. 5, pp. 591- 598, Sep. 1974." or the Newton's method "J. Selva, "ML Estimation and Detection of Multiple Frequencies Through Periodogram Estimate Refinement," IEEE Signal Process. Lett., vol. 24, no. 3, pp. 249-253, Mar. 2017."); b) an iterative method for binary search, known as the dichotomous search of the periodogram peak as in "E. Aboutanios, "A Modified Dichotomous Search Frequency Estimator," IEEE Signal Process. Lett., vol. 11, no. 2, pp. 186- 188, Feb. 2004."; c) interpolation methods amenable to iterative implementation like "E. Aboutanios and B. Mulgrew, "Iterative Frequency Estimation by Interpolation on Fourier Coefficients," IEEE Trans. Signal Process., vol. 53, no. 4, pp. 1237-1242, Apr. 2005. " ; d) the combination of the above mentioned dichotomous search with various interpolation techniques as "Y. V. Zakharov, V. M. Baronkin, and T. C. Tozer, "DFT-based frequency estimators
with narrow acquisition range," IEE Proc. - Commun., vol.148 , no. 1, pp. 1-7, Feb. 2001. e) the computation of the first derivative of the spectrum as in "C.-F. Huang, H.-P. Lu, and W.- H. Chieng, "Estimation of Single-Tone Signal Frequency with Special Reference to a Frequency-Modulated Continuous Wave System," Meas. Sci. Technol., vol. 23, no. 3, Mar. 2012.". Signal and System Models
In this session, we focus on the problem of estimating all the parameters of the complex sequence
and its real counterpart
with n = 0, 1, ..., N — 1 ; here, Ai and F( 6 0,1) denote the complex amplitude and the normalised frequency, respectively, of the l- th complex tone appearing in the right-hand side (RHS) of eq. (1), a( = |i4(| , xpi = A{Ai) , wc n is the n-th sample of an additive white Gaussian noise (AWGN) sequence (whose elements have zero mean and variance 2s2 ) , wr n = ¾{wcn), N is the overall number of samples, and ¾{x) and arg(x) denote the real part and the phase, respectively, of the complex quantity x. It is useful to point out that the signal model (2) can be rewritten as xr,n = åi=o [Qexp(/27mF;) + Q*exp(-y'27mF()] + wr n, ( 3 ) where
represents the complex amplitude of the real tone appearing in the RHS of eq. (2).
The signal models (1) and (2) appear in a number of problems concerning radar systems, wireless communications and biomedical applications. In the remaining part of this paragraph, we analyse
their meaning in the context of colocated and bistatic multiple- input multiple output (MIMO) radar systems operating in time division multiplexing (TDM) (see "S. M. Patole, M. Torlak, D. Wang, and M. Ali, "Automotive radars: A review of signal processing techniques," IEEE Signal Process. Mag., vol. 34, no. 2, pp. 22-35, Mar. 2017.") and at millimetre waves. These systems have the following important characteristics: a) They are equipped with a transmit (TX) antenna array and a receive (RX) antenna array, that consist of NT elements and NR elements, respectively; these allow to radiate orthogonal waveforms from different antennas and to receive distinct replicas of the electromagnetic echoes generated by multiple targets. b) Their TX and RX arrays are made of distinct antenna elements, placed at different positions. Moreover, TX antennas are close to the RX ones and, in particular, are usually placed on the same shield. c) Orthogonality of the transmitted waveforms is achieved by radiating them through distinct TX antennas over disjoint time intervals.
The architecture of the radar system considered in this work is illustrated in Fig. 1. In the remaining part of this section, two different models are described for the RF signal generated by the voltage controlled oscillator (VCO) of its transmitter and radiated by its TX array after power amplification. In the first case, corresponding to a frequency modulated continuous wave (FMCW) radar system, the VCO is fed by a periodic ramp generator; this produces a chirp FM signal, whose instantaneous frequency evolves periodically, as illustrated in Fig. 2. In this figure, the parameters T, TR and T0 represent the chirp interval, the reset time and the pulse period (or pulse
repetition interval), respectively, whereas the parameters f0 and B are the start frequency and the bandwidth, respectively, of the transmitted signal. If we assume that the radar system exploits all the available TX diversity (i.e., all its TX antennas) and that a time slot of T0 s is assigned to each TX antenna, transmission from the whole TX array is accomplished over an interval lasting TF = NTT0 s; this interval represents the duration of a single transmission frame.
Let us focus now on a single chirp interval and, in particular, on the time interval (0,G), and assume that, in that interval, the p-th TX antenna is employed by the considered radar system (with p 6 {0 , 1, ..., NT — 1}); the signal radiated by that antenna can be expressed as
S RF(0 = ARF ¾{s(t)}, (5) where ARF is its amplitude, s(t)= exp(/0(t)), (6)
and m = B/T is the chirp rate (i.e., the steepness of the generated frequency chirp). Let r RP(.t') denote the signal available at the output of the q-th RX antenna in the same time interval, with q = 0, 1, NR — 1 (see Fig. 1). This signal feeds a low noise amplifier (LNA), whose output undergoes downconversion, filtering and analog-to-digital conversion at a frequency fs = 1/Ts; here, Ts denotes the sampling period of the employed analog-to-digital converter (ADC). If the radiated signal sRF(t) (5) is reflected by L static point targets, the useful component of (t) consists of the superposition of L echoes, each originating from a distinct target. In this case,
if the propagation environment is static or undergoes slow variations, a simple mathematical model can be developed to represent the sequence of samples generated by the ADC in a single chirp interval. In deriving that model, the couple of the involved physical TX and RX antennas (namely, the p-th TX antenna and the q-th RX antenna) of the considered bistatic radar is usually replaced by a single virtual antenna of an equivalent monostatic radar. In particular, if it is assumed that all the TX and RX antennas belong to the same planar shield and a reference system lying on the physical antenna array is defined, the abscissa xv and the ordinate yv of the 17-th virtual antenna element associated with the p-th TX antenna and the q-th RX antenna can be computed as (e.g., see "J. Li and P. Stoica, Eds., MIMO Radar Signal Processing. Hoboken, NJ: J. Wiley & Sons, 2009. ")
and
_ yp+yg yp 2 (9) respectively, with v = 0, 1, NVR — 1; here,
are the coordinates of the TX (RX antenna) and NVR=NT-NR represents the overall size of the resulting virtual array. Based on these assumptions, the n-th received signal sample acquired through the 17-th virtual antenna element (with v = 0, 1, ..., NVR — 1) can be expressed as (e.g., see "J. Gamba, Radar Signal Processing for Autonomous Driving. Springer Singapore, 2020.")
with n = 0, 1, ..., N 1; here, N is the overall number of samples acquired over a chirp period, a( is the amplitude of the Z-th component of the useful signal,
is the normalized version of the frequency
characterizing the Z-th target detected on the v-th virtual RX antenna,
is the delay of the echo generated by the Z-th target and observed on the 17-th virtual channel, R\, 0; and fi denote the range of the Z-th target, its azimuth and its elevation, respectively, f\n) = 2nfQr\v) , ( 14 )
(nL and wrn, the n-th sample of the AWGN sequence affecting the received signal, is a Gaussian random variable having zero mean and variance s2 (assumed to be independent of v). It is important to point out that: a) the real signal model (10) is identical to that expressed by eq. (2) and can be adopted in all the FMCW radar systems that extract only the in-phase component of the signal captured by each RX antenna; b) some commercial MIMO radar devices provide both the in-phase and quadrature components of the received RF signals (e.g., see "M. A. Richards, Fundamentals of Radar Signal Processing. New York: McGraw-Hill, 2005."). In the last case, the complex model
equivalent to that expressed by eq. (1), must be adopted in place of its real counterpart (10) for any n; here,
for any v and Z.
The second case we consider in the generation of the radiated waveform corresponds to a stepped frequency continuous wave (SFCW) radar system represented in Fig. 3. Its name is
motivated by the fact that the VCO of its transmitter is fed by a staircase generator. For this reason, the instantaneous frequency of the resulting RF signal takes on N distinct and uniformly spaced values in an interval lasting T s for each TX antenna (see Fig. 4); the n-th value of the instantaneous frequency is fn=fo+nAf, (17) with n = 0, 1, , N — 1; here, f0 is the minimum radiated frequency, D/ is the frequency step size and N is the overall number of transmitted frequencies. It can be shown that, under the same assumptions made in the derivation of eq. (15), the measurement acquired through the v-th virtual element at the n- th frequency can be expressed as with v = 0, 1, ..., NV
R — 1; here, the amplitude ipt and the phase
(V)
At are still expressed by eqs. (14) and (16), respectively,
is the normalised delay characterizing the Z-th target and observed on the v-th virtual antenna, and the parameters (a;,
have exactly the same meaning as the one illustrated for the received signal model (15).
Then, in both the considered FMCW and SFCW radar systems, the signal observed on each virtual channel can be represented a superposition of L real or complex oscillations; moreover, the value of the parameter L has to be considered unknown. In the following derivations, the real samples
n = 0, 1, ..., N — 1} or their complex counterpart
n = 0, 1, ..., N — 1} acquired on the l?-th virtual channel are collected in the iV-dimensional vector
with z = r or c. This vector is processed by the next stages of the radar receiver for target detection and estimation. As it can be easily inferred from eqs. (10)—(12) ((15)— (16) and (18)— (19)), in the considered radar system, the problem of target detection and range estimation on the v-th virtual channel is equivalent to the classic problem of detecting multiple overlapped sinusoids (multiple overlapped complex exponentials) in the presence of AWGN and estimating their frequencies. In fact, if, in a FMCW radar system, the Z-th tone is found at the frequency ft , the presence of a target at the range (see eqs. (12) and (13)) p — lc f(v)
Kv ( ·1 ~ 2 21) ~m h is detected. Similarly, in a SFCW radar system, the normalised
—ft?) delay F; estimated on the the v -th virtual channel is associated with a target whose range is (see eqs. (13) and (19))
Kp _ hd p( v,l v) 1 (22)
2 —Af G, .
Information about the angular coordinates (namely, the azimuth and the elevation) of this target, instead, can be acquired
(V\ through the estimation of the set of NVR phases {ipt ; v = 0, 1, ..., NVR — 1} observed over the available virtual antennas. In fact, since (see eqs. (13) and (14)) where
is the wavelength associated with the frequency /0, the sequence
exhibits a periodic behavior characterized by the normalised horizontal spatial frequency
if the considered virtual elements form an horizontal uniform linear array (ULA), whose adjacent elements are spaced dH m apart. Dually, if a virtual vertical ULA is assumed, the periodic variations observed in the same sequence of phases are characterized by the normalised vertical spatial frequency
Fv,i A2^sin(0;), (26) where dv denotes the distance between adjacent elements of the virtual array itself. Consequently, angle finding can be easily accomplished by digital beamforming, i.e. by performing FFT processing on the estimated phases taken across multiple elements of the virtual array in a single frame interval.
From the considerations illustrated above, the following conclusions can be easily inferred: a) on the one hand, achieving precise estimation of the range of a given target requires the availability of an accurate estimate of the normalised frequency (or delay) of the real or complex tone associated with the target itself over at least one RX antenna; b) on the other hand, the quality of the estimate of the direction of arrival (DOA) characterizing the radar echo generated by a given target depends on the accuracy of the phase estimated over multiple virtual channels. The last consideration motivates the importance of accurately estimating this parameter over multiple antennas. Finally, in analysing the suitability of multiple tone estimators to colocated MIMO radar systems operating at millimetre waves, the following additional technical issues need to be taken carefully into account:
1) These radar systems often operate at short ranges and in the presence of extended targets. Each of resulting radar images is a cloud of point targets whose mutual spacing can be very small. For this reason, the accuracy of these images depends on
the frequency resolution (delay resolution) achieved on each virtual antenna in a FMCW (SFCW) radar system. In fact, this makes the receiver able to separate point targets characterized by similar ranges.
2) Information about the RCS of each point targets can be exploited in the classification of extended targets; for this reason, the availability of an accurate estimate of the amplitude of each radar echo can be very useful in a number of applications (e.g., in SAR imaging).
3) Distinct radar echoes can be characterized by substantially different signal-to-noise ratios (SNRs).
4) The number N of samples acquired over each virtual channel usually ranges from few hundreds to few thousands.
The last two issues explain why significant attention must be paid to the accuracy achieved by the adopted estimation algorithms at low SNRs and/or for relatively small values of N, since this can appreciably influence the quality of the generated radar image.
Summary of the invention
Scope of the present invention is to provide for methods for estimation of a single tone, which has a low computation complexity and that can be easily implemented in low computation power such has radars and in the base stations of cellular communication networks.
The basic idea of the present invention is based on the approximate maximum likelihood estimation of the single tone. The present invention proposes a new method that can be declined into different ways according to the reference domain, such as real or complex frequency or time domain. This is the reason for which the same inventive concept is described with the aid of the claims 1 and 10.
More in detail, when the single tone represents a frequency, the present invention proposes a method for computer estimation of the frequency, amplitude and phase of the single tone defining a signal acquired by means of a receiver, wherein the tone is real or complex; the method including
- estimation of a spectrum (X0) of said signal and its first NFFT — 1 derivatives ({Xfc;k = 1,...,NFFT — 1}) with NFFT = 3 when said signal is real and NFFT = 4 when said signal is complex, through FFT calculation with an oversampling factor (M ),
- iterative calculation of an estimate of the parameters of the most dominant tone, including phase (y), amplitude (a) and frequency (/) on the basis of said spectrum, and its derivatives, through NFFT — 1 FFTs.
When, instead, the single tone represents a delay, the present invention proposes a method for computer estimation of the delay, amplitude and phase of a single complex tone defining a signal acquired by means of a receiver; the method including
- estimation of the impulse response (X0) characterizing the
— (v) communication channel and its first NFFT — 1 derivatives ({Xfc ;k = 1,...,NFFT — 1}) with NFFT = 4, through IFFT calculation, with an oversampling factor (M),
- iterative calculation of an estimate of the parameters of the most dominant tone, including phase (ip), amplitude (a) and delay (t) on the basis of said impulse response, and its derivatives, through NFFT — 1 IFFTs.
It is immediately apparent that the two implementations depicted above are dual concepts.
It is also clear that the receiver can be part of a radar or of any device capable to acquire a single tone signal. Thus, the
present invention can be implemented also in 5G base stations. The present invetion aims at providing a total of three approaches about the estimation of a tone and their exploitation in colocated frequency modulated continuous wave (FMCW) and stepped frequency continuous wave (SFCW) multiple-input multiple-output (MIMO) radar systems operating at millimeter waves. In the present invention, a novel ML-based iterative estimator of a single real and complex tone is presented. The derivation of this estimator is based on: a) expressing the dependence of the ML metric on the frequency residual in an approximate polynomial form through standard approximations of trigonometric functions; b) exploiting the alternating minimization technique for the maximization of this metric. Moreover, its most relevant feature is represented by the fact that it requires the evaluation of spectral coefficients that are not exploited by all the other related estimation methods available in the technical literature.
It should be clear that the FFT represents the practical implementation of the DFT through computer means. Therefore, DFT and FFT can be used and referenced in an interchangeable way. The same concept applies also to IFFT and IDFT.
The present invention will be made fully clear with the aid of the attached dependent claim describing preferred embodiment of the invention.
Brief description of the figures
Figure 1 shows an example of architecture of a MIMO FMCW radar system;
Figure 2 discloses the time diagram of the instantaneous frequency of the RF signal radiated in a FMCW radar system;
Figure 3 shows an example of architecture of a MIMO SFCW radar system; Figure 4 discloses the time evolution of the instantaneous frequency of the RF signal radiated in a SFCW radar system;
Figure 5 shows a block diagram of the SFE#1 algorithm;
Figure 6 shows a block diagram of the SFE#2 algorithm;
Figure 7 shows a block diagram of the CSFE#1 algorithm; Figure 8 shows a block diagram of the CSFE#2 algorithm;
Figure 9 shows a block diagram of the CSDE#1 algorithm;
Figure 10 shows a block diagram of the CSDE#2 algorithm;
Detailed description of exemplary embodiments of the invention
Approximate Maximum Likelihood Estimation of a Single Tone In this section, we derive a new method for estimating the parameters of a single real or complex tone.
Estimation of a single frequency
Let us focus first on the problem of estimating the parameters (namely, the frequency and complex amplitude) of a single tone contained in the real sequence {xr,n; n = 0, 1, ..., N — 1}, whose n- th sample is expressed by eq. (10) with L = 1, i.e. as
or, equivalently, as (see eq. (3))
where (see eq. (4))
and F are the complex amplitude and the normalised frequency, respectively, of the tone itself. It is well known that the ML
estimates FML and CML of the parameters F and C, respectively, represent the solution of the least square problem
(FML< CML) = argmine(R,C), ( 30 )
F,C where F and C represent trial values of F and Cr respectively,
is the mean square error (MSE) evaluated over the whole observation interval,
is the square error between the noisy sample xrn (28) and its useful component sn(F,C)= Cexp (33)
evaluated under the assumption that F = F and C = C.
The cost function e(R, C) (31) can be easily reformulated in a different way as follows. To begin, substituting the RHS of eq. (33) in that of eq. (32) produces, after some manipulation,
fh = 2phR, (35) and ¾{x} (3{x}) denotes the real (imaginary) part of the complex variable x . Then, substituting the RHS of eq. (34) in that of eq. (31) yields
where
XR{P)± M{X (F)}, X,(P)A 3{X(P)},
is, up to the scale factor 1/iV, the Fourier transform of the sequence {xr,n},
Based on eq. (36), it is not difficult to show that the optimization problem expressed by eq. (30) does not admit a closed form solution because of the nonlinear dependence of the function e(R,C) on its variable P . However, an approximate solution to this problem can be derived by a) Exploiting an iterative method, known as alternating minimization (AM) (e.g., see "I. Ziskind and M. Wax, "Maximum likelihood localization of multiple sources by alternating projection," IEEE Trans. Acoust. Speech Signal Process., vol. 36, no. 10, pp. 1553-1560, Oct. 1988."). This allows us to transform the two-dimensional (2D) optimization problem expressed by eq. (30) into a couple of interconnected one¬ dimensional (ID) problems, one involving the variable P only (conditioned on the knowledge of C), the other one involving the variable C only (conditioned on the knowledge of P). b) Expressing the dependence of the function a(F,C) on the variable F through the couple (Fc,5) such that
P = Fc + dFdft, (40) where Fc is a given coarse estimate of F, d is a real variable called residual,
is the normalized fundamental frequency associated with the N0- th order discrete Fourier transform (DFT)
of the zero padded version
of the vector
collecting all the elements of the sequence {xr,n}, M is a positive integer (dubbed oversampling factor), 0D is a D —dimensional (column) null vector and iV0 = M ·N . c) Expressing the dependence of the function a(F,C) (36) on the variable d through its powers {di; 0 < l< 3}; this result is achieved by approximating various trigonometric functions appearing in the expression of s(F,C) with their Taylor expansions truncated to a proper order.
Let us show now how these principles can be put into practice. First of all, the exploitation of the above mentioned AM approach requires solving the following two sub-problems: PI) minimizing the cost function a(F,C) (36) with respect to C, given F = F; P2) minimizing the same function with respect to F, given C = C . Sub-problem PI can be easily solved thanks to the polynomial dependence of the cost function s(F,C) on the variables CR and Cj . In fact, the minimum of the function s(F,C) with respect to the last variables is the solution of the equations
and C,- X,(F)- C,gR(P)+ CR g,(F)= 0, (46) that result from computing the partial derivative of the RHS of eq. (36) with respect to CR and CIr respectively. Solving the linear system of equations (45)-(46) in the unknowns CR and Cj produces the optimal values
and f-- -XR (F)gi(F)+Xl(F) [1+9R(F)]
(48) i-ls(F)l2 of CR and Cj , respectively. Putting together the last two formulas yields
Therefore, given F = F, the optimal value C of the variable C can be computed exactly through the last equation; this requires the evaluation of X(F) and g(F) (see eqs. (38) and (39), respectively). Note that, on the one hand, g(F) can be easily evaluated through its exact expression
1exp(-;'4reNF)-l 5(0 7 (50)
N7 exp(-;'4reF)-l’ which is easily derived from its definition (39). On the other hand, 7(F) can be computed exactly through its expression (38) or, in an approximate fashion, through a computationally efficient procedure based on the fact that the vector
Xs = M X0 (51) collects iV0 uniformly spaced samples of the function X(F) (38). In fact, the k-th element
of the vector X0 (42) can be expressed as
where
Ffc= k FDFT (54) with k = 0 , 1 , ..., iV0 — 1 . For this reason, an approximate evaluation of the quantity X(F) at a normalised frequency F different from any multiple of FDFT (41) can be accomplished by
interpolating the elements of the vector Xs (51); the last vector, in turn, can be easily computed after evaluating the N0- th fast Fourier transform (FFT) of CO,ZR (43), i.e. the vector X0 (42).
Let us take into consideration now sub-problem P2 . Unluckily, this sub-problem, unlike the previous one, does not admit a closed form solution. For this reason, an approximate solution is developed below. Such a solution is based on representing the normalized frequency F in the same form as F (see eq. (40)), i.e. as
F = Fc + d FDFT (55) and on a novel method for estimating the real residual d, i.e. for accomplishing the fine estimation of F . This method is derived as follows. Representing the trial normalized frequency F according to eq. (40) allows us to express the variable fh (35) as
Fti @ n FtiA, (56) where
D = 2pd FQ PT ( 57 ) is a new variable and
§n = 2 ph Fc. (58)
Then, substituting the RHS of eq. (56) in that of eq. (34) (with C = C) yields
that can be rewritten as
£ n(F, C) = X2 in + 2 [c2 + Cf]
+2(C| — C2)[cos(20n)cos(2nA) — sin(20n)sin(2nA)]
— 4CRC/[sin(20n)cos(2n7\) + cos(20n)sin(2nA)]. ( 60)
If the normalized frequency FDFT (41) is small enough (i.e., if the FFT order N0 is large enough) , the trigonometric functions cos (knK) and sin (knK) appearing in the RHS of eq. (60) (with k = 1 and 2) can be approximated as sin(nA) — nA, (61)
and
Substituting the RHSs of eqs . (61) -(64) in that of eq. (60) and, then, substituting the resulting approximate expression in the RHS of eq. (31) yields
£SFE(A C) =
here, p = FC/Fdft, 66
gp\n\= np for p = 1, 2 and 3, g0\n\= 1,
for any p and k = 1,2, and
Xk,n — 71 ·Cg-gi (69) with n = 0, 1, ..., N — 1. It is important to point out that: a) if p is an integer, the quantity Xkp (68) (with k = 1 and 2 ) represents the p-th element of the vector
generated by the N0-th order DFT of the zero padded version
of the vector
b) if p is not an integer, the quantity Xkp can be evaluated exactly on the basis of eq. (68) or, in an approximate fashion, by interpolating a subset of adjacent elements of the vector Xfc (70); c) the evaluation of the vectors {Xfc; k = 1,2} requires two additional FFTs.
Since the function e5RE(D,C) (65) is a polynomial of degree three in the variable D, an estimate D of D and, consequently, an estimate (see eq. (57))
of d, can be obtained by computing the derivative of this function with respect to D, setting it to zero and solving the resulting quadratic equation
a(r)D2 + b(p)D+ c(p)= 0, (74) in the variable D; here,
and
Note that only one of the two solutions of eq. (74), namely b(p)+J(b(p))Z®a(p)c(p)
D (78) 2a(p) has to be employed. A simpler estimate of D is obtained neglecting the contribution of the first term in the left-hand side of eq. (74), i.e. setting a(p)= 0. This leads to a first- degree equation, whose solution is
D = —c{p)/b{p). (79)
Given an estimate D of D (and, consequently, and estimate d of d; see eq. (73)), the fine estimate
of F can be evaluated on the basis of eq. (55).
The mathematical results derived above allow us to derive novel estimation algorithms. Our first estimation algorithm, called single frequency estimator #1, (SFE #1), for iteratively estimating the normalized frequency F and the complex amplitude C is initialised by
1) Evaluating: a) the vector X0 (42); b) the initial coarse estimate F^ of F as
Fc (0)= aFDFT, (81) where the integer ά is computed by means of the well known
periodogram method (e.g., see "E. Aboutanios and B. Mulgrew, "Iterative Frequency Estimation by Interpolation on Fourier Coefficients," IEEE Trans. Signal Process., vol. 53, no. 4, pp. 1237-1242, Apr. 2005."), i.e. as ά = are
c) the initial estimate C® of C on the basis of eq. (49) with F =
; d) the spectral coefficients
and X2p on the basis of eq. (68); e) the coefficients {Kp(2a); p = 1, 2} ({Kp(2a ); p = 1, 2, 3}) according to eq. (67), the coefficients (a(a), ό(a), c(a)} ( {&(<$) , c(a)} according to eqs. (75)-(77) and the first estimate D^0-1 of D on the basis of eq. (78) (eq. (79)); f) the first fine estimate F® of F on the basis of eqs. (55) and (73), i.e. as
2) Setting its iteration index i to 1.
Then, an iterative procedure is started. The i-th iteration is fed by the estimates
of F and C, respectively, and produces the new estimates F® and
of the same quantities (with i = 1 , 2 , . .. , iVSFE , where iVSFE is the overall number of iterations); the procedure employed for the evaluation of F® and C® consists of the two steps described below (the p-th step is denoted SFE#1-Sp). SFE#1-Sl - The new estimate D® of D is computed through eq. (78) (eq. (79)); in the evaluation of the coefficients (a(p), b(p), c(p)} ({b(p) , c(p)}) appearing in the RHS of these equations, is assumed. Then,
is evaluated.
SFE#1-S2 - The new estimate C® of C is evaluated through eq. (49); F = F® is assumed in this case. Moreover, the index i is incremented by one before starting the next iteration.
At the end of the last (i.e., of the iVSFE-th) iteration, the fine estimates F = F(WSFE) and C = C^Nsfe^ of F and C, respectively, become available.
It is important to point out that: a) The estimate 5® of the residual d computed by the SFE#1 in its i-th iteration is expected to become smaller as i increases, since F® should progressively approach F if our algorithm converges. b) The coefficients {Kp(2a); p = 1, 2, 3} can be computed exactly on the basis of eq. (67). However, since the definition (67) can be put in the equivalent form
where
the identities
and
+(3N3 - 3 N2 -3 N - 1 )qN+1 (89)
+( —37V3 + 6 N2 - 4 )qN+2 + (N - 1 )3qN+3 (90) holding for any q 6 (C, can be exploited for the efficient
computation of Kv(2d) for any p and a. c) The quantities
k = 1,2} required in the first step of the i-th iteration can be computed exactly on the basis of eq. (68). However, they can be also evaluated, in an approximate fashion, by interpolating / adjacent elements of the iV0 - dimensional vectors Xfc (70), where / denotes the selected interpolation order. d) The estimate D® evaluated according to eq. (79) is expected to be less accurate than that computed on the basis of eq. (78). However, our numerical results have evidenced that both solutions achieve similar accuracy. e) The SFE#1 can be employed even if the single tone appearing in the RHS of eq. (27) is replaced by the superposition of L distinct tones (see eq. (3)). In this case, the strongest (i.e., the dominant) tone is detected through the periodogram method (see eq. (82)) and the parameters of this tone are estimated in the presence of both Gaussian noise and the interference due to the remaining tones. Therefore, the estimation accuracy of the SFE#1 is affected by both the amplitudes and the frequencies of the other (L— 1) tones. f) A stopping criterion, based on the trend of the sequence {D® ; /= 1, 2, ...}, can be easily formulated for the SFE#1. For instance, the execution of its two steps can be stopped if, at the end of the i-th iteration, the condition
|D®|<eD (91) is satisfied; here, eD represents a proper threshold.
Our second estimation algorithm (called single frequency estimator #2, SFE#2) relies on the same mathematical results exploited in the derivation of the SFE#1 and has the same structure. However, it is also based on the idea of taking the
estimate generated in its (i— l)-th iteration as a coarse estimate of F in the following (i.e., in the i-th) iteration. This leads to representing the (unknown) normalised frequency F as
F = F(i_1) + 5Fdft, (92) where d is a residual whose estimate 5®, generated in the i-th iteration of SFE#2, is expected to tend to zero as / increases. The initialization phase of SFE#2 is very similar to that of SFE#1, the only difference being represented by the fact that the quantity
is also defined. The first step of SFE#2 is different from the corresponding step of SFE#1 and can be summarised as follows.
SFE#2-Sl) The new estimate D® of D is computed through eq. (78) (eq. (79)); in the evaluation of the coefficients (a(p), b(p), c(p)} ({b(p) , c(p)}) appearing in the RHS of these equations, C = CV-V and
are assumed. Then,
F(i) = F(i-i) +— (95)
2TC is evaluated.
The second step of SFE#2, instead, is the same as that of SFE#1.
Preferably, the method for the estimation of the frequency, the amplitude and the phase of a real sequence (see Steps 2.1.0 - 2.1.2 of Fig. 5) can be described as follows:
- (Step 2.1.0) representation of a tone normalized frequency F = f/fs r where fs is the sampling frequency, as sum of a main portion (QCFdft) and a residual portion (5Fdft), where FDFT = l/iV0, where
N0 is the FFT order, wherein said main portion is determined by searching the index (a) of the vector element Xoa corresponding to the maximum absolute value in a first half plus one elements of the same vector X0 and said residual frequency is found through an iterative procedure together with the corresponding tone complex amplitude C= aexp(jxp)/2, where a is said tone amplitude and y is said tone phase.
- (Step 2.1.1) said iterative procedure includes a preliminary initialization comprising: a) setting the iteration index i to 1, the initial estimate of tone complex amplitude C equal to said vector element (Xoa)/ the initial estimate of A = 2pdFDFT to zero and the initial estimate of F to said main portion (CCFdft); b) computing the complex coefficients (Kp(2a);p = 1,2,3} as the 2a- th element of the iV0-order FFT of the sequence np, with p = 1,2,3; then iteratively
- (Step 2.1.2) the i-th iteration is fed by the estimates of D, F and C computed in the previous iteration, to generate new estimates of D , F and C , for a predetermined number of iterations, wherein each iteration includes a) a first step where the residual D is calculated solving the first order equation Z¾D + c¾ = 0 OR the second order equation adD2 + ί¾D + CQ = 0 in the variable D with:
where ¾{x} denotes the real part and 3{x} denotes the imaginary part of the complex quantity x; discarding one of the two solutions, which does not satisfy the inequalities: if a = 0
otherwise if said second order equation is solved; then the new estimate of F is obtained by adding said residual (A Fdft) to said main portion (afDFT); b) a second step wherein said new estimate of F is employed to compute the new estimate of C as
where
is generated by interpolating the elements of said vector MX0 where M is said oversampling factor employed in said FFT calculations.
The above mentioned Steps 2.1.1 - 2.1.2 of Fig. 5 can be replaced by an alternative method including (see Steps 2.2.1 - 2.2.2 of Fig. 6):
- (Step 2.2.1) said iterative procedure includes a preliminary initialization comprising: a) setting the iteration index i to 1, the initial estimate of C equal to said vector element (^o,a)/ the initial estimate of A = 27T5 FDFT to zero and the initial estimate of F to said main portion (aFDFT); b) computing the complex coefficients (Kp(2a); p = 1,2,3} as the 2a- th element of the iV0-order FFT of the sequence np, with p = 1, 2, 3 ; c) computing the real quantity p = F/FDFT; then iteratively
- (Step 2.2.2) the i-th iteration is fed by the estimates of D, F and C and the real quantity p computed in the previous iteration, to generate new estimates of D, F, C and p , for a predetermined number of iterations, wherein each iteration includes a) a first step where the residual D is calculated solving the first order equation Z¾D + c¾ = 0 OR the second order equation adD2 + Z¾D + CQ = 0 in the variable D with:
where ¾{x} denotes the real part, 3{L:} denotes the imaginary part of the complex quan e computed as:
with p = 1, 2, 3 and ix k,p> k = 1-2} as:
where xrn is the n-th element of the said real signal; said coefficients ({Xkp;k = 1,2}) can be also computed by interpolating the elements of said spectrum (Xfc); discarding one of the two solutions, which does not satisfy the inequalities:
(0 < D < 7 TFDFT if a = 0
—TTFQPT < D < TTFQPT otherwise
if said second order equation is solved; then the new estimate of F is obtained by adding said residual (A Fdft ) to said estimate of F; b) a second step wherein the new estimate of p is computed by dividing said new estimate of F by Fdft ; c) a third step wherein said new estimate of F is employed to compute the new estimate of C as:
where is generated by
interpolating the elements of said vector MX0 where M is said oversampling factor employed in said FFT calculations.
All the results developed above refer to the real sequence {xr,n}, whose n-th element is expressed by eq. (27). However, similar estimation method - dubbed complex SFEs, CSFEs - can be developed for the complex counterpart, i.e. for a complex sequence {xc,n} such that (see eq. (15) with L = 1 ) xc,n = Aexp(j2nnF) + wc n, (96) with n = 0, 1 , N — l . In this case, the ML estimation of the parameters F and A can be formulated in a similar way as eq. (30), the only differences being represented by the fact that: a) the parameter C is replaced by A ; b) the term sn(F, C) appearing in the RHS of eq. (31) is replaced by
where
represents the useful component of xc n (96) under the assumption that F = F and A = A . Substituting the RHS of the last equation
in that of eq. (97) and, then, the resulting expression in the RHS of the MSE (see eq. (1))
yields, after some manipulation,
where
The procedure adopted to minimise the last function with respect to F and A is similar to that illustrated above for the SFE. For this reason, in the following we limit to show essential mathematical results only. First all, it is easy to show that, given F = F, the function s{F,A) (100) is minimised with respect to A selecting
A=A = X(F), (102) where X(F) can be computed exactly through its expression (see eq. (38)) replacing the real sequence {xr>n} with the complex sequence {xc,n} ) or, in an approximate fashion, through an interpolation technique following the same mathematical approach as that illustrated for the SFE. Note also that the mathematical result expressed by eq. (102) is exact. On the contrary, given A = A, a closed form expression for the value of F minimising the function s(F,i4) cannot be derived because of the nonlinear dependence of this function on F. However, following the same mathematical approach as that illustrated for the SFE, the approximate expression
can be obtained for the cost function s{F,A) (100); the quantities appearing in the last equation have the same meaning as the ones defined for eq. (65). It worth stressing that: a) The considerations expressed about the evaluation of the quantities {Xk,p} in our derivation of the SFEs apply to the CFSEs too. However, in this case, three additional spectral coefficients, namely {Xk p; k = 1,2,3}, need to be computed. b) The initial estimate p^ = a (see eq. (93)) of p is evaluated in a similar way as the SFEs (see eq. (82)), i.e., as a = arg„ <ZE{0m,1,.a..,xNQ —1J}*o,a|· (104)
The minimization of the function ^CSFEC^W4) (103) with respect to D is achieved by taking its partial derivative with respect to this variable and setting it to zero; this results again in the quadratic equation (74), whose coefficients are
b{p)= ¾{i4*X2,p} (106) and c(p)A -3{i¾ ,p}. (107)
Following the approach adopted in the development of the SFEs and exploiting the mathematical results expressed by eq. (102) and by eqs. (105)-(107) allow us to easily derive the CSFEs. Preferably, the method for the estimation of the frequency, the amplitude and the phase of a complex sequence (see Steps 2.3.0 - 2.3.2 of Fig. 7) can be described as follows:
- (Step 2.3.0) representation of a tone normalized frequency F = f/fs r where fs is the sampling frequency, as sum of a main portion
(OCFdft) and a residual portion (5Fdft), where FDFT = l/iV0, where iV0 is the FFT order, wherein said main portion is determined by searching the index (a) of the vector element X0Q corresponding to the maximum absolute value of the same vector X0 and said residual frequency is found through an iterative procedure together with the corresponding tone complex amplitude A = aexp(jxp), where a is said target amplitude and y is said target phase.
- (Step 2.3.1) said iterative procedure includes a preliminary initialization comprising: setting the iteration index i to 1, the initial estimate of A equal to said vector element (Xo,a)t the initial estimate of A = 2i FDFT to zero and the initial estimate of F to said main portion (CCFdft); then iteratively
- (Step 2.3.2) the i-th iteration is fed by the estimates of D, F and A computed in the previous iteration, to generate new estimates of D , F and A , for a predetermined number of iterations, wherein each iteration includes a) a first step where the residual D is calculated solving the first order equation Z¾D + % = 0 OR the second order equation
in the variable D with: aa = bji4 X3Qj,
where ¾{x} denotes the real part and 3{x} denotes the imaginary part of the complex quantity x; discarding one of the two solutions, which does not satisfy the inequalities: if a = 0
otherwise if said second order equation is solved; then the new estimate of F is obtained by adding said residual (A Fdft) to said main portion (afDFT); b) a second step wherein said new estimate of F is employed to compute the new estimate of A as A = Xint{P) and Xintif^) is generated by interpolating the elements of said vector MX0 where M is said oversampling factor employed in said FFT calculation. The above mentioned Steps 2.3.1 - 2.3.2 of Fig. 7 can be replaced by an alternative method including (see Steps 2.4.1 - 2.4.2 of Fig. 8):
- (Step 2.4.1) said iterative procedure includes a preliminary initialization comprising: a) setting the iteration index i to 1, the initial estimate of A equal to said vector element (X0jg)/ the initial estimate of A = 27T5 FDFT to zero and the initial estimate of F to said main portion (a FDFT ) ; b) computing the real quantity p = F/FDFT; then iteratively
- (Step 2.4.2) the i-th iteration is fed by the estimates of D, F and A and the real quantity p computed in the previous iteration, to generate new estimates of D, F, A and p, for a predetermined number of iterations, wherein each iteration includes:
a) a first step where the residual D is calculated solving the first order equation b^A + % = 0 OR the second order equation adD2 F ί¾D F CQ = 0 in the variable D with:
where ¾{x} denotes the real part, 3{x} denotes the imaginary part of the complex quantity x and {Xkp; k = 1,2,3} is computed as:
where xcn is the n -th element of the said complex; said coefficients ({Xkp; k = 1,2,3}) can be also computed by interpolating the elements of said spectrum (Xfc); discarding one of the two solutions, which does not satisfy the inequalities: ifa = 0
otherwise if said second order equation is solved; then the new estimate of F is obtained by adding said residual (A Fdft) to said estimate of F; b) a second step wherein the new estimate of p is computed by dividing said new estimate of F by Fdft; c) a third step wherein said new estimate of F is employed to compute the new estimate of A as A = Xint
generated by interpolating the elements of said vector MX0 where M is said oversampling factor employed in said FFT calculations.
The above-described methods for the estimation of the frequency, amplitude and phase of a complex sequence can be employed for range (R) calculation of a point target detected through a single-input single-output frequency modulated continuous wave radar system as
where m is the radar chirp rate, c is the speed of light, fs is the sampling frequency and F is said estimated normalized frequency calculated according to said methods.
Additionally, said methods can be employed for angle of arrival (Q) calculation of a point target having range (R) detected through a multiple-input multiple-output frequency modulated continuous wave radar system equipped with a plurality of transmitting (TX) and receiving (RX) antennas, wherein each couple of said TX and RX antennas is replaced with the equivalent virtual antenna; said virtual antennas forming a virtual array composed by one uniform linear array; said angle of arrival being computed as: q= arcsin
here, l is the radar wavelength, d is the distance between two adjacent virtual antennas of said uniform linear array, and F is said estimated normalized frequency according to said methods. Estimation of single delays
All the estimation techniques developed in the previous two paragraphs refer to the real (complex) sequence {xr,n} ({xc,n})r whose n-th element is expressed by eq. (10) (eq. (15)). As shown in the previous section, these models are suitable to represent
the sampled downconverted signal made available by the receiver of a colocated FMCW MIMO radar. However, similar estimation methods can be also developed for a colocated SFCW MIMO radar system, since the complex signal model (18) illustrated for this system is very similar to the model (15) referring to its complex FMCW counterpart, the only differences being represented by the fact that: a) the physical meaning of the parameter F; , since this represents a normalized delay in eq. (18) and a normalized frequency in eq. (15); b) the sign of the argument of all the complex exponentials appearing in the RHS of eq. (18), since this is the opposite of that of the corresponding functions contained in the RHS of eq. (15). The implications of this sign reversal can be summarised as follows. In the derivation of the single tone estimator, the function X(F) (38) is replaced by
where xcn is the n-th sample of the received signal sequence (see eq. (18)). Consequently, the vectors X0 (42) and X¾ (70) are generated by taking the iV0 order inverse DFT (IDFT) of the zero- padded vectors x0 ZP (43) and xfc ZP (71), respectively; note that the n-th element of the N0-th dimensional vector x0,zp (x /c,zp ) is equal to xcn (nk ■ xcn; see eq. (69)) for n = 0 , 1, ..., N — 1 and is equal to zero for n > N . Two algorithms for estimating a single delay can be easily developed by following the same line of reasoning as the CSFEs; these are called complex single delay estimators (CSDEs). More specifically, the CSDEs compute an estimate d of the delay residual d by solving eq. (74) in the variable D (57); in this case, the coefficients a(p) and c(p) are still expressed by eqs. (105) and (107), respectively, whereas b(p)A -¾{i*X2,p}, (109) where A is the estimate of the complex gain A computed on the
basis of eq. (102). Moreover, in the evaluation of the coefficients {a(p), b(p), c(p)} of eq. (74), the quantity Xkp(i-q (see eq. (68)) is replaced by
where pO-1-1 is still defined by eq. (94).
Preferably, the method for the estimation of the delay, the amplitude and the phase of a complex sequence (see Steps 2.5.0
- 2.5.2 of Fig. 9) can be described as follows:
- (Step 2.5.0) representation of a tone normalized delay T = t Af , where Af is the frequency step size of the employed radar, as sum of a main portion (QCFidft) and a residual portion (5Fidft), where FIDFT = l/iV0, where iV0 is the IFFT order, wherein said main portion is determined by searching the index (a) of the vector element Xoa corresponding to the maximum absolute value of the same vector X0 and said residual portion is found through an iterative procedure together with the corresponding target complex amplitude A = a exp(-jxp), where a is said target amplitude and y is said target phase.
- (Step 2.5.1) said iterative procedure includes a preliminary initialization comprising: setting the iteration index i to 1, the initial estimate of A equal to said vector element (Xo,a) r the initial estimate of D = 2i FIDFT to zero and the initial estimate of T to said main portion (aFIDFT); then iteratively
- (Step 2.5.2) the i-th iteration is fed by the estimates of D, T and A computed in the previous iteration, to generate new
estimates of D , T and A , for a predetermined number of iterations, wherein each iteration includes a) a first step where the residual D is calculated solving the first order equation Z¾D + c¾ = 0 OR the second order equation adD2 + Z¾D + CQ = 0 in the variable D with:
where ¾{x} denotes the real part and 3{x} denotes the imaginary part of the complex quantity x ; discarding one of the two solutions, which does not satisfy the inequalities: ifa = 0
otherwise if said second order equation is solved; then the new estimate of T is obtained by adding said residual (A Fidft) to said main portion (afIDFT); b) a second step wherein said new estimate of T is employed to compute the new estimate of A as A = Xint
generated by interpolating the elements of said vector MX0 where M is said oversampling factor employed in said IFFT calculation. The above mentioned Steps 2.5.1 - 2.5.2 of Fig. 9 can be replaced by an alternative method including (see Steps 2.6.1 - 2.6.2 of Fig. 10):
- (Step 2.6.1) said iterative procedure includes a preliminary initialization comprising: a) setting the iteration index i to 1, the initial estimate of A equal to said vector element X0¾, the initial estimate of D =
2ir<5FIDFT to zero and the initial estimate of T to said main portion (afIDFT); b) computing the real quantity p = T/FIDFT; then iteratively - (Step 2.6.2) the i-th iteration is fed by the estimates of D,
T, A and the real quantity p computed in the previous iteration, to generate new estimates of D, T, A and p, for a predetermined number of iterations, wherein each iteration includes: a) a first step where the residual D is calculated solving the first order equation b%A + c¾ = 0 OR the second order equation ¾D2 + ¾D + ¾ = 0 in the variable D with:
where ¾{x} denotes the real part, 3{x} denotes the imaginary part of the complex quant computed as:
where xc n is the n-th element of the said complex signal; said coefficients ({Xkp; k = 1,2,3}) can be also computed by interpolating the samples of said impulse response (X¾); discarding one of the two solutions, which does not satisfy the inequalities: if a = 0
otherwise if said second order equation is solved; then the new estimate of T is obtained by adding said residual (AFidft) to said estimate of T ;
b) a second step wherein the new estimate of p is computed by dividing said new estimate of T by Fidft; c) a third step wherein said new estimate of T is employed to compute the new estimate of is
generated by interpolating the elements of said vector MX0 where M is said oversampling factor employed in said IFFT calculations.
The above-described methods for the estimation of the delay, amplitude and phase of a complex sequence can be employed for range (R) calculation of a point target detected through a single-input single-output stepped frequency continuous wave radar system as
where D/ is the frequency step size of the employed radar system, c is the speed of light and T is said estimated normalized delay, calculated according to said methods.
Additionally, said methods can be employed for angle of arrival (Q) calculation of a point target having range (R) detected through a multiple-input multiple-output stepped frequency continuous wave radar system equipped with a plurality of transmitting (TX) and receiving (RX) antennas, wherein each couple of said TX and RX antennas is replaced with the equivalent virtual antenna; said virtual antennas forming a virtual array composed by one uniform linear array; said angle of arrival being computed as:
here, l is the radar wavelength, d is the distance between two
adjacent virtual antennas of said uniform linear array, and T is said estimated normalized delay according to said methods.
This invention can be implemented advantageously in a computer program comprising program code means for performing one or more steps of such method, when such program is run on a computer. For this reason, the patent shall also cover such computer program and the computer-readable medium that comprises a recorded message, such computer-readable medium comprising the program code means for performing one or more steps of such method, when such program is run on a computer, such as a computing device associated to a radar system or to a base station of cellular telecommunication, such as 5G and similar ones.
Many changes, modifications, variations and other uses and applications of the subject invention will become apparent to those skilled in the art after considering the specification and the accompanying drawings which disclose preferred embodiments thereof as described in the appended claims.
The features disclosed in the prior art background are introduced only in order to better understand the invention and not as a declaration about the existence of known prior art. In addition, said features define the context of the present invention, thus such features shall be considered in common with the detailed description.
Further implementation details will not be described, as the man skilled in the art is able to carry out the invention starting from the teaching of the above description.
Claims
1. Method for computer estimation of the frequency, amplitude and phase of a single tone defining a signal acquired by means of a receiver, wherein the tone is real or complex; said method including:
- (Step 1: FFT processing) estimation of a spectrum (X0) of said signal and its first NFFT — 1 derivatives ({Xfc;k= 1,.. NFFT — 1}) with NFFT = 3 when said signal is real and NFFT = 4 when said signal is complex, through FFT calculation with an oversampling factor (M),
- (Step 2) iterative calculation of an estimate of the parameters of the most dominant tone, including phase (y), amplitude (a) and frequency (/) on the basis of said spectrum, and its derivatives, through NFFT — 1 FFTs.
2. [SFE] Method according to claim 1, wherein said signal is real, said (Step 2) calculation includes:
- (Step 2.1.0) representation of a tone normalized frequency F = f/fsr where fs is the sampling frequency, as sum of a main portion (aFdft) and a residual portion (5Fdft), where FDFT = l/iV0, where iV0 is the FFT order, wherein said main portion is determined by searching the index (a) of the vector element Xoa corresponding to the maximum absolute value in a first half plus one elements of the same vector X0 and said residual frequency is found through said iterative calculation together with the corresponding tone complex amplitude C = aexp(jxp)/2, where a is said tone amplitude and y is said tone phase.
3. [SFE#1] Method according to claim 2, wherein said iterative procedure includes a preliminary initialization (Step 2.1.1)
comprising: a) setting the iteration index i to 1, the initial estimate of tone complex amplitude C equal to said vector element (X0,a)t the initial estimate of A = 2pd FDFT to zero and the initial estimate of F to said main portion (CCFdft); b) computing the complex coefficients (Kp(2a); p = 1,2,3} as the 2a- th element of the iV0-order FFT of the sequence np, with p = 1, 2, 3 ; then iteratively
- (Step 2.1.2) the i-th iteration is fed by the estimates of D, F and C computed in the previous iteration, to generate new estimates of D , F and C , for a predetermined number of iterations, wherein each iteration includes a) a first step where the residual D is calculated solving the first order equation b^K + c¾ = 0 OR the second order equation adD2 + ί¾D + CQ = 0 in the variable D with:
where ¾{x} denotes the real part and 3{x} denotes the imaginary part of the complex quantity x; discarding one of the two solutions, which does not satisfy the inequalities: if a = 0
otherwise
if said second order equation is solved; then the new estimate of F is obtained by adding said residual (AFdft) to said main portion (aFDFT); b) a second step wherein said new estimate of F is employed to compute the new estimate of C as
where
is generated by interpolating the elements of said vector MX0 where M is said oversampling factor employed in said FFT calculations.
4. [SFE#2] Method according to claim 2, wherein said iterative procedure includes a preliminary initialization (Step 2.2.1) comprising: a) setting the iteration index i to 1, the initial estimate of C equal to said vector element (X0,a)t the initial estimate of A =
2i FDFT to zero and the initial estimate of F to said main portion (afDFT); b) computing the complex coefficients (Kp(2a);p = 1,2,3} as the 2a- th element of the iV0-order FFT of the sequence np , with p = 1,2,3; c) computing the real quantity p=F/FDFT; then iteratively
- (Step 2.2.2) the i-th iteration is fed by the estimates of D, F and C and the real quantity p computed in the previous iteration, to generate new estimates of D, F, C and p, for a predetermined number of iterations, wherein each iteration includes
a) a first step where the residual D is calculated solving the first order equation Z¾D + c¾ = 0 OR the second order equation adD2 + Z¾D + CQ = 0 in the variable D with: a? - 23 {(r)2 K3(2p)j,
where ¾{x} denotes the real part, 3{x} denotes the imaginary part of the complex quantity x, {Kr(2b); p = 1,2,3} are computed as:
with p = 1, 2, 3
as:
where cLh is the n-th element of the said real signal; said coefficients ({Xk,p> k = 1,2}) are also computed by interpolating the elements of said spectrum (Xfc ); discarding one of the two solutions, which does not satisfy the inequalities: if a = 0
otherwise if said second order equation is solved; then the new estimate of F is obtained by adding said residual (A Fdft) to said estimate of F;
b) a second step wherein the new estimate of p is computed by dividing said new estimate of F by Fdft; c) a third step wherein said new estimate of F is employed to compute the new estimate of C as:
where
is generated by interpolating the elements of said vector MX0 where M is said oversampling factor employed in said FFT calculations.
5. [CSFE] Method according to claim 1, wherein said signal is complex, said (Step 2) calculation includes:
- (Step 2.3.0) representation of a tone normalized frequency F = f/fsr where fs is the sampling frequency, as sum of a main portion (aFdft) and a residual portion (5Fdft), where FDFT = l/iV0, where N0 is the FFT order, wherein said main portion is determined by searching the index (a) of the vector element Xoa corresponding to the maximum absolute value of the same vector X0 and said residual frequency is found through said iterative calculation together with the corresponding tone complex amplitude A = aexp(jxp), where a is said target amplitude and y is said target phase.
6. [CSFE#1] Method according to claim 5, wherein said iterative procedure includes a preliminary initialization (Step 2.3.1) comprising: setting the iteration index i to 1, the initial estimate of A equal to said vector element (Xo,a)t the initial estimate of A =
2ir<5FDFT to zero and the initial estimate of F to said main portion (afDFT); then iteratively
- (Step 2.3.2) the i-th iteration is fed by the estimates of D, F and A computed in the previous iteration, to generate new estimates of D , F and A , for a predetermined number of iterations, wherein each iteration includes a) a first step where the residual D is calculated solving the first order equation b^L·+ c¾ = 0 OR the second order equation adD2 + Z¾D + CQ = 0 in the variable D with:
where ¾{x} denotes the real part and 3{x} denotes the imaginary part of the complex quantity x; discarding one of the two solutions, which does not satisfy the inequalities: ifa = 0
otherwise if said second order equation is solved; then the new estimate of F is obtained by adding said residual (AFdft) to said main portion (afDFT); b) a second step wherein said new estimate of F is employed to compute the new estimate of A as A =
generated by interpolating the elements of said vector MX0 where M is said oversampling factor employed in said FFT calculation.
7. [CSFE#2] Method according to claim 5, wherein said iterative procedure includes a preliminary initialization (Step 2.4.1) comprising: a) setting the iteration index i to 1, the initial estimate of A equal to said vector element (X0@) , the initial estimate of A = 2i FDFT to zero and the initial estimate of F to said main portion (aFDFT); b) computing the real quantity p=F/FDFT; then iteratively
- (Step 2.4.2) the i-th iteration is fed by the estimates of D, F and A and the real quantity p computed in the previous iteration, to generate new estimates of D, F, A and p, for a predetermined number of iterations, wherein each iteration includes: a) a first step where the residual D is calculated solving the first order equation b^A+ CQ = 0 OR the second order equation
where ¾{x} denotes the real part, denotes the imaginary part of the complex quantity x and {Xkp;k = 1,2,3} is computed as:
where xcn is the n -th element of the said complex; said coefficients ({Xkp; k = 1,2,3}) can be also computed by interpolating the elements of said spectrum (Xfc); discarding one of the two solutions, which does not satisfy the inequalities: ifa = 0
otherwise if said second order equation is solved; then the new estimate of F is obtained by adding said residual (A Fdft) to said estimate of F; b) a second step wherein the new estimate of p is computed by dividing said new estimate of F by Fdft; c) a third step wherein said new estimate of F is employed to compute the new estimate of A as A = Xint
is generated by interpolating the elements of said vector MX0 where M is said oversampling factor employed in said FFT calculations.
8. Method for range (R) calculation of a point target detected through a single-input single-output frequency modulated continuous wave radar system as
where m is the radar chirp rate, c is the speed of light, fs is the sampling frequency and F is said estimated normalized frequency calculated according to any one of previous claims from 1 to 7.
9. Method for angle of arrival (Q) calculation of a point target having range (R) detected through a multiple-input multiple- output frequency modulated continuous wave radar system equipped with a plurality of transmitting (TX) and receiving (RX) antennas, wherein each couple of said TX and RX antennas is replaced with the equivalent virtual antenna; said virtual antennas forming a virtual array composed by one uniform linear array; said angle of arrival being computed as:
Q= arcsin
here, l is the radar wavelength, d is the distance between two adjacent virtual antennas of said uniform linear array, and F is said estimated normalized frequency according to any one of the claims from 1 to 7.
10. Method for computer estimation of the delay, amplitude and phase of a single complex tone defining a signal acquired by means of a receiver; the method including:
- (Step 1: IFFT processing) estimation of the impulse response (X0) characterizing the communication channel and its first
— (v) NFFT — 1 derivatives ({Xfc ;k = 1,.. NFFT — 1}) with NFFT = 4, through IFFT calculation, with an oversampling factor (M ),
- (Step 2) iterative calculation of an estimate of the parameters of the most dominant tone, including phase (y), amplitude (a) and delay (t) on the basis of said impulse response, and its derivatives, through NFFT — 1 IFFTs.
11. Method according to claim 10, wherein said (Step 2) calculation includes:
- (Step 2.5.0) representation of a tone normalized delay T = t Af , where Af is the frequency step size of the employed radar, as sum of a main portion (OCFidft) and a residual portion (5Fidft), where FIDFT = l/iV0, where iV0 is the IFFT order, wherein said main portion is determined by searching the index (a) of the vector element Xoa corresponding to the maximum absolute value of the same vector X0 and said residual portion is found through said iterative calculation together with the corresponding target complex amplitude A = a exp(-jxp), where a is said target amplitude and y is said target phase.
12. [CSDE#1] Method according to claim 11, wherein said iterative procedure includes a preliminary initialization (Step 2.5.1) comprising: setting the iteration index i to 1, the initial estimate of A equal to said vector element (Xo,a) / the initial estimate of A = 27Z1>FIDFT to zero and the initial estimate of T to said main portion (aFIDFT); then iteratively
- (Step 2.5.2) the i-th iteration is fed by the estimates of D, T and A computed in the previous iteration, to generate new estimates of D , T and A , for a predetermined number of iterations, wherein each iteration includes a) a first step where the residual D is calculated solving the first order equation OR the second order equation
in the variable D with:
where ¾{x} denotes the real part and 3{x} denotes the imaginary part of the complex quantity x ; discarding one of the two solutions, which does not satisfy the inequalities: if a = 0
otherwise if said second order equation is solved; then the new estimate of T is obtained by adding said residual to said main
portion (aFIDFT); b) a second step wherein said new estimate of T is employed to compute the new estimate of and is
generated by interpolating the elements of said vector where
M is said oversampling factor employed in said IFFT calculation.
13. [CSDE#2] Method according to claim 11, wherein said iterative procedure includes a preliminary initialization (Step 2.6.1) comprising: a) setting the iteration index i to 1, the initial estimate of A equal to said vector element the initial estimate of D =
2i FIDFT to zero and the initial estimate of T to said main portion
b) computing the real quantity
then iteratively
- (Step 2.6.2) the i-th iteration is fed by the estimates of D, T, A and the real quantity p computed in the previous iteration, to generate new estimates of D, T, A and p, for a predetermined number of iterations, wherein each iteration includes:
a) a first step where the residual D is calculated solving the first order equation OR the second order equation
:
where ¾{x} denotes the real part, 3{x} denotes the imaginary part of the complex quantity x and {Xkp;k = 1,2,3} is computed as:
where xcn is the n-th element of the said complex signal; said coefficients ({Xk,p>k = 1,2,3}) can be also computed by interpolating the samples of said impulse response (Xfc); discarding one of the two solutions, which does not satisfy the inequalities: ifa = 0
otherwise if said second order equation is solved; then the new estimate of T is obtained by adding said residual to said estimate
of T; b) a second step wherein the new estimate of p is computed by dividing said new estimate of T by Fidft; c) a third step wherein said new estimate of T is employed to compute the new estimate of A as A = Xint
generated by interpolating the elements of said vector where
M is said oversampling factor employed in said IFFT calculations.
14. Method for range calculation of a point target detected
through a single-input single-output stepped frequency continuous wave radar system as
where D/ is the frequency step size of the employed radar system, c is the speed of light and T is said estimated normalized delay, calculated according to any one of the previous claims from 10 to 13.
15. Method for angle of arrival (Q) calculation of a point target having range ( R) detected through a multiple-input multiple- output stepped frequency continuous wave radar system equipped with a plurality of transmitting (TX) and receiving (RX) antennas, wherein each couple of said TX and RX antennas is replaced with the equivalent virtual antenna; said virtual antennas forming a virtual array composed by one uniform linear array; said angle of arrival being computed as:
here, l is the radar wavelength, d is the distance between two adjacent virtual antennas of said uniform linear array, and T is said estimated normalized delay according to any one of the claims 10 - 13.
16. A MIMO FMCW radar system including:
- a MIMO FMCW radar equipped with a plurality of transmitting (TX) and receiving (RX) antennas, arranged to generate real or complex signals in response to a propagation scenario including a plurality of point targets, and
- computing means arranged to receive signals generated by said transmitting and receiving antennas and configured to replace each couple of said transmitting and receiving antennas with an equivalent model defining virtual antennas and to perform all the steps of any one of the previous claims from 1 to 9.
17. A SISO FMCW radar system including:
- a SISO FMCW radar equipped with a single transmitting (TX) and receiving (RX) antennas, arranged to generate real or complex signals in response to a propagation scenario including a plurality of point targets, and
- computing means arranged to receive signals generated by said receiving antenna and configured to perform all the steps of any one of the previous claims from 1 to 9.
18. A MIMO SFCW radar system including:
- a MIMO SFCW radar equipped with a plurality of transmitting (TX) and receiving (RX) antennas, arranged to generate complex signals in response to a propagation scenario including a plurality of point targets, and
- computing means arranged to receive signals generated by said transmitting and receiving antennas and configured to replace
each couple of said transmitting and receiving antennas with an equivalent model defining virtual antennas and to perform all the steps of any one of the previous claims from 10 to 15.
19. A SISO SFCW radar system including:
- a SISO SFCW radar equipped with a single transmitting (TX) and receiving (RX) antennas, arranged to generate complex signals in response to a propagation scenario including a plurality of point targets, and - computing means arranged to receive signals generated by said receiving antenna and configured to perform all the steps of any one of the previous claims from 10 to 15.
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
IT102021000011156A IT202100011156A1 (en) | 2021-05-03 | 2021-05-03 | Method for two- and three-dimensional imaging based on co-located multiple inputs and multiple outputs radar |
IT202100011243 | 2021-05-03 | ||
PCT/EP2022/061876 WO2022233889A1 (en) | 2021-05-03 | 2022-05-03 | Methods for computer estimation of a single tone and their application to radar systems |
Publications (1)
Publication Number | Publication Date |
---|---|
EP4334742A1 true EP4334742A1 (en) | 2024-03-13 |
Family
ID=81849297
Family Applications (3)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
EP22726140.1A Pending EP4334740A1 (en) | 2021-05-03 | 2022-05-03 | Method for two-dimensional and three-dimensional imaging based on colocated multiple-input multiple-output radars |
EP22726141.9A Pending EP4334741A1 (en) | 2021-05-03 | 2022-05-03 | Method for two-dimensional and three-dimensional imaging based on collocated multiple-input multiple-output radars |
EP22726142.7A Pending EP4334742A1 (en) | 2021-05-03 | 2022-05-03 | Methods for computer estimation of a single tone and their application to radar systems |
Family Applications Before (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
EP22726140.1A Pending EP4334740A1 (en) | 2021-05-03 | 2022-05-03 | Method for two-dimensional and three-dimensional imaging based on colocated multiple-input multiple-output radars |
EP22726141.9A Pending EP4334741A1 (en) | 2021-05-03 | 2022-05-03 | Method for two-dimensional and three-dimensional imaging based on collocated multiple-input multiple-output radars |
Country Status (2)
Country | Link |
---|---|
EP (3) | EP4334740A1 (en) |
WO (3) | WO2022233889A1 (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP4418000A1 (en) * | 2023-02-17 | 2024-08-21 | IMRA Europe S.A.S. | Radar data processing by dnn |
CN117991245B (en) * | 2024-02-04 | 2024-07-23 | 南京畅淼科技有限责任公司 | Ship azimuth estimation method based on combination of multi-channel secondary DFT and interpolation |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006079181A1 (en) * | 2005-01-31 | 2006-08-03 | Genesys Design Pty Ltd | Frequency estimation |
US11346936B2 (en) * | 2018-01-16 | 2022-05-31 | Infineon Technologies Ag | System and method for vital signal sensing using a millimeter-wave radar sensor |
DE102018119278B4 (en) * | 2018-08-08 | 2020-04-02 | Infineon Technologies Ag | METHOD AND DEVICES FOR PROCESSING AN OFDM RADAR SIGNAL |
-
2022
- 2022-05-03 EP EP22726140.1A patent/EP4334740A1/en active Pending
- 2022-05-03 WO PCT/EP2022/061876 patent/WO2022233889A1/en active Application Filing
- 2022-05-03 WO PCT/EP2022/061875 patent/WO2022233888A1/en active Application Filing
- 2022-05-03 EP EP22726141.9A patent/EP4334741A1/en active Pending
- 2022-05-03 EP EP22726142.7A patent/EP4334742A1/en active Pending
- 2022-05-03 WO PCT/EP2022/061874 patent/WO2022233887A1/en active Application Filing
Also Published As
Publication number | Publication date |
---|---|
WO2022233888A1 (en) | 2022-11-10 |
WO2022233887A1 (en) | 2022-11-10 |
EP4334741A1 (en) | 2024-03-13 |
EP4334740A1 (en) | 2024-03-13 |
WO2022233889A1 (en) | 2022-11-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5600866B2 (en) | Detecting and ranging apparatus and detecting and ranging method | |
Baral et al. | Joint Doppler frequency and direction of arrival estimation for TDM MIMO automotive radars | |
EP4334742A1 (en) | Methods for computer estimation of a single tone and their application to radar systems | |
Oh et al. | Low-complexity range-azimuth FMCW radar sensor using joint angle and delay estimation without SVD and EVD | |
Aubry et al. | Experimental analysis of block-sparsity-based spectrum sensing techniques for cognitive radar | |
Yin et al. | Direct localization of multiple stationary narrowband sources based on angle and Doppler | |
Qin et al. | A decoupled direct positioning algorithm for strictly noncircular sources based on Doppler shifts and angle of arrival | |
Demissie | Clutter cancellation in passive radar using GSM broadcast channels | |
Yin et al. | Single-step localization using multiple moving arrays in the presence of observer location errors | |
Sirignano et al. | A comparative analysis of deterministic detection and estimation techniques for MIMO SFCW radars | |
Di Viesti et al. | Novel methods for approximate maximum likelihood estimation of multiple superimposed undamped tones and their application to radar systems | |
Shoudha et al. | Reduced-complexity decimeter-level Bluetooth ranging in multipath environments | |
Di Viesti et al. | Novel deterministic detection and estimation algorithms for colocated multiple-input multiple-output radars | |
Zakariyya et al. | Joint DoA and CFO estimation scheme with received beam scanned leaky wave antenna for Industrial Internet of Things (IIoT) systems | |
Nuss et al. | 3D radar image fusion using OFDM-based MIMO radar | |
Gu et al. | System architecture and signal processing for frequency-modulated continuous-wave radar using active backscatter tags | |
Feger et al. | Sparse antenna array design and combined range and angle estimation for FMCW radar sensors | |
Ulrich et al. | Wavelength-diverse MIMO radar: parameter-coupling, array-carrier optimization and direction-of-arrival estimation | |
Storrer et al. | Experimental implementation of a multi-antenna 802.11 ax-based passive bistatic radar for human target detection | |
Zhang et al. | Joint range and velocity fast estimation for OFDM-based RadCom shared signal | |
Zeng et al. | Delay compensation for distributed MIMO radar with non-orthogonal waveforms | |
Hashempour et al. | Cyclic prefix-based OFDM ISAR imaging | |
Hussain et al. | LDL decomposition-based FGPA real-time implementation of DOA estimation | |
Ding et al. | Integrated Sensing and Communication: Prototype and Key Processing Algorithms | |
Wang et al. | Performance bounds for displaced sensor automotive radar imaging |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: UNKNOWN |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: THE INTERNATIONAL PUBLICATION HAS BEEN MADE |
|
PUAI | Public reference made under article 153(3) epc to a published international application that has entered the european phase |
Free format text: ORIGINAL CODE: 0009012 |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: REQUEST FOR EXAMINATION WAS MADE |
|
17P | Request for examination filed |
Effective date: 20231204 |
|
AK | Designated contracting states |
Kind code of ref document: A1 Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR |
|
DAV | Request for validation of the european patent (deleted) | ||
DAX | Request for extension of the european patent (deleted) |