[筆記][總結] MIT線性代數 Gilbert Strang 矩陣運算

作者水平有限,歡迎大家提出文中錯誤


本文主要討論除加法和數乘之外的矩陣常用運算

轉置

Aij=AjiT A_{ij}=A^T_{ji}

轉置的性質

  1. (AT)1(A^T)^{-1}=(A1)T(A^{-1})^T
  2. 對於任意矩陣RRRTRR^TR是對稱矩陣

求逆

使用高斯消元法求解逆矩陣

通過高斯消元法,可以將可逆矩陣化爲單位矩陣,但是此時左乘的矩陣不僅僅是消元矩陣,而是涉及到所有初等行變換。如果將合作用依然寫作EE
那麼有
EA=I EA=I

E=A1 E=A^{-1}
如果同時對矩陣AA和單位陣II做初等行變換,則可以構造增廣矩陣[AI][A|I]
E[AI]=[IA1] E[A|I]=[I|A^{-1}]

矩陣逆的代數表達式

這裏會用到本部分後面的內容
A1=1det ACT A^{-1}=\frac{1}{det\ A}C^T
矩陣CC被稱爲代數餘子式矩陣
C=[A11A12A1nA21A22A2nAn1An2Ann] C= \left[ \begin{matrix} A_{11}&A_{12}&\cdots &A_{1n}\\ A_{21}&A_{22}&\cdots &A_{2n}\\ \vdots&\vdots&\ddots &\vdots\\ A_{n1}&A_{n2}&\cdots &A_{nn}\\ \end{matrix} \right]
CTC^T稱爲伴隨矩陣

證明:
ACT=[a11a12a1na21a22a2nan1an2ann][A11A21An1A12A22An2A1nA2nAnn] AC^T= \left[ \begin{matrix} a_{11}&a_{12}&\cdots &a_{1n}\\ a_{21}&a_{22}&\cdots &a_{2n}\\ \vdots&\vdots&\ddots &\vdots\\ a_{n1}&a_{n2}&\cdots &a_{nn}\\ \end{matrix} \right] \left[ \begin{matrix} A_{11}&A_{21}&\cdots &A_{n1}\\ A_{12}&A_{22}&\cdots &A_{n2}\\ \vdots&\vdots&\ddots &\vdots\\ A_{1n}&A_{2n}&\cdots &A_{nn}\\ \end{matrix} \right]
得到
(ACT)ij=ai1Aj1+ai2Aj2++ainAjn (AC^T)_{ij}=a_{i1}A_{j1}+a_{i2}A_{j2}+\cdots+a_{in}A_{jn}
注意到i=ji=j時,(ACT)ii=det A(AC^T)_{ii}=det\ A
對於其他項,考慮其意義,是矩陣AA某一行的元素,乘以一個不同行對應的代數餘子式,如果將矩陣的第jj行用第ii行進行代替,其行列式的值正是上式
a11a12a1nai1ai2ainaj1aj2ajnan1an2anna11a12a1nai1ai2ainai1ai2ainan1an2ann] \left| \begin{matrix} a_{11}&a_{12}&&\cdots &&a_{1n}\\ \vdots&\vdots&&&&\vdots\\ a_{i1}&a_{i2}&&\cdots &&a_{in}\\ \vdots&\vdots&&&&\vdots\\ a_{j1}&a_{j2}&&\cdots &&a_{jn}\\ \vdots&\vdots&&&&\vdots\\ a_{n1}&a_{n2}&&\cdots &&a_{nn}\\ \end{matrix} \right| \rightarrow \left| \begin{matrix} a_{11}&a_{12}&&\cdots &&a_{1n}\\ \vdots&\vdots&&&&\vdots\\ a_{i1}&a_{i2}&&\cdots &&a_{in}\\ \vdots&\vdots&&&&\vdots\\ a_{i1}&a_{i2}&&\cdots &&a_{in}\\ \vdots&\vdots&&&&\vdots\\ a_{n1}&a_{n2}&&\cdots &&a_{nn}\\ \end{matrix} \right]
所以iji\not= j時,(ACT)ij=0(AC^T)_{ij}=0
所以ACT=det AIAC^T=det\ A\cdot I,證畢

這樣的公式對實際中矩陣逆的求解沒有什麼價值,因爲時間複雜度太高,但是它提供了逆矩陣的代數表達式

