# Adaptive clutter suppression for airborne random pulse repetition interval radar based on compressed sensing.

1. INTRODUCTIONRandom selection of the pulse repetition intervals (PRI) within fixed coherent process interval (CPI) as a means of smearing the ambiguity peaks has been alluded to in the past [1-3], which can achieve great ECCM capabilities. Traditional processing approach for the random PRI radar signal such as correlation is constrained by the uncertainty principle and suffers from high sidelobe pedestal, which reduces the performance of target detection and velocity measurement. Considering of the random modulation scheme and the inherent sparsity of the target echo signal, the compressed sensing (CS) theory [4-8] has been applied to successfully resolve the velocity ambiguity and suppress the sidelobe pedestal [9].

In some actual applications such as airborne systems, however, echoes form land clutter consisting of stable objects like mountains and buildings [10] can violate the assumptions for CS severely [11], which needs to be investigated. By now, several literatures have directly addressed to alleviate strong clutter interference [9, 11-13], where the approaches can be classified into two groups. One is pre-filtering before CS recovery that is mathematically represented as a projection to the clutter free subspace [11, 12], and the other is weighting the clutter support combined with CS optimization [9, 13], which is actually a sparsity enhancing technology [14]. Both of the methods have a common point that the clutter spectrum needs to be estimated as a prior knowledge for clutter cancellation. Therefore, the key requirement is the accurate estimation of power spectral density (PSD) of clutter.

Traditionally, statistical spectrum estimation methods such as autoregressive model [9] need sufficient statistically independent and identically distributed (i.i.d.) training data to obtain an effective estimation of the clutter. However, if the clutter scenario is not homogeneous, the range stationarity is destroyed, which results in an inaccurate estimation of the clutter distribution [15]. Referring to the direct data domain (D3) approach designed for space-time adaptive processor (STAP) [15-17], we exploit the priori sparsity of the clutter spectrum and propose a new approach to obtain the high-resolution spectrum with only the test cell.

In this paper (which expands upon [18]), the recently developed CS algorithm based on iterative grid optimization (CS-IGO) with slight grid mismatch is applied to recover the Doppler spectrum. Since the sensing matrix and spectrum obtained by CS-IGO are well adapted to the data, there is almost no clutter sidelobe in the spectrum and then the clutter echo can be obtained by using prior knowledge of its boundary. This guarantees the availability of clutter covariance matrix (CCM) and corresponding adaptive prewhitening filter, which can effectively cancel the mainlobe clutter. Moreover, referring to the sparsity enhancing technology, the iterative reweighted [l.sub.1] minimization is introduced to suppress the clutter residue and boost the target response. Simulation results show that the approach offers a high performance of output signal to clutter and noise ratio (SCNR) in the nonhomogeneous clutter scenario.

The remainder of this paper is organized as follows. In the next section, we mainly introduce the echo signal model of random PRI radar and exploit the sparsity of spectrum for estimation. Section 3 constructs the clutter suppression sketch by the pre-filter and iterative reweighted [l.sub.1] minimization. Section 4 makes some complementary discussions. In Section 5 simulations and numerical illustrations are used to test the performance of clutter estimation and suppression. Some analyses are given as well. Finally, the conclusions are drawn in the last section.

2. DATA MODEL AND SPECTRUM ESTIMATION WITH CS-IGO

2.1. Data Model

In random PRI radar we transmit M pulses at [T.sub.m], which is i.i.d. within the CPI (0, [MT.sub.r]) randomly where [T.sub.r] is the average PRI. Then the transmitted signal can be represented as

s(t) = [M - 1.summation over (m=0)]A rect (t - T/2 - [T.sub.m]/T)exp([j2[pi]f.sub.0]t) (1)

where rect [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], m = 0, 1, ..., M - 1, A is the pulse amplitude, [f.sub.0] the carrier frequency, and T the pulse width.

Suppose that there are K targets in the test range cell moving toward the radar with radial velocity [v.sub.k]. The received signal is generated as the sum of the Doppler-shifted replicas of the transmitted waveform, each of which is multiplied by the respective scattering coefficient. Sampling the received signal of mth pulse after down-conversion and low-pass filtering can obtain

