SpringerOpen Newsletter

Receive periodic news and updates relating to SpringerOpen.

This article is part of the series Radar and Sonar Networks.

Open Access Research

Compressive sensing in distributed radar sensor networks using pulse compression waveforms

Lei Xu1*, Qilian Liang1, Xiuzhen Cheng2 and Dechang Chen3

Author Affiliations

1 Department of Electrical Engineering, University of Texas at Arlington, Arlington, TX 76010, USA

2 Department of Computer Science, The George Washington University, Washington DC 20052, USA

3 Department of Preventive Medicine and Biometrics Uniformed Services, University of the Health Sciences Bethesda, Maryland 20814-4799, USA

For all author emails, please log on.

EURASIP Journal on Wireless Communications and Networking 2013, 2013:36  doi:10.1186/1687-1499-2013-36

The electronic version of this article is the complete one and can be found online at: http://jwcn.eurasipjournals.com/content/2013/1/36

Received:25 October 2012
Accepted:5 January 2013
Published:19 February 2013

© 2013 Xu et al.; licensee Springer.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License(http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


Inspired by recent advances in compressive sensing (CS), we introduce CS to the radar sensor network (RSN) using pulse compression technique. Our idea is to employ a set of stepped-frequency (SF) waveforms as pulse compression codes for transmit sensors, and to use the same SF waveforms as the sparse matrix to compress the signal in the receiving sensor. We obtain that the signal samples along the time domain could be largely compressed so that they could be recovered by a small number of measurements. A diversity gain could also be obtained at the output of the matched filters. In addition, we also develop a maximum likelihood (ML) algorithm for radar cross section (RCS) parameter estimation and provide the Cramer-Rao lower bound (CRLB) to validate the theoretical result. Simulation results show that the signal could be perfectly reconstructed if the number of measurements is equal to or larger than the number of transmit sensors. Even if the signal could not be completely recovered, the probability of miss detection of target could be kept zero. It is also illustrated that the actual variance of the RCS parameter estimation <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M1">View MathML</a> satisfies the CRLB and our ML estimator is an accurate estimator on the target RCS parameter.

Compressive sensing; Radar sensor networks; Pulse compression; Stepped-frequency waveform; Target RCS

1 Introduction

Current requirements in warfighting functionality result in obtaining accurate and timely information about battlespace objects and events so that the warfighters can make decision about reliable location, tracking, combat identification and targeting information. While massive amounts of data will be generated by a penetrating sensor, it is important for the warfighters to find technologies that not only integrate information from diverse sources but also provide indications of information significance in ways that help them to make tactical decision. The RCS is the property of a scattering object, or target, which represents the magnitude of the echo signal returned to the radar by the target. Hence, we could have different classes with different RCS values representing corresponding targets, such as bird, conventional unmanned winged missile, small single-engine aircraft and large flight aircraft. In this article, we will study the target RCS in a radar sensor network (RSN) by using compressive sensing techniques.

It is well known that wireless sensor networks (WSN) are a fast growing class of systems. In [1], the authors presented a new method that makes use of the properties of data of multiple sensors to enable reliable data collection. In [2], the authors adopted a mutual-information-based sensor selection (MISS) algorithm to help sensing devices collaborate among themselves to improve the target localization and tracking accuracies. Alike WSN, RSN has been recently considered to overcome the performance degradation of a single radar. In [3], the authors design a network of distributed radar sensors that work in an ad hoc fashion and the simulation results showed that proposed RSN can provide much better detection performance than that of single radar sensor. However, RSN is quite different from WSN. The waveform of each radar sensor has to be properly designed, otherwise, these radar sensors are likely to badly interfere with each other in the RSN. As a result, the design of radar waveforms has been under the study with the goal of optimizing the performance of the RSN. In [4], binary coded pulses using simulated annealing in RSN are proposed. Liang [5] studied constant frequency (CF) pulse waveform design and proposed Maximum-Likelihood (ML) automatic target recognition (ATR) approach for both nonfluctuating and fluctuating targets in a network of multiple radar sensors. Furthermore, Liang [6] studied the RSN design based on linear frequency modulation (LFM) waveforms and also applied the LFM waveforms to RSN in the context of ATR with delay-Doppler uncertainty. In addition, it is known that the pulse compression technique allows a radar to achieve both the energy of a long pulse and the resolution of a short pulse, without the high peak power which is required by a high energy short duration pulse [7]. Pulse compression waveforms are obtained by adding frequency or phase modulation to a simple pulse. A stepped-frequency (SF) waveform is a frequency modulation waveform for obtaining a large bandwidth, and thus a fine range resolution without requiring intrapulse frequency modulation. The most common SF waveform employs a linear frequency stepping pattern, where the RF frequency of each pulse is increased by ΔF. This representation motivates the applicability of the recently proposed compressive sensing (CS) theory [8,9] that refers to such signals as ‘sparse’ or ‘compressible’.

Recent results in compressive sensing have shown that the information from a signal may be captured with a small set of nonadaptive, linear measurements as long as the signal is sparse in some basis or frame [10-12]. If the signal is properly chosen, the number of measurements to recover the signal can be much smaller than the number of Nyquist-rate samples. The application of compressive sensing to radar imaging has been investigated in [13,14]. A CS-based data acquisition and imaging method was proposed to study a number of point-like targets for SF continuous wave ground penetrating radars (SFCW-GPRs) in [15]. In [16], the authors proposed the step-frequency with compressive sampling (SFR-CS), that achieves high target range and speed resolution using significantly smaller bandwidth than transitional step-frequency radar. The application of CS to multiple input and multiple output (MIMO) radar has also attracted a lot of attention in recent research. The study of [17] studied angle-Doppler estimation of multiple targets for MIMO radar system. In addition, the authors in [18] also considered the range estimation performance of their proposed method for MIMO radar systems that employ CS. Based on adaptive radar design, the authors in [19] studied MIMO radar with widely separated antennas in the context of sparse modeling for estimating the positions and velocities of multiple targets.

Due to the expansion of data introduced to RSN, the compression and reconstruction of the received data is a design challenge of future RSN. Unlike the above research, the study in this article explores how to exploit compressive sensing in RSN composed of a number of transmit sensors but only one receiving sensor. It is known that the signal must be ‘compressible’ for compressive sensing to have benefit. Recognizing that the SF train could act as the sparsity basis for the signal, we apply it as a pulse compression code to construct the ‘compressible signal’ for a transmit sensor. We choose the Gaussian matrix as the measurement matrix that satisfies the restricted isometry property (RIP) with this basis. However, there are still a number of challenges in most CS theory applied to radar which specifically mentions that it eliminates the need for matched filter in the radar receiving sensor. In [13], only the range parameter of the target was studied and the target reflectivity being probed must be compressible in some basis before their CS-based radar system could work. In [20], the transmitted signal must be sufficiently ‘incoherent’ and the targets have to be radially aligned with the transmitter and receiver. In this article, we propose and investigate a totally different model of CS-based RSN system. Hence, the matched filters are still used in the receiving sensor.

After applying compressive sensing to RSN, we also perform target RCS value estimation. We propose an ML algorithm to estimate the target RCS parameter and use the Cramer-Rao lower bound (CRLB) to validate our theoretical result. In the simulation parts, the performance of signal recovery and the performance of target detection are studied as well as the performance of target RCS value estimation. The simulation results show that the signal could be precisely recovered if the number of measurements is no less than the number of sensors in RSN. The target could be perfectly detected even if the signal could not be precisely recovered. As a result, much smaller measurement matrix could be used on the receive part for the purpose of target detection. Finally, the actual variance of the RCS parameter estimation <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M2">View MathML</a> satisfies the CRLB.

The article is organized as follows. In Section 2, we present our basic model by applying CS to RSN exploiting the pulse compression technique. In Section 3, we show that the range resolution is increased for RSN using Step-Frequency waveforms as the pulse compression codes. In Section 4, we propose a Maximum Likelihood (ML) algorithm for target RCS value estimation and derive the CRLB accordingly. In Section 5, the simulation results illustrate the performance of signal recovery after decompression, target detection and target RCS value estimation. Finally, conclusions are drawn in Section 6.

2 The basic model

2.1 The produced signal for compressed sensing

The SF pulse train is a pulse burst waveform which obtains large overall bandwidth while maintaining narrow instantaneous bandwidth. Each pulse in the burst is a simple, constant-frequency pulse; however, the RF frequency is added by a frequency step ΔF between consecutive pulses. The most common SF waveform employs a linear frequency stepping pattern, where the RF frequency of each pulse is increased by ΔF Herts from the preceding pulse. In addition, the frequency steps can be added to a train of unmodulated pulses, as well as to a train of modulated pulses. The SF waveform we study in this article is expressed as follows:

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M3">View MathML</a>


where Tp is the pulse repetition interval and I is the number of pulses in a burst (train) of pulses. Observing the right side of Equation (1), we see that <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M4">View MathML</a> can be expressed into the orthogonal Fourier basis <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M5">View MathML</a>, where <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M6">View MathML</a> and n is the index of a transmit sensor in a RSN. Here, we let the frequency step ΔF among each pulse burst waveform corresponding to each specific transmit sensor be different and related to the index n. The goal of CS is to perform good reconstruction of the signal by using a few sparse samples. Here, we construct the data in such a way that is compressed already. We establish a RSN which uses the SF waveforms as pulse compression codes. The RSN consists of N transmit sensors indexed by n and only one receiving sensor. Then <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M7">View MathML</a> denotes the transmit signal for the transmit sensor n. Let <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M8">View MathML</a> and <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M9">View MathML</a> denote the information signal and the pulse compression waveforms, where N is the number of transmit sensors and I is still the number of pulses in a burst. In order for the later use, we express the preceding equations in vector format that <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M10">View MathML</a> and <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M11">View MathML</a>. Then, S = ΨXT is the transmit matrix, where <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M12">View MathML</a> and <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M13">View MathML</a>.

The RCS is the property of a scattering target that is included in the radar equation to represent the magnitude of the echo signal returned to the radar by the target. As a result, the RCS returned to the receiving sensor by the target is assumed to have isotropic reflectivity modeled by zero-mean, unit-variance, independent and identically distributed (i.i.d.) Gaussian complex random variables λn. Assuming that it is the slow fading system, where the amplitude and phase change imposed by the channel can be considered roughly constant over the period of use, we suppose that the Rayleigh distortion does not fluctuate during the whole pulse burst. The RCS can be modeled by the diagonal matrix

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M14','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M14">View MathML</a>


where the normalization factor makes the target average <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M15">View MathML</a> independent of the number of transmit sensors in the model. A popular method for representing the fluctuations of targets are the four statistical models described by Swerling in [21]. The nonfluctuating target modeled using non-zero constants for λn=λ is identified as ‘Swerling 0’ or ‘Swerling 5’ model. For the fluctuating target, if |λn| is drawn from the Rayleigh P.D.F. and varies independently from path to path, the target model represents a classical ‘Swerling 2’ model.

Taken all the above parameters into account, we use the following block diagram illustrating as Figure 1 to show our basic model.

thumbnailFigure 1. The block diagram of the model.

A vector form of the received signal including the useful signal and the system noise is <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M16">View MathML</a> and we will process with this <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M17','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M17">View MathML</a> in the following section

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M18','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M18">View MathML</a>


<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M19','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M19">View MathML</a>


where <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M20','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M20">View MathML</a> is the noise vector.

2.2 Decomposition and recovery of the signal

In order to obtain some benefit form CS, we choose the same orthogonal basis <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M21','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M21">View MathML</a> as the sparse matrix Ψ to study the received signal <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M22','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M22">View MathML</a> including the useful signal and the noise, where

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M23','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M23">View MathML</a>


with <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M24','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M24">View MathML</a>. The CS can only work when there are K non-zero coefficients in the vector <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M25','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M25">View MathML</a>, where K<I. According to the concepts of CS, there should be only a few large coefficients and many small ones in <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M26','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M26">View MathML</a>.

The CS approach uses a set of measurement probes <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M27','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M27">View MathML</a>, where M is significantly less than the dimensionality I of each individual probe (I is just the size of receive signal <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M28','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M28">View MathML</a> in our case). In a different way, we would like to recover all the I coefficients of <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M29','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M29">View MathML</a> by observing or measuring a subset M of these

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M30','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M30">View MathML</a>


