112考資工研究所時寫的筆記,有錯歡迎留言指正
內文有引用他人筆記/教材/文章的地方
若有侵權十分抱歉,告知後將立刻撤除


1 矩陣與線性方程組

1-1 矩陣及矩陣運算

  1. 矩陣分割(滿簡單,略)

    • p1-20範例3
  2. 共軛轉置conjugate transpose

    • 定義:
      AT
      後,裡面element取共軛,
      aijaij
    • (AH)H=A
    • (AB)H=BHAH
    • (αA±βB)H=αAH±βBH
  3. 轉置transpose

    • 定義: 略
    • (AT)T=A
    • (AB)T=BTAT
    • (αA±βB)T=αAT±βBT
  4. 跡數trace

    • 定義: 對角線元素合
    • tr(αA±βB)=αtr(A)±βtr(B)
    • AFm×n,BFn×m,tr(AB)=tr(BA)
    • 證明: 用sigma把裡面元素操作一下就完成了
    • 其他:
      • tr(ATA)=0A=O

1-2 反矩陣

  • 反矩陣(inverse matrix)=非奇異矩陣(nonsingular matrix)
    • 只有
      n×n
      矩陣存在反矩陣
    • (AB)1=(B)1(A)1
    • (A1A2Ak)1=(Ak)1(Ak1)1(A1)1
      • 證明inverse=某串算式,不是要你求inverse,而是要你把它跟某人相乘得到
        I
      • I
        不是
        I
        ,而是
        AA1
  • n×n
    矩陣的inverse, 時間複雜度
    =O(n3)
  • 其他練習: 例17, 18, 精選範例3, 5

1-3 列基本運算

  1. elementary row operation

    • 對矩陣A的三種操作
    • 第一型列運算: A的i, j列交換,
      rij(A)
    • 第二型列運算: A的第i列乘k倍,k
      0,ri(k)(A)
    • 第三型列運算: A的第i列乘k倍加到第j列,
      rij(k)(A)
  2. row elementary matrix

    • 用法和例子看p1-35例22
  3. 其他練習: 精選範例4,5

1-4 線性方程組

  1. 擴增矩陣 Augmented matrix

    • [A|b]
      Ax=b
      的擴增矩陣
  2. Ax=b有解(unique, infinitely)稱作consistent,否則稱作inconsistent

  3. 齊次/非齊次 Homogeneous/nonhomogeneous

    • Ax=0
      為齊次線性系統,他的解叫做齊次解,且他必有trivial solution 0
    • Ax0
      為非齊次線性系統,他的解叫做非齊次解
  4. 列梯形型式 row-echelon form (ref)

    • 所有0列在非0列下方
    • 每一列最左邊的非0項愈下方愈靠右
    • 每一列最左邊的非0項為1
    • 每個矩陣都row equivalent某個唯一row echelon form
    • 用高斯消去法(Gaussian elimination)
  5. 簡化列梯形型式 reduced row echelon form (rref)

    • 為列梯形型式
    • 每一列最左邊非0項所在的行,其他元素為0
    • 每個矩陣都row equivalent某個唯一reduced row echelon form
    • 用Gauss-Jordan消去法(Gauss-Jordan elimination)
    • 考試若規定用某某消去法,一定要做到最後一步
  6. 這張有關rank的內容,直接看第四章

  7. 其他練習: 精選範例3,5

1-5 可逆矩陣充要條件

  • A,BFm×n,ABPFm×m使B=PA

  • Let

    AFn×n,下列敘述等價

    1. A is 可逆 (invertible)
    2. Ax=0
      只有0解
    3. xA=0
      只有0解
    4. A(行/列)等價於
      In
    5. A是若干(行/列)基本矩陣乘積
    6. ABFn×nA1=B
    7. ACFn×nA1=C
    8. bFn×1,Ax=b
      有解
    9. bF1×n,xA=b
      有解
    10. bFn×1,Ax=b
      有唯一解
    11. bF1×n,xA=b
      有唯一解
    12. det(A)0
    13. rank(A)=dim(RS(A))=dim(CS(A))=n
    14. adj(A)
      is invertible
    15. ker(A)=Lker(A)=0
    16. A為行/列獨立
  • 上面證明都滿經典的,最好都做一次

  • 其他題目: 精選範例3,4,5

1-6 LU分解

  • 定義

    AFm×n,LFm×m為下三角矩陣,
    UFm×n
    是樞元(pivot)未必為1的列梯形矩陣,
    A=LU

  • 步驟參考p1-76

  • LDU分解:

    • L, U跟上面一樣,
      DFm×m
      是對角矩陣
    • 步驟: 例p1-81例42
  • LU分解求

    Ax=b的解

    • let A=LU, so Ax=b>LUx=b
    • let Ux=y, then Ly=b
    • 解Ly=b的y
    • 解Ux=y的x
  • PTLU分解: 好像很少考,先不寫

  • Permutation matrix: 把

    I的行or列做任意排列,得到的矩陣

  • 補充: LU分解過程不能做列交換,因此並非每個矩陣皆可做LU分解

  • 其他題目: 精選範例1,2,3

