Sample covariance matrices are widely used in multivariate statistical analysis. The central limit theorems (CLTs) for linear spectral statistics of high-dimensional noncentralized sample covariance ...matrices have received considerable attention in random matrix theory and have been applied to many high-dimensional statistical problems. However, known population mean vectors are assumed for noncentralized sample covariance matrices, some of which even assume Gaussian-like moment conditions. In fact, there are still another two most frequently used sample covariance matrices: the ME (moment estimator, constructed by subtracting the sample mean vector from each sample vector) and the unbiased sample covariance matrix (by changing the denominator n as N = n — 1 in the ME) without depending on unknown population mean vectors. In this paper, we not only establish the new CLTs for noncentralized sample covariance matrices when the Gaussianlike moment conditions do not hold but also characterize the nonnegligible differences among the CLTs for the three classes of high-dimensional sample covariance matrices by establishing a substitution principle: by substituting the adjusted sample size N = n — 1 for the actual sample size n in the centering term of the new CLTs, we obtain the CLT of the unbiased sample covariance matrices. Moreover, it is found that the difference between the CLTs for the ME and unbiased sample covariance matrix is nonnegligible in the centering term although the only difference between two sample covariance matrices is a normalization by n and n — 1, respectively. The new results are applied to two testing problems for high-dimensional covariance matrices.
A fundamental challenge in cognitive radio is the detection of primary users in a licensed spectrum. Cooperative sensing, which utilizes multiple receivers distributed across different locations, ...offers the advantage of utilizing multiple antennas and achieving spatial diversity gain. However, successful implementation of cooperative sensing relies on the ideal exchange of information among cooperating receivers, which may not always be feasible in real-world scenarios. In this paper, we consider the cooperative sensing problem in a non-ideal communication scenario, where the raw data broadcasted from a receiving node can be received by only a subset of the nearby nodes. Existing multiantenna detectors can not deal with such a scenario. To tackle this issue, we propose a novel cooperative sensing scheme, where each node sends only its local correlation coefficients to the fusion center. A detection mechanism based on factorizing the partially received sample covariance matrix is developed. To achieve fast convergence and avoid exhaustive step size tuning, a Bregman proximal method, based on an alternating minimization algorithm (with convergence guarantees), is also developed. The advantages of our proposed cooperative scheme is demonstrated through numerical simulations.
Abstract
Studying the positive definiteness of the covariance matrix of discrete samples helps to determine whether the dimensionality of the samples can be reduced, which is beneficial for ...optimizing the number of samples and designing optimal plans for sampling surveys. This paper aims to provide a method to determine the variable numbers of the sample subjecting to Poisson distribution.
Methods
. It is based on the theory of
I
-linear combination and its properties which are the author’s previous studying results.
Results
. study shows the covariance matrix of multi-Poisson distribution is positively defined and the probability of the sample covariance matrix of multi-poisson distribution is about 1 when the sample capacity is very large.
Conclusion
. The dimension size of the sample data matrix of multi-poisson distribution can be reduced when the sample capacity
n
is no more than the dimension size
p
.
Frequency diverse array and multi-input and multi-output (FDA-MIMO) radar have been applied to many fields due to its angle-range-dependent beampattern. A part of the related array processing ...algorithms can be successfully used on the basis of the known covariance matrix (CM). The estimation accuracy of the CM directly influences the algorithm performance. In particular, the estimation performance of the sample covariance matrix (SCM) will be degraded shapely once the sample size is less than the channel number. Aiming to improve the estimation accuracy of the SCM with FDA-MIMO radar, we propose a novel shrinkage-to-tapering (ST)-based method combined with the block Toeplitz rectification (STT). Firstly, each block matrix in the CM is processed by using the Toeplitz rectification method, and then, three different estimation matrix structures based on the ST method are proposed. Next, by plugging the unbiased estimators into the optimal shrinkage coefficient and normalized mean square error (MSE), their closed forms of their estimators can be given. The numerical simulation results demonstrate that the proposed three different STT approaches are superior to a number of existing methods in terms of CM estimation performance, and the estimation performance is related with the choice of tapering matrix.
In this paper, we study the asymptotic behavior of the extreme eigenvalues and eigenvectors of the high-dimensional spiked sample covariance matrices, in the supercritical case when a reliable ...detection of spikes is possible. In particular, we derive the joint distribution of the extreme eigenvalues and the generalized components of the associated eigenvectors, that is, the projections of the eigenvectors onto arbitrary given direction, assuming that the dimension and sample size are comparably large. In general, the joint distribution is given in terms of linear combinations of finitely many Gaussian and Chi-square variables, with parameters depending on the projection direction and the spikes. Our assumption on the spikes is fully general. First, the strengths of spikes are only required to be slightly above the critical threshold and no upper bound on the strengths is needed. Second, multiple spikes, that is, spikes with the same strength, are allowed. Third, no structural assumption is imposed on the spikes. Thanks to the general setting, we can then apply the results to various high dimensional statistical hypothesis testing problems involving both the eigenvalues and eigenvectors. Specifically, we propose accurate and powerful statistics to conduct hypothesis testing on the principal components. These statistics are data-dependent and adaptive to the underlying true spikes. Numerical simulations also confirm the accuracy and powerfulness of our proposed statistics and illustrate significantly better performance compared to the existing methods in the literature. In particular, our methods are accurate and powerful even when either the spikes are small or the dimension is large.
Covariance matrix tapers have a long history in signal processing and related fields. Examples of applications include autoregressive models (promoting a banded structure) or beamforming (widening ...the spectral null width associated with an interferer). In this paper, the focus is on high-dimensional setting where the dimension <inline-formula><tex-math notation="LaTeX">p</tex-math></inline-formula> is high, while the data aspect ratio <inline-formula><tex-math notation="LaTeX">n/p</tex-math></inline-formula> is low. We propose an estimator called Tabasco (TApered or BAnded Shrinkage COvariance matrix) that shrinks the tapered sample covariance matrix towards a scaled identity matrix. We derive optimal and estimated (data adaptive) regularization parameters that are designed to minimize the mean squared error (MSE) between the proposed shrinkage estimator and the true covariance matrix. These parameters are derived under the general assumption that the data is sampled from an unspecified elliptically symmetric distribution with finite 4th order moments (both real- and complex-valued cases are addressed). Simulation studies show that the proposed Tabasco outperforms all competing tapering covariance matrix estimators in diverse setups. An application to space-time adaptive processing (STAP) also illustrates the benefit of the proposed estimator in a practical signal processing setup.
We study the asymptotic distributions of the spiked eigenvalues and the largest nonspiked eigenvalue of the sample covariance matrix under a general covariance model with divergent spiked ...eigenvalues, while the other eigenvalues are bounded but otherwise arbitrary. The limiting normal distribution for the spiked sample eigenvalues is established. It has distinct features that the asymptotic mean relies on not only the population spikes but also the nonspikes and that the asymptotic variance in general depends on the population eigenvectors. In addition, the limiting Tracy–Widom law for the largest nonspiked sample eigenvalue is obtained.
Estimation of the number of spikes and the convergence of the leading eigenvectors are also considered. The results hold even when the number of the spikes diverges. As a key technical tool, we develop a central limit theorem for a type of random quadratic forms where the random vectors and random matrices involved are dependent. This result can be of independent interest.
A highly popular regularized (shrinkage) covariance matrix estimator is the shrinkage sample covariance matrix (SCM) which shares the same set of eigenvectors as the SCM but shrinks its eigenvalues ...toward the grand mean of the eigenvalues of the SCM. In this paper, a more general approach is considered in which the SCM is replaced by an M-estimator of scatter matrix and a fully automatic data adaptive method to compute the optimal shrinkage parameter with minimum mean squared error is proposed. Our approach permits the use of any weight function such as Gaussian, Huber's, Tyler's, or <inline-formula><tex-math notation="LaTeX">t</tex-math></inline-formula> weight functions, all of which are commonly used in M-estimation framework. Our simulation examples illustrate that shrinkage M-estimators based on the proposed optimal tuning combined with robust weight function do not loose in performance to shrinkage SCM estimator when the data is Gaussian, but provide significantly improved performance when the data is sampled from an unspecified heavy-tailed elliptically symmetric distribution. Also, real-world and synthetic stock market data validate the performance of the proposed method in practical applications.
•Covariance matrices are used in various signal processing tasks.•A cross-validation method for shrinkage covariance matrix estimation is proposed.•The method provides closed-form solutions for ...general covariance matrix estimators.•Applications to array signal processing applications are demonstrated.
Shrinkage can effectively improve the condition number and accuracy of covariance matrix estimation, especially for low-sample-support applications with the number of training samples smaller than the dimensionality. This paper investigates parameter choice for linear shrinkage estimators. We propose data-driven, leave-one-out cross-validation (LOOCV) methods for automatically choosing the shrinkage coefficients, aiming to minimize the Frobenius norm of the estimation error. A quadratic loss is used as the prediction error for LOOCV. The resulting solutions can be found analytically or by solving optimization problems of small sizes and thus have low complexities. Our proposed methods are compared with various existing techniques. We show that the LOOCV method achieves near-oracle performance for shrinkage designs using sample covariance matrix (SCM) and several typical shrinkage targets. Furthermore, the LOOCV method provides low-complexity solutions for estimators that use general shrinkage targets, multiple targets, and/or ordinary least squares (OLS)-based covariance matrix estimation. We also show applications of our proposed techniques to several different problems in array signal processing.