Could you please check my code?

Question

hrushikesh kyathari 2019-6-29

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/469482-could-you-please-check-my-code

评论： Walter Roberson 2019-6-30

R=[1,-0.5,-0.5];
for i=2:25
a = -0.50;
b = 0.50;
R = [R;i,(b-a).*rand(1,1) + a,(b-a).*rand(1,1) + a];
end
R=[R;26,0.5,0.5];
disp(R);
[D,I] = pdist2(R(:,(2:3)),R(:,(2:3)),'euclidean','Smallest',3);
D(1,:)=[];
I(1,:)=[];
disp(D);
disp(I);
Edge=zeros(52,3);
for i=1:26
  %intersect is a function written correctly!
  if intersect(R(i,2:3),R(I(1,i),2:3))==1
      continue
  else
      Edge(i,1)=i;
      Edge(i,2)=I(1,i);
      Edge(i,3)=D(1,i);
  end
end
for j=1:26
    
    if intersect(R(j,2:3),R(I(2,j),2:3))==1
      continue
  else
      Edge(j+26,1)=j;
      Edge(j+26,2)=I(2,j);
      Edge(j+26,3)=D(2,j);
  end
end
disp(Edge);
Q=zeros(1,3);
for k=1:51
    if Edge(k,:)==Q
        Edge(k,:)=[];
    end
end
 disp(Edge); 

The error I am receiving is:

Index in position 1 exceeds array bounds (must not exceed 42).

Error in ehe (line 39)

if Edge(k,:)==Q

3 个评论
显示 1更早的评论隐藏 1更早的评论

hrushikesh kyathari 2019-6-29

编辑：hrushikesh kyathari 2019-6-29

Whenever I run the code, I am getting a change in error. Sometimes it is

Index in position....exceed 42,36,40.

I feel something is wrong with the last for loop used to delete empty rows,cause when I remove the loop, the code runs perfectly . Any suggestions for that.

the cyclist 2019-6-29

I have run your code repeatedly, and get no error.

Do you have the problem on a fresh restart of MATLAB?

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Walter Roberson 2019-6-29

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/469482-could-you-please-check-my-code#answer_381321

在 MATLAB Online 中打开

mask = any(Edge,2);
Edge = Edge(mask,:);

or

mask = ~any(Edge,2);
Edge(mask,:) = [];

or

for k=51:-1:1
    if Edge(k,:)==Q
        Edge(k,:)=[];
    end
end

2 个评论
显示无隐藏无

hrushikesh kyathari 2019-6-29

编辑：hrushikesh kyathari 2019-6-29

Thank you, Mr Roberson, your answer has solved my problem. Could you please explain what is the difference between the two for loops, so that in future I can solve the bug myself.