Next Article in Journal
Unbiased 3D Semantic Scene Graph Prediction in Point Cloud Using Deep Learning
Previous Article in Journal
Design and Development of a Flexible 3D-Printed Endoscopic Grasping Instrument
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Clutter Suppression Algorithm with Joint Intrinsic Clutter Motion Errors Calibration and Off-Grid Effects Mitigation in Airborne Passive Radars

1
School of Information and Electronic Engineering, Hunan City University, No. 518, Yingbin East Road, Yiyang 413000, China
2
All-Solid-State Energy Storage Materials and Devices Key Laboratory of Hunan Province, Hunan City University, Yiyang 413000, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2023, 13(9), 5653; https://doi.org/10.3390/app13095653
Submission received: 27 March 2023 / Revised: 28 April 2023 / Accepted: 30 April 2023 / Published: 4 May 2023

Abstract

:
In an airborne passive radar, multipath (MP) clutter, which is caused by MP signals contained in the contaminated reference signal, degrades the space-time adaptive processing (STAP) performance. The MP clutter suppression algorithm before STAP can mitigate the influence of impure reference signals. However, the performances of the existing MP clutter suppression methods deteriorate when the intrinsic clutter motion (ICM) exists because the sparse model of MP clutter is disturbed. To eliminate the impacts of ICM on MP clutter suppression, a joint optimization algorithm is developed for airborne passive radar. Firstly, the sparse model of MP clutter is modified by taking ICM fluctuation into account. Subsequently, the joint optimization function of the ICM fluctuation and MP clutter profile is derived. Finally, based on the local search technique, MP clutter is suppressed with ICM error calibration and off-grid effects mitigation. A range of simulations verify the reliability and superiority of the proposed method.

1. Introduction

Airborne passive radars exploit the existing transmitters as emitters of opportunity, having attracted widespread attention in recent years [1,2,3,4,5,6,7]. Compared with airborne active radars, the lack of dedicated transmitters offers the advantages of low cost, strong survivability, and antijamming. However, the moving airborne platform leading to the received ground clutter is extended in terms of the Doppler frequencies and angle. Thus, it is difficult for a traditional spatial or temporal filter to detect targets.
Space-time adaptive processing (STAP), which uses multichannel spatial information across a sequence of temporal pulses to form a two-dimensional filter, is a crucial technology in airborne passive radar for clutter cancellation. However, the application of STAP is restricted by practical problems, such as excessively high computational complexity and a large number of the required training snapshots. Subsequently, extensive research on STAP mainly considered how to reduce the computational burden and enhance the convergence speed. However, the practical implementation of STAP in airborne passive radar is based on the assumption that the reference signal is uncontaminated. When this condition is not satisfied, the clutter caused by the direct-path (DP) signal (DP clutter) and multipath (MP) signal (MP clutter) affects the covariance matrix estimation. Consequently, the target self-nulling phenomenon may exist [8].
The MP clutter cancellation approach before STAP can eliminate the influence of impure reference signals in airborne passive radars. Existing MP clutter suppression techniques can be categorized into two groups. One is signal equalization, which relies heavily on precise prior information [9]. Another is the sparse reconstruction of MP clutter that does not require prior knowledge [8,10,11]. However, the pulse-to-pulse fluctuations in MP signals, also known as the intrinsic clutter motion (ICM) [12], degrade the reconstruction quality. Thus, it is a challenge to suppress MP clutter with ICM errors in airborne passive radar systems.
In the current paper, to enhance the MP clutter suppression performance with ICM errors, we present a novel suppression algorithm with ICM error calibration and off-grid effects mitigation. In the proposed algorithm, the ICMs of clutter patches (ICM-CP) and MP scatterers (ICM-MP) are considered. After using reference signals for matched filtering, the DP clutter snapshot involves the temporal fluctuation vector caused by ICM-CP, and the MP clutter snapshot involves the temporal fluctuation vector caused by both ICM-CP and ICM-MP. Consequently, the MP clutter is not only related to the DP clutter snapshot and the modified steering vector, which are mentioned in SRA and LSA, but also related to the ICM-MP fluctuation. Thus, the existing sparse model of MP clutter is disturbed in the presence of ICM-MP. Based on the above analyses, an MP clutter suppression algorithm based on joint optimization is designed. The procedure of the proposed algorithm is described as follows: Initially, the ICM-MP fluctuation is introduced, and the sparse model of MP clutter is revised. Subsequently, by adding a constraint to ICM-MP fluctuation, we modify a traditional reconstruction issue into a joint optimization issue. Thirdly, the ICM-MP fluctuation and MP clutter profile are estimated simultaneously based on the local search technique. Finally, the MP clutter is suppressed with ICM-MP error calibration and off-grid effects mitigation. After the MP clutter cancellation algorithm, the existing STAP is utilized to mitigate the impacts of ICM-DP and cancel DP clutter.
The main contributions of the paper are summarized as follows:
(i)
We develop a new sparse model of MP clutter and joint optimization problems by considering ICM-MP fluctuation.
(ii)
We propose a joint optimization algorithm for improving the MP clutter suppression performance when ICM-MP fluctuation exists, where the local search technique is incorporated to mitigate off-grid effects.
(iii)
The calibration step is efficiently performed, and we discuss the feasibility of the proposed algorithm from practical implementation and computational efficiency perspectives.
The remainder of this paper is organized as follows: We first discuss the related work in Section 2, and the signal model, including ICM for airborne passive radar, is introduced in Section 3. Subsequently, we review the existing optimization problems in SRA and LSA, and then the proposed joint optimization algorithm for MP clutter suppression is derived in Section 4. Section 5 presents some simulated experiments demonstrating the performance and validity of the proposed approach. Finally, Section 6 gives the relevant conclusions.
Notation: In this paper, the Kronecker product and Hadamard product are expressed by the symbols and , respectively. The symbol E indicates the mathematical expectation, and Re represents the real part of the argument. Τ , * and H denotes the transposition, complex conjugation, and conjugate transposition, respectively. diag denotes the diagonal matrix, I A is the A × A identity matrix, Toeplitz is the Toeplitz matrix, and 1 B × C is the B × C matrix with all elements being 1. In addition, F , 2 and 1 represent the Frobenius-norm, l2-norm, and l1-norm operations, respectively.

2. Related Work

STAP has been studied extensively in previous works. Initially, the fully adaptive STAP, where adaptive weight is calculated using all elements and pulses, was considered [12,13]. Afterwards, a variety of STAP approaches have been developed to reduce the computational burden and enhance the convergence speed. To improve the detection ability in heterogeneous environments, Mevin et al. proposed a knowledge-aided STAP technique in which the prior information was incorporated within the conventional training strategies [14]. A similar idea was used in other works [15,16,17,18]. The reduced-rank and reduced-dimension STAP, such as subspace projection and extended factored algorithms, have been widely considered since the training requirements were reduced to twice the degrees of freedom or the clutter rank [19,20,21,22,23,24,25,26]. To further reduce the training requirements, Yang et al. proposed the sparse representation STAP [27]. By utilizing intrinsic sparsity and compressive sensing techniques to recover the clutter, estimating the accurate clutter spectrum with minor training support is achievable. Similar conclusions were also obtained in other works [28,29,30,31,32]. The deep learning technique provides a new idea for STAP design. Authors in [33] and [34] used deep convolution networks to offer the high-resolution profile, and simulations demonstrated that deep-learning STAP reduced the computational load compared with sparse representation STAP.
These previous studies for STAP provide some important insights into clutter suppression and target detection for airborne passive radar; however, the common assumption is that the reference signal is pure. In a realistic environment, this assumption is violated. To enable the target detection ability of STAP, the question is how to develop the MP clutter suppression approach before STAP.
In [9], Guo et al. used a blind equalization approach to construct the constraint optimization model and restore the pure DP signal based on the prior information. To overcome the dependence on prior knowledge, some sparse reconstruction algorithms were proposed for MP clutter suppression [8,10,11]. In our previous work [8,10], we derived the sparse reconstruction algorithms (SRAs), in which the sparse model of MP clutter was designed and several optimization algorithms based on the L1-norm were exploited to estimate the profile. In our previous work [11], under the framework of sparse models in SRAs, we used the local search algorithm (LSA) to select optimal atoms to construct the range-Doppler dictionary and address the off-grid issue. However, the presence of ICM disturbs the sparse model of MP clutter and corrupts the MP clutter suppression performance of SRA and LSA. To the best of our knowledge, very little literature addresses the ICM problem in MP clutter suppression. Although the ICM problem also arose in STAP [35,36,37,38], these methods only focused on the DP clutter cancellation. In this paper, we propose a novel MP clutter suppression method to mitigate ICM effects while also taking off-grids into account.

