时间戳

首页 > 解决方案 > 如果向量相同，计算两个向量之间角度的方法会引发错误？

python - 如果向量相同，计算两个向量之间角度的方法会引发错误？

问题描述

我觉得这应该相当简单，但我一直盯着代码，现在很愚蠢。我有一个方法，它接受两个输入 3D 向量并以度为单位返回它们之间的角度。

方法如下：

def CalculateAngleBetweenVector(vector, vector2):

    dp = np.dot(vector, vector2)

    maga = math.sqrt((vector[0] ** 2) + vector[1] ** 2 + vector[2] ** 2)
    magb = math.sqrt((vector2[0] ** 2) + vector2[1] ** 2 + vector2[2] ** 2)
    magc = maga * magb

    dpmag = dp / magc

    angleindeg = ((math.acos(dpmag)) * 180) / math.pi

    return angleindeg

目前我有两个相同的向量来测试其他东西，并且当我这样做时已经意识到这个方法错误。这两个向量是：

[ 0.38154065 -0.38688849 -0.83949034]

我用来计算输入单位向量的方法是：

UnitVector = Vector / np.linalg.norm(Vector)

以防万一这里有问题。

非常感谢任何和所有帮助。

谢谢

编辑：

抱歉，输入的向量实际上是 (0.38154065, -0.38688849, -0.83949034) 我刚刚复制了终端打印行。该错误是数学域错误。

编辑2：

Error traceback:

Traceback (most recent call last):
  File "file path info that contains personal info Documents/com~apple~CloudDocs/WORK/Code1/CompareAnglesBetweenProteinAndMembane.py", line 261, in <module>
    angle = CalculateAngleBetweenVector(ProteinRotatedUV, MemRotatedUV)
  File "file path info that contains personal info Documents/com~apple~CloudDocs/WORK/Code1/CompareAnglesBetweenProteinAndMembane.py", line 167, in CalculateAngleBetweenVector
    angleindeg = ((math.acos(dpmag)) * 180) / math.pi
ValueError: math domain error

标签： pythonnumpymathvectorangle

解决方案

正如评论中所讨论的，问题源于浮点算术中的舍入错误（参数acos略微超出范围）。

我会考虑以下更改：

  1. 使用更简单的公式
angle = acos(dot(A,B) / (|A|* |B|))

  1. 如果 arccos 的参数略高于 1.0 或低于 -1.0，则将其四舍五入。如果参数在范围内，则不需要这样做。

  2. 规范化仅在向量接近 [0,0,0] 或具有非常大的条目时才有帮助 - 计算在数值上更稳定。如果您的应用程序都不是这种情况，请考虑删除规范化。

推荐阅读