Accurate detection of moving targets via random sensor arrays and Kerdock codes
Abstract.
The detection and parameter estimation of moving targets is one of the most important tasks in radar. Arrays of randomly distributed antennas have been popular for this purpose for about half a century. Yet, surprisingly little rigorous mathematical theory exists for random arrays that addresses fundamental question such as how many targets can be recovered, at what resolution, at which noise level, and with which algorithm. In a different line of research in radar, mathematicians and engineers have invested significant effort into the design of radar transmission waveforms which satisfy various desirable properties. In this paper we bring these two seemingly unrelated areas together. Using tools from compressive sensing we derive a theoretical framework for the recovery of targets in the azimuthrangeDoppler domain via random antennas arrays. In one manifestation of our theory we use Kerdock codes as transmission waveforms and exploit some of their peculiar properties in our analysis. Our paper provides two main contributions: (i) We derive the first rigorous mathematical theory for the detection of moving targets using random sensor arrays. (ii) The transmitted waveforms satisfy a variety of properties that are very desirable and important from a practical viewpoint. Thus our approach does not just lead to useful theoretical insights, but is also of practical importance. Various extensions of our results are derived and numerical simulations confirming our theory are presented.
Key words and phrases:
Sparsity, Radar, Compressive Sensing, Random Sensor Arrays, MIMO, Kerdock Codes1. introduction
The detection and parameter estimation of moving targets is one of the most important radar applications. The use of antenna arrays greatly improves our ability to perform this task. Antenna arrays make it possible to estimate not only the range and Doppler frequency, but also the azimuth of the target. Furthermore, using multiple antennas can significantly increase signal strength and thus in turn can greatly enhance accuracy and our ability to locate low contrast targets (“faint” or “weak” targets).
Therefore it does not come as a surprise that in recent years radar systems employing multiple antennas at the transmitter and the receiver (also referred to as MIMO radar, where MIMO stands for multipleinput multipleoutput) have attracted enormous attention in the engineering and signal processing community. Despite the significant resources that have been devoted to MIMO radar, there exists fairly little rigorous mathematical theory for MIMO radar that addresses fundamental questions, such as how many targets can be detected at which azimuthrangeDoppler resolution and at what signaltonoise ratio. Existing theory focuses mainly on the detection of a single target [12, 26]. Only very recently, in the footsteps of compressive sensing, do we see the emergence of a rigorous mathematical theory for MIMO radar that addresses the more realistic and more interesting case of multiple targets [33]. However, for the widely popular case of randomly spaced antennas^{1}^{1}1In this paper we only consider the case of colocated transmitters and receivers, which is the most relevant situation in practice. We do not discuss the case of widely separated antennas [15]., the mathematical theory is still in its infancy.
In an independent and seemingly disparate line of research in radar, mathematicians and engineers have devoted substantial efforts to the design of radar transmission waveforms that satisfy a variety of desirable properties. The vast majority of this research has focused on single antenna radar systems, and it is a priori not clear if and how these waveforms can be utilized for MIMO radar. In this paper we bring together these two independent areas of research, MIMO radar with random antenna arrays and radar waveform design, by developing a rigorous mathematical framework for accurate target detection via random arrays, which at the same time utilizes some of the most attractive radar waveforms, such as Kerdock codes.
A radar system illuminates a region of interest in order to detect the location, velocity, and reflectivity of the objects (targets) in its field of view. We consider the following standard (narrowband) radar model [32]. Suppose a target located at range is traveling with constant velocity and has reflection coefficient . Suppose further just for the moment that we have only one target, one transmitter and one receiver (in which case we cannot detect direction). After transmitting signal , the receiver observes the reflected signal
(1) 
where is the round trip time of flight, is the speed of light, is the Doppler shift, and is the carrier frequency. The basic idea is that the rangevelocity information of the target can be inferred from the observed time delayDoppler shift of in (1). For only one target this can be done conveniently by correlating the received signal with timefrequency shifted versions of the transmitted signal. Since we are dealing with bandlimited signals, it suffices to consider discrete signals sampled at a properly chosen rate . It is therefore common practice to compute
(2) 
and then locate the largest value of in order to detect the target in the rangeDoppler domain.
In the presence of multiple targets more sophisticated methods are necessary. In order to resolve azimuth in addition to range and Doppler, we need to employ an array of antennas. We assume an array of transmit and receiver antennas that are colocated (also known as monostatic radar) as illustrated in Figure 1. A more detailed description of the setup is postponed to Section 2.
The transmit antennas send simultaneously probing signals, which can differ from antenna to antenna and can be chosen to our specifications. It is convenient to divide the region of interest into rangeazimuthDoppler cells corresponding to distance, direction and velocity, respectively. Let be a measurement matrix whose columns correspond to the signal recorded at each receive antenna from a single unitstrength scatterer at a specific rangeazimuthDoppler cell. Let denote a vector whose elements represent the complex amplitudes of the scatterers. In many cases the radar scene is sparse in the sense that only a small fraction (often a very small fraction) of the cells is occupied by the objects of interest. In this case most of the entries of will be zero, but we do not know which ones, otherwise we would have located the targets already. With representing a noise vector, we are faced with the linear system of equations
(3) 
where is a vector of measurements collected by the receive antennas over an observation interval. Typically this system will be underdetermined, which implies that it will have infinitely many solutions. What comes to our rescue here is the sparsity of . While conventional radar processing techniques do not take full advantage of sparsity of the radar scene, the recent development of compressive sensing provides us with the possibility to optimally utilize this property [17, 31, 33]. The approach pursued in this paper to obtain a sparse solution of (3) is based on the lasso [36], which gained tremendous popularity in connection with compressive sensing. The lasso solves
(4) 
where the parameter trades off goodness of fit with sparsity.
However, one of the main challenges in bringing compressive sensing theory into radar is that in radar the sensing matrix cannot be freely chosen. Its structure is dictated by the laws of physics on which radar is based. The crux is to carefully balance the desired resolution in the azimuthrangeDoppler domain with the degrees of freedom at our disposal in the formation of , such as the antenna locations and the transmit waveforms.
A great deal of work has been devoted in the mathematical and engineering literature to the design of radar transmission waveforms, see for instance [1, 3, 9, 13, 14, 29, 38] for a small sample of references. The design criteria for radar waveforms can be roughly split into two categories: (i) properties that are important from the viewpoint of hardware implementation, and (ii) properties that are relevant for target detection. Waveforms that fall in the first category are for example polyphase sequences^{2}^{2}2A polyphase sequence is a sequence whose coefficients are of the form for some , see e.g. [13]., since they give rise to signals with low peaktoaverage power ratio^{3}^{3}3The peaktoaver power ratio of a signal is defined, up to different normalizations, as . (PAPR). A low PAPR is desirable in the digitaltoanalog conversion of signals, since signals with large PAPR would require expensive power amplifiers. Polyphase sequences also have the advantage that they can be very convientiently implemented in hardware via simple lookup tables. The second category usually includes waveforms with low autocorrelation and (nearly) ideal ambiguity function. Quite a number of polyphase sequences, such as Alltop sequences or Kerdock codes, fall in this category. With the exception of [17] a rigorous mathematical theory concerning the benefits (or even optimality) of such sequences has only existed for the detection of a single target. Common to all these carefully constructed sequences in both categories is that they have been designed for singleantenna radar systems and it is a priori not clear at all if any of these sequences are useful in exploiting the potential benefits of a MIMO radar system.
Our paper provides two main contributions: (i) We derive the first rigorous mathematical theory for the detection of moving targets in the azimuthrangeDoppler domain for random sensor arrays. (ii) The transmitted waveforms satisfy a variety of properties that are very desirable and important from a practical viewpoint. In particular, we show that Kerdock sequences, which would perform very poorly in singleantenna radar, are nearly ideally suited for MIMO radar with randomly spaced antennas. Since Kerdock codes are polyphase sequences, they have excellent PAPR and they are easy to implement in hardware via a simple lookup table. Thus, our framework does not just lead to useful theoretical insights, but also has a very strong practical appeal.
1.1. Connections with prior work and innovations
Random sensor arrays have been around for half a century. The pioneering work [27, 28] by Lo contains a mathematical analysis of important specific characteristics of random arrays, such as sidelope behavior and antenna gain. There is extensive engineering literature that deals with random arrays in connection with phased array radar technology, e.g. see [11]. Recently, Carin made an explicit connection between the areas of random sensor arrays and compressive sensing [6]. He has shown that algorithms developed in these two seemingly different areas are in fact highly interrelated. The setup in [6] is quite different from ours, since the paper is only concerned with angular resolution (thus transmission waveforms do not even explicitly enter into the model), while it is often crucial in practice to be able to estimate range and Doppler as well. Moreover, the theoretical analysis in [6] follows more an engineering style and places less emphasis on mathematical rigor. The paper [7] provides interesting results for the angular estimation of stationary targets. Its setup is similar to that in [6], and quite different from ours, as it does not deal with waveform design nor with moving targets.
Kerdock codes have been proposed for radar in [20]. However in the setting of a single transmit antenna. Kerdock codes are known to perform rather poorly^{4}^{4}4This poor performance is caused by Property (ii) in Theorem 3.1. even in the case of single targets as considered in [20]. Only in the setting of mulitple transmit antennas can Kerdock codes exhibit their enormous potential. Our paper utilizes some properties of Kerdock codes proved in [20], but otherwise there is no overlap. In our paper we also present an extension of the main result, that allows for instance the use of the socalled finite harmonic oscillator system as transmission waveforms. These sequences have been derived in [14], where the authors also briefly sketch their use in a singleantenna radar system for the simple case of a single target. Thus, while our framework allows one to employ the finite harmonic oscillator system, there is essentially no overlap of our results with those in [14].
The paper [33] (coauthored by one of the authors) is closest to this paper, but the setting is in a sense complementary. [33] considers a MIMO radar setting with a very specific (nonrandom) choice for the antenna locations, but random waveforms, while the current paper deals with randomly spaced antennas, but very specific, deterministic waveforms. At first glance, the difference may appear to be mainly semantic. But in practice, the second setting has many advantages. From an engineer’s viewpoint random waveforms have several drawback over properly designed deterministic waveforms: they are much harder to implement on a digital device (requiring more complicated hardware, more memory, …); and they exhibit a larger peaktoaveragepower ratio. On the other hand it makes no difference from the viewpoint of physics or hardware, if we place the antennas at random or at deterministic locations. In particular, the current paper yields some important insights, which cannot be inferred from [33]: We obtain a theoretical framework for radar operating with random antenna arrays, a technique which have been around for half a century; we show that Kerdock sequences, which are not useful for SISO or SIMO radar^{5}^{5}5SISO stands for singleinputsingleoutput radar, and SIMO for singleinputmultipleoutput radar (i.e., a radar with one transmit and multiple receive antennas)., are excellent for MIMO radar; our approach allows for waveforms that satisfy a number of properties which are very desirable in practice, and are not satisfied by random waveforms. Indeed, as mentioned above, we also show that the finite harmonic oscillator system “plays well” with random antenna arrays.
Our paper is organized as follows. Section 2 describes the problem setup and the radar model. We review key properties of Kerdock codes in Section 3. Our main theorem is presented in Section 4 and Section 5 is devoted to the proof of the main theorem. In Section 6 we extend our framework to other deterministic waveforms, such as the finite harmonic oscillator system. Numerical simulations are presented in Section 7.
1.2. Notation
For a matrix , we use to denote its adjoint matrix, which is its conjugate transpose. The operator norm of is the largest singular value of and is denoted by . We denote the th column of by and the element in the th row and th column by . The coherence of is defined as
(5) 
The Discrete Fourier Transform (DFT) matrix is written as and the identity matrix as . For , let denote the circulant translation operator, defined by
(6) 
where is understood modulo , and let be the modulation operator defined by
(7) 
Acknowledgements
The authors acknowledge generous support by the National Science Foundation under grant DTRADMS 1042939 and by DARPA under grant N660011114090.
2. Problem Setup
We consider a MIMO radar employing antennas at the transmitter and antennas at the receiver. We assume for convenience that transmitters and receivers are colocated, cf. Figure 1. Furthermore, we assume a coherent propagation scenario, i.e., the element spacing is sufficiently small so that the radar return from a given scatterer is fully correlated across the array. The arrays and all the scatterers are assumed to be in the same 2D plane. The extension to the 3D case is straightforward.
The array manifolds , with randomly spaced antennas are given by
(8) 
and
(9) 
where we assume that the relative antenna spacings ’s and ’s are i.i.d. uniformly on . The th transmit antenna repeatedly transmits the signal , which is assumed to be a periodic, continuoustime signal of periodduration seconds and bandwidth . We observe the backscattered signal over a duration , and since its bandwidth is , it suffices that each receive antennas takes samples^{6}^{6}6Actually the received signal will have a somwhat larger bandwidth due to the Doppler effect. However, in practice this increase in bandwidth is small, so we can assume ., where and . It is convenient to introduce the finitelength vector associated with , via .
Let be the noisefree received signal matrix from a unit strength target at direction , delay , and Doppler (corresponding to its radial velocity with respect to the radar). Then
where is a matrix whose columns are the circularly delayed and Doppler shifted signals .
We let be the noisefree vectorized received signal. We set up a discrete azimuthrangeDoppler grid for , and , where and denote the corresponding discretization stepsizes. Using vectors for all grid points we construct a complete response matrix whose columns are for and , . In other words, is a matrix with columns
(10) 
Assume that the radar illuminates a scene consisting of scatterers located on points of the grid. Let be a sparse vector whose nonzero elements are the complex amplitudes of the scatterers in the scene. The zero elements correspond to grid points which are not occupied by scatterers. We can then define the radar signal received from this scene by
(11) 
where is an vector, is an sparse vector and is an complex Gaussian noise vector. Our goal is to solve for , i.e., to locate the scatterers (and their reflection coefficients) in the azimuthdelayDoppler domain.
Remark: The assumption that the targets lie on the grid points, while common in compressive sensing, is certainly restrictive. A violation of this assumption will result in a model mismatch, sometimes dubbed gridding error, which can potentially be quite severe [18, 8]. Recently some interesting strategies have been proposed to overcome this gridding error [10, 35]. But these methods – at least in their current form – are not directly applicable to our setting. This model mismatch issue is beyond the scope of this paper and will be addressed in our future research.
3. Kerdock codes
In this section we introduce one particularly useful set of transmission waveforms. Due to the setup in Section 2 it suffices that we deal with discrete, finitelength signals as transmission waveforms. We briefly review the construction of Kerdock codes and some of their fundamental properties. There is a long list of properties that radar waveforms should satisfy. As we will see in this paper, Kerdock codes fulfill many of them. Kerdock codes over (i.e., binary Kerdock codes) were originally introduced in [23]. In the seminal paper [4] the authors extend Kerdock codes from to . By doing so, they uncover many fascinating properties of Kerdock codes and reveal numerous deep connections between coding theory, discrete geometry and group theory. In the same paper, the authors also extend Kerdock codes to the setting of , where is an odd prime.
Kerdock codes are an example of socalled mutually unbiased bases [39, 34]. Kerdock codes have also been proposed for use in communications engineering [16, 22]. In [20] the authors suggest the use of Kerdock codes for radar, based on the peculiar properties of the discrete ambiguity function associated with Kerdock codes. We emphasize however that for the single transmit antenna radar scenario Kerdock codes would actually perform rather badly, as discussed after Theorem 3.1 and shown by Figure 5 in Section 7. It is only in the setting of multiple transmit antennas that Kerdock codes become useful for radar.
For the remainder of this paper we will only be concerned with Kerdock codes over . Some of the Kerdock codes over , namely those corresponding to desarguesian planes in the language of [4], have also been derived earlier in [25] and [24]. A simple way to construct these Kerdock codes is the following, in which they arise as eigenvectors of timefrequency shift operators. Let be an odd prime number. For each we compute the eigenvector decomposition of (which always exists, since is a unitary matrix)
(12) 
where the unitary matrix contains the eigenvectors of and the diagonal matrix the associated eigenvalues^{7}^{7}7The attentive reader will have noticed that is just the DFT matrix .. Furthermore, we define . Now, let be the th column of . The set consisting of the vectors forms a Kerdock code. There are numerous equivalent ways to derive this Kerdock code, but, as pointed out earlier, not all Kerdock codes over are equivalent (see also the comment following Corollary 11.6 in [4]). But we will be a bit sloppy, and simply refer to the Kerdock code constructed above as the Kerdock code.
In the following theorem we collect those key properties of Kerdock codes that are most relevant for radar. These properties are either explicitly proved in [4, 20] or can be derived easily from properties stated in those papers.
Theorem 3.1.
Kerdock codes over , where is an odd prime, satisfy the following properties:

