線性代數 (下)

相似性

  • assume

    A,BFn×n,若有一可逆矩陣
    PFn×n
    使的
    P1AP=B
    則,
    AB

  • 等價關係
    AA

    AB
    BA

    AB
    BC
    AC

  • AαIA=αI,可知道
    αI
    只與自己相似

  • TL(V,V)
    β,γ
    V
    的基底,則
    [T]β[T]γ

Image Not Showing Possible Reasons
  • The image file may be corrupted
  • The server hosting the image is unavailable
  • The image path is incorrect
  • The image format is not supported
Learn More →
  • if

    AB, then
    (1)
    tr(A)=tr(B)

    (2)
    det(A)=det(B)

    (3)
    rank(A)=rank(B)

    (4)
    nullity(A)=nullity(B)

    (5)
    ATBT

    (6)
    AkBk,kN

    (7)
    αAαB,αF

    (8)
    A+αIB+αI

    (9)
    f(A)f(B),f(x)F[x]

    (10)
    A1B1
    ,if they are invertible
    (11)
    pA(x)=pB(x)
    (有相同的特徵根,但特徵向量不一定相同)
    上述1~4皆成立的話,也不保證
    AB
    ,但要是有任何一項不成立,則代表
    A
    B
    不相似

  • 相似矩陣有相同的特徵根,特徵向量不一定相同

  • Image Not Showing Possible Reasons
    • The image file may be corrupted
    • The server hosting the image is unavailable
    • The image path is incorrect
    • The image format is not supported
    Learn More →

不變子空間

  • 不變子空間:
    假設

    TL(V,V),W為V的子空間且
    T(W)W,wW,T(w)W

  • TL(V,V)

    1. {0},V,ker(T),Im(T)
      皆為T-不變子空間
    2. T在W上的限制函數,
      TW:WV
      定義為
      TW(v)=T(v),vW
      ,可視為
      TW:WW,TWL(W,W)
  • {0},V,ker(T),Im(T)都是T-invariant subspace

  • T:VV,所有的T-invariant subspace的交集也是T-invariant subspace,和空間也是

  • Image Not Showing Possible Reasons
    • The image file may be corrupted
    • The server hosting the image is unavailable
    • The image path is incorrect
    • The image format is not supported
    Learn More →

特徵根及特徵向量

  • eigen value與eigen vector就對應到不變子空間,投影後仍在同一個subspace上

  • 相對於

    λ的特徵向量的任意線性組合,若不為0向量,則仍為
    λ
    的特徵向量

  • det(T)=det([T]β)
    tr(T)=tr([T]β)

  • λ
    A
    的特徵根
    det(AλI)=0

  • pA(x)=charA(x)=det(AxI)

  • pA(0)=det(A),pA(x)
    (x)n1
    的係數為
    tr(A)

  • Image Not Showing Possible Reasons
    • The image file may be corrupted
    • The server hosting the image is unavailable
    • The image path is incorrect
    • The image format is not supported
    Learn More →
  • V(λ)={vV|T(v)=λv},稱為T相對於
    λ
    的特徵空間,
    V(λ)
    是T的不變子空間

  • V(λ)=ker(TλI)

  • Ax=λx,則若
    A=[T]β,x=[v]β
    ,則
    [T]β=λ[v]β

  • V(λ1),..,V(λk)為獨立子空間,
    viV(λi)
    v1,...,vk
    為線性獨立 (證明)

  • AB,則
    pA(x)=pB(x)
    ,具有相同特徵根(特徵向量不一定),反之未必成立

  • A,AT具有相同特徵根但未必具有相同特徵向量

  • A,A1,Am,f(A)有相同的特徵根

  • AB,BA具有相同特徵根,且連特徵多項式也相同

  • Image Not Showing Possible Reasons
    • The image file may be corrupted
    • The server hosting the image is unavailable
    • The image path is incorrect
    • The image format is not supported
    Learn More →
  • A為可逆
    0不為
    A
    的特徵根
    A
    可逆且
    Ax=λx
    ,則
    A1x=1λx

    mZ+,Amx=λmx

    f(x)P
    或馬式級數,則
    f(A)x=f(λ)x

