学步园

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