Cramers ruleCramer's\ rule

上一小節知道了矩陣逆的代數表達式,現在考慮矩陣方程Ax=bAx=b,可知
x=1det ACTb x=\frac{1}{det\ A}C^Tb
考慮xx的其中一項
det Axi=b1C1i+b2C2i++bnCni det\ A\cdot x_i=b_1C_{1i}+b_2C_{2i}+\cdots+b_nC_{ni}
構造矩陣BiB_i
Bi=[a11b1a1na21b2a2nan1bnann] B_i= \left[ \begin{matrix} a_{11}&\cdots &b_1&\cdots&a_{1n}\\ a_{21}&\cdots&b_2&\cdots&a_{2n}\\ \vdots&&\vdots&&\vdots\\ a_{n1}&\cdots&b_n &\cdots&a_{nn}\\ \end{matrix} \right]
矩陣AA用列向量bb替換了第ii列,形成了矩陣BiB_i,所以
det Axi=det Bi det\ A\cdot x_i=det\ B_i
克拉默法則
矩陣方程Ax=bAx=b若是有解且唯一的,那麼其解爲
xi=det Bidet A x_i=\frac{det\ B_i}{det\ A}

方陣逆矩陣的性質

如果AA是可逆矩陣,AA1=A1A=IAA^{-1}=A^{-1}A=I

矩陣乘法

矩陣乘法的五種形式

1. 行列內積

Cmn=AmpBpn C_{m*n}=A_{m*p}B_{p*n}

cij=row i(A)column j(A) c_{ij}=row\ i(A)\cdot column\ j(A)

2. 矩陣與列向量之積

column i(C)=A column i(B) column\ i(C)=A\ column\ i(B)

3. 行向量和矩陣之積

row i(C)=row iAB row\ i(C)=row\ iA\cdot B

4. 若干秩1矩陣之和

C=i=1p[a1ia2iami][bi1bi2bin] C=\sum\limits_{i=1}^p \left[ \begin{matrix} a_{1i}\\ a_{2i}\\ \vdots\\ a_{mi}\\ \end{matrix} \right] \left[ \begin{matrix} b_{i1}&b_{i2}&\cdots &b_{in} \end{matrix} \right]
比如
[273849][1600]=[234][16]3.[789][00] \left[ \begin{matrix} 2&7\\ 3&8\\ 4&9 \end{matrix} \right] \left[ \begin{matrix} 1&6\\ 0&0\\ \end{matrix} \right]= \left[ \begin{matrix} 2\\ 3\\ 4 \end{matrix} \right] \left[ \begin{matrix} 1&6 \end{matrix} \right] 3. \left[ \begin{matrix} 7\\ 8\\ 9 \end{matrix} \right] \left[ \begin{matrix} 0&0\\ \end{matrix} \right]

5. 分塊乘法

[A11A12A21A22][B11B12B21B22]=[A11B11+A12B21A11B12+A12B22A21B11+A22B21A21B12+A22B22] \left[ \begin{matrix} A_{11}&A_{12}\\ A_{21}&A_{22}\\ \end{matrix} \right] \left[ \begin{matrix} B_{11}&B_{12}\\ B_{21}&B_{22}\\ \end{matrix} \right]= \left[ \begin{matrix} A_{11}B_{11}+A_{12}B_{21}&A_{11}B_{12}+A_{12}B_{22}\\ A_{21}B_{11}+A_{22}B_{21}&A_{21}B_{12}+A_{22}B_{22}\\ \end{matrix} \right]
分塊必須匹配

矩陣乘法的性質

(AB)1=B1A1(AB)^{-1}=B^{-1}A^{-1}

行列式

行列式的性質

1.det I=1det\ I=1
2.行交換或列交換,會使行列式值取反

結合性質1,2可知det Pij=1det\ P_{ij}=-1(單次交換)
3.
tatbcd=tabcd \left| \begin{matrix} ta&tb\\ c&d\\ \end{matrix} \right|=t \left| \begin{matrix} a&b\\ c&d\\ \end{matrix} \right|
det kA=kndet A det\ kA=k^ndet\ A
4.

a+ab+bcd=abcd+abcd \left| \begin{matrix} a+a'&b+b'\\ c&d\\ \end{matrix} \right|= \left| \begin{matrix} a&b\\ c&d\\ \end{matrix} \right|+ \left| \begin{matrix} a'&b'\\ c&d\\ \end{matrix} \right|

