Space-Time Adaptive Processing Based on Modified Sparse Learning via Iterative Minimization for Conformal Array Radar

Space-time adaptive processing (STAP) is a well-known technique for slow-moving target detection in the clutter spreading environment. For an airborne conformal array radar, conventional STAP methods are unable to provide good performance in suppressing clutter because of the geometry-induced range-dependent clutter, non-uniform spatial steering vector, and polarization sensitivity. In this paper, a knowledge aided STAP method based on sparse learning via iterative minimization (SLIM) combined with Laplace distribution is proposed to improve the STAP performance for a conformal array. The proposed method can avoid selecting the user parameter. the proposed method constructs a dictionary matrix that is composed of the space-time steering vector by using the prior knowledge of the range cell under test (CUT) distributed in clutter ridge. Then, the estimated sparse parameters and noise power can be used to calculate a relatively accurate clutter plus noise covariance matrix (CNCM). This method could achieve superior performance of clutter suppression for a conformal array. Simulation results demonstrate the effectiveness of this method.


Introduction
The conformal array radar has many advantages, including minimal payload weight, the potential larger effective aperture, raised field of scan without the demand for cumbersome mechanical couplings, and easy integration with various sensor functions. The conformal array design also avoids signal modulation caused by antenna rotation [1]. An efficient model of the conformal array signal is a necessary prerequisite for conformal array signal processing. Compared with the traditional linear array or planar array, the conformal array has a different element pattern and normal because of the varying curvature. However, previous study assumed relatively simple element pattern and ignored the influence of the polarization and mutual coupling in clutter model.
Space-time adaptive processing (STAP) is a technique for detecting moving target in the presence of a spreading clutter background [2]. The optimal STAP weight vector is expressed as the product of the inverse of clutter plus noise covariance matrix (CNCM) and the target space-time steering vector [3]. When detecting whether the target exists in the range cell, the CNCM of the range cell is usually estimated by the training samples which are selected from the range cells adjacent to the cell under test (CUT) [4]. Traditional STAP techniques using a linear array have the acceptable performance, and the relationship between the spatial and Doppler frequencies is range-independent. Therefore, the training samples from adjacent range cells are independent and identically distributed (IID) and can estimate relatively accurate CNCM so that the STAP processor can be effectively produced to improve the output signal-interference-noise-ratio (SINR) in the test cell [5]. However, in a complex configuration, e.g., conformal array, the relationship between the spatial and Doppler frequencies is non-linear and range-dependent, and the clutter appears nonstationary over the range [6]. The Reed Mallett Brennan (RMB) rule [7] suggest that the (2) The proposed method uses the accurate information of array geometry and radar system parameters to construct the CNCM based on the new received data model. (3) The original SLIM method has a hyper-parameter that should be chosen by user. The KA-MSLIM-STAP method utilizes the Laplace distribution to avoid selecting the user parameter. It can provide an excellent performance by using the CUT data only. (4) The proposed method is an iterative algorithm. It finds a local minimum of the cost function, but it converges rapidly. This method has lower computational complexity compared with other SR-STAP methods.
The remainder of this paper is as follows. Section 2 introduces the signal and clutter model based on the conformal array. Section 3 discusses the key theory of the proposed method and illustrate the advantages over other methods. Numerical experiments with simulated data are carried out in Section 4 to illustrate the performance of the proposed method. A summary is provided in Section 5.
Notations: We use boldface lower case for vectors (a), upper case for matrices (A) and the italics for scalars (a). The transpose, conjugate transpose and inverse of a matrix are represented by [·] T [·] H and [·] −1 respectively. | · | stands for the absolute value. · is Euclidean norm for vector and Frobenius norm for matrix. · 0 represents the l 0 norm. · 1 is the l 1 norm. ⊗ represents the Kronecker product. denotes the Hadamard product. Diag(·) is Diagonal matrix with the vector on its diagonal. I N is an N × N identity matrix.

