0
$\begingroup$

I have a question regarding singular value decomposition of a matrix which I have not been able to answer. Why the matrices in the svd decomposition of a matrix, are one rank matrices?

$\endgroup$
  • $\begingroup$ They are not necessarily so. An SVD of $I$ is $I \cdot I \cdot I$, none of which are rank one (for $n >1$). $\endgroup$ – copper.hat Jun 11 at 20:35
  • $\begingroup$ A rank $r$ matrix can be written as a sum of $r$ rank $1$ matrices. $\endgroup$ – Upstart Jun 11 at 21:02
1
$\begingroup$

Further to @copper.hat comment, the Matrices of the SVD are not necessarily rank 1. In fact, matrices could be anywhere between 1 and $n$. Recall the definition of Singular Value Decomposition,

A matrix $A \in \mathbb{C}^{m \times n}$ has a singular value decomposition $A=U\Sigma V^*$, where $U\in \mathbb{C}^{m \times m}, V \in \mathbb{C}^{n \times n}$ are unitary and $\Sigma = \text{diag}(\sigma_1, \dots , \sigma_p)\in \mathbb{R}^{m \times n}, p = \min (m,n),$ where $\sigma_1 \ge \sigma_2 \ge \dots \sigma_p \ge 0$.

Since the matrix $\Sigma$ is defined by its leading diagonal, the rank of $A$ (and $\Sigma$) is equal to the number of non-zero $\sigma_i$ values ($1\le i \le p$).

So the maximum special case is when $A$ is a square matrix of $\text{rank}(A)=n$. In this case we actually have the special case of Schur's theorem. Where $A$ is a normal matrix ($A^* A=AA^* $) too.

Schur's Theorem is:

Let $A\in \mathbb{C}^{n \times n}$. Then there exists a unitary matrix $U$ and an upper triangular matrix $T$ such that $T = U^* A U$.

The fact that $A$ is normal in this case means that it is diagonalisable and $T$ is diagonal with values equal to the eigenvalues of $A$.

In summary, the SVD is a type of a similarity transformation on a Rectangular matrix. The SVD is possibly the most important decomposition of a matrix as it shows all characteristics of a matrix, its singular values and left and right eigenvectors. If you are interested in this further, I'd recommend Gilbert Strang's description of the SVD and it's power in the new MIT course Learning from Data - lectures available here: https://www.youtube.com/watch?v=rYz83XPxiZo

$\endgroup$
1
$\begingroup$

Why the matrices in the svd decomposition of a matrix, are one rank matrices?

I think you mean to ask why the SVD is a sum of rank-one matrices. That is why if $A = U \Sigma V^{T}$ then we have that

$$ A = \sum_{i=1}^{n} \sigma_{i} u_{i} v_{i}^{T}$$

You can see $u_{i} v_{i}^{T}$ as the matrix generated by the outer product of the vectors $u_{i}, v_{i}$ which are the left and right singular vectors respectively. Then $\sigma_{i}$ acts as a scaling factor along these new axes.

$\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.