企业官网建设流程全解析

简介：三维比例导引法仿真针对飞行器制导控制领域常见需求，面向导弹制导律设计、航向控制及仿真验证方向的初学者与研究人员，提供一套可直接运行的Matlab实现，支持对静止目标和机动目标进行攻击仿真。压缩包体积仅5KB，共包含6个文件，其中5个为.m脚本，覆盖龙格库塔积分器、比例导引指令计算、静止目标模型、机动目标模型等核心功能，另有1个txt说明文档梳理代码结构与使用要点。目前已吸引712人学习，适合作为三维制导律课程设计、毕设仿真或算法验证的参考起点。借助这份代码，读者可以快速理解三维比例导引法的几何关系与算法流程，通过修改目标运动参数或导引系数，观察不同场景下的弹道特性和脱靶量，进而加深对制导律设计原理的掌握。

1. 三维比例导引法仿真：先看一版能跑的弹道，再谈命中精度

三维比例导引法仿真，听起来是制导领域的老话题，但真正把代码跑起来会发现：二维平面里比例导引只是一条视线角速率公式走到底，三维里却要处理位置矢量、速度矢量、视线旋转方向之间的配合，坐标系稍错一步，弹道就直接螺旋发散。这份资源把五段 MATLAB 代码和一个说明文档打包在一起，四阶龙格库塔积分、比例导引律、静止目标、正弦机动目标和方波机动目标均有对应实现，适合正在做末制导算法验证、需要快速拿到一组可改可跑的弹道仿真基线的人。本文按文件拆解到参数调优的顺序，把每条能落地的操作和值得注意的坑位过一遍。

2. 仿真包文件拆解：五个 m 文件和一个说明文档怎么配合

2.1 文件分工：哪段负责积分、哪段负责制导、哪段负责目标运动

拿到资源后先别急着双击运行，把几个文件的职责对一遍。bili3dnew.m 是入口脚本，里面集中了仿真初始条件、总时长、积分步长、导航比这些用户级参数；proportional3dnew.m 是核心算法文件，负责在每个积分步里根据弹目相对位置和相对速度生成制导指令加速度；rungekutta3d.m 是四阶龙格库塔积分器，把导弹和目标的运动微分方程推进一步；sinmotorizeda3dnew.m 和 squaremotorizeda3dnew.m 是两类机动目标模型，一个给正弦机动，一个给方波机动；代码说明文档.txt 里写明了调用顺序、输出量和典型参数。理解了这份分工，后面替换目标模型或者改导引律就非常顺：哪个文件动哪几行，心里有数。

从文件命名还能看出一个重要特征：这份代码走的是模块化路线，弹目相对运动通常被整理成 12 维状态向量，前 6 维是导弹的位置和速度，后 6 维是目标的位置和速度。目标模型文件只负责给出目标加速度，不关心导弹怎么制导；比例导引文件只读取弹目相对状态，不关心目标内部是正弦还是方波。这种划分是弹道仿真里的常见做法，本质上是把"目标怎么动"和"导弹怎么追"解耦，改目标模型时制导律代码可以一行不动。

2.2 初始化与主流程：bili3dnew.m 里要改哪些参数

bili3dnew.m 是主程序，用户在仿真里要调整的东西基本都集中在文件靠前的参数区。我复现这类资源时习惯先把这些参数过一遍，再跑第一次：

% bili3dnew.m 起始位置的参数区（示例，字段名以包内实际为准） dt = 0.01; % 固定积分步长，单位秒 t_end = 15.0; % 弹道最大仿真时间，单位秒 N_nav = 4.0; % 比例导引导航比，范围通常取 3~5 a_max = 49.0; % 导弹可用过载上限，单位 m/s^2，对应 5g % 初始状态：位置与速度均为 3 维列向量 r_m0 = [0; 0; 0]; % 导弹初始位置 v_m0 = [320; 0; 0]; % 导弹初始速度，单位 m/s r_t0 = [2600; 500; 300]; % 目标初始位置 v_t0 = [0; 0; 0]; % 目标初始速度，静止目标写 0 向量