結合性質3,4可知det det\ \cdot對矩陣某一行或一列的線性算子
5.矩陣中出現兩個相等的行或列,會使行列式爲0
6.消去操作(從某一行加上零一行的k倍)不影響行列式的值
a+kcb+kdcd=abcd+kckccd=abcd \left| \begin{matrix} a+kc&b+kd\\ c&d\\ \end{matrix} \right|= \left| \begin{matrix} a&b\\ c&d\\ \end{matrix} \right|+ \left| \begin{matrix} kc&kc\\ c&d\\ \end{matrix} \right|= \left| \begin{matrix} a&b\\ c&d\\ \end{matrix} \right|
7.全零行(列)會導致行列式爲零
8.三角陣的行列式等於對角元之積
9.奇異矩陣行列式爲0
10.det AB=det Adet Bdet\ AB=det\ A\cdot det\ B
11.det A1=(det A)1det\ A^{-1}=(det\ A)^{-1}
12.det A=det ATdet\ A=det\ A^T

行列式的計算方法

行列式的定義式

接下來,利用上面的性質1,2,3,可以得到行列式的計算方法
以二階方陣爲例
abcd=a0cd+0bcd= \left| \begin{matrix} a&b\\ c&d\\ \end{matrix} \right|= \left| \begin{matrix} a&0\\ c&d\\ \end{matrix} \right|+ \left| \begin{matrix} 0&b\\ c&d\\ \end{matrix} \right|=

a0c0+a00d+0bc0+0b0d=adbc \left| \begin{matrix} a&0\\ c&0\\ \end{matrix} \right|+ \left| \begin{matrix} a&0\\ 0&d\\ \end{matrix} \right|+ \left| \begin{matrix} 0&b\\ c&0\\ \end{matrix} \right|+ \left| \begin{matrix} 0&b\\ 0&d\\ \end{matrix} \right|=ad-bc
可見,對於n階方陣,依次分解第一行,第二行……第nn行,共得到nnn^n項,其中只要有某項的一列或一行全爲零,則此項爲0,所有非零項的種類數問題,變成了一個排列組合問題,相當於對1n1\cdots n進行全排列,項總數爲n!n!,留下的項不僅沒有空行空列,而且每行或每列都只有一個元素,都可以使用多次置換矩陣PP左乘對角陣得到。

而這n!n!個項的值又是什麼,每一項都等於置換矩陣和某對角陣的乘積,這個對角陣的行列式就是這個項的所有非零元素之積

對於符號,已經知道,一次置換就會使行列式符號相反,這個總的置換矩陣PP,假設非零項分別爲p1a1p2a2...pnanp_{1a_1}p_{2a_2}...p_{na_n}{a1 a2 an}\{a_1\ a_2\cdots\ a_n\}1,2n1,2\cdots n的一種排列

所以,置換矩陣PP的行列式值是+1+1還是1-1,就取決於{1 2n}\{1\ 2\cdots n\}要經過多少次置換才能達到{a1 a2 an}\{a_1\ a_2\cdots\ a_n\},是偶數次還是奇數次。

考慮這個問題,需要先介紹逆序數的概念,在序列{a1 a2 an}\{a_1\ a_2\cdots\ a_n\}中,所有的二元對(ai,aj)(a_i,a_j)都有兩種情況,假設i<ji<j,一種是ai>aja_i>a_j,另一種是ai<aja_i<a_j(注意ai=aja_i=a_j是不可能的,因爲{a1 a2 an}\{a_1\ a_2\cdots\ a_n\}1,2n1,2\cdots n的一種排列),而逆序數就是所有ai>aja_i>a_j情況的總數,記作τ(a1a2...an)\tau(a_1a_2...a_n),顯然τ(1 2 n)=0\tau(1\ 2\ \cdots n)=0

置換操作會如何影響逆序數?
{a1 a2aiaj an}{a1 a2ajai an} \{a_1\ a_2\cdots a_i\cdots a_j\cdots\ a_n\}\rightarrow \{a_1\ a_2\cdots a_j\cdots a_i\cdots\ a_n\}