3. Signal Model

The system being considered is an airborne passive radar. We assume that the transmitter is a ground-based non-cooperative emitter and that the airborne platform moves in a straight line with constant altitude H R and velocity v R . The side-looking airborne passive radar geometry is shown in Figure 1. The measurement signals are received by a uniform linear array consisting of N elements. The reference antenna is directed toward the non-cooperative transmitter and receives the reference signal. The impure reference signal involves the DP signal, whose time bin is ld, and the NT MP signals, whose time bins are lp ( p = 1 , , N T ), respectively. Let the reference signal and measurement signals be segmented uniformly. Then, signals are equivalent to M pulses with a constant pulse repetition interval T r in a coherent integration time. L is the number of fast-time samples in T r . After using the segmented signals for matched filtering, the received data set comprises LMN samples. The radar data corresponding to a fast-time sample can be arranged as a M N × 1 space-time snapshot.
Considering an impure reference signal and ignoring the influence of range ambiguities, the clutter snapshot at the l-th range cell is composed of two parts. The first part is the DP clutter snapshot, that is, the matched result of the scatterer echoes whose time bin is (l + ld) with the DP signal. The DP clutter space-time snapshot, including ICM-CP, is given by:
x d c , l = i = 1 N c ξ l + l d , i v s ϑ i v t ω i ω d γ 1 = i = 1 N c ξ l + l d , i T v ω d , i , ϑ i ,
where N c is the number of statistically independent clutter patches; ξ l + l d , i = ε l + l d , i ε d , where ε d and ε l + l d , i are the complex amplitudes of the DP signal and the i-th clutter patch at the (l + ld)-th time bin, respectively; ϑ i = d cos φ i / λ and ω i = T r v R cos φ i / λ are spatial frequency and normalized Doppler frequency of the i-th clutter patch, respectively, where λ is the wavelength, d is the inter-channel spacing, and φ i is the spatial cone angle of the i-th clutter patch; ω d , i = ω i ω d indicates the relative (relative to the DP signal) normalized Doppler frequency, where ω d = T r ω d and ω d are the normalized Doppler frequency and Doppler frequency of the DP signal, respectively; γ 1 is the M × 1 temporal fluctuation vector corresponding to ICM-CP, which leads to the small pulse-to-pulse fluctuations in the scatterer echoes; the ICM-CP fluctuation matrix T = I N diag γ 1 ; v ω , ϑ = v s ϑ v t ω is the space-time steering vector with
v s ϑ = 1 , exp j 2 π ϑ , , exp j 2 π N 1 ϑ T ,
v t ω = 1 , exp j 2 π ω , , exp j 2 π M 1 ω T ,
denoting the spatial and temporal steering vectors, respectively.
The second part is the matched result of the scatterer echoes, whose time bin is (l + lp) ( p = 1 , , N T ) with the MP signals, that is, the MP clutter snapshot. In the presence of ICM-CP and ICM-MP, the MP clutter space-time snapshot is given by:
x p c , l = ε p i = 1 N c ξ l + l p , i v s ϑ i v t ω i ω p γ 1 γ 2 = ε p Γ i = 1 N c ξ l + l p , i T v ω p , i , ϑ i ,
where ε p = ε p / ε d represents the relative (relative to the DP signal) amplitude and ε p is the complex amplitude of the p-th MP signal [8]; ω p , i = ω i ω p is the relative (relative to the p-th MP signal) normalized Doppler frequency, where ω p = T r ω p and ω p are the normalized Doppler frequency and Doppler frequency corresponding to the p-th MP signal, respectively; γ 2 is the M × 1 temporal fluctuation vector corresponding to ICM-MP, which means the small pulse-to-pulse fluctuations exist in MP signals; the vector γ 2 is also called the ICM-MP fluctuation vector in what follows; the ICM-MP fluctuation matrix Γ = I N d i a g γ 2 .
By defining the modified steering vector f ω d ω p = 1 N × 1 v t ω d ω p , we obtain v ω p , i , ϑ i = v ω d , i , ϑ i f ω d ω p . Thus, Equation (4) is reformulated as [8]:
x p c , l = ε p Γ x d c , l + l p l d f ω p ω d = ε p Γ x l + l p l d f ω p ω d ε p Γ p ¯ = 1 N T ε p ¯ Γ x l + l p + l p ¯ 2 l d f ω ¯ p + ω p 2 ω d .
where the first term and second term of the second equation in Equation (5) are called the true MP clutter and the derived ghost MP clutter, respectively [8]. Thus, the MP clutter snapshot correlates with three parameters, which are the DP clutter snapshot, the ICM-MP fluctuation matrix, and the modified steering vector. Then, the received snapshot at the l-th range cell, consisting of the DP clutter, the MP clutter, and the receiver thermal noise, is expressed as
x l = x d c , l + p = 1 N T x p c , l + x n ,
where x n is the noise vector, which is assumed to be spatially and temporally uncorrelated. The complex amplitudes of clutter patches ξ l , i , i = 1 , 2 , , N c , are assumed to be independent and identically distributed random variables. Hence, the corresponding covariance matrix R x is given by:
R x = E x l x l H = R d + p = 1 N T R p + R n ,
where the DP clutter covariance matrix R d = E x d c , l x d c , l H ; the p-th MP clutter covariance matrix R p = E x p c , l x p c , l H ; R n = E x n x n H is the noise covariance matrix.
From the above analyses, it is evident that the presence of MP clutter severely affects the estimated covariance matrix. As a consequence, directly applying the STAP algorithms to suppress interference (DP clutter and MP clutter) requires more degrees of freedom than merely canceling DP clutter [8]. Additionally, the detection performance of STAP degrades seriously when targets fall within the MP clutter area [8].
The MP clutter cancellation method prior to STAP can overcome the influence of contaminated reference signals effectively. However, the presence of ICM-MP fluctuation disturbs the sparse model of MP clutter, and it transforms the reconstruction issue into a joint optimization problem. Consequently, both the existing SRA and LSA suffer from a performance loss in MP clutter estimation, thereby degrading the target detection ability of STAP. As a result, it is significant to devise an MP clutter suppression algorithm with ICM-MP fluctuation mitigation.

4. Proposed MP Clutter Suppression Algorithm Based on Joint Optimization

Here, the existing optimization problem in SRA and LSA is briefly reviewed. Subsequently, the motivations for the proposed joint optimization algorithm are presented, where a constraint to ICM-MP fluctuation is added based on the exiting optimization problem. Finally, the proposed joint optimization method is described.

