Spatial covariance matrix(SCM) is essential in many multi-antenna systems such as massive multiple-input multiple-output(MIMO). For multi-antenna systems operating at millimeter-wave bands, hybrid analog-digital struc...Spatial covariance matrix(SCM) is essential in many multi-antenna systems such as massive multiple-input multiple-output(MIMO). For multi-antenna systems operating at millimeter-wave bands, hybrid analog-digital structure has been widely adopted to reduce the cost of radio frequency chains.In this situation, signals received at the antennas are unavailable to the digital receiver, and as a consequence, traditional sample average approach cannot be used for SCM reconstruction in hybrid multi-antenna systems. To address this issue, beam sweeping algorithm(BSA) which can reconstruct the SCM effectively for a hybrid uniform linear array, has been proposed in our previous works. However, direct extension of BSA to a hybrid uniform circular array(UCA)will result in a huge computational burden. To this end, a low-complexity approach is proposed in this paper. By exploiting the symmetry features of SCM for the UCA, the number of unknowns can be reduced significantly and thus the complexity of reconstruction can be saved accordingly. Furthermore, an insightful analysis is also presented in this paper, showing that the reduction of the number of unknowns can also improve the accuracy of the reconstructed SCM. Simulation results are also shown to demonstrate the proposed approach.展开更多
The paper discusses the statistical inference problem of the compound Poisson vector process(CPVP)in the domain of attraction of normal law but with infinite covariance matrix.The empirical likelihood(EL)method to con...The paper discusses the statistical inference problem of the compound Poisson vector process(CPVP)in the domain of attraction of normal law but with infinite covariance matrix.The empirical likelihood(EL)method to construct confidence regions for the mean vector has been proposed.It is a generalization from the finite second-order moments to the infinite second-order moments in the domain of attraction of normal law.The log-empirical likelihood ratio statistic for the average number of the CPVP converges to F distribution in distribution when the population is in the domain of attraction of normal law but has infinite covariance matrix.Some simulation results are proposed to illustrate the method of the paper.展开更多
We normalize data from 43 Chinese observatories and select data from ten Chinese observatories with most continuous records to assess the secular variations(SVs)and geomagnetic jerks by calculating the deviations betw...We normalize data from 43 Chinese observatories and select data from ten Chinese observatories with most continuous records to assess the secular variations(SVs)and geomagnetic jerks by calculating the deviations between annual observed and CHAOS-6 model monthly means.The variations in the north,east,and vertical eigendirections are studied by using the covariance matrix of the residuals,and we find that the vertical direction is strongly affected by magnetospheric ring currents.To obtain noise-free data,we rely on the covariance matrix of the residuals to remove the noise contributions from the largest eigenvalue or vectors owing to ring currents.Finally,we compare the data from the ten Chinese observatories to seven European observatories.Clearly,the covariance matrix method can simulate the SVs of Dst,the jerk of the northward component in 2014 and that of the eastward component in 2003.5 in China are highly agree with that of Vertically downward component in Europe,compare to CHAOS-6,covariance matrix method can show more details of SVs.展开更多
Most of the reconstruction-based robust adaptive beamforming(RAB)algorithms require the covariance matrix reconstruction(CMR)by high-complexity integral computation.A Gauss-Legendre quadrature(GLQ)method with the high...Most of the reconstruction-based robust adaptive beamforming(RAB)algorithms require the covariance matrix reconstruction(CMR)by high-complexity integral computation.A Gauss-Legendre quadrature(GLQ)method with the highest algebraic precision in the interpolation-type quadrature is proposed to reduce the complexity.The interference angular sector in RAB is regarded as the GLQ integral range,and the zeros of the threeorder Legendre orthogonal polynomial is selected as the GLQ nodes.Consequently,the CMR can be efficiently obtained by simple summation with respect to the three GLQ nodes without integral.The new method has significantly reduced the complexity as compared to most state-of-the-art reconstruction-based RAB techniques,and it is able to provide the similar performance close to the optimal.These advantages are verified by numerical simulations.展开更多
An improved two-channel Synthetic Aperture Radar Ground Moving Target Indication (SAR-GMTI) method based on eigen-decomposition of the covariance matrix is investigated. Based on the joint Probability Density Function...An improved two-channel Synthetic Aperture Radar Ground Moving Target Indication (SAR-GMTI) method based on eigen-decomposition of the covariance matrix is investigated. Based on the joint Probability Density Function (PDF) of the Along-Track Interferometric (ATI) phase and the similarity between the two SAR complex images, a novel ellipse detector is presented and is applied to the indication of ground moving targets. We derive its statistics and analyze the performance of detection process in detail. Compared with the approach using the ATI phase, the ellipse detector has a better performance of detection in homogenous clutter. Numerical experiments on simulated data are presented to validate the improved performance of the ellipse detector with respect to the ATI phase approach. Finally, the detection capability of the proposed method is demonstrated by measured SAR data.展开更多
In this article, the problem of estimating the covariance matrix in general linear mixed models is considered. Two new classes of estimators obtained by shrinking the eigenvalues towards the origin and the arithmetic ...In this article, the problem of estimating the covariance matrix in general linear mixed models is considered. Two new classes of estimators obtained by shrinking the eigenvalues towards the origin and the arithmetic mean, respectively, are proposed. It is shown that these new estimators dominate the unbiased estimator under the squared error loss function. Finally, some simulation results to compare the performance of the proposed estimators with that of the unbiased estimator are reported. The simulation results indicate that these new shrinkage estimators provide a substantial improvement in risk under most situations.展开更多
Multivariate seemingly unrelated regression system is raised first and the two stage estimation and its covariance matrix are given. The results of the literatures[1-5] are extended in this paper.
Differential evolution algorithm based on the covariance matrix learning can adjust the coordinate system according to the characteristics of the population, which make<span style="font-family:Verdana;"&g...Differential evolution algorithm based on the covariance matrix learning can adjust the coordinate system according to the characteristics of the population, which make<span style="font-family:Verdana;">s</span><span style="font-family:Verdana;"> the search move in a more favorable direction. In order to obtain more accurate information about the function shape, this paper propose</span><span style="font-family:Verdana;">s</span><span style="font-family:;" "=""> <span style="font-family:Verdana;">covariance</span><span style="font-family:Verdana;"> matrix learning differential evolution algorithm based on correlation (denoted as RCLDE)</span></span><span style="font-family:;" "=""> </span><span style="font-family:Verdana;">to improve the search efficiency of the algorithm. First, a hybrid mutation strategy is designed to balance the diversity and convergence of the population;secondly, the covariance learning matrix is constructed by selecting the individual with the less correlation;then, a comprehensive learning mechanism is comprehensively designed by two covariance matrix learning mechanisms based on the principle of probability. Finally,</span><span style="font-family:;" "=""> </span><span style="font-family:;" "=""><span style="font-family:Verdana;">the algorithm is tested on the CEC2005, and the experimental results are compared with other effective differential evolution algorithms. The experimental results show that the algorithm proposed in this paper is </span><span style="font-family:Verdana;">an effective algorithm</span><span style="font-family:Verdana;">.</span></span>展开更多
Regularized system identification has become the research frontier of system identification in the past decade.One related core subject is to study the convergence properties of various hyper-parameter estimators as t...Regularized system identification has become the research frontier of system identification in the past decade.One related core subject is to study the convergence properties of various hyper-parameter estimators as the sample size goes to infinity.In this paper,we consider one commonly used hyper-parameter estimator,the empirical Bayes(EB).Its convergence in distribution has been studied,and the explicit expression of the covariance matrix of its limiting distribution has been given.However,what we are truly interested in are factors contained in the covariance matrix of the EB hyper-parameter estimator,and then,the convergence of its covariance matrix to that of its limiting distribution is required.In general,the convergence in distribution of a sequence of random variables does not necessarily guarantee the convergence of its covariance matrix.Thus,the derivation of such convergence is a necessary complement to our theoretical analysis about factors that influence the convergence properties of the EB hyper-parameter estimator.In this paper,we consider the regularized finite impulse response(FIR)model estimation with deterministic inputs,and show that the covariance matrix of the EB hyper-parameter estimator converges to that of its limiting distribution.Moreover,we run numerical simulations to demonstrate the efficacy of ourtheoretical results.展开更多
Multiphoton microscopy is the enabling tool for biomedical research,but the aberrations of biological tissues have limited its imaging performance.Adaptive optics(AO)has been developed to partially overcome aberration...Multiphoton microscopy is the enabling tool for biomedical research,but the aberrations of biological tissues have limited its imaging performance.Adaptive optics(AO)has been developed to partially overcome aberration to restore imaging performance.For indirect AO,algorithm is the key to its successful implementation.Here,based on the fact that indirect AO has an analogy to the black-box optimization problem,we successfully apply the covariance matrix adaptation evolution strategy(CMA-ES)used in the latter,to indirect AO in multiphoton microscopy(MPM).Compared with the traditional genetic algorithm(GA),our algorithm has a greater improvement in convergence speed and convergence accuracy,which provides the possibility of realizing real-time dynamic aberration correction for deep in vivo biological tissues.展开更多
The estimation of covariance matrices is very important in many fields, such as statistics. In real applications, data are frequently influenced by high dimensions and noise. However, most relevant studies are based o...The estimation of covariance matrices is very important in many fields, such as statistics. In real applications, data are frequently influenced by high dimensions and noise. However, most relevant studies are based on complete data. This paper studies the optimal estimation of high-dimensional covariance matrices based on missing and noisy sample under the norm. First, the model with sub-Gaussian additive noise is presented. The generalized sample covariance is then modified to define a hard thresholding estimator , and the minimax upper bound is derived. After that, the minimax lower bound is derived, and it is concluded that the estimator presented in this article is rate-optimal. Finally, numerical simulation analysis is performed. The result shows that for missing samples with sub-Gaussian noise, if the true covariance matrix is sparse, the hard thresholding estimator outperforms the traditional estimate method.展开更多
In this paper,a new correlated covariance matrix for Multi-Input Multi-Output(MIMO)radar is proposed,which has lower Side Lobe Levels(SLLs)compared to the new covariance matrix designs and the well-known multi-antenna...In this paper,a new correlated covariance matrix for Multi-Input Multi-Output(MIMO)radar is proposed,which has lower Side Lobe Levels(SLLs)compared to the new covariance matrix designs and the well-known multi-antenna radar designs including phased-array,MIMO radar and phased-MIMO radar schemes.It is shown that Binary Phased-Shift Keying(BPSK)waveforms that have constant envelope can be used in a closed-form to realize the proposed covariance matrix.Therefore,there is no need to deploy different types of radio amplifiers in the transmitter which will reduce the cost,considerably.The proposed design allows the same transmit power from each antenna in contrast to the phased-MIMO radar.Moreover,the proposed covariance matrix is full-rank and has the same capability as MIMO radar to identify more targets,simultaneously.Performance of the proposed transmit covariance matrix including receive beampattern and output Signal-to-Interference plus Noise Ratio(SINR)is simulated,which validates analytical results.展开更多
This paper aims at achieving a simultaneously sparse and low-rank estimator from the semidefinite population covariance matrices.We first benefit from a convex optimization which develops l1-norm penalty to encourage ...This paper aims at achieving a simultaneously sparse and low-rank estimator from the semidefinite population covariance matrices.We first benefit from a convex optimization which develops l1-norm penalty to encourage the sparsity and nuclear norm to favor the low-rank property.For the proposed estimator,we then prove that with high probability,the Frobenius norm of the estimation rate can be of order O(√((slgg p)/n))under a mild case,where s and p denote the number of nonzero entries and the dimension of the population covariance,respectively and n notes the sample capacity.Finally,an efficient alternating direction method of multipliers with global convergence is proposed to tackle this problem,and merits of the approach are also illustrated by practicing numerical simulations.展开更多
A two-mode entangled state was generated experimentally through mixing two squeezed lights from two optical parametric amplifiers on a 50/50 beam splitter.The entangled beams were measured by means of two pairs of bal...A two-mode entangled state was generated experimentally through mixing two squeezed lights from two optical parametric amplifiers on a 50/50 beam splitter.The entangled beams were measured by means of two pairs of balanced homodyne detection systems respectively.The relative phases between the local beams and the detected beams can be locked by using the optical phase modulation technique.The covariance matrix of the two-mode entangled state was obtained when the relative phase of the local beam and the detected beam in one homodyne detection system is locked and the other is scanned.This method provides a way by which one can extract the covariance matrix of any selected quadrature components of two-mode Gaussian state.展开更多
This paper is concerned with tile proOlenl or improving hue ~lma^u~ u~ under Stein＇s loss. By the partial Iwasawa coordinates of covariance matrix, the corresponding risk can be split into three parts. One can use th...This paper is concerned with tile proOlenl or improving hue ~lma^u~ u~ under Stein＇s loss. By the partial Iwasawa coordinates of covariance matrix, the corresponding risk can be split into three parts. One can use the information in the weighted matrix of weighted quadratic loss to improve one part of risk. However, this paper indirectly takes advantage of the information in the sample mean and reuses Iwasawa coordinates to improve the rest of risk. It is worth mentioning that the process above can be repeated. Finally, a Monte Carlo simulation study is carried out to verify the theoretical results.展开更多
In this paper, the problem of estimating the covariance matrix in general linear mixed models is considered. A new class of estimators is proposed. It is shown that this new estimator dominates the analysis of varianc...In this paper, the problem of estimating the covariance matrix in general linear mixed models is considered. A new class of estimators is proposed. It is shown that this new estimator dominates the analysis of variance estimate under two squared loss functions. Finally, some simulation results to compare the performance of the proposed estimator with that of the analysis of variance estimate are reported. The simulation results indicate that this new estimator provides a substantial improvement in risk under most situations.展开更多
Covariance matrix plays an important role in risk management, asset pricing, and portfolio allocation. Covariance matrix estimation becomes challenging when the dimensionality is comparable or much larger than the sam...Covariance matrix plays an important role in risk management, asset pricing, and portfolio allocation. Covariance matrix estimation becomes challenging when the dimensionality is comparable or much larger than the sample size. A widely used approach for reducing dimensionality is based on multi-factor models. Although it has been well studied and quite successful in many applications, the quality of the estimated covariance matrix is often degraded due to a nontrivial amount of missing data in the factor matrix for both technical and cost reasons. Since the factor matrix is only approximately low rank or even has full rank, existing matrix completion algorithms are not applicable. We consider a new matrix completion paradigm using the factor models directly and apply the alternating direction method of multipliers for the recovery. Numerical experiments show that the nuclear-norm matrix completion approaches are not suitable but our proposed models and algorithms are promising.展开更多
Consider the problems of frequency-invariant beampattern optimization and robustness in broadband beamforming.Firstly,a global optimization algorithm,which is based on phase compensation of the array manifolds,is used...Consider the problems of frequency-invariant beampattern optimization and robustness in broadband beamforming.Firstly,a global optimization algorithm,which is based on phase compensation of the array manifolds,is used to construct the frequency-invariant beampattern.Compared with some methods presented recently,the proposed algorithm is not only available to get the global optimal solution,but also simple for physical realization.Meanwhile,a robust adaptive broadband beamforming algorithm is also derived by reconstructing the covariance matrix.The essence of the proposed algorithm is to estimate the space-frequency spectrum using Capon estimator firstly,then integrate over a region separated from the desired signal direction to reconstruct the interference-plus-noise covariance matrix,and finally caleulate the adaptive beamformer weights with the reconstructed matrix.The design of beamformer is formulated as a convex optimization problem to be solved.Simulation results show that the performance of the proposed algorithm is almost always close to the optimal value across a wide range of signal to noise ratios.展开更多
This paper addresses the issue of testing sphericity and identity of high-dimensional population covariance matrix when the data dimension exceeds the sample size.The central limit theorem of the first four moments of...This paper addresses the issue of testing sphericity and identity of high-dimensional population covariance matrix when the data dimension exceeds the sample size.The central limit theorem of the first four moments of eigenvalues of sample covariance matrix is derived using random matrix theory for generally distributed populations.Further,some desirable asymptotic properties of the proposed test statistics are provided under the null hypothesis as data dimension and sample size both tend to infinity.Simulations show that the proposed tests have a greater power than existing methods for the spiked covariance model.展开更多
基金supported by National Key Research and Development Program of China under Grant 2020YFB1804901State Key Laboratory of Rail Traffic Control and Safety(Contract:No.RCS2022ZT 015)Special Key Project of Technological Innovation and Application Development of Chongqing Science and Technology Bureau(cstc2019jscx-fxydX0053).
文摘Spatial covariance matrix(SCM) is essential in many multi-antenna systems such as massive multiple-input multiple-output(MIMO). For multi-antenna systems operating at millimeter-wave bands, hybrid analog-digital structure has been widely adopted to reduce the cost of radio frequency chains.In this situation, signals received at the antennas are unavailable to the digital receiver, and as a consequence, traditional sample average approach cannot be used for SCM reconstruction in hybrid multi-antenna systems. To address this issue, beam sweeping algorithm(BSA) which can reconstruct the SCM effectively for a hybrid uniform linear array, has been proposed in our previous works. However, direct extension of BSA to a hybrid uniform circular array(UCA)will result in a huge computational burden. To this end, a low-complexity approach is proposed in this paper. By exploiting the symmetry features of SCM for the UCA, the number of unknowns can be reduced significantly and thus the complexity of reconstruction can be saved accordingly. Furthermore, an insightful analysis is also presented in this paper, showing that the reduction of the number of unknowns can also improve the accuracy of the reconstructed SCM. Simulation results are also shown to demonstrate the proposed approach.
基金Characteristic Innovation Projects of Ordinary Universities of Guangdong Province,China(No.2022KTSCX150)Zhaoqing Education Development Institute Project,China(No.ZQJYY2021144)Zhaoqing College Quality Project and Teaching Reform Project,China(Nos.zlgc202003 and zlgc202112)。
文摘The paper discusses the statistical inference problem of the compound Poisson vector process(CPVP)in the domain of attraction of normal law but with infinite covariance matrix.The empirical likelihood(EL)method to construct confidence regions for the mean vector has been proposed.It is a generalization from the finite second-order moments to the infinite second-order moments in the domain of attraction of normal law.The log-empirical likelihood ratio statistic for the average number of the CPVP converges to F distribution in distribution when the population is in the domain of attraction of normal law but has infinite covariance matrix.Some simulation results are proposed to illustrate the method of the paper.
基金supported by the National Natural Science Foundation of China(Grant No.41404053)Special Project for Meteo-Scientifi c Research in the Public Interest(No.GYHY201306073)
文摘We normalize data from 43 Chinese observatories and select data from ten Chinese observatories with most continuous records to assess the secular variations(SVs)and geomagnetic jerks by calculating the deviations between annual observed and CHAOS-6 model monthly means.The variations in the north,east,and vertical eigendirections are studied by using the covariance matrix of the residuals,and we find that the vertical direction is strongly affected by magnetospheric ring currents.To obtain noise-free data,we rely on the covariance matrix of the residuals to remove the noise contributions from the largest eigenvalue or vectors owing to ring currents.Finally,we compare the data from the ten Chinese observatories to seven European observatories.Clearly,the covariance matrix method can simulate the SVs of Dst,the jerk of the northward component in 2014 and that of the eastward component in 2003.5 in China are highly agree with that of Vertically downward component in Europe,compare to CHAOS-6,covariance matrix method can show more details of SVs.
基金supported by the National Natural Science Foundation of China(618711496197115962071144)。
文摘Most of the reconstruction-based robust adaptive beamforming(RAB)algorithms require the covariance matrix reconstruction(CMR)by high-complexity integral computation.A Gauss-Legendre quadrature(GLQ)method with the highest algebraic precision in the interpolation-type quadrature is proposed to reduce the complexity.The interference angular sector in RAB is regarded as the GLQ integral range,and the zeros of the threeorder Legendre orthogonal polynomial is selected as the GLQ nodes.Consequently,the CMR can be efficiently obtained by simple summation with respect to the three GLQ nodes without integral.The new method has significantly reduced the complexity as compared to most state-of-the-art reconstruction-based RAB techniques,and it is able to provide the similar performance close to the optimal.These advantages are verified by numerical simulations.
基金Supported by the Aviation Science Fund (No. 20080152004)China Postdoctoral Foundation (No. 20090461119)
文摘An improved two-channel Synthetic Aperture Radar Ground Moving Target Indication (SAR-GMTI) method based on eigen-decomposition of the covariance matrix is investigated. Based on the joint Probability Density Function (PDF) of the Along-Track Interferometric (ATI) phase and the similarity between the two SAR complex images, a novel ellipse detector is presented and is applied to the indication of ground moving targets. We derive its statistics and analyze the performance of detection process in detail. Compared with the approach using the ATI phase, the ellipse detector has a better performance of detection in homogenous clutter. Numerical experiments on simulated data are presented to validate the improved performance of the ellipse detector with respect to the ATI phase approach. Finally, the detection capability of the proposed method is demonstrated by measured SAR data.
基金supported by the Funding Project for Academic Human Resources Development in Institutions of Higher Learning Under the Jurisdiction of Beijing Municipality (0506011200702)National Natural Science Foundation of China+2 种基金Tian Yuan Special Foundation (10926059)Foundation of Zhejiang Educational Committee (Y200803920)Scientific Research Foundation of Hangzhou Dianzi University(KYS025608094)
文摘In this article, the problem of estimating the covariance matrix in general linear mixed models is considered. Two new classes of estimators obtained by shrinking the eigenvalues towards the origin and the arithmetic mean, respectively, are proposed. It is shown that these new estimators dominate the unbiased estimator under the squared error loss function. Finally, some simulation results to compare the performance of the proposed estimators with that of the unbiased estimator are reported. The simulation results indicate that these new shrinkage estimators provide a substantial improvement in risk under most situations.
基金Supported by the NSF of Henan Province(0611052600)
文摘Multivariate seemingly unrelated regression system is raised first and the two stage estimation and its covariance matrix are given. The results of the literatures[1-5] are extended in this paper.
文摘Differential evolution algorithm based on the covariance matrix learning can adjust the coordinate system according to the characteristics of the population, which make<span style="font-family:Verdana;">s</span><span style="font-family:Verdana;"> the search move in a more favorable direction. In order to obtain more accurate information about the function shape, this paper propose</span><span style="font-family:Verdana;">s</span><span style="font-family:;" "=""> <span style="font-family:Verdana;">covariance</span><span style="font-family:Verdana;"> matrix learning differential evolution algorithm based on correlation (denoted as RCLDE)</span></span><span style="font-family:;" "=""> </span><span style="font-family:Verdana;">to improve the search efficiency of the algorithm. First, a hybrid mutation strategy is designed to balance the diversity and convergence of the population;secondly, the covariance learning matrix is constructed by selecting the individual with the less correlation;then, a comprehensive learning mechanism is comprehensively designed by two covariance matrix learning mechanisms based on the principle of probability. Finally,</span><span style="font-family:;" "=""> </span><span style="font-family:;" "=""><span style="font-family:Verdana;">the algorithm is tested on the CEC2005, and the experimental results are compared with other effective differential evolution algorithms. The experimental results show that the algorithm proposed in this paper is </span><span style="font-family:Verdana;">an effective algorithm</span><span style="font-family:Verdana;">.</span></span>
基金supported in part by the National Natural Science Foundation of China(No.62273287)by the Shenzhen Science and Technology Innovation Council(Nos.JCYJ20220530143418040,JCY20170411102101881)the Thousand Youth Talents Plan funded by the central government of China.
文摘Regularized system identification has become the research frontier of system identification in the past decade.One related core subject is to study the convergence properties of various hyper-parameter estimators as the sample size goes to infinity.In this paper,we consider one commonly used hyper-parameter estimator,the empirical Bayes(EB).Its convergence in distribution has been studied,and the explicit expression of the covariance matrix of its limiting distribution has been given.However,what we are truly interested in are factors contained in the covariance matrix of the EB hyper-parameter estimator,and then,the convergence of its covariance matrix to that of its limiting distribution is required.In general,the convergence in distribution of a sequence of random variables does not necessarily guarantee the convergence of its covariance matrix.Thus,the derivation of such convergence is a necessary complement to our theoretical analysis about factors that influence the convergence properties of the EB hyper-parameter estimator.In this paper,we consider the regularized finite impulse response(FIR)model estimation with deterministic inputs,and show that the covariance matrix of the EB hyper-parameter estimator converges to that of its limiting distribution.Moreover,we run numerical simulations to demonstrate the efficacy of ourtheoretical results.
基金supported by the National Natural Science Foundation of China(Nos.62075135 and 61975126)the Science,Technology and Innovation Commission of Shenzhen Municipality(Nos.JCYJ20190808174819083 and JCYJ20190808175201640)。
文摘Multiphoton microscopy is the enabling tool for biomedical research,but the aberrations of biological tissues have limited its imaging performance.Adaptive optics(AO)has been developed to partially overcome aberration to restore imaging performance.For indirect AO,algorithm is the key to its successful implementation.Here,based on the fact that indirect AO has an analogy to the black-box optimization problem,we successfully apply the covariance matrix adaptation evolution strategy(CMA-ES)used in the latter,to indirect AO in multiphoton microscopy(MPM).Compared with the traditional genetic algorithm(GA),our algorithm has a greater improvement in convergence speed and convergence accuracy,which provides the possibility of realizing real-time dynamic aberration correction for deep in vivo biological tissues.
文摘The estimation of covariance matrices is very important in many fields, such as statistics. In real applications, data are frequently influenced by high dimensions and noise. However, most relevant studies are based on complete data. This paper studies the optimal estimation of high-dimensional covariance matrices based on missing and noisy sample under the norm. First, the model with sub-Gaussian additive noise is presented. The generalized sample covariance is then modified to define a hard thresholding estimator , and the minimax upper bound is derived. After that, the minimax lower bound is derived, and it is concluded that the estimator presented in this article is rate-optimal. Finally, numerical simulation analysis is performed. The result shows that for missing samples with sub-Gaussian noise, if the true covariance matrix is sparse, the hard thresholding estimator outperforms the traditional estimate method.
文摘In this paper,a new correlated covariance matrix for Multi-Input Multi-Output(MIMO)radar is proposed,which has lower Side Lobe Levels(SLLs)compared to the new covariance matrix designs and the well-known multi-antenna radar designs including phased-array,MIMO radar and phased-MIMO radar schemes.It is shown that Binary Phased-Shift Keying(BPSK)waveforms that have constant envelope can be used in a closed-form to realize the proposed covariance matrix.Therefore,there is no need to deploy different types of radio amplifiers in the transmitter which will reduce the cost,considerably.The proposed design allows the same transmit power from each antenna in contrast to the phased-MIMO radar.Moreover,the proposed covariance matrix is full-rank and has the same capability as MIMO radar to identify more targets,simultaneously.Performance of the proposed transmit covariance matrix including receive beampattern and output Signal-to-Interference plus Noise Ratio(SINR)is simulated,which validates analytical results.
基金The work was supported in part by the National Natural Science Foundation of China(Nos.11431002,11171018,71271021,11301022).
文摘This paper aims at achieving a simultaneously sparse and low-rank estimator from the semidefinite population covariance matrices.We first benefit from a convex optimization which develops l1-norm penalty to encourage the sparsity and nuclear norm to favor the low-rank property.For the proposed estimator,we then prove that with high probability,the Frobenius norm of the estimation rate can be of order O(√((slgg p)/n))under a mild case,where s and p denote the number of nonzero entries and the dimension of the population covariance,respectively and n notes the sample capacity.Finally,an efficient alternating direction method of multipliers with global convergence is proposed to tackle this problem,and merits of the approach are also illustrated by practicing numerical simulations.
基金supported by the National Basic Research Program of China(Grant No.2011CB921601)the National Natural Science Foundation of China(Grant No.11234008)+1 种基金the NSFC Project for Excellent Research Team(Grant Nos.61121064 and 11234008)Doctoral Program Foundation of the Ministry of Education China(Grant No.20111401130001)
文摘A two-mode entangled state was generated experimentally through mixing two squeezed lights from two optical parametric amplifiers on a 50/50 beam splitter.The entangled beams were measured by means of two pairs of balanced homodyne detection systems respectively.The relative phases between the local beams and the detected beams can be locked by using the optical phase modulation technique.The covariance matrix of the two-mode entangled state was obtained when the relative phase of the local beam and the detected beam in one homodyne detection system is locked and the other is scanned.This method provides a way by which one can extract the covariance matrix of any selected quadrature components of two-mode Gaussian state.
基金supported by the National Natural Science Foundation of China under Grant No.11371236the Graduate Student Innovation Foundation of Shanghai University of Finance and Economics(CXJJ-2015-440)
文摘This paper is concerned with tile proOlenl or improving hue ~lma^u~ u~ under Stein＇s loss. By the partial Iwasawa coordinates of covariance matrix, the corresponding risk can be split into three parts. One can use the information in the weighted matrix of weighted quadratic loss to improve one part of risk. However, this paper indirectly takes advantage of the information in the sample mean and reuses Iwasawa coordinates to improve the rest of risk. It is worth mentioning that the process above can be repeated. Finally, a Monte Carlo simulation study is carried out to verify the theoretical results.
基金This research is supported by National Natural Science Foundation of China, Tian Yuan Special Foundation under Grant No. 10926059 and Zhejiang Provincial Natural Science Foundation of China under Grant No. Y6100053.
文摘In this paper, the problem of estimating the covariance matrix in general linear mixed models is considered. A new class of estimators is proposed. It is shown that this new estimator dominates the analysis of variance estimate under two squared loss functions. Finally, some simulation results to compare the performance of the proposed estimator with that of the analysis of variance estimate are reported. The simulation results indicate that this new estimator provides a substantial improvement in risk under most situations.
基金supported by National Natural Science Foundation of China(Grant Nos.10971122,11101274 and 11322109)Scientific and Technological Projects of Shandong Province(Grant No.2009GG10001012)Excellent Young Scientist Foundation of Shandong Province(Grant No.BS2012SF025)
文摘Covariance matrix plays an important role in risk management, asset pricing, and portfolio allocation. Covariance matrix estimation becomes challenging when the dimensionality is comparable or much larger than the sample size. A widely used approach for reducing dimensionality is based on multi-factor models. Although it has been well studied and quite successful in many applications, the quality of the estimated covariance matrix is often degraded due to a nontrivial amount of missing data in the factor matrix for both technical and cost reasons. Since the factor matrix is only approximately low rank or even has full rank, existing matrix completion algorithms are not applicable. We consider a new matrix completion paradigm using the factor models directly and apply the alternating direction method of multipliers for the recovery. Numerical experiments show that the nuclear-norm matrix completion approaches are not suitable but our proposed models and algorithms are promising.
基金supported by the National Natural Science Foundation of China(51279043,61201411)the Fundamental Research Funds for the Central Universities(HEUCF120502)the National Key Laboratory on Underwater Acoustic Technology Foundation of China(9140C200203110C2001)
文摘Consider the problems of frequency-invariant beampattern optimization and robustness in broadband beamforming.Firstly,a global optimization algorithm,which is based on phase compensation of the array manifolds,is used to construct the frequency-invariant beampattern.Compared with some methods presented recently,the proposed algorithm is not only available to get the global optimal solution,but also simple for physical realization.Meanwhile,a robust adaptive broadband beamforming algorithm is also derived by reconstructing the covariance matrix.The essence of the proposed algorithm is to estimate the space-frequency spectrum using Capon estimator firstly,then integrate over a region separated from the desired signal direction to reconstruct the interference-plus-noise covariance matrix,and finally caleulate the adaptive beamformer weights with the reconstructed matrix.The design of beamformer is formulated as a convex optimization problem to be solved.Simulation results show that the performance of the proposed algorithm is almost always close to the optimal value across a wide range of signal to noise ratios.
基金supported by the National Natural Science Foundation of China(Nos.61374027,11871357)the Sichuan Science and Technology Program(Nos.2019YJ0122)。
文摘This paper addresses the issue of testing sphericity and identity of high-dimensional population covariance matrix when the data dimension exceeds the sample size.The central limit theorem of the first four moments of eigenvalues of sample covariance matrix is derived using random matrix theory for generally distributed populations.Further,some desirable asymptotic properties of the proposed test statistics are provided under the null hypothesis as data dimension and sample size both tend to infinity.Simulations show that the proposed tests have a greater power than existing methods for the spiked covariance model.
