### Python 计算 3D 空间三点垂直夹角方法 #### 向量表示与内积计算 对于给定的三维空间中的三个点 \( P_1(x_1, y_1, z_1) \)，\( P_2(x_2, y_2, z_2) \)，以及 \( P_3(x_3, y_3, z_3) \)，可以形成两个向量 \( \overrightarrow{P_1P_2} \) 和 \( \overrightarrow{P_1P_3} \)[^3]。 这两个向量可以通过坐标差来表达： \[ \overrightarrow{P_1P_2} = (x_2-x_1, y_2-y_1, z_2-z_1) \] \[ \overrightarrow{P_1P_3} = (x_3-x_1, y_3-y_1, z_3-z_1) \] 为了找到这两条线段所成的角度，可利用向量间的内积公式: \[ A·B=|A||B|\cos(\theta) \] 其中 \( |A| \) 表示向量 A 的模长，\(|B|\) 是 B 的模长，而 \( θ \) 则代表两者之间的角度[^1]。 因此， \[ \cos(\theta)=\frac{\overrightarrow{P_1P_2}\cdot\overrightarrow{P_1P_3}} {|\overrightarrow{P_1P_2}| |\overrightarrow{P_1P_3}| } \] 最后通过反余弦函数得到实际角度值。需要注意的是，在计算机程序里通常会返回弧度制的结果而不是度数形式；如果需要转换为度，则需乘以 \( \frac{180}{π} \). 下面是具体的Python代码实现： ```python import numpy as np def calculate_angle(p1, p2, p3): v1 = np.array([p2[0]-p1[0], p2[1]-p1[1], p2[2]-p1[2]]) v2 = np.array([p3[0]-p1[0], p3[1]-p1[1], p3[2]-p1[2]]) cos_theta = np.dot(v1, v2)/(np.linalg.norm(v1)*np.linalg.norm(v2)) theta_rad = np.arccos(np.clip(cos_theta, -1.0, 1.0)) # Ensure numerical stability return np.degrees(theta_rad) # Example usage: point1 = [0, 0, 0] point2 = [1, 0, 0] point3 = [0, 1, 0] angle_degrees = calculate_angle(point1, point2, point3) print(f"The angle between the vectors is: {angle_degrees:.2f} degrees") ``` 此代码片段定义了一个名为 `calculate_angle` 函数用于接收三个点作为输入参数并输出它们构成的平面内的夹角大小（单位：度）。这里还考虑到了浮点运算可能导致的小误差问题，采用了 `np.clip()` 来防止超出有效范围的情况发生。