33
$\begingroup$

A method called "Robust PCA" solves the matrix decomposition problem

$$L^*, S^* = \arg \min_{L, S} \|L\|_* + \|S\|_1 \quad \text{s.t. } L + S = X$$

as a surrogate for the actual problem

$$L^*, S^* = \arg \min_{L, S} rank(L) + \|S\|_0 \quad \text{s.t. } L + S = X,$$ i.e. the actual goal is to decompose the data matrix $X$ into a low-rank signal matrix $L$ and a sparse noise matrix $S$. In this context: why is the nuclear norm a good approximation for the rank of a matrix? I can think of matrices with low nuclear norm but high rank and vice-versa. Is there any intuition one can appeal to?

$\endgroup$
3
  • 5
    $\begingroup$ Consider this question: Why is the $\ell_1$ norm a surrogate for $\ell_0$? $\endgroup$ – p.s. Jun 28 '12 at 16:33
  • $\begingroup$ You write "low-rank" but the formula has $\arg \max$... $\endgroup$ – user31373 Jun 28 '12 at 23:01
  • $\begingroup$ @LeonidKovalev: Oops. Of course I meant $\min$. $\endgroup$ – blubb Jun 29 '12 at 7:26
21
$\begingroup$

Why does compressed sensing work? Because the $\ell_1$ ball in high dimensions is extremely "pointy" -- the extreme values of a linear function on this ball are very likely to be attained on the faces of low dimensions, those that consist of sparse vectors. When applied to matrices, the sparseness of the set of eigenvalues means low rank, as @mrig wrote before me.

$\endgroup$
2
  • 2
    $\begingroup$ Thank you, the reduction to sparse coding of eigenvalues was indeed the link I was looking for! $\endgroup$ – blubb Jun 29 '12 at 7:39
  • $\begingroup$ I think for linear objective functions, there is a proofs saying optimum happens on vertices. But what if the objective is nonlinear? Does minimizing a convex nonlinear loss of a param + nuclear norm of the param has anything beyond pure wish to be low rank? I think the answer is yes, but can not see the reason. $\endgroup$ – user25004 Aug 11 '16 at 19:36
12
$\begingroup$

The nuclear norm can be thought of as a convex relaxation of the number of non-zero eigenvalues (i.e. the rank).

$\endgroup$
9
$\begingroup$

A nuclear norm of a matrix is equivalent to the L1-norm of the vector of its eigenvalues. Thus, you are injecting sparsity to the vector of eigenvalues. Essentially, this sparsity means you are reducing the rank of the original matrix.

$\endgroup$
2
  • $\begingroup$ Why does a sparse eigenvector imply low rank? $\endgroup$ – the_dude Jan 14 '18 at 6:37
  • 3
    $\begingroup$ He means the vector of eigenvalues not an eigenvector. Fewer non-zero eigenvalues corresponds to lower rank. $\endgroup$ – Srinivas Eswar May 8 '18 at 19:39
8
$\begingroup$

To be accurate, it has been shown that the $\ell_1$ norm is the convex envelope of the $\| \cdot \|_0$ pseudo-norm while the nuclear norm is the convex envelope of the rank.

As a reminder, the convex envelope is the tightest convex surrogate of a function. An important property is that a function and its convex envelope have the same global minimizer.

$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.