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
matlab不用for循环怎么做?
- 写回答
- 好问题 0 提建议
- 追加酬金
- 关注问题
- 邀请回答
-
2条回答 默认 最新
- 技术专家团-Joel 2021-10-09 19:31关注
你好,不用循环的话,可以考虑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,:);
本回答被题主选为最佳回答 , 对您是否有帮助呢?解决 1无用
悬赏问题
- ¥15 matlab数字图像处理频率域滤波
- ¥15 在abaqus做了二维正交切削模型,给刀具添加了超声振动条件后输出切削力为什么比普通切削增大这么多
- ¥15 ELGamal和paillier计算效率谁快?
- ¥15 file converter 转换格式失败 报错 Error marking filters as finished,如何解决?
- ¥15 ubuntu系统下挂载磁盘上执行./提示权限不够
- ¥15 Arcgis相交分析无法绘制一个或多个图形
- ¥15 关于#r语言#的问题:差异分析前数据准备,报错Error in data[, sampleName1] : subscript out of bounds请问怎么解决呀以下是全部代码:
- ¥15 seatunnel-web使用SQL组件时候后台报错,无法找到表格
- ¥15 fpga自动售货机数码管(相关搜索:数字时钟)
- ¥15 用前端向数据库插入数据,通过debug发现数据能走到后端,但是放行之后就会提示错误