Multichannel adaptive signal detection in space-time colored compound-gaussian autoregressive processes

In this article, we consider the problem of adaptive detection for a multichannel signal in the presence of spatially and temporally colored compound-Gaussian disturbance. By modeling the disturbance as a multichannel autoregressive (AR) process, we first derive a parametric generalized likelihood ratio test against compound-Gaussian disturbance (CG-PGLRT) assuming that the true multichannel AR parameters are perfectly known. For the two-step GLRT design criterion, we combine the multichannel AR parameter estimation algorithm with three covariance matrix estimation strategies for compound-Gaussian environment, then obtain three adaptive CG-PGLRT detectors by replacing the ideal multichannel AR parameters with their estimates. Owing to treating the random texture components of disturbance as deterministic unknown parameters, all of the proposed detectors require no a priori knowledge about the disturbance statistics. The performance assessments are conducted by means of Monte Carlo trials. We focus on the issues of constant false alarm rate (CFAR) behavior, detection and false alarm probabilities. Numerical results show that the proposed adaptive CG-PGLRT detectors have dramatically ease the training and computational burden compared to the generalized likelihood ratio test-linear quadratic (GLRT-LQ) which is referred to as covariance matrix based detector and relies more heavily on training.


Introduction
In an airborne radar system, space-time adaptive processing (STAP) has been widely used in radar target detection; see [1][2][3][4] and references therein. Various wellknown STAP based detectors have been extensively investigated under Gaussian assumption [4][5][6][7]. However, with the support of measured data, the Gaussian model is no longer suitable for background disturbance in many situations of practical interest. The conventional STAP detectors may suffer severe performance degradation when the disturbance is non-Gaussian. Instead, a compound-Gaussian model can successfully describe the non-Gaussian disturbance as a product of a spatially and temporally "slowly varying" texture and a locally "rapidly varying" Gaussian speckle component [8,9]. The texture component accounts for random power variations over range cells. This model includes the socalled spherically invariant random processes (SIRPs). Working with the compound-Gaussian model, a multitude of adaptive detectors have been studied in the past few years, for instance, the generalized likelihood ratio test-linear quadratic (GLRT-LQ) which was independently derived in [10,11], the detectors with Rao and Wald tests [12], the Bayesian optimum radar detector (BORD) [13], and so forth. Notice that all of the aforementioned STAP detectors proposed in Gaussian and compound-Gaussian environment can be considered as covariance matrix based detectors [8,14]. Implementing these detectors involves estimating and inverting a space-time covariance matrix of the disturbance signal for each cell under test (CUT) utilizing independent and identically distributed (IID) target-free training data (or secondary data). Obviously, when the joint spatial-temporal dimension is large, the training and computational requirements will be quite onerous. Moreover, some practical situations may exacerbate the training data selection and collection problem and limit the amount of appropriate IID training data. The lack of training may lead to ill-conditioned covariance matrix estimate and significant degradation in the covariance matrix based detection procedure.
To overcome the above difficulty caused by large joint spatial-temporal dimension, the structural information about the disturbance space-time covariance matrix can be exploited. More precisely, a multichannel autoregressive (AR) process has been found to be able to model the spatial-temporal correlation of the disturbance efficiently [15][16][17][18]. In [15], based on approximating the disturbance spectrum with a multichannel AR model of low order, a parametric adaptive matched filter (PAMF) for STAP detection was presented for multichannel system in Gaussian environment. The PAMF detector, which has been proved to be equivalent to a parametric Rao detector in [16], has dramatically outperformed the conventional adaptive matched filter (AMF) [6,7] with small training size. Also in Gaussian background, a parametric GLRT [17] and a simplified parametric GLRT [18] have been successively developed by utilizing a parametric model (multichannel AR model) in the GLR principle. Experimental results on simulated and real data show that two parametric GLRT detectors work well with limited or even no range training data [18][19][20]. However, under such conditions, the AMF detector and Kelly's GLRT cannot be implemented. Moreover, application of the multichannel AR model in non-stationary Gaussian clutter for STAP are investigated in [21][22][23].
For the corresponding problem in compound-Gaussian environment, a non-Gaussian parametric adaptive matched filter (NG-PAMF) has been derived in [24]. However, this test involves explicit knowledge of the disturbance statistics, which are not always available. Unlike the NG-PAMF, the normalized parametric adaptive matched filter (NPAMF) reported in [8,14] requires no a priori knowledge about the disturbance statistics. This feature is rather important in real-time operation. However, Michels et al. [14] still combined the multichannel AR identification algorithm with the sample covariance matrix (SCM) even in compound-Gaussian background. The sample matrix is the maximum likelihood estimate (MLE) of the covariance matrix for Gaussian disturbance, but is no longer the MLE for compound-Gaussian disturbance. The covariance matrix estimation in compound-Gaussian environment is generally intractable. Conte et al. have advocated the use of a normalized sample covariance matrix (NSCM) in [25]. Also, considering the texture component as an unknown deterministic quantity, an approximate ML (AML) estimator has been derived by Gini et al. [26].
Motivated by the previous studies, the main purpose of this article is to derive a parametric GLRT (PGLRT) for detecting a multichannel signal in the presence of compound-Gaussian disturbance modeled as a multichannel AR process. Without any knowledge about the disturbance statistics, we resort to a suboptimal GLRT algorithm considering the texture components as unknown deterministic parameters. In further derivation, to get round the difficulty in performing the joint maximization for all the unknown parameters, a twostep GLRT design criterion is adopted in this article. We first derive the model-based parametric GLRT in compound-Gaussian environment (CG-PGLRT), which possesses the perfect knowledge about the multichannel AR parameters. We apply three covariance matrix estimation strategies, i.e., sample covariance matrix (SCM), normalized sample covariance matrix (NSCM) and approximate ML (AML) estimator, to the multichannel AR parameter estimation procedure for estimating the unknown AR coefficient matrices. Then three adaptive versions of the CG-PGLRT detector: CG-PGLRT-SCM, CG-PGLRT-NSCM, and CG-PGLRT-AML are obtained. Finally, the performance assessments are presented. Numerical results indicate that the CG-PGLRT-SCM detector has no texture CFAR property, while all of the CG-PGLRT, CG-PGLRT-NSCM, and CG-PGLRT-AML detectors ensure CFAR property with respect to the texture probability density function (PDF). Compared to the covariance matrix based detector, all of the proposed adaptive detectors can handle the training-limited case and significantly decreased the computation complexity in compound-Gaussian environment.
This article is organized as follows. The problem statement and the signal and disturbance models are presented in Section 2. The parametric GLRT detector for compound-Gaussian environment (CG-PGLRT) and adaptive CG-PGLRT detectors are derived in Section 3. Then the performance assessment of the proposed detectors is displayed in Section 4. Finally, conclusions are given in Section 5.
Notation: Vectors (matrices) are denoted by boldface lower (upper) case letters, all vectors are column ones, superscripts (•) T and (•) H denote transpose and complex conjugate transpose, respectively, C N (μ,R) denotes the multivariate complex Gaussian distribution with mean μ and covariance matrix R. C denotes the complex number field and det{•} takes the determinant of a matrix.
2 Problem statement and signal model

