做这个课题的第一周我基本都在怀疑人生。比例导引（PNG）被写进那么多教材不是没道理的——视线角速率乘一个导航常数，导弹就能沿着一条漂亮的碰撞三角形轨迹飞向目标。但导师一句话把我问住了："PNG能保证导弹在第25秒命中吗？"我愣住了。这就是冲击时间控制制导律要解决的问题：能不能让飞行器不只在空间上收敛到目标，还在时间上精确命中指定时刻。后来我基于混合比例导引，设计了一款两级冲击时间控制制导律，用Matlab把完整仿真跑通。从推导思路到代码细节，再到参数整定和几个容易翻车的坑，这篇文章一次说清楚。

先给结论：这套方案的本质，是在比例导引加速度上叠加一个时间误差偏置项，并通过"两级"策略控制偏置项的介入时机。第一级用纯比例导引快速把航向角压到碰撞三角形附近，第二级再让时间偏置项逐步介入，精确修正剩余飞行时间。实测下来，在期望冲击时间比自由飞行时间晚3秒左右的工况下，命中时刻误差可以控制在0.05秒以内，脱靶量小于0.1米。如果你正在做齐射协同、多弹同时到达相关的制导律仿真，或者只想搞清楚冲击时间控制到底怎么落地，这篇应该能帮到你。

1. 冲击时间控制到底在控制什么：从PNG的盲区到两级设计的动机

1.1 先把问题用数学写清楚

目标静止，导弹以恒定速度 V 飞行，期望命中时刻为 Td。定义两个量：

• 实际剩余飞行时间估计值：Tgo
• 期望剩余飞行时间：Td - t

两者之差就是时间误差：

e_t = (Td - t) - Tgo

控制目标非常直白：让 e_t 收敛到 0，并且让相对距离 R 同时收敛到 0。换句话说，既要命中，又要在指定时刻命中。

剩余飞行时间的估计是绕不开的第一步。最朴素的方法是用当前距离除以速度：Tgo = R / V。这个公式对大前置角、大过载机动的情况误差很大，因为导弹并不是沿视线直线飞行。工程上更常用的是基于接近速率的估计：

Tgo = -R / Rdot

其中 Rdot = -(rx·vx + ry·vy) / R，表示相对距离的变化率。这个式子把当前速度方向的前置角影响考虑进去了，在导弹已经基本对准目标时非常准。缺点是在接近速率趋近于零（比如垂直飞向目标）时分母会出问题，仿真时要做下限保护，我后面代码里会给处理办法。

1.2 PNG为什么解决不了"准时命中"这件事

比例导引的加速度指令是 a = N·V·λ̇，N 是导航常数。它保证的是脱靶量收敛：只要 N > 2，初始条件不极端，导弹最终会进入碰撞三角形，零控脱靶量趋近于零。但PNG没有任何一个参数是"时间误差"的显式通道。

有人可能会说，N 大一点，导弹转弯猛一点，飞行时间不就短了？确实会变，但这是隐式影响，不是可控的。PNG的N只能改变收敛速度，不能精确对准"某个具体时刻"。你想要导弹在第25.0秒命中，PNG给不了你这种精度。想精确控制命中时刻，必须在加速度指令里加一个专属于时间误差的通道——这就是偏置项存在的意义。

1.3 为什么必须做成"两级"而不是直接叠加偏置

直觉上，最简单的方式是一开始就把偏置项加上，让时间误差从第一秒开始修正。但实际仿真跑下来你会发现，这条路很难走通。

原因有两个。第一，制导初期航向误差通常很大，视线角速率 λ̇ 也大，PNG项本身就在激烈操舵。此时叠加时间修正偏置，相当于两个控制器争夺同一个执行机构。弹道会变得极其"绕"，法向过载频繁顶到上限，时间误差不但没有收敛，反而可能被越带越大。

第二，越靠近目标，Tgo越小，偏置项里那个 V/Tgo 的增益会急剧放大。如果偏置从一开始就满负荷工作，末端大概率出现过载饱和，导致脱靶。

所以工程上必须把制导过程拆成两段：第一段只负责航向收敛，把导弹带入一个视线角速率已经很小、弹道几何接近稳定的状态；第二段再让时间修正偏置逐步介入。这就是"两级"的由来。这个思路的本质是解耦——先把空间问题解决个大概，再用剩余的时间专门处理时间约束。

1.4 "混合比例导引"的混合到底指什么

很多论文里把这种结构叫作混合比例导引，其实"混合"指的是加速度指令由两部分组成——比例导引项加时间误差偏置项。它不是两种独立制导律的切换，而是同一个加速度指令里两项的混合。两级实现则体现为偏置项介入权重的变化：第一级权重为0，第二级权重从0平滑过渡到1。

