2D移动目标射击拦截求解 二次方程返回无穷值问题排查

用二次方程求解运动目标射击提前量的思路是正确的，代码存在多处错误导致计算始终返回无穷值：

• 目标速度取值完全错误（核心问题）
  代码中double victimSpeed = victim.Position.Length();取的是目标坐标到世界坐标系原点的距离，根本不是目标的移动速度。目标速度的大小应该取速度向量的模长：victim.LinearVelocity.Length()。
  只要目标离原点有一定距离，哪怕完全静止（LinearVelocity为0），算出的victimSpeed都会远大于子弹速度，直接导致二次方程判别式为负，返回无穷值。
• 相对向量方向存在隐患
  定义的vectorFromVictim = bullet.Position - victim.Position是目标指向发射点的向量，虽然计算向量长度、点乘时不会影响系数数值，但后续如果扩展逻辑（比如加入A的速度、子弹初速度叠加）很容易出现符号错误，建议统一使用发射点指向目标的初始向量deltaP = victim.Position - bullet.Position作为计算基准。
• 未考虑发射方A的移动速度
  需求中明确A沿特定路径移动，如果游戏物理规则是子弹发射后继承A的移动速度（绝大多数2D俯视角射击的默认逻辑），当前方程完全没有叠加A的瞬时速度，会导致命中位置计算偏差；如果子弹发射后速度完全独立于A的运动，可以忽略此项。
• 二次方程根的筛选逻辑有漏洞
  当子弹速率小于目标移动速率时，二次项系数a为负，函数开口向下，当前优先取较小根的逻辑会错误丢弃有效拦截解；同时没有显式过滤t<=0的无效根（负t代表碰撞发生在子弹发射前，没有实际意义）；当子弹速度与目标速度大小相等时a=0，会出现除以0的NaN错误，直接返回无穷。

拦截问题的核心约束是：子弹飞行t时间后的位置，与目标移动t时间后的位置完全重合。
设发射瞬间：

• 发射点A（子弹初始位置）坐标为P_a，子弹飞行速率为bulletSpeed，最大存活时间为bulletLifetime
• 目标B初始坐标为P_b，移动速度向量为v_v
  碰撞满足等式：P_a + v_b * t = P_b + v_v * t，其中v_b是子弹速度向量，模长固定为bulletSpeed。
  对等式两边取模长平方，整理为标准二次方程a*t² + b*t + c = 0，对应系数为：

Vector3 deltaP = victim.Position - bullet.Position; // 发射点指向目标的初始向量
float initialDist = deltaP.Length();
float victimSpeed = victim.LinearVelocity.Length();
// 如果子弹继承A的速度，需要将目标相对速度改为 victim.LinearVelocity - a.CurrentVelocity
double a = bulletSpeed * bulletSpeed - victimSpeed * victimSpeed;
double b = 2 * Vector3.Dot(deltaP, victim.LinearVelocity);
double c = -initialDist * initialDist;

求解得到有效t之后，拦截位置为victim.Position + victim.LinearVelocity * t，把A当前位置指向该拦截位置的单位向量作为发射方向即可。

double hitT = float.PositiveInfinity;
// 处理子弹与目标速度相等的一次方程边界情况
if (Math.Abs(a) < 0.0001f)
{
    double t = -c / b;
    if (t > 0.001f && t <= bulletLifetime) hitT = (float)t;
    return hitT;
}

double discriminant = b * b - 4 * a * c;
if (discriminant < 0)
{
    return float.PositiveInfinity; // 无可行拦截解
}

double sqrtD = Math.Sqrt(discriminant);
double t1 = (-b + sqrtD) / (2 * a);
double t2 = (-b - sqrtD) / (2 * a);

// 筛选所有t>0（碰撞在发射后发生）且t在子弹生命周期内的解，取最小t（最快命中）
if (t1 > 0.001f && t1 <= bulletLifetime)
{
    hitT = (float)t1;
}
if (t2 > 0.001f && t2 <= bulletLifetime)
{
    if (t2 < hitT)
    {
        hitT = (float)t2;
    }
}

return hitT;

补充：如果游戏中A是静止单位，只需要修正目标速度的取值、替换根筛选逻辑，就能解决返回无穷值的问题。

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