2.1 引言数学与物理基础的核心地位雷达仿真的本质是在数字世界中重构电磁波与物理世界相互作用的复杂过程。这一重构的基石并非编程技巧或软件工程而是对雷达系统背后深刻的数学原理和物理定律的精确把握。一个高保真的雷达仿真系统其灵魂在于其内在的数学模型。这些模型决定了仿真系统能否真实地反映雷达的探测能力、测量精度以及在复杂环境中的生存能力。因此在着手构建任何仿真模块之前必须首先夯实从雷达方程到信号模型的完整理论体系。本部分将系统性地阐述雷达仿真所依赖的核心数学与物理基础。我们将从最基本的雷达方程出发揭示决定雷达探测性能的根本因素然后深入探讨雷达信号的数学表征建立贯穿整个仿真链路的信号模型最后简要介绍随机过程与噪声模型为后续处理杂波、干扰和信号检测奠定基础。理解这些内容是读懂后续所有建模细节的前提。2.2 雷达方程探测能力的根本法则雷达方程Radar Range Equation是连接雷达系统参数、目标特性、环境因素与最终探测性能的桥梁。它是所有雷达仿真中评估信噪比SNR和探测距离的起点。2.2.1 基本雷达方程单基地、点目标对于一部发射天线和接收天线共址单基地的雷达探测一个远场点目标时其接收到的回波信号功率 Pr​可由以下经典方程描述Pr​4πR2Pt​Gt​​⋅4πR2σ​⋅LAe​​参数释义Pt​雷达发射机的峰值功率瓦W。Gt​发射天线在目标方向上的增益无量纲。它表示天线将辐射能量集中到特定方向的能力。R雷达与目标之间的斜距米m。方程中出现两次 R2分别对应发射路径的功率密度衰减和接收路径的功率密度衰减合起来就是 R4关系这直观地说明了距离对雷达探测能力的巨大影响。σ目标的雷达散射截面积Radar Cross Section, RCS单位平方米m²。它是一个等效面积表征目标截获入射功率并将其各向同性重辐射回雷达的能力。RCS是目标形状、尺寸、材料、表面处理、姿态角以及雷达波长的复杂函数是雷达仿真中目标建模的核心。Ae​接收天线的有效孔径面积平方米m²。对于同一个天线既发射又接收的情况通常有 Ae​4πGr​λ2​其中 Gr​是接收天线增益通常 Gr​Gt​λ是雷达工作波长。L系统总损耗因子1无量纲。这是一个综合性参数包含了发射馈线损耗、接收馈线损耗、大气传播损耗、信号处理损耗、波束形状损耗等所有非理想因素造成的功率损失。将 Ae​4πGλ2​假设 Gt​Gr​G代入上式得到更常见的形式Pr​(4π)3R4LPt​G2λ2σ​(1)2.2.2 雷达方程的扩展形式基本方程描述了单个脉冲照射点目标的情况。在实际仿真中我们需要考虑更多因素积累脉冲的影响雷达通常对同一目标发射多个脉冲并进行相干或非相干积累以提高SNR。对于 N个脉冲的非相干积累平均信噪比改善因子约为 N在检测概率较高时。方程可修正为(SNR)N​(4π)3R4kTs​Bn​Fn​LPt​G2λ2σN​其中分母部分引入了接收机噪声见下文。波束驻留与扫描的影响对于机械扫描雷达天线波束扫过目标需要一定时间驻留时间 Td​。在此期间内雷达发射的脉冲数 NTd​⋅PRFPRF为脉冲重复频率。波束形状也会造成信号幅度调制引入“波束形状损耗”。分布式目标与杂波对于体杂波如雨或面杂波如地、海目标不再是点而是分布在某个分辨率单元由距离门宽度 ΔR和波束宽度 θ决定内的无数散射体。此时雷达方程中的 σ需要替换为单位体积的杂波散射系数 η对于体杂波或单位面积的杂波散射系数 σ0对于面杂波并乘以分辨率单元的体/面积。这是杂波功率计算的基础。2.2.3 信噪比SNR与探测距离雷达的探测性能最终取决于接收机输出端的信噪比。接收机噪声功率 Pn​可表示为Pn​kTs​Bn​k玻尔兹曼常数1.38×10−23J/K。Ts​系统噪声温度K包括天线噪声和接收机内部噪声。Bn​接收机噪声带宽Hz通常近似为匹配滤波器的带宽。因此单个脉冲的信噪比为SNRPn​Pr​​(4π)3R4kTs​Bn​LPt​G2λ2σ​(2)雷达能够可靠检测目标所需的最小信噪比称为检测门限信噪比 (SNR)min​。它由期望的检测概率 Pd​、虚警概率 Pfa​以及信号处理方式决定。令 SNR(SNR)min​可以从方程(2)解出雷达的最大探测距离 Rmax​Rmax​[(4π)3kTs​Bn​L(SNR)min​Pt​G2λ2σ​]1/4(3)这就是著名的雷达作用距离方程。它清晰地揭示了“四次方根”关系欲使探测距离加倍需要将发射功率提高16倍或通过其他参数的组合实现等效的16倍改善。在仿真中我们通常利用方程(2)计算特定距离下的SNR进而评估检测性能。图2-1雷达方程中各参数关系示意图[发射机] | v Pt (功率) | v Gt (天线增益) | v 辐射功率密度: Pt*Gt/(4πR1²) | v [自由空间传播] | 距离 R1 v 照射到目标的功率密度 | v σ (目标RCS) | v 目标散射功率 (各向同性) | v [自由空间传播] | 距离 R2 (R2R1R) v 到达接收天线的功率密度 | v Ae (有效孔径) | v Pr (接收功率) | v ┌──────────────┴──────────────┐ v v Pn (噪声功率) L (系统损耗) | | v v SNR Pr / (Pn*L) --- 与检测门限比较 --- 检测判决2.3 雷达信号的数学表征雷达通过发射特定形式的电磁波信号并分析其回波来提取信息。因此对雷达信号的数学描述是仿真中生成、处理和解析信号的基石。2.3.1 复信号表示法实际雷达处理的信号是实信号但采用复信号解析信号表示法可以极大地简化数学分析特别是在处理幅度和相位信息时。一个实带通信号 sr​(t)中心频率为 f0​可以表示为sr​(t)a(t)cos[2πf0​tϕ(t)]其中 a(t)是时变包络ϕ(t)是时变相位。其对应的复解析信号 s(t)定义为s(t)a(t)ej[2πf0​tϕ(t)]u(t)ej2πf0​t其中 u(t)a(t)ejϕ(t)称为信号的复包络​ 或基带信号。它包含了信号的所有调制信息幅度和相位而 ej2πf0​t代表载频。在仿真和数字信号处理中我们通常直接对复包络 u(t)进行采样和运算这相当于将信号下变频到了基带零中频。接收机中的正交解调I/Q解调正是为了获得信号的复包络。2.3.2 常见雷达波形及其数学模型雷达波形决定了距离分辨率、多普勒分辨率、模糊特性、抗干扰能力等关键性能。以下是几种基本波形的数学模型简单脉冲Simple Pulse:u(t)A⋅rect(τt​)其中 τ是脉冲宽度rect(t/τ)是矩形函数在 ∣t∣≤τ/2时为1否则为0。A是幅度。其距离分辨率 ΔRcτ/2c为光速。线性调频脉冲LFM, Chirp:u(t)A⋅rect(τt​)⋅ejπβ(t/τ)tA⋅rect(τt​)⋅ejπkt2其中 β是调频带宽kβ/τ是调频斜率。LFM脉冲通过脉冲压缩技术可以在保持较长脉冲宽度高能量的同时获得由带宽 β决定的高距离分辨率 (ΔRc/(2β))。相位编码脉冲Phase-Coded Pulse:将一个长脉冲划分为 N个宽度为 τc​的子脉冲码片每个码片赋予不同的相位通常为0或π对应BPSK或更多相位如QPSK。u(t)An0∑N−1​cn​⋅p(t−nτc​)其中 cn​ejϕn​是相位编码序列如巴克码、m序列p(t)是码片波形通常为矩形。相位编码也用于脉冲压缩。连续波CW与调频连续波FMCW:CW: s(t)Aej2πf0​t用于测速多普勒效应。FMCW: 发射频率随时间线性变化的连续波。其瞬时频率 f(t)f0​kt。通过比较接收信号与发射信号的瞬时频率差差拍频率 fb​可以同时测距 (R∝fb​) 和测速。2.3.3 雷达回波信号模型考虑一个点目标雷达发射信号复包络为 ut​(t)。经过时延 τ2R/c、多普勒频移 fd​2vr​/λvr​为径向速度以及由目标RCS和传播路径决定的幅度衰减 A后接收到的基带回波信号复包络 ur​(t)可建模为ur​(t)A⋅ut​(t−τ)⋅ej2πfd​(t−τ)n(t)其中 n(t)是加性噪声包括热噪声、杂波等。幅度 A包含了雷达方程中除时延和多普勒外的所有因素A(4π)3R4LPt​G2λ2σ​​⋅Pn​​1​(归一化到噪声电平)在仿真中我们通常根据雷达方程计算接收功率 Pr​然后结合噪声功率 Pn​得到信噪比 SNR进而确定回波信号的幅度 ASNR​假设噪声功率归一化为1。同时时延 τ和多普勒 fd​是根据目标运动轨迹实时计算的。对于扩展目标或多个目标回波信号是多个上述信号的叠加。对于杂波则是大量散射体回波的相干或非相干叠加。2.4 随机过程与噪声模型雷达接收的信号总是伴随着噪声和杂波它们本质上是随机的。用概率统计的方法描述这些随机过程是仿真中生成逼真环境信号和进行统计性能评估的关键。2.4.1 热噪声模型接收机前端产生的热噪声通常建模为加性高斯白噪声AWGN。其特性如下加性Additive与信号线性叠加。白White功率谱密度在信号带宽内是平坦的。高斯Gaussian瞬时幅度服从零均值的高斯正态分布。复基带AWGN n(t)nI​(t)jnQ​(t)其中同相分量 nI​(t)和正交分量 nQ​(t)是相互独立的零均值高斯随机过程其方差功率为 N0​/2单边功率谱密度为 N0​。在离散时间仿真中我们通过生成两个独立的高斯随机数序列来模拟I、Q两路噪声。2.4.2 杂波统计模型杂波地、海、气象的幅度统计特性远非高斯它严重依赖于分辨率、擦地角、极化方式等因素。常用的统计模型有瑞利Rayleigh分布适用于低分辨率、低擦地角下的均匀杂波如远距离地面杂波。其幅度分布概率密度函数PDF为p(a)σ2a​exp(−2σ2a2​),a≥0其中 σ2是平均功率。瑞利分布是复高斯随机过程的包络分布。韦布尔Weibull分布比瑞利分布更灵活能描述中等非均匀性的杂波如起伏地形。p(a)αη​(αa​)η−1exp[−(αa​)η],a≥0其中 η是形状参数η2时退化为瑞利α是尺度参数。对数正态Log-Normal分布适用于非常不均匀的杂波如城市环境或高分辨率下的海杂波尖峰。p(a)as2π​1​exp[−2s2(lna−m)2​],a0其中 m和 s分别是 lna的均值和标准差。K分布目前公认的对高分辨率海杂波建模最有效的模型。它可以用一个乘积模型来描述zx⋅y​其中 x服从伽马分布表示局部平均功率的起伏y服从指数分布表示给定局部功率下的散斑分量。K分布能很好地模拟海杂波的尖峰和长拖尾特性。在仿真中需要根据场景选择合适的杂波模型并生成符合该幅度分布和特定功率谱相关特性的随机序列。2.4.3 目标起伏模型Swerling模型目标的RCS并非恒定而是随时间起伏的。Swerling模型将目标起伏分为四种经典类型Swerling I II适用于由大量独立散射体组成的目标如飞机。起伏服从瑞利分布。I型慢起伏在一次扫描期间RCS不变扫描间变化II型快起伏脉冲间独立变化。Swerling III IV适用于一个主散射体加多个小散射体的目标。起伏服从4自由度的卡方分布。III型慢起伏IV型快起伏。在仿真中需要根据目标类型在每个相干处理间隔CPI或每个脉冲按照相应的分布生成RCS的随机样本用于计算回波幅度 A。2.5 小结与下部分预告本部分建立了雷达仿真的核心理论框架。雷达方程提供了从系统参数到接收信噪比的全局映射是性能评估的标尺。信号模型复信号表示、波形数学描述、回波模型构成了仿真中信号生成、传播、接收和处理的全链路数学语言。随机过程模型噪声、杂波、目标起伏则为仿真注入了现实世界的不确定性和统计特性是评估雷达在复杂环境下性能的基础。掌握这些基础后我们便具备了深入雷达系统内部进行建模的“词汇”和“语法”。在下一部分我们将开始深入雷达系统的第一个关键子系统雷达天线与波束形成的建模与仿真。我们将详细讨论天线方向图函数、相控阵天线波束扫描与形成的原理、数字波束形成DBF算法以及它们在仿真中的实现方法。