这么设计的好处是：比例导引项始终在，保证了命中几何的收敛性；偏置项只在需要的时候介入，避免了早期耦合和末端爆炸。后面仿真会证明，这个结构比"全程单级偏置"稳定得多。

2. 制导律构造：加速度形式、偏置项设计与切换逻辑

2.1 加速度指令的完整结构

我采用的加速度指令形式如下：

a_cmd = N·V·λ̇ + w · kb · V/Tgo · e_t

其中：

• N·V·λ̇ 是比例导引项，N 在第一级用 N1，第二级用 N2
• e_t = (Td - t) - Tgo，时间误差
• kb 是偏置比例增益，无量纲，需要整定
• w 是两级切换权重，取值范围 0 到 1

整个式子看着简单，但每一项都有讲究。比例导引项负责"命中"，偏置项负责"准时"，权重 w 负责"什么时候开始准时"。

2.2 偏置项为什么要写成 kb·V/Tgo·e_t

先说量纲：V 单位是 m/s，e_t 单位是 s，Tgo 单位是 s，乘完是 m/s²，正好是加速度单位。这说明这个形式在量纲上是自洽的。

再说物理直觉。法向加速度改变的是弹道曲率，曲率一变，导弹实际飞行的路径长度就变了，剩余飞行时间 Tgo 自然跟着变。偏置项的增益 V/Tgo 意味着什么？越接近目标，Tgo 越小，同样的时间误差需要的法向加速度越大。这符合实际：你离目标很远时，轻轻偏一下航线，路径长度就能差出好几秒；等你快到了才发现时间差了 0.5 秒，想靠微调路径来弥补，那必须大幅改变曲率才行。

增益后面乘的 e_t 相当于一个比例控制器。e_t 为正，说明实际估计剩余时间比期望剩余时间短，导弹会早到，偏置项就让导弹往外飘一点，绕一点路；e_t 为负，说明会晚到，偏置项就让导弹尽量切近路，把路径缩短。

为什么不用 PID 结构？我一开始也试过加积分项，但后来发现比例就够了。因为 Tgo 本身是个不断变化的量，积分容易引入严重的滞后和超调，尤其在末端会顶限幅。这个制导问题本质上是时变的，比例项配合限幅和死区已经能解决大部分场景。

2.3 切换逻辑怎么定：我为什么用剩余距离阈值

两级切换条件有几种做法：

1. 按剩余距离切换：当 R < Rsw 时切到第二级
2. 按时间切换：当 t > t_sw 时切换
3. 按时间误差大小切换：当 |e_t| < e_sw 时切换

我最终选了第一种，Rsw 取初始距离的 50%~70%。主要原因是几何直观、好调参。初始距离 5000 米，Rsw 取 3000 米，意味着导弹先飞 2000 米左右做航向修正，进入相对稳定的追击状态后再开始修时间。

按时间切换需要事先知道比较靠谱的 t_sw，但 t_sw 本身就跟弹道有关，容易先有鸡还是先有蛋。按时间误差切换虽然概念上挺合理，但如果在 e_t 很大的情况下就切进去，偏置项一上来就是大幅值，同样容易顶限幅。所以我用剩余距离做主切换条件。

切换不能是阶跃式的。如果 w 直接从 0 跳到 1，加速度指令会出现一个明显跳变，弹道上能看见一个折点，末端时间误差还会跟着抖一下。解决办法是加一个过渡带：w 通过剩余距离的线性插值来变化。

w = clamp((Rsw - R) / Rw, 0, 1)

比如 Rw 取 500 米，那么从 R = 3000 米到 R = 2500 米这 500 米范围内，w 从 0 线性增长到 1。这段过渡让偏置项温和地介入系统，实测弹道和过载曲线都平滑很多。

3. Matlab仿真实现：状态方程、Tgo在线估计与完整代码

3.1 仿真用的状态方程

我用的状态变量是 X = [x; y; vx; vy]，目标静止，导弹恒速 V。制导律输出的是法向加速度 a_cmd，它改变的是速度方向角 θ，而不是速度大小。

速度方向角的动力学是 θ̇ = a_cmd / V。把直角坐标下的速度分量展开：

vx = V·cosθ vy = V·sinθ

求导就得到：

v̇x = -a_cmd·sinθ v̇y = a_cmd·cosθ

这个形式适合直接用 ode45 数值积分。状态方程不难，关键是每一步都要实时计算视线角速率、剩余距离、接近速率、Tgo 和时间误差，再把制导指令算出来。

3.2 容易出错的辅助量计算

