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

引言

脉冲雷达的工作原理在教材里通常以抽象公式呈现：雷达方程、虚警概率、检测概率、Swerling 起伏模型……但要把这些公式连成一个能跑起来、能看见的仿真并不容易。本文基于 AN/MPQ-64 "哨兵"雷达的公开规格，完整实现了一个脉冲雷达系统仿真——从雷达方程推导最大探测距离，到多目标运动与起伏回波的蒙特卡洛模拟，再到 PPI / A-Scope / 真实位置三屏联动可视化，所有代码打包为一个自包含的纯基础 MATLAB 仿真脚本，不依赖 Phased Array System Toolbox，可在任意 MATLAB 版本运行。

仿真的设计目标是回答下面三个问题：

  1. 给定雷达参数（功率、增益、波束宽度、虚警/检测概率），理论最大探测距离是多少？SNR 沿距离轴如何衰减？
  2. 当多个目标以随机 RCS、随机位置朝向雷达运动时，哪些目标会被检测到？哪些会落入盲区或发生距离模糊？
  3. 能否用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 模型。
  • 目标重置：当Rk>1.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