4.1. Review of the Existing Optimization Problem in SRA and LSA

In the existing SRA and LSA, the MP clutter is related to the DP clutter snapshot and the modified steering vector in the absence of ICM-MP fluctuation. Thus, the whole temporal plane is discretized uniformly into Q grid points. The corresponding normalized Doppler frequency and modified steering vector are ω q and f ω q = 1 N × 1 v t ω q , q = 1 , 2 , , Q , respectively. Supposing the MP clutter can be expressed by D received snapshots, then the sparse model of MP clutter can be approximated by [8]:
p = 1 N T x p c , l = S l α ,
where the range-Doppler dictionary S l = S D , l F 1 , S D , l F 2 , , S D , l F Q ; the M N × D matrix S D , l = x l + 1 , x l + 2 , , x l + D , the M N × D matrix F q = 1 1 × D f ω q = 1 N × D v t ω q ; α indicates the D Q × 1 profile, and there exists a degree of sparsity in α [8].
By denoting S D , n = x n 1 1 × D and S n = S D , n F 1 , S D , n F 2 , , S D , n F Q , we obtain [8]:
E p = 1 N T x p c , l S l α 2 2 = C + E x l S l α 2 2 + 2 Re E x n H S l α = C + E x l S l α 2 2 + 2 Re E x n H S n α ,
where E x n H S n α = E x n H d i a g x n Φ α = ζ α , ζ = σ 2 1 1 × M N Φ , the modified steering matrix Φ = F 1 , F 2 , , F Q and σ 2 is the noise power [8], C is a constant independent of parameters S l and α .
By minimizing the output power, the MP clutter estimation is described as the following optimization problem [8]:
α ^ = arg min α E x l S l α 2 2 + 2 κ α 1 + 2 Re ζ α ,
where κ is a sparse constant parameter.
Based on the above optimization problem, SRA estimates the profile α to suppress MP clutter, and LSA modifies the range-doppler dictionary to overcome the off-grid problem. However, the presence of ICM-MP fluctuations significantly affects the MP clutter suppression performance of SRA and LSA.

4.2. The Motivation of the Proposed Joint Optimization Algorithm

If we define a l p , p ¯ = Γ x l + l p + l p ¯ 2 l d f ω ¯ p + ω p 2 ω d , then the p , p ¯ -th ghost MP clutter in Equation (5) is x l p , p ¯ = Γ ε p ε p ¯ a l p , p ¯ = Γ x l p , p ¯ . Thus, we get [35]:
E x l p , p ¯ x l p , p ¯ H E x l p , p ¯ x l p , p ¯ H = ε p 2 ε p ¯ 2 Γ ˜ 1 M N × M N R a ,
where R a = E a l p , p ¯ a l p , p ¯ H , Γ ˜ = I N × N Toeplitz ζ 0 , ζ 1 , , ζ M 1 and the temporal autocorrelation of the ICM fluctuations is ζ m = exp 8 π 2 σ v 2 T r 2 m 2 / λ 2 , m = 0 , 1 , , M 1 , λ is wavelength and σ v is the velocity standard deviation [12]. Since ε p 2 1 , ε p ¯ 2 1 and each element in the matrix Γ ˜ 1 M N × M N is far less than 1 in general [12], the effect of the covariance matrix E x l p , p ¯ x l p , p ¯ H in STAP can be equivalent to the effect of E x l p , p ¯ x l p , p ¯ H . Thus, Equation (5) can be approximated by:
x p c , l = ε p Γ x l + l p l d f ω p ω d ε p Γ p ¯ = 1 N T ε p ¯ x l + l p + l p ¯ 2 l d f ω ¯ p + ω p 2 ω d .
Therefore, the sparse model of MP clutter in the presence of ICM-MP becomes:
p = 1 N T x p c , l = Γ S l α = S l + Λ S l α ,
where Λ = Γ I M N = I N d i a g h and h = γ 2 1 M × 1 . Here, Λ S l and h , which are the functions of the ICM-MP fluctuation vector γ 2 , are defined as the ICM-MP errors calibration matrix and the ICM-MP errors calibration vector, respectively.
Substituting the above equation into Equation (9), we can get:
E p = 1 N T x p c , l Γ S l α 2 2 = C + E x l Γ S l α 2 2 + 2 Re E x n H Γ S n α = C + E x l Γ S l α 2 2 + 2 Re E x n H Γ d i a g x n Φ α ,
where
E x n H Γ d i a g x n Φ α = E x n H d i a g x n Γ Φ α = σ 2 1 1 × M N Γ Φ α .
We define ψ = σ 2 1 1 × M N Γ Φ , and we can rewrite the Equation (14) as
E p = 1 N T x p c , l Γ S l α 2 2 = C + E x l Γ S l α 2 2 + 2 Re ψ α .
By minimizing the ICM-MP error calibration matrix and the output power, the optimization problem can be formulated to estimate the profile vector α and calibration vector h [39]
α ^ , h ^ = arg min α , h J α , h = arg min α , h E x l Γ S l α 2 2 + 2 κ α 1 + 2 Re ψ α + Λ S l F 2 ,
where Γ = I N d i a g 1 M × 1 + h .
Based on the above analyses, the MP clutter suppression problem is equivalent to the estimation problem of the profile vector α and calibration vector h . To overcome this issue, the calibration vector h must first be fixed. Thus, the minimization of Equation (17) with a profile α becomes:
J 1 α = E x l Γ S l α 2 2 + 2 κ α 1 + 2 Re ψ α .
We compute the gradient terms of the above equation with respect to α , and let the gradient terms equal zero. Thus, we get:
α = R S 1 r S x κ sign α + ψ H ,
where R S = E S l H Γ H Γ S l ; r S x = E S l H Γ H x l and sign denotes a component-wise function.
On the other hand, by fixing the profile vector α , the estimation of h is transformed into the following optimization problem:
J 2 h = E x l Γ S l α 2 2 + 2 Re ψ α + Λ S l F 2 .
Since Λ = I N d i a g h and Γ = I N d i a g 1 M × 1 + h are the diagonal matrices, we get Λ S l F 2 = C 1 Λ F 2 = C 1 N h 2 2 where C 1 is a constant, and
Re ψ α = Re σ 2 1 1 × M N I N d i a g 1 M × 1 + h Φ α = Re σ 2 1 1 × M N Φ α + Re σ 2 1 1 × M N I N d i a g h Φ α = C 2 + Re σ 2 Φ α Τ 1 N × 1 I M h .
For simplicity, we denote ς = σ 2 Φ α Τ 1 N × 1 I M . Then, Equation (20) can be reformulated as:
J 2 h = E x l S l α I N d i a g h S l α 2 2 + 2 Re ς h + h 2 2 = E e l G l h 2 2 + 2 Re ς h + h 2 2 ,
where e l = x l S l α and G l = d i a g S l α 1 N × 1 I M . We obtain the gradient terms of Equation (22) with respect to the calibration vector:
J h = J 2 h h * = R G h r G e + ς H ,
where R G = E G l H G l + I M and r G e = E G l H e l . Let the above equation be equal to zero, and then we get:
h = R G 1 r G e ς H .
Now, expressions (19) and (24), derived from the optimization problem presented in Equation (17), are obtained to estimate the profile vector α and the calibration vector h . However, it is difficult to obtain the analytical solutions of α and h by directly solving these two expressions. This is because Equations (19) and (24) are not closed-form expressions for α and h , that is, the variables R S , r S x and ψ in Equation (19) are functions of h , and the variables R G , r G e and ς in Equation (24) are functions of α . Therefore, estimating the vectors α and h jointly provides an alternative solution.
The above development shows that the ICM-MP fluctuation vector/matrix affects the existing sparse model of MP clutter and optimization problems. Consequently, the MP clutter suppression performance is degraded in SRA and LSA, which will be presented in Section 5. On the other hand, it is concluded from reference [11] that off-grid issues are difficult to avoid in practice. Thus, it is necessary to develop a joint optimization algorithm with fluctuation error calibration and off-grid effect mitigation.

