Source author record

Philippe Loubaton

Philippe Loubaton appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

15works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

15 published item(s)

preprint2021arXiv

On the asymptotic behaviour of the eigenvalue distribution of block correlation matrices of high-dimensional time series

We consider linear spectral statistics built from the block-normalized correlation matrix of a set of $M$ mutually independent scalar time series. This matrix is composed of $M \times M$ blocks that contain the sample cross correlation between pairs of time series. In particular, each block has size $L \times L$ and contains the sample cross-correlation measured at $L$ consecutive time lags between each pair of time series. Let $N$ denote the total number of consecutively observed windows that are used to estimate these correlation matrices. We analyze the asymptotic regime where $M,L,N \rightarrow +\infty$ while $ML/N \rightarrow c_\star$, $0<c_\star<\infty$. We study the behavior of linear statistics of the eigenvalues of this block correlation matrix under these asymptotic conditions and show that the empirical eigenvalue distribution converges to a Marcenko-Pastur distribution. Our results are potentially useful in order to address the problem of testing whether a large number of time series are uncorrelated or not.

preprint2015arXiv

A CLT for an improved subspace estimator with observations of increasing dimensions

This paper deals with subspace estimation in the small sample size regime, where the number of samples is comparable in magnitude with the observation dimension. The traditional estimators, mostly based on the sample correlation matrix, are known to perform well as long as the number of available samples is much larger than the observation dimension. However, in the small sample size regime, the performance degrades. Recently, based on random matrix theory results, a new subspace estimator was introduced, which was shown to be consistent in the asymptotic regime where the number of samples and the observation dimension converge to infinity at the same rate. In practice, this estimator outperforms the traditional ones even for certain scenarios where the observation dimension is small and of the same order of magnitude as the number of samples. In this paper, we address a performance analysis of this recent estimator, by proving a central limit theorem in the above asymptotic regime. We propose an accurate approximation of the mean square error, which can be evaluated numerically.

preprint2015arXiv

Large system analysis of a GLRT for detection with large sensor arrays in temporally white noise

This paper addresses the behaviour of a classical multi-antenna GLRT test that allows to detect the presence of a known signal corrupted by a multi-path propagation channel and by an additive white Gaussian noise with unknown spatial covariance matrix. The paper is focused on the case where the number of sensors M is large, and of the same order of magnitude as the sample size N, a context which is modeled by the large system asymptotic regime M goes to infinity, N goes to infinity in such a way that M/N goes to c for c in (0, infinity). The purpose of this paper is to study the behaviour of a GLRT test statistics in this regime, and to show that the corresponding theoretical analysis allows to accurately predict the performance of the test when M and N are of the same order of magnitude.

preprint2015arXiv

On the almost sure location of the singular values of certain Gaussian block-Hankel large random matrices

This paper studies the almost sure location of the eigenvalues of matrices ${\bf W}_N {\bf W}_N^{*}$ where ${\bf W}_N = ({\bf W}_N^{(1)T}, ..., {\bf W}_N^{(M)T})^{T}$ is a $ML \times N$ block-line matrix whose block-lines $({\bf W}_N^{(m)})_{m=1, ..., M}$ are independent identically distributed $L \times N$ Hankel matrices built from i.i.d. standard complex Gaussian sequences. It is shown that if $M \rightarrow +\infty$ and $\frac{ML}{N} \rightarrow c_*$ ($c_* \in (0, \infty)$), then the empirical eigenvalue distribution of ${\bf W}_N {\bf W}_N^{*}$ converges almost surely towards the Marcenko-Pastur distribution. More importantly, it is established that if $L = \mathcal{O}(N^α)$ with $α< 2/3$, then, almost surely, for $N$ large enough, the eigenvalues of ${\bf W}_N {\bf W}_N^{*}$ are located in the neighbourhood of the Marcenko-Pastur distribution.

preprint2015arXiv

Performance analysis of spatial smoothing schemes in the context of large arrays

This paper adresses the statistical behaviour of spatial smoothing subspace DoA estimation schemes using a sensor array in the case where the number of observations $N$ is significantly smaller than the number of sensors $M$, and that the smoothing parameter $L$ is such that $M$ and $NL$ are of the same order of magnitude. This context is modelled by an asymptotic regime in which $NL$ and $M$ both converge towards $\infty$ at the same rate. As in recent works devoted to the study of (unsmoothed) subspace methods in the case where $M$ and $N$ are of the same order of magnitude, it is shown that it is still possible to derive improved DoA estimators termed as Generalized-MUSIC with spatial smoothing (G-MUSIC SS). The key ingredient of this work is a technical result showing that the largest singular values and corresponding singular vectors of low rank deterministic perturbation of certain Gaussian block-Hankel large random matrices behave as if the entries of the latter random matrices were independent identically distributed. This allows to conclude that when the number of sources and their DoA do not scale with $M,N,L,$ a situation modelling widely spaced DoA scenarios, then both traditional and Generalized spatial smoothing subspace methods provide consistent DoA estimators whose convergence speed is faster than $\frac{1}{M}$. The case of DoA that are spaced of the order of a beamwidth, which models closely spaced sources, is also considered. It is shown that the convergence speed of G-MUSIC SS estimates is unchanged, but that it is no longer the case for MUSIC SS ones.

