Skip to main content
QUICK REVIEW

[论文解读] On the Power of Truncated SVD for General High-rank Matrix Estimation Problems

Simon S. Du, Yining Wang|arXiv (Cornell University)|Feb 22, 2017
Sparse and Compressive Sensing Techniques参考文献 20被引用 10
一句话总结

本文证明,对一个与高阶正定矩阵 $\mathbf{A}$ 的谱范数接近的估计 $\widehat{\mathbf{A}}$ 进行截断奇异值分解(SVD),可得到 Frobenius 范数的乘法近似,从而在高阶矩阵补全、去噪和协方差估计中实现最优样本复杂度,且无需依赖低秩结构或谱间隙假设。

ABSTRACT

We show that given an estimate $\widehat{A}$ that is close to a general high-rank positive semi-definite (PSD) matrix $A$ in spectral norm (i.e., $\|\widehat{A}-A\|_2 \leq δ$), the simple truncated SVD of $\widehat{A}$ produces a multiplicative approximation of $A$ in Frobenius norm. This observation leads to many interesting results on general high-rank matrix estimation problems, which we briefly summarize below ($A$ is an $n imes n$ high-rank PSD matrix and $A_k$ is the best rank-$k$ approximation of $A$): (1) High-rank matrix completion: By observing $Ω(\frac{n\max\{ε^{-4},k^2\}μ_0^2\|A\|_F^2\log n}{σ_{k+1}(A)^2})$ elements of $A$ where $σ_{k+1}\left(A ight)$ is the $\left(k+1 ight)$-th singular value of $A$ and $μ_0$ is the incoherence, the truncated SVD on a zero-filled matrix satisfies $\|\widehat{A}_k-A\|_F \leq (1+O(ε))\|A-A_k\|_F$ with high probability. (2)High-rank matrix de-noising: Let $\widehat{A}=A+E$ where $E$ is a Gaussian random noise matrix with zero mean and $ν^2/n$ variance on each entry. Then the truncated SVD of $\widehat{A}$ satisfies $\|\widehat{A}_k-A\|_F \leq (1+O(\sqrt{ν/σ_{k+1}(A)}))\|A-A_k\|_F + O(\sqrt{k}ν)$. (3) Low-rank Estimation of high-dimensional covariance: Given $N$ i.i.d.~samples $X_1,\cdots,X_N\sim\mathcal N_n(0,A)$, can we estimate $A$ with a relative-error Frobenius norm bound? We show that if $N = Ω\left(n\max\{ε^{-4},k^2\}γ_k(A)^2\log N ight)$ for $γ_k(A)=σ_1(A)/σ_{k+1}(A)$, then $\|\widehat{A}_k-A\|_F \leq (1+O(ε))\|A-A_k\|_F$ with high probability, where $\widehat{A}=\frac{1}{N}\sum_{i=1}^N{X_iX_i^ op}$ is the sample covariance.

研究动机与目标

  • 为解决谱范数近似在高阶矩阵估计中的局限性,即由于秩放大导致 Frobenius 范数误差可能较大的问题。
  • 证明即使 $\mathbf{A}$ 为高阶矩阵,对粗略的谱范数估计 $\widehat{\mathbf{A}}$ 进行截断 SVD 也能得到 $\mathbf{A}$ 的 Frobenius 范数乘法近似。
  • 建立高阶矩阵补全与去噪的样本复杂度界,使其与 $\mathbf{A}$ 的谱间隙无关,优于以往结果。
  • 通过证明 $N \approx n$ 个独立同分布样本即可实现 Frobenius 范数下的相对误差,从而改进低维协方差估计,优于经典所需的 $N \approx n^2$。

