スパース行列のグラフィカルな表現

有限要素メッシュの表示

まず、x および y を、範囲が $\left[0,1\right]$ に収まるように $2^{-32}$ でスケーリングします。次に、(i,j) 結合からスパース隣接行列を作成し、これを正定値にします。最後に、(x,y) を頂点 (メッシュ点) の座標として使用して、隣接行列をプロットします。

% Scaling x and y
x = pow2(x,-32); 
y = pow2(y,-32);

% Forming the sparse adjacency matrix and making it positive definite
n = max(max(i),max(j));
A = sparse(i,j,-1,n,n);
A = A + A';
d = abs(sum(A)) + 1;
A = A + diag(sparse(d));

% Plotting the finite element mesh
gplot(A,[x y])
title('Airfoil Cross-Section')

スパース パターンの可視化

spy を使用すると、行列内の非ゼロの要素のパターンを可視化できるため、スパース行列内のスパース性のパターンを確認するには特に便利な関数です。spy(A) は、行列 A のスパース性のパターンをプロットします。

spy(A)
title('Airfoil Adjacency Matrix')

対称的な並べ替え - 逆 Cuthill-McKee

symrcm は、隣接行列の並べ替えに逆 Cuthill-McKee 手法を使用します。r = symrcm(A) は、A より対角に近い対角要素をもつ傾向がある A(r,r) となるような置換ベクトル r を返します。この形式は、"細長い" 問題から生じる行列の LU 分解またはコレスキー分解に対して、適切な並べ替えを行います。対称と非対称の行列の両方に対して機能します。

r = symrcm(A);
spy(A(r,r))
title('Reverse Cuthill-McKee')

対称的な並べ替え - 列置換

非ゼロ数の非減少のスパース行列で A の列を並べ替えるために、置換ベクトルを返す j = COLPERM(A) を使用します。この形式は、lu(A(:,j)) のような LU 分解の前処理の並べ替えとして有効な場合があります。

j = colperm(A);
spy(A(j,j))
title('Column Count Reordering')

対称的な並べ替え - 対称な近似最小次数

symamd は、対称な近似最小次数置換を与えます。対称正定値行列 A に対して、コマンド p = symamd(S) は S(p,p) が S よりもスパースなコレスキー因子をもつような置換ベクトル p を返します。しばしば、symamd は、対称不定値行列に対しても、問題なく機能します。

m = symamd(A);
spy(A(m,m))
title('Approximate Minimum Degree')