where Y=[y1,y2,…,yM]T, M⊂1,…,I. Since random matrix Φ with independent identically distributed (i.i.d.) entries, e.g., Gaussian or binary entries, exhibit a very low coherence with any fixed sparse matrix Ψ[13], we take the Gaussian matrix as the measurement matrix.

The recovery of <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M31','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M31">View MathML</a> is done by solving a constrained l1-norm minimization problem and the reconstructed <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M32','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M32">View MathML</a> is given by <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M33','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M33">View MathML</a>, where <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M34','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M34">View MathML</a> is the solution to the convex optimization program

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M35','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M35">View MathML</a>


We choose the solution whose coefficient sequence has the minimal l1 norm.

We describe how our model recovers a particular <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M36','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M36">View MathML</a> in section 2.3. We first produce the transmit signal and receive it as <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M37','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M37">View MathML</a>. Then, we compute Y by compressing the received signal <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M38','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M38">View MathML</a>. Finally, we use the l1 minimization with relaxed constraints to reconstruct the original signal <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M39','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M39">View MathML</a>.

2.3 The basic model


<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M40','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M40">View MathML</a>

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M41','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M41">View MathML</a> for all mM

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M42','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M42">View MathML</a>, subject to

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M43','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M43">View MathML</a>

Output:<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M44','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M44">View MathML</a>