Problem statement
Consider the scene that the radar transmits a coherent train of N pulses and receives the signal with a uniform linear array with J sensors. The received data collected over K range cells is organized in a J × N × K data cube. For the range cell under test (CUT), a binary hypothesis test is applied to the JN-dimensional complex baseband space-time vector of primary data Typically, x 0 contains an unwanted additive disturbance signal with unknown space-time covariance matrix R ∈ C JN×JN and may contain a target signal as with deterministic but unknown complex amplitude, a, and known target space-time steering vector s = [s T (0), s T (1), . . . , s T (N − 1)] T ∈ C JN×1 . The spacetime steering vector takes the form of the Kronecker product of the normalized spatial and temporal steering vectors. K IID complex baseband space-time vectors of target-free training data {x k = d k } K k=1 ∈ C JN×1 exist for assisting the signal detection.
The detection problem at hand can be expressed as the following binary hypotheses test: where The disturbance signals {d k } K k=0 lump clutter, jamming, and thermal noise, and may be correlated in space and time.

Signal model
Herein, a multichannel AR process is applied to model the disturbance signal d k (n) at time n, k = 0, 1, . . . , K: where A H (p) P p=1 denotes the unknown J×J coefficient matrices for AR process d k (n) of known order P, and ε k (n) denotes the J × 1 spatial noise driving vectors that are temporally white but spatially colored.
Here the speckle z k (n) are J × 1 complex, zero-mean, Gaussian vectors with unknown covariance matrix Q. The texture component τ k is a positive random variable over range, but constant over time when it has long temporal coherent. The texture PDF f τ (τ) is defined to be the characteristic PDF of the complex SIRV. Given a specific value of τ k , we have ε k (n)|τ k ∼ C N (0, τ k Q) . The covariance matrix of the driving process is E{ε k (n)ε k (n) H } = E(τ k )Q where the mean value E(τ k ) is also the average disturbance power.
Utilizing the multichannel parametric model allows signal whitening adopting a vector one-step linear prediction error filter (PEF), involving temporal whitening via an inverse moving-average (MA) filter followed by spatial whitening [27]. Thus, the temporally whitened versions of the steering vector, the primary data and the secondary data are obtained as following: It is worth noting that the Pth-order linear prediction coefficients A H (p) P p=1 are identically equal to the AR (P) process coefficients.
For large data records, the first P values {x k (n)} P−1 n=0 , k = 0, 1, . . . , K can be ignored [28]. Thus, the joint conditional PDF or likelihood function fi(x k (0), x k (1), . . . , x k (N -1)|τ k ;a,A H ,Q) under hypothesis H i , i = 0 or 1, can be suitably approximated as where ε k (n), for k = 1, . . . , K, is a temporally whitened form of the training signals given by (6) and (7), whereas ε 0 (n) is given by (4), (5), and (8) with a = 0 under H 0 and a ≠ 0 under H 1 . Notice that the secondary data {x k } K k=1 and the primary data x 0 are independent. Define The final joint conditional PDF is given by where τ = [τ 0 , τ 1 , . . . , τ K ] T and is expressed in terms of the Hermitian operation for notational convenience, but is not necessarily a Hermitian matrix.

