3,423
社区成员
发帖
与我相关
我的任务
分享
function out=Sim(i,j) %i,j是两子矩阵存储的文件名
D1=dlmread(strcat(num2str(i),'.txt'));
D2=dlmread(strcat(num2str(j),'.txt'));
A1=spconvert(D1);
A2=spconvert(D2);
AF1=full(A1);
AF2=full(A2);
if size(AF1,2) < 531543 //稀疏矩阵存储,为0的未存储,将其补上
AF1=[AF1,zeros(size(AF1,1),531543-size(AF1,2))];
end
if size(AF2,2) < 531543
AF2=[AF2,zeros(size(AF2,1),531543-size(AF2,2))];
end
out=AF1*AF2';
Step1:
for i=0:179
A=Sim(i,0);
for j=1:179
A=[A,Sim(i,j)]; %横向拼接
end
dlmwrite(strcat('sim',num2str(i),'.txt'),A,'delimiter',' ');
end
Step2:
A=dmlread(strcat('sim0.txt'));
for i=1;179
A=[A;dmlread(strcat('sim',num2str(i),'.txt')); %纵向拼接
end
dlmwrite('sim.txt',A,'delimiter',' ');