matlab源代码为:
1、文件 test_bcsstk01.m:
% 绘制稀疏矩阵 bcsstk01.rsa 和 bcsstm01.rsa 的结构图
clear
% 读入rsa形式的刚度矩阵
Dpi = 150;
hb_to_mmm('matrix\bcsstk01.rsa','matrix\bcsstk01');
load matrix\bcsstk01;
ni=bcsstk01(1,1);
nj=bcsstk01(1,2);
nz=bcsstk01(1,3);
nl=nz+1;
KL=sparse( bcsstk01(2:nl,1), bcsstk01(2:nl,2), bcsstk01(2:nl,3),
ni, nj );
[dr,va]=get_diag(bcsstk01(2:nl,1),bcsstk01(2:nl,2),bcsstk01(2:nl,3),nz);
KD=sparse( dr, dr, va, ni, nj
);
K = KL + KL' -
KD;
clear KL;
clear KD;
clear dr;
clear va;
figure(1)
subplot(1, 2, 1)
spy(K)
title('矩阵 bcsstk01.rsa 的结构图');
% 读入rsa形式的质量矩阵
hb_to_mmm('matrix\bcsstm01.rsa','matrix\bcsstm01');
load matrix\bcsstm01;
ni=bcsstm01(1,1);
nj=bcsstm01(1,2);
nz=bcsstm01(1,3);
nl=nz+1;
ML=sparse( bcsstm01(2:nl,1), bcsstm01(2:nl,2), bcsstm01(2:nl,3),
ni, nj );
[dr,va]=get_diag(bcsstm01(2:nl,1),bcsstm01(2:nl,2),bcsstm01(2:nl,3),nz);
MD=sparse( dr, dr, va, ni, nj
);
M = ML + ML' -
MD;
clear ML;
clear MD;
clear dr;
clear va;
subplot(1, 2, 2)
spy(M)
title('矩阵 bcsstm01.rsa 的结构图');
SDpi=['-r',num2str(Dpi)];
print(1,'-dtiff', SDpi, 'bcsstk01_bcsstm01'
);
2、文件 get_diag.m:
function [dr,va]=get_diag(row,col,val,nz)
index=1;
sum =0;
for i=1:nz
end
if( sum > 0 )
end
end
3、文件 hb_to_mmm.m:
下载网址:
参考文献:
【总结】Harwell-Boeing格式矩阵转换为稀疏矩阵的方法
注:
【问题】昨天这个程序似乎有些问题:hb_to_mmm.m文件得到的非零元个数似乎比实际的要少。
例如:bcsstk09.rsa这个文件,如果只考虑下三角部分,实际上共有9760个非零元,而若先用hb_to_mmm.m文件得到一个matlab格式的矩阵文件bcsstk09,再读入用spconvert函数转化成稀疏矩阵之后,虽然仍有9760个非零元,但其中许多非零元的数值却被设置为0了(直接查看bcsstk09.rsa文件,发现这些位置处的非零元数值为e-7/e-8量级)!这样,当再调用spy函数画出刚度矩阵下三角的非零元个数时,就只有9039个了。
【解决方案】修改hb_to_mmm.m文件的msm_to_mm_coordinate_real_general函数,将第1445行:
fprintf ( fid, '
改为
fprintf ( fid, '
即可。
【结果】