我有以下数据,我想在for循环中使用前两列来检测没有接触的粒子。前两列显示了两个接触的粒子。我如何得到一个表格形式的输出。基本上,我想检测那些没有接触的粒子。粒子数为484,总接触数为956。
1324 0.720141865396207E+00 0.000000000000000 E+00 0.000000000000000 E+00
1 439 0.940050831127248e +00 0.000000000000000e +00 0 0.000000000000000e +00
1441 0.1316099985588E+01 0.000000000000000 E+00 0.000000000000000 E+00
2 25 0.363379567116406e +00 -0.000000000000000e +00 0 0.000000000000000e +00
289 0.75309039602638E+00 0.000000000000000 E+00 0.000000000000000 E+00
21090.851303474754872E+00.000000000000000 E+00.000000000000000 E+00
2200 0.898508917122121E+00-0.000000000000000 E+00 0.000000000000000 E+00
36 0.890341369366110E+00 0.000000000000000 E+00 0.000000000000000 E+00
326 0.283864328102457E+00 0.000000000000000 E+00 0.000000000000000 E+00
数据=导入数据(“dem22.txt”);
Particle1=数据(:,1);
Particle2 =数据(:,2);
粒子数=484
total_contacts =长度(Particle1);
= Particle1 + Particle2联系
表=单元(956,2);
对于i=1:1:ncont;
j=1:npa;
触点(ncont,1)=i;
i=i+1;
触点(ncont,2)=j;
j=j+1;
表(i,:)={ij};
fprintf(“%d%d\n”,表{i,:});
终止
T=表格(联系人,npa);