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表示σdBsm10log10(σm2)。5.2.2 RCS的频率特性与极化特性频率特性根据目标尺寸 L与波长 λ的相对关系散射机制可分为三个区域瑞利区L≪λRCS与频率的4次方成正比f4散射较弱。适用于低频对小型目标的探测分析。谐振区L∼λRCS随频率剧烈振荡是目标固有谐振模式的体现。复杂目标的RCS在此区域难以预测。光学区L≫λRCS趋于稳定与波长无关主要由目标的几何形状决定。大多数军用雷达和目标处于此区域。极化特性散射矩阵 S描述了入射波与散射波极化状态之间的关系[EhsEvs]Re−jkR[ShhSvhShvSvv][EhiEvi]其中下标 h,v分别代表水平极化和垂直极化。对于互易介质有 ShvSvh。极化信息是目标识别的重要特征。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)则 vx2y2服从均值为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 部件分解法CSCComponent Scattering Center该方法将复杂目标分解为多个简单的几何形状如平板、圆柱、球、角反射器等每个部件用一个“散射中心”来表征。原理每个散射中心 i具有以下属性位置 (xi,yi,zi)在目标坐标系中。频率依赖关系 σi∝fαi其中 αi是频率指数如平板 α2圆柱 α1尖锥 α0。姿态角依赖关系通常用一个方向图函数 Fi(ϕ,θ)来描述如平板的sinc函数圆柱的余弦函数。极化散射矩阵 Si。总RCS计算目标的总散射场是各散射中心散射场的矢量和考虑相位干涉Etotalsi1∑NEis⋅exp(−j2kk^⋅ri)其中 k^是雷达视线方向的单位矢量ri是第 i个散射中心的位置矢量k2π/λ。总RCS为 σtotallimR→∞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会产生多普勒频移 fd2vr/λ。在仿真中这体现为回波信号的载频偏移或脉冲间的相位变化。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处的目标双程大气衰减因子为 Latmexp(−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和高频仿真等多种建模方法的原理与适用场景。此外我们还讨论了目标微动产生的微多普勒效应以及大气传播对雷达波的影响。最后提出了一个集成化的目标与环境仿真模块设计方案。
雷达仿真（5）：雷达目标特性与电磁散射建模
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表示σdBsm10log10(σm2)。5.2.2 RCS的频率特性与极化特性频率特性根据目标尺寸 L与波长 λ的相对关系散射机制可分为三个区域瑞利区L≪λRCS与频率的4次方成正比f4散射较弱。适用于低频对小型目标的探测分析。谐振区L∼λRCS随频率剧烈振荡是目标固有谐振模式的体现。复杂目标的RCS在此区域难以预测。光学区L≫λRCS趋于稳定与波长无关主要由目标的几何形状决定。大多数军用雷达和目标处于此区域。极化特性散射矩阵 S描述了入射波与散射波极化状态之间的关系[EhsEvs]Re−jkR[ShhSvhShvSvv][EhiEvi]其中下标 h,v分别代表水平极化和垂直极化。对于互易介质有 ShvSvh。极化信息是目标识别的重要特征。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)则 vx2y2服从均值为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 部件分解法CSCComponent Scattering Center该方法将复杂目标分解为多个简单的几何形状如平板、圆柱、球、角反射器等每个部件用一个“散射中心”来表征。原理每个散射中心 i具有以下属性位置 (xi,yi,zi)在目标坐标系中。频率依赖关系 σi∝fαi其中 αi是频率指数如平板 α2圆柱 α1尖锥 α0。姿态角依赖关系通常用一个方向图函数 Fi(ϕ,θ)来描述如平板的sinc函数圆柱的余弦函数。极化散射矩阵 Si。总RCS计算目标的总散射场是各散射中心散射场的矢量和考虑相位干涉Etotalsi1∑NEis⋅exp(−j2kk^⋅ri)其中 k^是雷达视线方向的单位矢量ri是第 i个散射中心的位置矢量k2π/λ。总RCS为 σtotallimR→∞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会产生多普勒频移 fd2vr/λ。在仿真中这体现为回波信号的载频偏移或脉冲间的相位变化。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处的目标双程大气衰减因子为 Latmexp(−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和高频仿真等多种建模方法的原理与适用场景。此外我们还讨论了目标微动产生的微多普勒效应以及大气传播对雷达波的影响。最后提出了一个集成化的目标与环境仿真模块设计方案。