We explicitly tie together the parameters for the decompressive process, i.e.,

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M45','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M45">View MathML</a>


One of the important properties that A=ΦTΨ should satisfy is the restricted isometry property (RIP) [11]. A sufficient condition for the RIP is that the measurement vectors and the sparsity basis must be incoherent with each other. The coherence between the measurement matrix Φ and the sparsity basis Ψ here is defined as

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M46','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M46">View MathML</a>


If Φ and Ψ contain uncorrelated elements, the coherence is large. So the ‘incoherency’ here means that the inner products between the probes and the sparsity basis vectors are small, or, in other words, μ is small.

In this article, we express the SF coefficients in the method of Fourier basis as the sparse matrix

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M47','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M47">View MathML</a>


Here, the sparse matrix can be non-squared, i.e., the number of rows can be different from the number of columns in the sparse matrix. The number of rows is based on the number of transmit sensors N and the number of columns depends on the number of time samples I. The reason why we can use such a non-squared matrix is that we use the matrix same as the sparse matrix to produce the ‘compressible’ transmit signal. In addition, the random waveforms φm(t) with independent identically distributed (i.i.d) vectors are largely incoherent with the fixed basis Ψ. Hence, we choose the complex Gaussian random vectors as the measurement matrix in our study.

2.4 The output of the matched filter

