I have three vectors V1, V2, and V3. Their origin points are on the axes' origin. How could I determine whether V3 is between V1 and V2 when I move around counterclockwise from V1 to V2?
alt text http://www.freeimagehosting.net/uploads/1448ea8896.jpg
It can't be done with obtaining their angles and evaluating these kind of conditions (pseudo-code):
if angle(V3) > angle(V1) && angle(V3) < angle(V2)
printf("V3 is between V1 and V2")
else
printf("out of the interval")
To see its defect, suppose that the angle
function gives angles in the range of [-pi pi]. So, if angle(V1) = 120 (in degree), angle(V2) = -130 and angle(V3) = 150 then the answer (according to the above code) is "out of the interval" although if you move around counterclockwise from V1 to V2, it is between them.
You may suggest adding 2*pi to angle(V2) or something like that, but I've tried such things and it doesn't work.
I'm programming in MATLAB.
EDIT 1 : it is in 2D.
Since you are doing this in MATLAB, here is one solution that should work:
crossProds = [V1(1)*V2(2)-V1(2)*V2(1) ...
V1(1)*V3(2)-V1(2)*V3(1) ...
V3(1)*V2(2)-V3(2)*V2(1)];
if (all(crossProds >= 0) || ...
(crossProds(1) < 0) && ~all(crossProds(2:3) < 0)),
disp("V3 is between V1 and V2");
else
disp("out of the interval");
end
EXPLANATION:
The cross product between the 2-D vectors V1 and V2 is stored in the first element of crossProds. This value will be greater than or equal to zero if the counter-clockwise angle between V1 and V2 is between 0 and 180 degrees, inclusive. In this case, when V3 is between V1 and V2 in the counter-clockwise direction then the cross products (V1,V3) and (V3,V2) are also greater than or equal to zero. This explains the first logical check:
all(crossProds >= 0)
If the counter-clockwise angle between V1 and V2 is greater than 180 degrees, then the cross product of these two vectors will be less than zero. In this case, when V3 is between V1 and V2 in the clockwise direction then the cross products (V1,V3) and (V3,V2) are also less than zero. Therefore, if these cross products are not both less than zero then V3 must be between V1 and V2 in the counter-clockwise direction. This explains the next two logical checks:
(crossProds(1) < 0) && ~all(crossProds(2:3) < 0)
The above logical checks should cover all possible situations. The operators || and && are short circuit operators in MATLAB: they will skip the second statements if they are not necessary. For example, if the first statement in an OR is true, there is no reason to check the second statement since only one argument in an OR needs to be true for the result to be true.
Calculate angle(V1), angle(V2) and angle(v3) (a1, a2, a3).
Modify a2 and a3 (add 2*pi if needed) so that
a1 <= a2 < a1 + 2*pi
a1 <= a3 < a1 + 2*pi
Now you simply have to compare a2 and a3. V3 is between V1 and V2 is resulting a3 is inferior to a2.
If you love us? You can donate to us via Paypal or buy me a coffee so we can maintain and grow! Thank you!
Donate Us With