4.3. The Proposed Joint Optimization Algorithm

Here, the joint optimization algorithm for the ICM-MP fluctuation and profile estimation is derived. Subsequently, a two-step suppression method for MP clutter and DP clutter cancellation is presented, and the overall procedure of the proposed joint optimization algorithm is summarized.
In the joint optimization algorithm, two steps are performed in each external alternation: (a) the profile estimation step, where the local search technique presented in reference [11] is exploited to mitigate off-grid effects; and (b) the ICM-MP errors calibration step, where the calibration vector h is estimated.
In the profile estimation step of the k-th external alternation, the profile vector α is estimated with off-grid effect mitigation. In each external alternation, the optimal atoms representing MP clutter are selected by several iterative searches [11]. Each iterative search contains a global search and several local searches, and an optimal atom is determined when each iterative search is completed. For simplicity, let symbol k , m denote the m-th ( m = 1 , 2 , 3 , ) iterative search in the k-th external alternation, and let symbols k , m , 0 and k , m , n ( n = 1 , 2 , ) be the global search and the n-th local search of the k , m -th iterative search, respectively.
Initially, we discretize the temporal plane with a relatively large interval Δ ω , and build the following range-Doppler dictionary, termed the global dictionary.
S l = x l + 1 f ω 1 , , x l + D f ω 1 , , x l + 1 f ω Q , , x l + D f ω Q = g l , 1 , , g l , D , , g l , Q 1 D + 1 , , g l , D Q ,
where g l , η ( η = 1 , 2 , , D Q ) is the η -th global atom. Then, the global index in the k , m , 0 -th global search is picked up by [11]:
η k , m , 0 = arg max η 1 , 2 , , D Q E g l , η H Γ ^ H k 1 e ^ l k , m ,
where e ^ l k , m is the residual vector, which is initialized as e ^ l k , 1 = x l . Based on the selected global index, the range cell and normalized Doppler frequency corresponding to the optimal global atom g l , η k , m , 0 are denoted by τ ^ k , m and ω k , m , 0 , respectively.
Subsequently, we conducted several local searches. Motivated by the local search technique [11], the local dictionary in the k , m , n -th local search is expressed as:
S l k , m , n = x l + τ ^ k , m f ω k , m , n 1 Δ ω 2 n , x l + τ ^ k , m f ω k , m , n 1 , x l + τ ^ k , m f ω k , m , n 1 + Δ ω 2 n = g l , 1 k , m , n , g l , 2 k , m , n , g l , 3 k , m , n .
Correspondingly, the local index in the k , m , n -th local search is [11]:
η k , m , n = arg max η 1 , 2 , 3 a η k , m , n = arg max η 1 , 2 , 3 E g l , η H k , m , n Γ ^ H k 1 e ^ l k , m ,
where a η k , m , n denotes the local response. Thus, the optimal local atom is g l , η k , m , n and the corresponding normalized Doppler frequency is ω k , m , n .
When max a η k , m , n a η k , m , n 1 is less than the local threshold δ 1 , the local search terminates, which also means the k , m -th iterative search has been completed and an optimal atom is determined. Supposing the optimal normalized Doppler frequency is ω ^ k , m , then the optimal atom is g ^ l k m = x l + τ ^ k , m f ω ^ k , m . Let S ^ l and Φ ^ be the optimal dictionary and modified steering matrix, respectively. S ^ l and Φ ^ are initialized as S ^ l k , 0 = and Φ ^ l k , 0 = , and updated by S ^ l k , m = S ^ l k , m 1 , g ^ l k , m and Φ ^ k , m = Φ ^ k , m 1 , f ω ^ k , m , respectively. Therefore, based on Equation (19), the profile vector α at the k , m -th iterative search is calculated by using Γ ^ k 1 .
α ^ k , m = R ^ S 1 k , m r ^ S x k , m + ψ ^ H k , m ,
where R ^ S k , m = E S ^ l H k , m Γ ^ H k 1 Γ ^ k 1 S ^ l k , m ; r ^ S x k , m = E S ^ l H k , m Γ ^ H k 1 x l and ψ ^ k , m = σ 2 1 1 × M N Γ ^ k 1 Φ ^ k , m . Correspondingly, the residual vector is updated by e ^ l k , m + 1 = x l Γ ^ k 1 S ^ l k , m α ^ k , m .
The iterative search terminates until I k , m I k , m 1 is smaller than the search threshold δ 2 , where I k , m = S I N R k , m / S N R opt denotes the signal-to-interference- plus-noise-ratio (SINR) loss of the k , m -th iterative search, S I N R k , m and S N R opt are SINR of the k , m -th iterative search and the optimum output signal-to-noise ratio in an interference-free environment, respectively [12]. It also means the profile estimation step of the k-th external alternation has been completed, and all optimal atoms representing MP clutter have been selected. Consequently, we obtain the profile vector α ^ k in the k-th external alternation. The corresponding optimal dictionary and modified steering matrix are expressed as S ^ l k and Φ ^ k , respectively.
In the ICM-MP error calibration step of the k-th external alternation, the calibration vector h is calculated based on Equation (24)
h ^ k = R ^ G 1 k r ^ G e k ς ^ H k ,
where the covariance matrix R ^ G k = E G ^ l H k G ^ l k + I M ; the cross-correlation vector r ^ G e k = E G ^ l H k e ^ l k ; e ^ l k = x l S ^ l k α ^ k and G ^ l k = d i a g S ^ l k α ^ k 1 N × 1 I M ; and ς ^ k = σ 2 Φ ^ k α ^ k Τ 1 N × 1 I M .
To efficiently estimate h ^ k , let
B ^ k = diag S ^ l k α ^ k diag S ^ l k α ^ k = diag b ^ 1 k , b ^ 2 k , , b ^ N k ,
where the 1 × M vector b ^ i k is composed of i 1 M + 1 -th to i M -th diagonal elements of B ^ k . Then, we get:
G ^ l H k G ^ l k = 1 1 × N I M B ^ k 1 N × 1 I M = diag b ^ 1 k , diag b ^ 2 k , , diag b ^ N k I M , I M , , I M Τ = diag i = 1 N b ^ i k .
Thus, G ^ l H k G ^ l k is the diagonal matrix, and the diagonal matrix R ^ G k = d i a g R ^ k , 1 , R ^ k , 2 , , R ^ k , M . Then, Equation (30) is efficiently calculated by:
h ^ k = r ^ k , 1 / R ^ k , 1 , r ^ k , 2 / R ^ k , 2 , , r ^ k , M / R ^ k , M Τ ,
where r ^ k , i is the i-th element of r ^ G e k ς ^ H k . Correspondingly, the update of the ICM-MP fluctuation matrix is given by Γ ^ k = I N d i a g 1 M × 1 + h ^ k .
When the gain of the cost function expressed by Equation (17) is less than a predetermined alternation threshold δ 3 [39], the alternating step between Equations (29) and (33) terminates.
A summary of the proposed joint optimization algorithm is presented in Figure 2.
After obtaining the profile vector and the calibration vector, a two-step suppression process is performed. First, the MP clutter is canceled by utilizing the proposed joint optimization algorithm. By exploiting the estimated α ^ and h ^ , the output at the l-th range cell is:
y l = x l Γ ^ S ^ l α ^ .
Then, the STAP algorithm is used for DP clutter cancellation. Supposing the estimated ICM-CP fluctuation matrix T ^ is obtained by prior knowledge or the alternating direction method [40], then the snapshot after MP clutter suppression and ICM-CP calibration can be expressed as y l = T ^ 1 y l = T ^ 1 x l Γ ^ S ^ l α ^ . The STAP weight vector is calculated by:
w y = R ^ u 1 v ϖ t , ϑ t v H ϖ t , ϑ t R ^ u 1 v ϖ t , ϑ t ,
where R ^ u = 1 / K k = 1 , k l K + 1 y k y k H and v ϖ t , ϑ t is the space-time vector of the hypothetic target. Then, the scalar output is:
z l = w y H y l = w y H T ^ 1 x l Γ ^ S ^ l α ^ .