preprint2012arXiv

A Subspace Estimator for Fixed Rank Perturbations of Large Random Matrices

This paper deals with the problem of parameter estimation based on certain eigenspaces of the empirical covariance matrix of an observed multidimensional time series, in the case where the time series dimension and the observation window grow to infinity at the same pace. In the area of large random matrix theory, recent contributions studied the behavior of the extreme eigenvalues of a random matrix and their associated eigenspaces when this matrix is subject to a fixed-rank perturbation. The present work is concerned with the situation where the parameters to be estimated determine the eigenspace structure of a certain fixed-rank perturbation of the empirical covariance matrix. An estimation algorithm in the spirit of the well-known MUSIC algorithm for parameter estimation is developed. It relies on an approach recently developed by Benaych-Georges and Nadakuditi, relating the eigenspaces of extreme eigenvalues of the empirical covariance matrix with eigenspaces of the perturbation matrix. First and second order analyses of the new algorithm are performed.

preprint2011arXiv

Almost sure localization of the eigenvalues in a gaussian information plus noise model. Applications to the spiked models

Let $\boldsymbolΣ_N$ be a $M \times N$ random matrix defined by $\boldsymbolΣ_N = \mathbf{B}_N + σ\mathbf{W}_N$ where $\mathbf{B}_N$ is a uniformly bounded deterministic matrix and where $\mathbf{W}_N$ is an independent identically distributed complex Gaussian matrix with zero mean and variance $\frac{1}{N}$ entries. The purpose of this paper is to study the almost sure location of the eigenvalues $\hatλ_{1,N} \geq ... \geq \hatλ_{M,N}$ of the Gram matrix ${\boldsymbol Σ}_N {\boldsymbol Σ}_N^*$ when $M$ and $N$ converge to $+\infty$ such that the ratio $c_N = \frac{M}{N}$ converges towards a constant $c > 0$. The results are used in order to derive, using an alernative approach, known results concerning the behaviour of the largest eigenvalues of ${\boldsymbol Σ}_N {\boldsymbol Σ}_N^*$ when the rank of $\mathbf{B}_N$ remains fixed when $M$ and $N$ converge to $+\infty$.

preprint2011arXiv

Diversity of the MMSE receiver in flat fading and frequency selective MIMO channels at fixed rate

In this contribution, the evaluation of the diversity of the MIMO MMSE receiver is addressed for finite rates in both flat fading channels and frequency selective fading channels with cyclic prefix. It has been observed recently that in contrast with the other MIMO receivers, the MMSE receiver has a diversity depending on the aimed finite rate, and that for sufficiently low rates the MMSE receiver reaches the full diversity - that is, the diversity of the ML receiver. This behavior has so far only been partially explained. The purpose of this paper is to provide complete proofs for flat fading MIMO channels, and to improve the partial existing results in frequency selective MIMO channels with cyclic prefix.

preprint2011arXiv

Improved subspace estimation for multivariate observations of high dimension: the deterministic signals case

We consider the problem of subspace estimation in situations where the number of available snapshots and the observation dimension are comparable in magnitude. In this context, traditional subspace methods tend to fail because the eigenvectors of the sample correlation matrix are heavily biased with respect to the true ones. It has recently been suggested that this situation (where the sample size is small compared to the observation dimension) can be very accurately modeled by considering the asymptotic regime where the observation dimension $M$ and the number of snapshots $N$ converge to $+\infty$ at the same rate. Using large random matrix theory results, it can be shown that traditional subspace estimates are not consistent in this asymptotic regime. Furthermore, new consistent subspace estimate can be proposed, which outperform the standard subspace methods for realistic values of $M$ and $N$. The work carried out so far in this area has always been based on the assumption that the observations are random, independent and identically distributed in the time domain. The goal of this paper is to propose new consistent subspace estimators for the case where the source signals are modelled as unknown deterministic signals. In practice, this allows to use the proposed approach regardless of the statistical properties of the source signals. In order to construct the proposed estimators, new technical results concerning the almost sure location of the eigenvalues of sample covariance matrices of Information plus Noise complex Gaussian models are established. These results are believed to be of independent interest.

preprint2011arXiv

Large information plus noise random matrix models and consistent subspace estimation in large sensor networks

In array processing, a common problem is to estimate the angles of arrival of $K$ deterministic sources impinging on an array of $M$ antennas, from $N$ observations of the source signal, corrupted by gaussian noise. The problem reduces to estimate a quadratic form (called "localization function") of a certain projection matrix related to the source signal empirical covariance matrix. Recently, a new subspace estimation method (called "G-MUSIC") has been proposed, in the context where the number of available samples $N$ is of the same order of magnitude than the number of sensors $M$. In this context, the traditional subspace methods tend to fail because the empirical covariance matrix of the observations is a poor estimate of the source signal covariance matrix. The G-MUSIC method is based on a new consistent estimator of the localization function in the regime where $M$ and $N$ tend to $+\infty$ at the same rate. However, the consistency of the angles estimator was not adressed. The purpose of this paper is to prove the consistency of the angles of arrival estimator in the previous asymptotic regime. To prove this result, we show the property that the singular values of M x N Gaussian information plus noise matrix escape from certain intervals is an event of probability decreasing at rate O(1/N^p) for all p. A regularization trick is also introduced, which allows to confine these singular values into certain intervals and to use standard tools as Poincaré inequality to characterize any moments of the estimator. These results are believed to be of independent interest.