對角化

  • 假設

    TL(V,V)為一個線性映射且存在
    V
    的一組基底
    β
    使得
    [T]β=D
    為對角矩陣,則稱
    T
    為可對角化

  • 假設

    AFn×n,若存在一可逆矩陣
    PFn×n
    使得
    P1AP=D
    為對角矩陣,則稱
    A
    為可對角化

  • 有些人會誤以為方形矩陣存在rref就代表可對角化(例如我) 但其實不是這樣

  • T可對角化
    T
    dim(V)
    個線性獨立特徵向量

  • T可對角化
    pT(x)
    在F中可分解且
    am(λi)=gm(λi),i=1,2,..,r
    V=V(λ1)V(λ2)...V(λr)

  • β中的元素為
    T
    的n個線性獨立特徵向量且
    D
    的對角項元素依序為這n個特徵向量對應的根,用矩陣
    A
    來看的話則是
    P
    的行向量為
    A
    的n個線性獨立特徵向量

  • λ
    pT(x)
    中的重根數
    λ
    的代數重數
    am(λ)

  • dim(V(λ))為幾何重數
    gm(λ)

  • gm(λ)=dim(V(λ))=dim(ker(AλI))=nullity(AλI)=nrank(AλI)

  • 1gm(λ)am(λ)n,n是
    pT(x)
    為n次多項式的意思

  • 可分解: 假設

    f(x)F[x]
    f(x)
    可完全分解成佈於
    F
    一次因式乘積,則稱
    f(x)
    is split over
    F

  • 假設

    AFn×n
    pA(x)
    is split over
    F
    pA(x)=(λ1x)...(λnx)
    ,則
    det(A)=λ1...λn

    tr(A)=λ1+...+λn

  • abb...bbab...bbba...b......bbb...a
    有兩個eigen value,
    λ1=ab,λ2=a+(n1)b

    tr(A)
    來解

  • AFn×n,如果
    A
    有n個相異eigen value,則A可對角化

  • 假設

    A,B可同步對角化,則
    AB=BA
    ,逆命題未必成立
    假設
    A,B
    皆可對角化,則上述逆命題成立

  • Image Not Showing Possible Reasons
    • The image file may be corrupted
    • The server hosting the image is unavailable
    • The image path is incorrect
    • The image format is not supported
    Learn More →

對角化的應用

  • f(A)時,先做
    A=PDP1
    ,之後得
    f(A)=Pf(D)P1
    ,但
    A
    的所有eigen value都要落在
    f(x)
    的收斂區內

  • 如果

    XFn×n,Xn=A,則定義
    A1n=X
    ,先將
    A
    對角化得
    D
    再求出
    D1n

  • 可以解微分方程系統,

    xˊ=Ax,且
    A
    可對角化,
    v1,...,vn
    A
    的n個線性獨立eigen vector,此時
    x=i=1ncieλ1tvi

  • 可以解線性遞迴關係式

冪等算子與矩陣

  • TL(V,V),若
    T2=T
    ,則
    T
    V
    上的Idempotent operator或稱為projection operator,同理套用在矩陣上

  • T(A)=A+AT2是一個idempotent operator

  • 假設

    TL(V,V),則
    ker(T)ker(T2)

    Im(T2)Im(T)

    且下列敘述等價
    (1)
    Im(T2)=Im(T)

    (2)
    rank(T2)=rank(T)

    (3)
    nullity(T2)=nullity(T)

    (4)
    ker(T2)=ker(T)

  • Sylvester第二定理
    假設

    TL(V,V),則下列敘述等價
    (1)
    V=ker(T)Im(T)

    (2)
    V=ker(T)+Im(T)

    (3)
    ker(T)Im(T)={0}

    要證明(1)只需要證明(2)或(3)

  • 假設

    T
    V
    上的idempotent operator,則
    V=ker(T)Im(T)

    假設
    AFn×n
    為idempotent matrix,則
    Fn×1=ker(A)CS(A)

  • eigen value不是0就是1

    V(0)=ker(T),V(1)=Im(T)
    可對角化

  • 對於冪等算子而言,

    rank(T)=tr(T)

冪零算子

  • 假設

    T=L(V,V)且存在正整數k使得
    Tk=O
    ,則
    T
    為nilpotent operator,且最小的k為
    T
    的index

  • 零矩陣

    O為冪零算子,單位矩陣
    I
    不為冪零算子

  • 假設

    T為冪零算子,且
    V=span(β)
    ,則
    [T]β
    為冪零算子

  • 嚴格上三角或下三角矩陣都是index = n冪零矩陣

  • 下列敘述等價
    (1)

    T為index=
    k
    的冪零算子
    (2)
    Tk=O
    Tk1O

    (3)
    ker(Tk)=V
    ker(Tk1)V

    (4)
    vV,Tk(v)=0
    且存在
    uV,Tk1(u)0

  • k階下移矩陣:

    (Sk)ij=1,i=j+1
    (Sk)ij=0,ij+1

    (Sk)ijt=1,i=j+t

    (Sk)k=O

  • ker(Tk)
    Im(Tk)
    都為T-不變子空間

  • 局部冪零:

    TW為冪零算子,局部可逆:
    TW
    為局部可逆
    W
    為使
    TW
    為冪零算子的最大不變子空間,則
    W
    為最大冪零區
    W
    為使
    TW
    為可逆算子的最大不變子空間,則
    W
    為最大可逆區

  • ker(Tk)=ker(Tk+1),則
    ker(Tk+1)=ker(Tk+2)

    Im(Tk)=Im(Tk+1)
    ,則
    Im(Tk+1)=Im(Tk+2)

  • ker(Tk)ker(Tk+1),Im(Tk)Im(Tk+1)
    因此
    {0}ker(T)...ker(Ti)...V

    Im(T)Im(T2)⊇⊇...Im(Ti)...

  • ker(Tk)=ker(Tk+1)Im(Tk)=Im(Tk+1)

  • 存在最小正整數k使得

    ker(Tk)=ker(Tk+1)=...,也存在這樣的k使得
    Im(Tk)=Im(Tk+1)=...

    W=i=1ker(Ti)=ker(Tk)
    且W為最大冪零區且
    TW
    的指標為k
    W=i=1Im(Ti)=Im(Tk)
    且W為最大冪零區且
    TW
    的指標為k
    上述定理告訴我們,T的次方愈大,則它的核集愈大,但不會一直變大
    次方愈大,像集愈小,但不會一直變小

  • TL(V,V),dim(V)=n<,則存在一正整數k使得
    V=ker(Tk)Im(Tk)

    ker(T)=ker(T2)V=ker(T)Im(T)

  • 所以要是T為冪零算子,則

    V=ker(T)Im(T)
    det(T)=0,Tn=O

    T的eigen value
    λ=0
    所以T的矩陣A保證為不可逆矩陣
    反之不一定成立,但要是A為可逆,則A為不為可冪矩陣

  • T為具有指標k的冪零算子,則存在

    vker(Tk)ker(Tk1)使得
    {v,T(v),...,Tk1(v)}
    為線性獨立集