这段代码的作用是把仿真时间、步长、过载限制和初始弹目几何全部收敛到同一个文件里，后续所有弹道对比都以这组参数为基线。需要强调的是：整套代码按 SI 单位制走，长度是米、速度是米每秒、加速度是米每二次方秒，单位混用是后续最常翻车的地方。参数区的变量名不一定和包内完全一致，但含义是固定的，先核对这些名字，就能快速定位到 bili3dnew.m 的核心设置。

主流程部分，bili3dnew.m 里通常是循环推进：每个时间步先调目标模型函数拿到目标当前加速度，再调 proportional3dnew.m 生成制导指令，把指令写进导弹状态导数，最后交给 rungekutta3d.m 做一步积分。这一层不需要每次都改，但理解了它，才知道"把目标从静止切到机动"具体是在哪个环节下手。仿真推进过程中，代码一般会把每个时间步的状态量存进数组，仿真结束后在工作区里能看到完整的弹道时间序列。

2.3 四阶龙格库塔积分器：rungekutta3d.m 的输入输出

rungekutta3d.m 本质上是标准 RK4 的封装，接口通常是传入当前时刻、当前状态向量、状态导数函数和积分步长，返回推进一步之后的状态向量。很多仿真项目喜欢用变步长 ode45，但弹目拦截这类需要固定时间间隔输出曲线的场景，固定步长 RK4 反而更好用：输出数组尺寸干净，脱靶量统计、视线角速率曲线做差分都方便。RK4 的核心是取四个斜率再按 1/6 权重合成，单步精度和稳定性对大多数三维拦截场景是够用的。

% rungekutta3d.m 的典型接口结构 function x_next = rungekutta3d(f, t, x, h) % f: 状态导数函数，输入(t, x)输出 dx/dt % t: 当前时刻；x: 当前状态向量；h: 积分步长 k1 = f(t, x); k2 = f(t + 0.5*h, x + 0.5*h*k1); k3 = f(t + 0.5*h, x + 0.5*h*k2); k4 = f(t + h, x + h*k3); x_next = x + (h/6) * (k1 + 2*k2 + 2*k3 + k4); end

这段代码是纯粹的数学模板，没有业务逻辑，但工程上真正的坑恰好不在模板本身，而在于它默认输入的一切都是对的。我看这类三维仿真代码时有个习惯，先把 k1 到 k4 四个中间量里有没有 NaN 或极大值检查一遍，如果有问题，优先怀疑状态向量的组织方式，而不是积分器本身。这里的导数函数 f 必须保持脚本风格一致，返回的列向量尺寸和 x 完全一致。如果导弹状态和目标状态的定义顺序在两个文件里不一致，RK4 会把不同物理量的导数叠加到错误位置，弹道虽然还能画出来，但物理含义已经彻底错了。

2.4 从零复现的调用顺序：先跑谁、后看谁

代码说明文档.txt 里一般会写一行"直接运行 bili3dnew.m"。但第一次跑我建议做一个干净的基线对照：先保持静止目标，把导航比 N 设成 3，看一条平滑的拦截弹道；再把目标模型切成正弦机动，导航比加到 4，对比脱靶量变化。这样能确认仿真链路是通的，之后才去改速度、位置这些场景参数。如果一上来同时改目标模型、过载限幅和初始位置，出了问题根本说不清是哪个环节引起的。

# 推荐的复现顺序（在 MATLAB 命令行里依次执行） # 1. 用编辑器打开 bili3dnew.m # 2. 确认 dt=0.01, t_end=10 # 3. 目标模型函数名处填写 sinmotorizeda3dnew 或 squaremotorizeda3dnew # 4. 运行 bili3dnew # 5. 在工作区检查输出量，画弹目距离曲线 plot(t_out, dist_out); xlabel('t (s)'); ylabel('miss distance (m)');

