求两向量超过180度的夹角

Question

楼主的公式错了：角度余弦 = 向量1 点乘 向量2 / (|向量1|*|向量2|);如果没错，那么函数里面的公式也写错了：float c=(float)sqr(dx1*dx1+dx2*dx2+dy1*dy1+dy2*dy2);
|向量1|*|向量2|=sqr(dx1*dx1+dy1*dy1)*sqr(dx2*dx2+dy2*dy2)
sqr(dx1*dx1+dx2*dx2+dy1*dy1+dy2*dy2)不等于sqr(dx1*dx1+dy1*dy1)*sqr(dx2*dx2+dy2*dy2)
至于求大于180度的问题我也遇到了，希望高手解答！

Answer 1

1 判断从向量1旋转到向量2是 要顺时针还是逆时针 （近距离旋转）
2 算角度
3 考虑是否2PI-角度 若我没有记错的话 上面是可行的
另外我算角度是arctg

Answer 2

通过叉积判断是否大于180度nA = normalize(A)
nB = normalize(B)c1 = cross(nA, nB)
c2 = cross(-nA, nB)d = dot(nA, nB)
r = acos(d)如果c1和c2相等则大于180度，否则c1和c2相反
根据这个条件对r进行相应调整即可

Answer 3

    算角度用arctg是不行的，因为arctg的值域是(-pi/2,pi/2)，这里没有等号。那么-pi/2和pi/2这些角度就永远求不出来了。
     

Answer 4

double ans =  atan2(Dir1*Dir2, Dir1%Dir2);    晕倒，我还是参考一个高手的代码

Answer 5

     兄弟，我对你这种判断方法表示怀疑。现在顺时针、逆时针的判断很好判断，可以通过
((dx1*dy2-dy1*dx2)>0) 就是顺时针，反之逆时针。只是夹角大于180度的情况不好判断。

Answer 6

    夹角大于180度的判断办法想到了，通过叉乘判断
||a*b|| = ||a||*||b||*sin&  (&为夹角)来判断。

Answer 7

struct Point 
{
    float x,y;
};float vectorMag(const Point &start, const Point &end)
{ 
 return sqrtf((start.x - end.x) * (start.x - end.x) +
               (start.y - end.y) * (start.y - end.y));
}float DotProduct(const Point &cen, const Point &first, const Point &second)
{
   dx1 = first.x - cen.x;
   dx2 = second.x - cen.x;
   dy1 = first.y - cen.y;
   dy2 = second.y - cen.y;   return dx1 * dx2 + dy1 * dy2; 
}//假定逆时针为正向
bool SignByCrossProduct(const Point &cen, const Point &first, const Point &second)
{
   dx1 = first.x - cen.x;
   dx2 = second.x - cen.x;
   dy1 = first.y - cen.y;
   dy2 = second.y - cen.y;
   
   if(dx1 * dy2 - dy1 * dx2 > 0) return true;
   return false;   
}//夹角范围[-PI, PI]
float Angle(Point cen,Point first,Point second)
{
   float dx1,dx2,dy1,dy2;
   float angle;
   
   float mag = vectorMag(cen, first) * vectorMag(cen, second);
   angle = (float)acos(DotProduct(cen, first, second) / mag);   bool sign = SignByCrossProduct(cen, first, second);
   
   if (!sign) angle = -angle;   return angle;
}

Answer 8

噢 没看清楼主是要求大于180度的夹角...//夹角范围[-PI, PI]
float Angle_180(Point cen,Point first,Point second)
{
   float dx1,dx2,dy1,dy2;
   float angle;
   
   float mag = vectorMag(cen, first) * vectorMag(cen, second);
   angle = (float)acos(DotProduct(cen, first, second) / mag);   bool sign = SignByCrossProduct(cen, first, second);
   
   if (sign) 
   {
     angle = 2 * PI - angle;
   }
   else
   {
     angle = angle - 2 * PI;
   }    return angle;
}

Answer 9

       if (sign) 
   {
     angle = 2 * PI - angle;
   }
   else
   {
     angle = angle - 2 * PI;
   }      兄弟，你这里判断有问题的。假设sign为true时，只有 angle 大于180度时,才需要angle = 2 * PI - angle; 否则应该直接返回angle.

Answer 10

楼主，我看你前面的要求是要求大于180的夹角。angle从acos算出来的时候都是在[0,PI]内的，你不是要求对应的那个大于180度的夹角么，所以用2*PI - angle

Answer 11

http://topic.csdn.net/u/20110706/12/61af0aaf-0df1-4a17-88c4-41dce79fd2ac.html
回个帖吧 我结贴

Answer 12

我搞工业cad的 技术相通的地方较多 加个qq  ？

Answer 13

楼主做出来了速度去写blog分享下

Answer 14

   二维: GDI、GDI+
   三维:  osg

Answer 15

http://bbs.csdn.net/topics/80026626 我觉得15楼正解！

Answer 16

方法一，计算方位角
a = atan2(x,y) //0~90度 的弧度表示
其他根据象限处理
两角相减，得到角度差 d
统一化为 0～360  或者 -180 ～180
方法2
求向量点积，可以求夹角出余弦。和360 的差是较大的角。