Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Efficiency and Security Evaluation of Lightweight Cryptographic Algorithms for Resource-Constrained IoT Devices
Previous Article in Journal
Channel-Blind Joint Source–Channel Coding for Wireless Image Transmission
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Fast Iterative STAP Method with a Coprime Sampling Structure

1
School of Aeronautics and Astronautics, University of Electronic Science and Technology of China, Chengdu 611731, China
2
Chengdu Textile College, Chengdu 611731, China
*
Author to whom correspondence should be addressed.
Sensors 2024, 24(12), 4007; https://doi.org/10.3390/s24124007
Submission received: 22 April 2024 / Revised: 4 June 2024 / Accepted: 18 June 2024 / Published: 20 June 2024
(This article belongs to the Section Intelligent Sensors)

Abstract

:
In space-time adaptive processing (STAP), the coprime sampling structure can obtain better clutter suppression capabilities at a lower hardware cost than the uniform linear sampling structure. However, in practical applications, the performance of the algorithm is often limited by the number of training samples. To solve this problem, this paper proposes a fast iterative coprime STAP algorithm based on truncated kernel norm minimization (TKNM). This method establishes a virtual clutter covariance matrix (CCM), introduces truncated kernel norm regularization technology to ensure the low rank of the CCM, and transforms the non-convex problem into a convex optimization problem. Finally, a fast iterative solution method based on the alternating direction method is presented. The effectiveness and accuracy of the proposed algorithm are verified through simulation experiments.

1. Introduction

In airborne radar systems, space-time adaptive processing (STAP) is crucial for clutter suppression and moving target detection [1,2,3,4,5,6,7,8,9,10,11,12,13]. However, the performance of the STAP algorithm is largely affected by the accuracy of clutter covariance matrix (CCM) estimation. In practical applications, radar systems face complex clutter environments generated by factors such as terrain, buildings, and weather conditions. These clutters not only vary significantly in intensity but also possess distinct spatial and temporal characteristics. The received radar signal can be expressed as x = s + c + n , where s represents the target echo signal, c represents the clutter signal, and n represents the noise signal. The goal of the STAP algorithm is to extract the target echo signal s from the received signal x while suppressing the clutter signal c and noise signal n . The clutter covariance matrix is a crucial parameter that describes the statistical properties of the clutter signal, including its spatial and temporal correlations. Ideally, the clutter covariance matrix can be estimated using independently and identically distributed training samples from adjacent cells (without targets). Typically, the CCM is estimated using independent and identically distributed training samples (without targets) of the neighboring test cells. But, in a heterogeneous environment, it is difficult to meet this requirement. In addition, in order to ensure that the signal-to-interference-to-noise ratio (SINR) loss does not exceed 3 dB, the number of training samples should be at least twice the system degrees of freedom (DOF) [14]. However, the full vector filters require significant computational complexity and storage space. Therefore, algorithms with low samples and low computational complexity are particularly important in practical applications.
For decades, researchers have proposed various algorithms to estimate the CCM with limited training samples. These algorithms include reduced dimension [14,15], reduction-rank methods [16,17,18,19], and persymmetry methods [20,21], which can reduce the need for training samples. In addition, the direct data domain method [22] only uses test data and avoids dependence on training samples. However, these methods usually sacrifice system DOF, resulting in performance degradation. In recent years, the emergence of sparse arrays has provided new ideas for solving these problems. Through difference operations, sparse arrays can obtain virtual arrays similar to physical arrays, allowing the system to have more virtual sensors than physical sensors. In addition, the larger element spacing in the sparse array also significantly reduces the mutual coupling effect of the array, further improving the performance and stability of the system.
Common sparse arrays include minimal redundancy arrays [23,24], nested arrays [5,25,26,27], and coprime arrays [5,28,29,30,31]. Minimally redundant arrays provide the largest contiguous array for a given number of sensors, but lack closed expressions for sensor locations and DOF. Compared with nested arrays, coprime arrays have relatively weak mutual couple effects due to the large spacing between array sensors. However, a large number of holes appear in the difference coarray, which affects the full utilization of virtual elements. Due to the low rank of the CCM, it can be estimated by the compressed sensing algorithm. The literature [32] uses the symmetric structure [33] and its positive definiteness [34] of the covariance matrix to construct the kernel norm and proposes a sparse sensing algorithm that uses the cyclic rank minimization method to reconstruct the covariance matrix. Meanwhile, the literature [35] proposed an array interpolation algorithm using atomic norm minimization. However, the research results in [36] show that the rank kernel norm relaxation cannot guarantee the low rank of the CCM of the multi-level Toplitz matrix. In [37], Cai et al. demonstrated that kernel norm minimization based on low-rank matrix estimation problems can be solved through singular value contraction operations. In addition, Lin et al. [38] proposed an exact iterative method called augmented Lagrange multiplier to solve this problem. Although these algorithms perform well in terms of estimation accuracy, a large number of sensors are required to achieve the performance of these algorithms in practical applications. It has relatively high computational complexity, which poses challenges to the real-time performance and efficiency of the system. In order to overcome these problems, this paper proposes a fast iterative coprime STAP algorithm (TNNM-FIC-STAP) based on truncated kernel norm minimization.
First, a detailed analysis of the coprime sampling structure of the STAP is presented, along with the corresponding CCM determined in the virtual domain. Then, the concept of truncated kernel norm is introduced, and a fast iterative coprime STAP algorithm based on this theory is proposed. The final solution comes from the convex relaxation process. Finally, the effectiveness of the algorithm is verified through numerous simulation experiments.
The content of this article is arranged as follows: Section 2 introduces the coprime sampling structure model, difference operation, and filter weight vector estimation process. In Section 3, the truncated nuclear norm and its fast iteration process are explained. In Section 4, simulation results are used to verify the performance advantages of the algorithm. Finally, Section 5 summarizes the main contents of this article and presents conclusions. The list of abbreviations used in the text is presented in Table 1.