提出的方法

  • 对满足 $\|\widehat{\mathbf{A}} - \mathbf{A}\|_2 \leq \delta$ 的矩阵 $\widehat{\mathbf{A}}$ 应用截断 SVD,得到低秩近似 $\widehat{\mathbf{A}}_k$。
  • 证明 $\|\widehat{\mathbf{A}}_k - \mathbf{A}\|_F \leq C(k,\delta,\sigma_{k+1}(\mathbf{A})) \|\mathbf{A} - \mathbf{A}_k\|_F$,即在 Frobenius 范数下实现乘法近似。
  • 利用 Weyl 不等式与谱分解控制 $\widehat{\mathbf{A}}_k$ 与 $\mathbf{A}$ 之间的误差,利用 $\widehat{\mathbf{A}}$ 与 $\mathbf{A}$ 的谱范数接近性。
  • 分析稀疏性条件 $n\|\mathbf{A}\|_{\max} \leq \mu_0\|\mathbf{A}\|_F$,以确保矩阵补全中的样本复杂度界。
  • 建立多项式衰减谱下的极小极大率,推测 $R_1 \asymp (n/N)^{(2\beta-1)/(2\beta)}$ 与 $R_2 \asymp \nu^{2(2\beta-1)/(2\beta)}$,其中 $\beta > 1/2$。
  • 使用自助法或参数化特征值衰减模型,在 $\sigma_{k+1}(\mathbf{A})$ 未知时估计最优截断秩 $k$。

实验结果

研究问题

  • RQ1对谱范数接近的估计 $\widehat{\mathbf{A}}$ 进行截断 SVD,能否对一般高阶正定矩阵 $\mathbf{A}$ 实现 Frobenius 范数的乘法近似?
  • RQ2高阶矩阵补全是否能在样本复杂度与谱间隙无关的前提下,实现 Frobenius 范数下的 $(1+\varepsilon)$ 相对误差?
  • RQ3在不假设真实矩阵为精确低秩的前提下,能否仅利用噪声观测的谱范数接近性完成矩阵去噪?
  • RQ4低维协方差估计是否可在 $N \approx n$ 个样本下实现 Frobenius 范数下的相对误差,优于经典的 $N \approx n^2$ 要求?
  • RQ5在多项式谱衰减 $\sigma_j \propto j^{-\beta}$ 条件下,矩阵补全与去噪的极小极大率为何?

主要发现

  • 对满足 $\|\widehat{\mathbf{A}} - \mathbf{A}\|_2 \leq \delta$ 的谱范数估计 $\widehat{\mathbf{A}}$ 进行截断 SVD,可得到乘法 Frobenius 近似:$\|\widehat{\mathbf{A}}_k - \mathbf{A}\|_F \leq C(k,\delta,\sigma_{k+1}(\mathbf{A})) \|\mathbf{A} - \mathbf{A}_k\|_F$。
  • 在 $\mu_0$-稀疏性条件下,高阶矩阵补全的截断 SVD 估计器可实现 Frobenius 范数下 $(1+O(\varepsilon))$ 的相对误差,且样本复杂度与谱间隙无关。
  • 在高阶矩阵去噪中,该方法可在不假设 $\mathbf{A}$ 为精确低秩的前提下,恢复出 $\mathbf{A}$,且 Frobenius 范数误差可控。
  • 对于低维协方差估计,$N \approx n$ 个独立同分布样本即可实现对总体协方差 $\mathbf{A}$ 的 Frobenius 范数相对误差近似,优于经典要求的 $N \approx n^2$。
  • 本文推测,在多项式谱衰减 $\sigma_j \propto j^{-\beta}$ 条件下,矩阵补全与去噪的极小极大率分别为 $R_1 \asymp (n/N)^{(2\beta-1)/(2\beta)}$ 与 $R_2 \asymp \nu^{2(2\beta-1)/(2\beta)}$,其中 $\beta > 1/2$。
  • 该方法避免了核范数最小化,且样本复杂度优于以往工作,包括 Achlioptas 和 McSherry(2007)以及 Hardt 和 Wootters(2014)的结果。

更好的研究,从现在开始

从阅读论文到最终审阅,大幅缩短您的研究时间。

无需绑定信用卡

本解读由 AI 生成,并经人工编辑审核。