MATLAB Answers

2D matrix indexing for element wise operation

15 views (last 30 days)
Raj Raj
Raj Raj on 19 Feb 2016
Edited: Star Strider on 19 Feb 2016
Hello all, I am using following code to estimate pairwise angle between for a vector. Where the angle should be between 0 to 2*pi. I am able to predict the angle. But I need your advise on implementing this two for loops. This is slowing the code. Looking forward for your attention. In the code, XYdata is a n x 2 matrix representing the Cartesian coordinates of individual point [xi, yi]. I presume this slow is because of two FOR loops.
for i2 =1:length(XYdata_new)
for j2 = 1:length(XYdata_new)
angle_ij(i2,j2)= Angle360(XYdata_new(i2,1), XYdata_new(i2,2), XYdata_new(j2,1),XYdata_new(j2,2));
end
end
Thanks in advance.

  3 Comments

Raj Raj
Raj Raj on 19 Feb 2016
function angleout = Angle360(x1,y1,x2,y2)
if x1 == x2 && y1 ==y2
angleout = Inf;
end
angle = atan((y2-y1)/(x2-x1));
if angle > 0
if (y2 < y1) % y2 > y1, x2 > x1 % dy > 0, dx > 0 %first quadrant
angleout=pi+ angle;
else %
angleout= angle; % 3rd quadrant
end
elseif angle < 0
if x2 < x1 && y2 > y1
angleout= pi + angle; % 2nd quadrant
elseif y2 < y1 && x2 > x1
angleout= 3*pi/2 + (pi/2+angle); % 4th quadrant
elseif y2 < y1 && x2-x1 == 0
angleout= pi - angle; % 4th quadrant
end
elseif angle == 0
if x2 < x1
angleout = pi +angle;
else
angleout = angle;
end
end

Sign in to comment.

Accepted Answer

Star Strider
Star Strider on 19 Feb 2016
I would use the atan2 or atan2d functions. They are vectorised, so you can likely give your matrix columns to them directly and then do the angle difference calculations.

  8 Comments

Show 5 older comments
Star Strider
Star Strider on 19 Feb 2016
In the matrix bsxfun calculates, the first row are the angle differences between the first angle an all the others (so ‘AngleDifferences(1,2)’ is ‘Angles360(1)-Angles360(2)’ and so for the rest. This is the angular equivalent of the pdist function.
If you have two different sets of angles, run ‘Angles’ and ‘Angles360’ on each set separately, producing for example ‘Angles360A’ and ‘Angles360B’. The bsxfun call then becomes:
AngleDifferences = bsxfun(@minus, Angles360A(:), Angles360B(:)'); % Pairwise Angle Differences
The ‘AngleDifferences(1,2)’ matrix element would then become for example ‘Angles360A(1)-Angles360B(2)’, and so for the rest. (Note that in this line I forced them to both be column vectors.) (Make the appropriate changes to my code if you are using radian measure.)
This the angle equivalent of pdist2, so it should do what you want.
Raj Raj
Raj Raj on 19 Feb 2016
Thanks Star, I have marked it as answered.
Star Strider
Star Strider on 19 Feb 2016
My pleasure!
Thank you.
EDIT More efficient versions of the angle conversion lines are:
Angles360 = rem(360+a, 360);
Angles2pi = rem(2*pi+a, 2*pi);

Sign in to comment.

More Answers (0)

Tags

Community Treasure Hunt

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

Start Hunting!