4.4. Analyses of Computational Complexity

In this subsection, we discuss the computational complexity of the proposed algorithm in terms of the number of additions and multiplications. In the k-th external alternation, L 2 + 1 M N D Q multiplications and L 2 D Q M N 1 additions are required to compute η k , m , 0 in global search, where L 2 denotes the number of samples used to compute mathematical expectation, 3 M N K n L 2 + 2 multiplications and 3 L 2 K n M N 1 additions are required to compute η k , m , n in local search, where K n is the local search number, ο m 3 + m 2 M N L 2 + 1 + M N 3 m L 2 + m + 1 multiplications and ο m 3 + m 2 M N L 2 L 2 + 1 + m M N 1 L 2 + 1 additions are required to compute α ^ k , m , 2 m M N multiplications and m M N additions are required to compute e ^ l k , m . Thus, in the step of updating α ^ , m = N 1 N 2 1 N 1 ο m 3 + m 2 + M N m 2 L 2 + 3 m L 2 + 3 m + 1   + N 2 M N L 2 D Q + 3 L 2 K n + D Q + 6 K n multiplications and N 2 L 2 M N 1 3 K n + D Q   + m = N 1 N 2 1 N 1 ο m 3 + L 2 M N 1 m 2 + m + m m + 2 M N 1 additions are required, where N 1 is the iterative search number, and N 2 is the number of updating S ^ l k , m . In the step of updating h ^ , L 2 M N N 1 + 1 multiplications and L 2 M N N 1 1 + M additions are required to compute R ^ G k , L 2 M N 2 M + 1 multiplications and L 2 M 2 M N + N 2 additions are required to compute r ^ G e k , and M + M N M + N 1 multiplications and M N M + N 1 1 additions are required to compute h ^ k . Thus, L 2 M N 2 M + N 1 + 2 + M + M N M + N 1 multiplications and L 2 M 2 M N + N N 1 + N 3   + M + M N M + N 1 1 additions are required. From the above analyses, the added ICM-MP error calibration step slightly increases the computational cost of the proposed algorithm.

5. Simulations and Performance Analyses

In this section, some examples with simulated data are used to verify the suppression performance of the proposed joint optimization algorithm. In the simulation, the transmitted digital video broadcasting signal is approximated by the Gaussian random signal [8]. The important parameters of the simulated scenarios are given in Table 1. The reference channel receives the DP signal and three MP signals. The normalized Doppler frequency of the DP signal is 0.5, and the DP signal to noise ratio is set to 60 dB where the noise power is normalized. The relative (relative to DP signal) time bins-relative (relative to DP signal) normalized Doppler frequencies-MP signal to DP signal ratios (MDRs) of three MP signals are (3, 0.23, −18 dB), (4, −0.08, −18 dB), and (9, −0.22, −25 dB), respectively. The velocity standard deviations of ICM-CP fluctuation and ICM-MP fluctuation are σ 1 , v = 0.5 and σ 2 , v = 0.7 , respectively. A total of 20 training samples are used in global and local searches to select atoms. Further, 200 training samples are used to estimate the covariance matrix. In addition, we set Q = 11 , local threshold δ 1 = 0.00001 and a search threshold δ 2 = 0.1 for the proposed joint optimization algorithm.

5.1. Selection of the Number of Snapshots D

First, we present the SINR loss performance of the proposed joint optimization algorithm versus the number of snapshots in the range-Doppler dictionary (the value of D). Figure 3 gives the simulation results. It is seen that the SINR loss performance depends on the value of D. As D increases, more atoms are used to construct a global dictionary and match MP clutter. Hence, the MP clutter estimation error is decreased and the SINR loss performance is enhanced. On the other hand, an increasing dimension in the range-Doppler dictionary leads to higher computational complexity. Therefore, the value of D is selected by making a trade-off between the computational cost and SINR loss performance. Based on the experimental results, D is set at 9.

5.2. Setting of the Alternation Threshold δ 3

Figure 4 gives the gain in performance in the cost function versus the external alternation number in the proposed method. It is seen that the gain in cost function tends to zero gradually with the increasing number of external alternations, which means the profile estimation and the ICM-MP error calibration are close to steady-state performance. However, the value of the cost function will not improve much when the external alternation number is too large and the computational load is heavy. Thus, for this reason and to keep the complexity low, we adopt δ 3 = 0.1 in the following experiments. In Figure 4, the alternation threshold is expressed as the red horizontal line. As a result, the number of external alternations equals eight.

5.3. Comparison with the Exiting Methods