Signal Steering Vector Model
Consider a conformal array which consists of N radiating elements. As shown in Figure 1, ϕ and θ represent the azimuth and elevation angle of a clutter scattering point. v stands for the radar velocity.
(3) The original SLIM method has a hyper-parameter that should be chosen by us KA-MSLIM-STAP method utilizes the Laplace distribution to avoid selecting t parameter. It can provide an excellent performance by using the CUT data on (4) The proposed method is an iterative algorithm. It finds a local minimum of t function, but it converges rapidly. This method has lower computational com compared with other SR-STAP methods.
The remainder of this paper is as follows. Section 2 introduces the signal and model based on the conformal array. Section 3 discusses the key theory of the pr method and illustrate the advantages over other methods. Numerical experimen simulated data are carried out in Section 4 to illustrate the performance of the pr method. A summary is provided in Section 5. Notations

 
Diag  is Diagonal matrix with the vector on its diagonal. N I is an identity matrix.

Signal Steering Vector Model
Consider a conformal array which consists of N radiating elements. As sh Figure 1,  and  represent the azimuth and elevation angle of a clutter sca point. v stands for the radar velocity. The v is written as where  is crab angle.
The unit vector pointing to propagating direction of   ,   can be represent The v is written as where ψ is crab angle. The unit vector pointing to propagating direction of (ϕ, θ) can be represented as A circular arc array is used here, and the configuration of conformal array is shown in Figure 2. The blue dots denote the position of the elements and red arrows represent the normal direction of elements. λ stands for the radar wavelength. The array is composed of N elements and array inter-element spacing is λ/2. The position pointing vector of the ith antenna element is The conventional normalized array manifold vector of the clutter patch from (ϕ, θ) can be represented as a N × 1 vector where ϕ is the azimuth angle and θ denotes elevation angle.
We design a tangential individually-polarized dipole conformal array which is composed of dipoles that are placed tangentially to circumference. Due to the different orientations, their diverse patterns sample polarization diversely, making the array polarization-sensitive.
We denote that f i is the ith element vector pattern. It can be decomposed into two components of orthogonal polarization, and it is shown as where f i ϕ (ϕ, θ) and f i θ (ϕ, θ) are the ϕ polarized component and θ polarized component, relatively.
Because the installation direction and polarization of each element in the conformal array are different, the element vector patterns are also different in the global coordinate system.
In general, each element vector pattern is determined in their local coordinate, where z is the normal to the surface, x is the horizontal tangent of the carrier surface and y is the vertical tangent of the carrier surface. Thus, we can obtain the element vector patterñ f i ( x, y, z) in the local coordinate system. In this paper,f i ( x, y, z) = [1, 0, 0]. This suggests that the array elements are placed along the horizontal tangent of the carrier surface. We should transform the element pattern from local coordinate to global coordinate to acquire the expression of the element pattern in the global coordinate, and the exact expression of the conformal array spatial steering vector can be further obtained. The Euler rotation matrix [30,31] is a helpful tool for the transformation. As for the coordinate transformation (x, y, z) ⇒ ( x, y, z) , the corresponding rotation matrix can be written as where The above represents the rotation matrix with the x-axis as the rotating axis. α x is the rotation angle. Counterclockwise rotation is positive. Similarly, With the Euler rotation matrix, we can get So, we can transform the element pattern from the local coordinate to the global coordinate by With regard to the polarized components transformation from the Cartesian coordinate system to the spherical coordinate system, the following mathematical manipulations are required.
Then, we can obtain the f i ϕ (ϕ, θ) and f i θ (ϕ, θ). In general, for the common conformal array (e.g., cylindrical, circular, conical, and spherical array), two successive Euler rotation matrices are usually adequate. For the generalization of the conformal array signal model, three Euler rotation matrices mentioned would be used here. Now, conformal array spatial steering vector can be expressed in the following form: where Given the influence of mutual coupling, a mutual coupling matrix Z will be leftmultiplied to form a new spatial steering vector matrix. For simplicity, we ignore the influence of mutual coupling in following section.
We assumed that array transmits a series of M pulse at a invariant pulse repetition frequency (PRF) f r in a coherent processing interval (CPI). The normalized Doppler steering vector can be written as a M × 1 vector. The conformal array signal space-time steering vector can be expressed as