[s.sub.r][m] = [K.summation over (k=1)][A.sub.k] exp([j2[pi]u.sub.k][T.sub.m) (2)

where [A.sub.k] is the complex amplitude proportional to the square root of target radar cross section (RCS), and [v.sub.k] = [2f.sub.0][u.sub.k]/c is the Doppler frequency of the k-th target. Obviously, (2) can be regarded as a typical underdetermined linear system

[S.sub.r] = [THETA](xi)[alpha] + w (3)

when we consider N (> M) Doppler 'grids' at [xi][n][[delta].sub.u], where [s.sub.r] = [[s.sub.r][0], [s.sub.r][1], ..., [s.sub.r][[M - 1]].sup.T], [xi] = [[xi][0], [xi][1], ..., [xi] [[N - 1]].sup.T], [THETA]([xi]) is an M x N variable random sensing matrix, whose elements are [THETA][m, n] = exp(j2n[pi][xi][n][[delta].sub.u][T.sub.m]), the perturbation w = c + n represents clutter c plus noise n, [alpha] is the 'scene' vector with [[alpha].sub.n] = [A.sub.k] if target k locates at [xi][n][[delta].sub.u] and zero else, n = 0, ..., N - 1, [[delta].sup.u] = [l/MT.sub.r] and

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (4)

In this paper, we further assume that the clutter and noise are all Gaussian distributed and mutually independent. Clutter c is zero-mean with covariance matrix [R.sub.c] and noise n is zero-mean with covariance matrix [[sigma].sup.2][I.sub.M] ([I.sub.M] is the identity matrix of size M).

If the signal-to-clutter ratio (SCR) and signal-to-noise ratio (SNR) are sufficiently high, the random matrix [THETA] can be regarded to satisfy the RIP approximately. Unfortunately, the land clutter is always much stronger than the targets, which can violate the assumptions for CS severely [19]. The approaches of clutter cancellation by pre-processing [11] and weighting [14] are two practical ways to preserve the favorable properties of CS provided that the clutter spectrum is estimated as a prior knowledge. Therefore, the key requirement is the accurate estimation of clutter PSD.

2.2. Spectrum Estimation with CS-IGO

Traditionally, autoregressive model together with Burg's algorithm can obtain an effective estimation of the clutter by using sufficient i.i.d. training data [9]. However, it doesn't work well if the range stationarity is destroyed. So we propose a new approach similar to D3 in STAP [15-17] with only the test cell by exploiting the priori sparsity of the clutter spectrum.

[FIGURE 1 OMITTED]

Assume there are only mainlobe clutters in the echo which are caused by stable objects like mountains and buildings [10]. As shown in Fig. 1, it is well known that the Doppler frequency component of mainlobe clutter is approximately limited within

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (5)

where [u.sub.pf] is the platform velocity, [lambda] the radar wavelength, [[theta].sub.a] azimuth beamwidth, and [??] the angle between the velocity vector and the beam axis. In order to enlarge the non-clutter area in Doppler domain, the number of Doppler cells can be increased as large as the matrix RIP holds owe to the unambiguity, and also narrow beamwidth is always formed which leads to small number cells of clutter compared with the whole Doppler cells. Thus, the significant elements only exist in the area of the mainlobe as well as several discrete target positions, which guarantees sparsity of the spectrum.

Assume the clutter c in the test range cell can also be represented by S statistical independent scatter patches with complex amplitudes [C.sub.s] and corresponding Doppler frequencies [u.sub.s], that is

c[m] = [S.summation over (s=1)] [C.sub.s] exp([j2[pi]u.sub.s][T.sub.m]) (6)

where [u.sub.s] = [2.sub.upf] cos([??] + [[theta].sub.s])/[lambda] [member of] [[DELTA].sub.u] with [[theta].sub.s] as the corn angle between the beam axis and the look direction, then the clutter model can be represented as

c = [[THETA].sub.c]([[xi].sub.c])[[alpha].sub.c] (7)

where [[THETA].sub.c]([[xi].sub.c]) and [[alpha].sub.c] represent the sensing matrix and 'scene' vector of the clutter, which are similar to that of the targets. Finally, the signal model can be modified as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (8)

where [THETA]'([xi]') and [alpha]' correspond to the new Doppler 'grids' [[xi]'[n][delta].sub.u] and complex amplitudes containing both targets and clutter.

If only we knew the exact [xi]' we could directly apply the [l.sub.1] minimization method to obtain the parameters by solving the following convex optimization problem

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (9)

where [epsilon] [approximately equal to] [square root of M[[sigma].sup.2]] is the noise level determined by the noise variance [[sigma].sup.2], which plays a significant role on the recovery performance of (9). For example, if it is set to be too small or too large, in the reconstruction, either a significant part of the strong noise is treated as signal or the acquired signal components are not fully reconstructed, resulting in distortion due to too many false components or too few strong true components. In our situation, the noise level can be almost exactly estimated by using the nearby range cells with no clutter following the method proposed in [13].

In practice, however, we always have no priori information of the perfect Doppler 'grids', which may lead to unacceptable results with direct CS method suffering from serious mismatch effect. Therefore, in our situation the CS-IGO algorithm is used by deeming the clutter as signals with high SNR. By updating the 'grids' with estimated mismatch error iteratively to achieve better signal model fit, the parameter estimation mean-squared error is close to the Cramer-Rao lower bound (CRLB). And then the final sensing matrix [??]'([xi]') and the spectrum estimation [??]' with clutters and targets can be well obtained, which is summarized in Table 1.

