特征值+SVD分解+伪逆(广义逆)

  • Post author:
  • Post category:其他

方阵的特征值与矩阵的相似

A

n

,

n

A

x

=

λ

x

使

d

e

t

(

λ

I

A

)

=

0

λ

A_{n,n}有分解Ax=\lambda x\\ 使用det(\lambda I-A)=0求 \lambda

An,nAx=λx使det(λIA)=0λ

(

)

(

)

d

e

t

几何重数(特征值的特征子空间的维数)\leq 代数重数(根的重数)\\ 特征值和为迹,积为det(利用多项式根与系数的关系)

()()det

A

,

B

:

B

=

M

1

A

M

A

,

B

线

A

B

B

A

A

B

x

=

λ

x

,

B

A

B

x

=

λ

B

x

A,B矩阵相似:B=M^{-1}AM\\ A,B有相同的特征方程与特征值,同一个线性变换在不同基下的矩阵相似\\ AB与BA有相同的非零特征值若ABx=\lambda x,BABx=\lambda B x

A,B:B=M1AMA,B线ABBAABx=λx,BABx=λBx

A

X

=

A

[

x

1

,

x

2

,

x

3

]

=

X

Λ

A

=

X

Λ

X

1

使

Λ

,

Z

(

λ

)

AX=A[x_1,x_2,x_3]=X\Lambda \\ A=X\Lambda X^{-1}\\ 但是即使固定\Lambda ,Z也不唯一(如果一个特征值\lambda 对应的特征子空间不是一维的)

AX=A[x1,x2,x3]=XΛA=XΛX1使Λ,Z(λ)

E

V

D

[

F

k

+

2

F

k

+

1

]

=

[

1

1

1

]

[

F

k

+

1

F

k

]

n

EVD分解斐波那契数列通项公式: \begin{bmatrix}F_{k+2}\\F_{k+1}\end{bmatrix}=\begin{bmatrix}1&1\\&1\end{bmatrix}\begin{bmatrix}F_{k+1}\\F_{k}\end{bmatrix}\\本质是求矩阵的n次方

EVD[Fk+2Fk+1]=[111][Fk+1Fk]n

谱定理

A

n

,

A

Q

=

Q

Λ

Q

A

=

Q

Λ

Q

t

A

A特征值实数有n个正交的特征向量,AQ=Q\Lambda,Q是可逆矩阵\\ A=Q\Lambda Q^t\\ A是实对称方阵,满足以上条件

An,AQ=QΛQA=QΛQtA

SVD分解(链接为实例)

A

A

m

,

n

=

U

Σ

V

t

U

m

,

m

     

Σ

m

,

n

    

V

n

,

n

A

t

A

n

,

n

(

A

t

A

)

v

i

=

λ

i

v

i

A

t

A

n

,

n

U

A

A

m

,

m

t

V

,

(

A

A

t

)

u

i

=

λ

i

u

i

(

A

t

A

A

A

t

)

A

m

,

n

=

U

Σ

V

t

A

V

=

U

Σ

A

v

i

=

σ

i

u

i

σ

i

=

A

v

i

u

i

A

=

U

Σ

V

t

A

t

=

V

Σ

U

t

,

A

t

A

=

V

Σ

2

V

t

A

t

A

S

V

D

V

σ

i

=

λ

i

实数矩阵A不为方阵时,设:A_{m,n}=U\Sigma V^t\\ 其中U_{m,m}为酉矩阵 \ \ \ \ \ \Sigma _{m,n}为主对角非零矩阵 \ \ \ \ V_{n,n}为酉矩阵\\ ①A^tA_{n,n}为方阵,可求特征值(A^tA)v_i=\lambda_i v_i\\ 将A^tA_{n,n}的特征向量张成的矩阵记为U\\ ②将AA^t_{m,m}张成的记为V,(AA^t)u_i=\lambda_i u_i\tiny ( \color{red}A^tA和AA^t有相同的非零特征值\color{black})\normalsize\\\\ ③A_{m,n}=U\Sigma V^t \Rightarrow AV=U\Sigma \Rightarrow Av_i=\sigma_iu_i \\ \sigma_i= \frac{Av_i}{u_i} \\ ④A=U\Sigma V^t \Rightarrow A^t=V\Sigma U^t,则A^tA=V\Sigma^2 V^t \\可见A^tA特征向量组成的矩阵的确为SVD中的V矩阵,且\sigma_i=\sqrt \lambda_i