Different from previous CS-based radar system, we proposed a totally new model for CS-based RSN. Since we use the matched filters in the receiving sensor, let us assume that the received signal could be well recovered as <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M48','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M48">View MathML</a>

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M49','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M49">View MathML</a>

where <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M50','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M50">View MathML</a> is the noise vector. Actually, the power of noise <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M51','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M51">View MathML</a> in the recovered signal is less than the system noise in the received signal in Equation (3). Since the <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M52','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M52">View MathML</a> may not have good sparsity property in the basis Ψ, some information of <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M53','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M53">View MathML</a> could not be exactly recovered by <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M54','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M54">View MathML</a> Here, we choose the matched filters corresponding to the transmit signal S=ΨXT. Processing by a bank of matched filters, we express the output of the matched filter as the following matrix operations

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M55','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M55">View MathML</a>


Since Ψ is the Fourier basis in this article, and X is information bits expressed as pulsed signals, it is easy to obtain that

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M56','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M56">View MathML</a>


where E is the energy of the burst of pulses. Therefore,

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M57','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M57">View MathML</a>


where λn is the radar cross section for nth transmit sensor. According to the Equation (12), all the N radar cross section parameters could be taken into account to detect or recognize the target, so that the diversity gain of transmit side can be obtained as well.

3 Increased range resolution

The range resolution improvement is one of the important properties of the linear SF waveforms (LSFWs). The details of the Doppler response and range resolution can be expressed by the ambiguity function (AF) of the LSFW [22]. The AF is defined as

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M58','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M58">View MathML</a>


We designate the complex envelope of the SF pulse train by u(t) given by

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M59','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M59">View MathML</a>


where u1(t) is a radar pulse waveform. Figure 2 presents a sample AF of a SF train of unmodulated pulses calculated by using the MATLAB tool.

thumbnailFigure 2. Ambiguity function of SF train of unmodulated pulses.

It is obvious to see that the range resolution of the signal is improved, but there are still prominent sidelobes in delay and ambiguity in Doppler. As a result, LFMs and SF can be combined to mitigate the raging lobes, i.e., SF train of LFM pulses could be used. An example of SF train of LFM is shown in Figure 3.

