Next Article in Journal
Large Aerial Image Tie Point Matching in Real and Difficult Survey Areas via Deep Learning Method
Next Article in Special Issue
On the Efficient Implementation of Sparse Bayesian Learning-Based STAP Algorithms
Previous Article in Journal
MFCD-Net: Cross Attention Based Multimodal Fusion Network for DPC Imagery Cloud Detection
Previous Article in Special Issue
Transmit Antenna Selection and Power Allocation for Joint Multi-Target Localization and Discrimination in MIMO Radar with Distributed Antennas under Deception Jamming
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Sparse Bayesian Space-Time Adaptive Processing Algorithm to Mitigate Off-Grid Effects

National Laboratory of Radar Signal Processing, Xidian University, Xi’an 710071, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(16), 3906; https://doi.org/10.3390/rs14163906
Submission received: 14 July 2022 / Revised: 3 August 2022 / Accepted: 9 August 2022 / Published: 11 August 2022
(This article belongs to the Special Issue Small or Moving Target Detection with Advanced Radar System)

Abstract

:
Space-time adaptive processing (STAP) algorithms based on sparse recovery (SR) have been researched because of their low requirement of training snapshots. However, once some portion of clutter is not located on the grids, i.e., off-grid problems, the performances of most SR-STAP algorithms degrade significantly. Reducing the grid interval can mitigate off-grid effects, but brings strong column coherence of the dictionary, heavy computational load, and heavy storage load. A sparse Bayesian learning approach is proposed to mitigate the off-grid effects in the paper. The algorithm employs an efficient sequential addition and deletion of dictionary atoms to estimate the clutter subspace, which means that strong column coherence has no effect on the performance of the proposed algorithm. Besides, the proposed algorithm does not require much computational load and storage load. Off-grid effects can be mitigated with the proposed algorithm when the grid-interval is sufficiently small. The excellent performance of the novel algorithm is demonstrated on the simulated data.

1. Introduction

Space-time adaptive processing (STAP) exhibits great potential in suppressing clutter and detecting slow-moving targets [1]. The Reed–Mallet–Brennan (RMB) rule [2] states that signal-to-clutter-plus-noise ratio (SCNR) loss can be within 3 dB only if the number of independent and identically distributed (i.i.d) training snapshots is at least twice the system degrees of freedom (DOFs). Unfortunately, due to array configurations and the complex clutter environment, sufficient training snapshots cannot be guaranteed, and therefore the performance of STAP degrades significantly.
Motivated by compressed sensing (CS) theory, STAP algorithms have been enriched over the past 15 years with the application of sparse recovery (SR) techniques. These algorithms are termed SR-STAP algorithms [3,4,5,6,7,8,9]. The clutter-plus-noise covariance matrix (CNCM) can be estimated with SR-STAP algorithms when the training snapshots are limited. Most of these algorithms are exploited to solve an 1 -norm minimization problem instead of an 0 -norm minimization problem [10]. However, large coefficients in the 1 -norm are penalized more heavily, and their performances are greatly affected by the choice of regularization parameters. Compared with the above SR-STAP algorithms, sparse Bayesian learning (SBL) [11,12,13,14,15,16,17,18,19,20,21] has received much attention and has been applied to the STAP framework [22,23,24] because of its robustness and excellent performance. However, the performance of all SR-STAP algorithms relies on the match between clutter and the dictionary matrix.
The dictionary used in SR-STAP algorithms consists of the space-time steering vectors corresponding to all the discrete grids in the angular-Doppler plane. Based on the assumption that the clutter ridge lies exactly on the sampling grids, most SR-STAP algorithms aim to improve the accuracy of signal recovery. However, when this assumption does not hold, some portion of clutter is not located on the grid, i.e., off-grid problems, which leads to inaccuracy of the CNCM estimated by SR-STAP algorithms and performance degradation [25]. Unfortunately, off-grid problems are unavoidable and common in practical applications.
Some algorithms have been proposed to overcome the off-grid problems in direction-of-arrival (DOA) estimation [26,27,28,29] but the off-grid problems in SR-STAP have not received attention. The first-order approximation method, which is common and effective in the off-grid DOA estimation, cannot be directly applied to off-grid SR-STAP because STAP is two-dimensional signal processing. In [25], to mitigate the off-grid effects, the dictionary is constructed by exploiting the knowledge of clutter ridge. Regrettably, precise environmental knowledge is hard to obtain and changes with time. In [30,31], the global atoms are selected from the global STAP dictionary and then the optimal atoms are searched from the local STAP dictionary. However, the convergence of these algorithms is not guaranteed.
It is obvious that an effective way to deal with the off-grid problems is to reduce the grid interval to obtain a denser grid set. However, traditional SR-STAP algorithms suffer from heavy computational load and storage load when utilizing the dense grid set. Moreover, if the grid-interval is small enough, the dictionary matrix cannot satisfy the restricted isometry property (RIP) condition [27], which also results in the inaccuracy of the CNCM and the performance degradation of SR-STAP algorithms.
In this paper, a novel SR-STAP algorithm is proposed to mitigate the effects of off-grid problems with sufficiently dense grids. The specific steps are as follows: Firstly, the angular-Doppler plane is discretized into sufficiently dense grids to construct an overcomplete dictionary; secondly, based on the marginal likelihood maximization criterion, we select the atoms that are in the clutter subspace and reserve these atoms in the new dictionary; thirdly, we estimate the CNCM and calculate the STAP weight vector.
Although sufficiently dense grids lead to strong column coherence of atoms in the dictionary, the core idea is to select the most accurate atoms to represent the clutter subspace based on the marginal likelihood maximization criterion, which indicates that the performance of the proposed cannot be affected by strong column coherence of atoms. Besides, due to the sparsity of clutter, the number of the reserved atoms is small, and the new dictionary is a low-dimensional matrix, which indicates that the proposed algorithm does not require a heavy computational load and storage load.
In conclusion, compared with the traditional SR-STAP algorithm, there are two main advantages of the proposed algorithm: (i) Compared with traditional SR-STAP algorithms, the run time and storage load of the proposed algorithm are both less in the absence and presence of off-grid problems; furthermore, (ii) the proposed algorithm can mitigate off-grid effects with sufficiently dense grids.
Notation: Scalar quantities, vectors, and matrices are denoted by italic typeface, boldface small letters, and boldface capital letters, respectively. The i -th entry of x is denoted by x i . The i -th column and ( i , j ) -th element of A are denoted by a i and A i j , respectively. The expectation is denoted by Ε ( · ) . diag ( · ) denotes a diagonal matrix formed by a vector or diagonal elements of a matrix. The matrix transpose, conjugate transpose, and inverse are denoted by ( · ) T , ( · ) H , and ( · ) 1 , respectively.   ·   F denotes the Frobenius norm and   ·   2 , 0 denotes the 0 -norm of the vector formed by the 2 -norm of each row. |   ·   | denotes the determinant. I denotes the identity matrix.