Parametric GLRT derivation in compoundgaussian environment
In the Neyman-Pearson sense, the optimum solution for the composite hypothesis testing problem (1) is the likelihood ratio test (LRT). But for the case at hand, it cannot be implemented due to total ignorance of the signal parameter a, the multichannel AR parameters Q and A H and the texture PDF f τ (τ). Hence, we resort to a suboptimal GLRT algorithm where the τ k s are modeled as unknown deterministic parameters, and perform the maximum likelihood estimation (MLE) for all the unknown parameters under each hypothesis. Unfortunately, the exact maximization with respect to the unknown parameters is rather difficult and does not exist a close-form expression. Therefore, to get round the above difficulty, the two-step GLRT design criterion is adopted. We first assume that the AR parameters Q and A H are perfectly known, and derive the CG-PGLRT detector based on the primary data. The adaptive versions of CG-PGLRT detector are then obtained by substituting the unknown AR parameters with their estimates based on the training signals only.

CG-PGLRT detector for known AR parameters
The GLRT for the case at hand, under the assumption that Q and A H are known as well as τ k s are modeled as unknown deterministic parameters, is given by where h is the threshold value to be set in order to ensure the desired probability of false alarm (P fa ). Under H 1 and H 0 , the joint conditional PDFs of the primary data are and respectively. The ML estimates of τ 0 are easily derived, from the definition wherẽ X 0 = [x 0 (P),x 0 (P + 1), . . . ,x 0 (N − 1)] ∈ C J×(N−P) (17) and tr{·} denotes the trace operator. Utilizing the two estimators in (14) and (15), the joint conditional PDFs under H 1 and H 0 , respectively, are: Performing the maximization of (18) with respect to a is tantamount to minimizing the following expression: When the positive factor containing a is made to vanish, the minimum is clearly attained Substitutingα as well asτ 0 0 andτ 0 1 into the decision rule (11), after some algebra manipulations, the detection statistic for parametric GLRT in compound-Gaussian environment (CG-PGLRT) can be recast as where h CG-PGLRT is the appropriate modification of the original threshold in (11).
Note that the normalized parametric adaptive matched filter (NPAMF), originally developed in [8,14] for compound-Gaussian environment, is closely related to the CG-PGLRT detector but replaces true Q and A H with their estimates obtained from the training signals. In [8,14,29], several multichannel parameter estimation algorithms are considered in the NPAMF detector, including the Nuttall-Strand and the multichannel least squares methods [30]. However, we obtain herein adaptive version of the CG-PGLRT detector relying on maximum likelihood (ML) parameter estimation criterion. Specifically, our multichannel parameter estimation approach is similar to that proposed in [16,17] for Gaussian parametric model-based STAP detectors, whereas we utilize only training signals for parameter estimation.

