清空筛选条件
清空筛选条件

Find Angle Between Two Line

5 次查看（过去 30 天）
mnhafizuddin
mnhafizuddin 2012-5-18
How i want to find the angle between two line.. Example my code:
%%%Angle Between Two Line%%%
% Line 1: point(0,3) to (4,3)
% Line 2: point(0,2) to (3,0)
v1=[4,3]-[0,3];
v2=[0,2]-[3,0];
angle = mod( atan2( det([v1,v2]) , dot(v1,v2) ) , 2*pi );
what i get was
Error using ==> det
Matrix must be square..
Why I cannot get the angle.???

请先登录，再回答此问题。

回答（2 个）

Andrei Bobrov
Andrei Bobrov 2012-5-18
v1=[4,3]-[0,3];
v2=[0,2]-[3,0];
angle = mod( atan2( det([v1;v2]) , dot(v1,v2) ) , 2*pi );
EDIT
a1 = mod(atan2( det([v1;v2;]) , dot(v1,v2) ), 2*pi );
angleout = abs((a1>pi/2)*pi-a1)
  3 个评论
显示 1更早的评论
Jan
Jan 2012-5-18
编辑：Jan 2012-12-22
@mnhafizuddin: I don't understand the question. The problem of your code is the comma in "[v1,v2]", bcause it creates a [1 x 4] vector, but DET needs a [2x2] matrix. See some of the most popular thread in CSSM: http://www.mathworks.com/matlabcentral/newsreader/view_thread/151925#870338
mnhafizuddin
mnhafizuddin 2012-5-19
thanks for your suggestion website..
It's really helps me alot..
thanks!!!

请先登录，再进行评论。

Geoff
Geoff 2012-5-18
When I want to find the angle between two vectors, I take the dot product of the unit vectors:
clamp = @(val, low, high) min(max(val,low), high);
angle = acos( clamp(dot(v1,v2) / norm(v1) / norm(v2), -1, 1);
If you want to know which direction that's in, you can do a trick with the cross-product. Since it's 2D, you do this:
vc = cross([v1,0], [v2,0]);
anticlock = vc(3) > 0;
  3 个评论
显示 1更早的评论
Jan
Jan 2012-5-18
ACOS(DOT(u,v)) is instable if u and v are near to (anti)-parallel. ASIN(CROSS()) is not sufficient near to perpendicular vectors. Therefore the ATAN2 approach is recommended.
Geoff
Geoff 2012-5-21
Oh righto. I never had issues with stability. But that's what the clamp is for. It prevents the value from exceeding 1 (or -1) due to floating-point precision errors.

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Creating and Concatenating Matrices 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by