2. Signal Model

Without a loss of generality, an airborne pulsed-Doppler radar under consideration is equipped with a uniform linear array (ULA). N antenna elements are placed with a half-wavelength inner spacing. M pulses are transmitted at a constant pulse repetition frequency (PRF). The received signal x can be formulated as the superposition of clutter, thermal noise, and might targets. Ignoring the range ambiguity, the clutter x c M N × 1 can be modelled as the superposition of echoes from all the clutter patches in the cell under test (CUT).
x c = i = 1 N c α i υ st ( f d , i , f s , i )
υ st ( f d , i , f s , i ) = υ t ( f d , i ) υ s ( f s , i )
υ t ( f d , i ) = [ 1 , exp ( j 2 π f d , i ) , , exp ( j 2 π ( M 1 ) f d , i ) ] T
υ s ( f s , i ) = [ 1 , exp ( j 2 π f s , i ) , , exp ( j 2 π ( N 1 ) f s , i ) ] T
where N c is the number of clutter patches in the CUT; α i , f d , i , and f s , i are the complex amplitude, the normalized Doppler frequency, and the spatial frequency of the i -th clutter patch, respectively; the space-time steering vector υ st can be obtained by the Kronecker product operation of the temporal steering vector υ t and the spatial steering vector υ s ; denotes the Kronecker product operation.
In the sparsity-based STAP, we need to construct a dictionary matrix to represent the angular-Doppler plane. The angular-Doppler plane is discretized into K = N s N d grids, where N s = φ s N ( φ s > 1 ) and N d = φ d M ( φ d > 1 ) are the number of normalized spatial frequency bins and Doppler frequency bins along the spatial frequency axes and the Doppler frequency axes, respectively. Each grid corresponds to a space-time steering vector, and the dictionary consists of all the space-time steering vectors. Assuming that clutter is located exactly on the grids, the received clutter plus noise snapshots X M N × L from L range cells in the multiple measurement vectors (MMV) case can be expressed as
X = D A + N
where D M N × K , A K × L , and N M N × L denote the dictionary matrix, the sparse coefficient matrix, and the noise matrix, respectively. D can be expressed as
D = [ υ 1 , υ 2 ,   , υ K ]
where
υ k = υ t ( f d , k ) υ s ( f s , k )
In (7), f d , k and f s , k are the normalized Doppler frequency and the spatial frequency of the k -th grid, respectively.
In the SR-STAP algorithms, the noise-contaminated X is required to be denoted with as few atoms as possible.
A = arg min A A 2 , 0 , s . t . X D A F 2 ε
where ε is an error tolerance parameter related to the noise power.

3. Off-Grid Problems

In most SR-STAP algorithms, the grids are uniformly sampled in the angular-Doppler plane along the spatial frequency axes and the Doppler frequency axes, respectively. The performance of SR-STAP algorithms can be achieved on the assumption that the clutter ridge lies on the grids. Once some portion of the clutter ridge is not sampled by the grids, i.e., off-grid problems, the accuracy of the CNCM estimated by SR-STAP algorithms cannot be guaranteed, which leads to significant performance degradation. For example, (i) the slope of the clutter ridge is not equal to the ratio of N s and N d in the side-looking radar case and (ii) the clutter ridge of a non-side-looking radar is a non-linear curve.
Although a reduction in the grid interval can mitigate the effects of the off-grid problems, traditional SR-STAP algorithms suffer from heavy computational load and storage load. Moreover, if the grid- interval is small enough, the dictionary matrix cannot satisfy the restricted isometry property (RIP) condition. In the next section, a novel SR-STAP algorithm is proposed to overcome the aforementioned problems when the grid- interval is reduced.

4. The Proposed Algorithm to Mitigate Off-Grid Effects

4.1. Construction of the Dictionary

The values of φ s and φ d can be set to an integer between 2 and 5 in the absence of off-grid problems. When off-grid problems occur, we need to reduce the grid interval by increasing the values of φ s and φ d in the proposed algorithm. In the presence of off-grid problems, the values of φ s and φ d can be set to an integer more than 10 or even higher.
Although precise environmental knowledge is hard to obtain, we still can know the approximate location of the clutter. To speed up the proposed algorithm, the grids that are definitely not in the approximate location of clutter can be removed from the dictionary at the beginning. However, the dictionary used in traditional SR-STAP needs to be over-completed. For a fair comparison, all SR-STAP algorithms exploit all grids in the whole angular-Doppler plane to construct the dictionary in the simulations.

4.2. Estimation of the Clutter Subspace