3. CLUTTER SUPPRESSION WITH PRE-FILTER AND REWEIGHTED [L.sub.1] MINIMIZATION

Because some of the target energy may be losing during the procedure of linear optimization, further adaptive processing for clutter suppression is necessary for subsequent target detection. Considering that the mainlobe clutter is localized and strong in the spectrum whose boundary can be determined by the prior knowledge, we improve the optimization problem from two different aspects: The first is to utilize a pre-filter to cancel the mainlobe clutter, and the second is to construct a reweighted optimization penalty to suppress the residue. As we deem all the components within the boundary as clutter, further assumption should be made that the spectrum of clutter does not overlap with that of targets' spectrum.

3.1. Suppression of Mainlobe Clutter with Pre-filter

In classical processing, cancellation for mainlobe clutter is often performed simply by designing a notch filter, which is mathematically represented as a projection to the clutter free subspace. The same procedure can be performed for the CS processing. Generally, the pre-filter can be represented by an M x M nonsingular matrix F. The sensing matrix is modified to F[THETA]', and the filtered perturbation Fw contains only noise if the clutter contribution is totally removed. Then the target spectrum [[alpha].sub.T] can be obtained by solving the following transformed optimization problem:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (10)

where [eta] = [[parallel]Fw[parallel].sub.2] is the noise level after projection, which can also be estimated by the method proposed in [13] if F is obtained. As [THETA]'([xi]') has already been obtained in (9), the only unknown part of (10) is the pre-filter matrix F. Therefore, we will focus on the problem of how to form F as follows.

In our specific case of airborne systems, because the location of mainlobe clutter in the estimated Doppler spectrum [alpha]' is determined by the platform velocity, the radar wavelength, the look direction and the beamwidth as in (5), all values of which are a prior knowledge when the radar works [20]. Therefore the boundary of clutter can be obtained approximately. Then the estimated clutter can be effectively obtained as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (11)

where [SIGMA] represents a diagonal matrix with elements ones corresponding to the locations with clutter on the diagonal and zeros else. Then the estimated CCM is given by

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (12)

which is proved to be an positive semi-definite Hermite matrix where (x)H denotes conjugate transpose.

Proof: The CCM is obviously an Hermite matrix by testifying [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] . In order to prove its property of positive semi-definite, we rewriter the CCM as the matrix form

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (13)

where [[??].sub.p] denotes the pth element in C, (x)* and denotes conjugate. For any given vector z = [[[z.sub.0],..., [z.sub.M-1]].sup.T] [member of] [C.sup.Mx1], we have

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (14)

Therefore, the CCM is positive semi-definite according to the definition, thus proves [R.sub.c] the positive semi-definite Hermite matrix.

Notate [[sigma].sup.2] as the estimated noise variance and then the estimated covariance matrix of the perturbation w can be estimated as

[R.sub.w] = [R.sub.c] + [[sigma].sup.2][I.sub.M] (15)

which is obviously positive definite and can be expressed by Cholesky factorization [[??].sub.w] = [D.sup.H]D with a nonsingular matrix D. Finally the adaptive prewhitening filter can be derived as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (16)

In order to see the prewhitening effect, we can testify the covariance matrix of w' = Fw as follows

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (17)

which leads to the noise level [eta] = [[parallel]Fw[paralle].sup.2] [approximately equal to] [square root of M].

3.2. Further Suppression of Clutter Residue by Reweighted [L.sub.i] Minimization

By the pre-filter, the strong clutters are almost suppressed without doing harm to the reconstruction process. However, as the prewhitening filter is not perfect, we construct a weighted optimization penalty in (10) to suppress the residue. In terms of forming target spectrum with high quality, one usually needs extraction of the strongest signal components and rejection of others such as clutter/noise residue. Considering that the discrimination between targets and clutter/noise in the spectrum

corresponds to their amplitude difference, we make a distinction in the penalty function. Instinctively, a weighted formulation of [l.sub.1] minimization can be designed as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (18)

where W is the diagonal matrix with the corresponding weight [w.sub.i] for the ith component of on the diagonal and zeros elsewhere. As illuminated in [14], small weights could be used to encourage the recovery of strong components, while the large weights used potentially suppress the small ones. In our special case, as the target spectrum [alpha] with small clutter residue has been obtained by (10), we consider to assign weights directly according to the spectrum. For example, the weight for the ith coefficient of [alpha] can be calculated as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (19)

