第 7 章 線性代數:矩陣,向量,行 列式,線性方程組
7.1 矩陣,向量:加法與純量乘積 7.2 矩陣乘法
7.3 線性方程組,高斯消去法
7.4 線性獨立,矩陣的秩,向量空間 7.5 線性系統的解: 存在性,唯一性 7.6 參考用: 二階與三階行列式
7.7 行列式,柯拉瑪法則
7.8 反矩陣,高斯—喬丹消去法
反矩陣,高斯喬丹消去法
將 n × n 矩陣 A = [ajk
]
的反矩陣以符號記為 A-1,反矩陣為 n × n 矩陣且滿足(1)
其中,I 為 n × n 單位矩陣(見 6.2 節)
若 A 有反矩陣,則 A 稱為非奇異矩陣(nonsingular matrix)。
若 A 沒有反矩陣,則 A 稱為奇異矩陣(singular matrix)。若 A 有 反矩陣,則此反矩陣為唯一。
事實上,若 B 與 C 均為 A 的反矩陣,則 AB
=
I 與CA=
I,導 致由下式證得唯一性定理 1 反矩陣的存在性
定理 1 證明
令 A 為已知 n × n 矩陣,並考慮線性系統 (2)
若反矩陣 A
-1存在,則在上式兩邊由左邊乘上 A
-1,且利用 (1) 式可得
相反地,令 rank A= n。那麼根據相同定理,系統 (2) 對任 意 b 有唯一解 x。接著依照高斯消去法(6.3 節)的反向代入,
顯示 x 的分量 x
j為 b 向量的線性組合。故可寫成
(3)
定理 1 證明
對任意 b 而言。故 C = AB = I,單位矩陣。同理,若將 (2) 式代入 (3) 式得到
對任意 x 而言(且b = Ax),故 BA = I。即存在 B =
A-1。
高斯喬丹法求反矩陣
對於實際求解非奇異 n × n 矩陣 A 的反矩陣 A-1 而言,我們可 使用高斯消去法(6.3 節),事實上是它的變種,稱為高斯—喬丹 消去法(Gauss-Jordan elimina
tion)。此解法觀念如下:
使用 A 來構成 n 個線性系統
其中 e
(1),…,e
(n)為 n × n 單位矩陣 I 的行;因此, e
(1)= [1 0 … 0]
Τ, e
(2)= [0 1 0 …]
Τ等。這些是以未知向量 x
(1), …, x
(n)
表示的 n 個向量方程式。再合併成單一矩陣方程式 AX = I,
其中此未知矩陣 x 具有 x
(1), …, x
(n)等行。然後,將對應的 n
個擴大矩陣[A e
(1)],…, [A e
(n)] 合併成一個 n × 2n 的擴大矩陣
à =[A I]。
那麼,將 AX = I 由左邊乘以
A-1而得 X =
A-1 I = A-1。 由AX = I 求解 X,我們可對 Ã =[A I] 使用高斯消去法而得 形式為 [U H] 的矩陣,其中 U 為上三角,因高斯消去法使系 統成為三角化。
高斯
—喬丹消去法進一步使用基本列運算,將 U 簡化成對 角形式,事實上成為單位矩陣 I。此做法為消除 U 主對角線 上的項,並以乘法使對角線項變為全部是 1(見下面範例)。
當然,此方法在整個矩陣 [U H]上運算,同時將 H 轉換成 某一矩陣 K,因此[U H] 轉換成 [I K]。這也是 IX = K 的擴 大矩陣。那麼如前所示 IX = X =
A-1。比較而得 K =
A-1
,如此一來,可由 [I K] 讀出
A-1。
求
的反矩陣
A-1。
解
應用高斯消去法(6.3 節)到下列 n × 2n = 3 × 6 的矩陣,
其中藍色恆指向先前的矩陣。
範例 1 反矩陣,高斯喬丹消去法
此即 [U H],如同依高斯消去法所產生一樣。那麼再依高 斯喬丹的步驟,將 U 減少至 I,亦即成為主對角線上均為 1 的對角矩陣。
範例 1 (續)
範例 1 (續)
最後三行構成
A-1,驗算
因此,A
A-1= I。同理,
A-1A = I。範例 1 (續)
定理 2 反矩陣
定理 2 證明
將 (4) 式右邊表為 B,並證明 BA = I。首先寫出 (5)
然後再證出 G = I。那麼,根據矩陣乘法的定義以及 (4) 式 右邊 B 的形式,得到
(6)
(注意!是 C
sk而非 C
ks)
定理 2 證明
那麼由 6.7 節中 (9) 與 (10) 式知,當 l = k 時,上式右邊的 和 (…) 為 D = det A;而當 l ≠ k ,此和為零。因此
尤其當 n = 2 時,由 (4) 式的第 1 列得到 C
11= a
22、C
21=
-a
12,以及第 2 列得 C
12= -a
21、C
22= a
11。如此得證 (4*)
式。
範例 2 2 × 2 矩陣的反矩陣
範例 3 定理2 的進一步說明
試用 (4) 式計算
解
計算得到 det A = -1(-7)-1.13+2.8 =10 以及 (4)
式中
範例 3 (續)
使得(4) 式的結果與範例 1 一致。
對角矩陣
對角矩陣(diagonal matrix)A = [ajk