第一次跑完"静止目标加 N=3"的基线场景，看到的效果应该是这样的：导弹从初始位置出发，前段基本平稳飞行，进入末段时轨迹出现一个平滑弯曲，把指向目标的初始偏差逐步修正回来。弹目距离曲线应单调下降并在终端附近触及最小值。如果发现轨迹中途出现大幅折线，或者距离曲线不是单调下降而是先降后升，先回头检查 t_end 是否足够长、积分步长是否过大，再考虑制导指令是否在初始阶段就饱和。这个基线场景是整个调参往下走的前提，之后再换机动目标模型时，手里已经有了一条可靠的对照曲线。

3. 比例导引律实现：从视线角速率到过载指令

3.1 比例导引的本质：视线旋转角速度的 N 倍放大

经典比例导引在二维里一句话就能讲完：导弹加速度指令与视线角速度成正比，比例系数就是导航比 N。三维里的核心结论没有变，视线角速率从二维平面的标量变成了空间中的角速度矢量，制导指令同样取 N 倍的角速度乘上特征速度，再作用在垂直于视线的方向上。这个逻辑在工程里长盛不衰，是因为它不需要目标预知模型，只依赖当前弹目相对运动：视线在转，导弹就加过载把视线转率压下去。

proportional3dnew.m 就是沿着这个思路实现。实现上有几种表达方式，常见做法是用向量叉积直接计算视线角速率，省去显式定义俯仰和偏航两个视线角的步骤，在任意三维几何下都成立，也避开角度跳变问题。这份资源里的代码应该走的是这类纯向量路线：拿到弹目相对位置和相对速度，叉积除以距离平方得到视线角速度矢量，再乘导航比和接近速度得到指令。整个过程不需要三角函数分支，写起来短，跑起来也稳。

经典结论是，对静止目标或匀速直线目标，导航比略大于 2 就能保证运动收敛；对机动目标，导航比通常要取到 4 以上才有足够的跟踪刚性。这个资源支持静止目标和两类机动目标，正好把这两个区间都覆盖到了。

3.2 三维视线坐标系：弹目相对量与视线角速率的向量算法

三维仿真里最容易绕晕的是坐标系，但比例导引本质只需要惯性系下的弹目相对信息，不需要把运动全部转进弹体坐标系。设目标相对导弹的位置矢量为 R = r_t - r_m，相对速度矢量为 V = v_t - v_m，弹目连线方向单位向量为 R_hat。视线角速度矢量的标准表达式是 R 与 V 的叉积除以 R 模的平方，因为叉积同时包含了视线在两个正交方向上的旋转信息，除以距离平方后量纲就是角速度。

