I have two matrices.
Bondsbystep is a 64X30000 matrix listing bond lengths at each step of a calculation.
a is a 1x30000 matrix listing the important bond number in Bondsbystep.
I'm trying to write a script so that the bond length (from Bondsbystep) corresponding to the right bond (listed in a) at each step is listed in a 1x30000 matrix. I have done this with a set of data which is all real integers using:
for i= 1:30000;
x=bondsbystep(a(1,i),:);
end
However using a set of data where some a(1,i) values are real, and some are NaN produces the error:
Subscript indices must either be real positive integers or logicals.
Ideally I would have a 1x30000 matrix where some values are real bond lengths, and some are 0 (or nan). This will enable me to plot bond change over time.
I've tried to use an if function of the following:
x2=nan(1,numel(a));
for i=1:numel(a);
if a(1,i)==0;
x2(i)=Bondsbystep(a(1,i),:);
end
end
But I think I missing something in the way I'm relating the two matrices. Any help or pointers on where I should be look is much appreciated.
Edit:
I have also tried using a if/else loop:
for i=1:30000;
if a(1,i) > 0;
x2=Bondsbystep(c(1,i),:);
else
if a(1,i)==NaN;
x2(1,i)=0;
end
end
end
Which runs but doesn't give me the correct output.
I've tried using the continue function as well, which again runs but not the correct output.
Perhaps this is a clearer example of what I need:
for i=1:30000;
if a(1,i) > 0;
x2(1,i)=Bondsbystep(c(1,i),:);
elseif a(1,i)==0;
x2(1,i)=0;
end
end
I have tried the above but got the error:
Assignment has more non-singleton rhs dimensions than non-singleton subscripts
Many thanks for any input.