2. Coprime Sampling Structured STAP

2.1. Coprime Sampling Structured Model

Assuming a side-looking airborne phased array radar with a coprime sampling structure for both the array and pulse repetition interval (PRI), which emits M pulses in a coherent processing interval (CPI) and has N array sensors, as shown in Figure 1. The array consists of two sub-uniform linear arrays (ULAs). The denser sub-ULA contains N 2 array elements at position { N 1 p 2 d , 0 p 2 N 2 1 } , the sparser sub-ULA contains 2 N 1 1 array sensors at position { N 2 p 1 d , 1 p 1 2 N 1 1 } , where N 1 and N 2 are mutually prime integers satisfying N 1 < N 2 , d is the minimum distance between array sensors, as shown in Figure 1a. Similar to the coprime array structure, Figure 1b illustrates the coprime PRI structure. The impulse consists of two uniform sub-impulse, containing M 2 and 2 M 1 1 impulses, respectively. The position coordinates of the sub-impulses are { M 1 q 2 T r , 0 q 2 M 2 1 } and { M 2 q 1 T r , 1 q 1 2 M 1 1 } , where M 1 and M 2 are mutually prime integers satisfying M 1 < M 2 . The minimum PRI is T r . Without the ranger ambiguity, the space-time snapshot received from a range bin can be expressed as
x = a t v ( φ t , f t ) + x u
where a t represents the target complex amplitude, and v ( φ t , f t ) = v ( φ t ) v ( f t ) is the target space-time steering vector. And the spatial steering vector is
v ( φ t ) = [ 1 , e 2 π j n 1 φ t , , e 2 π j n N 1 φ t ] T
and the time steering vector is
v ( f t ) = [ 1 , e 2 π j m 1 f t , , e 2 π j m M 1 f t ] T
where φ t = d cos ( θ ) / λ , f t = 2 v r T r cos ( θ ) / λ , λ is the radar’s wavelength, v r is the radar velocity, and θ is the spatial cone angle where the target is located. Assuming that N c is the number of independent clutter patches representing in a range ring, then the clutter plus noise data x u is
x u = i = 1 N c a c , i v ( φ c , i , f c , i ) + n = i = 1 N c a c , i v ( φ c , i ) v ( f c , i ) + n
The normalized spatial and time frequencies of the i th clutter patch are φ c , i and f c , i , respectively. a c , i is its complex amplitude, n is the noise vector, and the space steering vector is
v ( φ c , i ) = [ 1 , e 2 π j n 1 φ c , i , , e 2 π j n N 1 φ c , i ] T .
The time steering vector is
v ( f c , i ) = [ 1 , e 2 π j m 1 f c , i , , e 2 π j m M 1 f c , i ] T
and the corresponding space-time steering vector is
v ( φ c , i , f c , i ) = v ( φ c , i ) v ( f c , i ) = [ 1 e 2 π j n 1 φ c , i e 2 π j n N 1 φ c , i ] [ 1 e 2 π j m 1 f c , i e 2 π j m M 1 f c , i ] = [ v 0 , i v 1 , i v N M 1 , i ]
where v l M + r 1 , i = e 2 π j ( n l φ c , i + m r 1 f c , i ) , l = 0 , , N 1 , r = 1 , , M , i = 1 , , N c . From formula (4), the clutter plus noise covariance matrix (CNCM) R u can be calculated by
R u = E [ x u x u H ] = i = 1 N c E ( | a c , i | 2 ) v ( φ c , i , f c , i ) v H ( φ c , i , f c , i ) + Q N M = V P V H + Q N M = R c + Q N M ,
where V = [ v ( f c , 1 , φ c , 1 ) , v ( f c , 2 , φ c , 2 ) , , v ( f c , N c , φ c , N c ) ] is the clutter space-time steering matrix, P = d i a g ( [ p 1 , p 2 , , p N c ] T ) ,   p k = E ( | a c , k | 2 ) represents the clutter power matrix, the CCM can be expressed as R c , and Q N M = d i a g ( σ 1 2 , σ 2 2 , , σ M N 2 ) is the noise covariance matrix. The background noise is the Gaussian white noise when the ideal condition Q N M = σ n 2 I M N is satisfied.

2.2. Difference Operation

To clearly explain the virtual array and impulse train, define the difference coarray D of the array A as
D = { d n 1 d n 2 | d n 1 , d n 2 A }
where d n 1 is the position of the n 1 array sensor. As shown in Figure 2a, the de-duplication set D l from set D is
D l = { ( 2 N 1 1 ) N 2 , ( 2 N 1 2 ) N 2 , , ( 2 N 1 1 ) N 2 } .
Holes in the arrays appear after ( N 1 N 2 + N 1 1 ) and ( N 1 N 2 + N 1 1 ) , and two discontinuous difference sub-arrays are symmetric. For a sparse impulse train , the difference operations set is
= { q m 1 q m 2 | q m 1 , q m 2 }
with its de-duplication set l
l = { ( 2 M 1 1 ) M 2 , ( 2 M 1 2 ) M 2 , , ( 2 M 1 1 ) M 2 }
As shown in Figure 2b, the structure and properties of impulses are similar to those of different arrays. To enlarge the DOF in the space and time domains, virtual array sensors and impulses can be increased by filling the holes.
A filled circle stands for physical sensors and an empty circle represents empty space. A filled rectangle stands for transmitting pulses and an empty rectangle represents empty space.