Conformal Array Clutter Model and Received Data Model
In clutter simulation, we make this assumption that the radar coherent processing interval (CPI) is short enough so that we can think that the velocity, RCS and polarization characteristics of the target and the RCS, polarization characteristics of the clutter patch are constant in a CPI. In this paper, we assumed that the wave received by radar array is the completely polarized wave whose polarized state does not change with time. From the above discussion, we can obtain that the conformal array can receive the polarized information of the wave. The polarization angle is µ ∈ [0, π/2], and the polarization phase difference is β ∈ [−π, π], where the definition of the polarization parameters is given in Ref. [32].
The Jones vector which can be used to express the completely polarized wave is showed as With the above description, the data of the clutter patch from azimuth ϕ and elevation θ angle c ∈ C N M×1 can be expressed as where ς stands for a complex scalar with ς ∼ CN(0, 1), α represents the complex amplitude of the clutter patch, and e p denotes the polarized state of the echo wave from clutter patch. Therefore, the received data of the kth range cell can be represented as ς l,c α l,c s l,c + n (18) where N l is the number of ambiguous range gates, N c indicates the number of statistically independent clutter patches of each range cell, and s l,c = S m (ϕ c , θ l ) · e p,(l,c) , where e p,(l,c) is the polarization vector of the echo of the l, cth clutter patch, and n ∈ C N M×1 is the Gaussian white thermal noise vector. The detection can be represented by following [25].
where hypotheses H 1 and H 0 correspond to target presence and absence, respectively, α t is the complex amplitude of the target and s t denotes the target space-time steering vector. The precise formula of the target vector can be written as where ϕ 0 and θ 0 is the azimuth angle and elevation angle of the target, respectively, e pt denotes the polarized state of the target.
Considering that clutter patches are mutually independent, clutter and noise are independent of each other. The CNCM R k ∈ C N M×N M can be expressed as where σ l,c is the complex amplitude of the l, cth clutter patch and σ 2 n is the noise power. The STAP processor can alleviate the clutter component and provide coherent gain on target by performing an inner product between the weight vector w and the data from the CUT. The well-known optimal weight vector, which maximizes the output SINR, can be written as Because the statistics of interference environment are unknown, the CNCM is usually unknown a priori. It should be estimated by some algorithms in order to obtain an excellent performance of STAP.

Dictionary Matrix of Conformal Array
The range cell can be separated into N k cells. In this way, the snapshot can be rewritten as where D is a dictionary matrix, which contains steering vectors. D could be represented by can be expressed as We note that the space-time steering vectors which is obtained from the clutter ridge are related with the array geometry and the radar system parameters. Because the elevation angle θ is constant in a certain range gate, the space-time steering vectors of clutter are determined by azimuth angle ϕ in essence. Therefore, we could get enough values of azimuth angle to produce the ideal space-time steering vectors corresponding to the clutter ridge of the CUT by the prior knowledge.
Here, we choose N k > N M, and then D is an over-complete dictionary matrix which contains enough space-time steering vectors to recover the clutter accurately. y denotes sparse coefficient vector where non-zero elements represent the presence of clutter on the space-time profile. It also contains the complex amplitude and polarization parameter of the clutter. n ∈ C N M×1 is zero-mean Gaussian noise.
When given the range ambiguity situation, the dictionary matrix D ∈ C N M×2N l N k can be rewritten as The received data from lth range gate can be represented in the following form.
In the following subsection, we remark that K = 2N l N k .

