5.1 引言：目标——雷达信息的来源
在雷达系统中，目标不仅是探测的对象，更是信息的载体。雷达通过发射电磁波并接收目标散射的回波，从中提取目标的距离、速度、角度以及形状、尺寸、材料等特征信息。因此，对目标电磁散射特性的准确建模，是雷达仿真能否真实反映探测性能、识别能力和抗干扰效果的关键。一个高保真的目标模型，需要能够模拟雷达波与复杂物体相互作用产生的幅度、相位、极化和时频特性变化。
现实世界中的雷达目标千差万别，从简单的金属球体到复杂的飞机、舰船、车辆，乃至自然地貌（如山体、海面）。它们的散射特性随频率、极化、姿态角（视角）剧烈变化。在仿真中，我们无法也无需对每个目标的微观电磁过程进行全波仿真，而是需要建立不同层次、不同保真度的模型，以平衡计算复杂度和物理真实性。本部分将系统阐述雷达目标散射的基本原理、雷达截面积（RCS）的定义与统计特性、经典的目标起伏模型，并深入探讨复杂目标RCS的建模方法与仿真实现，最后简要介绍电磁波传播的大气效应。
5.2 雷达截面积（RCS）基础
5.2.1 RCS的定义与物理意义
雷达截面积（Radar Cross Section, RCS），通常用符号 σ表示，是衡量目标向雷达方向散射电磁波能力的物理量。其严格定义基于平面波照射下的远场条件：
σ=R→∞lim4πR2∣Ei∣2∣Es∣2
其中：
R是雷达与目标之间的距离。
Ei是目标处入射电场的强度。
Es是雷达接收机处散射电场的强度。
物理意义：RCS可以理解为目标等效的一个各向同性散射体的截面积，该散射体将入射功率均匀地散射到全空间，并在雷达处产生与实际目标相同的功率密度。RCS的单位是平方米（m²），但常用分贝平方米（dBsm）表示：σdBsm=10log10(σm2)。
5.2.2 RCS的频率特性与极化特性
频率特性：根据目标尺寸 L与波长 λ的相对关系，散射机制可分为三个区域：
瑞利区（L≪λ）：RCS与频率的4次方成正比（f4），散射较弱。适用于低频对小型目标的探测分析。
谐振区（L∼λ）：RCS随频率剧烈振荡，是目标固有谐振模式的体现。复杂目标的RCS在此区域难以预测。
光学区（L≫λ）：RCS趋于稳定，与波长无关，主要由目标的几何形状决定。大多数军用雷达和目标处于此区域。
极化特性：散射矩阵 S描述了入射波与散射波极化状态之间的关系：
[EhsEvs]=Re−jkR[ShhSvhShvSvv][EhiEvi]
其中下标 h,v分别代表水平极化和垂直极化。对于互易介质，有 Shv=Svh。极化信息是目标识别的重要特征。
5.2.3 简单形状目标的RCS
一些简单形状目标的RCS有解析解，是理解和验证模型的基础。
理想导电球体：在光学区，其RCS为 σ=πa2，其中 a为球半径，与视角无关。这是各向同性散射体。
平板：面积为 A的平板，当其法线方向对准雷达时，RCS最大：σmax=λ24πA2。当偏离法线时，RCS迅速下降。
圆柱体：半径为 a，长度为 L的圆柱体，当其轴线垂直于雷达视线时，RCS为 σ=λ2πaL2。
角反射器：由三个相互垂直的金属平板构成，能将入射波沿原路反射回去，产生巨大的RCS。二面角反射器的最大RCS为 σ=λ28πa2b2，三面角反射器则更大。常用于校准或作为诱饵。
5.3 目标起伏的统计模型（Swerling模型）
实际目标的RCS并非定值，由于目标相对雷达的运动（如翻滚、摆动）导致其各散射中心的相对相位变化，使得接收到的回波幅度随机起伏。Swerling模型是描述这种起伏统计特性的经典模型。
5.3.1 Swerling模型分类
Swerling模型根据起伏快慢（扫描间起伏/扫描内起伏）和概率密度函数（PDF）的形状，分为四类：
模型类型 | 起伏快慢 | RCS概率密度函数 p(σ) | 幅度（电压）分布 | 典型目标 |
|---|---|---|---|---|
Swerling I | 慢起伏（扫描间） | 指数分布：p(σ)=σav1e−σ/σav | 瑞利分布 | 由大量相似散射中心组成的复杂目标（如飞机） |
Swerling II | 快起伏（脉冲间） | 指数分布：p(σ)=σav1e−σ/σav | 瑞利分布 | 同上，但起伏更快 |
Swerling III | 慢起伏（扫描间） | χ2分布（4自由度）：p(σ)=σav24σe−2σ/σav | 一个主导散射中心加多个小散射中心的目标 | |
Swerling IV | 快起伏（脉冲间） | χ2分布（4自由度）：p(σ)=σav24σe−2σ/σav | 同上，但起伏更快 |
其中，σav是平均RCS。“扫描间”指在一次天线波束驻留时间内（多个脉冲），RCS保持不变，但扫描到扫描之间变化；“脉冲间”指每个脉冲的RCS都独立变化。
5.3.2 Swerling模型在仿真中的实现
在雷达信号仿真中，我们需要为每个目标（或每个脉冲）生成一个随机的RCS值 σ，使其符合指定的Swerling分布。
生成瑞利分布幅度：对于Swerling I/II型，回波电压幅度 v服从瑞利分布。可以先生成两个独立的标准正态分布随机数 x,y∼N(0,1)，则 v=x2+y2服从均值为0、方差为1的瑞利分布。为了得到平均RCS为 σav的幅度，需要设置比例因子，因为RCS正比于幅度平方 σ∝v2。具体地，令 v=−σavln(u)，其中 u是[0,1]上的均匀分布随机数，则 σ=v2服从指数分布。
生成Swerling III/IV型RCS：可以生成两个独立的指数分布随机变量 σ1,σ2（均值为 σav/2），然后求和 σ=σ1+σ2，则 σ服从自由度为4的 χ2分布。
关联性处理：对于慢起伏模型（I, III），在一次相干处理间隔（CPI）或一次扫描内，所有脉冲的RCS应使用相同的随机值。对于快起伏模型（II, IV），每个脉冲的RCS应独立生成。
5.3.3 起伏模型对检测性能的影响
Swerling起伏目标的检测概率 Pd通常低于非起伏（Marcum）目标。在给定虚警概率 Pfa和信噪比（SNR）下，起伏目标的检测需要更高的SNR才能达到相同的 Pd。雷达信号处理中的恒虚警率（CFAR）检测器设计必须考虑目标的起伏特性。
5.4 复杂目标RCS的建模方法
对于飞机、舰船等复杂目标，其RCS是姿态角（方位角 ϕ、俯仰角 θ）的复杂函数，呈现出剧烈的闪烁特性。在仿真中，主要有以下几种建模方法：
5.4.1 查找表法（RCS Pattern）
这是最直接和常用的方法。
原理：通过电磁计算软件（如FEKO, HFSS）或外场实测，获取目标在全方位角、全俯仰角以及不同频率、极化下的RCS数据，形成一个多维数据表（查找表）。
数据结构：通常是一个以方位角 ϕ、俯仰角 θ为索引的二维数组 σ(ϕ,θ)。对于宽带或全极化仿真，可能需要增加频率维和极化维。
仿真使用：在仿真运行时，根据目标相对于雷达的瞬时视角 (ϕ,θ)，通过双线性插值从查找表中获取对应的RCS值 σ。
优点：精度高，能反映复杂的RCS闪烁特性。
缺点：数据获取成本高（计算或测量），数据量大，且只能用于已建模的特定目标。
5.4.2 部件分解法（CSC：Component Scattering Center）
该方法将复杂目标分解为多个简单的几何形状（如平板、圆柱、球、角反射器等），每个部件用一个“散射中心”来表征。
原理：每个散射中心 i具有以下属性：
位置 (xi,yi,zi)（在目标坐标系中）。
频率依赖关系 σi∝fαi，其中 αi是频率指数（如平板 α=2，圆柱 α=1，尖锥 α=0）。
姿态角依赖关系：通常用一个方向图函数 Fi(ϕ,θ)来描述（如平板的sinc函数，圆柱的余弦函数）。
极化散射矩阵 Si。
总RCS计算：目标的总散射场是各散射中心散射场的矢量和（考虑相位干涉）：
Etotals=i=1∑NEis⋅exp(−j2kk^⋅ri)
其中 k^是雷达视线方向的单位矢量，ri是第 i个散射中心的位置矢量，k=2π/λ。总RCS为 σtotal=limR→∞4πR2∣Etotals∣2/∣Ei∣2。
优点：物理意义清晰，数据量小，易于参数化修改（如改变部件尺寸、位置），能模拟目标运动时RCS的闪烁。
缺点：忽略了部件间的多次散射和遮挡效应，精度低于全波计算。
5.4.3 几何绕射理论/物理绕射理论（GTD/PTD）
这是一种高频近似方法，适用于电大尺寸目标（L≫λ）。
原理：将目标的散射归结为几种典型散射机理的贡献之和，包括：
镜面反射：来自平滑表面的法向反射。
边缘绕射：来自边缘的不连续处，如机翼前缘、尾翼边缘。
尖顶绕射：来自锥尖、角点。
曲面爬行波：沿光滑曲面传播的波。
仿真实现：需要目标的几何模型（CAD模型），通过射线追踪算法确定雷达波照射到目标后产生的反射点、绕射点等，然后根据GTD/PTD公式计算每个散射源的贡献，并相干叠加。
优点：计算速度比全波法快得多，能处理电大尺寸目标，且能提供散射中心的位置信息，适用于高分辨率雷达（如ISAR）成像仿真。
缺点：在焦散区（如边缘的法线方向）失效，需要过渡函数修正。
5.4.4 计算机电磁学（CEM）全波仿真
包括矩量法（MoM）、有限元法（FEM）、时域有限差分法（FDTD）等。它们直接求解麦克斯韦方程组，理论上最精确。
应用：主要用于生成高精度的RCS查找表，或验证简化模型的准确性。由于计算量巨大，通常不用于实时或大规模场景的雷达系统仿真。
图5-1：复杂目标RCS建模方法对比与流程
┌─────────────────────────────────────────────────────────────┐ │ 复杂目标RCS建模方法选择与流程 │ ├─────────────────────────────────────────────────────────────┤ │ │ │ ┌─────────────┐ ┌─────────────┐ ┌─────────────┐ │ │ │ 建模需求分析 │ │ 目标几何信息 │ │ 可用计算资源 │ │ │ │ (精度/速度/ │ │ (CAD模型/ │ │ (CPU/内存/ │ │ │ │ 动态性) │ │ 简单描述) │ │ 时间预算) │ │ │ └──────┬──────┘ └──────┬──────┘ └──────┬──────┘ │ │ │ │ │ │ │ └──────────┬────────┴──────────┬───────┘ │ │ ▼ ▼ │ │ ┌─────────────────┐ ┌─────────────────┐ │ │ │ 高精度需求 │ │ 快速/参数化 │ │ │ │ (成像/识别) │ │ 需求 (系统仿真)│ │ │ └─────────┬───────┘ └────────┬────────┘ │ │ │ │ │ │ ┌───────▼───────┐ ┌───────▼───────┐ │ │ │ 全波仿真或 │ │ 部件分解法 │ │ │ │ 高精度测量 │ │ (CSC模型) │ │ │ │ (生成查找表) │ │ │ │ │ └───────┬───────┘ └───────┬───────┘ │ │ │ │ │ │ └──────────┬────────┘ │ │ ▼ │ │ ┌──────────────────┐ │ │ │ RCS数据库/模型 │ │ │ │ (σ(θ,φ,f,pol)) │ │ │ └─────────┬────────┘ │ │ │ │ │ ┌─────────▼────────┐ │ │ │ 仿真运行时调用 │ │ │ │ 1. 根据目标姿态 │ │ │ │ (θ,φ)插值查询 │ │ │ │ 2. 考虑频率/极化 │ │ │ │ 3. 应用起伏模型 │ │ │ │ (Swerling) │ │ │ │ 4. 计算回波幅度 │ │ │ └──────────────────┘ │ │ │ └─────────────────────────────────────────────────────────────┘5.5 目标动态特性与微多普勒效应
除了RCS起伏，目标的宏观运动（平动）和微观运动（振动、旋转）都会对回波信号产生调制，这些信息是目标识别的重要特征。
5.5.1 平动与宏观多普勒
目标相对于雷达的径向速度 vr会产生多普勒频移 fd=2vr/λ。在仿真中，这体现为回波信号的载频偏移或脉冲间的相位变化。
5.5.2 微动与微多普勒
目标或其部件（如螺旋桨、旋翼、履带）的周期性微动（振动、旋转、摆动）会对回波产生频率调制，在时频域（如短时傅里叶变换）形成独特的“微多普勒谱线”。
旋转部件：如直升机的旋翼。叶片尖端线速度很高，会产生很强的多普勒调制。假设叶片长度为 L，旋转角速度为 Ω，则叶片尖端的最大多普勒频移为 fd,max=±(2ΩL/λ)。回波信号中会包含以旋转频率为间隔的边带。
振动目标：如车辆发动机引起的车体振动。设振动幅度为 Dv，振动频率为 fv，则引入的相位调制为 ϕ(t)=(4πDv/λ)sin(2πfvt)，在频谱上会产生以 fv为间隔的边带。
仿真建模：需要在计算目标上每个散射中心的回波时，将其位置视为时间的函数 ri(t)，然后计算随时间变化的时延和多普勒频移。
5.6 电磁波传播的大气效应
雷达波在大气中传播时，会受到衰减、折射和湍流等效应的影响，在仿真中需要加以考虑。
5.6.1 大气衰减
主要由大气分子（氧气、水蒸气）吸收和散射（雨、雾、云）引起。衰减系数 γ(dB/km) 是频率、温度、压力和湿度的函数。
氧气吸收峰：在60 GHz附近。
水蒸气吸收峰：在22.2 GHz和183 GHz附近。
降雨衰减：与降雨率 R(mm/hr) 密切相关，经验模型如ITU-R P.838建议书。
仿真实现：对于距离 R处的目标，双程大气衰减因子为 Latm=exp(−2∫0Rγ(s)ds)。在简化模型中，可使用平均衰减系数 γˉ，则 Latm≈e−2γˉR。该因子将乘入雷达方程。
5.6.2 大气折射与多路径效应
折射：由于大气折射率随高度变化，雷达波传播路径发生弯曲，导致测高误差。可使用等效地球半径模型（4/3地球半径）进行近似修正。
多路径效应：低仰角探测时，雷达波可能经由直射路径和目标反射到地面的间接路径到达目标，两者在接收端干涉，导致仰角测量误差和信号起伏（闪烁）。建模较为复杂，需要计算直射波和反射波的路径差及反射系数。
5.6.3 大气湍流
导致信号幅度和相位的随机起伏（闪烁），影响跟踪精度。通常用对数正态分布或伽马分布来描述强度起伏。
5.7 目标与环境仿真模块设计
在雷达系统仿真框架中，目标与环境模块负责管理场景中所有目标的运动、散射特性，并计算它们对雷达信号的贡献。
图5-2：目标与环境仿真模块结构框图
┌─────────────────────────────────────────────────────────────┐ │ 目标与环境仿真模块接口定义 │ ├─────────────────────────────────────────────────────────────┤ │ 输入接口： │ │ 1. 场景配置：目标列表（类型、初始状态）、环境参数（大气、地面）│ │ 2. 雷达状态：雷达平台位置、速度、波束指向、发射信号参数 │ │ 3. 时序控制：仿真时间、脉冲重复间隔（PRI） │ ├─────────────────────────────────────────────────────────────┤ │ 内部核心模型： │ │ │ │ ┌────────────────────────────────────────────────────┐ │ │ │ 目标运动模型 │ │ │ │ 根据目标类型（匀速、机动、飞行器动力学）更新每个 │ │ │ │ 目标在全局坐标系中的位置、速度、姿态(θ,φ,ψ) │ │ │ └────────────────────────┬───────────────────────────┘ │ │ │ │ │ ┌────────────────────────▼───────────────────────────┐ │ │ │ 目标RCS模型 │ │ │ │ ┌──────────┐ ┌──────────┐ ┌──────────┐ │ │ │ │ │ RCS数据库 │ │ 部件分解 │ │ 简单形状 │ │ │ │ │ │ (查找表) │ │ (CSC) │ │ (解析式) │ │ │ │ │ └─────┬────┘ └────┬─────┘ └────┬─────┘ │ │ │ │ │ │ │ │ │ │ │ └─────┬──────┴──────┬──────┘ │ │ │ │ ▼ ▼ │ │ │ │ ┌─────────────────┐┌─────────────────┐ │ │ │ │ │ 姿态角计算 ││ RCS值获取/计算 │ │ │ │ │ │ (目标->雷达) ││ (考虑频率/极化) │ │ │ │ │ └─────────────────┘└─────────────────┘ │ │ │ └────────────────────────┬───────────────────────────┘ │ │ │ │ │ ┌────────────────────────▼───────────────────────────┐ │ │ │ 传播效应模型 │ │ │ │ 计算：1. 自由空间路径损耗 1/(4πR^2)^2 │ │ │ │ 2. 大气衰减 L_atm(R) │ │ │ │ 3. 多路径效应（如需要） │ │ │ └────────────────────────┬───────────────────────────┘ │ │ │ │ │ ┌────────────────────────▼───────────────────────────┐ │ │ │ 回波信号生成引擎 │ │ │ │ 对于每个目标i，每个阵元m（考虑相控阵）： │ │ │ │ 1. 计算时延 τ_i = 2R_i/c │ │ │ │ 2. 计算多普勒 f_di = 2(v_ri)/λ │ │ │ │ 3. 计算空间相位差（阵元间） Δφ_mi = k * (r_m·k_i) │ │ │ │ 4. 获取RCS起伏值 σ_i(t) (Swerling模型) │ │ │ │ 5. 计算回波幅度 A_i ∝ sqrt(σ_i) * L_prop │ │ │ │ 6. 生成阵元m接收到的目标i信号： │ │ │ │ s_mi(t) = A_i * s_tx(t-τ_i) * exp(j2πf_di t) * │ │ │ │ exp(jΔφ_mi) │ │ │ │ 7. 对所有目标求和：s_m(t) = Σ_i s_mi(t) │ │ │ └────────────────────────┬───────────────────────────┘ │ │ │ │ │ ┌────────────────────────▼───────────────────────────┐ │ │ │ 杂波与干扰模型 │ │ │ │ (将在后续章节详细展开) │ │ │ └────────────────────────────────────────────────────┘ │ ├─────────────────────────────────────────────────────────────┤ │ 输出接口： │ │ 每个接收阵元m的基带回波信号 s_m(t) （包含所有目标、杂波、│ │ 干扰的贡献） │ └─────────────────────────────────────────────────────────────┘关键实现考虑：
计算效率：当场景中目标数量众多时，逐个目标、逐个阵元的计算量巨大。需要采用空间分区（如网格化）、可见性判断（遮挡剔除）和近似算法来加速。
模型保真度选择：根据仿真目的（系统性能评估、信号级算法验证、成像仿真）选择合适的RCS模型（查找表、CSC、简单模型）。
动态RCS与姿态管理：对于机动目标，其姿态角连续变化，需要实时计算或插值RCS。对于螺旋桨飞机等目标，需要集成微动模型。
相干性保持：在整个相干处理间隔（CPI）内，需要保持目标回波相位的连续性，以正确模拟多普勒处理。
5.8 小结
本部分深入探讨了雷达目标特性与电磁散射建模。我们从雷达截面积（RCS）的基本定义出发，阐述了其频率、极化和统计特性（Swerling模型）。针对复杂目标，系统介绍了查找表法、部件分解法、GTD/PTD和高频仿真等多种建模方法的原理与适用场景。此外，我们还讨论了目标微动产生的微多普勒效应以及大气传播对雷达波的影响。最后，提出了一个集成化的目标与环境仿真模块设计方案。