where [alpha](i) denotes the ith component of [alpha]. Because the spectrum [alpha] is recovered by [l.sub.1] minimization, it is always sparse with many zero elements. To ensures the existence of [W.sup.-1], we introduce a small constant [??] which can be preset as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (20)

Then, these weights can be applied in (21) to encourage the signal support and suppress the residual clutter support.

To solve the problem more efficiently, a modification of (18) is given by

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (21)

and we have [[??]'.sub.T] = [W.sup.-1]x. Herein, the weights are applied to the matrix but not to the coefficients to estimate. Thus, we can also directly use [l.sub.1] minimization method.

Moreover, in order to achieve successively better estimation of the nonzero coefficient locations, an iterative reweighted [l.sub.1] minimization algorithm [14] is used to construct the weights, which is shown in Table 2.

4. COMPLEMENTARY ISSUES AND DISCUSSIONS

After a series of processing, the targets' velocities and amplitudes can be finally obtained with considerable accuracy. However, it should be noted that there are several issues to be further explained.

Firstly, as the number of statistical independent clutter scatters distributed within the boundary may be more than that of the clutter Doppler cells, what we recovered by CS-IGO is not the exact clutter characteristics but the main components of the coherent summation of all the scatters. Fortunately, due to the nearly orthogonal columns of the sensing matrix, the mainlobe clutter can be well represented by the Doppler cells within the boundary because the signal can be approximated sufficiently by a linear combination of the neighboring grids [11], and more precisely with out optimized grids. Moreover, when we transform the clutter spectrum back into the time domain, most of the clutter components and energy can be reserved, which hardly affects the design of pre-filter.

Secondly, apart from the [l.sub.1] minimization there are two additional iteration techniques in our approach, whose convergence properties are still not characterized completely. Because the sensing matrix is only approximately known at any intermediate iteration and hence the traditional proof techniques do not apply. In principle, given access to a sufficient number of measurements and high SNR, we may expect similar convergence behavior for the two algorithms as conventional CS algorithm. Empirically, they can be shown to be stable to small amounts of noise in the sparse signal, which can be demonstrated with the help of various numerical experiments in Section 5.

Finally, the approach's computational complexity is another concerned issue. It can be seen that the primary runtime cost is incurred in solving CS-IGO and iterative reweighted [l.sub.i] minimization, which are both determined by the performance of solving a linear system and the iteration number respectively. The overall runtime of the approach is relatively larger than traditional methods, which requires further investigation on fast [l.sub.i] minimization. However, these fast methods are not covered in this paper and our attention is mainly focused on sound clutter suppression with only the test range cell.

5. NUMERICAL SIMULATIONS

In this section, we conduct some simulations to demonstrate the effectiveness and the feasibility of the proposed method. Assume that the radar works on X-band whose wavelength is 0.03 m, the coherent processing interval is 6.4 ms and the width of each pulse in the train is 2 (is. We transmit M = 64 pulses and consider N = 128 Doppler cells. What's more, there are some preconditions to be demonstrated.

1) For the random PRI signal, the pulses are transmitted at time deviations which are selected in the coherent processing interval with discrete uniform distribution.

2) In the following simulations the SNR, SCR or clutter-to-noise ratio (CNR) represents for the sampling echo signal after down-conversion and low-pass filtering, which is added by the complex white Gaussian noise.

3) For the clutter scenario, we directly assume that in the 3140th Doppler cells there are 50 independent clutter scatters, each of which has a Gaussian distributed complex amplitude and a uniform distributed Doppler frequency.

4) All the convex optimization problems are solved by employing cvx [21] as an effective tool.

5) Our simulations are performed in MATLAB7 environment using a Pentium (R) 4 CPU 3.00 GHz processor with 1 GB of memory, and under Microsoft Windows XP operating system.

[FIGURE 2 OMITTED]

5.1. Clutter Spectrum Estimation by CS and CS-IGO

In this experiment, to characterize the clutter spectrum estimation performance of the proposed approach, we add white Gaussian noise to the echo signal with only clutter and no target, and both direct CS and CS-IGO approaches are used for clutter recovery.

Firstly, we perform several experimental trials under CNR = 10 dB and the typical estimated clutter spectrums with direct CS and CS-IGO approaches are shown in Fig. 2. It indicates that the 50 independent clutter scatters can be well represented by the 10 Doppler cells, which are regarded as main components of the coherent summation of all the scatters. What's more, the clutter spectrum with CS-IGO approach shows smaller false components than that of direct CS method, thus maintaining much more energy when we transform the clutter spectrum back into the time domain.