SR-STAP Model and Principle
In sparse recovery methods, we use as few atoms from dictionary matrix D as possible to represent the measurement x l contaminated by noise. Hence, the objective function is written as [24]ŷ where ε 1 is a parameter that controls the error. As all known, we cannot solve (28) because it is a NP-hard problem. Hence, instead of minimizing the l 0 norm, l 1 norm has been used in the sparse represent problem [23].
where ε 2 is a user parameter. It can be equivalent to the following optimization problem [20].
where η > 0 is the user parameter. This SR method is also known as LASSO, whose disadvantages are that the user parameter is difficult to choose a compatible value and it is easily affected by the number of steering vectors. The computational complexity increases rapidly with the improvement of the size of the dictionary and the number of measurements.

KA-MSLIM-STAP
In this framework, to utilize the above sparsity structure, we adopt the hierarchical Bayesian model [33]. We suppose that the noise vector n l is a complex Gaussian random vector. It has zero mean and covariance matrix ηI. η denotes an unknown parameter. From the assumption, we obtain the following probability density functions (PDF) which can be written as In Ref. [27], where q is the hyper-parameter chosen by user. In order to making user parameter free, Ref. [27] proposed incorporating the Bayesian information criterion (BIC) to estimate q automatically. This will introduce additional computational burden. To overcome this problem, we should consider other methods. It is reasonable that we assumed that the signal coefficients follow complex Gaussian distribution. In Ref. [28], the proposed SLIM-0 method assumed that the signal coefficients follow a complex Gaussian distribution. However, we can see from Ref. [28] that the final derived result in SLIM-0 is equivalent to the SLIM with q = 0. It does not always get the optimal performance in STAP. So, we consider another distribution like Laplace distribution.
In Ref. [29], we know that Laplace prior cannot be directly applied to coefficients. To apply Laplace distribution, a hierarchical framework is modeled for the coefficient y l .
The following prior is employed on y l : where γ = (γ 1 , γ 2 , · · · , γ K ). Then, model the γ n by using the following hyper priors: Combining the hierarchical Bayesian model, we can obtain: p(y l |υ ) = p(y l |γ )p(γ|υ )dγ = K ∏ n=1 p(y l,n |γ n )p(γ n |υ )dγ n = υ K/2 It can be seen that the p(y l |υ ) is Laplace distribution. In this paper, we assume the signal coefficients present a Laplace distribution, that is Then, we model υ as the following form The estimates of y l , υ, and η are acquired by solving a maximum a posteriori (MAP) problem as follow: min where From (31), (36), (37), and later discarding irrelevant constants, g MSLI M (y l , υ, η) can be indicated as The minimization problem of (38) could be resolved by utilizing the cyclic minimization (CM) method [34]. First, suppose we have obtained y l in the tth iteration, and the minimization problem can be written as We set the (∂/∂υ)g MSLI M to zero, and can get Then, suppose we have obtained y l , υ and η in the tth iteration, we aim to get the y l and η in the t + 1th iteration and the minimization problem can simplified as follows: We can set the complex derivative to zero and solve it for y l , this results in: where P = diag 2 y n,l /υ , n = 1, 2, · · · , K. Then, solving (44), we can get the solution Similarly, setting the (∂/∂η)g MSLI M to zero leads to The solution of (38) can be obtained via repeating above iteration (42), (45), and (46) until this convergence criterion is satisfied where ∆ is a small, positive number. If the maximum number of iterations is attained in advance, the process suspends. We initialize the iterative process by applying a match filter, so that y n,l (0) = d H n x l / d n 2 , for n = 1, 2, · · · , K, where d n is the nth column of D, and η (0) = (1/NM) x l − Dy l (0) 2 . By the above iterative process, we can obtain the relatively accurate estimate of y l and η. The CNCM can be calculated as follows: The weight vector could be acquired by (22). The KA-MSLIM-STAP method is summarized in Algorithm 1. Step 1 Initialize values of parameters as: Step 2 Repeat the following for t = 0, 1, 2, · · · υ (t+1) = K−1

