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^5105 | W | 发射脉冲峰值功率 |
| 天线增益GantG_{ant}Gant | 30 | dB（等效 1000） | 天线对信号的放大能力 |
| 工作频率fff | 9 | GHz | X 波段 |
| 脉冲宽度τ\tauτ | 1 | μs | 单个脉冲持续时间 |
| 系统损耗LsL_sLs | 5 | dB | 系统各环节信号衰减 |
| 噪声系数FnF_nFn | 4 | dB | 接收机内部噪声水平 |
| 系统温度T0T_0T0 | 290 | K | 标准噪声温度 |
| 平均 RCSσ\sigmaσ | 1 | m² | 目标平均雷达截面积 |
| 虚警概率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 m | c/fc/fc/f |
| 距离分辨率RresR_{res}Rres | 149.90 m | cτ/2c\tau/2cτ/2 |
| 最小可测距离RminR_{min}Rmin | 149.90 m | 与RresR_{res}Rres相等 |
| 所需最小 SNR | 20.59 dB | Albersheim 近似 |
| 噪声功率N0N_0N0 | 1.01×10−141.01\times10^{-14}1.01×10−14W | kT0FnBnkT_0F_nB_nkT0FnBn |
| 最大探测距离RmaxR_{max}Rmax | 11.13 km | SNR 达到阈值 |
| 无模糊距离Rmax,naR_{max,na}Rmax,na | 14.84 km | (4/3)Rmax(4/3)R_{max}(4/3)Rmax |
| PRF | 10.10 kHz | c/(2Rmax,na)c/(2R_{max,na})c/(2Rmax,na) |
| PRI | 0.099 ms | 1/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π)3LsN0SNRreqPdgGant2λ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