Secondly, in order to see the influence of estimated noise level to the convex optimization problem, we set [epsilon] = [kappa] [square root of M[[sigma].sup.2] with [kappa] = 1,1.2,1.4,1.6 respectively. Defining the clutter relative mean square errors relative MSE (Relative MSE) as

Relative MSE = [[parallel]c - c[parallel].sup.2.sub.2]/[[parallel]c[parallel].sup.2.sub.3] (22)

the distributions of relative MSEs under CNR=10dB are shown in Fig. 3 with 100 Monte Carlo simulations. It can be seen from the histograms that the relative MSEs of clutters recovered by the CS-IGO approach are much smaller, which means more precise estimation. And we can observe another phenomenon that in both direct CS and CS-IGO, as the estimated noise level deviates form the true value, the estimation error of clutter increases obviously, which reflects the significance of exact estimation of noise level.

Thirdly, we will testify how noise influences the approach performance. The relative MSE results of estimated clutters are calculated by averaging 50 Monte Carlo realizations under each CNR, which increases from 0 to 20dB. As shown in Fig. 4, it is obvious that the performance meliorates with the increase of CNR. Especially, when CNR is above 10 dB the estimated clutter with CS-IGO changes slowly and the precision is always acceptable in practice, however, when CNR is below 10 dB the relative MSE increases sharply with the decrease of CNR, which is due to the sensitivity of CS theory to low CNR.

[FIGURE 3 OMITTED]

[FIGURE 4 OMITTED]

5.2. Comparison of Whitening Algorithms

In our first experiment, we will analyze the realization of weighting algorithm for whitening proposed in [9], the basic idea of which is to solve the following minimization problem

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (23)

where [[??].sub.c] is the estimated clutter spectrum, [zeta] a constant scale factor, and [??] the Hadamard product. It is indicated in [9] that by the clutter whiten weighting clutter is suppressed to the noise level and [[alpha].sub."T] only denotes the unambiguous Doppler spectrum of targets. In our particular situation, however, as the clutter spectrum [[??].sub.c] = [SIGMA][??]' is a sparse vector, those zero elements in the vector should be replaced with a small constant to guarantee the solvability of (23), and then the clutter spectrum for weighting can be represented as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (24)

where [[??].sub.c](i) denotes the ith component of [[??].sub.c]. In practice, weighting is actually a sparse enhancement technique and weighting only the clutter supports is a locally sparse enhancement technique. Therefore, on the one hand, when the constant [zeta] is set too large, the clutter will not be so well whitened that large mainlobe clutter residue exists. On the other hand, when the constant [zeta] is set too small, the signal energy will weaken and even the joint optimization problem becomes irresolvable so that the target spectrum cannot be recovered.

Assume that there are three targets in the test range cell located randomly among the non-clutter Doppler area and SCR = -10dB, SNR = 10 dB, we present the simulation results of a typical experiment. The original spectrum with CS-IGO and the target spectrums after weighting with [zeta] = 1,[10.sup.-2],[10.sup.-4] are shown in Fig. 5 respectively. It can be seen that the weighting algorithm can well suppress the mainlobe clutter with proper value of constant [zeta] (e.g., [zeta] = [10.sup.-2],[10.sup.-4]) and the clutter residue will be decreased with the reduction of (. However, the signal energy will decrease simultaneously, and what's more, when the constant is set to [zeta] = [10.sup.-6], the optimization problem (23) becomes irresolvable and the target spectrum cannot be recovered. Therefore, it can be concluded that the value of [zeta] is crucial to the performance of weighting algorithm and it is hard to decide the value of [zeta] in a time varying clutter environment.

[FIGURE 5 OMITTED]

In the next experiment, we will compare the pre-filter whitening algorithm with the weighting algorithm. For justice, we consider the same scenario as in last subsection. By directly following the algorithm steps, the CCM and corresponding adaptive filter can be formed with the estimated clutter spectrum and noise variance. Then the target spectrum after pre-whitening is shown in Fig. 6(a), from which it can be seen that mainlobe clutter can also be well suppressed with our method. Further comparison to the spectrum by weighting algorithm (with [zeta] = [10.sup.-4]) as shown in Fig. 6(b) indicates that, using our approach, the mainlobe clutter residue is more trivial and there are less energy losses of the targets. Most importantly, there is no parameter needed to be manually setting in our algorithm, which is more robust in practice. Although some sidelobe clutter residues appear in the spectrum, they can be suppressed by the iterative reweighted [l.sub.1] minimization algorithm, which is presented in Fig. 6(c). And finally the targets can be recovered with almost no clutter residue, which provides much high precision of amplitude and velocity estimation of the targets.

In order to compare the computational performance of the algorithms, we present the consumed time and memory cost in Table 3 for data size of M = 64 pulses and N = 128 Doppler cells. It can be seen from the table that pre-filtering is a litter more time-consuming and needs a larger memory cost than weighting due to calculation and storage of the prefilter. For the pre-filtering and reweighted [l.sub.1] minimization method, as there are additional convex optimization problems, more time is consumed, however, the memory cost remains unchangeable. Although the cost time seems relatively much larger than traditional methods, some fast algorithms for the [l.sub.1]-minimization problems have already been proposed [22-25], which is not covered in this paper.

[FIGURE 6 OMITTED]

5.3. Performance of Clutter Suppression

In order to characterize the clutter suppression performance of the proposed algorithm more quantitatively, the output SCNR defined in (25) and the normalized square errors (NSE) of the target amplitude modulus defined in (26) are calculated for a single target under each input SNR and SCR.

SCNR = [10log.sub.10][[A.sup.2]/([[parallel][[alpha].sub.T][parallel].sup.2.sub.2] - [A.sup.2])] (25)

NSE = [(A - A).sup.2]/[A.sup.2] (26)

where A is the estimated amplitude of the target, A the actual amplitude of the target, and [[??].sub.T] the recovered target spectrum. For comparison, the following four different clutter suppression schemes are used: a) weighting with [zeta] = [10.sup.-2]; b) weighting with [zeta] = [10.sup.-4]; c) pre-filtering; d) proposed algorithm by pre-filtering and reweighted [l.sub.1] minimization.