Until convergence
Step 3 Calculate CNCM Compute KA-MSLIM-STAP weight w Due to the cyclic minimization character, the cost function is monotonically nonincreasing in each iteration. This will indicate that the proposed method is convergent. We note that the cost function is not convex. Global optimum is not guaranteed. However, as we test numerically, excellent estimation results can be obtained, and the proposed method converges quickly. In the following simulations, this method shows excellent performance after five iterations.
The computational complexity of the methods is measured in terms of the number of multiplications. Based on the analysis above, the proposed method consists of the estimate of y l , υ, P and η in each iteration. We can obtain that the computational complexity of the proposed method is o (NM) 3 +2K(NM) 2 + 4KNM + 2K + 2 . The computational complexity of the KASPICE-STAP method is o 2(NM) 3 +2(NM) 2 +2K(NM) 2 + 4KNM + NM + 3K . We can see the proposed method has lower computational complexity.
The proposed KA-MSLIM-STAP method could acquire better performance by using the snapshot of the testing range cell only. The proposed method also possesses advantages, including low computational complexity and no need to select user parameters compared with other SR-STAP methods.

Simulation Experiments
In this section, some simulation results will show the performance of different methods with conformal array. In Figure 2, the conformal array contains 12 elements in a circular arc. We use a uniform grid here and the size is determined as N k = 300. We will verify the performance of the sample matrix inversion (SMI) method, the angle-Doppler compensation (ADC) method, the registration-based compensation with sparse representation (SR-RBC) method, the STAP method based on knowledge aided sparse iterative covariance-based estimation (KASPICE-STAP) method, the least absolute shrinkage and selection operator (LASSO) method and the proposed KA-MSLIM-STAP method with the simulated data. In this section, we consider forward-looking conformal array in the ideal condition and nonideal condition. The range ambiguity condition is taken into account here. Then, we use the improvement factor (IF) as the measurement of performance, which can be expressed as where R denotes the exact CNCM of the CUT. The relevant radar parameters are given in Table 1.

Ideal Conditoin
The clutter data of the 300th range cell will be the input data for the spectrums. In Figure 3, there are four clutter capon spectrums figures in the ideal condition. Figure 3a-d illustrate the clutter capon spectrums estimated by the ideal CNCM, the SMI method, the KASPICE-STAP method and the proposed method. We can get that the clutter capon spectra gained by the KASPICE-STAP and the proposed method are closer to the accurate spectrum. where R denotes the exact CNCM of the CUT. The relevant radar parameters are given in Table 1.

Ideal Conditoin
The clutter data of the 300th range cell will be the input data for the spectrums. In Figure 3, there are four clutter capon spectrums figures in the ideal condition. Figure 3a-d illustrate the clutter capon spectrums estimated by the ideal CNCM, the SMI method, the KASPICE-STAP method and the proposed method. We can get that the clutter capon spectra gained by the KASPICE-STAP and the proposed method are closer to the accurate spectrum.   Compared with the ADC, SMI, and SR-RBC methods, the excellent performance could be seen to emerged at the 300th and 500th range cells with the proposed method, when the number of IID training samples is inadequate. The proposed method only utilizes the testing range cell and knowledge of the clutter ridge to construct the dictionary matrix and estimate the CNCM, which can avoid the problem of inadequate IID training samples.
This subsection has exhibited the performance analysis of various methods in ideal conditions. In following subsections, we are going to probe the factors that impair the performance of these methods.     Compared with the ADC, SMI, and SR-RBC methods, the excellent performance could be seen to emerged at the 300th and 500th range cells with the proposed method, when the number of IID training samples is inadequate. The proposed method only utilizes the testing range cell and knowledge of the clutter ridge to construct the dictionary matrix and estimate the CNCM, which can avoid the problem of inadequate IID training samples.
This subsection has exhibited the performance analysis of various methods in ideal conditions. In following subsections, we are going to probe the factors that impair the performance of these methods.  Compared with the ADC, SMI, and SR-RBC methods, the excellent performance could be seen to emerged at the 300th and 500th range cells with the proposed method, when the number of IID training samples is inadequate. The proposed method only utilizes the testing range cell and knowledge of the clutter ridge to construct the dictionary matrix and estimate the CNCM, which can avoid the problem of inadequate IID training samples.
This subsection has exhibited the performance analysis of various methods in ideal conditions. In following subsections, we are going to probe the factors that impair the performance of these methods.

