k=0;c=0;
for i=1:l-1
for j=i+1:l
k=k+1;
distance(k)=sqrt((positionplus(i,1)-positionplus(j,1))^2+(positionplus(i,2)-positionplus(j,2))^2+(positionplus(i,3)-positionplus(j,3))^2);
if(distance(k)>range1 && distance(k)<range2)
c=c+1;
a=[positionplus(i,1),positionplus(i,2),positionplus(i,3)];
b=[positionplus(j,1),positionplus(j,2),positionplus(j,3)];
position1(c,1)=a(1,1);position1(c,2)=a(1,2);position1(c,3)=a(1,3);
position2(c,1)=b(1,1);position2(c,2)=b(1,2);position2(c,3)=b(1,3);
end
end
end
你好,不用循环的话,可以考虑pdist来计算。
positionplus = rand(100,3); range1 = 0.1; range2=0.2;%因为题主没给数据我这里随机生成一些,便于调试
D = squareform(pdist(positionplus, 'squaredeuclidean'));%以距离平方和来计算距离矩阵
D = triu(D);%获取上三角
[i,j] = find(D>range1^2 & D<range2^2);%获取满足条件的两个点的编号i和j
position1 = positionplus(i,:);
position2 = positionplus(j,:);
一是考虑用其它循环语句;二是考虑用递归嵌套