Graph Sampling for Matrix Completion Using Recurrent Gershgorin Disc Shift

Fen Wang,Yongchao Wang,Gene Cheung,Cheng Yang
DOI: https://doi.org/10.1109/tsp.2020.2988784
IF: 4.875
2020-01-01
IEEE Transactions on Signal Processing
Abstract:Matrix completion algorithms fill missing entries in a large matrix given a subset of observed samples. The problem of how to pre-select a subset of entries for sampling to maximize the reconstructed matrix fidelity is largely unaddressed. In this paper, we propose two sampling algorithms to tackle this problem: (i) a fast base sampling algorithm on general single graphs, and (ii) an extended sampling algorithm from the base algorithm for active matrix completion. Assuming the double graph smoothness prior where rows / columns of the matrix signal are smooth with respect to row / column factor graphs, we first show that a matrix signal can be recovered from partial observations by solving a system of linear equations ${\mathbf Q}{\mathbf x}={\mathbf y}$ , where the sought matrix is vectorized as ${\mathbf x}$ and interpreted as a signal residing on a single graph. On this single graph, to promote small reconstruction error and stability of the linear system, we maximize the smallest eigenvalue $\lambda _{\min }$ of ${\mathbf Q}$ by greedily selecting samples corresponding to the largest magnitude entries in the first eigenvectors of ${\mathbf Q}$ , based on a corollary of the Gershgorin circle theorem. Extending to the active matrix completion case, we alternately choose entries in the rows and columns corresponding to the largest magnitude entries in the first eigenvectors of the row / column coefficient matrices ${\mathbf Q}_r$ and ${\mathbf Q}_c$ . Our algorithm benefits computationally from warm start as the first eigenvectors are computed recurrently for more samples. Extensive experiments on both synthetic and real-world datasets show that our extended graph sampling algorithm outperforms existing sampling schemes for matrix completion, when combined with a range of modern matrix completion algorithms.
What problem does this paper attempt to address?