如何通过编程计算两条相交直线的夹角？（求PHP/JS实现思路）

没问题，我来一步步帮你搞定这个问题~

要计算两条直线的夹角，首先得把直线转换成方向向量——用直线上的终点坐标减去起点坐标就能得到。对应你的场景：

• 直线LineA的方向向量a：P2A - P1A → (11-10, 10-9) = (1, 1)
• 直线LineB的方向向量b：P2B - P1B → (11-10, 9-10) = (1, -1)

有了向量之后，就可以代入你给出的向量夹角公式计算了。

公式：cosΘ = (a₁×b₁ + a₂×b₂) / (√(a₁²+a₂²) × √(b₁²+b₂²))
代入向量a(1,1)和b(1,-1)：

• 分子（点积）：1*1 + 1*(-1) = 0
• 分母（模长乘积）：√(1²+1²) * √(1²+(-1)²) = √2 * √2 = 2
• 计算得cosΘ = 0/2 = 0，对应的角度就是90度，和预期一致。

<?php
// 定义四个点的坐标
$p1a = ['x' => 10, 'y' => 9];
$p1b = ['x' => 10, 'y' => 10];
$p2a = ['x' => 11, 'y' => 10];
$p2b = ['x' => 11, 'y' => 9];

// 计算两条直线的方向向量
$vecA_x = $p2a['x'] - $p1a['x'];
$vecA_y = $p2a['y'] - $p1a['y'];

$vecB_x = $p2b['x'] - $p1b['x'];
$vecB_y = $p2b['y'] - $p1b['y'];

// 计算向量点积
$dotProduct = $vecA_x * $vecB_x + $vecA_y * $vecB_y;

// 计算两个向量的模长
$magA = sqrt($vecA_x ** 2 + $vecA_y ** 2);
$magB = sqrt($vecB_x ** 2 + $vecB_y ** 2);

$denominator = $magA * $magB;

// 处理边界情况（避免除以0）
if ($denominator == 0) {
    echo "无法计算夹角：其中一条直线的两个点重合了";
} else {
    // 计算cosθ，修正浮点精度问题
    $cosTheta = $dotProduct / $denominator;
    $cosTheta = max(min($cosTheta, 1), -1);
    
    // 弧度转角度
    $thetaRad = acos($cosTheta);
    $thetaDeg = rad2deg($thetaRad);
    
    echo "两直线夹角为：" . number_format($thetaDeg, 2) . "度"; // 输出90.00度
}
?>

// 定义四个点的坐标
const p1a = { x: 10, y: 9 };
const p1b = { x: 10, y: 10 };
const p2a = { x: 11, y: 10 };
const p2b = { x: 11, y: 9 };

// 计算方向向量
const vecA = {
  x: p2a.x - p1a.x,
  y: p2a.y - p1a.y
};
const vecB = {
  x: p2b.x - p1b.x,
  y: p2b.y - p1b.y
};

// 计算点积
const dotProduct = vecA.x * vecB.x + vecA.y * vecB.y;

// 计算模长
const magA = Math.sqrt(vecA.x ** 2 + vecA.y ** 2);
const magB = Math.sqrt(vecB.x ** 2 + vecB.y ** 2);

const denominator = magA * magB;

// 边界处理与计算
if (denominator === 0) {
  console.log("无法计算夹角：其中一条直线的两个点重合了");
} else {
  let cosTheta = dotProduct / denominator;
  // 修正浮点精度问题，避免acos返回NaN
  cosTheta = Math.max(Math.min(cosTheta, 1), -1);
  
  // 弧度转角度
  const thetaRad = Math.acos(cosTheta);
  const thetaDeg = thetaRad * (180 / Math.PI);
  
  console.log(`两直线夹角为：${thetaDeg.toFixed(2)}度`); // 输出90.00度
}

• 方向向量的选择：直线的方向向量可以是终点减起点，也可以是起点减终点，不影响最终夹角结果（因为cosθ和cos(-θ)相等）
• 浮点精度问题：计算机计算浮点数时可能出现微小误差，导致cosθ超出[-1,1]范围，所以需要用max/min修正，否则acos会返回NaN
• 边界情况处理：如果其中一条直线的两个点重合，向量模长为0，此时无法计算夹角，需要提前判断

内容的提问来源于stack exchange，提问作者whitebear