這裏要證明一次置換會導致逆序數的奇偶性發生改變,爲了避免對序列其他部分的影響,通過相鄰元素的置換完成置換
{a1 a2ai ai+1aj an}{a1 a2ai+1 aiaj an}} \{a_1\ a_2\cdots a_i\ a_{i+1}\cdots a_j\cdots\ a_n\}\rightarrow \{a_1\ a_2\cdots a_{i+1}\ a_i\cdots a_j\cdots\ a_n\}\}
這樣的一次置換使得序列的逆序數只加11或減11,奇偶性一定會改變
假設原序列中aia_iaja_j之間有nn個元素
{a1 a2aiaj an}{a1 a2ai+1aj ai an} \{a_1\ a_2\cdots a_i\cdots a_j\cdots\ a_n\}\rightarrow \{a_1\ a_2\cdots a_{i+1}\cdots a_j\ a_i\cdots\ a_n\}

需要n+1n+1

{a1 a2ai+1aj ai an}{a1 a2aj ai an} \{a_1\ a_2\cdots a_{i+1}\cdots a_j\ a_i\cdots\ a_n\}\rightarrow \{a_1\ a_2\cdots a_j\cdots \ a_i\cdots\ a_n\}

需要nn
總步數是2n+12n+1則不管置換髮生在什麼位置,逆序數奇偶性一定會改變

那麼逆序數的奇偶性就能代表置換次數的奇偶性,這正是所求的,τ(1 2...n)=0\tau(1\ 2...n)=0,之後的每次置換都會同步改變det Pdet\ P的極性,以及序列逆序數的奇偶性。至此置換矩陣PP的行列式就已經清晰了
det P=(1)τ(p1p2...pn) det\ P=(-1)^{\tau(p_1p_2...p_n)}
再次說明,p1p2...pnp_1p_2...p_n是置換矩陣每行非零元所在列數(或每列非零元所在行數)

最後有行列式的定義式
det(Ann)=all of permutations(1)τ(p1p2...pn)a1p1a2p2...anpn det(A_{n*n})=\sum\limits_{all\ of\ permutations}(-1)^{\tau(p_1p_2...p_n)}a_{1p_1}a_{2p_2}...a_{np_n}

τ(p1p2...pn) is the inversion number of permutation p1p2...pn. {\tau(p_1p_2...p_n)}\ is\ the\ inversion\ number\ of\ permutation\ p_1p_2...p_n.

代數餘子式法

在剛纔的定義式推導的過程中,所有n!n!項其實可以進行如下劃分
{n! terms}=\{n!\ terms\}=
{(n1)! terms with nonzero a11}+\{(n-1)!\ terms\ with\ nonzero\ a_{11}\}+
{(n1)! terms with nonzero a12}+\{(n-1)!\ terms\ with\ nonzero\ a_{12}\}+
\vdots
{(n1)! terms with nonzero a1n}\{(n-1)!\ terms\ with\ nonzero\ a_{1n}\}
而這每一個劃分之和,又好像滿足一個n1n-1階行列式的定義式,比如
a1100000a1200000a1300 \left| \begin{matrix} a_{11}&0&0\\ 0&*&*\\ 0&*&*\\ \end{matrix} \right| \left| \begin{matrix} 0&a_{12}&0\\ *&0&*\\ *&0&*\\ \end{matrix} \right| \left| \begin{matrix} 0&0&a_{13}\\ *&*&0\\ *&*&0\\ \end{matrix} \right|
*是可以非零的位置

所以這個nn階的行列式,可以分裂成nnn1n-1類似行列式的量之和。
注意這裏的劃分不是唯一的,可以對第一行展開,也可以對任意行展開,甚至可以對任意列展開,不過是何種展開,都能得到這個低一階的類似行列式的量
這個類似行列式的量被稱爲代數餘子式,記作AijA_{ij}
det(Ann)=j=1naijAij,(i=1,2,...,n)det(A_{n*n}) = \sum\limits_{j=1}^na_{ij}A_{ij},(i=1,2,...,n)

det(Ann)=i=1naijAij,(j=1,2,...,n)det(A_{n*n}) = \sum\limits_{i=1}^na_{ij}A_{ij},(j=1,2,...,n)

代數餘子式

爲什麼說代數餘子式是類似行列式的量,它和行列式差在哪裏?
333*3矩陣爲例
a11a12a13a21a22a23a31a32a33= \left| \begin{matrix} a_{11}&a_{12}&a_{13}\\ a_{21}&a_{22}&a_{23}\\ a_{31}&a_{32}&a_{33}\\ \end{matrix} \right|=