循環子空間及循環分解

  • T-cyclic subspace generated by

    v,
    Cv(T)=span{v,T(v),T2(v),..,Tk(v),...}

  • v=0Cv(T)={0}dim(Cv(T))=0

  • Cv(T)
    T
    不變子空間

  • TL(V,V),vker(Tk)ker(Tk1),β={v,T(v),...,Tk1(v)}
    W=Cv(T)
    的一組基底且
    [TW]β=Sk
    β
    稱為
    W
    的循環基底
    如果
    β={Tk1(v),...,v}
    ,則
    [TW]β=SkT

  • 只要

    dim(W)=k就可保證
    β={v,T(v),...,Tk1(v)}
    W
    的基底,但此時
    Tk(v)
    未必為0

循環分解定理 cyclic decomposition theorem

TL(V,V)為具有指標k的冪零算子,則存在
v1,...,vrV
使得
V=Cv1(T)Cvr(T)

其中
dim(Cvi(T))=ni,n1n2nr
,且
viker(Tni)ker(Tni1)

dim(ker(T))=r,k=n1

Jordan form

  • 假設

    TL(V,V),dim(V)=n<,dim(i=1ker(Ti))=m

    1. m=0
      ,0不為
      pT(x)=0
      的根
    2. m1
      ,0為
      pT(x)
      的m重根
  • 假設

    TL(V,V),dim(V)=n<

    1. λ
      不為
      pT(x)=0
      的根
      0
      不為
      pTλI(x)=0
      的根
    2. λ
      pT(x)=0
      的m重根
      0
      pTλI(x)=0
      的m重根
  • TL(V,V),dim(V)=n<,λ
    T
    的特徵根,
    W=i=1ker((TλI)i)
    ,若
    dim(W)=n

    1. 因為0為
      pTλI(x)
      的m重根,所以
      λ
      pT(x)
      的m重根,
      am(λ)=m
    2. (TλI)W
      為冪零算子
  • 廣義特徵空間

    K(λ)=i=1ker((TλI)i),T相對於
    λ
    的廣義特徵空間
    也可定義為
    K(λ)={vV|(TλI)p(v)=0}
    ,for some
    pZ+

    廣義特徵向量
    vK(λ),v0
    ,v為T相對於
    λ
    的廣義特徵向量
    也可定義為
    (TλI)p(v)=0

  • K(λ)
    V
    的子空間也為T-不變子空間

  • Jk(λ)=Sk+λIk為k階Jordan矩陣

  • u,λ
    T
    的相異特徵根

    1. K(u)
      (TλI)i
      -不變子空間 (證明)
    2. (TuI)K(λ)
      為一對一函數 (證明)
    3. ((TuI)i)K(λ)
      為一對一函數

Cayley Hamilton定理及應用

  • Cayley-Hamilton定理

    AFn×n,f(x)=pA(x)=det(AxI),則
    f(A)=O

  • TL(V,V),W為T-invariant,則
    Tw
    的特徵多項式整除
    T
    的特徵多項式

  • pA(x)可分解,則,
    A
    為idempotent matrix
    A
    的所有特徵根都為0 (證明)

  • AFn×n,g(x)P,若
    g(x)
    除以
    pA(x)
    的餘式為
    r(x)
    ,則
    g(A)=r(A)

  • A為可逆,可以證明
    A1span{I,A,A2,,An1}

內積

  • <,>:V×VF

    1. u,v,wV,<u+v,w>=<u,w>+<v,w>
    2. u,vV,αF,<αu,v>=α<u,v>
    3. u,vV,<u,v>=<v,u>
    4. vV,v0,<v,v>>0

      則稱
      <,>
      為V上的內積
  • 要證明內積,只須證明以下三個

    1. <αu+βv,w>=α<u,w>+β<v,w>
    2. <u,v>=<v,u>
      ,若
      F=R
      ,則證明
      <u,v>=<v,u>
    3. <v,v>=0v=0

