MATLAB程式设计进阶章线性代数.pptVIP

* * 本章指令彙整 與線性代數相關的函式，彙整如下： 本章指令彙整 (2) 與線性代數相關的函式，彙整如下： 本章指令彙整 (3) 與線性代數相關的函式，彙整如下： * * * * * * * * * * * * * * * * * * * * * * * * * * * * * MATLAB 程式設計進階篇：線性代數 MATLAB 程式設計：進階篇線性代數 張智星 (Roger Jang) jang@ /jang 台大資工系 多媒體檢索實驗室 反矩陣 反矩陣： 一個矩陣 A 的反矩陣可表示成 ，滿足下列恆等式： 只有在 A 為方陣時， 才存在。 若 不存在，則 A 稱為 Singular  Quiz! 反矩陣 (2) inv： MATLAB 的 inv 指令可用於計算反矩陣 範例6-1：inv01.m 結果： B = 4.0000 -6.0000 4.0000 -1.0000 -6.0000 14.0000 -11.0000 3.0000 4.0000 -11.0000 10.0000 -3.0000 -1.0000 3.0000 -3.0000 1.0000 Note that det(pascal(n)) is always equal to 1. A = pascal(4); % 產生 4x4 的 Pascal 方陣 B = inv(A) 反矩陣 (3) inv： 若矩陣 A 為 Singular （即其反矩陣不存在），則在使用 inv 指令時，會產生警告訊息 範例6-2：inv02.m 結果： A = [1 2 3; 4 5 6; 7 8 9]; B = inv(A) 行列式 det： 欲計算矩陣 A 的行列式，可用 det 指令 範例6-3：det01.m 結果：d = 29 A = [1 3 4; -3 -4 -1; 2 2 5]; d = det(A) 反矩陣公式 由 Crammer Rule 可知矩陣 A 的行列式和反矩陣有下列關係式： 其中 代表 Ａ 的行列式， 代表 Ａ 的 Adjoint Matrix Equivalent statements is singular does not exist Quiz! Quiz! Cofactor and Adjoint Take a 3x3 matrix A for example: Quiz! How to Verify the Formula? Quiz! Quiz Candidates If A is a 3x3 matrix, what is cofactor(A)? What is adj(A)? And can you verify that Quiz! 反矩陣與行列式之驗證 (1/2) 若 A 為整數矩陣，則 乘上 必為整數矩陣，可驗証如下： 範例6-4：det02.m 結果： A = [1 3 4; -3 -4 -1; 2 2 5]; det(A)*inv(A) 反矩陣與行列式之驗證 (2/2) 若 A 為整數矩陣，將 inv(A) 以有理形式（Rational Format，即分子和分母都是整數的分數）來表示，即可察覺出它和行列式的關係 從這裡可以很明顯的看出，inv(A) 中每個元素的分母值，就是 det(A) 範例6-5：det03.m  結果： A = [1 3 4; -3 -4 -1; 2 2 5]; format rat % 以有理形式表示數值 inv(A) format short % 改回預設的數值表示形式 固有值與固有向量 一個方陣 A 的固有向量（Eigenvector） 與固有值 （Eigenvalue） 的關係式如下：上式可化簡成 由於 x 不是一個零向量，因此 必須是 Singular，上式才會有解。當 是 Singular 時，其行列式為零： Quiz! 固有值分解 (1/2) 若 A 為 n×n 的矩陣，則上式為 的 n 次多項式 ，代表 將有 n 個解 ，滿足下列關係式：或可寫成矩陣形式： 其中 如果 存在，則由上式可得矩陣A的固有值分解（Eigenvalue decomposition）： Quiz! 固有值分解 (2/2) If A is symmetric, then Suppose A is