% proportional3dnew.m 中视线角速率与指令计算的核心逻辑 function a_cmd = proportional3dnew(r_rel, v_rel, N, a_max) % 输入 r_rel: 目标相对导弹的位置矢量 [x;y;z], 单位 m % v_rel: 目标相对导弹的速度矢量 [vx;vy;vz], 单位 m/s % 输出 a_cmd: 导弹在惯性系下的指令加速度, 单位 m/s^2 r_norm = norm(r_rel); if r_norm < 1e-4 a_cmd = [0;0;0]; % 认为已进入命中区，停止生成新指令 return; end omega_los = cross(r_rel, v_rel) / (r_norm^2); % 视线角速率矢量，单位 rad/s v_close = -(r_rel' * v_rel) / r_norm; % 接近速率，接近时为正值 a_cmd = N * v_close * omega_los; % 比例导引指令，量纲 m/s^2 a_cmd_l = norm(a_cmd); if a_cmd_l > a_max a_cmd = a_cmd * (a_max / a_cmd_l); % 过载限幅，保持指令方向不变 end end

这段代码是比例导引函数的核心骨架，逻辑上做了三件事：先算视线角速率，再算接近速率，最后生成指令并限幅。重点看三个细节。第一，r_rel 的方向约定要全局统一，这里采用"目标相对导弹"，若采用相反约定，omega_los 会反向，N 前面的符号也要跟着变，三维里这种符号错误不像二维那样容易一眼看出来。第二，v_close 取的是负的视线方向速度内积，接近时为正，指令方向才会在物理上表现为"压视线角速率"。第三，限幅作用于合成后的指令矢量，不是分别对三个通道做饱和，分开限幅会改变指令方向，这是很多实现容易写错的地方。

接近速率的含义值得单独说一句：R_hat 与 V 的内积为正，说明目标沿视线方向远离，接近速率为负；取负号后接近时为正值。工程上有另一类实现用相对速度模值代替接近速率，那就是 PPN 的表达形式，和这里的 TPN 在末端大角度场景下略有差异。这份资源用接近速率的版本，数学上更贴合 TPN 的定义，末端稳定性会好一点。

3.3 制导指令限幅与终端条件：别让导弹过载爆掉

限幅环节在代码量上只有两行，却往往是弹道是否物理可用的分水岭。不设限幅时，比例导引在大初始误差或强机动目标条件下，指令加速度可以冲到几百米每二次方秒，仿真曲线依然"漂亮"地命中目标，但这不是导弹能飞出来的轨迹，拿到六自由度平台里根本没法复现。因此 a_max 必须根据飞行器可用过载来设，典型值从 3g 到 10g 不等。对这个资源里的弹道级仿真，a_max = 49 米每二次方秒即 5g 是比较均衡的起点，既能体现拦截机动目标需要的过载，又不会让曲线失真。

终端条件同样重要。比例导引在距离很小时，视线角速率公式里的 r² 项会快速放大，产生指令数值尖峰。因此距离小于某个阈值后停用制导指令是常见保护做法，否则末端一小段轨迹反而出现不应有的弯曲。主脚本里通常对两类终端条件做或运算：弹目距离小于阈值，或者仿真时间到达 t_end，谁先到谁结束。

% 主循环终端判定的常见写法（示意，与包内 bili3dnew.m 对应） t = 0; idx = 1; while t < t_end % 计算导弹与目标状态，得到当前弹目距离 dist if dist < 1.0 % 脱靶量阈值，命中即停 break; end % 推进一个步长，记录数组 t = t + dt; idx = idx + 1; end miss_dist = dist; % 仿真结束后的弹目距离即脱靶量

这段循环体现了一个工程习惯：脱靶量的取值不在积分器输出的最后一刻，而是在整个时间序列的最小值。真实场景里弹道可能因为过载限幅出现"穿过目标"的情况，弹目距离先降到 0.5 米又上升到 30 米，用最后一步距离作评价会严重高估性能。统计口径放到整个距离序列的最小值上，是这类仿真评价的默认动作。

4. 目标模型的三类选择：静止、正弦机动与方波机动

4.1 静止目标：视线角速率天然平稳，适合做链路基线

目标静止时，目标速度恒为零，状态在仿真过程中完全不动。这类场景的价值不在"目标真实"，而在它是最干净的链路自检手段。导弹从初始位置出发，只要接近速度为正、引导方向合理，视线角速率会随距离减小逐渐被压小，脱靶量应单调收敛到很小值。如果静止目标下弹道都走形、脱靶量大于几十米，问题几乎一定出在制导律实现、积分步长或状态定义上，而不是目标模型本身。

从参数角度看，静止目标几乎不挑导航比，N=3 时就能得到平滑弹道。这个场景里值得观察的是初始几何偏差对弹道弯曲程度的影响：初始横向偏差越大，前段轨迹弯曲越明显，限幅越容易触发。有的初学者把初始位置偏差设在几公里外，期望导弹从很远就开始修正，结果初始时刻限幅一直顶满，弹道呈现一条折线。这正是对"比例导引是渐进修正"理解不到位。静止目标场景中把初始横向距离控制在几百米量级，看到的才是比例导引的典型平滑弧形。

4.2 正弦机动目标：sinmotorizeda3dnew.m 的模型与参数

sinmotorizeda3dnew.m 建立的是横向加速度随时间正弦变化的目标，模拟飞行器做周期性蛇形机动时的横向运动规律。目标在某个方向保持匀速或匀加速，在另一个垂直通道上加幅值为 a_t、频率为 f 的正弦加速度。这两个参数直接决定目标机动强度：幅值常用 1g 到 3g，频率常用 0.1 到 0.5 赫兹。幅度太大时导弹需要在很晚时刻仍保持高机动输出，频率太高时视线角速率更新率跟不上，制导滞后会被明显放大。

正弦机动目标对比例导引的影响来自持续加速度。对静止目标 N=3 够用，但对持续横向摆动的目标，N 取 3 时末端会有可观察的跟踪误差，N 取 4 到 5 时误差明显压低，代价是末端指令噪声变大。我在用这类资源时会固定 a_t=2g、f=0.2 赫兹，把 N 从 3 到 5 各跑一遍，画出脱靶量随导航比的变化曲线，这比反复盯轨迹直观得多。

需要多说一句：正弦机动目标的状态更新应该写在微分方程右侧，而不是在主循环里用解析位置开环更新。常见误用是直接在循环里写"目标位置等于某个正弦函数"，让目标凭空移动，这会破坏 RK4 对系统状态的统一积分。正确做法是目标模型只返回目标加速度，位置和速度由 rungekutta3d.m 与导弹一起推进。

% sinmotorizeda3dnew.m 的目标加速度函数（示意） function at = sinmotorizeda3dnew(t, params) % params: 目标模型参数，一般为结构体 % params.a_t: 横向机动幅值，单位 m/s^2 % params.f_t: 机动频率，单位 Hz % 示意：目标沿 x 向匀速，y 向作正弦摆动，z 向无机动 at = [0; params.a_t * sin(2*pi*params.f_t*t); 0]; end

这段代码表明目标模型只是把加速度算出来，剩下交给积分器。想让机动方向从 y 轴换成 z 轴，把零分量和正弦项位置调换即可；想让目标在正弦摆动的同时向前匀速冲，就在目标初始速度里把纵向分量设成巡航速度。改变机动幅度时，a_t 直接决定比例导引的压力水平，建议从小幅值往上加，观察脱靶量的变化趋势。

4.3 方波机动目标：squaremotorizeda3dnew.m 的切换逻辑

squaremotorizeda3dnew.m 建模的是方波机动，目标横向加速度在两个极限值之间周期性跳变，典型场景对应飞行器短时间内切换大幅度横摆或做规避性阶跃加速。方波机动比正弦机动的仿真难度高一个档次，因为切换点上目标加速度不连续，视线角速率出现跳变，比例导引输出随之出现一个尖峰。尖峰超过限幅时，导弹在那一瞬间进入饱和；如果接近速度很高，留给制导律修正的时间窗口非常小。

从评价角度讲，方波机动的"杀伤力"比正弦模型更强。正弦机动属于能量受限的连续机动，方波机动在切换瞬间含有高频成分，对制导律响应速度、导航比和积分步长都更敏感。因此仿真设计里通常用正弦模型做常规性能测试，用方波模型做边界能力测试。同一个 a_t 幅值下，方波目标造成的脱靶量往往显著大于正弦目标，这不是代码 bug，而是目标"变向更快"的物理体现。

% squaremotorizeda3dnew.m 的目标加速度函数（示意） function at = squaremotorizeda3dnew(t, params) % 方波机动：横向加速度在两个值之间周期切换 % 切换半周期 T_switch，单位 s；a_sw 为幅值，单位 m/s^2 T_switch = params.T_switch; a_sw = params.a_sw; if mod(t, 2*T_switch) < T_switch a_lat = a_sw; else a_lat = -a_sw; end at = [0; a_lat; 0]; end

这段逻辑用 mod 函数把时间映射到切换周期，小于半周期为正、大于半周期为负。做对比实验时有个实用技巧：把 T_switch 和正弦模型的周期设为同一数值、幅值也相同，就能在相近机动周期下比较连续机动和阶跃机动的差异。建议的正弦-方波实验矩阵是：正弦模型下跑 N=3、4、5 三组，方波模型下跑 N=4、5、6 三组，把脱靶量记录成一张表，基本就能看出比例导引对两类机动的适应边界。方波条件下缩小积分步长意义更大，dt 从 0.01 降到 0.005 时，方波场景的脱靶量通常会继续下降一个台阶。

5. 避坑与常见问题：参数飞了、轨迹发散、脱靶量下不去

这一章集中写复现这类三维比例导引资源时最常遇到的五类问题，每一条都能直接对上现象。这些坑不是推导出来的，是实际跑代码踩出来的，很多问题在静止目标下不出现，一换到机动目标就暴露。

5.1 步长太大：固定步长 RK4 在方波目标下发散

现象：静止目标跑得好好的，换到方波机动目标后，弹道最后几米处出现剧烈抖动，弹目距离曲线末端上下跳动，脱靶量统计出来反而比正弦目标大几十倍。 原因：目标加速度在切换点不连续，RK4 对不连续激励的积分误差和步长成正比。dt=0.05 时，切换点引入的大误差被后续几步放大，视觉上就是弹道末端在"抽搐"。这不是比例导引算法的问题，是数值积分精度不足。 解决：把 dt 从 0.05 降到 0.01，如果现象依旧，继续降到 0.005，观察脱靶量是否随步长减半而明显改善。固定步长 RK4 在步长逐级减半时，脱靶量应当逐步趋近同一稳定值；若仍反复跳动，则目标函数在切换点需要做平滑处理，或者干脆改用变步长积分器。这个坑的实操价值在于：先降步长再改算法，不要一看到发散就去调导航比 N，那是拿错误的积分误差去补偿另一个错误。

5.2 初始视线角速率爆表：开头就顶满限幅

现象：仿真第 0 秒开始，导弹指令加速度就顶着 a_max 不放，轨迹不是平滑弧形而是折线，有时甚至绕一个大圈才回来。 原因：初始弹目几何设置不合理。目标初始位置在弹目连线横向方向分量过大，或者目标初始速度方向与弹目连线接近垂直，都会导致初始视线角速率非常大。视线角速率经过 N 倍放大后远超可用过载，从第一拍起就饱和。注意比例导引是比例控制，不是 bang-bang 控制，指令被顶到限幅时实际已经变成开环制导，命中性能无从谈起。 解决：先检查初始横向距离和横向相对速度，把初始条件调成视线角速率在合理范围的场景。一般经验是初始视线角速率不超过每秒几度到十几度，换算成弧度后乘以接近速度和导航比，应明显低于 a_max。如果要测大初始误差场景，应该分级加大初始偏差，而不是一步到位。

5.3 限幅单位混乱：g 与 m/s² 混用

现象：同一份代码在一台机器上运行脱靶量 3 米，换台机器只改一处参数就变成 300 米，代码逻辑看起来完全没动。 原因：a_max 在代码内部使用米每二次方秒，注释里却写"过载 5g"。赋值时写成 a_max = 5，实际只有约 0.5g，导弹根本没有能力完成末端机动；另一次写成 9.8 乘以 5 等于 49 米每二次方秒，又恢复正常。单位不一致在弹道仿真里不会像编译错误那样报错，只会让你看到异常弹道，然后开始怀疑制导律写错了。 解决：在所有代码文件里统一使用米每二次方秒，只在绘图或输出时除以 9.8 得到过载值。在参数区写注释时明确标注当前数值对应几个 g，例如 a_max = 49; 百分号后写等于 5g。每次改这个参数前先换算，这是成本最低的防呆习惯。

5.4 状态向量顺序不一致：RK4 把导数叠加错了位置

现象：单独跑静止目标正常，把目标模型换成机动目标后，导弹和目标轨迹出现怪异的分段现象，甚至目标做正弦机动时，导弹轨迹也跟着同步摆动。 原因：12 维状态向量里导弹状态与目标状态的排列顺序在两个文件里不一致。目标模型按自己的顺序返回加速度，主脚本组装状态向量时却按另一套顺序排列，导入加速度时索引错位。RK4 不理解物理含义，只会忠实地把导数按顺序加到状态上，顺序错一个块，后面每一步都叠加错误。 解决：在 bili3dnew.m 顶部明确写一段状态向量定义注释：1 到 3 为导弹位置，4 到 6 为导弹速度，7 到 9 为目标位置，10 到 12 为目标速度。目标模型函数的返回接口固定为只返回目标加速度，由主脚本负责映射到对应位置。这样换目标模型时只动目标文件，不碰主流程的状态索引。

5.5 命中判定过早：距离阈值太大掩盖性能问题

现象：仿真结束条件写成"弹目距离小于 10 米即命中"，脱靶量报告永远不超过 10 米。再调什么参数结果都一样，看起来性能优越，实际算法水平已经无法区分了。 原因：把命中阈值当成了最大允许脱靶量。阈值设得越大，评价越宽松，越容易掩盖制导律本身的问题。比如真实最小脱靶量是 2 米，阈值 10 米时仿真在第 10 米处就停止，后续更近的距离根本不会被记录。 解决：命中阈值按实际引信和战斗部要求来，弹道级验证通常设 1 米甚至 0.5 米。记录最终距离后，再单独计算全时间序列的最小弹目距离，两个指标同时看。主循环里结束条件是阈值判断，脱靶量统计用整体最小值函数，而不是取跳出循环那一刻的距离。这也是第 3 章末那段循环代码的核心用意。

提示：所有避坑条目都适用于包内代码的复现调试。如果某条现象对不上，先检查自己的 MATLAB 版本是否把浮点运算的细节改了，再回来对照状态向量定义。

6. 命中判定与导航比配法：脱靶量怎么统计、N 怎么取

6.1 脱靶量统计口径：梯度下降的曲线才是好曲线

弹道仿真里最容易误用的指标是"仿真结束时刻的弹目距离"。当导弹穿过目标或者在限幅饱和下往复振荡时，结束时刻可能是距离反弹后的值，远大于真实的最小脱靶量。正确的统计方式是对整个时间序列求最小值：

% 从输出数组中统计脱靶量（示意） dist_ser = sqrt(sum((r_m - r_t).^2, 2)); % 每个时间步的弹目距离 miss = min(dist_ser); % 全弧段最小值即脱靶量

统计之外更要看曲线的形状。健康的弹目距离曲线是单调下降的，即使末端有起伏，也应该在命中的瞬间触及全局最小值。如果距离曲线像波浪一样先降后升再降，说明制导过程在某个时段失去了对视线角速率的控制力。我从这类资源调参的习惯是：每次改完参数，先看距离曲线形态，再看脱靶量数值，形态不对时数值没有参考意义。

6.2 导航比 N 的三档配法

导航比 N 是比例导引唯一真正的算法参数。取低了跟踪刚性不足，取高了放大视线角速率噪声，末端指令毛刺增多。按这个资源支持的三类目标，推荐配法如下：

对机动目标来说，N 从 4 加到 5 的效果在方波场景下比正弦场景更明显，因为方波切换瞬间的视线角速率跳变需要更强的指令响应。但 N 超过 5 之后收益递减，反而会让末段指令抖动变大。一般我不建议在这个弹道级仿真里把 N 拉到 6 以上，除非同时降低积分步长并给视线角速率做滤波。

6.3 一个能带走的调试习惯：先静止、再正弦、后方波

这类三维比例导引仿真包里最容易改动的地方是比例导引函数和目标模型之间的接口。我自己的调试习惯是固定一套流程：先静止目标下把 N 和 a_max 调到脱靶量小于 1 米，再切正弦目标确认跟踪能力，最后换方波目标看边界表现。每次只动一个变量，记录脱靶量和距离曲线，状态向量定义每次先在注释里核对一遍。这套流程是一个做飞行器仿真的前辈反复强调的，从那以后我每次拿到新的制导仿真资源，都强制先走一遍静止到机动的完整链条，再谈改参数或加新功能。这份资源虽然只有弹道级的五段代码，但把基线跑稳之后，往六自由度模型或末端制导律扩展时，接口和统计口径都不用推倒重来。希望帮到你。

本文还有配套的精品资源，点击获取