Gram-Schmidt正交化及QR分解

  • S={v1,...,vk}V為一個不含零向量的正交集,若
    vspan(S),v=i=1kαivi,αj=<v,vj><vj,vj>=<v,vj>||vj||2

  • aj
    v
    相對於
    S
    的傅立葉級數

  • SV為不含零向量的正交集,則
    S
    為線性獨立集

  • SV為線性獨立集,令
    u1=v1,uk=vki=1k1<vk,ui><ui,ui>ui
    ,則
    {u1,..,un}
    為不含零向量的正交集且
    span{v1,...,vn}=span{u1,...,un}

  • AFm×n,rank(A)=n,A=QR,QFm×n,他的行向量形成單範正交集,
    R
    為可逆的上三角矩陣

  • 任何矩陣皆可做QR分解

正交投影

  • W
    V
    的子空間,
    vV
    ,若存在
    v0W
    使
    <vv0,w>=0,wW
    ,則稱
    v0
    v
    W
    上的正交投影向量

  • β={v1,...,vk}
    W
    的一組基底,
    <v,w>=0,wW⇔<v,vi>=0,viβ

  • β={v1,...,vk}
    W
    的一組正交基底,
    v
    W
    上存在唯一的正交投影向量
    projWv=v0=i=1k<v,vi><vi,vi>vi
    (第一種找正交投影的方法)

  • W
    V
    的子空間,
    vWprojWv=v

    projW0=0

    vW,projWv=0

  • P:VV,P(v)=projWv,vV

    1. P
      為線性映射
    2. vW,P(v)=v
    3. Im(P)=W
    4. P2=P
  • W
    V
    的子空間,
    β={v1,..,vk}
    W
    的一組正交基底

    1. ||P(v)||2=i=1k|<v,vi>|2||vi||2,vV
      (證明)
    2. ||P(v)||||v||,vV
      (Bessel不等式)
    3. ||P(v)||=||v||,vW
  • ||vP(v)||||vw||,wW
    v
    P(v)
    的距離永遠小於等於v與w的距離
    P(v)
    W
    上與
    v
    最靠近的向量
    ,也是
    v
    W
    上的最佳逼近best approximation

  • ker(AHA)=ker(A)
    rank(AHA)=rank(A)

    ARm×n,AH=AT,rank(ATA)=rank(AAT)

  • 假設

    AFm×n
    A
    為行獨立
    AHA
    為可逆矩陣
    W=R(A),bFm×1,projWb=A(AHA)1AHb

  • 假設

    AFm×n,W=R(A),bFm×1
    xFn×1
    使得
    ||Axb||
    為最小
    AHAx=AHb
    (證明)

  • AHAx=AHb稱為正規方程式,必定有解,當A不為行獨立時,
    rank(AHA)=rank(A)<n
    ,則有無限多組解,但
    projWb=Ax
    必定唯一,也就是可能解出很多
    x
    (最小平方解),但
    Ax
    必定唯一,這個
    Ax
    就是
    b
    CS(A)
    的正交投影,

  • Image Not Showing Possible Reasons
    • The image file may be corrupted
    • The server hosting the image is unavailable
    • The image path is incorrect
    • The image format is not supported
    Learn More →
  • Ax就是
    projWb
    了,為甚麼還要用那個麻煩的公式? 因為你不知道
    x
    ,只知道b,要從b得到x就要解正規方程式,直接套那個公式比較快

  • Ax=b為inconsistent的時候,可以用最小平方解來得到最近的答案

  • 求最小平方直線也是用相同方法

  • AFm×n,rank(A)=n,bFm×1,A=QR,則
    ||Axb||
    為最小
    Rx=QHb
    (證明)

  • AFm×n為行獨立,
    P=A(AHA)1AH
    為正交投影矩陣,則
    P2=P,PH=P

    若反過來,
    PFn×n,P2=P,PH=P
    ,則
    P
    為投影在
    W=R(P)
    上的正交投影矩陣 (證明)

  • QFm×n,其中
    Q
    的行向量為單範正交集,
    W=P(Q)
    ,則
    projWb=QQHb
    ,投影在
    R(Q)
    的正交投影矩陣為
    QQH

  • Image Not Showing Possible Reasons
    • The image file may be corrupted
    • The server hosting the image is unavailable
    • The image path is incorrect
    • The image format is not supported
    Learn More →

正交補空間

  • V
    為佈於
    F
    的內積空間,
    SV
    ,定義
    S={vV|<v,s>=0,sS}
    ,稱為S的正交補空間

注意!!! 沒有規定