In the third experiment, ICM-CP and ICM-MP are considered, and we compare the performances of the proposed algorithm with those of the existing others, i.e., STAP, LSA, and SRA. In STAP, the existing STAP algorithm is used to suppress DP and MP clutter simultaneously with the ICM-CP mitigation. For a fair and reasonable comparison, LSA, SRA, and the proposed algorithm use the same STAP algorithm to remove DP clutter. In the existing LSA and SRA, we separate the temporal plane into 11 cells, and D = 9 . For SRA, the forgetting factor is set to 0.95 and the constraint parameter is 70. For LSA, we set the local threshold at 0.0001 and the iteration threshold at 0.1, which are the same as the local threshold and search threshold in the proposed algorithm, respectively.
First, the distributions of optimal atoms selected by the proposed joint optimization algorithm and LSA are discussed. Table 2 gives the range-Doppler results of the selected atoms for different methods. The atoms corresponding to the real MP clutter component are called as real atoms, and the range-Doppler parameters of real atoms are also presented in Table 2. In the simulation, the number of iterative searches in LSA and the proposed algorithm are both six. Based on the simulated MDRs of three MP signals, the six real atoms are composed of three true atoms and three ghost atoms. Based on Equation (5), the range-Doppler parameters of a true atom and a ghost atom are denoted by l p l d , ω p ω d and l p + l p ¯ 2 l d , ω ¯ p + ω p 2 ω d , respectively. As can be seen from Table 2, the selected atoms of LSA, especially the sixth selected atom, deviate far from the real atoms. This is because the presence of ICM-MP fluctuations affects the sparse model and disturbs the selection criteria. Compared with LSA, the selected atoms in the proposed algorithm are closer to the real atoms owing to the added ICM-MP error calibration step. This also indicates that the proposed algorithm can eliminate the impact of ICM-MP fluctuation on local search techniques, providing an accurate estimation for the MP clutter profile.
Figure 5 presents the SINR loss performances of the proposed algorithm and three existing methods along the main beam direction. The normalized Doppler frequencies span from −0.5 to 0.5. From the simulation results, we find that STAP has nulls both in the DP clutter area (normalized Doppler frequencies of −0.5 and 0.5) and the MP clutter area (normalized Doppler frequencies of −0.23, 0.08, 0.22, −0.46, −0.15 and 0.16), while the other methods exhibit nulls only in the DP clutter area. This is because SRA, LSA, and the proposed approach effectively cancel the MP clutter before DP clutter suppression. As a consequence, the SINR loss performance in the non-DP clutter region reflects the MP clutter suppression capability. Furthermore, it is also shown from Figure 5 that the proposed algorithm provides a better SINR loss performance than LSA and SRA. This is because the proposed algorithm adds the ICM-MP error calibration step to mitigate its effect on MP profile estimation compared with LSA and addresses the off-grid and ICM problems compared with SRA.
Figure 6 presents the target detection results for different methods. Snapshots from the 992th to 1008th range cells are filtered in each method. It is found that STAP yields the worst target detection performance because the target located in the MP clutter region is also canceled. Moreover, SRA and LSA suffer from some performance loss in the presence of ICM-MP fluctuation when compared with the proposed algorithm. Therefore, it is concluded from the third experiment that the proposed joint optimization approach can provide satisfactory cancellation capability and detection performance in the presence of ICM-MP fluctuation and off-grid issues.
In the fourth experiment, we compare the performance of the proposed joint optimization algorithm with other existing methods in different ICM-MP scenarios. First, we consider three ICM-MP scenarios, namely no ICM-MP fluctuation, that is, σ 2 , v = 0 , σ 2 , v = 0.8 , and σ 2 , v = 1.6 . In the simulation, the MDR of three MP signals is −20 dB, and the remaining simulation parameters are the same as in the third experiment.
Figure 7 presents the SINR loss performance against the normalized Doppler frequency in different ICM-MP scenarios. From the simulation results, it is found that SRA, LSA, and the proposed algorithm outperform STAP regardless of the presence or absence of the ICM-MP fluctuation. Moreover, the SINR loss performance of SRA suffers a lot in these three cases. The reason is that the presence of ICM-MP fluctuations and off-grid problems disturbs the MP clutter estimation in SRA, leading to a severe SINR performance loss. Moreover, the SINR loss performance of LSA changes obviously with respect to ICM-MP fluctuations. Especially as the velocity standard deviation σ 2 , v increases, the performance of LSA becomes worse. This is because the accuracy of the selected atoms and estimated profile in LSA is sensitive to the ICM-MP fluctuation. Conversely, the proposed joint optimization algorithm always features acceptable SINR loss performance, whether the ICM-MP fluctuation is present or not.
In Figure 8, we present the corresponding detection results of the proposed joint optimization approach and the other three methods. The results illustrate that STAP yields substandard detection performance, and SRA no longer performs properly even in the absence of CIM-MP fluctuation. Additionally, as σ 2 , v increases, the detection performance of LSA improves. However, the proposed algorithm still provides satisfactory detection performance.
Then, we make a quantitative comparison between the proposed algorithm, SRA, and LSA with different ICM-MP fluctuations, i.e., σ 2 , v = 0.2 , 0.6, 1.0, 1.4, and 1.8. The SINR loss value at the target Doppler frequency considering different σ 2 , v is listed in Table 3. It is clear that the proposed joint optimization method can cope with the off-grid issues and is robust to the ICM-MP fluctuation.
In the fifth experiment, the SINR loss performances of different approaches with different MDR values are discussed. Figure 9 presents the SINR loss performances of different algorithms. In the simulation, we consider three MP signals having the same MDR values and set the MDR values as −22 dB and −17 dB, respectively. Other simulation parameters are the same as the third experiment.
It is seen from Figure 9 that the SINR loss performances of STAP and LSA vary obviously with MDR, and the SINR loss performance deteriorates with an increasing MDR value. The reason is that the high MDR causes deep nulls in the MP clutter area for STAP and residual MP clutter in profile estimation for LSA. Additionally, SRA provides substandard performance, and the performance degradation is due to the off-grid issue and the ICM problem. Conversely, the proposed joint optimization approach provides superior SINR loss performance in these two MDR scenarios.
Then, the proposed algorithms, SRA and LSA, are compared quantitatively. Table 4 represents the SINR loss values at the target Doppler frequency considering different MDR values, i.e., MDR = −15 dB, −19 dB, −23 dB, −27 dB, and −31 dB, respectively. It is found from Table 4 that the proposed approach provides enhanced SINR loss performance compared to SRA and LSA when the MDR value changes. Therefore, it is concluded that regardless of the MDR values, the proposed algorithm can resolve the off-grid issue and is robust to the ICM-MP fluctuation.
Finally, we evaluate the computational complexities of different algorithms. The computational complexity of the proposed algorithm is discussed in Section 4.4, and the computational complexity of SRA and LSA has been analyzed in references [10,11]. The simulation results are presented in Figure 10. The curves indicate that the proposed algorithm has lower computational complexity than SRA. This is because the inverse operation of the high-dimensional matrix is computationally intensive for SRA. In addition, the proposed algorithm has mildly higher complexities than LSA due to the dictionary updating and the added calibration step. However, the proposed algorithm outperforms LSA in MP clutter suppression, which was presented in the previous experiments. Therefore, it is concluded that the proposed algorithm has satisfactory suppression performance and moderate complexity in the presence of ICM-MP errors.

6. Conclusions

In the presence of ICM-MP fluctuations, the existing sparse model of MP clutter is disturbed, and the performance of the existing MP clutter suppression methods drops. In this paper, a MP clutter cancellation approach is presented to eliminate the influence of ICM-MP fluctuation. The proposed algorithm modifies the sparse model of MP clutter, derives a joint optimization problem, and alternately estimates the ICM-MP fluctuations and MP clutter profile with off-grid effects mitigation. The experimental results indicate the proposed joint optimization algorithm outperforms the existing STAP, LSA, and RSA in MP clutter cancellation when ICM-MP fluctuation exists. Moreover, the proposed approach can deal with the off-grid issue and is robust to the ICM-MP fluctuation, regardless of the MDR values.
By ignoring the influence of range ambiguities, the proposed algorithm can suppress MP clutter with ICM-MP calibration. However, the assumption may be violated in practice, and the performance of the proposed algorithm may degrade. Consequently, extending our model by considering range ambiguities is the research direction of our future work. In addition, possible future studies include robust MP clutter suppression in nonhomogeneous and nonstationary clutter environments for airborne passive radar.

Author Contributions