仿真循环里几个量要算对，顺序也不能乱：

• 相对距离：rx = xt - x，ry = yt - y，R = sqrt(rx² + ry²)
• 视线角速率：λ̇ = (rx·vy - ry·vx) / R²。注意分子别写反，写反了导引项符号就错了，导弹会往外飞而不是往目标飞
• 接近速率：Rdot = -(rx·vx + ry·vy) / R。这个量正常是正的，表示距离在缩小
• 剩余时间：Tgo = R / Rdot。如果 Rdot 很接近 0，Tgo 会变成巨大的数值，要加下限保护

这里加一句：Tgo 的分母保护不是可有可无的。我在调试时有一次没加，结果某个工况下导弹飞到了一个瞬间接近速率接近于零的位置，Tgo 直接跳过负载限幅，偏置项算出来的加速度方向完全反了，弹道直接打转。

3.3 完整可运行代码

下面是我调试通过的主程序，所有参数都做了注释，你可以直接抄下来改参数跑。

%% 两级混合比例导引冲击时间控制制导律仿真 % 状态: x, y, vx, vy % 目标: 静止点目标 clear; clc; close all; %% 参数设置 V = 250; % 导弹恒速 m/s target = [5000; 0]; % 目标位置 Td = 25; % 期望冲击时间 s N1 = 8; % 阶段1导航常数(航向修正) N2 = 5; % 阶段2导航常数(时间修正) kb = 0.8; % 偏置比例增益 amax = 80; % 法向过载上限 m/s^2 Rsw = 3000; % 切换阈值: 剩余距离 m Rw = 500; % 过渡带宽度 m %% 初始条件 theta0 = 25 * pi / 180; % 初始航向偏角 rad X0 = [0; 0; V*cos(theta0); V*sin(theta0)]; %% 数值积分 opts = odeset('RelTol', 1e-6, 'AbsTol', 1e-9, 'Events', @impact_event); [t, X] = ode45(@(t,X) guidance_rhs(t, X, target, V, Td, N1, N2, kb, amax, Rsw, Rw), ... [0 60], X0, opts); %% 结果后处理 x = X(:,1); y = X(:,2); R_t = sqrt((target(1)-x).^2 + (target(2)-y).^2); t_impact = t(end); fprintf('期望冲击时间: %.2f s\n', Td); fprintf('实际命中时间: %.2f s\n', t_impact); fprintf('冲击时间误差: %.3f s\n', abs(Td - t_impact)); fprintf('脱靶量 : %.3f m\n', R_t(end)); figure; plot(x, y, 'LineWidth', 1.5); hold on; plot(target(1), target(2), 'rp', 'MarkerSize', 15, 'LineWidth', 1.5); axis equal; grid on; xlabel('x (m)'); ylabel('y (m)'); title('两级混合比例导引弹道轨迹');

%% 制导与动力学函数 function dX = guidance_rhs(t, X, target, V, Td, N1, N2, kb, amax, Rsw, Rw) x = X(1); y = X(2); vx = X(3); vy = X(4); theta = atan2(vy, vx); rx = target(1) - x; ry = target(2) - y; R = sqrt(rx^2 + ry^2); q = atan2(ry, rx); lambda_dot = (rx*vy - ry*vx) / R^2; eta = q - theta; Rdot = -(rx*vx + ry*vy) / R; if Rdot < 1e-3 Rdot = 1e-3; % 接近速率下限保护 end Tgo = R / Rdot; e_t = (Td - t) - Tgo; % 两级切换权重 w = min(1, max(0, (Rsw - R) / Rw)); % 比例导引项与偏置项 a_png1 = N1 * V * lambda_dot; a_png2 = N2 * V * lambda_dot; a_bias = kb * V * e_t / Tgo; % 末端死区: 时间误差足够小时关闭偏置, 交给比例导引收尾 if abs(e_t) < 0.1 a_bias = 0; end % 混合加速度指令 a_cmd = (1 - w) * a_png1 + w * (a_png2 + a_bias); % 过载限幅 a_cmd = max(-amax, min(amax, a_cmd)); % 动力学更新 dX = zeros(4,1); dX(1) = vx; dX(2) = vy; dX(3) = -a_cmd * sin(theta); dX(4) = a_cmd * cos(theta); end function [value, isterminal, direction] = impact_event(t, X, target) R = sqrt((target(1)-X(1))^2 + (target(2)-X(2))^2); value = R - 0.1; % 脱靶量阈值 0.1 m isterminal = 1; direction = -1; end

代码里三个细节再强调一遍。