[FIGURE 7 OMITTED]

Firstly, by setting SCR = -10 dB, the output SCNR and NSE results are averaged over 50 independent Monte Carlo realizations under each SNR, which increases from 0 to 20 dB. It can be seen in Fig. 7 that, for the pre-filtering algorithm, although its output SCNR is a litter lower than the weighting algorithm, its amplitude precision is higher. Moreover, for the proposed algorithm by pre-filtering and reweighted l1 minimization, its output SCNR and amplitude precision are both higher, which is more effective than merely weighting the clutter supports. Furthermore, as shown in Fig. 7(a), when the SNR is not very high, the targets can be recovered with almost no clutter residue by the proposed algorithm which leads to high output SCNR. However, when the SNR is higher than 14 dB, the output SCNR performance of the proposed algorithm is almost the same as that of weighting with [zeta] = [10.sup.-4]. This is because that when the noise level is too low, the energy of sidelobe clutter residues after pre-filtering is comparable to the noise, which leads to small improvement of output SCNR by reweighted [l.sub.1] minimization and there is still some sidelobe clutter residues in the final spectrum.

Secondly, by setting SNR= 10 dB, the output SCNR and NSE results are averaged over 50 independent Monte Carlo runs under each SCR, which decreases from 0 to -20 dB. As shown in Fig. 8, it can also be seen that under high SCRs (e.g. higher than -14 dB), the proposed clutter suppression scheme is much more effective than merely weighting the clutter supports. Similarly, it can be seen form Fig. 8(a) that when the SCR is lower than -14 dB, the output SCNR of the proposed algorithm is sometimes lower than that of the weighting algorithm, which is also caused by the large energy of sidelobe clutter residues after pre-filtering.

[FIGURE 8 OMITTED]

Finally, it can be concluded that the proposed algorithm performs well in amplitude estimation at all times and can improve the SCNR dramatically under proper situations, which is much more effective than merely weighting the clutter supports.

6. CONCLUSION

In this paper, an adaptive clutter suppression method is proposed for airborne random pulse repetition interval radar. After recovering the clutter spectrum with only the test range cell by exploiting its intrinsic sparsity, the optimization problem is improved by utilizing a pre-filter and an iterative weighted penalty, offering a high performance of amplitude precision and output SCNR in the nonhomogeneous clutter scenario. Simulation results demonstrate that the CS-IGO approach is more precise than direct CS method in estimating the clutter spectrum and the proposed clutter suppression scheme is also more effective than merely weighting the clutter supports. How to decrease the overall runtime of the approach and how to extend the approach to wideband imaging radar are underway.

Received 20 February 2012, Accepted 21 May 2012, Scheduled 1 June 2012

ACKNOWLEDGMENT

This work is supported by the innovation project for excellent postgraduates of National University of Defense Technology under Grant B110404, and Hunan Provincial Innovation Foundation for Postgraduate under Grant CX2011B019. The authors would also like to thank the Editor and the anonymous reviewers for their helpful comments and suggestions to improve the quality of this paper.

REFERENCES

[1.] Kaveh, M. and G. R. Cooper, "Average ambiguity function for a randomly staggered pulse sequence," IEEE Trans. Aerosp. Electron. Syst., Vol. 12, No. 3, 410-413, May 1976.

[2.] Vergara-Dominguez, L., "Analysis of the digital MTI filter with random PRI," IEE Proceedings-F, Vol. 140, No. 2, 129-137, Apr. 1993.

