一个简单的MATLAB读写文本文件程序--稀疏矩阵向量乘法

最新推荐文章于 2024-09-25 20:52:58 发布

原创最新推荐文章于 2024-09-25 20:52:58 发布 · 6.5k 阅读

5 ·

CC 4.0 BY-SA版权

文章标签：

#matlab #matrix

MATLAB 专栏收录该内容

9 篇文章

订阅专栏

本文介绍了一个使用MATLAB实现的稀疏矩阵与稀疏向量相乘的方法，并将稀疏矩阵及其乘积结果以三元组形式写入文件的过程。该程序首先生成一个150,000 x 150,000的稀疏矩阵和一个同样长度的稀疏向量，然后进行了矩阵向量乘法运算。

%矩阵向量乘，写入文件 %双重for循环处，运行了很久，还需要改进，请多关照 clear clc n=150000;%矩阵的维数 a = sprand(n, n, 0.0001);%产生一个稀疏矩阵 fid1 = fopen('matrix.mtx', 'w');%要写入的矩阵文件 fid2 = fopen('matrix_b.mtx', 'w');%向量文件 disp('1') len = nnz(a);%稀疏矩阵中非零元个数 disp('2') fprintf(fid1, '%d %d %d/n', n, n, len);%行、列、非零元个数 disp('3') for j = 1:n for i = 1:n if a(i,j) ~= 0 fprintf(fid1, '%d %d %f/n', i, j, a(i,j));%三元组形式写入文件 end end end disp('4') fprintf(fid2, '%d/n', n); b = sprand(n, 1, 0.9);%产生一个稀疏向量 disp('6') for i = 1:n fprintf(fid2, '%d/n', b(i,1));%写入向量 end disp('5') s = a *b; full_s = full(s);%将结果转换成全矩阵 save s.txt full_s -ascii %保存矩阵向量乘的结果 fclose(fid1); fclose(fid2);