thumbnailFigure 3. Ambiguity function of SF train of LFM.

We compare the AF of a SF train of LFM pulses, as shown in Figure 3, with the AF of a SF train of unmodulated pulses as seen in Figure 2. Clearly, by adding the LFM, the range and the Doppler resolutions are improved by canceling the sidelobes along the delay and Doppler axes. As a result, we choose to use the SF train of LFM pulses as the transmit waveforms in our model to obtain both the range and Doppler resolution gain.

4 Target RCS value estimation

In this section, we use the maximum-likelihood (ML) estimation algorithm to perform target radar cross section (RCS) parameter estimation [23] in the proposed RSN model. For the ‘Swerling 2’ model, the RCS voltage |λ(u)| follows a Rayleigh distribution and the I and Q subchannels of λ(u) follow zero-mean complex Gaussian distribution with a variance γ2 (the RCS average power value)

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M60','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M60">View MathML</a>


In addition, n(u)=nI(u)+jnQ(u) follows a zero-mean complex Gaussian distribution with a variance σ2 for each I and Q subchannel. We express Equation (12) as following

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M61','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M61">View MathML</a>


Here, we let

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M62','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M62">View MathML</a>


Since λn(u) and <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M63','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M63">View MathML</a> are zero-mean complex Gaussian random variables, <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M64','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M64">View MathML</a> is zero-mean complex Gaussian random variable with a variance E2γ2+σ2.

Assuming that yn≅|Zn(u)| follows a Rayleigh distribution:

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M65','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M65">View MathML</a>


Let <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M66','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M66">View MathML</a>, <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M67','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M67">View MathML</a> and yn=|Zn(u)|, we can obtain that y≅ [y1,y2,…,yN]. Assuming that yn are independent of each other, then the P.D.F. of Y is

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M68','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M68">View MathML</a>


If θγ2, we can express (18) as

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M69','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M69">View MathML</a>


Therefore, we represent the ML algorithm to estimate the RCS average value θ as

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M70','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M70">View MathML</a>


It is equivalent to maximize log f(y) (natural logarithm),

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M71','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M71">View MathML</a>


It is a continuous function for yn>0 and θ>0, hence, a necessary condition for the ML estimation is

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M72','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M72">View MathML</a>


Equation (24) has the unique solution

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M73','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M73">View MathML</a>


Since θ>0,

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M74','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M74">View MathML</a>



<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M75','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M75">View MathML</a>


this solution gives the unique maximum of logf(y). The expectation of <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M76','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M76">View MathML</a> is then

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M77','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M77">View MathML</a>


The mean value of yn is <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M78','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M78">View MathML</a>, and its variance is (4−Π)(E2θ+σ2)/2. Since yn are independent of each other, it is

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M79','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M79">View MathML</a>



<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M80','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M80">View MathML</a>


As a result, it is an unbiased estimator.

Fisher’s information [24] in this case can be obtained as

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M81','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M81">View MathML</a>


Taking Equation (28) into account, we can obtain the Cramer-Rao lower bound (CRLB) [24]

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M82','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M82">View MathML</a>


From (31), we observe that CRLB is inversely proportional to the number of radars N in the RSN, which means that the RSN with large N will have a low CRLB. We draw this conclusion by assuming that the radar pulses are independent (in time and space) and follow a Rayleigh distribution, according to the ‘Swerling 2’ model.

5 Simulation results

5.1 Signal recovery

In this section, we study the scenario where we have N transmit sensors but only one receiving sensor. We assume that the number of samples is set to 500. We apply the SF train of LFM pulses as pulse compression codes on the transmit side and use the CS technique in the receiving side. The classical ‘Swerling model 2’ and Gaussian noise are also considered in the simulation. The reconstructed signal is compared with the original signal by calculating the mean square error (MSE) in order to evaluate the reconstruction ability. The Mean Square Error between the original signal and the reconstructed signal is shown in Figures 4 and 5 as a function of different number of measurements M and SNR values. The noise considered here is introduced by the propagation in the air but not by compressing and decompressing process. We use the Monte-Carlo simulation model here and the results are averaged by 105 runs/iterations. The cases of N=50 and N=100, where N is the number of transmit sensors are illustrated in Figures 4 and 5, separately.

