Multi-Carrier Agile Phased Array Radar
Tianyao Huang, Nir Shlezinger, Xingyu Xu, Dingyou Ma, Yimin Liu, and, Yonina C. Eldar

TL;DR
This paper introduces two multi-carrier frequency agile radar schemes, WMAR and CAESAR, which enhance target recovery in complex environments by leveraging wideband and narrowband waveforms, respectively, with CAESAR offering practical implementation advantages.
Contribution
The paper proposes two novel multi-carrier radar schemes, WMAR and CAESAR, extending FAR with improved performance and practical implementation considerations.
Findings
Both schemes outperform traditional FAR in simulations.
CAESAR's performance is close to wideband radar with simpler hardware.
The algorithms reliably recover target parameters under certain conditions.
Abstract
Modern radar systems are expected to operate reliably in congested environments. A candidate technology for meeting these demands is frequency agile radar (FAR), which randomly changes its carrier frequencies. FAR is known to improve the electronic counter-countermeasures (ECCM) performance while facilitating operation in congested setups. To enhance the target recovery performance of FAR in complex electromagnetic environments, we propose two radar schemes extending FAR to multi-carrier waveforms. The first is Wideband Multi-carrier Agile Radar (WMAR), which transmits/receives wideband waveforms simultaneously with every antenna. To mitigate the demanding hardware requirements associated with wideband waveforms used by WMAR, we next propose multi-Carrier AgilE phaSed Array Radar (CAESAR). CAESAR uses narrowband monotone waveforms, thus facilitating ease of implementation of the system,…
| Characters | Frequency agility | Beam pattern | # of observation | Transmit bandwidth | Receive bandwidth |
|---|---|---|---|---|---|
| WMAR | Yes | Focused | Moderate | Large | Large |
| CAESAR | Yes | Moderately focused | Moderate | Small | Small |
| FAR | Yes | Focused | Small | Small | Small |
| FDMA-MIMO | No | Omnidirectional | Large | Small | Large |
| SUMMeR | Yes | Omnidirectional | Moderate | Small | Large |
| FDA | No | Omnidirectional | Large | Small | Small |
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Multi-Carrier Agile Phased Array Radar
Tianyao Huang, Nir Shlezinger, Xingyu Xu, Dingyou Ma, Yimin Liu, and Yonina C. Eldar Parts of this work were presented in the 2018 IEEE International Workshop on Compressed Sensing applied to Radar, Multimodal Sensing, and Imaging (CoSeRa). This work received funding from the National Natural Science Foundation of China under Grants 61571260 and 61801258, from the European Union’s Horizon 2020 research and innovation program under grant No. 646804-ERC-COG-BNYQ, and from the Air Force Office of Scientific Research under grant No. FA9550-18-1-0208. T. Huang, Y. Liu, X. Xu, and D. Ma are with the EE Department, Tsinghua University, Beijing, China (e-mail: {huangtianyao, yiminliu}@tsinghua.edu.cn; {xy-xu15, mdy16}@mails.tsinghua.edu.cn).
N. Shlezinger and Y. C. Eldar are with the Faculty of Math and CS, Weizmann Institute of Science, Rehovot, Israel (e-mail: [email protected]; [email protected]).
Abstract
Modern radar systems are expected to operate reliably in congested environments. A candidate technology for meeting these demands is frequency agile radar (FAR), which randomly changes its carrier frequencies. FAR is known to improve the electronic counter-countermeasures (ECCM) performance while facilitating operation in congested setups. To enhance the target recovery performance of FAR in complex electromagnetic environments, we propose two radar schemes extending FAR to multi-carrier waveforms. The first is Wideband Multi-carrier Agile Radar (WMAR), which transmits/receives wideband waveforms simultaneously with every antenna. To mitigate the demanding hardware requirements associated with wideband waveforms used by WMAR, we next propose multi-Carrier AgilE phaSed Array Radar (CAESAR). CAESAR uses narrowband monotone waveforms, thus facilitating ease of implementation of the system, while introducing spatial agility. We characterize the transmitted and received signals of the proposed schemes, and develop an algorithm for recovering the targets, based on concepts from compressed sensing to estimate the range-Doppler parameters of the targets. We then derive conditions which guarantee their accurate reconstruction. Our numerical study demonstrates that both multi-carrier schemes improve performance compared to FAR while maintaining its practical benefits. We also demonstrate that the performance of CAESAR, which uses monotone waveforms, is within a small gap from the wideband radar.
Index Terms:
Frequency agile radar, multi-carrier agility, compressed sensing
I Introduction
Modern radars must be reliable, but at the same time compact, flexible, robust, and efficient in terms of cost and power usage [1, 2, 3, 4, 5]. A possible approach to meet these requirements is by exploiting frequency agility [1], namely, to utilize narrowband waveforms, while allowing the carrier frequencies to vary between different radar pulses. Among the main advantages of frequency agile radar (FAR) are its excellent electronic counter-countermeasures (ECCM) and electromagnetic compatibility (EMC) performance [1], and the fact that it has the flexibility of supporting spectrum sharing [4]. Furthermore, FAR is compatible with phased array antennas. Finally, by utilizing narrowband signals with varying frequencies, FAR systems can synthesize a large bandwidth with narrowband waveforms [6, 2], which simplifies the implementation of the waveform generator, facilitates the receiver operation, and allows the usage of non-linear amplifiers without limiting their power efficiency.
A major drawback of FAR compared to wideband radar is its reduced range-Doppler reconstruction performance of targets. This reduced performance is a byproduct of the relatively small number of radar measurements processed by FAR, which stems from its usage of a single narrowband waveform for each pulse. The performance reduction can be relieved by using compressed sensing (CS) algorithms that exploit sparsity of the target scheme [7]. However, the degradation becomes notable in extremely congested or contested electromagnetic environments [8], where there may be no vacant bands in some pulses or some radar returns of the transmitted pulses may be discarded due to strong interference [9, 10].
The performance degradation of FAR can be mitigated by using multi-carrier transmissions. When multiple carriers are transmitted simultaneously in a single pulse, the number of radar measurements is increased, and the target reconstruction performance is improved. Various multi-carrier radar schemes have been studied in the literature, including (FDMA-MIMO) [11, 12], sub-Nyquist radar (SUMMeR) [13], and frequency diversity array (FDA) radar [14, 15]. In the aforementioned schemes, different array elements transmit waveforms at different frequencies, usually forming an omnidirectional beam and illuminating a large field-of-view[16]. This degrades radar performance, especially in track mode, where a highly directional beam focusing on the target is preferred[16]. In addition, frequency agility is not exploited in FDMA-MIMO and FDA. The derivation of frequency agile multi-carrier schemes for phased array radar, which leads to a focused beam with high gain, is the focus of this work.
Here, we propose two multi-carrier agile phased array radar schemes. The first uses all the antenna elements to transmit a single waveform consisting of multiple carriers simultaneously in each pulse. Frequency agility is induced by randomly selecting the carriers utilized, resulting in a wideband multi-carrier agile radar (WMAR) scheme. While the increased number of carriers is shown to achieve improved reconstruction performance compared to conventional FAR [8], WMAR utilizes multiband signals of large instantaneous bandwidth. Therefore, its implementation does not benefit from the simplifications associated with utilizing conventional narrowband monotone signals, and may suffer from envelope fluctuation [17].
To overcome the use of instantaneously wideband waveforms, we next develop multi-Carrier AgilE phaSed Array Radar (CAESAR), which combines frequency agility and spatial agility. Specifically, CAESAR selects a small number of carrier frequencies on each pulse and randomly allocates different carrier frequencies among its antenna elements, such that each array element transmits a narrowband constant modulus waveform, facilitating system implementation. An illustration of this transmission scheme is depicted in Fig. 1.
For each carrier frequency, dedicated phase shifts on the corresponding sub-array elements are used to yield a directional transmit beam, allowing to illuminate the tracked target in a similar manner as phased array radar. Despite the fact that only a sub-array antenna is utilized for each frequency, the antenna-frequency hopping strategy of CAESAR results in array antenna gain loss and a relatively small performance gap compared to wideband radar equipped with the same antenna array. Furthermore, the combined randomization of frequency and antenna allocation can be exploited to realize a dual function radar-communications (DFRC) system [18, 19, 20, 21] by embedding digital information into the selection of these parameters. We study the application of CAESAR as a DFRC system in a companion paper [22], and focus here on the radar and its performance.
To present WMAR and CAESAR, we characterize the signal model for each approach, based on which we develop a recovery algorithm for high range resolution (HRR), Doppler, and angle estimation of radar targets. Our proposed algorithm utilizes CS methods for range-Doppler reconstruction, exploiting its underlying sparsity, and applies matched filtering to detect the angles of the targets. We provide a detailed theoretical analysis of the range-Doppler recovery performance of our proposed algorithm under complex electromagnetic environments. In particular, we prove that CAESAR and WMAR are guaranteed to recover with high probability a number of scattering points which grows proportionally to the square root of the number of different narrowband signals used, i.e., the number of carrier frequencies that are simultaneously transmitted in each pulse. This theoretical result verifies that increasing the number of carriers improves target recovery, and reduces performance degradation due to intense interference in complex electromagnetic environments. WMAR and CAESAR are evaluated in a numerical study, where it is shown that their range-Doppler reconstruction performance as well as robustness to interference are substantially improved compared to FAR. Additionally, it is demonstrated that the performance of CAESAR is only within a small gap from that achievable using wideband WMAR.
The remainder of the paper is structured as follows: Sections II and III present WMAR and CAESAR, respectively. Section IV introduces the recovery algorithm to estimate the range, Doppler, and angle of the targets. In Section V we discuss the pros and cons of each scheme compared to related radar methods. Section VI derives theoretical performance measures of the recovery method. Simulation results are presented in Section VII, and Section VIII concludes the paper.
Throughout the paper, we use , to denote the sets of complex, real numbers, respectively, and use for the magnitude or cardinality of a scalar number, or a set, respectively. Given , denotes the largest integer less than or equal to , and represents the binomial coefficient. Uppercase boldface letters denote matrices (e.g., ), and boldface lowercase letters denote vectors (e.g., ). The -th element of a matrix is denoted as , and similarly is the -th entry of the vector . Given a matrix , and a number (or a set of integers, ), () is the -th column of (the sub-matrix consisting of the columns of indexed by ). Similarly, is the sub-vector consisting of the elements of indexed by . The complex conjugate, transpose, and the complex conjugate-transpose are denoted , , , respectively. We denote as the norm, is the number of non-zero entries, and is the Frobenius norm. The probability measure is , while and are the expectation and variance of a random argument, respectively.
II WMAR
In this section we present the proposed WMAR scheme, which originates from FAR [1], aiming to increase the number of radar measurements and improve the range-Doppler recovery performance. We first briefly review FAR in Subsection II-A. Then, we detail the proposed WMAR in Subsection II-B, and present the resulting radar signal model in Subsection II-C.
II-A Preliminaries of FAR
FAR [1] is a technique for enhancing the ECCM and EMC performance of radar systems by using randomized carrier frequencies. In the following we consider a radar system equipped with antenna elements, uniformly located on an antenna array with distance between two adjacent elements. Let be the number of radar pulses transmitted in each coherent processing interval (CPI). Radar pulses are repeatedly transmitted, starting from time instance to , , where and represent the pulse repetition interval and pulse duration, respectively, and . Let be the set of available carrier frequencies, given by , where is the initial carrier frequency, , is the number of available frequencies, and is the frequency step.
In the -th radar pulse, FAR randomly selects a carrier frequency from . The waveform sent from each antenna for the -th pulse at time instance is , where
[TABLE]
and for and zero otherwise, representing rectangular envelope baseband signals.
In order to direct the antenna beam pointing towards a desired angle , the signal transmitted by each antenna is weighted by a phase shift [23], given by
[TABLE]
where denotes the speed of light. Define the vector whose -th entry is . The transmitted signal can be written as
[TABLE]
The vector in (3) denotes the transmission vector of the full array for the -th pulse at time instance .
The fact that FAR transmits monotone waveform facilitates its realization. Furthermore, the frequency agility achieved by randomizing the frequencies between pulses enhances survivability in complex electromagnetic environments. However, this comes at the cost of reduced number of radar measurements, which degrades the target recovery performance, particularly in the presence of interference, where some of the radar returns are missed [8]. To overcome these drawbacks, in the following we propose WMAR, which extends FAR to multi-carrier transmissions.
II-B WMAR Transmit Signal Model
WMAR extends FAR to multi-carrier signalling. Broadly speaking, WMAR transmits a single multiband waveform from all its antennas, maintaining frequency agility by randomizing a subset of the available frequencies on each pulse.
Specifically, in the -th radar pulse, WMAR randomly selects a set of carrier frequencies from , . We assume that the cardinality of is constant, i.e., for each , and write the elements of this set as , . The portion of the -th pulse of WMAR in the -th frequency is given by , and the overall transmitted vector is , i.e.,
[TABLE]
where the factor guarantees that (4) has the same total power as the FAR signal (3).
FAR is a special case of WMAR under the setting . By using multiple carriers simultaneously via wideband signalling, i.e., , WMAR transmits a highly directional beam, while improving the robustness to missed pulses compared to FAR. The improved performance stems from the use of multi-carrier transmission, which increases the number of radar measurements. To see this, we detail the received signal model of WMAR in the following subsection.
II-C WMAR Received Signal Model
We next model the received signal processed by WMAR for target identification. To that aim, we focus on the time interval after the -th pulse is transmitted, i.e., . In this period, the radar receives echoes of the pulse, which are sampled and processed in discrete-time.
To formulate the radar returns, we assume an ideal scattering point, representing either target or clutter, with scattering coefficient located in the transmit beam of the radar with direction angle , i.e., . Denote by the range between the target/clutter and the first radar antenna array element at time . The scattering point is moving at a constant velocity radially along with the radar line of sight, i.e., . Under the “stop and hop” assumption [24, Page 99, Ch. 2], which assumes that the target hops to a new location when the radar transmits a pulse and stays there until another pulse is emitted, the range in the -pulse is approximated as
[TABLE]
To model the received signal, we first consider the -th radar pulse that reaches the target, denoted by . Let be its component at frequency , i.e., . Note that is a summation of delayed transmissions from the corresponding antenna elements. The delay for the -th array element is . Under the narrowband, far-field assumption, using (2), we have that
[TABLE]
Substituting (5) and the definition of into (6) yields
[TABLE]
where is the relative direction sine with respect to the transmit beam, and is the transmit gain, expressed as
[TABLE]
Note that approaches when .
Having modeled the signal which reaches the target, we now derive the radar returns observed by the antenna array. After being reflected by the scattering point, the signal at the -th frequency propagates back to the -th radar array element with an extra delay of , resulting in
[TABLE]
The echoes vector can be written as
[TABLE]
where follows from (7).
The received echoes at all frequencies are then separated and sampled independently by each array element. The signal is sampled with a rate of at time instants , , such that each pulse is sampled once. Every sample time instant corresponds to a coarse range cell (CRC), . The division to CRCs indicates coarse range information of scattering points. We focus on an arbitrary -th CRC, assuming that the scattering point does not move between CRCs during a CPI, i.e., there exists some integer such that
[TABLE]
Collecting radar returns from pulses and elements at the same CRC yields a data cube with entries
[TABLE]
where is the CRC index. The data cube is processed to estimate the refined range information, Doppler, and angle of the scattering point. Data cubes from different CRCs are processed identically and separately.
Finally, we formulate how the unknown parameters of the targets are embedded in the processed data cube . To that aim, define as the high-range resolution distance, as the carrier frequency index, and as the relative frequency factor. Then, denoting by , and the generalized scattering intensity, and the normalized range and velocity, respectively, and substituting (10) into (12), we have that
[TABLE]
The unknown parameters in (13) are , , and , which are used to reveal the scattering intensity , HRR range , velocity and angle of the target.
The above model can be naturally extended to noisy multiple scatterers. When there are scattering points inside the CRC instead of a single one as assumed previously, the received signal is a summation of returns from all these points corrupted by additive noise, denoted by . Following (13), the entries of the data matrix are
[TABLE]
where , , and represent the sets of factors of scattering coefficients, ranges, velocities, and angles of the scattering points, respectively, which are unknown and should be estimated. A method for recovering these parameters from the data cube is detailed in Section IV.
WMAR has several notable advantages: First, as an extension of FAR, it preserves its frequency agility and is suitable for implementation with phased array antennas. Furthermore, as we discuss in Section V, its number of radar measurements for each CRC is increased by a factor of compared to FAR, thus yielding increased robustness to interference. However, WMAR transmitters simultaneously send multiple carriers instead of a monotone as in FAR, which requires large instantaneous bandwidth, leading to envelope fluctuation and low amplifier efficiency. To overcome these issues, we introduce CAESAR in the following section, which utilizes narrowband radar transceivers while introducing spatial agility, enabling multi-carrier transmission using monotone signals at a cost of a minimal array antenna gain loss.
III CAESAR
CAESAR, similarly to WMAR, extends FAR to multi-carrier transmission. However, unlike WMAR, CAESAR utilizes monotone signalling and reception, and is thus more suitable for implementation. We detail the transmit and receive models of CAESAR in Subsections III-A and III-B, respectively.
III-A CAESAR Transmit Signal Model
Broadly speaking, CAESAR extends FAR to multi-carrier signalling by transmitting monotone waveforms with varying frequencies from different antenna elements. The selection of the frequencies, as well as their allocation among the antenna elements, is randomized anew in each pulse, thus inducing both frequency and spatial agility.
To formulate CAESAR, we consider the same pulse radar formulation detailed in Section II. Similarly to WMAR detailed in Subsection II-B, in the -th radar pulse, CAESAR randomly selects a set of carrier frequencies from . While WMAR uses the set of selected frequencies to generate wideband waveforms, CAESAR allocates a sub-array for each frequency, such that all the antenna array elements are utilized for transmission, each at a single carrier frequency. Denote by the frequency used by the -th antenna array element, . After phase shifting the waveform to direct the beam, the -th array element transmission can be written as
[TABLE]
The vector in (15) denotes the full array transmission vector for the -th pulse at time . Here, unlike FAR which transmits a single frequency from the full array (3), CAESAR assigns diverse frequencies to different sub-array antennas, as illustrated in Fig. 1.
The transmitted signal (15) can also be expressed by grouping the array elements which use the same frequency . Let with zero padding represent the portion of which utilizes , i.e.,
[TABLE]
where is a diagonal selection matrix with diagonal , whose -th entry is one if the -th array element transmits at frequency and zero otherwise, i.e., and when . The transmitted signal is thus , namely
[TABLE]
Comparing (17) with (4), we find that each array element of CAESAR transmits a single frequency with unit amplitude while in WMAR all frequencies with amplitudes scaled by a factor are sent by each element.
The diagonal selection matrices uniquely describe the allocation of antenna elements for the -th pulse. CAESAR transmission scheme implies that , i.e., all the antenna elements are utilized for the transmission of the -th pulse. The trace of represents the number of antennas using the -th frequency. Without loss of generality, we assume that is an integer and , for each and .
Phased array FAR and FDA [15] are special cases of CAESAR with and , respectively. A fundamental difference between these radar schemes is the transmit beam pattern. In FAR, the same carrier frequency is utilized by all the elements, i.e., and are identical for each and , respectively, resulting in highly directional beam. In FDA, all available frequencies are transmitted simultaneously and one frequency corresponds to a single antenna element, leading to an omnidirectional beam which degrades radar performance and is not suitable for target tracking [16]. The proposed CAESAR uses only a subset of the available frequencies in each pulse and multiple antenna elements share the same frequency, thus achieving a compromise radiation beam that only illuminates the desired angle. Despite the gain loss in comparison with FAR discussed in Section V, CAESAR achieves improved range-Doppler reconstruction performance and increased robustness to interference, as numerically demonstrated in Section VII.
III-B CAESAR Received Signal Model
We next model the received signal processed by CAESAR. Unlike WMAR, in which each antenna receives and separates different frequency components, in CAESAR, the -th antenna element only receives radar returns at frequency , and abandons other frequencies. This enables the use of narrowband receivers, simplifying the hardware requirements.
Note that the derivation of the signal component received at the -th frequency in (6), , does not depend on the specific radar scheme. Here, substituting (16) into (6) yields
[TABLE]
where is the transmit gain of the selected sub-array antenna, expressed as
[TABLE]
Note that, in contrast to the transmit gain of WMAR in (8) which tends to , approaches when . By repeating the arguments in the derivation of (10), the echo vector can be written as
[TABLE]
CAESAR receives and processes impinging signals by the corresponding elements of the antenna array. In particular, only a sub-array, whose elements are indicated by , receives the impinging signal ; the other array elements are tuned to other frequencies. The zero-padded received signal at the -th frequency, denoted by , is thus . The full array received signal is given by .
The observed signal is sampled in a similar manner as detailed in Subsection II-C. Since CAESAR processes a single frequency component per antenna element, the measurements from each CRC are collected together as a data matrix , as opposed to a cube processed by WMAR. By repeating the arguments used for obtaining (13), it holds that
[TABLE]
which can be extended to account for multiple targets and noisy measurements as in (14), i.e.,
[TABLE]
where is the additive noise. In order to recover the unknown parameters from the acquired data matrix (22), in the following section we present a dedicated recovery scheme.
IV Target Recovery Method
Here, we present an algorithm for reconstructing the unknown HRR range, velocity, angle, and scattering intensity parameters of the scattering points from the radar measurements of both WMAR and CAESAR. Detection is performed based on the estimated scattering intensities. The detected scattering points may belong to either target or clutter, and they are identified by their Doppler estimates. The motivation for this approach is that in many ground-based radar systems, fast moving targets are of interest, while static or slow moving scatterers with zero or nearly zero Doppler are regarded as clutter. We henceforth model the Doppler of both targets and clutter as unknown parameters, which are simultaneously estimated. In specific applications where the clutter Doppler is a-priori known, one can apply clutter mitigation [25] in advance to the target recovery method. A similar procedure is also applied in pulse Doppler radars [24, Ch. 5.5.1], where the moving target indication filtering for gross clutter removal is placed prior to the pulse Doppler filter bank.
In order to maintain feasible computational complexity, we do not estimate all the parameters simultaneously: our proposed algorithm first jointly recovers the range-Doppler parameters followed by estimation of the unknown angles. When performing joint range-Doppler estimation, we assume that all the scattering points are located within the mainlobe of the transmit beam, and that the difference of the angle sine is negligible, i.e., . We then estimate the direction angles of scattering points based on their range-Doppler estimates.
We divide the target recovery method into three stages: 1) apply receive beamforming such that the magnitude of the received signal is enhanced, facilitating range-Doppler recovery; 2) apply CS methods for joint reconstruction of range and Doppler, followed by a target detection procedure; and 3) angle and scattering intensity estimation. These steps are discussed in Subsections IV-A-IV-C, respectively. A theoretical analysis of the range-Doppler estimation performance of our algorithm is provided in Section VI, where we quantify how using multiple carriers improves the range-Doppler reconstruction performance.
IV-A Receive Beamforming
The first step in processing the radar measurements is to beamform the received signal in order to facilitate recovery of the range-Doppler parameters. This receive beamforming is applied to radar returns at different frequencies separately. To formulate the beamforming technique, we henceforth focus on the -th frequency of the -th pulse, . For both CAESAR and WMAR, a total of measurements correspond to , and are denoted by . For CAESAR, is given by , of which only elements corresponding to the selected sub-array are nonzero. For WMAR, consists of the entries for each . These measurements are integrated with the weights such that the receive beam is pointed towards , resulting in
[TABLE]
Define for CAESAR, and for WMAR. When , i.e., the beam direction is close to the true angle of the target, the resulting beam pattern can be simplified as stated in the following lemma:
Lemma 1**.**
If the difference of the angle sine satisfies , then in (23) can be approximated as
[TABLE]
Proof.
See Appendix -A. ∎
The receive beamforming produces the matrix whose entries are , for each , . Under the approximation (24), the obtained is used for range-Doppler reconstruction, as discussed in the next subsection.
IV-B Range-Doppler Reconstruction Method
To reconstruct the range-Doppler parameters and detect targets in the presence of noise and/or clutter, we first recast the beamformed signal model of Lemma 1 in matrix form, and then apply CS methods to recover the unknown parameters, exploiting the underlying sparsity of the resulting model. The targets of interest are then identified based on the estimated parameters.
To obtain a sparse recovery problem, we start by discretizing the range and Doppler domains. Recall that and denote the normalized range and Doppler parameters, with resolutions and , corresponding to the numbers of available frequencies and pulses, respectively. Both parameters belong to continuous domains in the unambiguous region . We discretize and into HRR and Doppler grids, denoted by grid sets and , with grid intervals and , respectively, and assume that the targets are located precisely on the grids. The target scene can now be represented by the matrix with entries
[TABLE]
We can now use the sparse structure of (25) to formulate the range-Doppler reconstruction as a sparse recovery problem. To that aim, let and be the vectorized representations of and , respectively, i.e., and . From (24), it holds that
[TABLE]
where the entries of are given by
[TABLE]
, , and . The matrix is determined by the frequencies utilized in each pulse. Consequently, is a random matrix, as these parameters are randomized by the radar transmitters, whose realization is known to the receiver.
In the presence of noisy radar returns, (26) becomes
[TABLE]
where the entries of the noise vector are the beamformed noise, e.g., for CAESAR these are given by .
Since in each pulse only a subset of the available frequencies are transmitted, i.e., , the sensing matrix in (28) has more columns than rows, , indicating that solving (28) is naturally an under-determined problem. When is -sparse, which means that there are at most non-zeroes in and , CS algorithms can be used to solve (28), yielding the estimate .
Particularly, CS methods aim to solve under-determined problems such as (26) by seeking the sparsest solution, i.e.,
[TABLE]
The optimization in (29) is generally NP-hard. To reduce computational complexity, many alternatives including optimization and greedy approaches have been suggested to approximate (29), see [26, 27].
We take optimization as an example, under which we provide a theoretical analysis and numerically evaluate the performance in Sections VI and VII, respectively. In particular, in the absence of noise, we use the basis pursuit algorithm, which solves
[TABLE]
instead of (29). In noisy cases, recovering can be formulated as minimizing the norm under a constraint on the fidelity:
[TABLE]
Problem (31) can be solved using the Lasso method [26], which finds the solution to the regularized least squares as
[TABLE]
where and are predefined parameters.
Having obtained the estimate using CS methods, we can use it to identify which of these estimated parameters correspond to a true target of interest. Elements with significant amplitudes in are detected as dominant scattering points. Denote by the support set indexing these dominant scattering points, whose range-Doppler parameters are recovered from the corresponding indices in . For example, we may use some threshold to determine whether the amplitude is significant[24, Page 295, Ch. 6]. 5 In this case, the support set is given by . According to their Doppler estimates, these dominant scattering points are categorized into target of interests or clutter individually.
With the recovered range-Doppler values, we can estimate the angle and refine the scattering intensity, as detailed in the following subsection. The procedure for detecting which parameters correspond to targets of interest detailed above is based on . While this vector here represents the coarsely estimated scattering intensity, the procedure can be repeated for further separating targets from clutter based on the refined estimate obtained in the sequel.
IV-C Angle and Scattering Intensity Estimation
In this part, we refine the angle estimation of the scattering points, which are coarsely assumed within the transmit beam in the receive beamforming step, i.e., . While the following formulation focuses on CAESAR, the resulting algorithm is also applicable for WMAR as well as FAR.
We estimate the directions of the scatterers individually, as different points may have different direction angles. Since after receive beamforming some directional information is lost in , we recover the angles from the original data matrix (22). Using the obtained range and Doppler estimates, we first isolate echoes for each scattering point with an orthogonal projection, and apply a matched filter to estimate the direction angle of each scattering point. Finally, we use least squares to infer the scattering intensities.
IV-C1 Echo isolation using orthogonal projection
In order to accurately estimate the angle of each scattering point, it is necessary to mitigate the interference between scattering points. To that aim, we use an orthogonal projection to isolate echoes from each scatterer.
Let be the support set of , and infer the normalized range and Doppler parameters from . According to (22), given these parameters, the original data vector from the -th array element can be written as
[TABLE]
where has entries , and denotes the effective scattering intensities corresponding to all discrete range-Doppler grids, with -th entry . The intensities, containing unknown phase shifts and antenna gains due to angles , are estimated as
[TABLE]
where and we assume that \big{|}\widehat{\mathcal{S}}\big{|}<N and is invertible. The received radar echo from the -th scattering point, , , is then reconstructed by setting the -th row as
[TABLE]
IV-C2 Angle estimation using matched filter
With the isolated echoes of the -th scattering point, we use a matched filter to refine the unknown angle , which is coarsely assumed within the beam in the previous receive beamforming procedure, i.e., . Using (22), we write the isolated echo as , where denotes the noise matrix corresponding to the -th scattering point. The entries of the steering matrix are
[TABLE]
which can be computed using (19) with given and the estimates of the range-Doppler parameters. Note that , and are unknown, and is of interest. The value of the intensity recovered next is refined in the sequel to improve accuracy. Here, we apply least squares estimation, i.e.,
[TABLE]
Substituting \hat{\tilde{\beta}}_{s}=\left({\rm vec}\left({\bm{Y}}_{s}(\vartheta_{s})\right)\right)^{{\dagger}}{\rm vec}\big{(}\widehat{\bm{Y}}_{s}\big{)}=\frac{{\rm tr}\left({\bm{Y}}_{s}^{H}(\vartheta_{s}){\widehat{\bm{Y}}}_{s}\right)}{\left\|{\bm{Y}}_{s}(\vartheta_{s})\right\|_{F}^{2}} into (35) yields a matched filter
[TABLE]
The angle is estimated for each via (36) separately.
IV-C3 Scattering intensity estimation using least squares
When , there exist approximation errors in (24) and the resultant intensity estimate . We thus propose to refine the estimation of from the original data matrix once the range-Doppler and angle parameters are acquired. Given estimated angles , we concatenate the steering vectors into
[TABLE]
. The model (22) is rewritten as , and can be re-estimated via least squares as
[TABLE]
The overall parameter estimation method is summarized as Algorithm 1.
V Comparison to Related Radar Schemes
We next compare our proposed WMAR and CAESAR schemes, and discuss their relationship with relevant previously proposed radar methods.
V-A Comparison of WMAR, CAESAR, and FAR
We compare our proposed techniques to each other, as well as to FAR, which is a special case of both WMAR and CAESAR obtained by setting . We focus on the following aspects: 1) instantaneous bandwidth; 2) the number of measurements in a CPI; and 3) signal-to-noise ratio (SNR). A numerical comparison of the target recovery performance of the considered radar schemes is provided in Section VII.
In terms of instantaneous bandwidth, recall that CAESAR and FAR use narrowband transceivers, and only a monotone signal is transmitted or received by each element. In WMAR, multi-tone signals are sent and received simultaneously in each pulse, thus it requires instantaneous wideband components.
To compare the number of obtained measurements, we note that for each CRC, WMAR acquires a data cube with samples, while FAR and CAESAR collect samples in the data matrix. After receive beamforming, the number of observations become , and , for FAR, CAESAR, and WMAR, respectively, via (23). This indicates that the multi-carrier waveforms of CAESAR and WMAR increase the number of measurements after receive beamforming.
The aforementioned radar schemes also differ in their SNR, as the transmitted power and antenna gains differ. Here, as in [24, Page 304, Ch. 6], SNR refers to the ratio of the power of the signal component to the power of the noise component after coherently accumulating the radar returns. To see this difference, we consider the case when there exists a target with range-Doppler-angle , scattering coefficient and the noise elements in are i.i.d. zero-mean proper-complex Gaussian with variance . In this case, coherent accumulation of radar returns reduces to in WMAR and in CAESAR (FAR can be regarded as a special case of WMAR/CAESAR with ), and the antenna gains are and , respectively. It follows from (14) that the signal amplitude in radar returns of WMAR is . After coherent accumulation, the amplitude becomes , leading to a signal power of , while the power of the noise component becomes . Hence, the SNR of WMAR is . In CAESAR, the amplitude of the signal component in is , which becomes after accumulation. Since there are only noise elements in CAESAR, the noise power after accumulation is and the resultant SNR is . Letting implies that the SNR of FAR is also . The SNR calculation indicates that CAESAR has an SNR loss by a factor of compared to WMAR and FAR. This loss stems from the fact that CAESAR uses a subset of the antenna array for each carrier, and thus has lower antenna gain than FAR and WMAR. This SNR reduction can affect the performance of CAESAR in the presence of noise, as numerically demonstrated in Section VII.
The above comparison reveals the tradeoff between instantaneous bandwidth requirement, number of observations, and post-accumulation SNR. Among these three factors, the number of observations is crucial to the target recovery performance especially in complex electromagnetic environments, where some observations may be discarded due to strong interference [9, 10]. The proposed multi-carrier schemes, WMAR and CAESAR, are numerically shown to outperform FAR in Section VII, despite the gain loss of CAESAR. CAESAR also achieves performance within a relatively small gap compared to WMAR, while avoiding the usage of instantaneous wideband components. The resulting tradeoff between number of beamformed observations and SNR, induced by the selection of , is not the only aspect which must be accounted for when setting the value of , as it also affects the frequency agility profile. In particular, smaller values result in increased spectral flexibility, as different pulses are more likely to use non-overlapping frequency sets. Consequently, in our numerical analysis in Section VII we use small values of , for which the gain loss between CAESAR and WMAR is less significant, and increased frequency agility is maintained. In addition, CAESAR can also exploit its spatial agility character, which is not present in FAR or WMAR, to realize a DFRC system, as discussed in our companion paper [22].
V-B Comparison to Previously Proposed Schemes
Similarly to CAESAR, previously proposed FDMA-MIMO radar [11], SUMMeR [13], and FDA radar [14, 15] schemes also transmit a monotone waveform from each antenna element while different elements simultaneously transmit multiple carrier frequencies. The main differences between our approaches and these previous methods are beam pattern and frequency agility. Due to the transmission of diverse carrier frequencies from different array elements of FDMA-MIMO/SUMMeR/FDA, the array antenna does not form a focused transmit beam and usually illuminates a large field-of-view [16]. This results in a transmit gain loss which degrades the performance, especially for track mode, where a high-gain directional beam is preferred [16]. By transmitting each selected frequency with an antenna array (the full array in WMAR and a sub-array in CAESAR), our methods achieve a focused beam pattern that facilitates accurate target recovery.
Furthermore, FDMA-MIMO and FDA transmit all available frequencies simultaneously, and thus do not share the advantages of frequency agility, e.g., improved ECCM and EMC performance, as the multi-carrier version of SUMMeR and the proposed WMAR/CAESAR. In addition, FDMA-MIMO, SUMMeR and WMAR receive instantaneous wideband signals with every single antenna, as opposed to FDA [15] and CAESAR, which use narrowband receivers.
To summarize, we compare in Table I the main characteristics of these radar schemes. Unlike previously proposed radar methods, our proposed techniques are based on phased array antenna and frequency agile waveforms to achieve directional transmit beam and high resistance against interference. In terms of instantaneous bandwidth, CAESAR is preferred for its usage of monotone waveforms and simple instantaneous narrowband receiver.
VI Performance Analysis of Range-Doppler Reconstruction
Range-Doppler reconstruction plays a crucial role in target recovery. This section presents a theoretical analysis of range-Doppler recovery using CS methods. Since both WMAR and CAESAR are generalizations of FAR, the following analysis is inspired by the study of CS-based FAR recovery in [7]. In particular, we extend the results of [7] to multi-carrier waveforms, as well as to extremely complex electromagnetic environments, where some transmitted pulses are interfered by intentional or unintentional interference. In the presence of such interference, only partial observations in the beamformed matrix remain effective for range-Doppler reconstruction. To present the analysis, we first briefly review some basic analysis techniques of CS in Subsection VI-A, followed by the range-Doppler recovery performance analysis in Subsection VI-B.
VI-A Preliminaries
There have been extensive studies on theoretical conditions that guarantee unique recovery for noiseless models or robust recovery for noisy models [26, 27]. The majority of these studies characterize conditions and properties of the measurement matrix , including spark, mutual incoherence property (MIP) and restricted isometry property (RIP).
Following [7], we focus on the MIP. A sensing matrix is said to satisfy the MIP when its coherence, defined as
[TABLE]
is not larger than some predefined threshold. Bounded coherence ensures unique or robust recovery using a variety of computationally efficient CS methods. We take optimization as an example to explain the bounds on matrix coherence. In the absence of noise, the uniqueness of the solution to (30) is guaranteed by the following theorem:
Theorem 2** ([28]).**
Suppose the sensing matrix has coherence . If solves (30) and has support size at most , then is the unique solution to (30).
In noisy cases, the following result shows that also guarantees stable recovery.
Theorem 3** ([29]).**
Consider the model (28) with and , and let sensing matrix have coherence . If solves (31), then
[TABLE]
Based on Theorems 2 and 3, we next analyze the coherence measure of the sensing matrix in (27) for CAESAR and WMAR (whose sensing matrices are identical), and establish the corresponding performance guarantees.
VI-B Performance Analysis
Here, we analyze the range-Doppler reconstruction of WMAR and CAESAR. Since the sensing matrix is random, we start by analyzing its statistics, and then derive conditions that ensure unique recovery by invoking Theorem 2.
We assume that the frequency set is uniformly i.i.d. over . For mathematical convenience, in our analysis we adopt the narrow relative bandwidth assumption from [7], i.e., , such that (27) becomes
[TABLE]
Numerical results in [7] indicate that large relative bandwidth has negligible effect on the MIP of . In addition, recall that all the targets precisely lie on the predefined grid points, as assumed in Subsection IV-B. Here, we adopt the on-the-grid assumption for mathematical convenience. Consequently, the accuracy and actual resolutions of range and Doppler reconstruction results are restricted by the grid intervals, i.e., and , respectively. In practical scenarios, we may use denser grid points, as will be discussed by simulations in Subsection VII-B. Denser grid enhances the obtainable accuracy/resolution and alleviates the performance loss when the real parameters are off the grid points. However, the density of grid points cannot go to infinity, because denser grid affects the incoherence property of the observation matrix while increasing the memory requirements and the computational burden. An alternative approach to overcome the need to specify a range-Doppler grid is to utilize off-the-grid CS methods, see [30, 31].
In complex electromagnetic environments, some of the radar echoes may be corrupted due to jamming or interference. Heavily corrupted echoes are unwanted and should be removed before processing in order to avoid their influence on the estimation of target parameters [9, 10]. In this case, the corrupted radar returns are identified, as such echoes typically have distinct characters, e.g., extremely large amplitudes. These interfered observations are regarded as missing, where we consider two kinds of missing patterns: 1) pulse selective, i.e., all observations in certain pulses are missing, which happens when the interference in these pulses is intense over all sub-bands; 2) observation selective, namely, only parts of the observations are missed when the corresponding pulse is interfered. We consider the first case, assuming that the radar receiver knows which pulses are corrupted, and leave the analysis under the second case for future investigation. In particular, we adopt the missing-or-not approach [32], in which each pulse in has a probability of , , to be corrupted, and the missing-or-not status of the pulses are statistically independent of each other.
After removing the corrupted returns, only part of the observations in the beamformed vector (28) are used for range-Doppler recovery. Equivalently, corresponding rows in can be regarded as missing, affecting the coherence of the matrix and thus the reconstruction performance. Denote by the random set of available pulse indexes and by the corresponding index set of available observations. The signal model (28) is now
[TABLE]
where , , and .
Consider the inner product of two columns in , denoted and , corresponding to grid points and , respectively, , , . While there are different pairs of , the magnitude of the inner product \big{|}\left[{\bm{\Phi}}_{\ast}\right]_{l_{1}}^{H}\left[{\bm{\Phi}}_{\ast}\right]_{l_{2}}\big{|}, which determines the coherence of , takes at most distinct random values. To see this, note that
[TABLE]
indicating that the inner product depends only on the difference of the grid points, i.e., and , and not on the individual values of the column indices and . It follows from (42) that the MIP of can be written as
[TABLE]
where , take values in the sets and , respectively.
Next, we define
[TABLE]
where the random variable satisfies when and 0 otherwise. In addition, let
[TABLE]
Some of the magnitudes are duplicated since and . To eliminate the duplication and remove the trivial nonrandom value , we restrict the values of and to and , respectively, and define the set
[TABLE]
with cardinality , such that each value of \big{|}\left[{\bm{\Phi}}_{\ast}\right]_{l_{1}}^{H}\left[{\bm{\Phi}}_{\ast}\right]_{l_{2}}\big{|} (except the trivial case ) corresponds to a single element of the set . We can now write (43) as
[TABLE]
The coherence in (47) is a function of the dependent random variables and . To bound , we derive bounds on and , respectively. To this aim, we first characterize the statistical moments of for some fixed , which we denote henceforth as , in the following lemma:
Lemma 4**.**
The sequence of random variables satisfies
[TABLE]
[TABLE]
Furthermore, for each ,
[TABLE]
Proof.
See Appendix -B. ∎
Using Lemma 4, the probability that the magnitude is bounded can be derived as in the following Corollary:
Corollary 5**.**
Let . For any and it holds that
[TABLE]
Proof.
Based on the definition (44) and the independence assumption on the frequency selection and missing-or-not status of each pulse, it holds that are independent zero-mean complex-valued random variables. Now, since
[TABLE]
according to (48) and equals 0 for , recalling that , we have . Then, it holds that . Combining Bernstein’s inequality [33, Thm. 12] with the fact that by (50), , results in (51). ∎
We next derive a bound on the number of effective pulses in the following lemma:
Lemma 6**.**
For any , it holds that
[TABLE]
Proof.
Since, by its definition, obeys a binomial distribution, (53) is a direct consequence of [34, Thm. 1]. ∎
Based on the requirement in Theorem 2, we now use Corollary 5 and Lemma 6 to derive a sufficient condition on the radar parameters , , , as well as the intensity of interference , guaranteeing that the measurement matrix meets the requirement with high probability. This condition is stated in the following theorem:
Theorem 7**.**
For any constant , the coherence of satisfies when
[TABLE]
Proof.
See Appendix -C. ∎
Recall that the value of depends on the quantities , and . When is reasonably large such that , we have . When there is no noise in radar returns, a number of scattering points (on the grid) in the scale of guarantees a unique reconstruction of range-Doppler parameters with high probability according to Theorems 2 and 7. Note that this rather simple asymptotic condition assumes that , i.e., that the overall number of available frequencies is substantially larger than the number of frequencies utilized in each pulse , thus ensuring the agile character in frequency domain. Compared to the asymptotic condition of FAR with full observations [7], we find that the presence of corrupted observations, i.e. when , leads to degraded range-Doppler reconstruction performance. However, by increasing the number of transmitted frequencies in each pulse while maintaining , the performance deterioration due to missing observations can be mitigated, enhancing the interference immunity of the radar in extreme electromagnetic environments. In the special case that and , i.e., FAR in an interference free environment, the two conditions coincide as .
The above condition is proposed for the noiseless case, indicating that the inherent target/clutter reconstruction capacity increases with . In practical noisy cases, the reconstruction performance does not monotonically increase with , because the transmit power of each frequency decreases with , thus degrading the SNR in both CAESAR and WMAR. Particularly, in CAESAR, larger means that less antennas () are allocated to each frequency, which affects the radiation beam and enlarges the gain loss. In addition, a small maintains the practical advantages of frequency agility in terms of, e.g., ECCM and EMC performance.
VII Simulation Results
In this section, we numerically compare the performance of WMAR, CAESAR, and FAR in noiseless/noisy, clutter, and/or jamming environments. The performance is evaluated in terms of target detection probability, accuracy and resolution, probability of correct reconstruction, and mutual interference, as presented in Subsections VII-A - VII-D, respectively.
We consider a frequency band starting from GHz, with 4 available carriers and carrier spacing of MHz. The radar system is equipped with an antenna array of elements with spacing of , and utilizes pulses focusing on . CAESAR and WMAR use frequencies at each pulse. In noisy scenarios, we use the term SNR for the post-accumulation SNR of WMAR/FAR, i.e., as derived in Subsection V-A. To guarantee fair comparison, we use the same definition for all the radar schemes. In the presence of jamming, we test the pulse selective missing pattern. In order to implement target recovery via Algorithm 1 for the three radar schemes, we use the convex optimization toolbox [35] to implement basis pursuit (30) in noiseless cases or the Lasso algorithm (32) with in noisy setups for range-Doppler reconstruction.
VII-A Target Detection in Clutter Environment
We compare the detection performance of the proposed WMAR and CAESAR schemes with FAR and conventional fixed frequency radars, which can be regarded as a special case of WMAR, CAESAR and FAR, with and . To that aim, we evaluate the detection probabilities, , of a moving target under ground clutter environment. The ground clutter is modeled as radar returns from many static scattering points with velocity being zero. Thus, moving targets and ground clutter are distinguishable by observing their Doppler values. We denote by the index set corresponding to zero Doppler, i.e., , and by its complementary set, i.e., . The number of target scattering points is . The normalized range and Doppler parameter of the target is determined by randomly selecting the index from . The target intensity is selected to match the desired SNR. The number of clutter scattering points is set as , and the intensity of each scattering point is unity, with a random phase uniformly distributed over . The normalized range parameters of these scatterers are uniformly distributed over , and are not assumed to lie on the predefined grid points. Since these scatterers have identical velocity, the superposition of their echoes, i.e., the clutter signal, can be represented by radar returns from the grid points indexed by . Echoes from both target and clutter are treated as unknown, and are reconstructed simultaneously by CS. Particularly, we apply Lasso (32) for range-Doppler reconstruction, yielding , where elements indexed in are regarded as equivalent clutter intensities and the remaining elements are regarded as intensities of moving targets. The curves are plotted versus SNR, and we choose .
The simulations are carried out in two parts. In the first part, there are only clutter and additive noises in the radar returns, without returns from moving targets. The resulting radar measurements are used to determine the detection thresholds for all radar schemes, respectively, under a given probability of false alarm, denoted . In the second part, these thresholds are used to detect the existence of a target from the received echoes, which now include noise as well as returns from both clutter and moving target. The target detection procedure is based on the estimated parameters , as detailed in Subsection IV-B. In the first part, we set and perform Monte Carlo trials. While radar systems typically operate at lower values of , we use this value for computational reasons. It can be faithfully simulated under the given number of trials, and the selected value of provides a characterization and understanding of the behavior of the considered radar schemes. A false alarm is proclaimed if any nonzero-Doppler element in exceeds certain threshold , i.e., . In the second part, we execute Monte Carlo trials. A successful detection is proclaimed if the estimated intensity of the moving target is larger than the threshold, , . To evaluate the influence of missing observation caused by jamming, we set the survival rate and compare the resulting probability of detection, curves with those of the full observation cases in Fig. 2. Note that the detection thresholds for these curves are calculated individually.
As shown in Fig. 2, the fixed frequency radar has higher detection probabilities than the counterparts of frequency agile schemes. The advantage of fixed frequency radar stems from the property of its observation matrix , where and becomes an orthogonal matrix, benefiting the Doppler reconstruction performance of CS methods. While in the frequency agile schemes, generally it holds that , resulting in an incomplete observation matrix and degradation of clutter/target recovery performance. However, we note that the fixed frequency scheme is vulnerable in a jamming environment. In the full observation cases, WMAR outperforms FAR because of the increased number of observations. Though CAESAR has identical number of observations with WMAR after receive beamforming, the values of CAESAR are less than those achieved by WMAR with a SNR gap of approximately 6 dB. This follows since CAESAR suffers from an antenna gain loss of as discussed in Subsection V-A (i.e., 6 dB since ). When some of the observations are missing due to jamming, the detection probabilities are affected. Both WMAR and CAESAR suffer from an SNR loss of approximately 3 dB, while FAR fails to detect any moving target in the scenarios under test. When FAR is lacking in radar observations, the mutual coherence property of its observation matrix becomes degraded, leading to many spurious peaks of high intensities in the recovery results . These spurious peaks significantly increase the detection threshold, thus operating under a fixed of results in notably reduced detection probability .
To summarize, we find from the simulation results that 1) the proposed multi-tones schemes (WMAR and CAESAR) enhance the immunity against missing data over the single-tone FAR, and 2) WMAR outperforms CAESAR due to its higher antenna gain, which comes at the cost of increased instantaneous bandwidth.
VII-B Accuracy and Resolution
Here, we compare the range, Doppler and angle estimate results under different range-Doppler grid points. To this aim, two sets of range-Doppler grid points are tested: one uses the standard grid as mentioned in Subsection IV-B, where the intervals of consecutive range and Doppler grid points are and , respectively; The latter uses a denser grid, setting and , and the consequent simulation results are denoted with label “-2”, e.g., CAESAR-2. The number of scattering points is without clutter. The normalized range-Doppler parameter, , of the scattering point is uniformly, randomly set over , and the angle is randomly set within the beam . Under this setting, the ground truth of the range-Doppler parameter may be off the grid, which leads to inevitable estimation error. The CS method applied for range-Doppler reconstruction is based on (32), and we estimate range-Doppler, denoted by , from the index of the element with maximum magnitude in . We then use root mean squared error (RMSE) as the metric of accuracy, defined by , taking normalized range as an example. The remaining settings are the same as those used in Subsection VII-A. We run 500 Monte Carlo trials and the range, Doppler and angle accuracy results are shown in Figs. 3, 4 and 5, respectively.
As expected, the RMSEs become lower when we increase SNR, while we observe in Figs. 3 and 4 that the RMSEs of range and Doppler estimates reach error floors as the SNR increases. The error floors depend on the grid intervals or , and denser grid points lead to lower error floors. The results also reveal that WMAR and FAR have similar accuracy performance, while CAESAR has an SNR loss of 6 dB in moderate SNR levels because of its lower antenna gain. In high SNR scenarios, the RMSEs of CAESAR also reach the error floor.
We next examine the ability of CAESAR, FAR and WMAR in separating closely spaced scattering points, i.e., obtainable resolution. In the simulations, we use dense grid points with intervals and . We consider two closely spaced scattering points, of which angles are set and range-Doppler parameters are on the grid. Particularly, we fix the range-Doppler parameter of one scattering point and change the counterpart of the other scattering from , such that the range/Doppler separation (/) between two points are changed. In this experiment, we disregard noise and the scattering intensities are both with random phase. We use (30) for range-Doppler recovery, and the two most dominant elements in the estimate are regarded as scattering points. The indices of these two elements are compared with the corresponding ground truth, and a successful recovery (also referred to a hit) is proclaimed if both indices are correct. The number of Monte Carlo trails are 500. The achievable hit rate results versus separation between scattering points are shown in Fig. 6. The results demonstrate that all the frequency agile schemes, CAESAR, FAR and WMAR have close performance in resolution, while the hit rates of CAESAR and WMAR are slightly higher than those of FAR, because they have more observations and the measurement matrix has better coherence property.
VII-C Reconstruction of Multiple Scattering Points
In this subsection, we evaluate the proposed radar schemes in recovering a set of scattering points in noiseless and noisy setups. In both scenarios, we use the standard grid points with grid intervals . The range-Doppler parameters of scattering points are randomly selected from the grid points, angle parameters are randomly set from the continuous set , and scattering intensities are all set to unity. We apply CS methods for range-Doppler recovery, and the indices of most significant entries in are regarded as elements of the estimated support set. Hit rates are applied as performance metric, and a hit is proclaimed if the obtained support set is identical to the ground truth, which means all the range-Doppler parameters are reconstructed correctly.
In the noiseless experiment, we simulate different numbers of recoverable scattering points, . We set the survival rate for jamming environments. The resulting hit rates versus are depicted in Fig. 7. As expected, the hit rates decrease as increases. The performance of CAESAR is within a very small gap of that achievable using WMAR, because CAESAR and WMAR use the same amount of transmitted frequencies , and the number of beamformed measurements is also the same. Hit rates of CAESAR and WMAR exceed that of FAR significantly. This gain stems from the fact that transmitting multi-carriers in each pulse of CAESAR and WMAR increases the number of observations, and thus raises the number of recoverable scattering points.
We then consider the noisy case, and compare the range-Doppler recovery performance versus SNR, which is changed by varying . We set , and we let for the jamming environment. The hit rates of the range-Doppler parameters are depicted in Fig. 8.
Observing Fig. 8, we note that, as expected, WMAR achieves the best performance in range-Doppler reconstruction. While WMAR and CAESAR have the same number of observations, CAESAR has a lower antenna gain as noted in Subsection V-A, which results in its degraded performance compared to WMAR. In the full observation case with high SNRs, i.e., SNR dB, CAESAR has higher hit rates than FAR due to the advantage of increased number of transmitted frequencies, while in low SNRs of less than 20 dB, FAR exceeds CAESAR owing to its higher antenna gain. In the jamming scenario, CAESAR outperforms FAR, and that FAR almost fails to reconstruct scattering points (with hit rates around 0.25). The superiority of WMAR/CAESAR over FAR demonstrates the advantage of the proposed multi-carrier waveforms.
From the experimental results in Subsections VII-A - VII-C, we find that the multi-carrier signals used by CAESAR and WMAR significantly enhance range-Doppler reconstruction performance over the monotone waveform in FAR. The advantage becomes more distinct in jamming environments, where some radar measurements are invalid. In reasonably high SNR scenarios, the reconstruction performance of CAESAR, which uses narrowband constant modulus waveforms for each antenna element, approach those of WMAR, which uses instantaneously wideband waveforms.
VII-D Mutual Interference
One of the main advantages of frequency agile transmission is its relatively low level of mutual interference, which implies that multiple transmitters can coexist in dense environments. To demonstrate this property of the proposed radar schemes, which all utilize some level of frequency agility, we next evaluate the unintended mutual interference of closely placed radars transmitting the same waveform pattern. We compare the frequency agile schemes of FAR, WMAR and CAESAR, with an instantaneous wideband radar which transmits all subbands simultaneously. In the simulation, we consider a scenario with radars operating independently. Mutual interference occurs if a reference radar is receiving echoes while another radar is transmitting at the same subcarriers with their antenna beams directed towards each other. In this case the echoes of the reference radar at the conflicted subcarriers are corrupted. The level of mutual interference is measured by the average number of uncorrupted subcarriers, denoted .
We use to represent the probability that one radar may interfere the reference radar, i.e., that it is radiating during the reception period of the reference radar and their beams are directed towards each other. The number of subcarriers transmitted in each pulse varies from to , where represents the FAR while represents the radar using full bandwidth. As WMAR and CAESAR transmit the same number of subcarriers for a specific , their performance on mutual interference is the same. To demonstrate the mutual interference intensity versus interference probabilities, we simulate Monte Carlo trials for each interference probability and calculate the average number of uncorrupted subcarriers as shown in Fig. 9. From the results, we observe that, as expected, when the interference probability is small, e.g., less than , radar systems transmitting more subcarriers are capable of effectively utilizing their bandwidth reliably. However, as the probability of interference grows, wideband radar induce severe mutual interference, resulting in a negligible average number of uncorrupted subcarriers for . The frequency agile schemes, such as FAR and WMAR/CAESAR operating with , are still capable of reliably utilizing a notable portion of their bandwidth in the presence of such high interference. These results indicate that the less frequency agile the scheme is, the severer the mutual interference becomes in the high interference probability regime.
VIII Conclusion
In this work we developed two multi-carrier frequency agile schemes for phase array radars: WMAR, which uses wideband waveforms; and CAESAR, which transmits monontone signals and introduces spatial agility. We modeled the received radar signal, and proposed an algorithm for target recovery. We then characterized theoretical recovery guarantees. Our numerical results demonstrate that our proposed schemes achieve enhanced survivability in extreme electromagnetic environments. Furthermore, it is shown that CAESAR is capable of achieving performance which approaches that of wideband radar, while utilizing narrowband transceivers. An additional benefit which follows from the introduction of frequency and spatial agility is the natural implementation of CAESAR as a DFRC system, studied in a companion paper.
-A Proof of Lemma 1
In the following we prove (24) for CAESAR. The proof for WMAR follows similar arguments and is omitted for brevity.
Substituting the definitions of and into (23) yields
[TABLE]
Recall that when it holds that . Then, (-A.1) reduces to (24), proving the lemma. ∎
-B Proof of Lemma 4
We first prove (48) and (49), after which we address (50).
-B1 Proof of (48) and (49)
For brevity, let , , , and . We set when or , and when .
We first compute {\rm E}\left[\chi_{n}\right]=\frac{1}{K}{\rm E}\big{[}I_{n}e^{jqn}\sum_{k=0}^{K-1}e^{jpc_{n,k}}\big{]}. The expectation is taken over the indicator and frequency codes . Since they are independent and , it holds that
[TABLE]
Since frequencies are selected uniformly (but not independently), it follows that
[TABLE]
where denotes the -th frequency in the -th combination. Out of these combinations, there are that contain a given selection . Thus, we have that
[TABLE]
Substituting (-B.3) into (-B.2) yields
[TABLE]
As , it holds that {\rm E}\big{[}\sum_{k=0}^{K-1}e^{jpc_{n,k}}\big{]}=K if and zero otherwise. Substituting this into (-B.1), we have
[TABLE]
which proves (48).
To obtain \!{\rm D}\left[\chi_{n}\right]\!:=\!{\rm E}\big{[}\!\left|\chi_{n}\!-\!{\rm E}\!\left[\chi_{n}\right]\right|^{2}\big{]}, we consider two cases, and . When , we have and
[TABLE]
where holds since and in we apply .
When , the random variable has zero mean, and its variance is given by
[TABLE]
where we use . To compute (-B.7), we note that
[TABLE]
where follows since in the second term can be replaced by . From the derivation of (-B.5), it holds that for the second summand in (-B.8) vanishes, resulting in
[TABLE]
Plugging (-B.9) into (-B.7), we obtain
[TABLE]
Finally, to prove (49), we calculate for and . When , from (-B.6), we have that
[TABLE]
When , it follows from (-B.10) that
[TABLE]
Combining (-B.11) and (-B.12) proves (49). ∎
-B2 Proof of (50)
We again consider the two cases and separately: When , it follows from (-B.5) that
[TABLE]
When , , which is not larger than 1 by definition of (44), thus proving (50). ∎
-C Proof of Theorem 7
By fixing some positive , setting and
[TABLE]
we have that for any
[TABLE]
Here holds since the event implies that at least one of the conditions and is satisfied; and follows from Corollary 5 and Lemma 6.
Using the bound (-C.2) on the magnitude of the normalized correlation, we next bound the probability of to exceed some constant. By applying the union bound to (47), we have
[TABLE]
According to (-C.3), for any , it holds that
[TABLE]
where is obtained from via (-C.1). The right hand side of (-C.4) is not smaller than when , implying that when satisfies
[TABLE]
Finally, by (-C.1), fixing implies that
[TABLE]
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1] S. R. J. Axelsson, “Analysis of random step frequency radar and comparison with experiments,” IEEE Trans. Geosci. Remote Sens. , vol. 45, no. 4, pp. 890–904, April 2007.
- 2[2] J. Yang, J. Thompson, X. Huang, T. Jin, and Z. Zhou, “Random-frequency SAR imaging based on compressed sensing,” IEEE Trans. Geosci. Remote Sens. , vol. 51, no. 2, pp. 983–994, 2013.
- 3[3] Y. Liu, G. Liao, J. Xu, Z. Yang, and Y. Zhang, “Adaptive OFDM integrated radar and communications waveform design based on information theory,” IEEE Commun. Let. , vol. 21, no. 10, pp. 2174–2177, Oct 2017.
- 4[4] D. Cohen, K. V. Mishra, and Y. C. Eldar, “Spectrum sharing radar: Coexistence via Xampling,” IEEE Trans. Aerosp. Electron. Syst. , vol. 54, no. 3, pp. 1279–1296, 2018.
- 5[5] D. Cohen and Y. C. Eldar, “Sub-nyquist radar systems: Temporal, spectral, and spatial compression,” IEEE Signal Processing Magazine , vol. 35, no. 6, pp. 35–58, Nov 2018.
- 6[6] T. Huang, Y. Liu, G. Li, and X. Wang, “Randomized stepped frequency ISAR imaging,” in Proc. IEEE Radar Conf , May 2012, pp. 0553–0557.
- 7[7] T. Huang, Y. Liu, X. Xu, Y. C. Eldar, and X. Wang, “Analysis of Frequency Agile Radar via Compressed Sensing,” IEEE Trans. Signal Process. , vol. 66, no. 23, pp. 6228–6240, Dec 2018.
- 8[8] T. Huang, Y. Liu, Y. C. Eldar, and X. Wang, “Multiple carrier agile radar via compressed sensing,” in International Workshop on Compressed Sensing applied to Radar, Multimodal Sensing, and Imaging (Co Se Ra) 2018 , 2018.