[3.] Cook, C. E. and M. Bernfeld, Radar Signals: An Introduction to Theory and Application, Academic Press, New York, 1967.

[4.] Donoho, D., "Compressed sensing," IEEE Trans. Inf. Theory, Vol. 52, No. 4, 1289-1306, Apr. 2006.

[5.] Candes, E. and T. Tao, "Near optimal signal recovery from random projections: Universal encoding strategies?" IEEE Trans. Inf. Theory, Vol. 52, No. 12, 5406-5425, Dec. 2006.

[6.] Candes, E., J. Romberg, and T. Tao, "Robust uncertainty principles: Exact signal reconstruction from highly incomplete frequency information," IEEE Trans. Inf. Theory, Vol. 52, No. 2, 489-509, Feb. 2006.

[7.] Wei, S. J., X. L. Zhang, J. Shi, and G. Xiang, "Sparse reconstruction for SAR imaging based on compressed sensing," Progress In Electromagnetics Research, Vol. 109, 63-81, 2010.

[8.] Wei, S. J., X. L. Zhang, and J. Shi, "Linear array SAR imaging via compressed sensing," Progress In Electromagnetics Research, Vol. 117, 299-319, Jun. 2011.

[9.] Quan, Y. H., L. Zhang, M. D. Xing, and Z. Bao, "Velocity ambiguity resolving for moving target indication by compressed sensing," Electronics Letters, Vol. 47, No. 22, Oct. 2011.

[10.] Barton, D. K. and S. A. Leonov, Radar Technology Encyclopedia, Artech House, Boston, London, 1998.

[11.] Ender, J. H. G., "On compressive sensing applied to radar," Signal Processing, No. 90, 1402-1414, 2010.

[12.] Khwaja, A. S. and J. Ma, "Applications of compressed sensing for SAR moving-target velocity estimation and image compression," IEEE Transactions on Instrumentation and Measurement, Vol. 60, No. 8, 2848-2860, 2011.

[13.] Zhang, L., M. Xing, C.-W. Qiu, J. Li, J. Sheng, Y. Li, etal., "Resolution enhancement for inversed synthetic aperture radar imaging under low SNR via improved compressive sensing," IEEE Trans. Geosci. Remote Sens., Vol. 48, No. 10, 3824-3838, Oct. 2010.

[14.] Candes, E., M. Wakin, and S. Boyd, "Enhancing sparsity by reweighted [l.sub.1] minimization," J. Fourier Anal. Appl., Vol. 14, No. 5, 877-905, Dec. 2008.

[15.] Sun, K., H. Meng, Y. Wang, and X. Wang, "Direct data domain STAP using sparse representation of clutter spectrum," Signal Processing, Vol. 91, No. 9, 2222-2236, 2011.

[16.] Choi, W., T. K. Sarkar, W. Hong, and E. L. Mokole, "Adaptive processing using real weights based on a direct data domain least squares approach," IEEE Transactions on Antennas and Propagation, Vol. 54, No. 1, 182-191, 2006.

[17.] Burintramart, S., T. K. Sarkar, Y. Zhang, and M. C. Wicks, "Performance comparison between statistical-based and direct data domain STAPs," Digital Signal Processing, Vol. 17, 737-755, 2007.

[18.] Liu, Z., H. Wang, Y. Qin, and X. Li, "Adaptive clutter suppression for airborne random PRI radar based on improved compressed sensing," Proc. CoSeRa2012, 2012.

[19.] Levanon, N. and E. Mozeson, Radar Signals, Wiley, New York, 2004.

[20.] Zhang, M. and X. Wang, Radar Systems, Publishing House of Electronics Industry, Beijing, 2006.

[21.] Grant, M. and S. Boyd, "CVX: Matlab software for disciplined convex programming," http://stanford.edu/_boyd/cvxCVX, 2008.

[22.] Tropp, J. A. and S. J. Wright, "Computational methods for sparse solution of linear inverse problems," Proceedings of the IEEE, Vol. 98, No. 6, 948-958, Jun. 2010.

[23.] Yang, A. Y., A. Ganesh, Z. Zhou, S. S. Sastry, and Y. Ma,

"A review of fast li-minimization algorithms for robust face recognition," http://arXiv:1007.3753v2 [cs.CV] 29 Jul 2010, 2010.

[24.] Mailhe, B., R. Gribonval, P. Vandergheynst, and F. Bimbot, "Fast orthogonal sparse approximation algorithms over local dictionaries," Signal Processing, doi:10.1016/j.sigpro.2011.01.004, 2011.

[25.] Donoho, D. L. and Y. Tsaig, "Fast solution of li-norm minimization problems when the solution may be sparse," IEEE Trans. Inf. Theory, Vol. 54, No. 11, 4789-4811, Nov. 2008.

