爱问基本MATLAB指令介绍.docVIP

基本MATLAB指令介紹 (I) I. 建立各種矩陣的指令 1. 直接輸入矩陣的方法：將數字置於[ ]中，並以分號” ; ”做為換行的記號。如果沒有分號，或以逗號分隔元素，均會產生一個列向量(row vector)。 試比較：a=[1 2 3;4 5 6;7 8 9]，a=[1 2 3 4 5 6 7 8 9]，a=[1;2;3;4;5;6;7;8;9]， a=[1, 2, 3, 4, 5, 6, 7, 8, 9] 2. 矩陣中可以放矩陣，這便是分割矩陣(partition matrix)的概念： 試比較：a=[1 2;3 4]，b=[5 6;7 8]，c=[a b]，d=[a; b] 3. 利用冒號” : ”可以從一個矩陣中擷取一子矩陣(submatrix)： 例如，若輸入a=[1 2 3;4 5 6;7 8 9]，則試比較：b=a(1:2,: )，c=a(:,1:2)，d=a(1:2,1:2) 4. 利用冒號” : ”可以使矩陣變成一行向量(column vector)： 試比較：a=[1 2;3 4;5 6]，b=a(:)，c=[1 2 3 4]，d=c(:) 5. 利用冒號” : ”可以等差方式建立列向量： 試比較：a=-1:0.5:3，x=1:5，z=6:-1:1 6. 建立元素均為1的矩陣：例如，ones(3,3)，ones(5,2) 7. 建立元素均為0的矩陣：例如，zeros(3,3)，zeros(5,2) 8. 建立單位矩陣：例如，eye(4)代表一個4╳4的單位矩陣 9. 建立對角線元素為的對角線矩陣：diag() 10. 擷取任意矩陣a的對角線元素：diag(a) 11. 擷取任意矩陣a的上三角部分：triu(a)；擷取任意矩陣a的下三角部分：tril(a) 12. 建立n×n之隨機矩陣：rand(n)；建立m×n之隨機矩陣：rand(m,n) II. 矩陣運算的指令 1. 基本算數指令：+、－、×、\ 2. 任意矩陣a的轉置矩陣(transpose)：a′ 3. 非奇異矩陣a的逆矩陣(inverse)：inv(a) 4. 任意方陣a的行列式(determinant)：det(a) 5. 兩個矩陣之元素間的對應運算：.* ./ .^ 例如，a=[1 2;3 4]，b=[5 6;7 8]，試求a.*b、a./b、a.^2 6. 將矩陣a化為Reduced Row Echelon Form：rref(a) 7. 若解線性系統Ax=b，且A為一非奇異方陣，則x = A\b = inv(A)*b 8. 查詢矩陣a的列數與行數：size(a) 或 [m,n]=size(a) 9. 查詢列向量或行向量X的元素數目：length(X) III. 向量空間之運算 1. 任意矩陣a的rank：rank(a) 2. 任意矩陣a之null space的基底向量：null(a) 3. 任意兩個向量a，b之內積：dot(a,b) 4. 任意兩個向量a，b之外積：cross(a,b) IV. 特徵值與特徵向量的運算 1. 任意矩陣a的特徵值：eig(a) 2. 任意矩陣a的對角化矩陣(D)及其特徵向量所構成之矩陣(P)：[P,D]=eig(a) (注意相似矩陣的定義：D=P-1AP) 3. 任意矩陣a之特徵多項式的係數，即det(a-λI)=0展開後形成λ之多項式的係數(依降冪方式排列)：poly(a) 注意事項： 1. 在MATLAB指令視窗(command window)中執行指令時，若在指令尾端打上分號，則執行結果不會顯示在視窗上，若不打分號，才可見到執行結果。 2. MATLAB設有help指令供使用者查詢指令的功能及使用方法。 3. 若欲建立可執行之MATLAB程式，應使用NOTEPAD建立附檔名為M之檔案，此工作可在指令視窗上方的File選單中選擇New內的M-file進行。 建議養成建立副程式的習慣，以使程式結構化。在MATLAB中是以function代表副程式，其寫法為：function [ 輸出參數 ] = 副程式名稱( 輸入參數 )。 參考書籍： 1. MATLAB使用入門手冊，歐崇明，高立圖書公司，1998。 2. MATLAB程式設計 基礎篇，鄭錦聰，全華圖書公司，2000。 空格号