科学计算实务matlab 08 稀疏矩阵与结构阵列 - 明新科技大学[ 电子 .pdfVIP

科學計算實務 MATLAB08 稀疏矩陣與結構陣列 主 講 ：呂 明 峰 副 教 授 明新科技大學電子工程系 E- mail : mflu@must.edu.tw Outline 稀疏矩陣 結構陣列 sparse指令 struct指令 spdiags指令 取用及改變結構陣列的資料 spconvert指令 Bracket Operation 稀疏矩陣的圖示 取用及改變結構陣列的欄位 稀疏矩陣的相關指令表 其他相關指令 科學計算實務 2 1 稀疏矩陣 一個矩陣中 ，如果包含很多零元素 ，此矩陣即可稱為稀疏 矩陣 （Sparse Matrix ）。稀疏矩陣在工程上的用途相當廣 泛 ，舉凡電路 、圖學、有限元素法及偏微分方程等 ，都會 用到稀疏矩陣。 但在稀疏矩陣中，由於大部分的元素都是0 ，因此只須儲 存 「非零元素的位置」及其 「元素值」即可 ，這種特殊的 內部儲存方式有下列兩大優點： 1. 節省記憶體儲存空間 2. 節省許多不必要的運算 （因為很多元素都是零） MATLAB 用三個內部陣列來儲存稀疏矩陣的相關資訊： 1. 第一個陣列：以double 方式儲存了所有的非零元素 。 2. 第二個陣列：以整數方式儲存了每個元素的列索引 。 3. 第三個陣列：以整數方式儲存了每個直行的起始指標 。 科學計算實務 3 sparse指令 MATLAB 的sparse 指令可將一個完全矩陣轉換成稀疏矩陣 ： A = [2 0 0 0; 0 0 0 1; 0 4 0 0] S = sparse(A) %將一個完全矩陣轉換成稀疏矩陣 S = (1,1) 2 (3,2) 4 (2,4) 1 由此可看出 ，S 是一個稀疏矩陣 ，MATLAB 只儲存其各個非零 元素的位置 （即其二維下標 「(1,1) 、(3,2) 、(2,4) 」）和元素值 （「2 、4 、1」）。若要比較矩陣A 和S 所佔用的記憶體大 小，可輸入如下 ： whos Name Size Bytes Class A 3x4 96 double array S 3x4 56 sparse array 科學計算實務 4 2 可以使用sparse 指令來直接產生稀疏矩陣 ，其使用語 法如下 ： S = sparse(i, j, s, m, n); 其中i 是列索引 ，j 是行索引 ，s 是非零元素所形成的向量 ，m 是 s 的列維度 ，n 是 s 的行維度。i 、j 、s 都是長度相同的向量 ，s(k) 的二維下標即是i(k)及j(k ) 。 例如，以下語法可以產生與前例相同的稀疏矩陣S ： S = sparse([1 3 2], [1 2 4], [2 4 1], 3, 4) S = (1,1) 2 (3,2) 4 (2,4) 1 科學計算實務 5 sp