thumbnailFigure 4. Normalized MSE between reconstructed signal and original signal for fixed N = 100. (a) Normalized MSE versus M. Each line corresponds to a different SNR (−5,4). (b) Normalized MSE versus SNR. Each line corresponds to a different number of measurements M.

thumbnailFigure 5. Normalized MSE between reconstructed signal and original signal for fixed N = 50. (a) Normalized MSE versus M. Each line corresponds to a different SNR (−5,4). (b) Normalized MSE versus SNR. Each line corresponds to a different number of measurements M.

According to both Figures 4a and 5a, MSE is reduced as the number of measurements M is increased. The system can perfectly reconstruct the signal which includes the received signal and the system noise when the number of measurements M is equal to the number of transmit sensors N. In addition, the slope of MSE versus the number of measurements M is almost a consistent for each SNR value. From Figures 4b and 5b, we draw the same conclusion that the closer the number of measurements approaches N (MN), the better performance of signal recovery is achieved. In addition, we also discover that the MSE does not depend much on the SNR, especially when M is large. As a result, the proposed model can be used under a low SNR if the number of measurements M could be properly chosen according to the number of transmit sensors N.

On the basis of the simulation results, we can draw a brief conclusion that the number of measurements M of our model only depends on the size of RSN even when the number of samples is fixed as large as 500 here. Another important result emerging from the simulations is that the probability of target miss detection is zero no matter how small a number of measurements we use in the recovery process. That is to say, less measurements can be used to detect the target in the system, since a kind of diversity gain is achieved at the output of the matched filters on the receiving sensors.

5.2 RCS parameter estimation

In this section, we will consider the fluctuating target with an RCS parameter θ (following Rayleigh distribution) investigated in the Section 4. We will apply the ML estimation algorithm to estimate the parameter <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M83','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M83">View MathML</a>. The scenario is similar to the one in the section above, but the number of samples in time domain is reduced to 100 for complexity reasons. We ran Monte Carlo simulations for 105 iterations at each SNR value. We have considered the fluctuating target with RCS parameter θ=2 (Small flighter aircraft or 4 passenger jet) in Figure 6. We plotted the variance of the RCS ML estimator with different number of radars in RSN.

thumbnailFigure 6. Variance of RCS ML estimator with different number of radars in RSN. The actual RCS parameter θ=2.

According to Figure 6, the variance of <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M84','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M84">View MathML</a> closely approaches the CRLB but doesn’t exactly match it. The reason why the variance of <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M85','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M85">View MathML</a> is not exactly the same as the CRLB is that the noise included in Equation (10) cannot be termed as sparse as the transmit signal in the same sparsity matrix. Therefore, the noise which is considered in CRLB can not be perfectly reconstructed after the decompression. The power of noise σ2 reduces so that the calculated CRLB might be lower than the practical one. As a result, a gap between the variance of <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M86','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M86">View MathML</a> and the CRLB. However, the gap between variance of <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M87','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M87">View MathML</a> and the CRLB is reduced by increasing the number of radars. It is also easy to see that the actual variance of <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M88','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M88">View MathML</a> reduces as the number of radars increases from N=10 to N=20. Hence, the actual variance of <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M89','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M89">View MathML</a> is inversely proportional to N, as we have shown in the theoretical result <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M90','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M90">View MathML</a>.

It is easy to see that the actual variance of <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M91','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/36/mathml/M91">View MathML</a> and the CRLB do not change much as the SNR increases. Stating differently, our ML estimator performs well even for low SNR ratios. In all, the simulation results validate the theoretical results. The variance of the RCS parameter estimation satisfies the CRLB and our ML estimator on the RCS parameter is an accurate estimator.

6 Conclusions

Motivated by the representation of SF waveforms, we introduced CS to the RSN exploiting the pulse compression technique. A set of SF waveforms were applied as pulse compression codes at the transmit sensors, and the sparse matrix is also constructed based on the same SF waveforms. We observed that the signal samples along the time domain can be significantly compressed and recovered by using a small number of measurements which depend on the number of transmit sensors. A diversity gain is also achieved after the matched filters in the proposed model, so the probability of target miss detection can be zero even if the signal could not be perfectly recovered. In addition, we propose a ML algorithm to estimate the target RCS parameter and use the CRLB to successfully verify our theoretical result.

Competing interests

The authors declare that they have no competing interests.


