Calculate angles correctly between two vectors using the dot product.
31 views (last 30 days)
Show older comments
Hi all,
As shown below in the figure, I have two vectors,
and
, and I want to calculate the angle with the norm,
and
, between each of these vectors with respect to vector L using the dot product.
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/1582606/image.png)
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/1582611/image.png)
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/1582616/image.png)
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/1582621/image.png)
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/1582626/image.png)
Here is a simple code to calculate these angles:
close all;
clear all;
%define the reference vector
L = [-0.5 -0.5];
%v1
v1 = [-0.6788 0.3214];
%theta1 with the x axis
theta1 = acos((v1(1)*L(1)+v1(2)*L(2))/(sqrt(v1(1)^2+v1(2)^2)*sqrt(L(1)^2+L(2)^2)));
%theta1 with the norm
theta1norm = pi/2 - theta1;
%v2
v2 = [0.3214 -0.6788];
%theta2 with the x axis
theta2 = acos((v2(1)*L(1)+v2(2)*L(2))/(sqrt(v2(1)^2+v2(2)^2)*sqrt(L(1)^2+L(2)^2)));
%theta1 with the norm
theta2norm = pi/2 - theta2;
My question here is that the product produces both angles to be equal, i.e.,
. However, and as shwon in the figure, it is clear that
, i.e.,
and
. How can I caulcate these angles correctly using the dot product?
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/1582631/image.png)
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/1582636/image.png)
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/1582641/image.png)
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/1582646/image.png)
Thanks.
Answers (3)
Hassaan
on 4 Jan 2024
% Close all figures and clear variables
close all;
clear all;
% Define the reference vector
L = [-0.5 -0.5];
% Define the first vector v1 and calculate the angle with L
v1 = [-0.6788 0.3214];
% Dot product of v1 and L
dot_v1_L = dot(v1, L);
% Norms of v1 and L
norm_v1 = norm(v1);
norm_L = norm(L);
% Angle theta1 with the x axis
theta1 = acos(dot_v1_L / (norm_v1 * norm_L));
% Angle theta1 with the norm
theta1norm = pi/2 - theta1;
% Define the second vector v2 and calculate the angle with L
v2 = [0.3214 -0.6788];
% Dot product of v2 and L
dot_v2_L = dot(v2, L);
% Norms of v2
norm_v2 = norm(v2);
% Angle theta2 with the x axis
theta2 = acos(dot_v2_L / (norm_v2 * norm_L));
% Angle theta2 with the norm
theta2norm = pi/2 - theta2;
% Convert angles to degrees
theta1_degree = radtodeg(theta1);
theta2_degree = radtodeg(theta2);
theta1norm_degree = radtodeg(theta1norm);
theta2norm_degree = radtodeg(theta2norm);
% Display the results
disp(['Theta1: ', num2str(theta1_degree), ' degrees']);
disp(['Theta2: ', num2str(theta2_degree), ' degrees']);
disp(['Theta1 from the norm: ', num2str(theta1norm_degree), ' degrees']);
disp(['Theta2 from the norm: ', num2str(theta2norm_degree), ' degrees']);
The functions dot and norm are used to calculate the dot product and the magnitude of the vectors, respectively, and acos computes the arccosine of the given value to obtain the angle in radians. The function radtodeg converts the angle from radians to degrees.
------------------------------------------------------------------------------------------------------------------------------------------------
If you find the solution helpful and it resolves your issue, it would be greatly appreciated if you could accept the answer. Also, leaving an upvote and a comment are also wonderful ways to provide feedback.
Professional Interests
- Technical Services and Consulting
- Embedded Systems | Firmware Developement | Simulations
- Electrical and Electronics Engineering
0 Comments
James Tursa
on 4 Jan 2024
Edited: James Tursa
on 4 Jan 2024
Also see this related link for a robust method using atan2 that can recover small angles:
0 Comments
Torsten
on 4 Jan 2024
I think an intuitive way is to compute the angles between the positive x-axis and the respective vector counterclockwise first (the result will be between 0 and 360) and then make the necessary subtractions.
theta1 = cart2pol(L(1)/norm(L),L(2)/norm(L))*180/pi
if theta1 <=0
theta1 = 360 + theta1;
end
theta2 = cart2pol(v1(1)/norm(v1),v1(2)/norm(v1))*180/pi
if theta2 <=0
theta2 = 360 + theta2;
end
theta3 = cart2pol(v2(1)/norm(v2),v2(2)/norm(v2))*180/pi
if theta3 <=0
theta3 = 360 + theta3;
end
0 Comments
See Also
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!