无预测寻的导弹对航路移动目标命中时间O(1)计算方法咨询

你描述的是经典的纯追踪（狗追）曲线问题，不存在覆盖任意目标轨迹的全局严格O(1)闭式解，但针对你给出的waypoint分段固定速率运动的目标场景，可以做到和航段数量线性相关的复杂度——单段匀速直航目标场景下是严格O(1)计算，完全可以替代小步长离散仿真，误差轻松控制在1ms以内。

核心求解思路 #

你的目标运动天然是分段直线的：两个相邻waypoint之间，目标速率固定、方向沿两点连线，属于匀速直航段，完全可以逐段解析计算，不需要逐tick迭代：

• 预处理所有航段：提前算好每一段的长度、目标在段内的速度矢量、段总飞行时长（段长/目标固定速率）
• 从初始状态开始逐段遍历：用单段纯追踪闭式公式，计算两个核心结果：① 导弹在当前段的时长范围内能不能追上目标；② 如果追不上，直接算出当前段结束时刻的导弹位置，作为下一段计算的初始状态
• 定位到命中发生的具体航段后，直接在段内解出精确命中时间，累加前面所有段的总时长，就是全局命中时间

单段匀速直航场景的O(1)闭式计算 #

先定义单段初始时刻（记为段内t=0）的参数：

• 导弹初始位置：M0 = (x_m0, y_m0)，导弹固定飞行速率：v_m
• 目标初始位置：T0 = (x_t0, y_t0)，段内目标固定速度矢量：v_t，目标速率：v_t_mag = norm(v_t)，段总时长：t_seg，段终点位置：T1 = T0 + v_t * t_seg
• 速度比：k = v_m / v_t_mag，你提到导弹必然命中，因此默认k>1

为了简化计算，先做坐标变换：把坐标系平移到t=0时的目标位置，再旋转坐标系让目标运动方向沿x轴正方向。变换后导弹的初始相对坐标为(x0, y0)，初始弹目距离r0 = sqrt(x0^2 + y0^2)，初始视线与目标运动方向的夹角θ0 = atan2(y0, -x0)。

这个场景下的纯追踪微分方程有成熟解析解，段内命中时间的计算公式为：

t_hit = (k * r0 * (1 - cosθ0)) / (v_m * (k*k - 1)) + (r0 * cosθ0) / (v_m + v_t_mag)

• 如果计算得到的t_hit落在[0, t_seg]区间内，说明命中就发生在当前段，直接返回结果即可
• 如果t_hit > t_seg，说明导弹在当前段追不上目标，代入纯追踪轨迹的位置闭式公式，算出t=t_seg时刻的导弹位置，带入下一段的初始条件继续计算即可

方案优势 #

• 不存在步长选择问题：单段计算仅涉及基础浮点运算，不管段长是100ms还是100s，计算量完全一致，没有离散仿真的线性化误差
• 整体复杂度为O(N)，N是目标waypoint的数量，和初始弹目距离、精度要求完全无关，比小步长仿真快几个数量级
• 精度仅受双精度浮点运算误差影响，不需要额外做碰撞检测的插值修正

边界处理说明 #

• 航段切换时目标速度方向突变的场景不需要特殊处理：直接把上一阶段结束时的导弹位置、新段的目标初始位置和速度矢量代入公式即可
• 若遇到k=1（导弹与目标速率相等）的特殊情况，只有初始时导弹位于目标正后方且视线与目标运动方向完全重合才能命中，否则无法追上，加个分支判定即可

内容的提问来源于stack exchange，提问作者Anıl Akkaya