This study was supported in part by National Science Foundation under Grants CNS-0964713, CNS-1017662, CNS-0963957, CNS-0964060, and Office of Naval Research under Grant N00014-11-1-0071.


  1. S Mukhopadhyay, C Schurgers, D Panigrahi, S Dey, Model-based techniques for data reliability in wireless sensor networks. IEEE Trans. Mob. Comput 8(4), 528–543 (2008)

  2. T Onel, C Ersoy, H Delic, Information content-based sensor selection and transmission power adjustment for collaborative target tracking. IEEE Trans. Mob. Comput 8(4), 1103–1116 (2009)

  3. J Liang, Q Liang, Design and analysis of distributed radar sensor networks. IEEE Trans. Parallel Distrib. Process 22(11), 1926–1933 (2011)

  4. H Deng, Synthesis of binary sequences with good correlation and cross-correlation properties by simulated annealing. IEEE Trans. Aerosp. Electron. Systs 8(8), 684–689 (2009)

  5. Q Liang, Waveform design and diversity in radar sensor networks: theoretical analysis and application to automatic target recognition. IEEE Sensor Ad Hoc Commun. Netws. Conf 2(28), 684–689 (2006)

  6. Q Liang, Radar sensor networks for automatic target recognition with Delay-Doppler uncertainty. IEEE Military Commun. Conf 23–25, 1–7 (2006)

  7. MA Richards, Fundamentals of Radar Signal Processing (McGraw-Hill, 2005)

  8. E Candes, M Walkin, An introduction to compressive sampling. IEEE Signal Process. Mag 25(2), 21–30 (2008)

  9. R Baraniuk, Compressive sensing. IEEE Signal Process. Mag 24(4), 118–121 (2007)

  10. E Candes, J Romberg, T Tao, Robust uncertainty principles: exact signal reconstruction from highly incomplete frequency information. IEEE Trans. Inf. Theory 52(2), 489–509 (2006)

  11. D Donoho, Compressed sensing. IEEE Trans. Inf. Theory 52, 1289–1306 (2006)

  12. E Candes, T Tao, Near-optimal signal recovery from random projections: universal encoding strategies. IEEE Trans. Inf. Theory 52(12), 5406–5425 (2006)

  13. R Baraniuk, P Steeghs, Compressive radar imaging. Radar Conf., 2007 IEEE, 128–133 (2007)

  14. AC Gurbuz, JH McClellan, WR Scott, Compressive sensing for GPR imaging. in Proc, ed. by . 41th Asilomar Conf. Signals, Syst. Comput (Pacofoc Grove, CA, 2007), pp. 2223–2227

  15. AC Gurbuz, JH McClellan, WR Scott, A compressive sensing data acqisition and imaging method for stepped frequency GPRs. IEEE Trans. Signal Process 57, 2640–2650 (2009)

  16. S Shah, Y Yu, A Petropulu, Step-frequency radar with compressive sampling (SFR-CS). 2010 IEEE International Conference on Acoustics Speech and Signal Processing (ICASSP), 1686–1689 (2010)

  17. Y Yu, AP Petropulu, HV Poor, MIMO radar using compressive sampling. IEEE J. Sel. Top. Signal Process 4(1), 146–163 (2010)

  18. Y Yu, AP Petropulu, HV Poor, Range estimation for MIMO step-frequency radar with compressive sensing. 2010 4th International Symposium on Communications, Control and Signal Processing (ISCCSP) (2010)

  19. S Gogineni, A Nehorai, Adaptive design for distributed MIMO radar using sparse modeling. 2010 International Waveform Diversity and Design Conference (WDD) (2010)

  20. MA Herman, T Strohmer, High-resolution radar via compressed sensing. IEEE Trans. Signal Process 57(6), 2275–2284 (2009)

  21. P Swerling, Probability of detection for fluctuating targets. IRE Trans. Inf. Theory 6, 269–308 (1960). Publisher Full Text OpenURL

  22. N Levanon, E Mozeson, Radar Signals (New York: Wiley, 2004)

  23. Q Liang, X Cheng, KUPS: knowledge-based ubiquitous and persistent sensor networks for threat assessment. IEEE Trans. Aerosp. Electron. Syst 44(3), 1060–1069 (2008)

  24. JM Mendel, Lessons in Estimation Theory for Signal Processing, Communications, and Control (Upper Saddle River, NJ: Prentice-Hall, 1995)