1-7基本行運算

  • elementary column operation

    • 對矩陣A的三種操作
    • 第一型行運算: A的i, j行交換,
      cij(A)
    • 第二型行運算: A的第i行乘k倍,k
      0,ci(k)(A)
    • 第三型行運算: A的第i行乘k倍加到第j行,
      cij(k)(A)
  • column elementary matrix

    • 用法和例子看p1-91例46
  • 基本運算elementary operation指基本列運算 或 基本行運算

  • 基本矩陣elementary matrix指列基本矩陣 或 行基本矩陣

  • 其他練習: 精選範例1,2


2 行列式

2-1 二階行列式

  • 定義:

    A=(abcd)F2×2,定義A的行列式Determinant為
    adbc

  • 計做

    det(A)|A|

2-2 高階行列式

  1. 定義: 略,反正都會寫。只須注意

    AF1×1det(A)=a11

  2. 餘因子: 略,會求就好,他就是算古典伴隨矩陣(adj)或反矩陣時會用到的那個咚咚

  3. det(A)=det(AT)

  4. 其他練習: 精選範例2,4

2-3 行列式的性質

  1. Let

    AFn×n,n2, Aor,det(A)=0

  2. Let

    AFn×n

    • detrij(A)orcij(A), det(A)det(A)
    • det某行/列乘
      k,det(A)kdet(A)
    • 某行/列乘k加上另一行/列,det不變:
      det(A)det(A)
  3. A,BFn×n,det(AB)=det(A)det(B)

  4. det(αA)=αndet(A)

  5. skew-symmetric:

    AT=A

  6. If A is invertible,

    det(A1)=1det(A)

  7. Vandermonde matrix

  8. A,B,CFn×n,(ACOB)=(AOCB)=det(A)det(B)

  9. Upper/Lower triangular matrix's

    det(A)=i=1naii (對角線元素相乘)

  10. 其他練習: 精選範例 1,5,11

2-4 古典伴隨矩陣

  1. 略,會求就好

  2. AFn×n is invertible,
    A1=1det(A)adj(A)

  3. AFn×n is invertible,
    adj(A)1=Adet(A)

  4. A is invertible
    adj(A)
    is invertible

  5. 克拉瑪公式Cramer's rule (略)

  6. 其他練習: p2-46例26, 精選範例5,7