Assuming that the noise in (5) is complex Gaussian white noise, the likelihood function of X can be expressed as
p ( X | A , σ 2 ) = ( π σ 2 ) M N L exp [ σ 2 l = 1 L x l D a l 2 ]
where A = [ a 1 , a 2 , , a L ] and the noise variance σ 2 can be calculated with the knowledge of the radar system parameters.
Assuming that a l ,   l are submitted to the same zero-mean complex Gaussian prior distribution, the prior of A can be expressed as
p ( A | Γ ) = π K L | Γ | L exp ( l = 1 L a l H Γ 1 a l )
where Γ = diag ( γ ) represents the prior variance of a l , l and γ = [ γ 1 , γ 2 , , γ K ] T is a vector of hyper-parameters corresponding to all grids. Since Γ represents the prior variance of a l , l , γ i 0 .
Based on the Bayes rule, the posterior density of A can be expressed as
p ( A | X , Γ , σ 2 ) = p ( X | A , σ 2 ) p ( A | Γ ) p ( X | A , σ 2 ) p ( A | Γ ) d A = π K L | Σ | L exp [ l = 1 L ( a l μ l ) H Σ 1 ( a l μ l ) ] = l = 1 L C N ( a l | μ l , Σ )
where
Σ = ( σ 2 D H D + Γ 1 ) 1
μ l = σ 2 Σ D H x l
The next step is to estimate γ . If we adopt the expectation maximization (EM) approach to estimate γ , then we can obtain the SBL algorithm with MMV (M-SBL) in [22]. Another approach is that γ can be point estimated by maximizing the marginal likelihood function. The marginal likelihood function p ( X | Γ , σ 2 ) is expressed as
p ( X | Γ , σ 2 ) = p ( X | A , σ 2 ) p ( A | Γ ) d A = π N M L | D Γ D H + σ 2 I | L exp [ l = 1 L x l H ( D Γ D H + σ 2 I ) 1 x l ]
Define
C D Γ D H + σ 2 I
and γ can be point estimated by
γ = arg max γ ln p ( X | Γ , σ 2 ) = arg max γ ( L ln | C | l = 1 L x l H C 1 x l )
Define
( Γ ) L ln | C | l = 1 L x l H C 1 x l
and (16) can then be expressed as
γ = arg max γ ( Γ )
An effective method to solve Equation (18) is to update a single hyper-parameter γ j at a time. The update of γ j makes the biggest contribution to the maximization of ( Γ ) . Next, we introduce how to select the serial number j .
C in (15) can be rewritten as
C = σ 2 I M N + 1 k K , k i γ k ν k ν k H + γ i ν i ν i H = C i + γ i ν i ν i H
where C i = σ 2 I M N + 1 k K , k i γ k ν k ν k H contains all terms that are independent of γ i .
Using the Woodbury Matrix Identity
C 1 = C i 1 C i 1 ν i ( γ i 1 + ν i H C i 1 ν i ) 1 ν i H C i 1
| C | = | C i | | 1 + γ i ν i H C i 1 ν i |
Equation (17) can then be expressed as
( Γ ) = l = 1 L [ x l H C i 1 x l + ln | C i | ] + l = 1 L [ x l H C i 1 ν i ( γ i 1 + ν i H C i 1 ν i ) 1 ν i H C i 1 x l ln | 1 + γ i ν i H C i 1 ν i | ] = i + ( γ i )
where
i = l = 1 L [ x l H C i 1 x l + ln | C i | ]
( γ i ) = l = 1 L [ x l H C i 1 ν i ( γ i 1 + ν i H C i 1 ν i ) 1 ν i H C i 1 x l ln | 1 + γ i ν i H C i 1 ν i | ]
Equation (22) has been divided into two parts: The part independent of γ i is denoted as i and the other part related to γ i is denoted as ( γ i ) .
Define
s ^ i ν i H C i 1 ν i , q ^ i , l ν i H C i 1 x l
and ( γ i ) can then be simplified to
( γ i ) = l = 1 L [ q ^ i , l H ( γ i 1 + s i ) 1 q ^ i , l ln | 1 + γ i s i | ]
Differentiate ( γ i ) with respect to γ i
( γ i ) γ i = l = 1 L [ q ^ i , l H ( 1 + γ i s i ) 2 q ^ i , l s i ( 1 + γ i s i ) 1 ] = 0
and we can obtain the optimal γ i by
γ i = { 1 L s i 2 l = 1 L ( q ^ i , l H q ^ i , l s i ) i f   l = 1 L ( q ^ i , l H q ^ i , l s i ) > 0 i 0 i f   l = 1 L ( q ^ i , l H q ^ i , l s i ) 0
Let γ = [ γ 1 , γ 2 , , γ K ] T . In the ( t + 1 ) -th iteration, we need to identify the update of which hyper-parameter makes the biggest contribution to the maximization of ( Γ ) . The serial number of the corresponding hyper-parameter is j . Compared with other hyper-parameters, the update of the hyper-parameter γ j makes the biggest contribution to the maximization of ( Γ ) . Therefore, j can be selected by the following equation
j = arg max i , 1 i K ( Γ ) = arg max i , 1 i K i + ( γ i )
where, herein Γ = diag ( γ ) and γ = [ γ 1 ( t ) , , γ i , , γ K ( t ) ] T .
In (29), i , i need to be calculated in each iteration. To improve efficiency, we define
Δ ( i ) i + ( γ i ) ( Γ ( t ) ) = i + ( γ i ) i ( γ i ( t ) ) = ( γ i ) ( γ i ( t ) )
Equation (29) can be simplified as the following equation since ( Γ ( t ) ) is a constant in the ( t + 1 ) -th iteration.
j = argmax i , 1 i K ( Γ ( t ) ) + Δ ( i ) = argmax i , 1 i K Δ ( i ) = argmax i , 1 i K ( γ i ) ( γ i ( t ) )
Compared with (29), it is more convenient to select j with (31) and ( γ i ( t ) ) ,   i been calculated in the last iteration. If j has been identified, replace γ j with γ j while fixing { γ i | 1 i K , i j } .
γ j ( t + 1 ) = γ j
γ i ( t + 1 ) = γ i ( t ) ,   1 i K , i j
As mentioned in the introduction, we select the atoms that are in the clutter subspace and reserve these atoms in the new dictionary. Therefore, the non-zero values in γ are reserved in Ω and the corresponding atoms are reserved in the new dictionary ψ .
Ω diag ( γ ω 1 , γ ω 2 , , γ ω J )
ψ [ ν ω 1 , ν ω 2 , , ν ω J ]
where J is the number of atoms in ψ and
{ ω 1 , ω 2 , , ω J } = { i | 1 i K ,   γ i 0 }
We initialize γ ( 0 ) = 0 , namely, Ω ( 0 ) = and ψ ( 0 ) = at the beginning of the proposed algorithm. With the knowledge of γ j ( t ) and γ j ( t + 1 ) , an efficient sequential addition and deletion of dictionary atoms can be taken to estimate ψ ( t + 1 ) and Ω ( t + 1 ) . If γ j ( t ) = 0 and γ j ( t + 1 ) > 0 , ψ ( t + 1 ) = [ ψ ( t ) , ν j ] and Ω ( t + 1 ) = [ Ω ( t ) γ j ( t + 1 ) ] ; if γ j ( t ) > 0 and γ j ( t + 1 ) > 0 , ψ ( t + 1 ) = ψ ( t ) and replace γ j ( t ) with γ j ( t + 1 ) in Ω ( t ) ; if γ j ( t ) > 0 and γ j ( t + 1 ) = 0 , delete ν j from ψ ( t ) and delete γ j ( t ) from Ω ( t ) ; if γ j ( t ) = 0 and γ j ( t + 1 ) = 0 , stop iteration because ( Γ ) has already converged.

4.3. Fast Computation of { s ^ i , q ^ i , l }

Calculating the matrix inversion of C i ,   i brings a heavy computational load when { s ^ i , q ^ i , l } are updated. Fast computation of { s ^ i , q ^ i , l } is introduced as follows.
Define
s i ν i H C 1 ν i   ,   q i , l ν i H C 1 x l
It is computationally efficient to calculate { s ^ i , q ^ i , l } with { s i , q i , l } because it is more convenient to calculate the matrix inversion of only C than that of all C i ,   i .
With (20), we can obtain
s i = ν i H C 1 ν i = ν i H ( C i 1 C i 1 ν i ( γ i 1 + ν i H C i 1 ν i ) 1 ν i H C i 1 ) ν i = s ^ i s ^ i ( γ i 1 + s ^ i ) 1 s ^ i = s ^ i 1 + γ i s ^ i
Thus, s ^ i can be calculated with s i .
s ^ i = s i 1 γ i s i
Similarly, q ^ i , l can be also calculated with q i , l .
q ^ i , l = ( 1 + γ i s ^ i ) q i , l = q i , l 1 γ i s i
To further improve the computational efficiency, we introduce the following approach to reduce computational complexities of { s i , q i , l } . With { Ω , ψ } , (15) can be also expressed as
C = σ 2 I M N + ψ Ω ψ H
Define α σ 2 , and calculate { s i , q i , l } with matrix inversion lemmas.
s i = α ν i H ν i α 2 ν i H ψ Σ ψ H ν i
q i , l = α ν i H x l α 2 ν i H ψ Σ ψ H x l
where Σ herein represents the covariance of χ l , l . χ l , l is defined as the elements in a l , l whose corresponding hyper-parameters are non-zeros in the current iteration.
Σ = ( σ 2 ψ H ψ + Ω 1 ) 1
The mean of χ l , l is expressed as
μ l = σ 2 Σ ψ H x l , l
With (42)~(45), we utilize { s i , q i , l , Σ , μ l } ( t ) calculated in the t -th iteration to update { s i , q i , l , Σ , μ l } ( t + 1 ) in the ( t + 1 ) -th iteration, which improves the computational efficiency. The updated formulas are listed in Appendix A.

4.4. Calculation of the STAP Filter Weight Vector

The CNCM R can be estimated by
R = 1 L l = 1 L i = 1 J μ i , l 2 υ ω i υ ω i H + β σ 2 Ι
where β 1 is a load factor.
The optimal STAP weight vector can be given by
w = R 1 s t s t H R 1 s t
where s t is the steering vector of the target.
The proposed algorithm is shown in Algorithm 1. To reduce storage load, we do not store D . When traversing the i -th grid, we use Formula (7) to generate the corresponding space-time steering vector.
Algorithm 1. Pseudocode for the proposed algorithm.
Step 1: Input: the data X, σ 2 .
Step 2: Initialize: γ ( 0 ) = 0 , C = C i = σ 2 I M N , and Ω ( 0 ) = , ψ ( 0 ) = .
Step 3: While not converged do
    Obtain all γ i , i by (28), and exploit (31) to find j -th hyper-parameter which needs to be updated in the current iteration.
    If γ j ( t + 1 ) > 0 and γ j ( t ) = 0 ,
       ψ ( t + 1 ) = [ ψ ( t ) , ν j ] , and Ω ( t + 1 ) = [ Ω ( t ) γ j ( t + 1 ) ] .
    If γ j ( t + 1 ) > 0 and γ j ( t ) > 0 ,
       ψ ( t + 1 ) = ψ ( t ) , and replace γ j ( t ) with γ j ( t + 1 ) .
    If γ j ( t + 1 ) > 0 and γ j ( t ) > 0 ,
      delete ν j from ψ ( t ) , and delete γ j ( t ) from Ω ( t ) .
    end
    Update { Σ   ,   μ   ,   s i   ,   q i , l } ,   i referring to Appendix A.
   end while
Step 4: Estimate the CNCM R by (46)
Step 5: Compute the space-time adaptive weight w using (47).
Step 6: The output of the space-time filter is w H X .

5. Analysis of Complexity, Storage and Convergence

5.1. Complexity Analysis

To reflect the low computational complexity of the proposed algorithm, some MMV SR-STAP algorithms, including the multiple focal underdetermined system solver (M-FOCUSS) and the M-SBL, are utilized as a comparison with the proposed algorithm. The number of multiplications for a single iteration is utilized as a measurement of the computational complexity. The dimension of ψ is M N × J and J is not fixed in each iteration. Because of the sparsity of clutter, J < M N . Suppose that J is equal to twice the rank of clutter when comparing the computational complexity of different algorithms in the simulations. The computational complexities of different algorithms are summarized in Table 1. For convenience, we set φ s = φ d , and define
φ φ s = φ d
Figure 1 illustrates the function of the computational complexities to φ in the side-looking radar case. In the simulations, M = N = 10 , L = 6 and β = 4 v / λ f PRF = 1 . From Figure 1, the computational complexity of the proposed algorithm is far less than that of the M-FOCUSS and the M-SBL algorithm.

5.2. Storage Analysis

In order to reflect the low storage requirement of the proposed algorithm more intuitively, we compare the proposed algorithm with the M-SBL algorithms. We need to store D M N × K and Γ K × K in the M-SBL algorithm while we store ψ M N × J , Ω J × J in the proposed algorithm. In the M-SBL algorithm, the dimension of Σ in (12) is K × K and the dimension of μ l in (13) is K × 1 . However, in the proposed algorithm, the dimension of Σ in (44) is J × J and the dimension of μ l in (45) is J × 1 . For example, when M = N = 10 and φ s = φ d = 20 , K = φ d M φ s N = 40,000 . In this case, the storage load in the M-SBL algorithm is heavy. Because of the sparsity of clutter, we can observe that J < M N = 100 . The storage load in the proposed algorithm is far less than the M-SBL algorithm.

5.3. Convergence Analysis

According to [32], ( Γ ) has an upper bound. We can conclude from (27) and (28) that ( γ j ( t + 1 ) ) = ( γ j ) ( γ j ( t ) ) and ( Γ ( t + 1 ) ) ( Γ ( t ) ) , t . ( Γ ) has an upper bound while it is also a monotonically increasing function, which means the proposed algorithm converges.

6. Performance Assessment

In this section, we verify the performance of the proposed algorithm and the other SR-STAP algorithms with the simulated data. The parameters of the airborne radar system are listed in Table 2. The SR-STAP algorithms for comparison are the M-FOCUSS and the M-SBL. We utilize the improvement factor (IF) as the measurement of performance.
IF = | w H s | 2 w H R w tr ( R ) s H s

6.1. Comparison of Clutter Spectrums Estimated by SR-STAP Algorithms

Clutter spectrums estimated by the proposed algorithm and other SR-STAP algorithms are compared in the absence and presence of off-grid problems. We consider three different cases: (i) A side-looking radar without off-grid problems (the platform velocity is 150 m/s and β = 1 ); (ii) a side-looking radar with off-grid problems (the platform velocity is 180 m/s and β = 1.2 ); (iii) a forward-looking radar (the platform velocity is 180 m/s). When off-grid problems occur, the effects of off-grid problems are mitigated by increasing the value of φ . We also set φ = 15 and φ = 25 in the M-FOCUSS and M-SBL algorithms. However, either their running time is far beyond our acceptable range, the computer crashes, or their results are wrong because of the RIP condition. Therefore, their results are not shown when φ = 15 and φ = 25 .
(i)
A side-looking radar without off-grid problems
We compare the high-resolution spectrum estimated by the M-FOCUSS, the M-SBL, and the proposed algorithm. In the absence of off-grid problems, the value of φ can be set to an integer between 2 and 5. In the experiment, φ = 5 . From the below figures in Figure 2, the clutter spectrums estimated by the M-SBL and the proposed algorithms are close to the ideal clutter spectrum in terms of the power and location of the clutter, which indicates the exact clutter power can be estimated by the proposed algorithm.
(ii)
A side-looking radar with off-grid problems
Considering that the radar is a side-looking radar with off-grid problems, we set φ equal to 5, 15, and 25, respectively. When φ = 15 or φ = 25 in the M-FOCUSS and M-SBL algorithms, either their running time is far beyond our acceptable range, or the computer crashes, or their results are wrong because of the RIP condition. Therefore, their results are not shown when φ = 15 and φ = 25 . From the below figures in Figure 3, the clutter spectrum estimated by the proposed algorithms is close to the ideal clutter spectrum in terms of the power and location of the clutter when φ = 25 , which indicates that the exact clutter power can be estimated by the proposed algorithm with sufficiently dense grids. The experiment shows that the proposed algorithm can effectively mitigate the off-grid effects.
(iii)
A forward-looking radar.
Considering that the radar is a forward-looking radar, we set φ equal to 5, 15, and 25, respectively. From the below figures in Figure 4, the clutter spectrum estimated by the proposed algorithms is close to the ideal clutter spectrum in terms of the power and location of the clutter when φ = 25 , which indicates that the exact clutter power can be estimated by the proposed algorithm with sufficiently dense grids. The experiment shows that the proposed algorithm can effectively mitigate the off-grid effects.

6.2. Comparison of IF Curves with SR-STAP Algorithms

In this experiment, we compare the clutter suppression performance of the proposed method with the M-FOCUSS and the M-SBL algorithms in the presence and absence of off-grid problems. Furthermore, we also consider the ideal case and the non-ideal case. Amplitude Gaussian error (standard deviation 0.03) and phase random error (standard deviation 2°) are different in all directions in the non-ideal case.
(i)
A side-looking radar without off-grid problems
In the absence of off-grid problems, the value of φ can be set equal to 5. From Figure 5, we note that the IF curves obtained by the M-SBL and the proposed algorithms are both close to the optimal IF curves obtained by the exact CNCM, which indicates that two algorithms can recover more exact clutter sources. The experiment shows that the proposed algorithm can also obtain good performance in the absence of off-grid problems.
(ii)
A side-looking radar with off-grid problems
Considering that the radar is a side-looking radar with off-grid problems, we set φ equal to 5, 15, and 25, respectively. From Figure 6, we find that increasing φ can effectively improve the clutter suppression performance in the presence of off-grid problems. When φ = 25 , the IF curves obtained by the proposed algorithm are higher and narrower than the other IF curves. The experiment shows that the proposed algorithm can effectively mitigate the off-grid effects.
(iii)
A forward-looking radar
Considering that the radar is a forward-looking radar, we set φ equal to 5, 15, and 25, respectively. From Figure 7, we find that increasing φ can effectively improve the clutter suppression performance in the presence of off-grid problems. When φ = 25 , the IF curves obtained by the proposed algorithm are higher and narrower than the other IF curves. The experiment shows that the proposed algorithm can effectively mitigate the off-grid effects. Theoretically, as long as the grid points are dense enough, the off-grid problems can be solved with the proposed algorithm.

6.3. Comparison of Running Time with SR-STAP Algorithms

The symbol represents the running time is far beyond our acceptable range (1 hour) or beyond the capacity of the computer. We consider three different cases: (i) A side-looking radar without off-grid problems, (ii) a side-looking radar with off-grid problems, and (iii) a forward-looking radar. The running time of the algorithm is not only affected by the convergence performance of the algorithm itself, but also the computational complexity and the burden of the storage on the computer. The computer in the simulations is equipped with Intel(R) Xeon(R) Gold 6140 CPU @ 2.30GHz 2.29GHz. The running time of the proposed algorithm in Table 3, Table 4 and Table 5 is far less than the traditional SR-STAP algorithms, especially when dense grids are exploited. One hundred Monte Carlo trials were performed to obtain the average running time.
(i)
A side-looking radar without off-grid problems
Table 3. Running time of different algorithms for a side-looking radar without off-grid problems.
Table 3. Running time of different algorithms for a side-looking radar without off-grid problems.
AlgorithmThe Average Running Time (s)
M-FOCUSS1.05
M-SBL26.32
the proposed algorithm0.88
(ii)
A side-looking radar with off-grid problems
Table 4. Running time of different algorithms for a side-looking radar with off-grid problems.
Table 4. Running time of different algorithms for a side-looking radar with off-grid problems.
AlgorithmThe Average Running Time (s)
φ = 5 φ = 15 φ = 25
M-FOCUSS1.1372.23
M-SBL29.79
the proposed algorithm1.063.3612.14
(iii)
A forward-looking radar
Table 5. Running time of different algorithms for a forward-looking radar.
Table 5. Running time of different algorithms for a forward-looking radar.
AlgorithmThe Average Running Time (s)
φ = 5 φ = 15 φ = 25
M-FOCUSS1.0264.14
M-SBL30.62
the proposed algorithm1.054.0312.94

7. Conclusions

A novel SR-STAP algorithm is proposed to mitigate off-grid effects. The proposed algorithm is based on the Bayes criterion and mitigates off-grid effects when the grid interval is sufficiently small. We pick out which atoms are in the subspace of clutter, which means the proposed algorithm does not need to satisfy the restricted isometry property (RIP) condition. Moreover, the complexity and storage requirement of the proposed algorithm is low, and its convergence can be promised.

Author Contributions

Conceptualization, T.W.; methodology, C.L.; validation, C.L., K.L. and X.Z.; writing—original draft preparation, C.L.; writing—review and editing, K.L. and X.Z.; visualization, C.L., K.L. and X.Z.; supervision, T.W.; project administration, T.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key R&D Program of China, grant number 2021YFA1000400.

Data Availability Statement

In this section, please provide details regarding where data supporting reported results can be found, including links to publicly archived datasets analyzed or generated during the study. Please refer to suggested Data Availability Statements in section “MDPI Research Data Policies” at https://www.mdpi.com/ethics. You might choose to exclude this statement if the study did not report any data.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

In order to make the following formulas more intuitive and concise, we remove ( t ) from quantities in the t -th iteration, and updated quantities in the ( t + 1 ) -th iteration are denoted by a tilde ( ~ ). As stated in the text, j represents the serial number of the hyper-parameter that needs to be updated in the ( t + 1 ) -th iteration. j also represents the position index of ν j in D . When ν j is already in ψ , g is defined as the position index of ν j in ψ .
  • When γ ˜ j > 0   & &   γ j = 0 ,
    Δ ( j ) = ( γ ˜ j ) = l = 1 L [ q ^ j , l H ( γ ˜ j 1 + s j ) 1 q ^ j , l ln | 1 + γ ˜ j s j | ]
    Σ ˜ = [ Σ + α 2 z Σ ψ H ν j ν j H ψ Σ α z Σ ψ H ν j α z ( Σ ψ H ν j ) H z ]
    μ ˜ l = [ μ l α z q j , l Σ ψ H ν j z q j , l ]
    s ˜ i = s i α 2 z ν i H e j e j H ν i , i
    q ˜ i , l = q i , l α z q i , l ν i H e j , i
    where z = ( γ ˜ j 1 + s j ) 1 and e j = ν j α ψ Σ ψ H ν j .
  • When γ ˜ j > 0   & &   γ j > 0 ,
    Δ L j = γ ˜ j γ j = l = 1 L q ^ j , l H γ ˜ j 1 + s j 1 q ^ j , l q ^ j , l H γ j 1 + s j 1 q ^ j , l ln 1 + γ ˜ j s j + ln 1 + γ j s j
    Σ ˜ = Σ κ Σ g Σ g H
    μ ˜ l = μ κ μ g , l Σ g
    s ˜ i = s i + α 2 κ ν i H ψ Σ g Σ g H ψ H ν i , i
    q ˜ i , l = q i , l + α κ μ g , l v i H ψ Σ g , i
    where κ = γ j γ ˜ j Δ γ i + Σ g g 1 , Δ γ j = γ ˜ j γ j and Σ g = Σ : , g .
  • When γ ˜ j = 0   & &   γ j > 0 ,
    Δ L j = γ j = l = 1 L q ^ j , l H γ j 1 + s j 1 q ^ j , l + ln 1 + γ j s j
    Σ ˜ = Σ Σ g g 1 Σ g Σ g H
    μ ˜ l = μ Σ g g 1 μ g , l Σ g
    s ˜ i = s i + α 2 Σ g g 1 ν i H ψ Σ g Σ g H ψ H ν i , i
    q ˜ i , l = q i , l + α Σ g g 1 μ g , l ν i H ψ Σ g , i

References

  1. Ward, J. Space-Time Adaptive Processing for Airborne Radar; Technical Report; MIT Lincoln Laboratory: Lexington, KY, USA, 1998. [Google Scholar]
  2. Reed, I.S.; Mallet, J.D.; Brennan, L.E. Rapid convergence rate in adaptive arrays. IEEE Trans. Aerosp. Electron. Syst. 1974, 10, 853–863. [Google Scholar] [CrossRef]
  3. Donoho, D.L.; Elad, M.; Temlyakov, V.N. Stable recovery of sparse overcomplete representations in the presence of noise. IEEE Trans. Inf. Theory 2006, 52, 6–18. [Google Scholar] [CrossRef]
  4. Sun, K.; Zhang, H.; Li, G.; Meng, H.D.; Wang, X.Q. A novel STAP algorithm using sparse recovery technique. IEEE Int. Geosci. Remote Sens. Symp. 2009, 1, 3761–3764. [Google Scholar]
  5. Yang, Z.C.; Li, X.; Wang, H.Q.; Jiang, W.D. On clutter sparsity analysis in space-time adaptive processing airborne radar. IEEE Geosci. Remote Sens. Lett. 2013, 10, 1214–1218. [Google Scholar] [CrossRef]
  6. Sen, S. Low-rank matrix decomposition and spatial-temporal sparse recovery for STAP radar. IEEE J. Sel. Top. Signal Process. 2015, 9, 1510–1523. [Google Scholar] [CrossRef]
  7. Yang, Z.; Wang, Z.; Liu, W.; de Lamare, R.C. Reduced-dimension space-time adaptive processing with sparse constraints on beam-Doppler selection. Signal Process. 2019, 157, 78–87. [Google Scholar] [CrossRef]
  8. Zhang, W.; An, R.; He, N.; He, Z.; Li, H. Reduced dimension STAP based on sparse recovery in heterogeneous clutter environments. IEEE Trans. Aerosp. Electron. Syst. 2019, 56, 785–795. [Google Scholar] [CrossRef]
  9. Liu, C.; Wang, T.; Zhang, S.G.; Ren, B. Clutter suppression based on iterative reweighted methods with multiple measurement vectors for airborne radar. IET Radar Sonar Navig. 2022, 16, 1–14. [Google Scholar] [CrossRef]
  10. Candes, M.; Wakin, M.; Boyd, S. Enhancing sparsity by reweighted 1 minimization. J. Fourier Anal. Appl. 2008, 5, 877–905. [Google Scholar] [CrossRef]
  11. Tipping, M.E. Sparse Bayesian learning and the relevance vector machine. J. Mach. Learn. 2001, 1, 211–244. [Google Scholar]
  12. Wipf, D.P.; Rao, B.D. Sparse Bayesian learning for basis selection. IEEE Trans. Signal Process. 2004, 52, 2153–2164. [Google Scholar] [CrossRef]
  13. Wipf, D.P.; Rao, B.D. An empirical Bayesian strategy for solving the simultaneous sparse approximation problem. IEEE Trans. Signal Process. 2007, 55, 3704–3716. [Google Scholar] [CrossRef]
  14. Tipping, M.E.; Faul, A.C. Fast marginal likelihood maximization for sparse Bayesian models. In Proceedings of the Ninth International Workshop on Artificial Intelligence and Statistics, Key West, FL, USA, 3–6 January 2003; Volume 1, pp. 276–283. [Google Scholar]
  15. Ji, S.H.; Xue, Y.; Carin, L. Bayesian compressive sensing. IEEE Trans. Signal Process. 2008, 56, 2346–2356. [Google Scholar] [CrossRef]
  16. Ji, S.H.; Dunson, D.; Carin, L. Multi-task compressive sensing. IEEE Trans. Signal Process. 2009, 57, 92–106. [Google Scholar] [CrossRef]
  17. Babacan, S.D.; Molina, R.; Katsaggelos, A.K. Bayesian compressive sensing using Laplace priors. IEEE Trans. Image Process. 2010, 19, 53–63. [Google Scholar] [CrossRef]
  18. Wu, Q.S.; Zhang, Y.M.; Amin, M.G.; Himed, B. Complex multitask Bayesian compressive sensing. In Proceedings of the 2014 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), Florence, Italy, 4–9 May 2014. [Google Scholar]
  19. Serra, J.G.; Testa, M.; Katsaggelos, A.K. Bayesian K-SVD using fast variational inference. IEEE Trans. Image Process. 2017, 26, 3344–3359. [Google Scholar] [CrossRef]
  20. Ma, Z.Q.; Dai, W.; Liu, Y.M.; Wang, X.Q. Group sparse Bayesian learning via exact and fast marginal likelihood maximization. IEEE Trans. Signal Process. 2017, 65, 2741–2753. [Google Scholar] [CrossRef]
  21. Liu, C.; Wang, T.; Zhang, S.G.; Ren, B. A fast space-time adaptive processing algorithm based on sparse Bayesian learning for airborne radar. Sensors 2022, 22, 2664. [Google Scholar] [CrossRef]
  22. Duan, K.Q.; Wang, Z.T.; Xie, W.C.; Chen, H.; Wang, Y.L. Sparsity-based STAP algorithm with multiple measurement vectors via sparse Bayesian learning strategy for airborne radar. IET Signal Process. 2017, 11, 544–553. [Google Scholar] [CrossRef]
  23. Wang, Z.T.; Xie, W.C.; Duan, K.Q. Clutter suppression algorithm base on fast converging sparse Bayesian learning for airborne radar. Signal Process. 2017, 130, 159–168. [Google Scholar] [CrossRef]
  24. Yang, X.P.; Sun, Y.Z.; Yang, J.; Long, T.; Sarkar, T.K. Discrete Interference suppression method based on robust sparse Bayesian learning for STAP. IEEE Access 2019, 10, 26740–26751. [Google Scholar] [CrossRef]
  25. Duan, K.Q.; Liu, W.J.; Duan, G.Q.; Wang, Y.L. Off-grid effects mitigation exploiting knowledge of the clutter ridge for sparse recovery STAP. IET Radar Sonar Navig. 2018, 12, 557–564. [Google Scholar] [CrossRef]
  26. You, K.T.; Guo, W.B.; Liu, Y.L.; Wang, W.B.; Sun, Z. Grid evolution: Joint dictionary learning and sparse Bayesian recovery for multiple off-grid targets localization. IEEE Commun. Lett. 2018, 22, 2068–2071. [Google Scholar] [CrossRef]
  27. Dai, J.S.; Bao, X.; Xu, W.C.; Chang, C.Q. Root sparse Bayesian learning for off-grid DOA estimation. IEEE Signal Process. Lett. 2017, 24, 46–50. [Google Scholar] [CrossRef]
  28. Fang, J.; Wang, F.Y.; Shen, Y.N.; Li, H.B.; Blum, R.S. Super-resolution compressed sensing for line spectral estimation: An iterative reweighted approach. IEEE Trans. Signal Process. 2016, 64, 4649–4662. [Google Scholar] [CrossRef]
  29. Fang, J.; Shen, Y.N.; Li, H.B.; Li, S.Q. Super-resolution compressed sensing: An iterative reweighted algorithm for joint parameter learning and sparse signal recovery. IEEE Signal Process. Lett. 2014, 21, 761–766. [Google Scholar]
  30. Li, Z.H.; Zhang, Y.S.; He, X.Y.; Guo, Y.D. Low-complexity off-grid STAP algorithm based on local search clutter subspace estimation. IEEE Geosci. Remote Sens. Lett. 2018, 15, 1862–1865. [Google Scholar] [CrossRef]
  31. Yuan, H.D.; Xu, H.; Duan, K.Q.; Xie, W.C.; Liu, W.J.; Wang, Y.L. Sparse Bayesian learning-based space-time adaptive processing with off-grid self-calibration for airborne radar. IEEE Access 2018, 6, 47296–47307. [Google Scholar] [CrossRef]
  32. Wipf, D.; Nagarajan, S. A new view of automatic relevance determination. In Advances in Neural Information Processing Systems 20; MIT Press: New York, NY, USA, 2008. [Google Scholar]
