The present invention relates generally to signal processing and analog to digital conversion, and particularly to methods and systems for sub-Nyquist sampling and processing of pulsed signals.
Various methods and systems for sub-Nyquist signal processing are known in the art. For example, U.S. Pat. No. 8,457,579, whose disclosure is incorporated herein by reference, describes a method for signal processing that includes distributing an analog input signal to a plurality of processing channels. In each processing channel, the input signal is mixed with a respective periodic waveform including multiple spectral lines, so as to produce a respective baseband signal in which multiple spectral slices of the input signal are superimposed on one another. The baseband signal produced in each of the processing channels is digitized, to produce a set of digital sample sequences that represent the input signal.
U.S. Pat. No. 8,032,085, whose disclosure is incorporated herein by reference, describes a signal processing method that includes sampling an analog signal, which has a spectral density defining one or more bands, to produce a digitized signal. A spectral transform of the digitized signal is expressed as a matrix multiplied by a vector, wherein the vector includes multiple elements that represent the spectral density of the analog signal as a function of frequency within respective spectral slices. Indices of a subset of the elements of the vector, in which the spectral density is concentrated, is determined. The analog signal is reconstructed from the digitized signal using the subset of the elements of the vector and a corresponding subset of columns of the matrix having the indices.
U.S. Patent Application Publication 2011/0225218, whose disclosure is incorporated herein by reference, describes a method that includes accepting an analog input signal that includes a sequence of pulses. The analog input signal is filtered so as to produce a filter output, using a filter whose time-domain response is confined to a finite time period and whose frequency-domain response is non-zero at a finite set of integer multiples of a frequency shift Δω, and is zero at all other integer multiples of Δω. The filter output is sampled so as to produce digital samples. Respective amplitudes and time positions of the pulses in the sequence are calculated based on the digital samples.
U.S. Patent Application Publication 2013/0038479, whose disclosure is incorporated herein by reference, describes a method that includes accepting an analog input signal including a sequence of pulses of a given pulse shape. The analog input signal is distributed to multiple processing channels operating in parallel. The analog input signal is sampled by performing, in each of the multiple processing channels, the operations of: mixing the analog input signal with a different, respective modulating waveform to produce a mixed signal; filtering the mixed signal; and digitizing the filtered mixed signal to produce a respective digital channel output.
U.S. Patent Application Publication 2013/0187682, whose disclosure is incorporated herein by reference, describes a method for signal processing that includes accepting an analog signal, which consists of a sequence of pulses confined to a finite time interval. The analog signal is sampled at a sampling rate that is lower than a Nyquist rate of the analog signal and with samples taken at sample times that are independent of respective pulse shapes of the pulses and respective time positions of the pulses in the time interval.
An embodiment of the present invention that is described herein provides a method including receiving a signal, which includes reflections of multiple pulses from one or more targets. A Doppler focusing function, in which the reflections of the multiple pulses from each target accumulate in-phase to produce a respective peak associated with a respective delay and a respective Doppler frequency of the target, is evaluated based on the received signal. Respective delays and Doppler frequencies of the targets are estimated based on the Doppler focusing function.
In some embodiments, receiving the signal includes sampling the signal with a sampling rate that is lower than a Nyquist rate of the pulses, and evaluating the Doppler focusing function includes processing the sampled signal at a processing rate that is lower than the Nyquist rate. In a disclosed embodiment, estimating the delays and the Doppler frequencies includes defining multiple Doppler-frequency bands and, for each band in at least a subset of the bands, estimating the delays for one or more of the targets whose Doppler frequency falls in the band.
In an embodiment, estimating the delays and the Doppler frequencies includes identifying peaks of the Doppler focusing function, and deriving the delays and the Doppler frequencies from the identified peaks. In another embodiment, receiving the signal includes sampling the signal, and estimating the delays and the Doppler frequencies includes identifying a given peak in the Doppler focusing function, deriving a delay and a Doppler frequency of a target from the given peak, adjusting the sampled signal to remove a contribution of the target, re-evaluating the Doppler focusing function based on the adjusted sampled signal, and searching for a next peak in the re-evaluated Doppler focusing function.
In yet another embodiment, evaluating the Doppler focusing function includes calculating a Fourier transform coefficient matrix over a uniform grid of the Doppler frequencies. In still another embodiment, evaluating the Doppler focusing function includes windowing the Doppler focusing function using a window function.
There is additionally provided, in accordance with an embodiment of the present invention, apparatus including a sampling unit and digital processing circuitry. The sampling unit is configured to sample a received signal that includes reflections of multiple pulses from one or more targets. The digital processing circuitry is configured to evaluate, based on the sampled signal, a Doppler focusing function in which the reflections of the multiple pulses from each target accumulate in-phase to produce a respective peak associated with a respective delay and a respective Doppler frequency of the target, and to estimate respective delays and Doppler frequencies of the targets based on the Doppler focusing function.
There is further provided, in accordance with an embodiment of the present invention, apparatus including a receiver and a processor. The receiver is configured to receive a signal that includes reflections of multiple pulses from one or more targets. The processor is configured to evaluate, based on the received signal, a Doppler focusing function in which the reflections of the multiple pulses from each target accumulate in-phase to produce a respective peak associated with a respective delay and a respective Doppler frequency of the target, and to estimate respective delays and Doppler frequencies of the targets based on the Doppler focusing function.
The present invention will be more fully understood from the following detailed description of the embodiments thereof, taken together with the drawings in which:
Embodiments of the present invention that are described herein provide improved methods and systems for processing radar signals. In the disclosed embodiments, a radar system transmits a signal comprising a train of pulses, and receives a reflected signal comprising reflections of the pulses from one or more targets. The system then samples and processes the received signal, so as to reconstruct the amplitude, delay and Doppler frequency associated with each target. In some embodiments, all sampling and reconstruction are performed at a sub-Nyquist rate, i.e., below the Nyquist rate of the pulses that is typically given by twice the highest signal frequency.
In some embodiments, the system uses the samples of the received signal to evaluate a Doppler focusing function. The Doppler focusing function is a two-dimensional function, in which reflections of the pulses from each target accumulate in phase to produce a peak that is indicative of the target's Doppler frequency and delay. The Doppler focusing function may be evaluated either in the time domain or in the frequency domain. The system then derives the amplitude, delay and Doppler frequency of each target from the maxima Doppler focusing function. The derivation of target parameters is performed either directly or in an iterative fashion repeating the process one target at a time.
The use of the Doppler focusing function provides high immunity to noise and clutter, as well as fine Doppler resolution. Since the entire sampling and parameter reconstruction process is performed at a low, sub-Nyquist rate, the cost, size and power consumption of the radar system can be significantly reduced.
In the embodiment of
Both sampling and subsequent signal processing are performed at a sub-Nyquist rate, i.e., at a rate that is considerably lower than the bandwidth of the radar signal. In the present example, processor 36 comprises a sub-Nyquist sampling unit 40, a Doppler focusing unit 44 and a parameter reconstruction unit 48. Units 44 and 48 are sometimes referred to herein jointly as digital processing circuitry, which carries out the digital signal processing functions described herein. The functions of the various units of processor 36 are explained in detail below.
The system configuration of
In some embodiments, some or all of the functions of processor 36 may be implemented using a general-purpose computer, which is programmed in software to carry out the functions described herein. The software may be downloaded to the processor in electronic form, over a network, for example, or it may, alternatively or additionally, be provided and/or stored on non-transitory tangible media, such as magnetic, optical, or electronic memory.
In the embodiments described herein, the radar signal transmitted by transceiver 32 comprises a series of P equally-spaced pulses h(t). The transmitted signal can thus be written as
wherein τ denotes the pulse-to-pulse delay, or Pulse Repetition Interval (PRI).
Assuming the transmitted signal is reflected by a set of L targets 24, the received signal can be written as
wherein αl, τl, νl denote the amplitude, delay and Doppler frequency of the lth target, respectively. The 3L parameters αl, τl, νl l=0 . . . L−1 are referred to herein collectively as target parameters. Generally, αl, τl, νl are indicative of the target Radar Cross Section (RCS), range and radial velocity, respectively.
It is possible to express the received signal as a sum of frames
The radar's unambiguous time-frequency region is [0,τ]×[−π/τ,π/τ].
In some embodiments, processor 36 evaluates a Doppler focusing function of the received signal, which is defined as
Consider the sum-of-exponents Σp=0P-1ej(ν-ν
Thus, it is possible to approximate Φ(t; ν) by the expression
wherein Λ(ν)={l:|ν−νl|<π/Pτ}, i.e., the set of targets inside the focus zone.
A similar Doppler focusing process can be defined and evaluated in the frequency domain. Transforming Equation [4] above to the frequency domain yields
wherein Xp(ω) denotes the continuous-time Fourier transform (CTFT) of xp(t+pτ).
Taking the CTFT of Φ(t; ν) as a function of t yields
Similarly to Φ(t; ν), Ψ(ω; ν) is also expressed in terms of the 3L target parameters, namely the amplitudes, delays and Doppler frequencies of the L targets. The expression of Ψ(ω; ν) also comprises the sum-of-exponents term Σp=0P-1ej(ν-ν
The Doppler focusing function can be evaluated in practice from properly chosen sub-Nyquist samples. The description that follows shows how this function can be used for target estimation and detection. These operations are typically performed on the sub-Nyquist samples of the signal. Sub-Nyquist sampling schemes that can be used for obtaining such samples are described, for example, in U.S. Patent Application Publication 2013/0038479, cited above; and in “A Sub-Nyquist Radar Prototype: Hardware and Algorithms,” Baransky et al., arXiv:1208.2515 [cs.IT], Aug. 13, 2012, which is incorporated herein by reference.
The Doppler focusing function comprises a two-dimensional function, which has local maxima that are indicative of targets. In other words, detecting a peak in Φ(τT; νT) is indicative of a target having a delay is τT and Doppler frequency νT. This representation enables processor 36 to reduce the delay-Doppler estimation problem (estimating delays and Doppler frequencies of multiple targets over the system's unambiguous time-Doppler range) into a set of delay-only estimation problems (estimating the delays of multiple targets in a selected narrow slice of Doppler frequencies). This reduction can be carried out either in the time domain using Φ(t; ν), or in the frequency domain using Ψ(ω; ν).
In some embodiments, the target parameters are extracted directly from the Doppler focusing function. In alternative embodiments, the target parameters are extracted in an iterative process, one target at a time, while adapting the Doppler focusing function to remove the contribution of each target that has already been detected. These two alternatives are described in detail below. Further alternatively, the target parameters can be extracted from the Doppler focusing function in any other suitable way.
After reducing the problem dimensionality, processor 36 solves each delay-only estimation problem at a sub-Nyquist rate. Processor 36 may use various techniques for estimating the target delays for a selected Doppler frequency. Example techniques for reconstructing the 2L amplitudes and delays in a signal of the form
from sub-Nyquist samples are given, for example, in U.S. Patent Application Publication 2013/0038479, cited above. Such techniques are also described by Gedalyahu et al., in “Multichannel Sampling of Pulse Streams at the Rate of Innovation,” IEEE Transactions on Signal Processing, volume 59, number 4, April, 2011, pages 1491-1504; by Vetterli et al., in “Sampling Signals with Finite Rate of Innovation,” IEEE Transactions on Signal Processing, volume 50, number 6, 2002, pages 1417-1428; by Tur et al., in “Innovation Rate Sampling of Pulse Streams with Application to Ultrasound Imaging,” IEEE Transactions on Signal Processing, volume 59, number 4, April, 2011, pages 1827-1842; and by Wagner et al., in “Compressed Beamforming in Ultrasound Imaging,” IEEE Transactions on Signal Processing, volume 60, issue 9, September, 2012, pages 4643-4657, which are all incorporated herein by reference, as well as in the article by Baransky et al., cited above. Processor 36 may use any such technique, or any other suitable technique, for estimating the target delays for a given Doppler frequency at a sub-Nyquist rate.
When using these techniques, processor 36 typically pre-filters the signal ϕ(t) and then samples the pre-filtered signal, so as to obtain Fourier coefficients c [k] of the Fourier expansion of ϕ(t):
In some embodiments, processor 36 (e.g., sampling unit 40) produces a set of coefficients c[k] by splitting the received signal into multiple processing channels, modulating the signal in each channel by a respective complex sinusoid, and integrating the result over [0, τ]. A detailed description of such sampling schemes is given in U.S. Patent Application Publication 2013/0038479, cited above.
In alternative embodiments, processor 36 filters the signal with several band-pass filters, and then samples the output of each filter with a low sampling rate. Sampling of this sort is described in the article by Baransky et al., cited above. Further alternatively, processor 36 may perform sub-Nyquist sampling using multi-channel filters or multi-channel modulators. All of these techniques result in low-rate samples from which the Fourier coefficients c[k] can be obtained over a desired set of frequencies that is lower than the signal bandwidth.
From Equation [11] it can be seen that the target parameters αl, τl are embedded in the Fourier coefficients c[k] in the form of a complex sinusoid problem. Without noise, it can be shown that for such a problem, 2L coefficients are sufficient for reconstructing αl, τl l=0 . . . L−1.
Many reconstruction techniques for solving the sinusoid problem and reconstructing the coefficients are known in the art, for example annihilating filters, matrix pencil, ESPRIT and MUSIC, among others. MUSIC is described by Schmidt, in “Multiple Emitter Location and Signal Parameter Estimation,” IEEE Transactions on Antennas and Propagation, volume 34, no. 3, March, 1986, pages 276-280, and by Bienvenu and Kopp, in “Adaptivity to Background Noise Spatial Coherence for High Resolution Passive Methods,” IEEE International Conference on Acoustics, Speech and Signal Processing, volume 5, April, 1980, pages 307-310, which are incorporated herein by reference. ESPRIT, and is described by Roy and Kailath, in “ESPRIT—Estimation of Signal Parameters via Rotational Invariance Techniques,” IEEE Transactions on Acoustics, Speech and Signal Processing, volume 37, no. 7, July, 1989, pages 984-995, which is incorporated herein by reference.
In other embodiments, an alternative compressed sensing scheme can be used for recovering the target delays from coefficients c[k] with improved performance in the presence of noise. Assume that the time delays are aligned to a grid
τl=nlΔτ 0≤nl≤Nτ [12]
wherein Δτ is chosen such that Nτ=τ/Δτ is integer. If the time delays do not lie exactly on a grid, local interpolation may be used to reduce quantization errors. For a set of indices κ={k0, . . . , k|κ|-1}, we define the corresponding vector of Fourier coefficients
c=[c[k0] . . . c[k|K|-1]]T [13]
Equation [11] above can be written in vector form as
wherein H is a |κ|×|κ| diagonal matrix with diagonal elements H (2πk/τ), and V is a |κ|×Nτ Vandermonde matrix with elements Vmq=e−j2πk
The vector of target delays x is L-sparse, with the nth vector element gives the amplitude of a target with delay nΔτ if such a target exists, or zero otherwise. Defining the Compressed Sensing (CS) dictionary
we get
c=Ax [15]
The target delays can be reconstructed by solving Equation [15] and finding the support of x, i.e., the set of indices of the non-zero elements in x. Any suitable CS technique can be used for solving Equation [15]. Example CS techniques are described by Eldar and Kutyniok, in “Compressed Sensing: Theory and Applications,” Cambridge University Press, 2012, which is incorporated herein by reference.
The description now returns to the problem of reconstructing the target delays and amplitudes, for a given Doppler frequency, from the sub-Nyquist-sampled coefficients c[k]. Let cp[k] denote the Fourier coefficients in the Fourier expansion of the component of the received signal contributed by the pth pulse:
Similarly to Equation [14], we can write for each pulse p:
From Equation [16] it can be seen that all 3L target parameters αl, τl, νl are embedded in the Fourier coefficients cp[k] in the form of a complex sinusoid problem. The number of Fourier coefficients sampled for each pulse, |Λ|, sets the trade-off between sampling rate and robustness to noise.
Thus, in some embodiments sampling unit 40 acquires Fourier coefficients cp[k] with a sub-Nyquist sampling rate. Based on this set of coefficients, Doppler focusing unit 44 evaluates the Doppler focusing function in the frequency domain:
As discussed above, for each target 1 whose Doppler frequency νl is within the focus zone of the selected Doppler frequency ν (i.e., |νl−ν|<π/Pτ), we get
Therefore, Doppler focusing can be performed on the low rate sub-Nyquist samples:
Thus, for any Doppler frequency ν we get a delay-only estimation problem of the form
wherein xν is L-sparse and
Ψν=[Ψν[k0] . . . Ψν[k|k|-1]]T [22]
This problem is identical to the CS problem of Equations [14]-[15] above. Note that, unlike other CS schemes that estimate delay and Doppler frequency, in the disclosed technique the CS dictionary does not grow with the number of pulses.
The Doppler focusing function of Equation [18] above is a continuous function of ν, and can be evaluated for any Doppler frequency up to the Pulse Repetition Frequency (PRF, defined as 1/PRI). The Doppler focusing function has no “blind speeds,” i.e., target velocities that are undetectable.
Let
denote the set of above-described Fourier coefficients, and let Ψν(C) denote the vector of Doppler-focused coefficients of equation [22] that were obtained from C using Equation [18]. Therefore, xν(C) can be recovered from Ψν(C) for any ν. Typically, processor 36 searches for large values of |xν(C)[n]|, which are indicative of targets. For each identified peak, processor 36 estimates the target delays and Doppler frequencies as nΔτ and ν, respectively.
In some embodiments, processor 36 runs an iterative process that finds a peak associated with a target, estimates the target parameters, adjusts the set of Fourier coefficients so as to remove the contribution of this target, and re-evaluates the Doppler focusing function using the adjusted Fourier coefficients. This sort of iterative process improves the capability to detect weak targets and removes spurious targets that may be generated by processing side-lobes.
The iterative process continues until some termination criterion is met. For example, if the number of targets L is known, processor 36 may continue the iterations until L targets are found. If the number of targets is unknown, the process may terminate when the identified peak drops below some amplitude threshold. The iterative process is summarized by the following example pseudo-code:
As noted above, this process is depicted purely by way of example. In alternative embodiments, processor 36 may recover the target parameters from the Doppler focusing function using any other suitable process, whether iterative or non-iterative.
In the above process, although the total number of targets is known, the distribution of the targets in the delay-Doppler plane is typically unknown. Therefore, when solving Equation [21] inside the above loop, processor 36 may estimate the model order 0≤Lν≤L, or make a worst-case assignment Lν=L. Processor 36 may alternatively stop the iterations upon detecting that the progress of the iterative process (e.g., reduction of error between successive iterations) falls below some predefined threshold.
In various embodiments, processor 36 may solve Equation [21] using any suitable method, e.g., one of many known methods for estimating the number of sinusoids in a noisy sequence. Solving Equation [21] with a accurate model order may decrease computation time, and may reduce detection of false targets. This option may be preferable in high Signal-to-Noise Ratio (SNR) scenarios when the Lν estimates are accurate. When SNR is poor, the worst case approach that assumes Lν=L may perform better.
Since V is a partial Fourier matrix, the problem defined in Equation [21], after normalization by H−1, or any other equalizer such as a matched filter, Wiener filter etc., is equivalent to the problem of recovering the individual frequencies from a sum of complex exponentials. Many techniques are known for solving this problem, and processor 36 may use any such method instead of CS. CS may be preferable for low-SNR scenarios. However, since Doppler focusing is independent of the underlying delay estimation, in various scenarios CS can be replaced with alternative delay estimation methods. For example, when noise is not a dominant factor, processor 36 may use an annihilating filter or MUSIC method instead of CS.
In some embodiments, processor 36 evaluates the argmax term above over a uniform grid of M Doppler frequencies. In these embodiments, it is possible to evaluate the Doppler focusing function efficiently using a length-M DFT or Fast Fourier Transform (FFT) of a series of length P:
Ψm[k]=DFTM{cp[k]} [23]
The following pseudo-code carries out Doppler focusing-based parameter estimation for the case of Doppler frequencies lying on a uniform grid:
In some embodiments, processor 36 applies windowing when evaluating the Doppler focusing function. In an example embodiment, multiplies the expression of Equation [18] above by a suitable window function w[p], p=0, . . . , P−1, such as a Hann or Blackman window. Windowing helps to reduce the effect of targets that are nearby one another in the delay-Doppler plane. On the other hand, windowing increases the size of the focus zone. Windowing is important, for example, when the dynamic range of target amplitudes is large.
In processor 36, sampling unit 40 samples the received signal with a sub-Nyquist rate, at a sampling step 58. Doppler focusing unit 44 evaluates the Doppler focusing function based on the sub-Nyquist sampled signal, at a Doppler focusing step 62. Parameter reconstruction unit 48 reconstructs the amplitudes, delays and Doppler frequencies (corresponding to the RCS, range and radial velocity) of the targets from the Doppler focusing function, at a parameter reconstruction step 66. The reconstructed target parameters are provided as output of system 20, e.g., displayed to an operator using a suitable display.
Although the embodiments described herein mainly address radar applications, the methods and systems described herein can also be used in other applications that involve estimating signal delays and Doppler frequencies, such as in estimation of multipath mobile communication channels, ultrasound imaging, sonar, radar imaging (e.g., Synthetic Aperture Radar—SAR), or in any other type of receiver that coherently integrates signals having different Doppler frequencies.
It will thus be appreciated that the embodiments described above are cited by way of example, and that the present invention is not limited to what has been particularly shown and described hereinabove. Rather, the scope of the present invention includes both combinations and sub-combinations of the various features described hereinabove, as well as variations and modifications thereof which would occur to persons skilled in the art upon reading the foregoing description and which are not disclosed in the prior art. Documents incorporated by reference in the present patent application are to be considered an integral part of the application except that to the extent any terms are defined in these incorporated documents in a manner that conflicts with the definitions made explicitly or implicitly in the present specification, only the definitions in the present specification should be considered.
This application claims the benefit of U.S. Provisional Patent Application 61/722,453, filed Nov. 5, 2012, whose disclosure is incorporated herein by reference.
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/IB2013/059602 | 10/24/2013 | WO | 00 |
Publishing Document | Publishing Date | Country | Kind |
---|---|---|---|
WO2014/080303 | 5/30/2014 | WO | A |
Number | Name | Date | Kind |
---|---|---|---|
4106019 | Alexander | Aug 1978 | A |
5416488 | Grover | May 1995 | A |
6188348 | Raney | Feb 2001 | B1 |
6313785 | Mattox | Nov 2001 | B1 |
6809682 | Madewell | Oct 2004 | B1 |
6822606 | Ponsford | Nov 2004 | B2 |
7151483 | Dizaji | Dec 2006 | B2 |
7667639 | Cheng | Feb 2010 | B2 |
8013781 | Stockmann | Sep 2011 | B2 |
8032085 | Mishali et al. | Oct 2011 | B2 |
8457579 | Mishali et al. | Jun 2013 | B2 |
8842037 | Emery | Sep 2014 | B2 |
20030174088 | Dizaji | Sep 2003 | A1 |
20050179579 | Pinder | Aug 2005 | A1 |
20070040729 | Ohnishi | Feb 2007 | A1 |
20100049470 | Szajnowski | Feb 2010 | A1 |
20110225218 | Eldar et al. | Sep 2011 | A1 |
20120105276 | Ryland | May 2012 | A1 |
20120313808 | Frohling | Dec 2012 | A1 |
20130038479 | Eldar et al. | Feb 2013 | A1 |
20130187682 | Eldar et al. | Jul 2013 | A1 |
20130257645 | Penney | Oct 2013 | A1 |
Number | Date | Country |
---|---|---|
2012020530 | Feb 2012 | WO |
2012056357 | May 2012 | WO |
Entry |
---|
Parker, Michael, “Radar Basics—Part 2: Pulse Doppler Radar” Altera Corporation. May 28, 2011. Accessed Sep. 19, 2018. <https:// www.eetimes.com/document.asp?doc_id=1278808> (Year: 2011). |
T. L. Marzetta, E. A. Martinsen and C. P. Plum, “Fast pulse Doppler radar processing accounting for range bin migration,” The Record of the 1993 IEEE National Radar Conference, Lynnfield, MA, USA, 1993, pp. 264-268.doi: 10.1109/NRC.1993.270453 (Year: 1993). |
Tur et al., “Innovation Rate Sampling of Pulse Streams with Application to Ultrasound Imaging”, IEEE Transactions on Signal Processing, vol. 59, No. 4, pp. 1827-1842, Apr. 2011. |
Gedalyahu et al., “Multichannel Sampling of Pulse Streams at the Rate of Innovation,” IEEE Transactions on Signal Processing, vol. 59, No. 4, pp. 1491-1504, Apr. 2011. |
Baransky et al., “A Sub-Nyquist Radar Prototype: Hardware and Algorithms”, Department of Electrical Engineering, Technion-Israel Institute of Technology, Submitted to IEEE Transactions on Aerospace and Electronic Systems, Special Issue on Compressed Sensing for Radar, 14 pages, Aug. 2012. |
Vetterli et al., “Sampling signals with finite rate of innovation”, IEEE Transactions on Signal Processing, vol. 50, No. 6, pp. 1417-1428, Jun. 2002. |
Mishali et al., “Xampling: Signal Acquisition and Processing in Union of Subspaces”, IEEE Transactions on Signal Processing, vol. 59, No. 10, pp. 4719-4734, Oct. 2011. |
Mishali et al., “Xampling: Analog to Digital at Sub Nyquist Rates”, Technion—Israel Institute of Technology, Department of Electrical Engineering, Published in IET Circuits, Devices & Systems, vol. 5, No. 1, pp. 8-20, Jan. 2011. |
Sarkar et al., “Using the Matrix Pencil Method to Estimate the Parameters of a Sum of Complex Exponentials”, IEEE Antennas and Propagation Magazine, vol. 37, No. 1, pp. 48-55, Feb. 1995. |
Rudelson et al., “On Sparse Reconstruction from Fourier and Gaussian Measurements”, Communications on Pure and Applied Mathematics, vol. LXI, pp. 1025-1045, 2008. |
Herman et al., “High-resolution radar via compressed sensing”, IEEE Transactions on Signal Processing, vol. 57, No. 6, pp. 2275-2284, Jun. 2009. |
Baraniuk et al., “Compressive Radar Imaging”, IEEE Radar Conference, 6 pages, Apr. 2007. |
Bajwa et al., “Identification of Parametric Underspread Linear Systems and Super-Resolution Radar”, IEEE Transactions on Signal Processing, vol. 59, No. 6, pp. 2548-2561, Jun. 2011. |
Zhang et al., Adaptive Compressed Sensing Radar Oriented Toward Cognitive Detection in Dynamic Sparse Target Scene, IEEE Transactions on Signal Processing, vol. 60, No. 4, Apr. 2012. |
Wagner et al., “Compressed Beamforming in Ultrasound Imaging”, IEEE Transactions on Signal Processing, vol. 60, No. 9, pp. 4643-4657, Sep. 2012. |
Skolnik., “Radar Handbook”, 3rd Edition, The McGraw-Hill Companies, 1352 pages, year 2008. |
Bienvenu et al., “Adaptivity to Background Noise Spatial Coherence for High Resolution Passive Methods”,IEEE International Conference on Acoustics, Speech and Signal Processing, vol. 5, pp. 307-310, Apr. 1980. |
Schmidt et al., “Multiple Emitter Location and Signal Parameter—Estimation”, IEEE Transactions on Antennas and Propagation, vol. AP-34, No. 3, pp. 276-280, Mar. 1986. |
Roy et al., “ESPRIT—Estimation of Signal Parameters via Rotational Invariance Techniques”, IEEE Transactions on Acoustics Speech and Signal Processing, vol. 37, No. 7 pp. 984-995, Jul. 1989. |
Mallat et al., “Matching pursuits with time-frequency dictionaries”, IEEE Transactions on Signal Processing, vol. 41, No. 12, pp. 3397, 3415, Dec. 1993. |
Demissie, B., “High-Resolution Range-Doppler Imaging by Coherent Block-Sparse Estimation,” International Workshop on Compressed Sensing Applied to Radar, 1 page, May 2012—Abstract. |
Blumensath et al., “Iterative hard thresholding for compressed sensing”, Applied and Computational Harmonic Analysis, vol. 27, pp. 265-274, year 2009. |
Smith et al., “Compressed Sampling for Pulse Doppler Radar”, IEEE Radar Conference, 6 pages, May 10-14, year 2010. |
Schwarz, G., “Estimating the Dimension of a Model”, Institute of Mathematical Statistics,The Annals of Statistics, vol. 6, No. 2 , pp. 461-464, Mar. 1978. |
Wax et al., “Detection of Signals by Information Theoretic Criteria”, IEEE Transaction on Acoustics, Speech and Signal Processing, vol. ASSP-33, No. 2, pp. 387-392, Apr. 1985. |
Fuchs et al., “Estimating the Number of Sinusoids in Additive White Noise”, IEEE Transaction on Acoustics, Speech and Signal Processing, vol. 36, No. 12, pp. 1846-1853, Dec. 1988. |
Harris, F., “On the Use of Windows for Harmonic Analysis with the Discrete Fourier Transform”, Processing of the IEEE, vol. 66, No. 1, pp. 51-83, Jan. 1978. |
Duarte et al., “Structured Compressed Sensing: From Theory to Applications”., IEEE Transactions on Signal Processing, vol. 59, No. 9, pp. 4053-4085, Sep. 2011. |
Bar-Ilan et al., “Sub-Nyquist Radar via Doppler Focusing”, Department of Electrical Engineering, Technion—Israel Institute of Technology, 11 pages, Nov. 4, 2012. |
International application # PCT/IB2013/059602 Search Report dated Mar. 5, 2014. |
Stoica et al., “Introduction to Spectral Analysis”, Prentice-Hall, 447 pages, year 2000—Abstract. |
Number | Date | Country | |
---|---|---|---|
20150285905 A1 | Oct 2015 | US |
Number | Date | Country | |
---|---|---|---|
61722453 | Nov 2012 | US |