Model Errors
In Figure 5, we utilize four methods in the non-ideal condition with 0.5% amplitude error and 0.5 • phase error. Figure 5a-d exhibit the clutter capon spectra estimated by the ideal CNCM, the SMI method, the KASPICE-STAP method, and the proposed method. We can also get that the clutter capon spectra gained by the KASPICE-STAP and the proposed method are closer to the accurate spectrum.

Model Errors
In Figure 5, we utilize four methods in the non-ideal condition with 0.5% amplitude error and 0.5° phase error. Figure 5a-d exhibit the clutter capon spectra estimated by the ideal CNCM, the SMI method, the KASPICE-STAP method, and the proposed method. We can also get that the clutter capon spectra gained by the KASPICE-STAP and the proposed method are closer to the accurate spectrum. Because of the presence of model errors, there is a mismatch between the steering vector in dictionary matrix and the real steering vector of the clutter. This will result in the inaccurate estimation of CNCM. The model errors consist of array amplitude and phase errors. Figure 6 shows the IF curves with various methods with model errors, where the model errors are 0.5% amplitude error and 0.5° phase error in Figure 6a and the model errors are 1.2% amplitude error and 1.2° phase error in Figure 6b. Compared with ideal condition, the proposed method demonstrates performance degradation. The results show that the performance decrease more when model errors are large. However, we can know that the effectiveness of the proposed method is still better than other conventional methods and is excellent in terms of SR STAP methods. Because of the presence of model errors, there is a mismatch between the steering vector in dictionary matrix and the real steering vector of the clutter. This will result in the inaccurate estimation of CNCM. The model errors consist of array amplitude and phase errors. Figure 6 shows the IF curves with various methods with model errors, where the model errors are 0.5% amplitude error and 0.5 • phase error in Figure 6a and the model errors are 1.2% amplitude error and 1.2 • phase error in Figure 6b. Compared with ideal condition, the proposed method demonstrates performance degradation. The results show that the performance decrease more when model errors are large. However, we can know that the effectiveness of the proposed method is still better than other conventional methods and is excellent in terms of SR STAP methods.

ICM
Intrinsic clutter motion (ICM) is an issue that will exist in practice. Many factors, due to both the radar system design and the environment, may in practice cause small pulseto-pulse fluctuations in the clutter return. The presence of ICM requires a wider clutter

ICM
Intrinsic clutter motion (ICM) is an issue that will exist in practice. Many factors, due to both the radar system design and the environment, may in practice cause small pulse-to-pulse fluctuations in the clutter return. The presence of ICM requires a wider clutter notch and more adaptive degrees of freedom for effective cancellation. According to Ref. [35], considering the ICM influence, the received data can be rewritten as where γ t is the temporal covariance weight matrix. In Figure 7, we utilize four methods in the non-ideal condition with ICM and the standard variance of clutter spectrum spread is 0.22 m/s. Figure 7a-d exhibit the clutter capon spectra estimated by the ideal CNCM, the SMI method, the KASPICE-STAP method, and the proposed method. We can get that the clutter capon spectra gained by the KASPICE-STAP and the proposed method are closer to the accurate spectrum.