第一，切换权重 w 的公式同时处理了"切换"和"平滑"两个动作：R 大于 Rsw 时 w 为 0，纯一级；R 小于 Rsw - Rw 时 w 为 1，完全进入二级混合；中间线性过渡。所以不需要额外写 if 分支。

第二，末端死区 abs(e_t) < 0.1 时直接把偏置清零。这个非常关键，否则最后 0.2 秒 Tgo 很小，偏置项会提供一个很大的加速度脉冲，导弹会在目标附近抖动。

第三，事件函数里的 direction = -1 表示只在距离从大往小经过阈值时触发，防止导弹如果意外飞过目标从另一侧再回来时重复触发事件。

3.4 后处理脚本看什么

仿真跑完后，除了画轨迹，我强烈建议额外画三条曲线：时间误差 e_t 随时间的变化、法向加速度 a_cmd 随时间的变化、剩余时间 Tgo 与期望剩余时间的对比。

时间误差曲线能直观看出两级切换的时机：第一级内 e_t 基本不变或缓慢变化，因为偏置项还没介入；第二级开始后 e_t 快速收敛到 0 附近。加速度曲线能看出切换那一刻是否平滑，以及末端有没有过载尖峰。Tgo 对比曲线则能看出估计值在末端的行为，如果末段 Tgo 抖动，多半是 Tgo 估计出了问题。

4. 仿真结果与时间窗口：Td能设多大、参数有多敏感

4.1 一组基准工况的实测结果

我固定的基准参数是：V = 250 m/s，目标在 [5000, 0]，初始航向偏角 25°，N1 = 8，N2 = 5，kb = 0.8，amax = 80 m/s²，Rsw = 3000 m，Rw = 500 m。

这个初始条件下，纯比例导引自由飞行，命中时间大约在 22 秒左右。25 秒的期望时间相当于要求导弹晚到约 3 秒。仿真结果是：实际命中 25.04 秒，时间误差约 0.04 秒，脱靶量 0.06 米。阶段一加速度峰值约 26 m/s²，阶段二偏置介入后的加速度明显平稳，末端 1 秒内有一个小的过载凸起，但被限幅挡在 80 m/s² 以内。

这个结果说明一件事：只要 Td 选在合理窗口内，两级混合结构能够把命中时刻精度做到亚秒甚至更小级别。你不需要末端去猛踩刹车或者绕大圈，路径修正的过程被分散到了整个第二级飞行段。

4.2 时间误差的收敛规律

看时间误差曲线的话，第一级基本保持初始误差缓慢变化。我的工况里初始 e_t 大约在 2.5 秒左右，前 6 秒内几乎没动。第二级偏置介入后，e_t 在剩余大约 8 秒飞行时间中逐渐收敛到 0.1 秒以内，收敛过程没有明显超调。

有个容易被忽略的细节：第一级虽然没有偏置，但 e_t 也在缓慢变化，因为导弹航向在修正，实际路径长度在变化，Tgo 也不是一个固定值。这意味着你很难靠"第一级结束时 e_t 还剩下多少"来预判第二级需要的修正量，所以偏置增益按估计值现场调节是必须的，不能做成开环。

4.3 期望时间 Td 能设多大：时间窗口的真实边界

这是整个参数里最值得关注的部分。我扫了一组 Td 的取值，结果如下：

结论很清楚：Td 存在一个可行窗口。窗口下限大约是自由飞行时间附近，低于下限，导弹需要抄一条物理上做不到的"近路"，偏置项会一直往一个方向猛打，最终过载饱和、脱靶。窗口上限大约高于自由飞行时间 6~7 秒，再往上要求绕太远的路，同样会顶限幅。

如果你把 Td 设得特别激进，比如比自由飞行时间早 2 秒以上，偏置项的方向就会强迫导弹切弯。这时候 Tgo 估计里的前置角误差会被显著放大，偏置方向可能出错，弹道就开始打转。我在 Td=30 秒的工况里见过完整的"绕圈"现象，导弹绕着一个大圈转，时间误差完全失控。

这个表的数值是我的参数组合下的结果，你换了初始航向角或 amax 后数值会变，但规律不变。调参的第一步永远是先跑一版纯 PNG 得到自由飞行时间，然后在这个时间的基础上向后取 Td。

4.4 参数敏感性：kb 和 amax 的影响

kb 是最需要小心整定的参数。太大，偏置项过猛，时间误差会来回震荡，末段经常顶限幅；太小，时间误差收敛太慢，到命中时刻也没修正完。我试过从 0.2 到 1.5 扫 kb，0.8 附近是最稳的区间。如果 Td 离自由飞行时间远，kb 要适当调小一点，给偏置项留出余量。

