SPARSE函数
S = sparse(A)
S = sparse(i,j,s,m,n,nzmax)
S = sparse(i,j,s,m,n), nzmax=length(s)
S = sparse(i,j,s), m = max(i), n = max(j)
S = sparse(m,n) == sparse([],[],m,n,0)
————————————————————————————————
备注: s 通常为一维向量 如 s = ones(N,1) 或 s = ones(1,N);
例如:
_______________________________________________
e =
D =
full(D)
_____________________________________________________
sparse用法: 在有限差分计算中,考虑4阶差分算法,其稀疏矩阵为
以6*6矩阵为例:
ans =
观察此矩阵,对角线系数为0,下三角第二对角线 -2/3,第三对角线 1/12, 末尾两位-1/12,2/3
算法:
N =6;
e = ones(N,1);D = sparse(1:N,[2:N 1],2*e/3,N,N)...
- sparse(1:N,[3:N 1 2],e/12,N,N);
D = (D-D');
rats(D)
差分算法 W = D*u/h