尧图网站开发 
简介这是一套面向飞行器制导控制研究者的攻击角约束制导律Matlab仿真程序包覆盖比例导引与变结构导引两类典型算法适用于本科毕业设计、研究生课题或科研人员的算法验证与性能对比。资源共474个文件压缩包大小6.65MB主要包含mat数据文件、c/h源码文件、obj编译中间文件、mexw64编译接口、slx模型文件及m脚本等既有完整Simulink模型也有便于直接运行的编译版本支持快速仿真测试与二次开发。目前已有222人学习下载。程序包提供多种攻击角约束策略实现如变结构导引律、比例导引以及带落角约束的中末制导模型并附带参数配置脚本与批处理文件读者可参照模型结构复现仿真结果分析导引参数对攻击角、过载等指标的影响对深入理解制导律设计具有较强工程参考价值。1. 攻击角约束制导律仿真能做、也能复现的最小闭环攻击角约束制导律说到底就是把“命中目标”这一个要求拆成两个指标脱靶量趋近于零终端攻击角同时落在期望范围内。常规比例导引只处理前者所以它能把导弹导到目标附近却无法保证以预定的角度命中而带攻击角约束的比例导引和变结构制导律是在原有制导指令之上再叠加一个与终端角度偏差相关的修正项。导弹制导仿真里最有价值的地方也在这里同一套相对运动方程换一段制导指令就能在 MATLAB 里清楚看到导弹飞行轨迹如何弯曲、攻击角误差如何被逐步修正、脱靶量和角度约束又是怎么互相换来的。本文把这套方案拆成四个部分攻击角约束的对象和几何定义、比例制导律与滑模变结构制导律的指令构造逻辑、可直接在 MATLAB 里跑通的最小仿真程序以及参数标定、发散排查和批量验证。适合做末段拦截方案对比、制导律课程设计或者在答辩前临时补仿真数据的工程师和研究生。下文所有代码只用到基本 MATLAB 语法不依赖 Simulink也不依赖任何工具箱R2020 以上版本都能直接跑。2. 攻击角约束的两条设计路线比例导引与滑模变结构制导律2.1 攻击角到底约束的是哪个角制导律里的“攻击角”在不同论文里有不同叫法终端角度约束、落角约束、碰撞角约束。它们说的其实是同一件事——命中瞬间导弹速度矢量与参考方向之间的夹角。二维平面里用航迹角 θm 描述导弹速度方向θt 描述目标速度方向攻击角约束就是让 θm(tf) 在命中时刻 tf 逼近期望值 θd。对于固定目标终端几何关系会进一步简化。导弹指向目标的视线角 q 与航迹角 θm 之间满足命中时刻若有 r→0 且视线角速率 qdot→0则 θm(tf)q(tf)。所以固定目标场景里约束攻击角可以等效为约束终端视线角。这也是后面比例导引指令里可以直接用 (q - θd) 做修正项的原因。目标做匀速直线运动时需要把 θd 换成“满足碰撞三角关系的期望视线角”公式会多一个由目标速度方向组成的偏差项机动目标则要基于目标加速度估计值在线补偿。常见做法是先按固定目标或匀速目标把主回路跑通再把目标机动当作扰动由变结构项去抑制。2.2 带攻击角项的比例导引怎么改写纯比例导引的指令加速度形式是a_cmd N1 · Vc · qdot其中 Vc 是接近速度qdot 是视线角速率N1 是导航比。这个形式只能把 qdot 压到零终端攻击角完全取决于初始几何。要加攻击角约束最常用的做法是增加一个“角度偏置项”a_cmd N1 · Vc · qdot N2 · Vc · (q - θd) / tgo这里 tgo 是剩余飞行时间估计最粗糙但实用的估计是 tgo r / Vc。第二项的作用是这样的当视线角 q 还没有收敛到期望攻击角 θd 时它会生成一个与角度误差成正比的过载指令并且误差项被 1/tgo 放大。飞行末段 tgo 不断变小这个偏置项的影响力越来越大从而把终端航迹方向拧到期望值附近。习惯上把这种结构叫“带攻击角约束的扩展比例导引”或“偏置比例导引”实现成本最低工程上也最容易定性解释第一项管脱靶量第二项管攻击角。参数 N2 的取值很关键。N2 太小角度修正能力弱终段只剩比例项在起作用N2 太大末段视线角会被强行拉离目标方向脱靶量反而增加。一般场景下 N1 取 3 到 5N2 取 1 到 4 比较稳妥具体值要靠仿真扫描得到一个兼顾两个指标的区间。2.3 变结构滑模制导律的优势与边界变结构制导律也就是滑模变结构制导律设计思路不是直接拼一个角度误差项而是构造一个能同时反映“视线角速率收敛”和“终端角度收敛”的滑动面。常见取法如下s qdot λ · (q - θd) / tgo当滑动面 s 收敛到 0 时qdot 与 (q - θd)/tgo 保持比例关系距离还很远时q 可以缓慢向 θd 靠拢接近目标时tgo 减小迫使 qdot 也变小最终实现 q→θd 且 qdot→0。制导指令用趋近律构造a_cmd N1 · Vc · qdot ε₁ · tanh(s / φ)其中 tanh(s/φ) 是符号函数 sign(s) 的连续化形式φ 称为边界层厚度。变结构体现在这一项当外界扰动或目标机动使 s 偏离零点时tanh 项会提供一个与 s 符号相反的修正过载把系统状态拉回滑动面。变结构方案对目标机动和模型不确定性的鲁棒性明显强于单纯的比例偏置方案。代价是需要多调两个参数λ 决定滑动面收敛速度ε₁ 决定切换强度的上限φ 决定抖振抑制程度。φ 取得过大就退化成一个线性反馈滑模特性不明显φ 取得过小切换作用接近理想的 sign(s)数字仿真里指令加速度会高频翻转表现为曲线上的锯齿抖振。在 MATLAB 里跑变结构制导律最常见的问题就是这个后面第 4 章专门展开。2.4 选用比例还是变结构的判据两条路线并不是替代关系而是应用场景不同。固定目标、目标机动幅度小、算力预算紧张的场景带攻击角约束的比例导引完全够用参数少仿真和半实物都容易调试。目标存在机动、或者传感器噪声明显希望制导指令不被目标加速度扰动带偏优先选变结构方案。还有一个工程判据是看可用过载余量。比例偏置项在末段会因 1/tgo 放大产生很大的指令加速度如果导弹可用过载有限偏置项先饱和攻击角误差只能靠剩余航程慢慢修正变结构方案虽然也有切换项但可以用边界层和趋近律限幅限制行为更稳定。实际项目里两种都会写进同一套 MATLAB 仿真架构用配置开关切换和标题里“比例变结构均有”的说法一致。3. 在 MATLAB 中搭建攻击角约束制导仿真状态方程与最小主程序3.1 仿真回路与坐标定义二维平面仿真里状态向量取 [xm; ym; θm; xt; yt; θt]前三个是导弹位置和航迹角后三个是目标位置和航迹角。制导指令是垂直于导弹速度方向的法向加速度 a_cmd单位 m/s²它改变航迹角速率θm_dot a_cmd / Vm导弹位置随时间的变化满足运动学关系xm_dot Vm · cos(θm) ym_dot Vm · sin(θm)目标侧同理。固定目标时 Vt0目标位置不变。相对距离 r、视线角 q、视线角速率 qdot 从状态量里实时算出来作为制导律的输入。整个仿真回路是“由状态求 r、q、qdot → 代入制导指令 → 更新导弹航迹角 → 积分一步 → 判断脱靶”这是所有平面制导仿真共同的骨架。仿真步长 dt 取 0.001 秒积分用四阶 Runge-Kutta。步长再大一些欧拉法也能跑但变结构末段的抖振会被数值积分方法放大所以这里直接上 RK4代码量只多几行稳定性却好得多。3.2 制导律指令函数与状态方程代码把制导指令的计算单独写成函数比例和变结构两种方案通过参数切换。下面是完整的可运行版本场景是导弹从原点出发拦截位于 (6000, 0) 的固定目标期望攻击角 θd -20°。function a_cmd guidance_command(q, qdot, r, Vc, theta_d, tgo, mode, param) % 攻击角约束制导指令 % mode PN_ANG 带攻击角偏置的比例导引 % mode VSS 滑模变结构制导律 switch mode case PN_ANG a_cmd param.N1 * Vc * qdot ... param.N2 * Vc * (q - theta_d) / tgo; case VSS s qdot param.lambda_s * (q - theta_d) / tgo; a_cmd param.N1 * Vc * qdot ... param.eps_s * tanh(s / param.phi); otherwise error(unknown guidance mode); end end函数的输入含义如下q 是当前视线角qdot 是视线角速率Vc 是接近速度tgo 是剩余飞行时间theta_d 是期望攻击角。param 这个结构体里放着各制导律自己的参数。代码里有一个细节需要注意VSS 分支里 tanh(s/phi) 完全取代了 sign(s)这就是边界层方法。s 的物理单位是 1/stanh 的输入无量纲所以 phi 要和 s 保持同一个量级一般取 0.01 到 0.1。下面是状态方程函数。它接收当前状态 X、指令加速度 a_cmd返回状态导数function dX plant(X, a_cmd, Vm, Vt) thm X(3); tht X(6); dX [Vm * cos(thm); Vm * sin(thm); a_cmd / Vm; % 航迹角速率由法向过载给出 Vt * cos(tht); Vt * sin(tht); 0]; end固定目标 Vt0所以第 4、5、6 个分量都是 0。追逐状态目标时把 Vt 和 θt 设成目标的速度和航向即可同一个 plant 函数不需要改。3.3 主循环RK4 积分与脱靶判定主程序把制导指令计算和内层积分串起来。每个仿真步长内先根据当前状态算出 a_cmd然后以这个 a_cmd 为输入做四步 RK4。这个处理方式在工程上等价于“制导指令零阶保持”即认为一个计算周期内指令不变符合数字制导计算机的实际工作方式。% 攻击角约束制导律仿真主程序 clear; clc; close all; Vm 300; Vt 0; X0 [0; 0; deg2rad(25); 6000; 0; deg2rad(0)]; theta_d deg2rad(-20); param.N1 4; % 比例项导航比 param.N2 2.5; % 偏置项增益 param.lambda_s 4; % 滑模面系数 param.eps_s 30; % 切换增益 param.phi 0.05; % 边界层厚度 mode VSS; % 或 PN_ANG dt 0.001; t_max 30; X X0; t 0; r_min 1e6; hist.t[]; hist.X[]; hist.a[]; hist.q[]; while t t_max xmX(1); ymX(2); thmX(3); xtX(4); ytX(5); thtX(6); dx xt - xm; dy yt - ym; r sqrt(dx^2 dy^2); if r 0.5, break; end q atan2(dy, dx); qdot (dx*(Vt*sin(tht)-Vm*sin(thm)) ... - dy*(Vt*cos(tht)-Vm*cos(thm))) / r^2; rdot Vt*cos(tht-q) - Vm*cos(thm-q); Vc -rdot; tgo r / max(Vc, 1); a_cmd guidance_command(q, qdot, r, Vc, theta_d, tgo, mode, param); % RK4 更新 k1 plant(X, a_cmd, Vm, Vt); k2 plant(X0.5*dt*k1, a_cmd, Vm, Vt); k3 plant(X0.5*dt*k2, a_cmd, Vm, Vt); k4 plant(Xdt*k3, a_cmd, Vm, Vt); X X dt/6*(k1 2*k2 2*k3 k4); t t dt; hist.t(end1)t; hist.X(:,end1)X; hist.a(end1)a_cmd; hist.q(end1)q; end theta_m_final X(3); err_angle rad2deg(theta_m_final - theta_d); fprintf(脱靶量 r_min%.4f m, 终端攻角误差%.4f deg\n, r, err_angle);qdot 的计算式用的是相对坐标求导的直接展开不需要额外定义视线坐标系代码简单且不容易出现符号错误。tgo 用 max(Vc, 1) 是防除零Vc 接近 0 时仿真早已经在脱靶判定之前跳出了。r 0.5 作为命中判定条件意味着脱靶量小于半米就认为命中。3.4 结果可视化与参数对应关系仿真结束后至少要看两条曲线导弹飞行轨迹以及终端攻击角误差随时间的变化。轨迹用 plot 画 Y 随 X 变化攻击角误差用 rad2deg(theta_d) 归一化后看末段是否收敛到 0。还有一个很重要的检查对象是 a_cmd 曲线变结构制导律如果出现沿零线附近高频抖振大多能在 a_cmd 历史曲线里一眼看出。固定目标场景里命中时刻指令加速度 a_cmd 大概率不为零。原因很简单制导律还在抑制 qdot 的残差直到最后一刻导弹都在微调航迹方向。这是正常现象不代表程序有 bug。4. 攻击角约束制导律仿真的参数设定与坑排查4.1 三个必须调的参数与判定依据带攻击角约束的比例导引里最需要调的是 N1 和 N2。N1 影响视线角速率收敛速度N1 太大指令饱和N1 太小飞行轨迹弯曲度不够N2 直接决定攻击角约束的强度。判断依据很简单看终端航迹角相对期望值的偏差以及脱靶量是否被第二项破坏。先固定 N14把 N2 从 1 扫到 4通常能在某个值附近同时满足脱靶量小于 0.5 米、攻击角偏差小于 0.5 度。变结构制导律要调三个量λ、ε₁、φ。λ 决定滑动面收敛快慢可以参考比例导引中 N2/N1 的比例关系起调。ε₁ 必须大于目标机动和模型不确定性造成的等效扰动上界否则滑动面无法建立攻击角误差无法消除。φ 控制抖振和稳态精度之间的平衡φ0.01 时指令加速度会出现明显的锯齿φ0.1 时曲线平滑但终端角度误差可能偏大。这三个参数没有固定解析公式工程上都是先定 λ再加 ε₁最后用 φ 磨平曲线。4.2 脱靶量与攻击角误差的取舍两个指标在同一个仿真里经常是矛盾的。偏置项过强导弹会在末端绕一个很大的弧线去“够”期望角度轨迹弯过头脱靶量放大偏置项太弱攻角约束满足不了。碰到这种情况先检查是不是指令饱和造成的而不是盲目调整 N2 或 ε₁。用 max(a_cmd, 可用过载) 做饱和限幅是更贴近实际的改法amax 100; % 最大法向加速度 m/s^2 a_cmd max(min(a_cmd, amax), -amax);限幅后终端攻击角误差增大是预期结果说明导弹可用过载不足以支撑期望弹道的弯曲。真实系统中这里对应的是气动舵面偏转限制不能因为仿真没限幅就忽略。限幅之后还可以把制导律切换为“先比例、末端角度修正”的分段策略保证飞行中段有足够的能量。4.3 仿真不收敛或曲线发散的常见原因第一个原因是单位不统一。atan2 返回弧度如果 theta_d 用角度直接代入两个量在数值上差了约 57.3 倍仿真结果会乱到无法解释。脚本里统一用 deg2rad 作为唯一入口绘图输出时才转回角度。第二个原因是 tgo 估计失真。接近速度 Vc 在固定目标场景末段未必单调趋近于 Vm如果目标有速度Vc 可能过零直接用 r/Vc 会算出无穷大的偏置项。用 max(Vc, 1) 或者加一个下限保护是好习惯同时要在循环开头检查 r 是否已经小于脱靶阈值。第三个原因是积分方法跟不上抖振。变结构制导律里 φ 取得过小时tanh(s/φ) 接近 sign(s)指令加速度在一个仿真步长内急剧翻转。此时 RK4 的表现只比欧拉好一点彻底解决要靠“保留滑模切换项但外加低通滤波”的工程处理。常见做法是给 a_cmd 加一阶惯性环节a_cmd_filt a_cmd_filt dt/tau * (a_cmd - a_cmd_filt);τ 取 0.02 到 0.05 秒相当于模拟制导系统本身的动态响应曲线会平滑很多而攻击角约束精度损失通常小于 0.1 度。这一步在论文仿真中可以省略但在考虑工程可实现性时必须保留。最后排查一种容易被忽略的现象仿真在 r 0.5 之前就提前跳出原因是目标初始位置或速度配置导致弹道越过目标但没有满足接近条件。这时不用急着调参数先画出轨迹看看是不是初始航向角与目标方向相差过大末端弹道绕过了目标后再折返。制导系统的捕获区有限初始航向不能偏离视线方向太远这是物理约束不是程序问题。5. 批量仿真验证把攻击角约束结果与参数矩阵同时沉淀下来验证一条制导律是否真的“带攻击角约束”单跑一条弹道说服力不够至少要扫一遍期望攻击角的覆盖范围。把上一章的主循环包成函数输出脱靶量、终端航迹角和指令加速度序列然后在主脚本里循环不同 θd、不同初始航向角。这个动作既能把整个数据变成表格也顺便暴露参数矩阵里那些单条弹道看不见的边界情况。theta_d_list deg2rad(-40:10:20); % 期望攻击角扫描 N numel(theta_d_list); results table(Size,[N,4], ... VariableTypes, {double,double,double,double}, ... VariableNames, {theta_d_deg,miss_m,angle_err_deg,max_a} ); for i 1:N [miss, err_deg, a_hist] run_guidance_scenario(theta_d_list(i)); results.theta_d_deg(i) rad2deg(theta_d_list(i)); results.miss_m(i) miss; results.angle_err_deg(i) err_deg; results.max_a(i) max(abs(a_hist)); end disp(results);run_guidance_scenario 把上一章从“参数初始化”到“脱靶判定”的全部代码收进去返回脱靶量、终端攻击角误差、指令加速度曲线。扫描结果里通常能看到明显边界某个期望攻击角和初始几何之间一旦超过一定角度差攻击角误差开始快速增大而脱靶量依然很小。这个边界值就是该制导律的可用捕获区写报告或答辩时比单条弹道光标图有用得多。更进一步扫描 N2 或 ε₁ 与 θd 的二维网格用 surf 或 imagesc 画“攻击角误差随参数和期望角变化”的等高线图可以直观标出参数可用区间。这个做法对不同版本的 MATLAB 都通用不涉及任何工具箱R2023b、R2026a、R2026b 安装包装好之后直接跑脚本即可。有人问现在带 AI 能力的代码生成工具能不能像执行 Python 一样直接操作 MATLAB 任务批量跑这种参数扫描自然可以但制导律的物理边界、指令限幅和积分步长设置仍然需要有人把关这正是写清楚主循环结构的意义所在。跑完批量扫描后建议把 results 表存成 CSV用 MATLAB 自带的 writetable 导出后续画图、写文档、做对比分析都可以直接复用同一份数据。本文还有配套的精品资源点击获取 
简介这是一套面向飞行器制导控制研究者的攻击角约束制导律Matlab仿真程序包覆盖比例导引与变结构导引两类典型算法适用于本科毕业设计、研究生课题或科研人员的算法验证与性能对比。资源共474个文件压缩包大小6.65MB主要包含mat数据文件、c/h源码文件、obj编译中间文件、mexw64编译接口、slx模型文件及m脚本等既有完整Simulink模型也有便于直接运行的编译版本支持快速仿真测试与二次开发。目前已有222人学习下载。程序包提供多种攻击角约束策略实现如变结构导引律、比例导引以及带落角约束的中末制导模型并附带参数配置脚本与批处理文件读者可参照模型结构复现仿真结果分析导引参数对攻击角、过载等指标的影响对深入理解制导律设计具有较强工程参考价值。1. 攻击角约束制导律仿真能做、也能复现的最小闭环攻击角约束制导律说到底就是把“命中目标”这一个要求拆成两个指标脱靶量趋近于零终端攻击角同时落在期望范围内。常规比例导引只处理前者所以它能把导弹导到目标附近却无法保证以预定的角度命中而带攻击角约束的比例导引和变结构制导律是在原有制导指令之上再叠加一个与终端角度偏差相关的修正项。导弹制导仿真里最有价值的地方也在这里同一套相对运动方程换一段制导指令就能在 MATLAB 里清楚看到导弹飞行轨迹如何弯曲、攻击角误差如何被逐步修正、脱靶量和角度约束又是怎么互相换来的。本文把这套方案拆成四个部分攻击角约束的对象和几何定义、比例制导律与滑模变结构制导律的指令构造逻辑、可直接在 MATLAB 里跑通的最小仿真程序以及参数标定、发散排查和批量验证。适合做末段拦截方案对比、制导律课程设计或者在答辩前临时补仿真数据的工程师和研究生。下文所有代码只用到基本 MATLAB 语法不依赖 Simulink也不依赖任何工具箱R2020 以上版本都能直接跑。2. 攻击角约束的两条设计路线比例导引与滑模变结构制导律2.1 攻击角到底约束的是哪个角制导律里的“攻击角”在不同论文里有不同叫法终端角度约束、落角约束、碰撞角约束。它们说的其实是同一件事——命中瞬间导弹速度矢量与参考方向之间的夹角。二维平面里用航迹角 θm 描述导弹速度方向θt 描述目标速度方向攻击角约束就是让 θm(tf) 在命中时刻 tf 逼近期望值 θd。对于固定目标终端几何关系会进一步简化。导弹指向目标的视线角 q 与航迹角 θm 之间满足命中时刻若有 r→0 且视线角速率 qdot→0则 θm(tf)q(tf)。所以固定目标场景里约束攻击角可以等效为约束终端视线角。这也是后面比例导引指令里可以直接用 (q - θd) 做修正项的原因。目标做匀速直线运动时需要把 θd 换成“满足碰撞三角关系的期望视线角”公式会多一个由目标速度方向组成的偏差项机动目标则要基于目标加速度估计值在线补偿。常见做法是先按固定目标或匀速目标把主回路跑通再把目标机动当作扰动由变结构项去抑制。2.2 带攻击角项的比例导引怎么改写纯比例导引的指令加速度形式是a_cmd N1 · Vc · qdot其中 Vc 是接近速度qdot 是视线角速率N1 是导航比。这个形式只能把 qdot 压到零终端攻击角完全取决于初始几何。要加攻击角约束最常用的做法是增加一个“角度偏置项”a_cmd N1 · Vc · qdot N2 · Vc · (q - θd) / tgo这里 tgo 是剩余飞行时间估计最粗糙但实用的估计是 tgo r / Vc。第二项的作用是这样的当视线角 q 还没有收敛到期望攻击角 θd 时它会生成一个与角度误差成正比的过载指令并且误差项被 1/tgo 放大。飞行末段 tgo 不断变小这个偏置项的影响力越来越大从而把终端航迹方向拧到期望值附近。习惯上把这种结构叫“带攻击角约束的扩展比例导引”或“偏置比例导引”实现成本最低工程上也最容易定性解释第一项管脱靶量第二项管攻击角。参数 N2 的取值很关键。N2 太小角度修正能力弱终段只剩比例项在起作用N2 太大末段视线角会被强行拉离目标方向脱靶量反而增加。一般场景下 N1 取 3 到 5N2 取 1 到 4 比较稳妥具体值要靠仿真扫描得到一个兼顾两个指标的区间。2.3 变结构滑模制导律的优势与边界变结构制导律也就是滑模变结构制导律设计思路不是直接拼一个角度误差项而是构造一个能同时反映“视线角速率收敛”和“终端角度收敛”的滑动面。常见取法如下s qdot λ · (q - θd) / tgo当滑动面 s 收敛到 0 时qdot 与 (q - θd)/tgo 保持比例关系距离还很远时q 可以缓慢向 θd 靠拢接近目标时tgo 减小迫使 qdot 也变小最终实现 q→θd 且 qdot→0。制导指令用趋近律构造a_cmd N1 · Vc · qdot ε₁ · tanh(s / φ)其中 tanh(s/φ) 是符号函数 sign(s) 的连续化形式φ 称为边界层厚度。变结构体现在这一项当外界扰动或目标机动使 s 偏离零点时tanh 项会提供一个与 s 符号相反的修正过载把系统状态拉回滑动面。变结构方案对目标机动和模型不确定性的鲁棒性明显强于单纯的比例偏置方案。代价是需要多调两个参数λ 决定滑动面收敛速度ε₁ 决定切换强度的上限φ 决定抖振抑制程度。φ 取得过大就退化成一个线性反馈滑模特性不明显φ 取得过小切换作用接近理想的 sign(s)数字仿真里指令加速度会高频翻转表现为曲线上的锯齿抖振。在 MATLAB 里跑变结构制导律最常见的问题就是这个后面第 4 章专门展开。2.4 选用比例还是变结构的判据两条路线并不是替代关系而是应用场景不同。固定目标、目标机动幅度小、算力预算紧张的场景带攻击角约束的比例导引完全够用参数少仿真和半实物都容易调试。目标存在机动、或者传感器噪声明显希望制导指令不被目标加速度扰动带偏优先选变结构方案。还有一个工程判据是看可用过载余量。比例偏置项在末段会因 1/tgo 放大产生很大的指令加速度如果导弹可用过载有限偏置项先饱和攻击角误差只能靠剩余航程慢慢修正变结构方案虽然也有切换项但可以用边界层和趋近律限幅限制行为更稳定。实际项目里两种都会写进同一套 MATLAB 仿真架构用配置开关切换和标题里“比例变结构均有”的说法一致。3. 在 MATLAB 中搭建攻击角约束制导仿真状态方程与最小主程序3.1 仿真回路与坐标定义二维平面仿真里状态向量取 [xm; ym; θm; xt; yt; θt]前三个是导弹位置和航迹角后三个是目标位置和航迹角。制导指令是垂直于导弹速度方向的法向加速度 a_cmd单位 m/s²它改变航迹角速率θm_dot a_cmd / Vm导弹位置随时间的变化满足运动学关系xm_dot Vm · cos(θm) ym_dot Vm · sin(θm)目标侧同理。固定目标时 Vt0目标位置不变。相对距离 r、视线角 q、视线角速率 qdot 从状态量里实时算出来作为制导律的输入。整个仿真回路是“由状态求 r、q、qdot → 代入制导指令 → 更新导弹航迹角 → 积分一步 → 判断脱靶”这是所有平面制导仿真共同的骨架。仿真步长 dt 取 0.001 秒积分用四阶 Runge-Kutta。步长再大一些欧拉法也能跑但变结构末段的抖振会被数值积分方法放大所以这里直接上 RK4代码量只多几行稳定性却好得多。3.2 制导律指令函数与状态方程代码把制导指令的计算单独写成函数比例和变结构两种方案通过参数切换。下面是完整的可运行版本场景是导弹从原点出发拦截位于 (6000, 0) 的固定目标期望攻击角 θd -20°。function a_cmd guidance_command(q, qdot, r, Vc, theta_d, tgo, mode, param) % 攻击角约束制导指令 % mode PN_ANG 带攻击角偏置的比例导引 % mode VSS 滑模变结构制导律 switch mode case PN_ANG a_cmd param.N1 * Vc * qdot ... param.N2 * Vc * (q - theta_d) / tgo; case VSS s qdot param.lambda_s * (q - theta_d) / tgo; a_cmd param.N1 * Vc * qdot ... param.eps_s * tanh(s / param.phi); otherwise error(unknown guidance mode); end end函数的输入含义如下q 是当前视线角qdot 是视线角速率Vc 是接近速度tgo 是剩余飞行时间theta_d 是期望攻击角。param 这个结构体里放着各制导律自己的参数。代码里有一个细节需要注意VSS 分支里 tanh(s/phi) 完全取代了 sign(s)这就是边界层方法。s 的物理单位是 1/stanh 的输入无量纲所以 phi 要和 s 保持同一个量级一般取 0.01 到 0.1。下面是状态方程函数。它接收当前状态 X、指令加速度 a_cmd返回状态导数function dX plant(X, a_cmd, Vm, Vt) thm X(3); tht X(6); dX [Vm * cos(thm); Vm * sin(thm); a_cmd / Vm; % 航迹角速率由法向过载给出 Vt * cos(tht); Vt * sin(tht); 0]; end固定目标 Vt0所以第 4、5、6 个分量都是 0。追逐状态目标时把 Vt 和 θt 设成目标的速度和航向即可同一个 plant 函数不需要改。3.3 主循环RK4 积分与脱靶判定主程序把制导指令计算和内层积分串起来。每个仿真步长内先根据当前状态算出 a_cmd然后以这个 a_cmd 为输入做四步 RK4。这个处理方式在工程上等价于“制导指令零阶保持”即认为一个计算周期内指令不变符合数字制导计算机的实际工作方式。% 攻击角约束制导律仿真主程序 clear; clc; close all; Vm 300; Vt 0; X0 [0; 0; deg2rad(25); 6000; 0; deg2rad(0)]; theta_d deg2rad(-20); param.N1 4; % 比例项导航比 param.N2 2.5; % 偏置项增益 param.lambda_s 4; % 滑模面系数 param.eps_s 30; % 切换增益 param.phi 0.05; % 边界层厚度 mode VSS; % 或 PN_ANG dt 0.001; t_max 30; X X0; t 0; r_min 1e6; hist.t[]; hist.X[]; hist.a[]; hist.q[]; while t t_max xmX(1); ymX(2); thmX(3); xtX(4); ytX(5); thtX(6); dx xt - xm; dy yt - ym; r sqrt(dx^2 dy^2); if r 0.5, break; end q atan2(dy, dx); qdot (dx*(Vt*sin(tht)-Vm*sin(thm)) ... - dy*(Vt*cos(tht)-Vm*cos(thm))) / r^2; rdot Vt*cos(tht-q) - Vm*cos(thm-q); Vc -rdot; tgo r / max(Vc, 1); a_cmd guidance_command(q, qdot, r, Vc, theta_d, tgo, mode, param); % RK4 更新 k1 plant(X, a_cmd, Vm, Vt); k2 plant(X0.5*dt*k1, a_cmd, Vm, Vt); k3 plant(X0.5*dt*k2, a_cmd, Vm, Vt); k4 plant(Xdt*k3, a_cmd, Vm, Vt); X X dt/6*(k1 2*k2 2*k3 k4); t t dt; hist.t(end1)t; hist.X(:,end1)X; hist.a(end1)a_cmd; hist.q(end1)q; end theta_m_final X(3); err_angle rad2deg(theta_m_final - theta_d); fprintf(脱靶量 r_min%.4f m, 终端攻角误差%.4f deg\n, r, err_angle);qdot 的计算式用的是相对坐标求导的直接展开不需要额外定义视线坐标系代码简单且不容易出现符号错误。tgo 用 max(Vc, 1) 是防除零Vc 接近 0 时仿真早已经在脱靶判定之前跳出了。r 0.5 作为命中判定条件意味着脱靶量小于半米就认为命中。3.4 结果可视化与参数对应关系仿真结束后至少要看两条曲线导弹飞行轨迹以及终端攻击角误差随时间的变化。轨迹用 plot 画 Y 随 X 变化攻击角误差用 rad2deg(theta_d) 归一化后看末段是否收敛到 0。还有一个很重要的检查对象是 a_cmd 曲线变结构制导律如果出现沿零线附近高频抖振大多能在 a_cmd 历史曲线里一眼看出。固定目标场景里命中时刻指令加速度 a_cmd 大概率不为零。原因很简单制导律还在抑制 qdot 的残差直到最后一刻导弹都在微调航迹方向。这是正常现象不代表程序有 bug。4. 攻击角约束制导律仿真的参数设定与坑排查4.1 三个必须调的参数与判定依据带攻击角约束的比例导引里最需要调的是 N1 和 N2。N1 影响视线角速率收敛速度N1 太大指令饱和N1 太小飞行轨迹弯曲度不够N2 直接决定攻击角约束的强度。判断依据很简单看终端航迹角相对期望值的偏差以及脱靶量是否被第二项破坏。先固定 N14把 N2 从 1 扫到 4通常能在某个值附近同时满足脱靶量小于 0.5 米、攻击角偏差小于 0.5 度。变结构制导律要调三个量λ、ε₁、φ。λ 决定滑动面收敛快慢可以参考比例导引中 N2/N1 的比例关系起调。ε₁ 必须大于目标机动和模型不确定性造成的等效扰动上界否则滑动面无法建立攻击角误差无法消除。φ 控制抖振和稳态精度之间的平衡φ0.01 时指令加速度会出现明显的锯齿φ0.1 时曲线平滑但终端角度误差可能偏大。这三个参数没有固定解析公式工程上都是先定 λ再加 ε₁最后用 φ 磨平曲线。4.2 脱靶量与攻击角误差的取舍两个指标在同一个仿真里经常是矛盾的。偏置项过强导弹会在末端绕一个很大的弧线去“够”期望角度轨迹弯过头脱靶量放大偏置项太弱攻角约束满足不了。碰到这种情况先检查是不是指令饱和造成的而不是盲目调整 N2 或 ε₁。用 max(a_cmd, 可用过载) 做饱和限幅是更贴近实际的改法amax 100; % 最大法向加速度 m/s^2 a_cmd max(min(a_cmd, amax), -amax);限幅后终端攻击角误差增大是预期结果说明导弹可用过载不足以支撑期望弹道的弯曲。真实系统中这里对应的是气动舵面偏转限制不能因为仿真没限幅就忽略。限幅之后还可以把制导律切换为“先比例、末端角度修正”的分段策略保证飞行中段有足够的能量。4.3 仿真不收敛或曲线发散的常见原因第一个原因是单位不统一。atan2 返回弧度如果 theta_d 用角度直接代入两个量在数值上差了约 57.3 倍仿真结果会乱到无法解释。脚本里统一用 deg2rad 作为唯一入口绘图输出时才转回角度。第二个原因是 tgo 估计失真。接近速度 Vc 在固定目标场景末段未必单调趋近于 Vm如果目标有速度Vc 可能过零直接用 r/Vc 会算出无穷大的偏置项。用 max(Vc, 1) 或者加一个下限保护是好习惯同时要在循环开头检查 r 是否已经小于脱靶阈值。第三个原因是积分方法跟不上抖振。变结构制导律里 φ 取得过小时tanh(s/φ) 接近 sign(s)指令加速度在一个仿真步长内急剧翻转。此时 RK4 的表现只比欧拉好一点彻底解决要靠“保留滑模切换项但外加低通滤波”的工程处理。常见做法是给 a_cmd 加一阶惯性环节a_cmd_filt a_cmd_filt dt/tau * (a_cmd - a_cmd_filt);τ 取 0.02 到 0.05 秒相当于模拟制导系统本身的动态响应曲线会平滑很多而攻击角约束精度损失通常小于 0.1 度。这一步在论文仿真中可以省略但在考虑工程可实现性时必须保留。最后排查一种容易被忽略的现象仿真在 r 0.5 之前就提前跳出原因是目标初始位置或速度配置导致弹道越过目标但没有满足接近条件。这时不用急着调参数先画出轨迹看看是不是初始航向角与目标方向相差过大末端弹道绕过了目标后再折返。制导系统的捕获区有限初始航向不能偏离视线方向太远这是物理约束不是程序问题。5. 批量仿真验证把攻击角约束结果与参数矩阵同时沉淀下来验证一条制导律是否真的“带攻击角约束”单跑一条弹道说服力不够至少要扫一遍期望攻击角的覆盖范围。把上一章的主循环包成函数输出脱靶量、终端航迹角和指令加速度序列然后在主脚本里循环不同 θd、不同初始航向角。这个动作既能把整个数据变成表格也顺便暴露参数矩阵里那些单条弹道看不见的边界情况。theta_d_list deg2rad(-40:10:20); % 期望攻击角扫描 N numel(theta_d_list); results table(Size,[N,4], ... VariableTypes, {double,double,double,double}, ... VariableNames, {theta_d_deg,miss_m,angle_err_deg,max_a} ); for i 1:N [miss, err_deg, a_hist] run_guidance_scenario(theta_d_list(i)); results.theta_d_deg(i) rad2deg(theta_d_list(i)); results.miss_m(i) miss; results.angle_err_deg(i) err_deg; results.max_a(i) max(abs(a_hist)); end disp(results);run_guidance_scenario 把上一章从“参数初始化”到“脱靶判定”的全部代码收进去返回脱靶量、终端攻击角误差、指令加速度曲线。扫描结果里通常能看到明显边界某个期望攻击角和初始几何之间一旦超过一定角度差攻击角误差开始快速增大而脱靶量依然很小。这个边界值就是该制导律的可用捕获区写报告或答辩时比单条弹道光标图有用得多。更进一步扫描 N2 或 ε₁ 与 θd 的二维网格用 surf 或 imagesc 画“攻击角误差随参数和期望角变化”的等高线图可以直观标出参数可用区间。这个做法对不同版本的 MATLAB 都通用不涉及任何工具箱R2023b、R2026a、R2026b 安装包装好之后直接跑脚本即可。有人问现在带 AI 能力的代码生成工具能不能像执行 Python 一样直接操作 MATLAB 任务批量跑这种参数扫描自然可以但制导律的物理边界、指令限幅和积分步长设置仍然需要有人把关这正是写清楚主循环结构的意义所在。跑完批量扫描后建议把 results 表存成 CSV用 MATLAB 自带的 writetable 导出后续画图、写文档、做对比分析都可以直接复用同一份数据。本文还有配套的精品资源点击获取