🤖 AI Summary
This paper investigates the fundamental statistical limits of Laplacian matrix estimation under general observation models. Addressing its intrinsic structural constraints—symmetry, zero-sum nullspace (i.e., $mathbf{1} in ker(L)$), and sparsity—we establish, for the first time, a structure-aware Cramér–Rao bound (CRB) framework. We propose two oracle CRBs incorporating support-set priors and extend them to settings with sparse support priors. Methodologically, we integrate linear reparameterization, constrained maximum likelihood estimation (CMLE), the Slepian–Bangs formula, and graph signal processing tools. The derived CRB provides a unified performance benchmark for Gaussian Markov random field modeling, graph filter design, and power system topology identification. Tightness is empirically validated in graphical LASSO and power grid topology recovery: both CMLE and oracle CMLE achieve asymptotic mean-squared error convergence to the bound as sample size increases.
📝 Abstract
In this paper, we analyze the performance of the estimation of Laplacian matrices under general observation models. Laplacian matrix estimation involves structural constraints, including symmetry and null-space properties, along with matrix sparsity. By exploiting a linear reparametrization that enforces the structural constraints, we derive closed-form matrix expressions for the Cramer-Rao Bound (CRB) specifically tailored to Laplacian matrix estimation. We further extend the derivation to the sparsity-constrained case, introducing two oracle CRBs that incorporate prior information of the support set, i.e. the locations of the nonzero entries in the Laplacian matrix. We examine the properties and order relations between the bounds, and provide the associated Slepian-Bangs formula for the Gaussian case. We demonstrate the use of the new CRBs in three representative applications: (i) topology identification in power systems, (ii) graph filter identification in diffused models, and (iii) precision matrix estimation in Gaussian Markov random fields under Laplacian constraints. The CRBs are evaluated and compared with the mean-squared-errors (MSEs) of the constrained maximum likelihood estimator (CMLE), which integrates both equality and inequality constraints along with sparsity constraints, and of the oracle CMLE, which knows the locations of the nonzero entries of the Laplacian matrix. We perform this analysis for the applications of power system topology identification and graphical LASSO, and demonstrate that the MSEs of the estimators converge to the CRB and oracle CRB, given a sufficient number of measurements.