%產生稀疏矩陣并轉變為全元素矩陣，得出的稀疏矩陣分布圖
a=eye(3);
a(4,:)=[-5 -2 -3]
b=sparse(a)  %創建稀疏矩陣
b=sparse([1 4 2 4 3 4],[1 1 2 2 3 3],[1 -5 1 -2 1 -3]) %創建與b相同的稀疏矩陣
full(b)
D=[3 2 9;2 4 9;1 1 4]
d=[0 1 2]
s=spdiags(D,d,4,3)
full(s)
 load spr.dat  			%裝載文本文件
b= spconvert(spr)  		%轉換成稀疏矩陣
 save spr b  				%保存spr.mat文件
 clear
 load spr b  				%裝載spr.mat文件
 b
 nnz(b)
 nonzeros(b)
 nzmax(b)
 spy(b)