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

 

%矩阵向量乘，写入文件%双重for循环处，运行了很久，还需要改进，请多关照clearclcn=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 endenddisp('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));%写入向量enddisp('5')s = a *b;full_s = full(s);%将结果转换成全矩阵save s.txt full_s -ascii %保存矩阵向量乘的结果fclose(fid1);fclose(fid2);