2.3. The STAP Method

Equation (8) can be written as
R u = R c + Q N M = [ R 0 , 0 R 1 , 0 R N M 1 , 0 R 0 , 1 R 1 , 1 R N M 1 , 1 R 0 , N M 1 R 1 , N M 1 R N M 1 , N M 1 ] + [ σ 1 2 0 0 0 σ 2 2 0 0 0 σ M N 2 ]
where
R ( l 1 M + r 1 1 ) , ( l 2 M + r 2 1 ) = k = 1 N c p k e 2 π j [ ( n l 1 n l 2 ) φ c , k + ( m r 1 1 m r 2 1 ) f c , k ]
and l 1 , l 2 = 0 , , N 1 ,   r 1 , r 2 = 1 , , M .
By comparing Equations (10), (12) and (13), it can be seen that the elements in the CCM correspond to the virtual sensors and the virtual impulse positions. To downplay the effects of the background noise on the covariance matrix elements, the covariance matrix elements corresponding to the same virtual array sensor and the virtual pulse are averaged. Let the i th element in R u corresponding to the d n virtual sensor and the q m virtual pulse be R u i ( d n , q m ) , then the mean value is
R ¯ u ( d n , q m ) = 1 w ( d n , q m ) i = 1 w ( d n , q m ) R u i ( d n , q m )
where w ( d n , q m ) denotes the number of times the array difference d n and impulse difference q m are repeated at the same time. Now, R u is
R u = R 1 + R 2  
where
[ R 1 ] i , j = { [ R u ] i , j , i j 0 , i = j
R 2 = d i a g ( [ R u ] 1 , 1 , [ R u ] 2 , 2 , , [ R u ] M N , M N )
and [ R u ] i , i = R i , i + σ i 2 .
To reduce the impact of non-uniform noise on the covariance matrix, all σ i 2 are replaced with σ min 2 = MIN ( σ i 2 ) to whiten it. In this way, the influence of non-uniform noise can be suppressed and the signal-to-noise ratio (SNR) can be improved. At this time, the estimated value of R 2 is denoted as R ^ 2 , with its corresponding covariance matrix denoted as R ^ u .
To increase the number of virtual array sensors and impulses, all holes are filled in the difference coarray, and the numbers of virtual sensors and virtual impulses are now 4 N 1 N 2 2 N 2 + 1 and 4 M 1 M 2 2 M 2 + 1 , respectively. Extending the covariance matrix R ^ u with zero padding results in an ( 4 N 1 N 2 2 N 2 + 1 ) × ( 4 M 1 M 2 2 M 2 + 1 ) matrix R E by the existing elements of wave range difference and impulse difference. It is assumed that R E = R 1 E + R 2 E , where R 1 E and R 2 E are the corresponding recovery matrices for R 1 and R 2   , respectively.

3. Fast Iterative Coprime STAP Method

3.1. Truncated Kernel Norm

The truncated kernel norm of the matrix R m × n ( m n ) is defined as the sum of the smaller n r singular values. Namely,
R r = i = 1 n σ i i = 1 r σ i = i = r + 1 n σ i
where σ i is the i th singular value. The singular value decomposition (SVD) of the matrix R is R = U Σ V T , U = [ u 1 , u 2 , u m ] , V = [ v 1 , v 2 , v n ] , Σ = diag ( σ 1 , σ 2 , , σ n ) . For R m × n ( m n ) and r < n , we have
t r ( X R Y T ) i = 1 r σ i , r N .
Here, X r × m and Y r × n satisfy X X T = I ,   Y Y T = I . When X = [ u 1 , u 2 , u r ] T = U r T and Y = [ v 1 , v 2 , v r ] T = V r T , the maximum value of t r ( X R Y T ) is obtained by
max   t r ( X R Y T ) = t r ( U r T R V r ) = i = 1 r σ i .

3.2. TNNM-STAP Method

From Equation (19), we have
min R 1 E , R 2 E   R 1 E r + λ R 2 E   1 ,   s . t .   R E = R 1 E + R 2 E
where R 1 E r is non-convex that cannot be directly solved, but it can be transformed into [39].
min R 1 E , R 2 E   R 1 E * max X X T = I , Y Y T = I t r ( X R 1 E Y T ) + λ R 2 E   1 ,   s . t .   R E = R 1 E + R 2 E
The first r larger singular values are assigned a weight of 0 and the remaining n-r smaller ones are assigned a weight of 1. This preserves the main component of the data and removes the noise. The tr ( X R 1 E Y T ) reaches its maximum value when X = U r T and Y = V r T , so its low-rank matrix can be recovered in two steps. In the k th iteration, the left singular vectors U r and right singular vectors V r obtained by singular value decomposition of R 1 E k can optimize Equation (22) to
min R 1 E , R 2 E   R 1 E * t r ( U r T R 1 E V r ) + λ R 2 E   1 ,   s . t .   R E = R 1 E + R 2 E .
The two steps are executed alternately until convergence. The specific steps are
Input: original matrix R E , termination error ε , and R 0 = 0 .
Output:  R ^ 1 E ,   R ^ 2 E .
Step 1: singular value decomposition of R 1 E k :
[ U k , k , V k ] = SVD ( R 1 E k )   U r = [ u 1 , u 2 , u r ] m × r ,   V r = [ v 1 , v 2 , v r ] m × r .
Step 2: solve the optimization problem:
( R 1 E ( k + 1 ) , R 2 E ( k + 1 ) ) = arg min R 1 E , R 2 E   R 1 E * t r ( U r T R 1 E V r ) + λ R 2 E   1 ,   s . t .   R E = R 1 E + R 2 E .
Repeat steps 1 and 2 for iterative operations.
When R E R 1 E ( k + 1 ) R 2 E ( k + 1 ) F / R E F ε , stop iteration and return to R ^ 1 E ,   R ^ 2 E .
Here, the alternating direction method is used to update R 1 E k and R 2 E k to solve Equation (24). The corresponding augmented Lagrangian function can be written as
L ( R 1 E , R 2 E , Y , μ ) = R 1 E * t r ( U r T R 1 E V r ) + λ R 2 E   1 = + t r ( Y T ( R E R 1 E R 2 E ) ) + μ 2 R E R 1 E R 2 E   F 2
where the Lagrange multiplier is Y , and the penalty coefficient is μ . The solution steps for the k th iteration are as follows:
(1)
Given R 1 E k , Y k , and μ k , we have
R 2 E ( k + 1 ) = arg min R 2 E   L ( R 1 E k , R 2 E , Y k , μ k ) = arg min R 2 E   { λ R 2 E   1 + t r ( Y k T ( R E R 1 E R 2 E ) ) + μ k 2 R E R 1 E k R 2 E   F 2 } = arg min R 2 E   { λ R 2 E   1 + μ k 2 R 2 E ( R E R 1 E k + μ k 1 Y k )   F 2 } .
According to the soft-thresholding method [40], we can know that
R 2 E ( k + 1 ) = S λ μ k 1 ( R E R 1 E k + μ k 1 Y k ) , S δ ( Z ) = max ( z i , j δ , 0 ) z i , j z i , j 2
where z i , j is the element of the i th row and j th column of the matrix Z = R E R 1 E k + μ k 1 Y , δ = λ μ k 1 .
(2)
Given R 2 E k , Y k , and μ k , we have
R 1 E ( k + 1 ) = arg min R 1 E   L ( R 1 E , R 2 E ( k + 1 ) , Y k , μ k ) = arg min R 2 E   { R 1 E   * t r ( ( U r T ) k R 1 E ( V r ) k ) + t r ( Y k T ( R E R 1 E R 2 E ( k + 1 ) ) ) + μ k 2 R E R 1 E R 2 E ( k + 1 )   F 2 }
which is equivalent to
R 1 E ( k + 1 ) = arg min R 2 E   { R 1 E   * + μ k 2 R 1 E [ R E R 2 E ( k + 1 ) + μ k 1 ( Y k + ( U r ) k ( V r T ) k ) ]   F 2 } .
According to the Singular Value Threshold (SVT) method [37], the above problem turns into
R 1 E ( k + 1 ) = μ k 1 ( R E R 2 E ( k + 1 ) + μ k 1 ( Y k + ( U r ) k ( V r T ) k ) ) = U S μ k 1 ( Σ ) V T
S μ k 1 ( Σ ) = d i a g { max [ ( σ i μ k 1 ) , 0 ] }
where vectors U and V are the left and right singular vectors of matrix Q = R E R 2 E ( k + 1 ) + μ k 1 ( Y k + ( U r ) k ( V r T ) k ) , and σ i being the i th singular value.
(3)
Given R 1 E ( k + 1 ) , R 2 E ( k + 1 ) , and μ k , we have
Y ( k + 1 ) = Y k + μ k ( R E R 1 E ( k + 1 ) R 2 E ( k + 1 ) ) .
(4)
Given R 1 E ( k + 1 ) , R 2 E ( k + 1 ) , and Y ( k + 1 ) , we have
μ k + 1 = min ( μ k ρ , μ max )
where ρ > 1 is a constant. The larger μ is, the faster the algorithm converges.
The above iterative algorithm updates { R 1 E ( k + 1 ) , R 2 E ( k + 1 ) } , which in turn updates { U r , V r } in step 1, until the iteration stops. The rank r of the reconstructed matrix is needed in the above process, but in practice, it is usually impossible to predict. If r is set too small, there would be reduced low-rank parts and inaccurate reconstruction. And if r is set too large, there would be slow convergence and poor timeliness. To reduce the rank dependence, the rank estimation is performed by the soft-thresholding method. Meanwhile, the SVD operation is replaced by the fast iteration to reduce the algorithm’s computational workload and complexity.
By introducing a soft thresholding operator into the calculation of singular values in Equation (30), Equation (31) is abbreviated as
S μ k 1 ( Σ i , j ) = { σ i μ k 1 , i f   σ i > μ k 1 0 , otherwise
Now, R 1 E ( k + 1 ) = U S μ k 1 ( Σ i , j ) V T . The above equation tells that the value of non-zero elements in S μ k 1 ( Σ i , j ) come by threshold   μ k 1 , where the rank of the estimated matrix is the number of the non-zero elements, denoted as r ^ . So, the algorithm steps become
Input: original matrix R E and termination error ε . R 0 = 0 , Y 0 = D / J ( D ) , J ( D ) = max ( D 2 , λ 1 D ) , ( U r ) 0 = 0 , ( V r ) 0 = 0 , r 0 = 0 .
Output:  R ^ 1 E , R ^ 2 E .
Step 1: R 2 E ( k + 1 ) = S λ μ k 1 ( R E R 1 E k + μ k 1 Y k ) .
Step 2:  R 1 E ( k + 1 ) = μ k 1 ( R E R 2 E ( k + 1 ) + μ k 1 ( Y k + ( U r ) k ( V r T ) k ) ) .
Step 3: sum ( U r ) k + 1 and ( V r ) k + 1 , and estimate the rank r k + 1 .
Step 4:  Y ( k + 1 ) = Y k + μ k ( R E R 1 E ( k + 1 ) R 2 E ( k + 1 ) ) .
Step 5:  μ k + 1 = min ( μ k ρ , μ max ) .
Repeat steps 1–5, until R E R 1 E k R 2 E k F / R E F ε .
The penalty parameter μ is automatically updated to increase the number of iterations and speed up convergence. A constant ρ greater than 1 must be chosen appropriately since a value too large would cause the algorithm to diverge. To demonstrate the speed of the proposed algorithm compared to ordinary iterative algorithms, Table 2 provides the specific complexity of each step of the algorithm.

4. Simulation Experiments

This section verifies the performance of the TNNM-FIC-STAP by comparing it with traditional STAP (T-STAP) and traditional coprime STAP (C-STAP) in terms of root mean square error (RMSE), system DOF, Beampatterns, and SCNR. The main parameters of the airborne radar system are listed in Table 3, and β = 1 . We set N 1 = M 1 = 2 , N 2 = M 2 = 3   ρ = 1.5 , μ max = 10 10 , μ 0 = 1.25 / R E 2 . All simulation results are averaged over 500 Monte Carlo experiments.

4.1. The Root Mean Square Error (RMSE)

The root mean square error (RMSE) of matrix estimation is defined as
RMSE = E { R ^ c R c 2 } R c 2
Figure 3 shows the relationships between the RMSE and the number of samples for the sample covariance matrix (SCM) method, the minimum nuclear norm (NNM) method, and the proposed method (TNNM-FIC).
It shows that the performance of NNM and TNNM-FIC is superior to that of SCM. The reconstruction errors of TNNM-FIC and NNM are also relatively stable. The overall relative error of TNNM-FIC reconstruction decreases rapidly, and the convergence speed is much faster than NNM and SCM. The convergence speed of the NNM is also fast, and as the number of iterations increases, the overall relative error of the NNM remains basically unchanged, while the TNNM-FIC is still decreasing. Although the overall relative error of the NNM can be achieved, it requires a larger sample size.

4.2. DOF

Compare the performance of the three algorithms from the system DOF. The ratio of DOF is defined as
γ ( K ) = L ( K ) / L T N N M F I C ( K )
where L T N N M F I C is the maximum DOF of the TNNM-FIC-STAP when N = M = K .
It can be seen from Equation (36) that the size of DOF is proportional to γ ( K ) . As shown in Figure 4, the C-STAP acquires a larger DOF in the virtual domain than the T-STAP due to the different structure. The TNNM-FIC-STAP uses the low-rank characteristics of the matrix to fill its holes, fully utilizes the virtual array sensors and virtual pulses in the difference structure, and its system DOF is relatively maximized. When K is any value, the DOF of the TNNM-FIC-STAP is higher than that of the other two algorithms. Especially when the value of K is small, the superiority of the TNNM-FIC-STAP is most significant. As its value grows, the DOF of the C-STAP becomes closer and closer to it until it is the same. However, the TNNM-FIC-STAP is more complex than the C-STAP algorithm, so its value must be selected reasonably in practical applications.

4.3. Beampatterns

Figure 5 shows the space-time Beampatterns of T-STAP, C-STAP and TNNM-FIC-STAP. The optimal Beampatterns are formed at the target location (normalized spatial frequency of 0.1 and normalized temporal frequency of −0.2) for the three algorithms. However, the TNNM-FIC-STAP shows a narrow main beam and lower resolution than the other two algorithms in the Doppler and spatial Beampatterns. Of all the three algorithms, there are notches as deep as 30 dB on the main clutter ridges, which can be suppressed effectively. As the system DOF increases, the notches become deeper, and the main clutter is easier to suppress.
To better compare the notch-forming capabilities of various algorithms, Figure 6 presents the spatial and temporal channel direction diagrams where the target is located. As can be seen from the figure, the proposed TNNM-FIC-STAP algorithm achieves the optimal spatial and temporal resolution at the target’s position. The closer to the grid boundary, the relatively poorer the resolution becomes.

4.4. SCNR

Figure 7 shows that the SCNR performance of the TNNM-FIC-STAP is superior to the T-STAP and C-STAP, indicating that the estimation of CCM is the most accurate. Due to the relatively small number of sensors and pulses in the C-STAP, there is a significant estimation error in the clutter subspace, resulting in a decrease in system performance. Both C-STAP and TNNM-FIC-STAP adopt the coprime structure and difference operation to extend their virtual uniform linear array and virtual uniform linear impulse train, thereby increasing the virtual space DOF and virtual time DOF, resulting in a higher virtual DOF for the system filter. The TNNM-FIC-STAP performs better than C-STAP, because the truncated kernel norm imposes a more stringent constraint on the clutter rank, and its iterative estimation of the CCM leads to smaller errors.

5. Conclusions

A fast iterative coprime STAP method based on truncated kernel norm minimization has been proposed. This method uses limited physical array elements and physical pulses, based on difference operations, to generate difference coarrays and copulse, and provides a covariance matrix with missing data. Then, the truncated kernel norm is used to constrain the clutter rank to solve the CCM estimation problem, ensuring the low-rank characteristics of the CCM. The low-rank matrix recovery model is built on the low-rank constraint by truncated kernel norm minimization, which transforms the non-convex optimization problem into a convex one, along with soft-threshold rank estimation of each matrix produced at each iteration. This results in a fast iterative algorithm that does not depend on the matrix rank, which reduces the computational complexity of the algorithm and improves the computational efficiency. The data of the virtual CCM in the difference structure are completed, which improves the system space DOF and time DOF and improves the target temporal resolution and spatial resolution of STAP detection under limited physical array sensors and physical pulses. That is, to meet the scheduled performance of the early warning radar system, the hardware equipment and power loss are relatively reduced, which can effectively increase the UAV equipment and endurance time. While meeting the predetermined performance of the early warning radar system, the hardware equipment and power loss required by the system are relatively reduced, which can effectively increase the unmanned aerial vehicles equipment and battery life. Compared to traditional STAP, the proposed method retains the coprime sampling structure to reduce the mutual coupling effect. Theoretical analysis shows the improved DOF of the system filter under a limited number of N physical sensors and M physical pulses. In simulation experiments, this method has been proven to be able to achieve higher target detection and reduce mutual coupling effects under limited physical array sensors and pulses, which is better than typical STAP methods such as T-STAP and C-STAP. Practical experiments show that the TNNM-FIC-STAP provides more accurate clutter suppression and less environmental influence when detecting dynamic targets. In fact, radar may lose some data. These missing data will affect the robustness of the algorithm, so it is necessary to study their robustness and the percentage of allowed data loss in the next step.

Author Contributions

Methodology, M.L.; Software, H.L.; Formal analysis, M.L. and H.L.; Writing—original draft, M.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Melvin, W.L. A STAP overview. IEEE Aerosp. Electron. Syst. Mag. 2004, 19, 19–35. [Google Scholar] [CrossRef]
  2. Cao, J.; Wang, T.; Wang, D. Beam-Space Post-Doppler Reduced-Dimension STAP Based on Sparse Bayesian Learning. Remote Sens. 2024, 16, 307. [Google Scholar] [CrossRef]
  3. Sun, Y.; Xie, J.; Gong, Y.; Zhang, Z.; Wang, L. Adaptive interference mitigation space-time array reconfiguration by joint selection of antenna and delay tap. IET Radar Sonar Navig. 2024, 18, 448–462. [Google Scholar] [CrossRef]
  4. Chen, Z.; Zhang, L.; Zhou, Y.; Lin, C.; Tang, S.; Wan, J. A non-adaptive space-time clutter canceller for multi-channel synthetic aperture radar. Signal Process. IET 2019, 13, 472–479. [Google Scholar] [CrossRef]
  5. Wang, W.; Zou, L.; Wang, X. A Novel Two-Level Nested STAP Strategy for Clutter Suppression in Airborne Radar. Math. Probl. Eng. 2019, 2019, 2540858. [Google Scholar] [CrossRef]
  6. Zhao, X.; He, Z.; Wang, Y.; Sun, G. Reduced-dimension STAP using a modified generalised sidelobe canceller for collocated MIMO radars. IET Radar Sonar Navig. 2018, 12, 1476–1483. [Google Scholar] [CrossRef]
  7. Wang, W.; Zou, L.; Wang, X.; Yang, Y. Deterministic-aided single dataset STAP method based on sparse recovery in heterogeneous clutter environments. EURASIP J. Adv. Signal Process. 2018, 2018, 24. [Google Scholar] [CrossRef]
  8. 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]
  9. Wang, W.-Q. Frequency Diverse Array SAR With Fast-Time Space-Time Adaptive Processing for Subsurface Imaging. IEEE Access 2023, 11, 122128–122140. [Google Scholar] [CrossRef]
  10. Ren, B.; Wang, T. Space-Time Adaptive Processing Based on Modified Sparse Learning via Iterative Minimization for Conformal Array Radar. Sensors 2022, 22, 6917. [Google Scholar] [CrossRef]
  11. Zhang, S.; Wang, T.; Liu, C.; Wang, D. A Space-Time Adaptive Processing Method Based on Sparse Bayesian Learning for Maneuvering Airborne Radar. Sensors 2022, 22, 5479. [Google Scholar] [CrossRef]
  12. Fu, D.; Liao, G.; Xu, J. Clutter Subspace Characteristics-Aided Space-Time Adaptive Outlier Sample Selection Method. Sensors 2021, 21, 3108. [Google Scholar] [CrossRef]
  13. Chen, N.; Wei, P.; Ai, X. A Robust STAP Method Based on Alternating Projection and Clutter Cancellation. IEEE Access 2019, 7, 179880–179890. [Google Scholar] [CrossRef]
  14. 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]
  15. Zhang, T.; Wang, Z.; Xing, M.; Zhang, S.; Wang, Y. Research on Multi-Domain Dimensionality Reduction Joint Adaptive Processing Method for Range Ambiguous Clutter of FDA-Phase-MIMO Space-Based Early Warning Radar. Remote Sens. 2022, 14, 5536. [Google Scholar] [CrossRef]
  16. Fa, R.; de Lamare, R.C.; Wang, L. Reduced-Rank STAP Schemes for Airborne Radar Based on Switched Joint Interpolation, Decimation and Filtering Algorithm. IEEE Trans. Signal Process. 2010, 58, 4182–4194. [Google Scholar] [CrossRef]
  17. Wong, M.; Aboutanios, E.; Rosenberg, L. A STAP Detection Scheme for Low Sample Support Maritime Environments. IEEE Trans. Aerosp. Electron. Syst. 2023, 59, 5671–5683. [Google Scholar] [CrossRef]
  18. Zhang, W.; Han, M.; He, Z.; Li, H. Data-dependent reduced-dimension STAP. IET Radar Sonar Navig. 2019, 13, 1287–1294. [Google Scholar] [CrossRef]
  19. Shen, M.; Wang, J.; Wu, D.; Zhu, D. An Efficient Moving Target Detection Algorithm Based on Sparsity-Aware Spectrum Estimation. Sensors 2014, 14, 17055–17067. [Google Scholar] [CrossRef]
  20. Tong, Y.; Wang, T.; Wu, J. Improving EFA-STAP performance using persymmetric covariance matrix estimation. IEEE Trans. Aerosp. Electron. Syst. 2015, 51, 924–936. [Google Scholar] [CrossRef]
  21. Ciuonzo, D.; Orlando, D.; Pallotta, L. On the maximal invariant statistic for adaptive radar detection in partially homogeneous disturbance with persymmetric covariance. IEEE Signal Process. Lett. 2016, 23, 1830–1834. [Google Scholar] [CrossRef]
  22. Sun, K.; Meng, H.; Wang, Y.; Wang, X. Direct data domain STAP using sparse representation of clutter spectrum. arXiv 2010, arXiv:1008.4184. [Google Scholar] [CrossRef]
  23. Rivera Velázquez, J.M.; Latorre, L.; Mailly, F.; Nouet, P. A New Algorithm for Fault Tolerance in Redundant Sensor Systems Based on Real-Time Variance Estimation. IEEE Sens. J. 2022, 22, 15410–15418. [Google Scholar] [CrossRef]
  24. Ma, H.; Mao, X.; Wang, X.; Gao, Y. A new array extension configuration method based on nesting and minimum redundancy. IET Radar Sonar Navig. 2023, 17, 748–758. [Google Scholar] [CrossRef]
  25. Wei, S.; Zhu, G.; Su, Y. A Novel Sparse Array Configuration for Direction of Arrival Estimation with Increased Uniform Degrees of Freedom and Reduced Mutual Coupling. Sensors 2024, 24, 808. [Google Scholar] [CrossRef] [PubMed]
  26. Wang, X.; Zhao, L.; Jiang, Y. Super Augmented Nested Arrays: A New Sparse Array for Improved DOA Estimation Accuracy. IEEE Signal Process. Lett. 2024, 31, 26–30. [Google Scholar] [CrossRef]
  27. Li, N.; Zhang, X.-K.; Zong, B.; Lv, F.; Xu, J.; Wang, Z. An off-grid direction-of-arrival estimator based on sparse Bayesian learning with three-stage hierarchical Laplace priors. Signal Process. 2024, 218, 109371. [Google Scholar] [CrossRef]
  28. Liu, M.; Wang, X.; Zou, L. Robust STAP With Reduced Mutual Coupling and Enhanced DOF Based on Super Nested Sampling Structure. IEEE Access 2019, 7, 175420–175428. [Google Scholar] [CrossRef]
  29. Wang, X.; Yang, Z.; Huang, J. Robust space-time adaptive processing for airborne radar with coprime arrays in presence of gain and phase errors. IET Radar Sonar Navig. 2020, 15, 75–84. [Google Scholar] [CrossRef]
  30. Li, B.; Lu, L.; Zhou, C. SADOE: Sequential-based angle-Doppler off-grid estimation with coprime sampling structures for space-time adaptive processing. IET Radar Sonar Navig. 2021, 15, 775–787. [Google Scholar] [CrossRef]
  31. Liu, M.; Tang, B.; Zheng, X.; Wang, Q.; Wang, S.; Wang, X.; Fang, M.; Zhu, J.; Feng, C. STAP for Airborne Radar Based on Dual-Frequency Space-Time Coprime Sampling Structure. IEEE Access 2020, 8, 138673–138681. [Google Scholar] [CrossRef]
  32. Liu, S.; Mao, Z.; Zhang, Y.D.; Huang, Y. Rank Minimization-Based Toeplitz Reconstruction for DoA Estimation Using Coprime Array. IEEE Commun. Lett. 2021, 25, 2265–2269. [Google Scholar] [CrossRef]
  33. Qiao, H.; Pal, P. Unified analysis of co-array interpolation for direction-of-arrival estimation. In Proceedings of the 2017 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), New Orleans, LA, USA, 5–9 March 2017; pp. 3056–3060. [Google Scholar]
  34. Zhou, C.; Gu, Y.; Shi, Z.; Zhang, Y.D. Off-Grid Direction-of-Arrival Estimation Using Coprime Array Interpolation. IEEE Signal Process. Lett. 2018, 25, 1710–1714. [Google Scholar] [CrossRef]
  35. Zhou, C.; Gu, Y.; Fan, X.; Shi, Z.; Mao, G.; Zhang, Y.D. Direction-of-Arrival Estimation for Coprime Array via Virtual Array Interpolation. IEEE Trans. Signal Process. 2018, 66, 5956–5971. [Google Scholar] [CrossRef]
  36. Stoica, P.; Wang, Z.; Li, J. Robust Capon beam map. IEEE Signal Process. Lett. 2003, 10, 172–175. [Google Scholar] [CrossRef]
  37. Cai, J.-F.; Candès, E.J.; Shen, Z. A Singular Value Thresholding Algorithm for Matrix Completion. SIAM J. Optim. 2010, 20, 1956–1982. [Google Scholar] [CrossRef]
  38. Lin, Z.; Liu, R.; Su, Z. Linearized Alternating Direction Method with Adaptive Penalty for Low-Rank Representation. In Proceedings of the Neural Information Processing Systems 2011, Granada, Spain, 12–15 December 2011. [Google Scholar]
  39. Hu, Y.; Zhang, D.; Ye, J.; Li, X.; He, X. Fast and Accurate Matrix Completion via Truncated Nuclear Norm Regularization. IEEE Trans. Pattern Anal. Mach. Intell. 2013, 35, 2117–2130. [Google Scholar] [CrossRef]
  40. Donoho, D. De-noising by soft-thresholding. IEEE Trans. Inf. Theory 1995, 41, 613–627. [Google Scholar] [CrossRef]
