vectorisation a for loop

1 view (last 30 days)
Johnny
Johnny on 8 Dec 2019
Commented: Stephen23 on 8 Dec 2019
c=5;
retint=0;
dist=log([10:-1:1]+retint);
for i=1:length(dist)
eta=exp(-c*abs(dist(i)-dist));
discrim(i)=1/sum(eta);
end
Does anyone know how to vectorise this for loop to make it more efficient?
  1 Comment
Stephen23
Stephen23 on 8 Dec 2019
Because you are not concatenating anything, square brackets are not needed here:
dist=log((10:-1:1)+retint);

Sign in to comment.

Accepted Answer

David Hill
David Hill on 8 Dec 2019
Although arrayfun is really a loop,
c=5;
retint=0;
dist=log([10:-1:1]+retint);
discrim=arrayfun(@(x)1/sum(exp(-c*abs(dist-x))),dist);

More Answers (1)

Stephen23
Stephen23 on 8 Dec 2019
Real vectorized code (no loop or arrayfun):
eta = exp(-c*abs(bsxfun(@minus,dist,dist(:))));
discrim = 1./sum(eta,1)
Or for MATLAB versions >=R2016b:
eta = exp(-c*abs(dist-dist(:)));
discrim2 = 1./sum(eta,1)

Categories

Find more on Creating and Concatenating Matrices in Help Center and File Exchange

Community Treasure Hunt

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

Start Hunting!