0001 0002 function [AA,ann] = findAngle(v1, v2) 0003 0004 normera = norm(v1)*norm(v2); 0005 0006 if normera == 0 0007 angle = 0; 0008 else 0009 angle = acos (dot(v1,v2)/normera)* (180/pi); 0010 end 0011 0012 a = angle; 0013 if v1(1) > 0 0014 if v2(2) < 0 0015 a = 360 - angle; 0016 else 0017 a = angle; 0018 end 0019 end 0020 0021 if v1(1) < 0 0022 if v2(2) > 0 0023 a = 360 - angle; 0024 else 0025 a = angle; 0026 end 0027 end 0028 0029 AA = a; 0030 ann = angle; 0031 AA; 0032 ann; 0033 0034