一、矩阵运算的范例.doc

PAGE  PAGE 4 一、矩陣運算的範例 輸入Matlab指令說明執行結果a=[1 2 3;4 5 6;7 8 9]輸入3 ×3矩陣a1 2 3 4 5 6 7 8 9b=[1 2 4;6 8 9]輸入2 ×3矩陣b1 2 4 6 8 9c=[1;4;5]輸入3 ×1矩陣c1 4 5d=[1 2; 4 6; 8 9]輸入3 ×2矩陣d1 2 4 6 8 9[c d]將c,d合併為一3 ×3的矩陣1 1 2 4 4 6 5 8 9a*b’計算17 49 38 118 59 187c*c’計算1 4 5 4 16 20 5 20 25a.^2將a矩陣的每個元素平方1 4 9 16 25 36 49 64 81b./2將b矩陣的每個元素除以23.5000 1.0000 2.0000 3.0000 4.0000 4.5000tril(a)取a矩陣的下三角矩陣1 0 0 4 5 0 7 8 9triu(a)取a矩陣的上三角矩陣1 2 3 0 5 6 0 0 9zeros(2,3)建立一2 ×3的零矩陣0 0 0 0 0 0ones(2,3)建立一2 ×3的矩陣且元素均為11 1 1 1 1 1eye(3)建立一3 ×3的單位矩陣1 0 0 0 1 0 0 0 1c*ones(1,3)建立一3 ×3的矩陣且各列元素相同1 1 1 4 4 4 5 5 5 二、解聯立方程式的範例 試解 方法一：以高斯-喬丹消去法求解 輸入Matlab指令執行結果a=[2 3 -4 ;5 -7 -3 ;1 5 11 ]2 3 -4 5 -7 -3 1 5 11b=[1;2;-3]1 2 -3augmtx=[a b]2 3 -4 1 5 -7 -3 2 1 5 11 -3c=rref(augmtx)1.0000 0 0 0.1338 0 1.0000 0 -0.0845 0 0 1.0000 -0.2465sol=c(:,4)0.1338 -0.0845 -0.2465 方法二：以逆矩陣法求解 輸入Matlab指令執行結果sol=inv(a)*b 0.1338 -0.0845 -0.2465 三、Matlab中冒號（：）的使用範例 1. 向量的產生 輸入Matlab指令說明執行結果x=1:5建立由1到5的列向量且增量為11 2 3 4 5y=0:pi/4:pi建立由0到π的列向量且增量為π/40 0.7854 1.5708 2.3562 3.1416z=6:-1:1建立由6到1的列向量且增量為-16 5 4 3 2 1a=[1 2 3];a(:)將列向量a變成行向量1 2 3 2. 子矩陣的產生 輸入Matlab指令說明執行結果a=[1 3 5 7 9;3 4 1 2 5;0 1 2 3 6]建立3 ×5的矩陣1 3 5 7 9 3 4 1 2 5 0 1 2 3 6b=-1*ones(3,3)建立3 ×3的矩陣且元素均為-1-1 -1 -1 -1 -1 -1 -1 -1 -1a(2,3)取a矩陣之第2列第3行的元素1a(3,4)取a矩陣之第3列第4行的元素3a(:,3)取a矩陣的第3行5 1 2a(2,:)取a矩陣的第2列3 4 1 2 5a(2:3,1:4)取a矩陣的第2列到第3列及第1行到第4行的子矩陣3 4 1 2 0 1 2 3a(:,[2 4 5])=b(:,