Adaptive CG-PGLRT detector for unknown AR parameters
To make the derived detector (22) fully adaptive, suitable estimates of the multichannel AR parameters Q and A H must be acquired. Employing the target-free training signals to estimate the multichannel AR parameters allows to decouple from estimation of target parameters, such as complex amplitude a and the texture component τ 0 . However, derivation of the MLEs of Q and A H with K unknown textures {τ k } K k=1 is still a more challenging task. The similar covariance estimation problem in compound-Gaussian environment has been solved in [25,26,31], where three different estimation strategies are introduced, including SCM, NSCM, and AML estimators.
In our opinion, the derivation of these three covariance estimators depends upon two different clutter models, that is, the dependent interference model [31] and the independent interference model [26], respectively. The former assumes that the textures of training signals are completely correlated, i.e., τ 1 = τ 2 = ··· = τ K , and adopts the SCM estimator. The latter, where the NSCM and AML estimators are usually applied, deals with the IID {τ k } K k=1 instead. For our problem, as suggested in [26], a two-step person-by-person maximization is followed. Specifically, we first assume that the textures of the training signals are perfectly known for the realization under observation.
. . , K , the log-likelihood function based on training signals is proportional to where ε k (n) =x k (n) and D τ t ∈ C K×K is a diagonal matrix whose diagonal entries are τ 1 , . . . , T K . We, next, use (23) to derive the MLEs of Q and A H .
By setting to zero the derivative of (23) with respect to Q, we get the MLE of Q given τ 1 = [τ 1 , . . . , τ K ] T and (23), we find that maximizing (23) with respect We further observe that [17] where the above correlation matrices conditioned on andR xy (τ t ) =R H yx (τ t ) . SinceR yy (τ t ) is non-negative definite and the second and third terms in (25) do not depend on A H , it follows that [17] wherê Inserting (30) into (25), we finally havê Obviously, once the MLE of τ t = [τ 1 , . . . , τ K ] T , namelŷ τ tML = τ 1ML , . . . ,τ KML T , is available, the exact MLEs of A H and Q can be readily obtained according to (30) and (31). However, so far the explicit solutions forτ tML are still inaccessible. Therefore, in the second step of the person-by-person maximization, we introduce the aforementioned two clutter models to motivate the use of the three different covariance estimation strategies.
For the dependent interference model, it is interesting that the MLE of Q with identical textures is proportional to the sample covariance matrix estimate of the temporally whitened training signals within a scalar multiplicative factor l [26,31]. More precisely, based on the classical sample covariance matrix (SCM) estimator Q ML ,Q ML for the dependent interference model is given aŝ Likewise, following the derivation of (25), (30), and (31), the MLEs of A H and Q for the dependent model can be straightforwardly expressed aŝ wherê Notice that an explicit solution for λ is not necessary since the AR coefficients estimateÂ Thus, we havễ S = [ŝ(P),ŝ(P + 1), . . . ,ŝ(N − 1)] ∈ C J×(N−P) and PluggingQ SCM ,Ŝ , andX 0 in place of Q,S, andX 0 into (22), respectively, leads to our adaptive CG-PGLRT-SCM detector given by Actually, the dependent interference model can be considered as a simple extension of the partially-homogeneous environment under the assumption of the deterministic texture component. Recently, two different parametric Rao tests, referred to as the normalized parametric Rao (NPRao) test and the scale-invariant parametric Rao (SI-PRao) test, respectively, have been developed in [32] for the partially-homogeneous environment. The above parametric Rao tests can be candidates for multichannel parametric detection in the dependent interference model. Different to the CG-PGLRT-SCM detector, the NPRao and SI-PRao tests use both training and test signals for parameter estimation. However, the above techniques are no longer applicable to the independent interference model because of the model mismatch.
In the case of the independent interference model, the texture PDF f τ (τ) is not always known a prior. To cope with this uncertainty, the random clutter texture component is treated as a deterministic and unknown parameter. The first candidate for the independent model is the normalized sample covariance matrix (NSCM) estimator [26,33]. Here, we normalize each term {x k (n)x H k (n)} in the sums in (24) by the data-dependent normalization factorτ k NSCM in place of true τ k to obtain is the sample estimate of the local disturbance power in the reference range cells.
According to (30) and (31), we can also obtainÂ where the covariance matrices are also normalized by From (41) and (42), it is worth observing that an interdependent relationship exists between the estimation ofτ k NSCM andÂ H NSCM . Therefore, it is rather difficult to perform joint estimation for them. The idea pursued here is to obtainτ k NSCM by utilizingÂ H SCM , namely, we getτ In like manner, the adaptive CG-PGLRT-NSCM detector is given by where the matrices S ∈ C J×(N−P) and X 0 ∈ C J×(N−P) consist of the temporally whitened versions of {x(n)} N−1 n=P and {s(n)} N−1 n=P , respectively, using the coefficient estimateÂ H NSCM . The second candidate for the independent model is the constrained approximate ML (AML) estimator, provided by the ML theory. Summarily, the constrained AML estimates of Q and A H are derived by jointly solving iteratively the following equations: for i = 0, 1, 2, . . . , N it , where N it is the number of iterations. Equation (52) guarantees the constraint tr{Q AML (i + 1)} = tr{ Q} to be satisfied at each iteration. Note that the covariance matricesR xx ,R yy , andR yx have the same structure as those in (44)-(46), the only difference being that the normalization factorτ k NSCM is replaced byτ k AML (i + 1) which is refreshed at each itera- In general, the adaptive versions of the CG-PGLRT detector can be defined as a parametric implementation of the following adaptive GLRT-LQ detector [34] GLRT -LQ = −JN ln where the original GLRT-LQ detector given in [10,11] has been reintroduced in (54). The space-time covariance matrix estimateR can be also obtained by exploiting the SCM, NSCM, and AML estimators, respectively. However, taking the SCM estimator as an example, we at least need K ≥ JN training signals to ensure a fullrank estimate of the JN × JN matrix R. Obviously, the large JN spatio-temporal product will impose excessive training and computational burdens to the detector.