Mutually unbiased bases: For all and all , there holds:

Timefrequency “autocorrelation”:
(a) For any fixed there exists a unique such that(13) (14) (b) For any fixed , there exist , such that
(15) 
Timefrequency crosscorrelation: For all and all and there holds:
(16) 
Polyphase property (Roots of unity property) in time and in frequency:
For any , there holds:(17) For any , there holds:
(18)
Proof.
Property (i) is proved for instance in Lemma 11.3 in [4]. Properties (ii) and (iii) appear in Theorem 3 of [20]. Statement (17) of property (iv) follows from the comment right after Corollary 11.6 in [4]. Finally, statement (18) of property (iv) follows from (12) together with property (3) and the wellknown fundamental relationships
∎
Kerdock codes have been proposed for adaptive radar in [20]. We emphasize again though that Kerdock codes would not be very effective for a radar system with a single transmit antenna (SISO or SIMO radar). This can be easily seen as follows: Assume we only have one antenna that transmits one waveform . Because of (15), is (up to a constant phase factor) equal to for some . In practice this ambiguity prevents us from determining the distance and the velocity of the object, when using Kerdock codes for SISO or SIMO.
As a consequence of the aforementioned ambiguity we will not use all of the Kerdock codes as transmission signals for our MIMO radar, instead we will choose one code for each index . The reason is that we need the waveforms to have low timefrequency crosscorrelation, while (16) only holds when and are different.
Definition 3.2 (Kerdock waveforms).
Let be a Kerdock code over . The Kerdock waveforms , where , are given by for some arbitrary . In other words, for each we pick an arbitrary vector from the orthonormal basis .
Note that Kerdock waveforms do not include any unit vectors, since only the first unitary matrices are considered and is strictly less than (recall that ).
4. The main result
As mentioned in the introduction, a standard approach to solve (11) when is sparse, is
(19) 
which is also known as lasso [36]. But instead of (19), we will use the debiased lasso. That means first we compute an approximation for the support of by solving (19). This is the detection step. Then, in the estimation step, we “debias” the solution by computing the amplitudes of via solving the reducedsize least squares problem , where is the submatrix of consisting of the columns corresponding to the index set , and similarly for .
We assume that the locations of the targets are random. To be precise, we assume that the nonzero coefficients of are selected uniformly at random and the phases of the nonzero entries of are random and uniformly distributed in . We will refer to this model as the generic sparse model.
We are now ready to state our main result.
Theorem 4.1.
Consider , where is defined as in (10) and . Assume that the positions of the transmit and receive antennas ’s and ’s are chosen i.i.d. uniformly on at random. Suppose further that each transmit antenna sends a different Kerdock waveform, i.e. the columns of the signal matrix are different Kerdock waveforms. Choose the discretization stepsizes to be ,, and suppose that
(20) 
and also
(21) 
If is drawn from the generic sparse scatterer model with
(22) 
for some constant , and if
(23) 
then the solution of the debiased lasso computed with satisfies
(24) 
with probability at least
(25) 
and
(26) 
with probability at least
(27) 
where
and
Remarks:

The condition in (21) is by no means necessary, but rather to make our computation a little cleaner. We could change it into , then the theorem would remain true with a slightly different probability of success.
5. Proof of the result
To prove Theorem 4.1, we use a theorem by Candès and Plan (Theorem 1.3 in [5]) which requires to estimate the operator norm of and the coherence of . The original theorem only treats the realvalued case, it can be extended to complexvalues case after some straightforward modifications (see Appendix B in [33]).
5.1. Auxiliary results
We first need the following Bernstein type lemma.
Lemma 5.1.
Suppose is an matrix, and are two joint independent random vectors in with zero means and for . If is a positive constant, then for any and ,

if for all , then
(28) and
(29) 
if for and , then
(30) and
(31)
Proof.
where denotes addition modulo .
Let us first assume that .
Since and are joint independent, then for any , the entries in are all joint independent and it is easy to check that and , then Theorem 4.5 in [21] will give,
(32) 
different from above, the entries in are no longer all jointly independent. But similar to the proof of Theorem 5.1 in [30] and Lemma 3 in [33], we observe that for any we can split the index set into two subsets , each of size , such that the variables are jointly independent for , and analogous for . (For convenience we assume here that is even, but with a negligible modification the argument also applies for odd .) In other words, each of the sums , contains only jointly independent terms.
So for each ,
(34) 
which implies that
(35) 
Now let us assume that for and .
Since and are joint independent and ,
(37) 
5.2. Estimation of the Operator Norm
Proof.
Since , we consider matrix as block matrix
where the blocks are matrices of size .
Via a simple permutation, we can turn into a matrix with blocks of size , where the th entry of the block is defined as
(41) 
Then it is easy to see that is blockdiagonal, and all the diagonalblocks are identical. So we only have to bound the first block .
Define , then
where is the matrixvalued random variable given by and therefore .
We apply the union bound over the possibilities associated with and get
which implies that
Then the fact that will give us the conclusion.
∎
5.3. Estimation of the Coherence
Lemma 5.3.
Proof.
We need to find an upper bound for
Recall the , it follows from the definition that
from which we readily compute
We use the discretization , , where , , with .
Since
for