Perturbation of linear forms of singular vectors under Gaussian noise
Abstract
Let A∈Rm× n be a matrix of rank r with singular value decomposition (SVD) A=Σk=1rσk (uk vk), where \σk, k=1,…,r\ are singular values of A (arranged in a non-increasing order) and uk∈ Rm, vk∈ Rn, k=1,…, r are the corresponding left and right orthonormal singular vectors. Let A=A+X be a noisy observation of A, where X∈Rm× n is a random matrix with i.i.d. Gaussian entries, Xij(0,τ2), and consider its SVD A=Σk=1m nσk(ukvk) with singular values σ1≥…≥σm n and singular vectors uk,vk,k=1,…, m n. The goal of this paper is to develop sharp concentration bounds for linear forms uk,x, x∈ Rm and vk,y, y∈ Rn of the perturbed (empirical) singular vectors in the case when the singular values of A are distinct and, more generally, concentration bounds for bilinear forms of projection operators associated with SVD. In particular, the results imply upper bounds of the order O((m+n)m n) (holding with a high probability) on 1≤ i≤ m|<uk-1+bkuk,eim>|\ \ and \ \ 1≤ j≤ n|<vk-1+bkvk,ejn>|, where bk are properly chosen constants characterizing the bias of empirical singular vectors uk, vk and \eim,i=1,…,m\, \ejn,j=1,…,n\ are the canonical bases of Rm, Rn, respectively.
Turn this paper into a lesson
ArcXiv compiles a structured reading guide from this paper's metadata: plain-English importance, contributions, prerequisite concepts, which sections to read first, flashcards, and a quiz. Grounded in the abstract, never invented.