Vectorization of a for loop

2 次查看(过去 30 天)
jpedro
jpedro 2015-11-9
编辑: Guillaume 2015-11-10
Hello,
I am trying to eliminate the loop and vectorise this piece of code. Any ideas?
for k = 1:(length(P)-2)
if ((P(k+1)-P(k)>0)&&(P(k+2)-P(k+1)<0));
B(k+1) = 1;
end
end
Thanks

回答(2 个)

Peter O
Peter O 2015-11-9
v = 1:length(P)-2;
ix = find( ((P(v+1) - P(v)) > 0) & ((P(v+2) - P(v+1)) < 0) ) + 1;
B2(ix) = 1;
  1 个评论
jpedro
jpedro 2015-11-10
Hello Peter,
Thanks for your answer, but that code doesn't give the same results as the for above

请先登录,再进行评论。


Guillaume
Guillaume 2015-11-10
编辑:Guillaume 2015-11-10
Well, P(k+1)-P(k) can be obtained in a vectorised way with the diff function, so, assuming that by default B(k) is zero:
B = diff(P(1:end-1)) > 0 & diff(P(2:end)) < 0 %no need to preinitialise B
If by default, B(k) is something else:
B(diff(P(1:end-1)) > 0 & diff(P(2:end)) < 0) = 1; %B must be preinitialised

类别

Help CenterFile Exchange 中查找有关 Loops and Conditional Statements 的更多信息

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by