ICM
Intrinsic clutter motion (ICM) is an issue that will exist in practice. Many factors, to both the radar system design and the environment, may in practice cause small pu to-pulse fluctuations in the clutter return. The presence of ICM requires a wider clu notch and more adaptive degrees of freedom for effective cancellation. According to R [35], considering the ICM influence, the received data can be rewritten as where t γ is the temporal covariance weight matrix.
In Figure 7, we utilize four methods in the non-ideal condition with ICM and standard variance of clutter spectrum spread is 0.22 m/s. Figure 7a-d exhibit the clu capon spectra estimated by the ideal CNCM, the SMI method, the KASPICE-ST method, and the proposed method. We can get that the clutter capon spectra gained the KASPICE-STAP and the proposed method are closer to the accurate spectrum.  In the presence of ICM, the IF curves with different methods have been shown in Figure 8. The standard variance of the clutter spectrum spread is 0.22 m/s in Figure 8a and 0.56 m/s in Figure 8b. We can see that the proposed method still preserves desirable performance. In the presence of ICM, the IF curves with different methods have been shown in Figure 8. The standard variance of the clutter spectrum spread is 0.22 m/s in Figure 8a and 0.56 m/s in Figure 8b. We can see that the proposed method still preserves desirable performance.

Target Detection
To further prove the performance of different methods, we use the probability of detection (PD) curves here. PD is a function of input signal-to-noise ratio (SNR) with a constant probability of false alarm fa P . The detection curves will be acquired by the cell average constant false alarm rate (CA-CFAR) detector. The CUTs in the area of the 8-12th Doppler bins and the 301-500th range cells are selected as the 1000 runs in the Monte-Carlo experiments. We set that fa P is 4 10  .
In Figure 9, we show PD versus SNR with ADC, SR-RBC, KASPICE, and the proposed method. We can see that, compared with the detection performance of the KASPICE, there is a slight advantage of about 1 dB in that corresponding to the proposed method. Compared with the SR-RBC and ADC methods, the proposed method demonstrates great detection performance. We can conclude that the result obtained from the probability of detection performance is similar to that obtained from the performance of STAP. This can further prove the superiority of the proposed method.

Target Detection
To further prove the performance of different methods, we use the probability of detection (PD) curves here. PD is a function of input signal-to-noise ratio (SNR) with a constant probability of false alarm P f a . The detection curves will be acquired by the cell average constant false alarm rate (CA-CFAR) detector. The CUTs in the area of the 8-12th Doppler bins and the 301-500th range cells are selected as the 1000 runs in the Monte-Carlo experiments. We set that P f a is 10 −4 .
In Figure 9, we show PD versus SNR with ADC, SR-RBC, KASPICE, and the proposed method. We can see that, compared with the detection performance of the KASPICE, there is a slight advantage of about 1 dB in that corresponding to the proposed method. Compared with the SR-RBC and ADC methods, the proposed method demonstrates great detection performance. We can conclude that the result obtained from the probability of detection performance is similar to that obtained from the performance of STAP. This can further prove the superiority of the proposed method.

Conclusions
In this paper, we have presented a novel KASTAP method based on modified sparse learning via iterative minimization named KA-MSLIM-STAP to suppress clutter for a conformal array radar. The new method reformulates the clutter model by introducing the influence of polarization and mutual coupling. Then, a new sparse learning via iterative minimization is derived with Laplace prior. The proposed method exploits the knowledge of array geometry and radar system parameters to construct the dictionary matrix and

Conclusions
In this paper, we have presented a novel KASTAP method based on modified sparse learning via iterative minimization named KA-MSLIM-STAP to suppress clutter for a conformal array radar. The new method reformulates the clutter model by introducing the influence of polarization and mutual coupling. Then, a new sparse learning via iterative minimization is derived with Laplace prior. The proposed method exploits the knowledge of array geometry and radar system parameters to construct the dictionary matrix and utilizes Laplace prior to avoid selecting user parameter. The effectiveness of the proposed method is verified through numerical simulations under ideal and non-ideal conditions (model errors and ICM). The results show the superiority of the proposed method in terms of clutter suppression and target detection. Moreover, we show that the proposed method has lower computational complexity compared with other SR-STAP.
For simplicity, we ignore the influence of the mutual coupling. Hence, this could be the subject of future work.