考古

  • 通靈題型: 2-59, 2-89, 2-90, 2-108
    • 還是有跡可循,觀察乘上的矩陣,都有類似的形狀(對角線為
      I
      的上or下三角矩陣

3 向量空間

3-1 向量空間

  • 假設V為非空集合,F為一個體,滿足十大公設,稱V為佈於F的向量空間(vector space),

    u,v,wV;α,βF

    1. 唯一存在
      u+vV
    2. 唯一存在
      αvV
    3. u+v=v+u
    4. (u+v)+w=u+(v+w)
    5. v+0=0+v=v
    6. vV,v+(v)=(v)+v=0
    7. α(u+v)=αu+αv
    8. (α+β)v=αv+βv
    9. (αβ)v=α(βv)
    10. 1v=v
  • 口訣:

    • 唯一存在(1,2)單位元(5,10)
    • 加法反(6)交換結合(3,4)
    • 乘分配交換結合(7,8,9)
  • 其他練習: 精選範例1,2,3

3-2 子空間

  1. 子空間(subspace)
    • 假設
      (V,+,)
      為一個佈於F的向量空間,若W為V的子集且
      (W,+,)
      仍為一個佈於F的向量空間,稱W為V的子空間
    • V為V的子空間
    • {0}稱為零空間(zero space),為任何向量空間的子空間
    • 證明某subset S是否為子空間(p3-15例5)
      1. 證明零矩陣(零空間)
        O
        在S中
      2. 已知
        A,BS
        ,證明
        αA+βB
        在S中
  2. 子空間充要條件
    • 假設
      (V,+,)
      為一個佈於F的向量空間,且
      WV,W
      ,則
    • (V,+,)
      (W,+,)
      的子空間
    • α,βF,u,vW,αu+βvW
  3. 子空間必要條件
    • 若W為V的子空間
    • 0W
    • If
      vWvW
  4. V是向量空間且
    W1,W2
    為的子空間,則
    W1W2V
  5. V是向量空間且
    W1,W2
    為的子空間,則
    • W1W2W1W2 or W2W1
  6. 和空間
    • 定義:
      Wi
      皆為V的子空間,定義
      W1,W2...Wk
      的和空間(sum space)
      W1+W2+...+Wk={w1+w2+...+wk|wiWi,i=1,2,...k}
    • W1,W2V
      的子空間,
      W1+W2
      為V的子空間
    • 四個基本子空間,
      AFn×m
      • CS(A)={Ax|xFn×1}
        , A的行空間
      • RS(A)={xA|xF1×m}
        , A的列空間
      • ker(A)={xFn×1,Ax=0}
        , A的核空間
      • Lker(A)={xF1×m,xA=0}
        , A的左核空間
    • A列等價於B
      • RS(A)=RS(B)
      • ker(A)=ker(B)
    • A行等價於B
      • CS(A)=CS(B)
      • Lker(A)=Lker(B)
  • 其他練習: p3-16例6, p3-17例7, p3-19定理3-4

3-3 生成與線性獨立

  1. 線性組合(略)
  2. 生成空間span
    • SV,span(S)={v|vS}
    • span(ϕ)={0}
  3. 線性獨立集
    • v1,v2,...vkV,if α1v1,α2v2,...αkvk=0
      α1=α2=...=αk=0
  4. 線性相依集
    • v1,v2,...vkV,
      存在有限個不全為0的純量
      α1,α2,...αkF
      使得
      α1v1,α2v2,...αkvk=0
  5. Wronskian 最好直接看課本
    • 算出來W後,求W(x)=det(W)
      • W(x)=0 x[a,b]
        linear dependent
      • W(0)0 x[a,b]
        linear independent
  • 其他練習: 例15,25,26,27, 精選範例4,8,9

3-4 基底與維度

  1. 基底basis
    • 定義: 滿足以下量條件,稱
      β
      為V的一組基底,
      β
      的元素個數稱為V的維度dimension
    • βVspan(β)=V
    • β
  2. S為V的基底
    S為最小生成集
    S為最大獨立集
  3. dim(W1+W2)=dim(W1)+dim(W2)dim(W1W2)
  • 其他練習: 例31,33,40,42 推廣3-4, 精選範例2,6,12,13

3-5 直和

  1. 獨立子空間
    • W1,W2,...WkV
      的子空間,若
      W1W2...Wk={0}
      ,稱
      W1,W2,...Wk
      為獨立子空間
  2. W1,W2,...Wk
    為獨立子空間
    • if w1+w2+...wk=0, then w1=w2=...wk=0,wiWi.
    • dim(W1+W2+...Wk)=dim(W1)+dim(W2)+...dim(Wk)
  3. 直和 direct sum
    • 滿足下列兩條件(
      W1,W2,...WkV
      )
    • V=W1+W2+...Wk
    • W1,W2,...Wk
      為獨立子空間
    • 記做
      V=W1W2...Wk

3-6 Lagrange內插法

偏少考,自己看課本


4 線性映射

4-1 線性映射

  1. 線性轉換/映射Linear transformation/mapping

    • u,vV,T(u+v)=T(u)+T(v)
    • αV,T(αv)=αT(v)
  2. 線性映射充要條件

    • α,βF,u,vV,T(αu+βv)=αT(u)+βT(v)
    • 證明是否為線性映射,用這個方法比較快
    • p4-6例5,
  3. 線性映射必要條件

    • T(0)=0
    • vV,T(v)=T(v)
    • u,vV,T(uv)=T(u)T(v)
    • 必要條件: 線性映射一定有這些性質,但有這些性質不一定線性映射!
    • p4-7例7, p4-22範例2
  4. 矩陣形狀變換
    以下都是乘在欲變換之矩陣的左邊!!!

    • 旋轉rotation
      沿著該軸逆時針轉
      θ
      度,沿著的軸為1
      ex. counterclockwise about z-axis
      (cosθsinθ0sinθcosθ0001)
    • 鏡射reflection
      沒提到的向量軸為-1,其他為1
      ex. reflection to/about yz-plane
      (100010001)
    • 投影projection
      沒提到的那個軸被銷掉,乘0,其他為1
      ex. project onto xy-plane
      (100010000)
    • 擴張/收縮dilation/Squeezing
      被收縮的向量乘k
      ex. dilation of factor 3
      (300030003)
    • 平移translation
      二維的平移需用三維矩陣,n維的平移需用n+1維矩陣。位移量在最右行
      ex. translation
      [2  3]T

      (102013001)
    • p4-23範例5, 範例10
  5. 其他練習: 精選範例1,3,9

4-2 座標化

  1. 座標向量

    • 定義: 把某向量轉到另一個向量空間後,用該向量空間的係數表示此向量
      x=[α1α2αn]
  2. 本節重點: 定理4-5

    • a,bF,u,vV,[au+bv]β=a[u]β+b[v]β
    • u=v[u]β=[v]β
    • xFn×1
      , 唯一存在
      vV
      使得
      [v]β=x
  3. 同構(isomorphic)

    • 定義: T: V> V'
      • T為線性映射
      • T為1-to-1 function
      • T為onto function
    • T稱作V至V'的線性同構函數,稱V與V'同構,
      VV
  4. (equivalence)

  5. β={u1,u2,un}γ={v1,v2,vn}的座標轉換矩陣(change of coordinate matrix)記做
    [I]βγ

  6. 定理4-8超重要,但自己看我懶得打

  • 其他練習: 例20,22,24,精選範例3,5

4-3 矩陣表示法與換底公式

  • 比較4-2, 4-3:
    • 4-2轉移矩陣(transition matrix)為
      [I]BA
      ,必定為
      n×n
      • 像形容詞,上例中形容B的單位矩陣,在A空間中長什麼樣子(A, B必同維)
    • 4-3相對於
      βγ
      的矩陣表示法(matrix representation of T relative to
      β
      and
      γ
      ),為
      [T]βγ
      ,可為
      n×m
      • 像動詞,把他乘在
        β
        空間的矩陣左邊,就是把原本
        β
        空間的向量變成
        γ
        空間的形狀
  1. 同空間的換底公式
    • [T]γγ=[I]βγ[T]ββ[I]γβ
    • 其中
      [T]γγ
      常被簡寫為
      [T]γ
  2. 不同空間的換底公式
    • [T]γγ=[I]βγ[T]ββ[I]γβ
  • 其他練習: 例28,30, 精選範例6,7,8,10

4-4 核空間與像集

  • 直像 direct image

    • T(S)={T(v)| vS}
  • 反像 inverse image

    • T1(S)={vV| T(v)S}
  • 核空間 kernel space = null space

    • ker(T)={vV| T(v)=0}
    • 可記做ker(T)或N(T)
  • range = image

    • R(T)=T(V)={T(v)| vV}
    • 可記做R(T), Im(T)或Range(T)
  • 核數 nullity

    • nullity=dim(ker(T))
  • 秩 rank

    • rank(T)=dim(R(T))
  • 維度定理

    • Let
      TL(V,V),
    • dim(ker(T))+dim(R(T))=nullity(T)+rank(T)=dim(V)
  • TL(V,V)

    • T 是一對一函數
      ker(T)=0
    • 若T是一對一函數
      dim(V)dim(V)
    • 若T是映成函數,
      dim(V)dim(V)
    • 若T是一對一且映成函數,
      dim(V)=dim(V)
  • A為一對一函數

    • N(T)=N(A)={0}
    • Ax=0
      只有零解 (或其他有關invertible的描述)
    • A為行獨立
    • A保獨立(某子集獨立,經過A轉換仍為獨立)
  • A為映成函數

    • R(T)=R(A)=Fm×1
    • 行生成
      Fm×1
    • A保生成(生成集,經過A轉換仍為生成集)
  • AFm×n

    • rank(A)+nullity(A)=n
    • rank(AT)+nullity(AT)=m
  • 面積,體積

    • A=[a1,a2], 為2*3矩陣,由a1, a2圍出來的平行四邊形面積為
      |det(A)|
    • A=[a1,a2,a3], 為3*3矩陣,由a1, a2, a3圍出來的平行六面體體積為
      |det(A)|
    • 原空間S上的面積/體積為r,經過A矩陣線性變換後
      • S為平行四邊形,T(S)面積為
        r|det(A)|
      • S為平行六面體,T(S)體積為
        r|det(A)|
  • 行獨立

    列獨立,行獨立
    列獨立

  • 行相依

    列相依,行相依
    列相依

  • 求行空間用行銷去;求列空間用列銷去,!!!用反會出事!!!

  • 寫題目時,需特別清楚以下答案寫法(表示法)

    • basis for ker/range of T
    • ker/range of T
    • 尤其是題目牽扯多項式時,很容易錯,需多寫題目孰悉
  • 其他練習: 例33,34,35, 精選範例1,4,6,7,13 (最好精選範例全寫)

4-5 矩陣的秩

  • rr(A)=dim(RS(A))cr(A)=dim(CS(A))
  • if
    rr(A)=cr(A)rank(A)=rr(A)=cr(A)
  • 維度定理(另一種表示法)
    • AFm×n
    • rank(A)+nullity(A)=n
  • rank(AB)min{rank(A),rank(B)}
  • rank(A+B)rank(A)+rank(B)
  • AFm×n
    • rank(A)=n
    • dim(RS(A))=n
    • dim(CS(A))=n
    • A
    • AF1×n
  • AFm×n
    • rank(A)=m
    • dim(RS(A))=m
    • dim(CS(A))=m
    • A
    • AFm×1

4-6 線性映射的合成與可逆

  • [T+U]βγ=[T]βγ+[U]βγ

  • [αT]βγ=α[T]βγ

  • [TU]βω=[T]γω[U]βγ

  • 都是早就知道的東西

  • 其他練習: 精選範例2

4-7 對偶空間與零化集

整張skip(資工系不用管)

5 對角化及其應用

5-1 相似性

  • A,BFn×n,PFn×n使P1AP=B,BAAB
    • A相似於B 跟 B相似於A,順序不是很重要,反正是iff的關係
  • 是一個等價關係
    • A~A
    • A~B -> B~A
    • A~B and A~C -> A~C
  • 換底公式、對角化有相似的關係
  • AB
    • tr(A)=tr(B)
    • det(A)=det(B)
    • rank(A)=rank(B)
    • nullity(A)=nullity(B)
  • 其他練習: 例2, 精選範例4

5-2 不變子空間

  • 不變子空間
    • TL(V,V)
      , 若
      wWT(w)W
      ,W為T-invariant subspace
  • TL(V,V)
    • {0}, V, ker(T), Im(T)皆為T不變子空間
  • 知道這樣就夠了,其他都看不太懂QQ(資工似乎少考)

5-3 特徵根及特徵向量

  • AFn×n,λF,若存在
    xFn×1,x0
    使得
    Ax=λx

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

  • 特徵多項式 characteristic polynomial

    • p(x)=det(AxI)=charA(x)
  • 不可先將A做列運算到上三角矩陣在計算特徵根,因為列運算不保證特徵根不變

  • If

    AB,

    • pA(x)=pB(x)det(AλI)=det(BλI)
    • A, B具有相同特徵根
    • !!!A, B未必具有相同特徵向量!!!
  • A,BFn×n,AB和BA有相同特徵根

  • A為可逆矩陣

    0不為A的特徵根

  • A為不可逆矩陣

    0是A的特徵根

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

  • If

    Ax=λx

    • A為可逆矩陣,則
      1λA1
      相對於x的特徵根
    • mN,λmAm
      相對於x的特徵根
  • 1+λ is an eigenvalue of
    I+A

    • (I+A)x=Ix+Ax=1x+Ax=1x+λx=(1+λ)x
    • A,BFn×n
      (定理5-13)
      • AB和BA有一樣的特徵根(特徵方程式)
      • A, B不為方陣時,AB, BA有相同非零特徵根
      • 應用:
        • 爆開
          ATA
          算特徵根>算到世界末日
        • 正解: 先求
          AAT
          特徵根,剩下的補0
          • (
            AAT
            最多只會求出兩個非零特徵根,但
            ATA
            有三個)
  • 其他練習: 例14, 17, 22, 23, 28, 精選範例4, 5, 10, 11

5-4 對角化

  • 可對角化 Diagonalizable

    • AFn×n
      ,若存在一可逆矩陣
      PFn×n
      使得
      P1AP=D
      為對角矩陣,稱A可對角化
  • Defective matrix: 不可對角化的矩陣(少於n個線性獨立特徵向量)

  • 代數重數 algebraic multiplicity

    • 特徵根
      λ
      的重根數,計做
      am(λ)
  • 幾何重數

    • dim(V(λ))=gm(λ)=nrank(AλI)
    • 詳細過程見p.5-60
  • gm(λ)am(λ)

  • If

    pA(x)pA(x)=(λ1x)(λ2x)...(λnx)

    • det(A)=λ1λ2...λn
      (含重根)
    • tr(A)=λ1+λ2+...+λn
      (含重根)
  • 例38,對角線為a,其他值皆為b的n*n矩陣

    • 超重要,直接背兩個特徵根和對應的特徵向量
    • 特徵根對應特徵向量:
      • λ1=ab
        ,特徵向量
        span{[1100][1010][1001]}
      • λ2==a+b(n1)
        ,特徵向量
        span{[1111]}
  • 可對角化充要條件

    • A有n個線性獨立特徵向量
    • pA(x)
      可分解且所有eigenvalue皆滿足
      gm(λ)=am(λ)
  • 可同步對角化 Simultaneously diagonalizable

    • If
      A,BFn×n
      ,A, B皆可對角化
      • A, B可同步對角化
        AB=BA
  • 其他練習: 例29, 30, 31, 35, 37, 39, 42, 45, 47, 範例3, 5, 7, 10, 11

5-5 冪等算子與矩陣

  • 冪等矩陣 idempotent matrix = 投影矩陣 projection matrix

    • AFn×n
      A2=A
  • 其他證明的部分沒有很懂,應該記定義就好了(吧

  • 其他練習: 精選範例4

5-6 對角化的應用

  • f(A)=Pf(D)P1

  • 二階矩陣的

    pA(x)=x2tr(A)x+det(A)

    • 推導: Let
      pA(x)=(xλ1)(xλ2)=x2(λ1+λ2)x+λ1λ2
  • 三階矩陣的

    pA(x)=x3tr(A)x2+(λ1λ2+λ2λ3+λ1λ3)xdet(A)
            
    =x3tr1(A)x2+tr2(A)xtr3(A)

    • where
      tr1(A)
      為一般的trace
    • tr2(A)
      : 對角線上任兩個元素,為2*2行列式的對角元素,把所有可能行列式相加
    • 圖為102台大數學
    • tr3(A)
      : 符合條件的3*3行列式的和,在上例中就是det(A),A為3*3
    • 推導類似二階
  • 題型:

    • 解微方(難!)
    • 解遞迴 (用離散的線性解比較快!)
    • 解高次矩陣
  • 一定要練過題目

  • 其他練習: 例57, 60, 62, 64

5-7 特徵根的近似解法

skip

5-8 Markov 鍊

  • 每行的機率必為1

  • 狀態向量合必為1

  • 所有

    λ1,且必存在
    λ=1

  • 穩定狀態向量: Px=x,可對角化且

    λ=1

    • 此時的
      v(1)=ker(λ1)=x
      即為穩定狀態
  • 定理5-42和下面的例79 一定要做

  • 難以言傳,務必寫過題目

  • 其他練習: 精選範例2, 6

6 Jordan型及其應用

6-1 冪零算子

  • 冪零算子/矩陣
    • 算子的部分是T,跳過
    • AFn×n
      , 若存在正整數k使得
      Ak=O
      ,則A稱為冪零矩陣,最小正整數k稱為A的指標index

6-2 循環子空間及循環分解

skip

6-3 Jordan型

  • p6-40前面放掉

  • 會用就好

  • 記得J從左上到右下,要按am(

    λ)來排序(大到小)

  • 其他練習: 例18, 精選範例4

6-4 Caylay Hamilton定理及其應用

  • 理論基礎薄弱,我也只會做題目,台大偶爾會考題10分的嚇人
  • 理論上都找的到不用這個定理的解法(爆破or觀察)
  • A1
    or 無法對角化的矩陣高次 很好用
  • 其他練習: 例22, 23, 29精選範例9, 10

6-5 Jordan型的應用

  • 用途: 無法對角化的矩陣A,求

    An

  • 課本寫的滿清楚的,我懶得打XD

  • 看到四階(含)以上的矩陣求

    An
    P
    ,不論是用對角化or用Jordan,需評估一下值不值得硬幹(不信的話挑戰精選範例4)

  • 其他練習: 例30, 31, 32

6-6 極小多項式

  • Minimal polynomial

    1. deg(f(x))1
      ,且f(x)最高項係數=1,稱作monic polynomial
    2. f(A)=O
    3. 對任意多項式g(x)滿足
      g(T)=O
      deg(f(x))deg(g(x))
  • Example

    • pA(x)=det(AxI)=(x2)(x3)3
    • mA(x){(x2)(x3),(x2)(x3)2,(x2)(x3)3}
    • 然後就從mA(x)中,次方最小的開始try,第一個為
      O
      的就是極小多項式
  • 長的很像Cayley-Hamilton的應用,極小多項式的次方通常比特徵多項式小,要慢慢try

  • 例39, 40, 41, 精選範例4, 10

7 內積空間

7-1 內積

  • 內積 Inner product

    • u,v,wV,<u+v,w> = <u,w>+<v,w>
    • u,vV,αF,<αu,v>=α<u,v>
    • u,vV,<u,v>=<v,u>
    • vV,v0,<v,v> > 0
      (大於0)
  • 長度 norm

    • 定義向量長度為
      ||v||=<v,v>
  • 柯西不等式 Cauchy-Schwarz inequality

    • |<u,v>|||u|| ||v||
  • 三角不等式 triangle inequality

    • ||u+v||||u||+||v||
  • 畢氏定理 pythagorean theorem

    • if
      uv,||u+v||2=||u||2+||v||2
  • 單範正交集 Orthonormal set

    • u,vS,uv,||u||=||v||=1,<u,v>=0
  • 其他練習: 例1, 3, 5, 10, 15, 精選範例2, 4, 6, 12

7-2 Gram-schmidt正交化及QR分解

  • S={v1,v2,v3...}
    為內積空間,且
    vS
    ,可以把v表示成S中向量的線性組合
    • v=<v,v1><v1,v1>v1+<v,v2><v2,v2>v2+<v,v3><v3,v3>v3 + ...
  • Gram-schmidt 正交化過程
    • 直接看例24
  • QR分解
    • 直接看例32
  • 題型沒啥變化,但極容易計算錯誤,務必小心寫
  • 其他練習: 例26, 28, 30, 31, 34

7-3 正交投影

  • V為內積空間,W為V的子空間,

    β={v1,v2,...vk}為W的一組正交基底
    vV

    • 正交向量唯一存在,記做
      projwv
    • projwv=i=1k<v,vi><vi,vi>vi
      (正交投影公式)
  • Hermitian matrix

    • AH
      : 對實係數矩陣來說就是
      AT
    • ker(AHA)=ker(A)
    • rank(AHA)=rank(A)
  • AFm×n行獨立,
    W=R(A),bFm×1
    ,則

    • projwb=A(AHA)1AHb
    • 其中
      A(AHA)1AH
      被稱作正交投影向量orthogonal projection matrix
  • xFn×1使
    ||Axb||
    最小
    AHAx=AHb
    (normal equation)

    • normal equation
      AHAx=AHb
      必定有解
    • A為行獨立,唯一解;否則x有無限多解
    • x稱為最小平方解 (Least square solution)
  • 比較兩個投影公式

    • projwv=i=1k<v,vi><vi,vi>vi
      • 可用時機: W需為正交基底
    • projAb=A(AHA)1AHb
      • 可用時機: A為行獨立,且只能用於歐氏空間和標準內積
  • 最小平方直線

    • 見例51
  • A=QR為A的QR分解

    • xFn×1
      使
      ||Axb||
      最小
      Rx=QHb
  • 其他練習: 例35, 37, 38, 39, 42, 46, 49, 精選範例16

7-4 正交補空間

  • SV,S={vV|<v,s>=0,sS}

    • 不懂看例56
  • V為內積空間,W為V的子空間

    • dim(V)=dim(W)+dim(W)
    • W⊥⊥=W
  • V為內積空間,W為V的子空間,P為V在W上的正交投影算子,則

    • vP(v)vW
      上的正交投影向量,即
      projWv=vprojWv
    • IP
      為V在
      W
      上的正交投影算子
      • projWv=vprojWv=vPv=(IP)v
  • R(AT)=N(A) ->可以想成T翻進去變倒的,N變R,R變N

    • 以下三個也成立
    • R(AT)=N(A)
    • R(A)=N(AT)
    • N(AT)=R(A)
  • 極小解 Minimal solution

    • Ax=b
      的極小解,令為s,則對其他滿足
      Ax=b
      的解u,
      ||s||||u||
    • 先解
      AATx=b
    • 取x為某一解,則
      ATx
      為極小解
  • 點到面距離(例64)

    • 另解:
      P(x0,y0,z0)E:ax+by+cz+d=0
      的距離
    • d=|ax0+by0+cz0+d|a2+b2+c2
  • 面到面距離(例65)

    • 另解:
      E1:ax+by+cz+d1E2:ax+by+cz+d2=0
      的距離
    • d=|d1d2|a2+b2+c2
  • 這張非常抽象,一定要練習

  • 其他練習: 例56, 63, 64, 65, 68, 69

投影公式整理

  • projAb=A(AHA)1AHb
  • 最小平方解
    • (
      AHAx=AHb
      )的解
    • (
      Ax=projAb
      )的解
  • 極小解
    • 先解
      AATu=b
    • 再解
      s=ATu

8 內積上的算子及應用

8-1 伴隨算子

  • 伴隨算子
    T
    adjoint operator
    • TL(V,V),
      若存在
      T:VV
      為線性映射使得
      <T(u),v>=<u,T(v)>,u,vV
    • 唯一存在
    • β={v1,v2...vn}
      為V的一組單範正交(orthonormal)基底
    • [T]β=[T]βH
      ,這邊H是矩陣取Hermitian matrix
  • A commute with B
    • AB=BA
  • 其他練習: 精選範例2,3

8-2 正規算子與矩陣

  • 正規算子 normal operator

    • TT=TT
  • 正規矩陣 normal matrix

    • AHA=AAH
    • λ1,λ2
      是兩個相異特徵根,且v1, v2為對應特徵向量,則
      v1v2
  • 自伴算子 self-adjoint operator

    • T=T
  • 斜自伴算子 skew self-adjoint operator

    • T=T
  • Conjugate transpose

    • AH=(A)T=AT
    • A的元素階為實數
      AH=AT
  • Hermitian matrix

    • AH=A
    • A為正規矩陣
    • A的特徵根階為實數
    • A的相異特徵根對應的特徵向量必正交
    • A的主對角線元素階為實數
  • skew-Hermitian matrix

    • AH=A
    • A為正規矩陣
    • A的特徵根階為0或虛數
    • A的相異特徵根對應的特徵向量必正交
    • A的主對角線元素階為0或虛數
  • symmetric matrix

    • AT=A
    • 適用hermitian matrix 的所有性質
  • skew-symmetric matrix

    • AT=A
    • 適用skew-hermitian matrix 的所有性質
  • 么正矩陣 unitary matrix

    • ACn×n,AHA=I
    • A為正規矩陣
    • A的特徵根滿足
      |λ|=1
    • A的相異特徵根對應的特徵向量必為正交
    • |det(A)|=1
  • 正交矩陣 orthogonal matrix

    • ARn×n,ATA=I
    • A為正規矩陣
    • A的特徵根滿足
      λ=±1
    • A的相異特徵根對應的特徵向量必為正交
    • det(A)=±1
  • ARm×n

    • AATATA
      都必為symmetric,而symmetric可正交對角化,因此
      AATATA
      都可正交對角化
    • AATATA
      皆必為半正定
  • 其他練習: 精選範例5

8-3 么正及正交算子的特性

  • 好像只有最後兩頁是資工會的

  • ACn×n(ARn×n),下列敘述等價

    • A為么正矩陣 (A為正交矩陣)
    • AHA=I
      ,
      (ATA=I)
    • A的行向量形成單範正交集
    • A的列向量形成單範正交集
  • 其他練習: 例12, 精選範例2,3,4

8-4 雙線性式與半雙線性式

skip

8-5 正定及正半定算子與矩陣

  • 也是有看沒有懂

  • 二次式 quadratic form

    • Q(x)=xHAx
  • 假設

    AFn×n,Q(x)=xHAx為A的二次式

    • 正定 (positive definite form):
      Q(x)>0,x0
    • 正伴定 (positive semidefinite form):
      Q(x)0,x
    • 負定 (negative definite form):
      Q(x)<0,x0
    • 負半定 (negative semidefinite form):
      Q(x)0,x
    • 不定型 (indefinite form):
      Q(x)
      可能的值有正也有負
  • 其他練習: 例19

8-6 么正及正交對角化

  • A,BCn×n,PHAP=B,PCn×n
    為么正矩陣
    • 稱為A與B么正相似 (unitarily similar)
    • 或稱A與B么正等價 (unitarily equivalent)
  • A,BRn×n,PHAP=B,PRn×n
    為正交矩陣
    • 稱為A與B正交相似 (orthogonally similar)
    • 或稱A與B正交等價 (orthogonally equivalent)
  • Schur's 定理
    • 假設
      ACn×n(Rn×n),
      存在么正矩陣(正交矩陣)
      PCn×n(Rn×n)使PTAP
      為上三角矩陣
  • A可么正對角化
    AAT=ATA
    (A為正規矩陣)
  • A可正交對角化
    A=AT
    (A為對稱矩陣)
  • 題型
    • Unitary (Orthogonal) matrix P s.t
      D=PTAP(D=P1AP)
      • 例26, 27, 精選範例2~7,十題裡面八題都這樣考!!!
    • 光譜分解 spectral decomposition
      • 例29, 30

8-7 正定及正半定矩陣的特性

  • 正定 positive definite
  • Δk(A)
    是A的主子行列式 Principal minors
    • 從左上角開始的正方形小矩陣取det(A'), A'從
      1×1,2×2...n×n
    • 自己看課本p8-83
  • ARn×n
    為對稱矩陣
    • A為正定
      A所有特徵根皆為正
      Δk(A)>0,k=1,2,...n
    • 若A不為對稱矩陣,不能直接從特徵根或
      Δk(A)
      判斷
    • Q(x)=xHAx=xHBx
      找到對稱的B,再從B的特徵根或
      Δk(A)
      判斷
    • 參見p8-51注意事項8-18
    • 例32,33,34, 精選範例2
  • ACn×n
    且A為Hermitian矩陣
    • A為正定(positive definite)
      A的所有特徵根皆為正
    • A為半正定(positive semidefinite)
      A的所有特徵根皆為非負
    • A為負定(negative definite)
      A的所有特徵根皆為負
    • A為半負定(negative semidefinite)
      A的所有特徵根皆為非正
    • 精選範例3,4
  • Cholesky分解
    • 分解
      A=LLT
      ,其中A為正定矩陣,L為下三角矩陣且對角元素皆為正
    • 兩種方法:
      1. A作LU分解後推得(例36)
      2. A作正交對角化,得
        A=PDPH
        後推得(例38)
    • 若A為正半定,能用法2算出
      A=BBH
      (例39)
    • LU分解做比較
      • LU分解的對角線沒有限制,Cholesky的限制對角線元素要為正
  • 其他練習: 精選範例5,6,7

8-8 二次式的應用

  • 兩種題型

    • 題型一: 二次式座標轉換
      • q(x,y)=xTAxxTPDPTx=(PTx)TD(PTx)=yTDy
      • 會發現其中y是轉換過的新座標(x',y'),且
        y=PTx
        , 其中
        PT
        是旋轉矩陣
      • PT=(cosθsinθsinθcosθ)
      • 一定要看例41, 42, 43
    • 題型二: 二次式的最大、最小值 (Rayleigh principle)
  • Rayleigh principle

    • 最大值為最大的eigenvalue,發生在對應的eigenvector
    • 最小值為最小的eigenvalue,發生在對應的eigenvector
    • 例45, 46, 47, 精選範例5, 12
  • 其他練習: 例41, 42

8-9 矩陣的長度及條件數

  • 矩陣長度

    • 滿足條件的函數都可稱該運算為矩陣長度 matrix norm
    • 4條件:
    • ||A||0||A||=0A=O
    • ||αA||=|α| ||A||
    • ||A+B||||A||+||B||
      三角不等式
    • ||AC||||A|| ||C||
      柯西不等式
  • Frobenius長度

    • 矩陣中每個元素的平方合開根號, 寫作
      ||A||F
    • ex.
      A=(1243)
      ,
      ||A||F=12+22+42+(3)2=30
  • 算子長度 operator norm

    • 1-norm
      ||A||1
      • 每行元素取絕對值的合,最大的即為
        ||A||1
    • 2-norm
      ||A||2
      • ATA
        的最大eigenvalue
        λmax
        開根號
        =λmax(ATA)
    • -norm
      ||A||
      • 每列元素取絕對值的合,最大的即為
        ||A||
  • 條件數 cond(A)

    • cond1(A)=||A||1||A1||1
    • cond2(A)=||A||2||A1||2=|λmax||λmin|
    • cond(A)=||A||||A1||
  • 其他練習: 例51, 53, 精選範例2

8-10 Householder轉換

  • Householder 矩陣

    • uRn×1, ||u||=1,H=I2uuT
    • H為對稱矩陣, 即
      HT=H
    • H為正交矩陣, 即
      HT=H1
    • 因此
      H=HT=H1
    • H只存在兩個特徵根:
      1, 1
    • gm(1)=1
    • gm(1)=n1
    • 因此
      tr(H)=sum of λ=11+1(n1)=n2
    • H可逆, 所以
      rank(H)=n
  • 單位法向量

    u

    • 直線
      ax+by=0 ()
      • u=1a2+b2[ab]
    • 平面
      ax+by+cz=0 ()
      • u=1a2+b2+c2[abc]
  • 定理8-51:

    x,yRn×1,xy,||x||=||y||,u=xy||xy||

    • H=I2uuT
      為相對u的householder矩陣且
      Hx=y
  • 解題目時盡量湊出

    uTu, 因為他運算後是常數 (when
    uRn×1
    )

  • 其他練習: 精選範例5, 6, 7

8-11 奇異值分解

  • Singular Value Decomposition (SVD)

    • A=UΣVT
    • ACm×n, UCm×m, VCn×n, ΣCm×n
    • σ1σ2...σs, s=min{m,n}
    • SVD一定存在,但U, V不唯一
    • rank(A)=r=A
      的非零signlar value個數
      • v1,...vrR(AH)
        的單範正交基底(orthonormal basis)
      • vr+1,...vnN(A)
        的單範正交基底(orthonormal basis)
      • u1,...urR(A)
        的單範正交基底(orthonormal basis)
      • ur+1,...umN(AH)
        的單範正交基底(orthonormal basis)
      • 很重要!!!
  • 虛反矩陣 pseudoinverse

    • 先求A的SVD:
      A=UΣVT
    • 再求
      Σ+
      : 先取
      ΣT
      ,再把所有
      σi
      1σi
    • 虛反矩陣
      A+=VΣ+UT
    • 若A行獨立
      rank(A)=n(AHA)1
      存在
      • A+=(AHA)1AH
        ,這樣算快多了
  • 這張的其他練習都一定要鍊

  • 其他練習: 例63, 64, 66, 68 精選範例4