Toeplitz Low-Rank Approximation with Sublinear Query Complexity

Abstract

We present a sublinear query algorithm for outputting a near-optimal low-rank approximation to any positive semidefinite Toeplitz matrix T ∈ Rd × d. In particular, for any integer rank k ≤ d and ε,δ > 0, our algorithm makes O (k2 · (1/δ) · poly(1/ε) ) queries to the entries of T and outputs a rank O (k · (1/δ)/ε ) matrix T ∈ Rd × d such that \| T - T\|F ≤ (1+ε) · \|T-Tk\|F + δ \|T\|F. Here, \|·\|F is the Frobenius norm and Tk is the optimal rank-k approximation to T, given by projection onto its top k eigenvectors. O(·) hides polylog(d) factors. Our algorithm is structure-preserving, in that the approximation T is also Toeplitz. A key technical contribution is a proof that any positive semidefinite Toeplitz matrix in fact has a near-optimal low-rank approximation which is itself Toeplitz. Surprisingly, this basic existence result was not previously known. Building on this result, along with the well-established off-grid Fourier structure of Toeplitz matrices [Cybenko'82], we show that Toeplitz T with near optimal error can be recovered with a small number of random queries via a leverage-score-based off-grid sparse Fourier sampling scheme.

0

Turn this paper into a full lesson

ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…