AAm,n=UΣVtUm,m     Σm,n    Vn,nAtAn,n(AtA)vi=λiviAtAn,nUAAm,mtV,(AAt)ui=λiui(AtAAAt)Am,n=UΣVtAV=UΣAvi=σiuiσi=uiAviA=UΣVtAt=VΣUt,AtA=VΣ2VtAtASVDVσi=λ
i

几何意义:
A是行空间中的, Ax 一定落在 A 的列空间中,AA^t是列空间的一组基组成的,A = USV’,V’ 的含义是把列空间中的向量投影到 r 维子空间中,\Sigma再进行旋转操作(也可能包括翻转),第三步 U的含义就是把这个旋转逆过来,或者说把中介空间中的向量旋转回左、右两个 r 维子空间中去,但原先跟 r 维子空间垂直的分量就恢复不回来了。

A

T

A

A

A

T

A

m

n

,

m

m

A

A

T

其实上边的A^TA和AA^T 写反了,从矩阵的维度可知,A和\sum都是m*n,分解左边应为m*m,只能是A*A^T

ATAAATAmn,mmAAT

伪逆 (广义逆)

,

A

L

A

=

E

,

A

A

L

E

,

A

L

=

(

A

t

A

)

1

A

t

,

A

A

R

=

E

,

A

R

A

E

,

A

R

=

A

t

(

A

t

A

)

1

列满秩,左逆矩阵:A^LA=E,AA^L\neq E,A^L=(A^tA)^{-1}A^t\\ 行满秩,右逆矩阵:AA^R=E,A^RA\neq E,A^R=A^t(A^tA)^{-1}\\

,ALA=E,AAL=E,AL=(AtA)1At,AAR=E,ARA=E,AR=At(AtA)1

A

B

使

A

B

A

=

A

,

B

A

B

=

B

,

(

A

B

)

t

=

A

B

,

(

B

A

)

t

=

B

A

,

B

A

A

=

U

[

Σ

0

0

0

]

V

t

B

=

V

[

Σ

1

0

0

0

]

U

t

对于矩阵A,如果存在B使得:ABA=A,BAB=B,(AB)^t=AB,(BA)^t=BA,则B为A的伪逆\\ 若A=U\begin{bmatrix}\Sigma &0\\0&0\end{bmatrix}V^t,则B=V\begin{bmatrix}\Sigma^{-1}&0\\0&0\end{bmatrix}U^t

AB使ABA=A,BAB=B,(AB)t=AB,(BA)t=BA,BAA=U[Σ000]VtB=V[Σ1000]Ut
添加链接描述

PCA降维

确定几个正交分主方向综合起来模拟原来的矩阵:

A

U

m

,

r

Σ

r

,

r

V

r

,

n

t

A\simeq U_{m,r}\Sigma_{r,r} V_{r,n}^t

AUm,rΣr,rVr,nt


MATLAB

程序:

A=[3,2,2;2,3,-2]
B=A'
C=A'*A
d = eig(C)%特征值
[V,D] = eig(C)%计算C的特征值对角阵D和特征向量V,使AV=VD成立。
[U,S,V] = svd (C)   %返回一个与C同大小的对角矩阵S,两个酉矩阵U和V,且满足C= U*S*V'。若C为m×n阵,则U为m×m阵,V为n×n阵。奇异值在S的对角线上,非负且按降序排列。

结果:

C=   13    12     2
     12    13    -2
      2    -2     8
d =  -0.0000   9.0000   25.0000
---------------------------------------------------
U =

   -0.7071   -0.7071
   -0.7071    0.7071
S =

    5.0000         0         0
         0    3.0000         0
V =

   -0.7071   -0.2357   -0.6667
   -0.7071    0.2357    0.6667
   -0.0000   -0.9428    0.3333

版权声明:本文为ResumeProject原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。