Complexity issues
We provide the floating-point operations (flops) involved in each step of the parametric and non-parametric STAP detector with different covariance estimators in Tables 1 and 2, respectively. Only the complex multiplication and division are considered but the complex addition and subtraction are ignored.
Note that the adaptive parametric implementations differ mainly in parameter estimation, they share identical steps in signal whitening and calculating the test statistic. The CG-PGLRT-NSCM detector is slightly more complex than the CG-PGLRT-SCM detector since it requires evaluatingτ k NSCM by utilizingÂ H SCM for normalization. However, the CG-PGLRT-AML detector is the most complex detector among the three. Likewise, similar conclusion can be made for the adaptive versions of the GLRT-LQ detector. For a quick comparison, suppose KN >JP for the parametric detectors in Table 1. Clearly, it can be seen from Tables 1 and 2 that the parametric detectors can reduce the computational complexity of their non-parametric counterparts, especially when the spatial-temporal dimension JN is large.

Performance assessment
This section is devoted to the performance analysis of the proposed detectors: CG-PGLRT, CG-PGLRT-SCM, CG-PGLRT-NSCM, and CG-PGLRT-AML in terms of the probability of false alarm (P fa ) and the probability of detection (P d ). For compound-Gaussian environment, the closed-form expressions for both P fa as well as P d are not available. Hence, we carry out the analysis via Monte-Carlo techniques based on 100/P fa and 100/P d independent trials, respectively. In order to limit the computational burden, we set the probability of false alarm P fa = 10 -3 and the number of sensors J = 4 throughout the section.
Moreover, the following statements hold for this section: (i) The spatial driving noise is K-distributed, which is obtained assuming that the texture component is Gamma distributed [35] where Г(·) is the gamma function, μ = E{τ} denotes the mean of the distribution, and ν is the shape parameter which provides a measure of noise spikiness. Lower ν means more spikes will appear in disturbance. Without loss of generality, we set μ = 1.
(ii) The disturbance signal is generated as a multichannel AR(2) compound-Gaussian process with the method in [36,37]. We assume that the order P of the multichannel AR process is known. If P is unknown, it can be selected using some model order selection criterion [38]. The following correlation matrices are utilized to synthesize simulated data [39]  where ρ relates to the disturbance power, ρ s decides the spatial correlation, ρ t controls the temporal correlation, and {(mn + l)ω} defines the phase of the correlation function (56). Here, we choose the value of ρ to satisfy tr{Q} = J. For reader's convenience, we give an overview of relevant parameter settings for the simulation processes in Table 3.
In Figure 6, the effects of the disturbance spikiness upon the detection performance of the CG-PGLRT detector are analyzed by varying the shape parameter ν. Figure 6 displays P d plots versus SINR for shape parameters ν = 0.1, 0.5, 1, 4.5, 10. The figure highlights that the detection performance is enhanced as ν decreases at the low SINR region. However, at the large SINR region (about > 10 dB), the CG-PGLRT detector performs worse with smaller ν. The lower the value of ν and the spikier the disturbance, for a given average disturbance  power μ, at low SINR (small value of a), the signal energy that can be accumulated from the cells under test increases, hence, the detection performance is improved; at large SINR (large value of a), more spikes in disturbance may have a negative effect on the signal detection. In addition, the detector exhibits close performance for ν = 4.5 and 10, since the K distribution collapses into the Gaussian distribution when ν ≥ 4.5. In Section 3.2, we combined the multichannel AR parameter estimation algorithm with three covariance matrix estimation strategies: SCM, NSCM, and AML estimators, and then gave three adaptive CG-PGLRT detectors, where the true multichannel AR parameters Q and A H are substituted with the estimated ones. Now various simulations are performed to evaluate the detection performance loss of the adaptive CG-PGLRT detectors with respect to the CG-PGLRT detector which possesses the perfect knowledge about the multichannel AR parameters Q and A H .
In Figures 7, 8, and 9, the limited-training (K = 2) case is considered for different numbers of pulses N (N = 16, 32) and shape parameters ν = 0.1, 0.5 and 4.5. In particular, Figure 7 refers to the shape parameter ν = 0.5 whereas Figures 8 and 9 to ν = 0.1 and ν = 4.5. All the detection curves in Figures 7, 8, and 9 show that the detection performance is improved as N increases. For instance, at P d = 0.9, the CG-PGLRT-AML is about 4 dB from the CG-PGLRT bound at N = 16; the gap reduces to about 1.5 dB at N = 32 for any value of ν in this case. This implies that we can remedy the performance degradation with respect to CG-PGLRT detector by increasing the temporal dimension N. In Figures 7  and 8, we find that the CG-PGLRT-NSCM bears almost the same performance as the CG-PGLRT-AML at N = 32 and slightly poorer than the CG-PGLRT-AML at N = 16. Figure 9 demonstrates that all the adaptive detectors: CG-PGLRT-SCM, CG-PGLRT-NSCM, and CG-PGLRT-AML perform almost the same in the homogeneous environment (ν ≥ 4.5), simply because the K distribution converges to the Gaussian distribution for increasing ν. Thus, the loss of the CG-PGLRT detector with the SCM estimator due to model mismatch decreases as ν increases. We also observe that the CG-PGLRT-AML detector outperforms the other two detectors for the heavy tailed disturbance (ν 0). Thus, the CG-PGLRT-AML detector is the most suitable detector to implement the adaptive detection in the realistic spiky disturbance environment.  In Figures 10 and 11, we compare the performance of the proposed adaptive CG-PGLRT detectors with that of the GLRT-LQ with R scm and with R nscm . Figures 10 and  11 refers to the plots of P d versus SINR of the adaptive CG-PGLRT detectors with training size K = 2, 3 and the GLRT-LQ with R scm (R nscm ) with K = 128. The large K chosen for the GLRT-LQ with estimated R is to ensure a nonsingular estimate of R and obtain performance within 3 dB from the optimum bound. It is observed that, with the limited-training (K = 2), the performance of the CG-PGLRT-SCM is close to that of the GLRT-LQ with R scm . And both CG-PGLRT-NSCM and CG-PGLRT-AML detectors can produce better detection performance with dramatically less training than the GLRT-LQ with R scm . Though all the proposed adaptive CG-PGLRT detectors with K = 2 are inferior to the GLRT-LQ with R nscm with K = 128, the performance degradation of the CG-PGLRT-NSCM and CG-PGLRT-AML detectors can be easily remedied only by increasing the training size to K = 3. Indeed, the parametric detectors are significantly less dependent on training than the covariance matrix based detector. So far the model order P of the multichannel AR process has been assumed to be known. In the following simulation example, we evaluate the performance of the proposed adaptive CG-PGLRT detectors when an model order estimation error occurs. Figure 12 depicts the performance of the proposed adaptive CG-PGLRT detectors using an overestimated model order. The curves indicate that some performance loss caused by mismatched model order P happens for N = 16 and K = 3. Since P samples have been discarded in the derivation of the proposed detectors as well as in calculating the test statistics, the detection performance is degraded as model order P is overestimated and the number of pulses N is not large enough.

