找到两点之间角度的最快方法

Question

为了提高找到角度的正弦/余弦的速度，我建立了一个参考表，而不是即时计算它们。我也有同样的想法来找到从一个点到另一个点的角度。

我创建了一个包含 3600 个归一化向量的表（3600 / 10 = 十分之一度的精度）。每当我需要知道从一个点到下一个点的角度时，我都会查看表格以找到最佳匹配。但是，我担心这可能比使用 math.atan2() 慢。

这是我正在使用的代码：

创建向量表：

// vector to angle table
var vectorToAngleTable = new Array();
for (i = 0; i < 3600; i += 1) {
    vectorToAngleTable[i] = new Vector2();
    vectorToAngleTable[i] = RotatePoint(forwardVector, i / 10);
}

从两点求角度：

function NormalizeVector(vector) {
    var toReturn = vector;
    var dist = Math.sqrt(vector.x * vector.x + vector.y * vector.y);
    toReturn.x /= dist.x;
    toReturn.y /= dist.y;
    return toReturn;
}

function PointDirection(position, target) {
    var vector = target;
    var toReturn = 0;
    var smallest = 1.0;
    vector.x -= position.x;
    vector.y -= position.y;
    vector = NormalizeVector(vector);
    for (i = 0; i < 3600; i += 1) {
        if (PointDistance(vectorToAngleTable[i], vector) < smallest) {
            smalllest = PointDistance(vectorToAngleTable[i], vector);
            toReturn = i;
        }
    }
    return toReturn;
}

function PointDistance(point1, point2) {
    return Math.sqrt(((point2.x - point1.x) * (point2.x - point1.x)) + ((point2.y - point1.y) * (point2.y - point1.y)));
}

如您所见，我关心的是它要遍历的所有代码行，以及它要遍历的表中有多少条目。我很想知道找到角度的最快方法，不管是什么方法。

Answer 1

作为angle(v1, v2) = acos( (v1x * v2x + v1y * v2y) / (sqrt(v1x^2+v1y^2) * sqrt(v2x^2+v2y^2)) )，我们知道v2 = [1, 0]

var v = {x: 0, y: 1},
    angleRad = Math.acos( v.x / Math.sqrt(v.x*v.x + v.y*v.y) ),
    angleDeg = angleRad * 180 / Math.PI;

其中v 是向量[point2.x - point1.x , point2.y - point1.y]

---

编辑 - 我刚刚意识到你可能意味着将每个点视为一个向量，在这种情况下它会是

var v1 = {x: 0, y: 1}, v2 = {x: 1, y: 0},
    angleRad = Math.acos( (v1.x * v2.x + v1.y * v2.y) / ( Math.sqrt(v1.x*v1.x + v1.y*v1.y) * Math.sqrt(v2.x*v2.x + v2.y*v2.y) ) ),
    angleDeg = angleRad * 180 / Math.PI;

其中v1 是向量[point1.x , point1.y] 而v2 是[point2.x , point2.y]

---

编辑 2
如果您多次使用向量长度，要加快速度，请将其保存为例如v.length = ... 这样你就可以得到它而无需再次重新计算。 如果您知道每个向量都需要多次计算其角度，请使用我编写的第一种方法并将其缓存，即v.angle = ...。你可以然后你可以做v2.angle - v1.angle来找到两者之间的角度，等等。
即有

function Vector(x, y){
    this.x = x;
    this.y = y;
    this.length = Math.sqrt(x*x + y*y);
    this.angle = Math.acos( x / this.length );
}

---

jsperf 的预计算和在 3601 项目的数组中查找与直接使用 acos

Answer 2

这肯定会比调用atan2 更小，因为它是一个平方根，然后是对 3600 种可能性的线性搜索。相反，许多处理器直接实现 atan2——它是英特尔领域的 FPATAN。