a11(a22a33a23a32)a12(a21a33a23a31)+a13(a21a32a22a31) a_{11}(a_{22}a_{33}-a_{23}a_{32})-a_{12}(a_{21}a_{33}-a_{23}a_{31})+a_{13}(a_{21}a_{32}-a_{22}a_{31})
從中可以觀察出規律aija_{ij}的代數餘子式Aij=(1)i+jMijA_{ij}=(-1)^{i+j}M_{ij},其中MijM_{ij}是除去iijj列的元素,組成的子行列式,被稱爲aija_{ij}餘子式

所以333*3矩陣的代數餘子式的極性爲
+++++ \left| \begin{matrix} +&-&+\\ -&+&-\\ +&-&+\\ \end{matrix} \right|

上三角陣對角線元素之積

以上兩種方法的時間複雜度都是O(n!)O(n!),其實根據行列式的倍加不變性質,可以知道det A=det Udet\ A=det\ U,而高斯消元法的時間複雜度只有O(n3)O(n^3),所以MATLAB在求解矩陣行列式時,都是先消元然後計算主元之積。
det(Ann)=i=1npivot elementidet(A_{n*n}) = \prod\limits_{i=1}^npivot\ element_i

例子

三對角矩陣AnA_n,比如
A4=[1111111111] A_4= \left[ \begin{matrix} 1&1&&\\ 1&1&1&\\ &1&1&1\\ &&1&1\\ \end{matrix} \right]
求一下det Andet\ A_n通項
可以直接看出A1=1,A2=0,A3=1A_1=1,A_2=0,A_3=-1,從A4A_4開始
det A4=1111111111= det\ A_4= \left| \begin{matrix} 1&1&&\\ 1&1&1&\\ &1&1&1\\ &&1&1\\ \end{matrix} \right|=
11111111+(1)111111=det A3det A2=1 1\cdot \left| \begin{matrix} 1&1&\\ 1&1&1\\ &1&1\\ \end{matrix} \right|+(-1)\cdot \left| \begin{matrix} 1&1&\\ &1&1\\ &1&1\\ \end{matrix} \right|=det\ A_3-det\ A_2=-1
不難發現det An=det An1det An2det\ A_n=det\ A_{n-1}-det\ A_{n-2}
所以{An}\{A_n\}是一個循環的數列,循環部分爲
{1 0 1 1 0 1} \{1\ 0\ -1\ -1\ 0\ 1\}

特徵值和特徵向量

矩陣AA的特徵向量是在左乘矩陣AA後,方向不變的向量
Ax=λ x Ax=\lambda\ x
λ\lambda是特徵向量xx的特徵值

根據不可逆矩陣的定義,存在非零向量xx,滿足Ax=0Ax=0,所以如果一個矩陣是奇異的,那麼它必有特徵值00

例子

投影矩陣P=A(ATA)1ATP=A(A^TA)^{-1}A^T
對於vC(A)v\in C(A)Pv=vPv=v,所以所有在AA列空間(投影到的子空間)的向量都是特徵向量,特徵值爲11
對於wN(AT)w\in N(A^T)Pw=0Pw=0,所以所有與AA列空間垂直的向量都是特徵向量,特徵值爲00

算法

Ax=λ x[AλI]x=0 Ax=\lambda\ x\rightarrow[A-\lambda I]x=0
可見xN(AλI)x\in N(A-\lambda I)
在計算中,常常首先計算det [AλI]det\ [A-\lambda I],先求出λs\lambda s
對於矩陣
A=[3113] A=\left[ \begin{matrix} 3&1\\ 1&3\\ \end{matrix} \right]
AλI=[3λ113λ] A-\lambda I=\left[ \begin{matrix} 3-\lambda&1\\ 1&3-\lambda\\ \end{matrix} \right]
det [AλI]=(3λ)21=0 det\ [A-\lambda I]=(3-\lambda)^2-1=0
得到λ1=2,λ2=4\lambda_1=2,\lambda_2=4
分別將λ1=2,λ2=4\lambda_1=2,\lambda_2=4代回方程
[3λ113λ] \left[ \begin{matrix} 3-\lambda&1\\ 1&3-\lambda\\ \end{matrix} \right]
λi\lambda_i對應的特徵向量就在這個矩陣的零空間中
λ1=2\lambda_1=2
[1111] \left[ \begin{matrix} 1&1\\ 1&1\\ \end{matrix} \right]
x1=[1 1]Tx_1=[1\ -1]^T