Formal analysis, Y.D. and Z.P.; methodology, Y.D. and Z.P.; validation, Y.D. and W.L.; writing—original draft, Y.D.; writing—review and editing, W.L. and D.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Hunan Provincial Natural Science Foundation of China, grant No. 2021JJ40029; Scientific Research Fund of the Hunan Provincial Education Department, grant No. 21B0710; and the 2021 Undergraduate Innovation and Entrepreneurship Training Program of Hunan Province, grant No. 3333.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing is not applicable to this article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Malanowski, M.; Rytel-Andrianik, R.; Kulpa, K.; Stasiak, K.; Ciesielski, M.; Kulpa, J. Geometric clutter analysis for airborne passive coherent location radar. IEEE Trans. Geosci. Remote Sens. 2022, 60, 5116314. [Google Scholar] [CrossRef]
  2. Tan, D.K.P.; Lesturgie, M.; Sun, H.; Lu, Y. Space–time interference analysis and suppression for airborne passive radar using transmissions of opportunity. IET Radar Sonar Navig. 2014, 8, 142–152. [Google Scholar] [CrossRef]
  3. Wojaczek, P.; Colone, F.; Cristallini, D.; Lombardo, P. Reciprocal-Filter-Based STAP for Passive Radar on Moving Platforms. IEEE Trans. Aerosp. Electron. Syst. 2019, 55, 967–988. [Google Scholar] [CrossRef]
  4. Rosenberg, L.; Duk, V. Land clutter statistics from an airborne passive bistatic radar. IEEE Trans. Geosci. Remote Sens. 2022, 60, 5104009. [Google Scholar] [CrossRef]
  5. Sui, J.; Wang, J.; Zuo, L.; Gao, J. Cascaded least square algorithm for strong clutter removal in airborne passive radar. IEEE Trans. Aerosp. Electron. Syst. 2022, 58, 679–696. [Google Scholar] [CrossRef]
  6. Brown, J.; Woodbridge, K.; Griffiths, H.; Stove, A.; Watts, S. Passive bistatic radar experiments from an airborne platform. IEEE Aerosp. Electron. Syst. Mag. 2012, 27, 50–55. [Google Scholar] [CrossRef]
  7. Sui, J.; Wang, J.; Zuo, L.; Gao, J. Multistage least squares algorithms for clutter suppression in airborne passive radar based on subband operation. IEEE Trans. Aerosp. Electron. Syst. 2023, 59, 1893–1909. [Google Scholar] [CrossRef]
  8. Deng, Y.; Wang, J.; Luo, Z.; Guo, S. Cascaded suppression method for airborne passive radar with contaminated reference signal. IEEE Access 2019, 7, 50317–50329. [Google Scholar] [CrossRef]
  9. Guo, S.; Wang, J.; Ma, H.; Wang, J. Modified blind equalization algorithm based on cyclostationarity for contaminated referent signal in airborne PBR. Sensors 2020, 20, 788. [Google Scholar] [CrossRef]
  10. Deng, Y.; Zhang, S.; Zhu, Q.; Zhang, L.; Li, W. Clutter suppression methods based on reduced-dimension transformation for airborne passive radar with impure reference signals. J. Appl. Remote Sens. 2021, 15, 016514. [Google Scholar] [CrossRef]
  11. Deng, Y.; Li, W.; Zhang, S.; Wang, F.; Xiao, W.; Cui, Z. Clutter suppression method for off-grid effects mitigation in airborne passive radars with contaminated reference signals. Sensors 2021, 21, 6339. [Google Scholar] [CrossRef]
  12. Ward, J. Space-Time Adaptive Processing for Airborne Radar. In Proceedings of the IEE Colloquium on Space-Time Adaptive Processing, London, UK, 6 April 1998. [Google Scholar]
  13. Klemm, R. Principles of Space-Time Adaptive Processing, 3rd ed.; Institute of Electical Engineering: London, UK, 2006. [Google Scholar]
  14. Mevin, W.L.; Showman, G.A. An approach to knowledge-aided covariance estimation. IEEE Trans. Aerosp. Electron. Syst. 2006, 42, 1021–1042. [Google Scholar] [CrossRef]
  15. Liu, M.; Zou, L.; Yu, X.; Zhou, Y.; Wang, X.; Tang, B. Knowledge aided covariance matrix estimation via gaussian Kernel function for airborne SR-STAP. IEEE Access 2020, 8, 5970–5978. [Google Scholar] [CrossRef]
  16. Riedl, M.; Potter, L.C. Multimodel shrinkage for knowledge-aided space-time adaptive processing. IEEE Trans. Aerosp. Electron. Syst. 2018, 54, 2601–2610. [Google Scholar] [CrossRef]
  17. Hu, J.; Li, J.; Li, H.; Li, K.; Liang, J. A novel covariance matrix estimation via cyclic characteristic for STAP. IEEE Geosci. Remote Sens. Lett. 2020, 17, 1871–1875. [Google Scholar] [CrossRef]
  18. Liu, C.; Wang, T.; Zhang, S.; Ren, B. A clutter suppression algorithm via weighted l2−norm penalty for airborne radar. IEEE Signal Process. Lett. 2022, 29, 1522–1525. [Google Scholar] [CrossRef]
  19. Zhou, Y.; Wang, L.; Chen, X.; Wen, C.; Jiang, B.; Fang, D. An improving EFA for clutter suppression by using the persymmetric covariance matrix estimation. Circuits Syst. Signal Process. 2018, 37, 4136–4149. [Google Scholar] [CrossRef]
  20. Cui, N.; Duan, K.; Xing, K.; Yu, Z. Beam-space reduced-dimension 3D-STAP for nonside-looking airborne radar. IEEE Geosci. Remote Sens. Lett. 2022, 19, 3506505. [Google Scholar] [CrossRef]
  21. Shen, S.; Tang, L.; Nie, X.; Bai, Y.; Zhang, X.; Li, P. Robust space time adaptive processing methods for synthetic aperture radar. Appl. Sci. 2020, 10, 3609. [Google Scholar] [CrossRef]
  22. Shen, M.; Zhu, D.; Zhu, Z. Reduced-rank space-time adaptive processing using a modified projection approximation subspace tracking deflation approach. IET Radar Sonar Navig. 2009, 3, 93–100. [Google Scholar] [CrossRef]
  23. Huang, P.; Zou, Z.; Xia, X.; Liu, X.; Liao, G. A novel dimension-reduced space–time adaptive processing algorithm for spaceborne multichannel surveillance radar systems based on spatial–temporal 2-D sliding window. IEEE Trans. Geosci. Remote Sens. 2022, 60, 5109721. [Google Scholar] [CrossRef]
  24. Shi, J.; Xie, L.; Cheng, Z.; He, Z.; Zhang, W. Angle-Doppler channel selection method for reduced-dimension STAP based on sequential convex programming. IEEE Comm. Lett. 2021, 25, 3080–3084. [Google Scholar] [CrossRef]
  25. Chen, W.; Xie, W.; Wang, Y. Short-range clutter suppression for airborne radar using sparse recovery and orthogonal projection. IEEE Geosci. Remote Sens. Lett. 2022, 19, 3500605. [Google Scholar] [CrossRef]
  26. Xie, L.; He, Z.; Tong, J.; Zhang, W. A recursive angle-Doppler channel selection method for reduced-dimension space-time adaptive processing. IEEE Trans. Aerosp. Electron. Syst. 2020, 56, 3985–4000. [Google Scholar] [CrossRef]
  27. Yang, Z.; de Lamare, R.C.; Li, X. L1-regularized STAP algorithms with a generalized sidelobe canceler architecture for airborne radar. IEEE Trans. Signal Process. 2012, 60, 674–686. [Google Scholar] [CrossRef]
  28. Wang, X.; Yang, Z.; Huang, J.; de Lamare, R.C. Robust two-stage reduced-dimension sparsity-aware STAP for airborne radar with coprime arrays. IEEE Trans. Signal Process. 2020, 68, 81–96. [Google Scholar] [CrossRef]
  29. Li, X.; Yang, X.; Wang, Y.; Duan, K. Gridless sparse clutter nulling STAP based on particle swarm optimization. IEEE Geosci. Remote Sens. Lett. 2022, 19, 4023205. [Google Scholar] [CrossRef]
  30. 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. 2020, 56, 785–795. [Google Scholar] [CrossRef]
  31. Hu, Z.; Wang, W.; Dong, F. A low-complexity MIMO radar STAP strategy for efficient sea clutter suppression. IEEE Geosci. Remote Sens. Lett. 2022, 19, 4024105. [Google Scholar] [CrossRef]
  32. Cui, N.; Xing, K.; Yu, Z.; Duan, K. Tensor-based sparse recovery space-time adaptive processing for large size data clutter suppression in airborne radar. IEEE Trans. Aerosp. Electron. Syst. 2023, 59, 907–922. [Google Scholar] [CrossRef]
  33. Duan, K.; Chen, H.; Xie, W.; Wang, Y. Deep learning for high-resolution estimation of clutter angle-Doppler spectrum in STAP. IET Radar Sonar Navig. 2022, 16, 193–207. [Google Scholar] [CrossRef]
  34. Zou, B.; Wang, X.; Feng, W.; Zhu, H.; Lu, F. DU-CG-STAP method based on sparse recovery and unsupervised learning for airborne radar clutter suppression. Remote Sens. 2022, 14, 3472. [Google Scholar] [CrossRef]
  35. Yang, Z.; Qin, Y.; de Lamare, R.C.; Wang, H.; Li, X. Sparsity-based direct data domain space-time adaptive processing with intrinsic clutter motion. Circuits Syst. Signal Process. 2017, 36, 219–246. [Google Scholar] [CrossRef]
  36. Yang, Z.; Li, X.; Wang, H.; Fa, R. Knowledge-aided STAP with sparse-recovery by exploiting spatio-temporal sparsity. IET Signal Process. 2016, 10, 150–161. [Google Scholar] [CrossRef]
  37. Gu, Y.; Wu, J.; Fang, Y.; Zhang, L.; Zhang, Q. End-to-end moving target indication for airborne radar using deep learning. Remote Sens. 2022, 14, 5354. [Google Scholar] [CrossRef]
  38. Duan, K.; Liu, W.; Duan, G.; Wang, Y. 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]
  39. Ma, Z.; Liu, Y.; Meng, H.; Wang, X. Sparse recovery-based space-time adaptive processing with array error self-calibration. Electron. Lett. 2014, 50, 952–954. [Google Scholar] [CrossRef]
  40. Yang, Z.; de Lamare, R.C.; Liu, W. Sparsity-based STAP using alternating direction method with gain/phase errors. IEEE Trans. Aerosp. Electron. Syst. 2017, 53, 2756–2768. [Google Scholar]