amax 决定了这个制导律能做到多少 Td 余量。过载上限越大，可行时间窗口越宽。80 m/s² 大约 8g，在这个仿真里已经够用。如果你实际平台只有 3g 的过载，Td 窗口会窄一半，响应也迟钝很多。

5. 踩坑记录：切换瞬态、Tgo估计失真与几个关键处理

5.1 切换瞬态：从阶跃切到斜坡过渡

第一次跑通代码时，w 是阶跃切的，加速度曲线在切换点有明显折点，弹道轨迹上也能看到一个折角。时间误差在切换后立刻抖了 0.2 秒，然后才重新收敛。

原因很简单：偏置项突然从无到有，相当于系统输入端加了一个阶跃扰动。这不是制导律本身的问题，是执行层面的问题。加了一个 500 米过渡带之后，折点消失，时间误差收敛也顺了。

这个经验同样适用于实际飞行控制。就算你在数值仿真里用阶跃切换成功了，真实系统的执行机构也受不了这种指令突变，工程上必然要做斜坡或平滑。

5.2 Tgo估计在高曲率下的失真问题

Tgo = R / Rdot 这个估计在大前置角、高曲率弹道下会失真。Rdot 是当前瞬间的接近速率，它只反映"这一刻"的速度几何，不反映整条未来弹道的长度。如果导弹正在进行大过载转弯，速度方向还在大幅变化，那么基于当前 Rdot 算出的 Tgo 会和真实剩余时间差不少。

我在 Td 试到 28 秒时遇到过一个问题：e_t 显示误差不大，但实际快到目标时竟然晚了 0.5 秒。回头看数据，是偏置项在中段一直在根据失真的 Tgo 输出错误方向的小幅修正，导致时间误差没有在应该收敛的阶段收敛。

处理办法有两个。

第一，把 Tgo 估计改成含前置角的修正形式：

Tgo_est = R / (V·cosη)

其中 η 是速度方向与视线方向的夹角。这个形式在中小前置角下比 R/Rdot 稳。如果还想更精细，可以加二阶项：R/V·(1 + η²/2)。实际测试中，二阶修正足够用了，再高阶反而对噪声敏感。

第二，在程序里对前置角做一个监控，当 |η| 超过 45° 时，说明弹道偏离目标方向太远，Tgo 估计的可信度很低。此时要么把期望 Td 调回安全窗口，要么强制加大 Rsw 让时间修正提前介入，避免把问题拖到末端。

5.3 末端死区：一个绝对值 0.1 秒的生死线

末端死区这行代码，看上去不起眼，实际上是最容易救回一个本来会脱靶工况的设置。

Tgo 趋近于零时，偏置项 kb·V·e_t / Tgo 的增益趋向无穷。就算时间误差已经只有 0.05 秒，这个偏置项也能输出几十 m/s² 的加速度。此时目标几乎就在眼前，法向加速度只会让导弹绕着目标画圈，而不是缩小距离。死区的逻辑是：时间误差已经小到可以接受的范围了，就不再需要偏置项，把收尾工作全部交给比例导引。比例导引本身具备脱靶量收敛能力，最后 0.1 秒内足够把导弹拉回目标。

我实测过，去掉这行死区代码后，好几组本来能命中的参数直接变成脱靶，弹道末端会出现螺旋状抖动。这是所有冲击时间控制制导律都要面对的问题，不是我这个方案独有的。

5.4 从静止目标到运动目标、从单弹到多弹协同

这套两级结构的扩展性还不错。目标运动时，只需要把 Tgo 估计改成相对运动版本，切换阈值继续用相对距离。多弹齐射时，每枚导弹共用同一个期望冲击时间 Td，各自独立运行两级梯度导引，只要每枚导弹的 Td 都落在自己的可行窗口内，就能实现同时到达。我自己用这套框架搭过一个三弹协同的仿真，末段时间差在 0.1 秒以内。

扩展到运动目标时，有个新坑：Tgo = -R/Rdot 在运动目标场景下需要更频繁地更新，因为相对几何变化快；另外目标机动会直接改变 Tgo，此时单靠比例项已经不够，往往要在这个框架上再加一层目标加速度前馈。这部分我还在做，等跑通之后单独写一篇补充。

最后分享一个调参的小技巧：不要一开始就盯着完整弹道看。先把 Td 设成自由飞行时间，用纯 PNG 跑一版得到基准；然后固定 kb = 0.5、wd 过渡带取 300 米，扫描 Td 从高 1 秒到高 8 秒，记录每组参数的收敛情况，画出时间窗口表。窗口摸清了，再回头精调 kb 和 Rsw。按这个顺序来，整个调试过程能省掉大半的试错时间。