λ2=4\lambda_2=4
[1111] \left[ \begin{matrix} -1&1\\ 1&-1\\ \end{matrix} \right]
x2=[1 1]Tx_2=[1\ 1]^T

AλIA-\lambda I的意義是什麼
矩陣
[0110] \left[ \begin{matrix} 0&1\\ 1&0\\ \end{matrix} \right]
其特徵值爲λ1=1,λ2=1\lambda_1=-1,\lambda_2=1,特徵向量爲x1=[1 1]Tx_1=[1\ -1]^Tx2=[1 1]Tx_2=[1\ 1]^T
對比上邊的例子,可以發現AλIA-\lambda I並不會影響特徵向量,而是將所有特徵值減去了λ\lambda

複數特徵值

90°90\degree旋轉矩陣
Q=[0110] Q= \left[ \begin{matrix} 0&-1\\ 1&0\\ \end{matrix} \right]
QλI=[λ11λ] Q-\lambda I= \left[ \begin{matrix} -\lambda&-1\\ 1&-\lambda\\ \end{matrix} \right]
det [QλI]=λ2+1=0 det\ [Q-\lambda I]=\lambda^2+1=0
可得λ1=i\lambda_1=-iλ2=i\lambda_2=i
注意特徵值是純虛數,而矩陣是反對稱的

虛數意味着某種旋轉

重複的特徵值

矩陣
A=[3103] A= \left[ \begin{matrix} 3&1\\ 0&3\\ \end{matrix} \right]
det [AλI]=(3λ)21=0 det\ [A-\lambda I]=(3-\lambda)^2-1=0
這個矩陣,有兩個相等的特徵值λ1=λ2=3\lambda_1=\lambda_2=3,也稱特徵值λ=3\lambda=3代數重度是2。
但是對於這個矩陣的特徵向量
A3I=[0100] A-3I= \left[ \begin{matrix} 0&1\\ 0&0\\ \end{matrix} \right]
發現只有一個特徵向量[1 0]T[1\ 0]^T
重複的特徵值可能會造成特徵向量的短缺,但是不是所有有重複特徵值的矩陣都會造成特徵向量的短缺,對於這樣矩陣的重複的特徵值,其對應的特徵空間維度等於其特徵值的代數重度。

特徵值的性質

nnn*n矩陣有nn個特徵值
特徵值之和等於矩陣的跡(trace)(trace)
特徵值之積等於矩陣行列式
復特徵值成對出現,互爲共軛

對角化和矩陣冪

假設矩陣AA存在nn個相互獨立的特徵向量,構造成特徵向量矩陣
S=[x1 x2xn] S=[x_1\ x_2\cdots x_n]
AA左乘SS
AS=[λ1x1 λ2x2λnxn] AS=[\lambda_1x_1\ \lambda_2x_2\cdots\lambda_nx_n]
分離特徵值和特徵向量
AS=[x1 x2xn][λ1λ2λn]=SΛ AS=[x_1\ x_2\cdots x_n] \left[ \begin{matrix} \lambda_1&&&\\ &\lambda_2&&\\ &&\ddots&\\ &&&\lambda_n\\ \end{matrix} \right]=S\Lambda
左乘S1S^{-1}
A=SΛS1 A=S\Lambda S^{-1}

對角化的條件

上一小節可以看到,不是所有矩陣都有數量足夠的線性無關的特徵向量,而可對角化的條件,正是要求nn階矩陣有nn個線性無關的特徵向量。(注意這裏不要求AA可逆)

矩陣冪

對角化提供了計算矩陣冪次的簡便算法,當然前提是矩陣可對角化(diagonalizable)(diagonalizable)
An=SΛS1SΛS1SΛS1=SΛnS1 A^n=S\Lambda S^{-1}S\Lambda S^{-1}\cdots S\Lambda S^{-1}=S\Lambda^n S^{-1}
上式說明,AnA^n相比於AA,特徵向量沒有變,但是特徵值變成λn\lambda^n

矩陣的穩定性

如果矩陣的所有特徵值絕對值都小於1,那麼
n,AnO n\rightarrow\infin,A^n\rightarrow O
稱矩陣AA是穩定的