Figure 1. Computational complexity versus φ for a single iteration.
Figure 1. Computational complexity versus φ for a single iteration.
Remotesensing 14 03906 g001
Figure 2. Comparison between the clutter angle-Doppler spectrum estimated by different algorithms for a side-looking radar without off-grid problems. (a) Ideal clutter spectrum; (b) M-FOCUSS, φ = 5 ; (c) M-SBL, φ = 5 ; (d) the proposed algorithm, φ = 5 .
Figure 2. Comparison between the clutter angle-Doppler spectrum estimated by different algorithms for a side-looking radar without off-grid problems. (a) Ideal clutter spectrum; (b) M-FOCUSS, φ = 5 ; (c) M-SBL, φ = 5 ; (d) the proposed algorithm, φ = 5 .
Remotesensing 14 03906 g002aRemotesensing 14 03906 g002b
Figure 3. Comparison between the clutter angle-Doppler spectrum estimated by different algorithms for a side-looking radar with off-grid problems. (a) Ideal clutter spectrum; (b) M-FOCUSS, φ = 5 ; (c) M-SBL, φ = 5 ; (d) the proposed algorithm, φ = 5 ; (e) the proposed algorithm, φ = 15 ; (f) the proposed algorithm, φ = 25 .
Figure 3. Comparison between the clutter angle-Doppler spectrum estimated by different algorithms for a side-looking radar with off-grid problems. (a) Ideal clutter spectrum; (b) M-FOCUSS, φ = 5 ; (c) M-SBL, φ = 5 ; (d) the proposed algorithm, φ = 5 ; (e) the proposed algorithm, φ = 15 ; (f) the proposed algorithm, φ = 25 .
Remotesensing 14 03906 g003aRemotesensing 14 03906 g003b
Figure 4. Comparison between the clutter angle-Doppler spectrum estimated by different algorithms for a forward-looking radar. (a) Ideal clutter spectrum; (b) M-FOCUSS, φ = 5 ; (c) M-SBL, φ = 5 ; (d) the proposed algorithm, φ = 5 ; (e) the proposed algorithm, φ = 15 ; (f) the proposed algorithm, φ = 25 .
Figure 4. Comparison between the clutter angle-Doppler spectrum estimated by different algorithms for a forward-looking radar. (a) Ideal clutter spectrum; (b) M-FOCUSS, φ = 5 ; (c) M-SBL, φ = 5 ; (d) the proposed algorithm, φ = 5 ; (e) the proposed algorithm, φ = 15 ; (f) the proposed algorithm, φ = 25 .
Remotesensing 14 03906 g004
Figure 5. IF versus normalized Doppler frequency for a side-looking radar without off-grid problems. (a) Ideal case; (b) non-ideal case with array errors.
Figure 5. IF versus normalized Doppler frequency for a side-looking radar without off-grid problems. (a) Ideal case; (b) non-ideal case with array errors.
Remotesensing 14 03906 g005
Figure 6. IF versus normalized Doppler frequency for a side-looking radar with off-grid problems. (a) Ideal case; (b) non-ideal case with array errors.
Figure 6. IF versus normalized Doppler frequency for a side-looking radar with off-grid problems. (a) Ideal case; (b) non-ideal case with array errors.
Remotesensing 14 03906 g006
Figure 7. IF versus normalized Doppler frequency for a forward-looking radar. (a) Ideal case; (b) non-ideal case with array errors.
Figure 7. IF versus normalized Doppler frequency for a forward-looking radar. (a) Ideal case; (b) non-ideal case with array errors.
Remotesensing 14 03906 g007
Table 1. The computational complexity for a single iteration.
Table 1. The computational complexity for a single iteration.
AlgorithmComputational Complexity
M-FOCUSS o ( M N K L + ( M N ) 3 + 2 ( M N ) 2 K )
M-SBL o ( M N K L + ( M N ) 3 + 2 ( M N ) 2 K + M N K 2 )
the proposed algorithm o ( ( M N + 2 M N J + J 2 + L ) K )
Table 2. Parameters of airborne radar system.
Table 2. Parameters of airborne radar system.
ParametersSymbolsValue
Distance between elements d 0.15 m
Wavelength λ 0.3 m
Platform height H 9000 m
Number of pulses M 8
Number of channels N 8
Pulse repetition frequency f PRF 2000 Hz
Range sampling frequency f s 2.5 MHz
Clutter to noise ratio CNR 40 dB
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liu, C.; Wang, T.; Liu, K.; Zhang, X. A Novel Sparse Bayesian Space-Time Adaptive Processing Algorithm to Mitigate Off-Grid Effects. Remote Sens. 2022, 14, 3906. https://doi.org/10.3390/rs14163906

AMA Style

Liu C, Wang T, Liu K, Zhang X. A Novel Sparse Bayesian Space-Time Adaptive Processing Algorithm to Mitigate Off-Grid Effects. Remote Sensing. 2022; 14(16):3906. https://doi.org/10.3390/rs14163906

Chicago/Turabian Style

Liu, Cheng, Tong Wang, Kun Liu, and Xinying Zhang. 2022. "A Novel Sparse Bayesian Space-Time Adaptive Processing Algorithm to Mitigate Off-Grid Effects" Remote Sensing 14, no. 16: 3906. https://doi.org/10.3390/rs14163906

APA Style

Liu, C., Wang, T., Liu, K., & Zhang, X. (2022). A Novel Sparse Bayesian Space-Time Adaptive Processing Algorithm to Mitigate Off-Grid Effects. Remote Sensing, 14(16), 3906. https://doi.org/10.3390/rs14163906

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop