### 计算两个平面的水平夹角和垂直夹角 在三维空间中，两个平面可以通过它们的法向量来定义。两个平面之间的夹角可以由它们法向量之间的夹角计算得出。假设两个平面的法向量分别为 \( \mathbf{n_1} = (a_1, b_1, c_1) \) 和 \( \mathbf{n_2} = (a_2, b_2, c_2) \)，则两个平面之间的夹角 \( \theta \) 可以通过以下公式计算： \[ \cos(\theta) = \frac{|\mathbf{n_1} \cdot \mathbf{n_2}|}{|\mathbf{n_1}| |\mathbf{n_2}|} \] 其中： - \( \mathbf{n_1} \cdot \mathbf{n_2} = a_1a_2 + b_1b_2 + c_1c_2 \) 是法向量的点积。 - \( |\mathbf{n_1}| = \sqrt{a_1^2 + b_1^2 + c_1^2} \) 和 \( |\mathbf{n_2}| = \sqrt{a_2^2 + b_2^2 + c_2^2} \) 分别是法向量的模。 水平夹角和垂直夹角的定义如下： - **水平夹角**：是指两个平面在水平面上投影后的夹角。 - **垂直夹角**：是指两个平面在垂直方向上的夹角。 以下是 Python 实现代码示例： ```python import numpy as np def calculate_plane_angles(n1, n2): # 将输入法向量转换为 NumPy 数组 n1 = np.array(n1) n2 = np.array(n2) # 计算法向量的模 magnitude_n1 = np.linalg.norm(n1) magnitude_n2 = np.linalg.norm(n2) # 计算法向量的点积 dot_product = np.dot(n1, n2) # 计算夹角的余弦值 cos_theta = abs(dot_product) / (magnitude_n1 * magnitude_n2) # 确保余弦值在有效范围内（避免浮点误差导致的问题） cos_theta = np.clip(cos_theta, -1.0, 1.0) # 计算夹角（弧度制） theta = np.arccos(cos_theta) # 转换为角度制 theta_degrees = np.degrees(theta) # 水平夹角：忽略 z 轴分量 n1_horizontal = n1[:2] n2_horizontal = n2[:2] magnitude_n1_horizontal = np.linalg.norm(n1_horizontal) magnitude_n2_horizontal = np.linalg.norm(n2_horizontal) dot_product_horizontal = np.dot(n1_horizontal, n2_horizontal) cos_theta_horizontal = abs(dot_product_horizontal) / (magnitude_n1_horizontal * magnitude_n2_horizontal) cos_theta_horizontal = np.clip(cos_theta_horizontal, -1.0, 1.0) theta_horizontal = np.arccos(cos_theta_horizontal) theta_horizontal_degrees = np.degrees(theta_horizontal) # 垂直夹角：仅考虑 z 轴分量 n1_vertical = n1[2] n2_vertical = n2[2] magnitude_n1_vertical = abs(n1_vertical) magnitude_n2_vertical = abs(n2_vertical) dot_product_vertical = n1_vertical * n2_vertical if magnitude_n1_vertical == 0 or magnitude_n2_vertical == 0: theta_vertical_degrees = 90 else: cos_theta_vertical = abs(dot_product_vertical) / (magnitude_n1_vertical * magnitude_n2_vertical) cos_theta_vertical = np.clip(cos_theta_vertical, -1.0, 1.0) theta_vertical = np.arccos(cos_theta_vertical) theta_vertical_degrees = np.degrees(theta_vertical) return theta_degrees, theta_horizontal_degrees, theta_vertical_degrees # 示例：计算两个平面的夹角、水平夹角和垂直夹角 n1 = [1, 2, 3] n2 = [4, 5, 6] angle, horizontal_angle, vertical_angle = calculate_plane_angles(n1, n2) print(f"两个平面的夹角: {angle:.2f} 度") print(f"两个平面的水平夹角: {horizontal_angle:.2f} 度") print(f"两个平面的垂直夹角: {vertical_angle:.2f} 度") ``` #### 结果解释 - **夹角**：通过 `np.arccos` 函数计算得到，表示两个平面之间的夹角。 - **水平夹角**：通过忽略法向量的 z 轴分量，仅考虑 x 和 y 分量的投影来计算。 - **垂直夹角**：通过仅考虑法向量的 z 轴分量来计算。 上述方法基于几何投影原理，确保了结果的正确性[^1]。 ### 注意事项 - 如果法向量为零向量，则无法定义夹角，需在实际应用中进行特殊处理。 - 使用 `np.clip` 避免由于浮点运算误差导致的余弦值超出 [-1, 1] 范围。