当前位置:   article > 正文

线性代数|机器学习-P5特征值和特征向量

线性代数|机器学习-P5特征值和特征向量

1. 特征值和特征向量

1.1 特征向量

假设有一个n行n列的方阵A,有 n 个不相同的特征值为 λ \lambda λ,特征向量为 x 1 , x 2 , ⋯   , x n x_1,x_2,\cdots,x_n x1,x2,,xn.等式如下:
A x i = λ i x i , i = 1 , ⋯   , n → A 2 x = λ 2 x

Axi=λixi,i=1,,nA2x=λ2x
Axi=λixi,i=1,,nA2x=λ2x

  • 特征向量的好处在于,对于向量x来说, A x = λ x Ax=\lambda x Ax=λx,通过左乘矩阵A,还是不改变向量的方向,只是按照 λ \lambda λ倍进行缩放。
    A k x = λ k x
    Akx=λkx
    Akx=λkx
  • 对于微分方程来说
    d u d t = A u , e A t = e λ t
    dudt=AueAt=eλt
    dtdu=AueAt=eλt
  • 通解表示如下:
    u ( t ) = S e Λ t S − 1 u ( 0 ) = e A t u ( 0 )
    u(t)=SeΛtS1u(0)=eAtu(0)
    u(t)=SeΛtS1u(0)=eAtu(0)

1.2 向量分解

假设矩阵A有n个线性无关的特征向量,那么对于任意矩阵v来说,可以分解为特征向量的线性组合
v = c 1 x 1 + c 2 x 2 + ⋯ + c n x n

v=c1x1+c2x2++cnxn
v=c1x1+c2x2++cnxn

  • 两边同时乘以 A k , A k x = λ k x A^k,A^{k}x=\lambda^kx Ak,Akx=λkx:
    A k v = c 1 λ 1 k x 1 + c 2 λ 2 k x 2 + ⋯ + c n λ n k x n
    Akv=c1λ1kx1+c2λ2kx2++cnλnkxn
    Akv=c1λ1kx1+c2λ2kx2++cnλnkxn
  • 特征向量在差分方程上的应用
    u k + 1 = A u k → u k = A k u 0 = λ k x u 0
    uk+1=Aukuk=Aku0=λkxu0
    uk+1=Aukuk=Aku0=λkxu0

2. 矩阵相似

2.1 特征值求解法-相似

假设我们有两个矩阵A,B如果存在一个可逆矩阵M,满足如下关系,可推出A相似于B
B = M − 1 A M → B ∼ A → A 和 B 有相同的特征值

B=M1AMBAAB
B=M1AMBAAB有相同的特征值

  • 假设矩阵A的特征值为 λ \lambda λ,特征向量为x,
    ∣ B − λ I ∣ = ∣ M − 1 A M − λ I ∣ = ∣ M − 1 A M − M − 1 λ M ∣ = ∣ M − 1 ∣ ∣ A − λ I ∣ ∣ M ∣ = ∣ A − λ I ∣
    |BλI|=|M1AMλI|=|M1AMM1λM|=|M1||AλI||M|=|AλI|
    BλI=M1AMλI=M1AMM1λM=M1∣∣AλI∣∣M=AλI
  • 所以可得如下:
    B ∼ A ⇒ λ A = λ B
    BAλA=λB
    BAλA=λB
  • Matlab中如何求解特征值
    对于给定的矩阵A来说,我们用一个可逆矩阵 M 1 M_1 M1右乘矩阵A,左乘 M 1 − 1 M_1^{-1} M11,使得矩阵A逐渐变成上三角矩阵,通过不断地左右乘 M 1 , M 2 M1,M2 M1,M2,最后得到一个上三角矩阵B,这样我们就通过相似的形式得到主对角线上的特征值了。
    B = ( M n ⋯ M 2 M 1 ) − 1 A ( M n ⋯ M 2 M 1 ) → B U p T r i a n g l e ∼ A → A 和 B 有相同的特征值
    B=(MnM2M1)1A(MnM2M1)BUpTriangleAAB
    B=(MnM2M1)1A(MnM2M1)BUpTriangleAAB有相同的特征值

2.2 特殊特征值

假设我们有两个矩阵A,B,令AB的特征值为 λ A B \lambda_{AB} λAB,特征向量为x,令BA的特征值为 λ B A \lambda_{BA} λBA,证明 λ A B = λ B A \lambda_{AB}=\lambda_{BA} λAB=λBA

  • 根据定义可得:
    A B x = λ A B x
    ABx=λABx
    ABx=λABx
  • 两边同时乘以B可得:
    B A B x = λ A B B x → ( B A ) ( B x ) = λ A B ( B x ) → λ A B = λ B A
    BABx=λABBx(BA)(Bx)=λAB(Bx)λAB=λBA
    BABx=λABBx(BA)(Bx)=λAB(Bx)λAB=λBA

2.3 反对称矩阵

假设我们有一个矩阵A表示如下:
A = [ 0 1 − 1 0 ] → A T = − A

A=[0110]AT=A
A= 0110 AT=A

  • 矩阵A实现的功能是将向量x顺时针旋转90°。
    在这里插入图片描述
  • 求矩阵A的特征值和特征向量如下:
    λ 1 = i , v 1 = [ 1 i ] ; λ 2 = − i , v 1 = [ 1 − i ] ; S = [ 1 1 i − i ] ; Λ = [ i 0 0 − i ] ;
    λ1=i,v1=[1i];λ2=i,v1=[1i];S=[11ii];Λ=[i00i];
    λ1=i,v1= 1i ;λ2=i,v1= 1i ;S= 1i1i ;Λ= i00i ;
  • 分解A如下:
    A = S Λ S − 1 ⇒ [ 0 1 − 1 0 ] = [ 1 1 i − i ] [ i 0 0 − i ] [ 1 1 i − i ] − 1 ;
    A=SΛS1[0110]=[11ii][i00i][11ii]1;
    A=SΛS1 0110 = 1i1i i00i 1i1i 1;

3.对称矩阵

对称矩阵具有实数特征值和正交的特征向量。我们定义矩阵A如下:
A = [ 0 1 1 0 ] → λ 1 = 1 , v 1 = [ 1 1 ] ; λ 2 = − 1 , v 2 = [ − 1 1 ] ;

A=[0110]λ1=1,v1=[11];λ2=1,v2=[11];
A= 0110 λ1=1,v1= 11 ;λ2=1,v2= 11 ;

  • 可得如下:
    A = S Λ S − 1 → [ 0 1 1 0 ] = [ 1 1 1 − 1 ] [ 1 0 0 − 1 ] [ 1 1 1 − 1 ] − 1
    A=SΛS1[0110]=[1111][1001][1111]1
    A=SΛS1 0110 = 1111 1001 1111 1
声明:本文内容由网友自发贡献,不代表【wpsshop博客】立场,版权归原作者所有,本站不承担相应法律责任。如您发现有侵权的内容,请联系我们。转载请注明出处:https://www.wpsshop.cn/w/Monodyee/article/detail/672247
推荐阅读
相关标签
  

闽ICP备14008679号