Figure 1. Airborne passive radar geometry.
Figure 1. Airborne passive radar geometry.
Applsci 13 05653 g001
Figure 2. Flow chart of the proposed joint optimization algorithm.
Figure 2. Flow chart of the proposed joint optimization algorithm.
Applsci 13 05653 g002
Figure 3. SINR loss performance of the proposed approach versus the value of D.
Figure 3. SINR loss performance of the proposed approach versus the value of D.
Applsci 13 05653 g003
Figure 4. Gain performance in cost function versus external alternation number.
Figure 4. Gain performance in cost function versus external alternation number.
Applsci 13 05653 g004
Figure 5. SINR loss performances of different methods.
Figure 5. SINR loss performances of different methods.
Applsci 13 05653 g005
Figure 6. Target detection performance versus range cell for different methods.
Figure 6. Target detection performance versus range cell for different methods.
Applsci 13 05653 g006
Figure 7. Impacts of ICM-MP fluctuation on SINR loss performance of the proposed algorithm and the existing methods. (a) σ 2 , v = 0 ; (b) σ 2 , v = 0.8 ; (c) σ 2 , v = 1.6 .
Figure 7. Impacts of ICM-MP fluctuation on SINR loss performance of the proposed algorithm and the existing methods. (a) σ 2 , v = 0 ; (b) σ 2 , v = 0.8 ; (c) σ 2 , v = 1.6 .
Applsci 13 05653 g007aApplsci 13 05653 g007b
Figure 8. Impacts of ICM-MP fluctuation on target detection performance of the proposed algorithm and the existing methods. (a) σ 2 , v = 0 ; (b) σ 2 , v = 0.8 ; (c) σ 2 , v = 1.6 .
Figure 8. Impacts of ICM-MP fluctuation on target detection performance of the proposed algorithm and the existing methods. (a) σ 2 , v = 0 ; (b) σ 2 , v = 0.8 ; (c) σ 2 , v = 1.6 .
Applsci 13 05653 g008
Figure 9. Impacts of MDR on SINR loss performance of the proposed algorithm and the existing methods. (a) MDR = −22 dB; (b) MDR = −17 dB.
Figure 9. Impacts of MDR on SINR loss performance of the proposed algorithm and the existing methods. (a) MDR = −22 dB; (b) MDR = −17 dB.
Applsci 13 05653 g009
Figure 10. Computational complexities of different algorithms. (a) number of multiplications; (b) number of additions.
Figure 10. Computational complexities of different algorithms. (a) number of multiplications; (b) number of additions.
Applsci 13 05653 g010
Table 1. Simulation parameters of an airborne passive radar system.
Table 1. Simulation parameters of an airborne passive radar system.
ParameterValue
Number of spatial elements10
Number of temporal pulses10
Antenna arraySide-looking array
inter-channel spacing λ / 2
Pulse repetition interval2.5 ms
Platform velocity100 m/s
Number of clutter patches in each range cell181
Range cell index of target1000
Normalized Doppler frequency of target−0.28
Signal-to-noise ratio of target0 dB
Table 2. Range-Doppler results of selected atoms for different methods.
Table 2. Range-Doppler results of selected atoms for different methods.
IndexTrue 1True 2True 3Ghost 1Ghost 2Ghost 3
Real atoms[3, 0.23][4, −0.08][9, −0.22][6, 0.46][7, 0.15][8, −0.16]
LSA[3, 0.234][4, −0.0795][9, −0.2164][6, 0.4738][7, 0.1773][1, 0.05]
Proposed method[3, 0.2328][4, −0.0805][9, −0.218][6, 0.4727][7, 0.1758][8, −0.1797]
Table 3. The SINR loss value at the target Doppler frequency with different σ 2 , v .
Table 3. The SINR loss value at the target Doppler frequency with different σ 2 , v .
σ 2 , v 0.20.61.01.41.8
SRA (dB)−5.046−5.032−5.097−5.104−5.096
LSA (dB)−4.451−4.827−5.115−5.481−5.668
Proposed algorithm (dB)−4.278−4.321−4.318−4.432−4.417
Table 4. The SINR loss value at the target Doppler frequency with different MDR values.
Table 4. The SINR loss value at the target Doppler frequency with different MDR values.
MDR Value (dB)−15−19−23−27−31
SRA (dB)−5.321−5.032−4.957−4.915−4.876
LSA (dB)−5.228−4.867−4.565−4.348−4.084
Proposed algorithm (dB)−4.656−4.444−4.268−4.032−4.081
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Deng, Y.; Pei, Z.; Li, W.; Jiang, D. Clutter Suppression Algorithm with Joint Intrinsic Clutter Motion Errors Calibration and Off-Grid Effects Mitigation in Airborne Passive Radars. Appl. Sci. 2023, 13, 5653. https://doi.org/10.3390/app13095653

AMA Style

Deng Y, Pei Z, Li W, Jiang D. Clutter Suppression Algorithm with Joint Intrinsic Clutter Motion Errors Calibration and Off-Grid Effects Mitigation in Airborne Passive Radars. Applied Sciences. 2023; 13(9):5653. https://doi.org/10.3390/app13095653

Chicago/Turabian Style

Deng, Yaqi, Zhengwang Pei, Wenguo Li, and Dongchu Jiang. 2023. "Clutter Suppression Algorithm with Joint Intrinsic Clutter Motion Errors Calibration and Off-Grid Effects Mitigation in Airborne Passive Radars" Applied Sciences 13, no. 9: 5653. https://doi.org/10.3390/app13095653

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