S必須是subspace,只需要是
V
的subset即可

  • SV

    1. {0}=V
      V={0}
    2. S
      V
      的子空間
    3. SS⊥⊥,S⊥⊥=(S)
      ,如果S為一個subspace,則
      S=S⊥⊥
    4. SS={{0} if 0S if 0S
      (證明)
  • Rn中,一條方程式決定的平面稱為hyperplane,若hyperplane通過原點時,它為一個子空間
    W
    ,,此hyperplane的法向量會是
    W
    的basis

  • P
    V
    W
    上的正交投影算子,
    ker(P)=W

  • V=WW (證明)
    V=WW⊥⊥

  • 根據以上定理可知

    1. dim(V)=dim(W)+dim(W)
    2. W⊥⊥=W
    3. vV
      ,存在唯一
      wW,uW
      使得
      v=w+u
  • P
    V
    W
    上的正交投影算子,
    projWv=vprojWv
    (IP)
    V
    W
    上的正交投影算子

  • WiWj,ij,稱
    W1,...,Wk
    為正交子空間

  • W1,..,Wk為正交子空間,則她們也為獨立子空間

  • ARm×n

    1. R(AT)=N(A)
      (證明)
    2. R(A)=N(AT)
    3. N(A)=R(AT)
    4. N(AT)=R(A)
  • 根據以上可以得出

    rank(A)=n
    A
    行獨立,所以
    R(A)
    為一個subspace
    A(ATA)1AT
    Rm
    投影在
    R(A)
    的正交投影矩陣
    IA(ATA)1AT
    Rm
    投影在
    N(AT)
    的正交投影矩陣
    rank(A)=m
    A
    列獨立
    AT(AAT)1A
    Rn
    投影在
    R(AT)
    的正交投影矩陣
    IAT(AAT)1A
    Rn
    投影在
    N(A)
    的正交投影矩陣

  • ARm×n,bRm×1
    Ax=b
    有解
    y,ATy=0
    ,則
    bTy=0

  • 極小解

    Ax=b若有解,極小解
    s
    是指
    As=b
    ||s||||u||,u,Au=b

    也就是所有解當中長度最小的

  • AFm×n,bFm×1,若
    Ax=b
    有解

    1. 唯一存在
      sR(AH)
      使得
      s
      為極小解
    2. (AAH)u=b
      ,則
      s=AHu

內積上的算子及其應用

伴隨算子

  • V為內積空間,
    TL(V,V)
    ,若存在
    T:VV
    為線性映射使得
    <T(u),v>=<u,T(v)>,u,vV
    ,稱
    T
    T
    的adjoint operator

  • f:VF為一個線性泛函,存在唯一
    vV
    使得
    f(u)=<u,v>,uV

  • T存在唯一

  • β={v1,,vn}
    V
    的單範正交基底,則
    [T]β=[T]βH
    (可用此方法求伴隨算子)

  • <Ax,y>=<x,AHy>,x,yFn×1

  1. (αT+βU)=αT+βU
  2. (T)=T
  3. (TU)=UT
  4. I=I
  5. T
    為可逆函數,則
    (T)1=(T1)
  6. λ
    T
    的eigen value,則
    λ
    T
    的eigen value,但不保證eigen vector相等

正規算子

  • TT=TT,T稱為normal operator
    AHA=AAH
    ,A稱為normal matrix

  • 旋轉矩陣為正規矩陣

  • 正規算子會有以下特性

    1. |T(v)|=|T(v)|,vV
    2. cF,TcI
      為正規算子
    3. Tv=λv
      ,則
      T(v)=λv
    4. Tv1=λ1v1,Tv2=λ2v2,λ1λ2
      ,則
      v1v2
  • T=T,稱
    T
    為自伴算子(self-adjoint operator)
    所以
    [T]β=[T]β
    ,為Hermitian matrix
    T=T
    ,稱
    T
    為斜自伴算子(skew self-adjoint operator)
    所以
    [T]β=[T]β
    ,為skew-Hermitian matrix

  • ACn×n為Hermitian matrix,若
    ARn×n
    ,以下也同樣成立

    1. A為normal matrix
    2. A的eigen value為實數
    3. A的相異eigen value對應的eigen vector必正交
    4. A的主對角項元素皆為實數
    5. det(A)R
  • ACn×n未必保證
    A
    為正規矩陣

  • ACn×n為skew-Hermitian matrix,若
    ARn×n
    ,以下也同樣成立

    1. A為normal matrix
    2. A的eigen value為0或純虛數
    3. A的相異eigen value對應的eigen vector必正交
    4. A的主對角項元素皆為0或純虛數
    5. 當n為偶數時,
      det(A)R
      ,當n為奇數時,
      det(A)
      為0或純虛數
  • TT=I
    F=C
    時稱為么正算子unitary operator
    F=R
    時稱為正交算子orthogonal operator
    正交矩陣,
    A1=AT

    么正矩陣,
    A1=AH

  • ACn×n,A1=AH

    1. A為normal matrix
    2. |λ|=1
    3. 相異eigen value對應的eigen vector必正交
    4. |det(A)|=1
  • ACn×n,A1=AT

    1. A為normal matrix
    2. λ=±1
    3. 相異eigen value對應的eigen vector必正交
    4. det(A)=±1

么正及正交算子的特性

  • TL(V,V),V為複內積空間

    1. <T(u),v>=0,u,vV
      ,則
      T=O
      (在實內積空間也成立)
    2. <T(v),v>=0,vV
      ,則
      T=O
      (證明) (在實內積空間不成立)
  • TL(V,V),V為複內積空間,以下等價 (實空間也成立)

    1. T
      為么正算子
    2. T
      保內積,
      <T(u),T(v)>=<u,v>,u,vV
    3. T
      保長度,
      |T(v)|=|v|,vV
  • polar identity

    V為佈於
    F
    的內積空間
    F=R,<u,v>=14|u+v|214|uv|2,u,vV

    F=C,<u,v>=14k=14ik|u+ikv|2,u,vV

  • A,BFn×n,下列敘述等價

    1. A=B
    2. yHAx=yHBx,x,yFn×1

    F=C

    1. A=B
    2. xHAx=xHBx,xCn×1
  • ACn×n,下列等價

    1. A
      為么正矩陣
    2. A
      保內積,
      <Ax,Ay>=<x,y>,x,yCn×1
    3. A
      保長度,
      |Ax|=|x|,xCn×1
  • ARn×n,下列等價

    1. A
      為正交矩陣
    2. A
      保內積,
      <Ax,Ay>=<x,y>,x,yRn×1
    3. A
      保長度,
      |Ax|=|x|,xRn×1
  • TL(V,V),V為複內積空間,下列等價

    1. T
      為么正算子
    2. β
      V
      的orthonormal basis,則
      T(β)
      也為
      V
      的orthonormal basis
    3. 存在
      V
      的一個orthonormal basis使得
      T(β)
      也為
      V
      的orthonormal basis
  • ACn×n,下列等價

    1. A
      為么正矩陣
    2. β={v1,,vn}
      Cn×1
      的orthonormal basis,則
      {Av1,,Avn}
      也為
      Cn×1
      的orthonormal basis
    3. 存在
      Cn×1
      的一個orthonormal basis
      β={v1,,vn}
      使得
      {Av1,,Avn}
      也為
      Cn×1
      的orthonormal basis
    4. A
      的行向量形成orthonormal set
    5. A
      的列向量形成orthonormal set
  • ARn×n,下列等價

    1. A
      為么正矩陣
    2. β={v1,,vn}
      Rn×1
      的orthonormal basis,則
      {Av1,,Avn}
      也為
      Rn×1
      的orthonormal basis
    3. 存在
      Rn×1
      的一個orthonormal basis
      β={v1,,vn}
      使得
      {Av1,,Avn}
      也為
      Rn×1
      的orthonormal basis
    4. A
      的行向量形成orthonormal set
    5. A
      的列向量形成orthonormal set

正定及正半定算子與矩陣

  • ACn×n
    A
    為Hermitian matrix
    xHAxR,xCn×1

    TL(V,V)
    V
    為佈於
    C
    的內積空間
    T
    為self-adjoint operator
    ⇔<T(v),v>∈R,vV

  • 二次式 quadratic form

    ACn×n
    Q(x)=xHAx,xCn×1

    ARn×n
    Q(x)=xTAx,xRn×1

  • 任何矩陣都有一個對應的對稱矩陣其二次式完全一樣 (這樣不就相等了?)

  • 以下定義需要滿足

    Q(x)R,也就是
    A
    為Hermitian
    positive definite matrix(operator):
    ACn×n,AH=A,Q(x)>0,x0
    T
    為自伴算子且
    <T(v),v> >0,v0

    positive semidefinite matrix(operator):
    ACn×n,AH=A,Q(x)0,x0
    T
    為自伴算子且
    <T(v),v> 0,v0

  • ARn×n,則要改為對稱矩陣

  • ACn×n為正定矩陣

    1. A
      為normal matrix,因為
      A
      為Hermitian matrix
    2. A
      為正半定矩陣
    3. 所有eigen value都為正
    4. 所有相異eigen value對應的eigen vectors都正交(normal matrix特性)
    5. A
      為可逆矩陣(證明)
    6. 主對角項元素皆為正(證明)
  • ACn×n為正半定矩陣

    1. A
      為normal matrix,因為
      A
      為Hermitian matrix
    2. 所有eigen value都為非負
    3. 所有相異eigen value對應的eigen vectors都正交(normal matrix特性)
    4. 主對角項元素皆為非負(證明)

么正及正交對角化

  • unitarily similar

    A,BCn×n,存在unitary matrix
    PHP=I
    使得
    PHAP=B

    orthogonally similar
    A,BRn×n
    ,存在orthogonal matrix
    PTP=I
    使得
    PTAP=B

  • 算子的么正對角化

    TL(V,V),若存在
    V
    的一組單範正交基底
    β
    使得
    [T]β
    為對角矩陣,則稱
    T
    為unitarily diagonalizable

  • Schur's定理

    TL(V,V),若
    pT(x)
    F
    中可分解,則存在
    V
    的一組單範正交基底
    β
    使得
    [T]β
    為上三角矩陣
    ACn×n
    ,則存在一個么正矩陣
    PCn×n
    使得
    PHAP
    為上三角矩陣
    ARn×n
    pA(x)
    R
    中可以分解,則存在一個正交矩陣
    PRn×n
    使得
    PTAP
    為上三角矩陣

  • V為佈於
    C
    的內積空間,若
    T
    為normal operator,則
    T
    可么正對角化
    V
    為佈於
    R
    的內積空間,若
    T
    為self-adjoint operator,則
    T
    可正交對角化
    ACn×n
    A
    可么正對角化
    A
    為normal operator
    ARn×n
    A
    可正交對角化
    A
    為對稱矩陣

  • ACn×n
    A
    為正規且上三角矩陣
    A
    為對角矩陣
    ARn×n
    A
    為對稱且上三角矩陣
    A
    為對角矩陣

  • ACn×n
    A
    為正規矩陣
    A
    可么正對角化(Hermitian矩陣一定可以么正對角化)
    ARn×n
    A
    為對稱矩陣
    A
    可正交對角化

  • 光譜定理spectral theorem

    TL(V,V),λ1,,λk為其相異eigen value,
    F=C
    時假設
    T
    為normal operator,
    F=R
    時假設
    T
    為self adjoint operator,
    Wi=V(λi),Ti
    V
    投影在
    Wi
    上的正交投影算子

    1. V=W1W2Wk
    2. Wi=jiWj
    3. TiTj=δijTi
      (證明)
    4. T1+T2++Tk=I
      (證明)
    5. T=λ1T1+λ2T2++λkTk
      (證明)

正定及正半定矩陣的特性

  • principal minors

    Image Not Showing Possible Reasons
    • The image file may be corrupted
    • The server hosting the image is unavailable
    • The image path is incorrect
    • The image format is not supported
    Learn More →

  • LDU分解:

    A=LDU=LDLT
    L
    為對角線皆為1的下三角矩陣,
    D
    為對角矩陣
    先做LU分解,再把U的對角項拿出來變成D

  • ARn×n為對稱矩陣(代表可以正交對角化),則以下等價

    1. A
      為正定
    2. A
      的所有eigen values都為正
    3. Δk(A)>0,k=1,,n
  • 可用以上定理檢查一個矩陣是否正定,若該矩陣不是對稱矩陣,可以先把他轉為對稱矩陣

  • ARn×n為對稱矩陣,若
    A
    為正定矩陣,則
    A=LDLT
    ,L為對角線皆為1的下三角矩陣且
    D
    為對角項元素皆為正的對角矩陣 (特徵根皆為正)

  • ACn×n,AH=A

    1. A
      為正定
      A
      的所有eigen values都為正
    2. A
      為正半定
      A
      的所有eigen values都為非負
  • Cholesky分解
    假設

    ARn×n為對稱矩陣,若
    A
    為正定矩陣,則
    A=LD(L)T=LDD(L)T=LDT(LDT)T=LLT
    ,其中
    L=LDT
    為下三角矩陣且對角項元素皆為正

  • ACn×n
    A
    為正定矩陣
    存在
    BCm×n,rank(B)=n,A=BHB

  • ARn×n
    A
    為對稱矩陣
    A
    為正定矩陣
    存在
    BRm×n,rank(B)=n,A=BTB

  • 假設

    ACn×n,則
    A
    為正半定矩陣
    存在
    BCm×n
    使得
    A=BHB

    假設
    ARn×n
    ,則
    A
    為正半定矩陣
    存在
    BRm×n
    使得
    A=BTB

  • 假設

    ACm×n,則
    AHA,AAH
    都是正半定矩陣
    假設
    ARm×n
    ,則
    ATA,AAT
    都是正半定矩陣

  • ACn×n,則
    A
    為正定矩陣
    存在
    B
    為正定矩陣使得
    A=B2

    A
    為正半定矩陣
    存在
    B
    為正半定矩陣使得
    A=B2

    ARn×n
    ,則
    A
    為正定矩陣
    存在
    B
    為正定矩陣使得
    A=B2

    A
    為正半定矩陣
    存在
    B
    為正半定矩陣使得
    A=B2

二次式的應用

  • 主軸定理principal axis theorem
    假設

    ACn×n為Hermitian matrix,
    Q(x)=xHAx=λ1|y1|2++λn|yn|2

    假設
    ARn×n
    為symmetric matrix,
    Q(x)=xHAx=λ1y12++λnyn2

  • Rayleigh商式

    ACn×n為Hermitian matrix,
    ARn×n
    為對稱矩陣,
    ρ(x)=Q(x)||x||2,x0
    稱為Rayleigh quotient

  • Rayleigh principle

    AFn×n
    F=C
    時為Hermitian matrix,
    F=R
    時為對稱矩陣,
    λ1λ2λn
    A
    的特徵根

    1. λ1ρ(x)λn
    2. maxx0ρ(x)=λn
    3. minx0ρ(x)=λ1
  • x
    A
    相對於
    λn
    的特徵向量
    xn
    ,則
    ρ(xn)=λn
    ,所以
    maxρ(x)=ρ(xn)

    x
    A
    相對於
    λ1
    的特徵向量
    x1
    ,則
    ρ(x1)=λ1
    ,所以
    maxρ(x)=ρ(x1)

Householder 轉換

  • uRn×1,||u||=1,稱
    H=I2uuT
    為相對於u的Householder matrix,也稱為elementary reflextor

    1. H
      為對稱矩陣
    2. H
      為正交矩陣,所以
      H1=HT=H
  • x,yRn×1,xy,||x||=||y||,u=xy||xy||,H=I2uuT,則

    1. ||xy||2=2(xy)Tx
    2. Hx=y (證明)
  • 對任意

    xRn×1,若
    α=||x||,y=αe1,u=xy||xy||,H=I2uuT
    ,則
    Hx=αe1

    就是x經過Householder轉換後,第一項之外元素都變成0

  • xRn×1,則存在一個Householder矩陣使得
    Hx=α1e1++αkek
    ,也就是
    Hx
    的後面n-k項為0

奇異值分解

  • σi,i=1,2,,s稱為
    A
    的奇異值singular value

  • ACm×n,s={m,n},若
    A=UΣVH,UCm×m,VCn×n,ΣCm×n

    其中
    UHU=I,VHV=I
    為么正矩陣,
    (Σ)ij=0,ij
    (Σ)ii=σi,σ1σ2σs

  • ARm×n,s={m,n},若
    A=UΣVH,URm×m,VRn×n,ΣRm×n

    其中
    UTU=I,VTV=I
    為正交矩陣,
    (Σ)ij=0,ij
    (Σ)ii=σi,σ1σ2σs

  • Σ可以不為方陣

  • A為方陣,
    U=V
    時,
    A=VΣVH
    ,此為么正對角化

  • A的奇異值分解步驟如下
    假設
    U=[u1  u2    um],V=[v1  v2    vn]

    1. 首先找

      AHA的eigen value
      λ1,,λn

      A
      的singular value為
      σi=λi
      ,記得要照大小排好

    2. V,因為
      V
      的行向量都是
      AHA
      的eigen vectors且為單範正交集,所以用
      λi
      V(λi)=ker(AHAλiI)
      的基底並做成單範正交集

    3. U分為兩部分,首先因為知道
      AV=UΣ
      ,所以
      Avi=σiui
      ,所以可先找
      u1,,un

      再來
      un+1,,um
      形成
      N(AH)
      的單範正交基底,所以解
      AH=0
      並化成單範正交基底就行

  • 經過證明可以知道任何矩陣都可以做SVD分解,並且有以下幾個特性

    • AHA
      的eigen values是
      λ1,,λn
      ,則
      A
      的singular value就是
      σi=λi
    • A
      的奇異值唯一,但
      U
      V
      未必唯一
    • V
      的行向量都是
      AHA
      的eigen vectors且形成一個單範正交集
    • U
      的行向量都是
      AAH
      的eigen vectors且形成一個單範正交集
    • Avi=σiui,i=1,,n
      ,因為
      AV=UΣ
    • AHui=σivi,i=1,,n
      並且
      AHui=0,i=n+1,,m
    • vi
      稱為
      A
      的right singular vector,
      ui
      稱為
      A
      的left singular vector
    • rank(A)=r
      ,則
      • v1,,vr
        形成
        R(AH)
        的單範正交基底
      • vr+1,,vn
        形成
        N(A)
        的單範正交基底
      • u1,,ur
        形成
        R(A)
        的單範正交基底
      • ur+1,,um
        形成
        N(AH)
        的單範正交基底
    • AHA
      AAH
      有相同的non zero eigenvalues

虛反矩陣

  • AFm×n,s=min{m,n}
    A=UΣVH

    定義
    A
    的pseudoinverse為
    A+=VΣ+UH

    其中
    Σ+Fn×m,(Σ+)ij=0,ij

    (Σ+)ij={1σi,  if  σi00,  if  σi=0

    又稱為Moore-Penrose廣義反矩陣

  • Σ+就是把
    Σ
    轉置之後把非零元素都變倒數

  • O+=O

統整算子

所有operator經過座標化都可以得到對應matrix,伴隨算子對應到

AH

名稱 定義 備註
Adjoint operator
<T(v),v>=<v,T(v)>
[T]β=[T]βH
Normal operator
TT=TT
AHA=AAH
self-adjoint operator
T=T
AH=A
也是normal operator,Hermitian matrix
skew self-adjoint operator
T=T,AH=A
也是normal operator,skew Hermitian matrix
unitary operator
C
中,
TT=I,AHA=I
orthogonal operator
R
中,
TT=I,ATA=I