MATLAB实战（23）：雷达脉冲系统仿真：从参数分析到多目标检测可视化

引言脉冲雷达的工作原理在教材里通常以抽象公式呈现：雷达方程、虚警概率、检测概率、Swerling 起伏模型……但要把这些公式连成一个能跑起来、能看见的仿真并不容易。本文基于 AN/MPQ-64 "哨兵"雷达的公开规格，完整实现了一个脉冲雷达系统仿真——从雷达方程推导最大探测距离，到多目标运动与起伏回波的蒙特卡洛模拟，再到 PPI / A-Scope / 真实位置三屏联动可视化，所有代码打包为一个自包含的纯基础 MATLAB 仿真脚本，不依赖 Phased Array System Toolbox，可在任意 MATLAB 版本运行。仿真的设计目标是回答下面三个问题：给定雷达参数（功率、增益、波束宽度、虚警/检测概率），理论最大探测距离是多少？SNR 沿距离轴如何衰减？当多个目标以随机 RCS、随机位置朝向雷达运动时，哪些目标会被检测到？哪些会落入盲区或发生距离模糊？能否用PPI / A-Scope / 真实位置三屏联动把检测过程直观呈现？本文按这三条主线展开。所有数值结果均来自仿真脚本sentinel_radar_animation.m的实际运行。一、系统整体设计1.1 参数配置仿真采用 AN/MPQ-64 "哨兵"雷达的公开规格（X 波段、3D 监视雷达），核心参数如下：参数数值单位说明峰值功率PdgP_{dg}Pdg​10510^5105W发射脉冲峰值功率天线增益GantG_{ant}Gant​30dB（等效 1000）天线对信号的放大能力工作频率fff9GHzX 波段脉冲宽度τ\tauτ1μs单个脉冲持续时间系统损耗LsL_sLs​5dB系统各环节信号衰减噪声系数FnF_nFn​4dB接收机内部噪声水平系统温度T0T_0T0​290K标准噪声温度平均 RCSσ\sigmaσ1m²目标平均雷达截面积虚警概率PfaP_{fa}Pfa​10−610^{-6}10−6—噪声误判为目标的概率上限检测概率PdP_dPd​0.9—正确探测到目标的概率波束宽度θbw\theta_{bw}θbw​3度天线瞬时覆盖角度仿真得到的关键派生参数（运行sentinel_radar_animation.m第 1–2 节可得）：派生参数数值说明波长λ\lambdaλ0.0333 mc/fc/fc/f距离分辨率RresR_{res}Rres​149.90 mcτ/2c\tau/2cτ/2最小可测距离RminR_{min}Rmin​149.90 m与RresR_{res}Rres​相等所需最小 SNR20.59 dBAlbersheim 近似噪声功率N0N_0N0​1.01×10−141.01\times10^{-14}1.01×10−14WkT0FnBnkT_0F_nB_nkT0​Fn​Bn​最大探测距离RmaxR_{max}Rmax​11.13 kmSNR 达到阈值无模糊距离Rmax,naR_{max,na}Rmax,na​14.84 km(4/3)Rmax(4/3)R_{max}(4/3)Rmax​PRF10.10 kHzc/(2Rmax,na)c/(2R_{max,na})c/(2Rmax,na​)PRI0.099 ms1/PRF1/\text{PRF}1/PRF注：shnidman函数（Phased Array Toolbox 提供）可直接用于计算 SNR 需求；本脚本在没有工具箱时改用 Albersheim 近似（N=1N=1N=1、非起伏目标），两者数值偏差 0.5 dB。1.2 检测场景仿真同时运行 5 个独立移动目标，初始参数随机：初始位置：极坐标均匀采样。径向距离在[Rmin+0.3Rmax,na, 1.3Rmax,na][R_{min}+0.3R_{max,na},\;1.3R_{max,na}][Rmin​+0.3Rmax,na​,1.3Rmax,na​]之间随机，部分在最大探测距离内、部分在外，同时体现"盲区""距离模糊"等典型现象。运动方向：均径直朝向雷达原点。速度：远距离目标 600–900 m/s，近距离目标 400–650 m/s，确保 15 s 仿真时长内有充分穿越。RCS 起伏：每个脉冲用exprnd(σ)\text{exprnd}(\sigma)exprnd(σ)采样瞬时 RCS，等效于Swerling I/II 模型。目标重置：当Rk1.5×Rmax,naR_k 1.5 \times R_{max,na}Rk​1.5×Rmax,na​时随机重置位置和速度，并清空该目标的航迹缓冲。二、核心算法2.1 雷达方程与最大探测距离由单基地雷达方程：Rmax=[Pdg Gant2 λ2 σ(4π)3 Ls N0 SNRreq]1/4 R_{max} = \left[\frac{P_{dg}\,G_{ant}^2\,\lambda^2\,\sigma}{(4\pi)^3\,L_s\,N_0\,\text{SNR}_{req}}\right]^{1/4}Rmax​=[(4π)3Ls​N0​SNRreq​Pdg​Gant2​λ2σ​]1/4代入参数得Rmax=11.13R_{max} = 11.13Rmax​=11.13km。MATLAB 代码片段：SNR_req_dB=A_alb+0.12*A_alb*B_alb+1.7*B_alb;SNR_req=10^(SNR_req_dB/10);N0=k_B*T0*Fn*Bn;Rmax=(Pdg*Gant^2*lambda^2