Consider a PSD matrix C=RRT. Let R=USVT be an SVD of R, we have C=UΛUT where Λ=SST=diag(σ12,…,σr2,0,…,0)
Hence, the eigenvalues of C are non-negative.
Conversely, if C=U∧UT with Λ=diag(λ1,…,λn),λi≥0 for every i, then C=RRT, with R=Udiag(λ1,…,λn)
Congruence transformation:
For any matrix A∈Rm,n it holds that:
ATA⪰0, and AAT⪰0;
ATA≻0 if and only if A is full-column rank, i.e., rankA=n;
AAT≻0 if and only if A is full-row rank, i.e., rankA=m
Matrix square-root
Let A∈Sn. Then A⪰0⇔∃B⪰0:A=B2A≻0⇔∃B≻0:A=B2.
Matrix B=A1/2 is called the matrix square-root of A.
Any A⪰0 admits the spectral factorization A=U∧UT, with U orthogonal and Λ=diag(λ1,…,λn),λi≥0,i=1,…,n.
Defining Λ1/2=diag(λ1,…,λ1) and B=UΛ1/2UT : A⪰0⇔∃B:A=BTBA≻0⇔∃B nonsingular : A=BTB
A is positive definite if and only if it is congruent to the identity
As Ellipsoids:
Positive-definite matrices are intimately related to geometrical objects called ellipsoids.
A full-dimensional, bounded ellipsoid with center in the origin can indeed be defined as the set E={x∈Rn:xTP−1x≤1},P≻0.
The eigenvalues λi and eigenvectors ui of P define the orientation and shape of the ellipsoid: ui the directions of the semi axes of the ellipsoid, while their lengths are given by λi.
Using the Cholesky DecompositionP−1=ATA, the previous definition of ellipsoid E is also equivalent to: E={x∈Rn:∥Ax∥2≤1}.