Figure 1. Coprime configuration: (a) coprime array; (b) coprime PRI.
Figure 1. Coprime configuration: (a) coprime array; (b) coprime PRI.
Sensors 24 04007 g001
Figure 2. Difference coarray and copulse: (a) difference array; (b) difference impulses.
Figure 2. Difference coarray and copulse: (a) difference array; (b) difference impulses.
Sensors 24 04007 g002
Figure 3. Relationship of the RMSE and the number of samples.
Figure 3. Relationship of the RMSE and the number of samples.
Sensors 24 04007 g003
Figure 4. The ratio of DOF.
Figure 4. The ratio of DOF.
Sensors 24 04007 g004
Figure 5. Beampatterns in angle Doppler: (a) T-STAP; (b) C-STAP; (c) TNNM-FIC-STAP.
Figure 5. Beampatterns in angle Doppler: (a) T-STAP; (b) C-STAP; (c) TNNM-FIC-STAP.
Sensors 24 04007 g005
Figure 6. Beampatterns: (a) spatial domain; (b) Doppler domain.
Figure 6. Beampatterns: (a) spatial domain; (b) Doppler domain.
Sensors 24 04007 g006
Figure 7. SCNR.
Figure 7. SCNR.
Sensors 24 04007 g007
Table 1. List of abbreviations.
Table 1. List of abbreviations.
AbbreviationFull Name
STAPspace-time adaptive processing
TKNMtruncated kernel norm minimization
CCMclutter covariance matrix
SINRsignal-to-interference-to-noise ratio
DOFdegrees of freedom
TNNM-FIC-STAPfast iterative coprime STAP algorithm
PRIpulse repetition interval
CPIcoherent processing interval
ULAsuniform linear arrays
CNCMclutter plus noise covariance matrix
SNRsignal-to-noise ratio
SVDsingular value decomposition
SVTsingular value threshold
T-STAPtraditional STAP
C-STAPtraditional coprime STAP
RMSEroot mean square error
SCMsample covariance matrix
NNMminimum nuclear norm
Table 2. Computational complexity.
Table 2. Computational complexity.
StepComputational Complexity
Step 1 O ( ( 4 N 1 N 2 2 N 2 + 1 ) 2 × ( 4 M 1 M 2 2 M 2 + 1 ) 2 )
Step 2 O ( ( 4 N 1 N 2 2 N 2 + 1 ) 3 × ( 4 M 1 M 2 2 M 2 + 1 ) 3 )
Step 3 O ( ( 4 N 1 N 2 2 N 2 + 1 ) 2 × ( 4 M 1 M 2 2 M 2 + 1 ) 2 )
Step 4 O ( ( 4 N 1 N 2 2 N 2 + 1 ) 2 × ( 4 M 1 M 2 2 M 2 + 1 ) 2 )
Step 5 O ( 1 )
Table 3. Main parameters of STAP radar.
Table 3. Main parameters of STAP radar.
SymbolQuantityValue
N the number of sensors6
M the number of pulses6
λ carrier wavelength0.2 m
T r minimal PRI0.5 ms
v r radar velocity100 m/s
σ n 2 noise power1 dB
N c number of independent clutter patches361
CNRclutter to noise ratio30 dB
SNRSignal-to-noise ratio0 dB
normalized angle frequency of target0.1
normalized Doppler frequency of target0.2
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

Li, M.; Li, H. A Novel Fast Iterative STAP Method with a Coprime Sampling Structure. Sensors 2024, 24, 4007. https://doi.org/10.3390/s24124007

AMA Style

Li M, Li H. A Novel Fast Iterative STAP Method with a Coprime Sampling Structure. Sensors. 2024; 24(12):4007. https://doi.org/10.3390/s24124007

Chicago/Turabian Style

Li, Mingfu, and Hui Li. 2024. "A Novel Fast Iterative STAP Method with a Coprime Sampling Structure" Sensors 24, no. 12: 4007. https://doi.org/10.3390/s24124007

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