preprint2011arXiv

On bilinear forms based on the resolvent of large random matrices

Consider a matrix $Σ_n$ with random independent entries, each non-centered with a separable variance profile. In this article, we study the limiting behavior of the random bilinear form $u_n^* Q_n(z) v_n$, where $u_n$ and $v_n$ are deterministic vectors, and Q_n(z) is the resolvent associated to $Σ_n Σ_n^*$ as the dimensions of matrix $Σ_n$ go to infinity at the same pace. Such quantities arise in the study of functionals of $Σ_n Σ_n^*$ which do not only depend on the eigenvalues of $Σ_n Σ_n^*$, and are pivotal in the study of problems related to non-centered Gram matrices such as central limit theorems, individual entries of the resolvent, and eigenvalue separation.

preprint2011arXiv

On the Capacity Achieving Covariance Matrix for Frequency Selective MIMO Channels Using the Asymptotic Approach

In this contribution, an algorithm for evaluating the capacity-achieving input covariance matrices for frequency selective Rayleigh MIMO channels is proposed. In contrast with the flat fading Rayleigh case, no closed-form expressions for the eigenvectors of the optimum input covariance matrix are available. Classically, both the eigenvectors and eigenvalues are computed numerically and the corresponding optimization algorithms remain computationally very demanding. In this paper, it is proposed to optimize (w.r.t. the input covariance matrix) a large system approximation of the average mutual information derived by Moustakas and Simon. The validity of this asymptotic approximation is clarified thanks to Gaussian large random matrices methods. It is shown that the approximation is a strictly concave function of the input covariance matrix and that the average mutual information evaluated at the argmax of the approximation is equal to the capacity of the channel up to a O(1/t) term, where t is the number of transmit antennas. An algorithm based on an iterative waterfilling scheme is proposed to maximize the average mutual information approximation, and its convergence studied. Numerical simulation results show that, even for a moderate number of transmit and receive antennas, the new approach provides the same results as direct maximization approaches of the average mutual information.

preprint2010arXiv

On the Capacity Achieving Covariance Matrix for Frequency Selective MIMO Channels Using the Asymptotic Approach

In this contribution, an algorithm for evaluating the capacity-achieving input covariance matrices for frequency selective Rayleigh MIMO channels is proposed. In contrast with the flat fading Rayleigh cases, no closed-form expressions for the eigenvectors of the optimum input covariance matrix are available. Classically, both the eigenvectors and eigenvalues are computed numerically and the corresponding optimization algorithms remain computationally very demanding. In this paper, it is proposed to optimize (w.r.t. the input covariance matrix) a large system approximation of the average mutual information derived by Moustakas and Simon. An algorithm based on an iterative water filling scheme is proposed, and its convergence is studied. Numerical simulation results show that, even for a moderate number of transmit and receive antennas, the new approach provides the same results as direct maximization approaches of the average mutual information.

preprint2010arXiv

On the capacity achieving covariance matrix for Rician MIMO channels: an asymptotic approach

The capacity-achieving input covariance matrices for coherent block-fading correlated MIMO Rician channels are determined. In this case, no closed-form expressions for the eigenvectors of the optimum input covariance matrix are available. An approximation of the average mutual information is evaluated in this paper in the asymptotic regime where the number of transmit and receive antennas converge to $+\infty$. New results related to the accuracy of the corresponding large system approximation are provided. An attractive optimization algorithm of this approximation is proposed and we establish that it yields an effective way to compute the capacity achieving covariance matrix for the average mutual information. Finally, numerical simulation results show that, even for a moderate number of transmit and receive antennas, the new approach provides the same results as direct maximization approaches of the average mutual information, while being much more computationally attractive.

preprint2009arXiv

On the precoder design of flat fading MIMO systems equipped with MMSE receivers: a large system approach

This paper is devoted to the design of precoders maximizing the ergodic mutual information (EMI) of bi-correlated flat fading MIMO systems equiped with MMSE receivers. The channel state information and the second order statistics of the channel are assumed available at the receiver side and at the transmitter side respectively. As the direct maximization of the EMI needs the use of non attractive algorithms, it is proposed to optimize an approximation of the EMI, introduced recently, obtained when the number of transmit and receive antennas $t$ and $r$ converge to $\infty$ at the same rate. It is established that the relative error between the actual EMI and its approximation is a $O(\frac{1}{t^{2}})$ term. It is shown that the left singular eigenvectors of the optimum precoder coincide with the eigenvectors of the transmit covariance matrix, and its singular values are solution of a certain maximization problem. Numerical experiments show that the mutual information provided by this precoder is close from what is obtained by maximizing the true EMI, but that the algorithm maximizing the approximation is much less computationally intensive.