Z. Liu *, X. Z. Wei, and X. Li

Research Institute of Space Electronics Information Technology, School of Electronic Science and Engineering, National University of Defense Technology, Changsha 410073, China

* Corresponding author: Zhen Liu (zhen_liu@nudt.edu.cn).

Table 1. The iterative grid optimization algorithm. Input: The signal sampling vector [s.sub.r] [member of] [R.sup.M], the random PRI [T.sub.m], the estimated noise bound [epsilon], the error threshold [??], the maximum number of iterations [L.sub.max], the Doppler resolution [[delta].sub.u] = [1/MT.sub.r], the number of Doppler cells N > M, initialization [[sigma].sup.(0).sub.[xi]] = [[sigma].sup.(0).sub.[n]] 1, l = 0, the elements of [[xi].sup.'(loop)] are set as [[xi].sup.'(loop)][n] = n, loop = 0, ..., [L.sub.max] - 1, n = 0, ..., N - 1. Output: The final sensing matrix [THETA]'([xi]') and the spectrum estimation [alpha]' with clutters and targets. Step1: Compute [[THETA].sup.'(0)] from [[xi].sup.'(0)] and obtain the coarse spectrum [[alpha].sup.'(0)] by using (9). Step2: Set l = l + 1 and estimate the support set vector [[LAMBDA].sup.(l)] of [[alpha].sup.'(l] - 1) by constant false alarm rate (CFAR) detection to form [K.sup.(l)], [[alpha].sup.'(l).sub.[LAMBDA]] and [[THETA].sup.[l].sub.[LAMBDA]], where [K.sup.(l)]\ is the vector length of [[LAMBDA].sup.(l)], [[alpha]'.sup.(l).sub.[LAMBDA]], [[THETA].sup.(l).sub.[LAMBDA] denote the elements and columns in [[alpha]'.sup.(l - 1)], [[THETA]'.sup.(0)] indexed within [[LAMBDA].sup.(l)]. Step3: Compute [A.sup.(l).sub.a]) as well as [A.sup.+(l).sub.a] by using [[LAMBDA].sup.(l)], [[alpha]'.sup.(l).sub.[LAMBDA]], [[sigma].sup.(l-1).sub.[xi]] and [[sigma].sup.(l - 1).sub.n] according to [A.sub.a] = [[sigma][xi][PHI][xi], [[sigma].sub.n][I.sub.M]], where [A.sup.+.sub.a] is the pseudo-inverse of Aa, [[PHI].sub.[beta]][m,k] = j2[pi][T.sub.m] exp(j2[pi][LAMBDA][k] [[delta].sub.u][T.sub.m])[a'.sub.[LAMBDA]][k], [a'.sub.[LAMBDA]][k] denotes the kth coefficient in [a'.sup.(l).sub.[LAMBDA]] and k = 0, ..., K(l) - 1. Step4: Update [a'.sup.(l).[LAMBDA]] according to [a.sub.[LAMBDA]] = [([H.sup.H]H).sup.-1][H.sup.H]y where H = [A+.sup.(l).sub.a][[THETA]].sup.(l).sub.[LAMBDA]] and y = [[??].sup.+(l).sub.a][s.sub.r]. Step5: [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] Step6: [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] and go to Step 2. Table 2. The iterative reweighted L1 minimization algorithm. Input: The obtained signal [F.sub.sr] and sensing matrix [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], the estimated noise bound [eta], the error threshold [??]', the maximum number of iterations [L'.sub.max], initialization l' = 0, the estimated target spectrum with clutter residue [[??].sup.(0).sub.T] = [[??].sub.T] obtained by (10), the weight matrix [W.sup.(0)] calculated by (19) and (20). Output: The final estimated target spectrum [[??].sub.TF]. Step1: Solve the weighted [l.sub.1] minimization problem [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. Step2: [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], stop; else update the weights by using [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] set l0 = l0 + 1 and go to Step 1. Table 3. Performance and cost time of different algorithms. Algorithms Consumed Memory Cost Time (s) (Double Complex) Weighting 1.578 8384 Pre-filtering 1.659 12480 Pre-filtering and reweighted [l.sub.1] 3.024 12480 minimization

Printer friendly Cite/link Email Feedback | |

Author: | Liu, Z.; Wei, X.Z.; Li, X. |
---|---|

Publication: | Progress In Electromagnetics Research |

Article Type: | Report |

Geographic Code: | 9CHIN |

Date: | Jun 1, 2012 |

Words: | 6350 |

Previous Article: | Multilayer model formulation and analysis of radar backscattering from sea ice. |

Next Article: | Printed wideband antenna with chip-capacitor-loaded inductive strip for LTE/GSM/UMTS WWAN wireless USB dongle applications. |

Topics: |