相似性

如果M\exists M是可逆矩陣,使得矩陣A,BA,B滿足,B=M1AMB=M^{-1}AM
相似性完成了對同階方陣的一個劃分

再議特徵值分解

S1AS=Λ S^{-1}AS=\Lambda
可對角化的矩陣,就是相似與某個對角陣的矩陣

相似矩陣的性質

  1. 矩陣相似是一種等價關係
    a) 相似矩陣與自身相似
    b) 如果矩陣AA相似於矩陣BB,則矩陣BB相似於矩陣AA
    c) 如果矩陣AA相似於矩陣BB,矩陣BB相似於矩陣CC,則矩陣AA相似於矩陣CC
  2. 若矩陣AA相似於矩陣BB,則AABB有相同的特徵值
  3. 若矩陣AA相似於矩陣BB,則矩陣f(A)f(A)相似於矩陣f(B)f(B),當AA可逆,矩陣A1A^{-1}相似於矩陣B1B^{-1}

相似矩陣的特徵值與特徵向量

相似矩陣具有相同的特徵值
證明
如果AA相似於BB
Bx=λx Bx=\lambda x
M1AMx=λxM^{-1}AMx=\lambda x
AMx=λMxAMx=\lambda Mx
證畢,同時,可以看出雖然相似矩陣的特徵值相等,但是特徵向量發生了變化,不過可以發現兩個相似的矩陣必有相同數目的線性無關特徵向量
但是這是相似矩陣的判據嗎?
考慮矩陣
[3103][3003] \left[ \begin{matrix} 3&1\\ 0&3\\ \end{matrix} \right] \left[ \begin{matrix} 3&0\\ 0&3\\ \end{matrix} \right]
先上結論,這個兩個矩陣不相似,但是他們具有相同的特徵值3,而且代數重度是2,左邊矩陣只能找到一個線性無關的特徵向量,而右邊的矩陣可以找到兩個線性無關的特徵向量。

其實右邊的對角陣,在利用相似性進行的矩陣劃分中,他所在的劃分只有他一個元素,換言之,對於特徵值爲λ\lambda,代數重度爲nn,具有nn個線性無關的特徵向量的nn階方陣,除自身外無其他矩陣與其相似。

左邊的矩陣中,其實在上一節已經知道,它是無法對角化的,當然它與右邊的對角陣不相似。但是在它所在的劃分中,它是最標準的形式,被稱爲Jordan formJordan\ form(諾爾當標準型)。

JordanJordan標準型

相似的矩陣都可以被表示成同一個JordanJordan標準型

再考慮矩陣
[11][11] \left[ \begin{matrix} &1&&\\ &&1&\\ &&&\\ &&&\\ \end{matrix} \right] \left[ \begin{matrix} &1&&\\ &&&\\ &&&1\\ &&&\\ \end{matrix} \right]
這兩個矩陣都有相同的特徵值00,和相同的代數重度44,而且都有兩個線性無關的特徵向量,但是其實兩者並不相似

JordanJordan

JordanJordan塊形如
[λi1λi1λi] \left[ \begin{matrix} \lambda_i&1&&\\ &\lambda_i&\ddots&\\ &&\ddots&1\\ &&&\lambda_i\\ \end{matrix} \right]
JordanJordan塊只有一個特徵向量
上述例子的兩個矩陣可以分塊爲
[J33O31O13J11][J22O22O22J22] \left[ \begin{matrix} J_{3*3}&O_{3*1}\\ O_{1*3}&J_{1*1}\\ \end{matrix} \right] \left[ \begin{matrix} J_{2*2}&O_{2*2}\\ O_{2*2}&J_{2*2}\\ \end{matrix} \right]
可見這兩個方陣無法得到相同的分塊,JordanJordan認爲這樣的矩陣不相似

JordanJordan定理

每個方陣AA都相似與一個JordanJordan矩陣
JordanJordan矩陣形如
J=[J1J2Ji] J= \left[ \begin{matrix} J_1&&&\\ &J_2&&\\ &&\ddots&\\ &&&J_i\\ \end{matrix} \right]
JJ_*JordanJordan
而且
#Jordan blocks=#eigenvectors \#Jordan\ blocks=\# eigenvectors
可對角化的矩陣的JordanJordan矩陣就是對角陣Λ\Lambda

發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章