Conclusions
In this article, multichannel signal detection problem in space-time colored compound-Gaussian environment is discussed. By exploiting the structural information about the disturbance space-time covariance matrix, we model the disturbance signal as a multichannel AR process to ease the training and computational burdens. Modeling the texture as an unknown deterministic parameter, we first derive the CG-PGLRT detector under the assumption that the multichannel AR parameters Q and A H are perfectly known. For the two-step GLRT design criterion, we combine the multichannel AR parameter estimation algorithm with three covariance matrix estimation strategies, i.e., sample covariance matrix (SCM), normalized sample covariance matrix (NSCM) and approximate ML (AML) estimator, and then obtain the adaptive versions of the CG-PGLRT detector by substituting the true multichannel AR parameters with their estimates. Finally, we show the CFAR behavior and detection performances of the proposed detectors: CG-PGLRT, CG-PGLRT-SCM, CG-PGLRT-NSCM, and CG-PGLRT-AML by Monte Carlo trials.
Conclusions of this article are summarized as follows: -The CG-PGLRT-SCM detector has no texture CFAR property, while the CG-PGLRT, CG-PGLRT-NSCM, and CG-PGLRT-AML detectors ensure texture CFAR property.
-The detection probability of the CG-PGLRT detector increases with increasing disturbance spikes (decreasing ν) at low SINR. However, at large SINR, we observe a degradation in performance as ν decreases.
-For the limited-training case, the detection performance loss of the adaptive CG-PGLRT detectors with respect to the CG-PGLRT detector can be remedied by increasing temporal dimension N.
-The CG-PGLRT-AML detector has the best detection performance in the heavy tailed disturbance (ν 0). Thus, the CG-PGLRT-AML detector is the most suitable detector to implement the adaptive detection in the realistic spiky disturbance. -Compared to the covariance matrix based detector, the proposed model-based adaptive detectors have significantly decreased the training requirements and the computation complexity.