3.1 引言天线——雷达系统的“眼睛”与“嘴巴”在雷达系统中天线扮演着无可替代的双重角色它既是将发射机产生的高频电能转换为定向辐射电磁波的“嘴巴”又是收集目标反射回波的微弱电磁能量的“眼睛”。天线的性能直接决定了雷达的探测距离、角度分辨率、抗干扰能力和空域覆盖范围。因此在雷达仿真中对天线及其波束形成网络进行高保真建模是准确评估雷达整体性能的前提。随着相控阵技术和数字波束形成DBF技术的成熟与发展现代雷达天线已从传统的机械扫描抛物面天线演进为能够实现快速电扫描、多波束、自适应调零等先进功能的复杂系统。这给仿真建模带来了新的挑战和机遇。本部分将系统阐述雷达天线的基本原理、关键参数、各类天线模型的构建方法并重点深入相控阵天线与数字波束形成的仿真实现为构建高保真雷达仿真系统奠定坚实的天线建模基础。3.2 天线基础参数与方向图函数建模天线模型的核心是其辐射方向图Antenna Pattern它描述了天线辐射能量或接收灵敏度在空间不同方向上的相对强度分布。3.2.1 关键天线参数方向图函数 F(θ,ϕ)在球坐标系中以天线相位中心为原点(θ,ϕ)分别表示俯仰角和方位角。方向图函数通常归一化最大值设为10 dB。主瓣Main Lobe与波束宽度Beamwidth方向图中辐射最强的区域称为主瓣。主瓣的宽度通常用半功率波束宽度HPBW即3 dB波束宽度θ3dB​和 ϕ3dB​来衡量它决定了雷达的角度分辨率。旁瓣Sidelobes与栅瓣Grating Lobes主瓣之外的其他辐射瓣。旁瓣电平SLL是衡量天线抗干扰和低截获概率性能的关键指标。对于阵列天线当阵元间距过大时会在空间其他方向形成与主瓣幅度相当的栅瓣这是必须避免的。天线增益GainG在最大辐射方向上的辐射强度与各向同性天线辐射强度的比值。它与方向性系数 D和辐射效率 η有关GηD。在雷达方程中增益是核心参数之一。有效孔径面积Effective Aperture AreaAe​表征天线接收电磁波的能力与增益的关系为 Gλ24πAe​​。3.2.2 解析方向图函数模型在仿真中我们常用解析函数来近似真实天线的方向图以平衡计算效率和精度。余弦模型适用于许多实际天线的主瓣近似。F(θ)cosn(θ)(对于∣θ∣≤π/2)指数 n决定了波束的宽窄。该模型简单但无法描述旁瓣。sinc函数模型适用于均匀照射的矩形或圆形孔径天线的主瓣及近旁瓣。矩形孔径尺寸 Lx​×Ly​F(θ,ϕ)sinc(λπLx​​sinθcosϕ)⋅sinc(λπLy​​sinθsinϕ)其中 sinc(x)sin(x)/x。圆形孔径直径 DF(θ)πDsinθ/λ2J1​(πDsinθ/λ)​其中 J1​是第一类一阶贝塞尔函数。其第一旁瓣电平约为 -17.6 dB。高斯模型一种光滑的、无旁瓣的近似常用于理论分析。F(θ)exp[−α(θ3dB​θ​)2]其中 α2ln2≈1.386使得 F(θ3dB​/2)1/2​。3.2.3 方向图插值与坐标变换在仿真中天线方向图通常以数据表查找表的形式存储包含在离散的 (θi​,ϕj​)角度网格上的增益值 G(θi​,ϕj​)。当需要计算任意方向 (θ,ϕ)的增益时需要进行二维插值如双线性插值。此外需要建立天线本体坐标系与全局坐标系如北东地坐标系之间的转换关系以计算目标相对于天线波束指向的方向余弦。3.3 常见雷达天线类型及其仿真建模3.3.1 机械扫描抛物面天线原理利用抛物面的几何特性将位于焦点的馈源辐射的球面波转换为平面波定向辐射。仿真建模要点方向图可采用圆形孔径sinc函数或实际测量的方向图数据。扫描建模通过随时间改变天线本体坐标系相对于全局坐标系的欧拉角来模拟方位和俯仰的机械旋转。扫描速度、加速度、扇区扫描模式需要精确模拟。波束形状损耗由于目标并不总是处于波束中心需要考虑波束形状对回波信号的幅度调制。信号幅度需乘以归一化方向图函数值 F(θoff​,ϕoff​)其中 (θoff​,ϕoff​)是目标偏离波束中心的角度。3.3.2 平板缝隙阵列天线原理在波导或微带板上开凿一系列辐射缝隙通过控制缝隙的激励位置、倾角、长度来形成所需的方向图。仿真建模要点通常作为相控阵天线的一种实现形式来建模其方向图由阵列因子和单元方向图的乘积决定。3.3.3 相控阵天线这是现代雷达的主流也是仿真建模的重点和难点我们将在下一节详细展开。3.4 相控阵天线原理与建模相控阵天线由大量按规则排列如平面阵、共形阵的辐射单元阵元组成通过控制每个阵元发射/接收信号的相位或时间延迟实现波束在空间的快速电扫描而无需机械转动。3.4.1 阵列因子Array Factor对于由 N个相同阵元组成的线性阵列假设阵元间距为 d其阵列因子 AF(θ)描述了阵列的干涉方向图忽略单元方向图AF(θ)n0∑N−1​wn​exp(jnψ)其中ψkdsinθβ是相邻阵元间的空间-相位差。k2π/λ是波数。β是相邻阵元间人为引入的相位偏移移相器设置。wn​是第 n个阵元的复加权幅度和相位用于控制波束形状和旁瓣。波束指向 θ0​由 β决定为了使所有阵元在 θ0​方向同相叠加需要满足 ψ0即β−kdsinθ0​这就是相控阵波束扫描的基本公式。通过电子方式快速改变 β即可实现波束的无惯性扫描。3.4.2 平面阵列与二维扫描对于 M×N个阵元组成的矩形平面阵阵元位于 (xmn​,ymn​)。其阵列因子为AF(θ,ϕ)m0∑M−1​n0∑N−1​wmn​exp[jk(xmn​uymn​v)]其中 usinθcosϕ, vsinθsinϕ是方向余弦。为了将波束指向 (θ0​,ϕ0​)需要给阵元 (m,n)设置相位补偿βmn​−k(xmn​u0​ymn​v0​),u0​sinθ0​cosϕ0​,v0​sinθ0​sinϕ0​3.4.3 单元方向图与总方向图相控阵天线的总方向图 Ftotal​(θ,ϕ)是阵列因子 AF(θ,ϕ)与单个阵元的方向图 Felement​(θ,ϕ)称为“单元方向图”或“有源单元方向图”的乘积Ftotal​(θ,ϕ)Felement​(θ,ϕ)⋅AF(θ,ϕ)单元方向图考虑了阵元间的互耦效应通常比孤立阵元方向图更宽。在仿真中单元方向图可以通过电磁仿真软件如HFSS、CST计算得到或以简化模型如余弦模型近似。3.4.4 量化误差与栅瓣抑制移相器量化误差数字移相器的相位调整是离散的如4位移相器有16种相位状态。这会导致波束指向误差和旁瓣电平升高需要在仿真中建模。栅瓣抑制为了避免出现栅瓣阵元间距 d必须满足d1∣sinθmax​∣λ​其中 θmax​是最大扫描角。在宽角扫描时需要更小的阵元间距。图3-1相控阵天线波束形成与扫描原理框图┌─────────────────────────────────────────────────────────────┐ │ 相控阵天线仿真模型 │ ├─────────────────────────────────────────────────────────────┤ │ ┌─────────────┐ ┌─────────────┐ ┌─────────────┐ │ │ │ 阵元1 │ │ 阵元2 │ │ 阵元N │ │ │ │ 坐标(x1,y1) │ │ 坐标(x2,y2) │ │ 坐标(xN,yN) │ │ │ │ 方向图Fe1 │ │ 方向图Fe2 │ │ 方向图FeN │ │ │ └──────┬──────┘ └──────┬──────┘ └──────┬──────┘ │ │ │ │ │ │ │ ┌──────▼──────┐ ┌──────▼──────┐ ┌──────▼──────┐ │ │ │ 复加权w1 │ │ 复加权w2 │ │ 复加权wN │ │ │ │ (幅度a1, │ │ (幅度a2, │ │ (幅度aN, │ │ │ │ 相位φ1) │ │ 相位φ2) │ │ 相位φN) │ │ │ └──────┬──────┘ └──────┬──────┘ └──────┬──────┘ │ │ │ │ │ │ │ └──────────────────┼──────────────────┘ │ │ ▼ │ │ ┌────────────────────┐ │ │ │ 信号合成器 │ │ │ │ (相干叠加) │ │ │ └─────────┬──────────┘ │ │ ▼ │ │ ┌────────────────────┐ │ │ │ 总输出信号y(t) │ │ │ └────────────────────┘ │ └───────────────────────────┬─────────────────────────────────┘ │ ┌───────────▼───────────┐ │ 波束控制算法 │ ├───────────────────────┤ │ 输入期望波束指向 │ │ (θ0, φ0) │ │ 输出各阵元加权值 │ │ wn an·exp(jφn)│ │ 其中φn -k·(xn u0 │ │ yn v0) Δφn│ │ (Δφn: 用于低旁瓣的附加相位)│ └───────────────────────┘3.5 数字波束形成DBF仿真数字波束形成是相控阵技术的进一步发展它在每个阵元通道后都进行数字化采样然后在数字域通过加权求和来形成波束。DBF提供了极大的灵活性可以同时形成多个独立波束并实现自适应抗干扰。3.5.1 DBF系统模型一个典型的接收DBF系统模型如下设有 M个阵元第 m个阵元在 k时刻接收到的复基带信号为 xm​(k)。则DBF的输出 y(k)为y(k)wHx(k)m1∑M​wm∗​xm​(k)其中x(k)[x1​(k),x2​(k),...,xM​(k)]T是 M×1的快拍数据向量。w[w1​,w2​,...,wM​]T是 M×1的复加权向量。H表示共轭转置∗表示复共轭。3.5.2 静态加权低旁瓣波束形成静态加权 w的设计目标是在期望方向形成主瓣同时抑制旁瓣。常用的加权函数有均匀加权wm​1主瓣最窄但第一旁瓣电平高-13.2 dB。切比雪夫Chebyshev加权在给定主瓣宽度下能获得等旁瓣电平是优化旁瓣的经典方法。泰勒Taylor加权近似切比雪夫加权但更易于实现旁瓣包络呈衰减趋势。余弦窗加权如汉宁窗、汉明窗旁瓣抑制效果好但主瓣会展宽。在仿真中可以根据所需的波束宽度和旁瓣电平要求计算相应的加权向量 w。3.5.3 自适应波束形成ADBF自适应波束形成能够根据外部干扰环境实时调整加权向量 w在期望信号方向形成主瓣的同时在干扰来向形成零陷。最著名的算法是最小方差无失真响应MVDR或线性约束最小方差LCMV波束形成器。MVDR准则在保证期望方向增益为1的约束下使阵列输出功率最小即抑制干扰和噪声。wmin​wHRws.t.wHa(θ0​)1其解为wopt​aH(θ0​)R−1a(θ0​)R−1a(θ0​)​其中RE[x(k)xH(k)]是接收数据的协方差矩阵需要通过采样数据估计R^K1​∑k1K​x(k)xH(k)。a(θ0​)是期望信号方向的导向矢量Steering Vector。对于位于 (θ0​,ϕ0​)的远场窄带信号第 m个阵元的导向矢量分量为am​(θ0​,ϕ0​)Felement,m​(θ0​,ϕ0​)⋅exp[jk(xm​u0​ym​v0​)]其中 Felement,m​是第 m个阵元在目标方向的单元方向图值。3.5.4 DBF仿真实现框架在雷达系统仿真中DBF模块通常位于接收通道数字化之后信号处理如脉冲压缩、CFAR之前。其仿真流程如下阵元信号生成根据目标、杂波、干扰模型计算每个阵元接收到的基带信号 xm​(k)其中包含了空间相位差。协方差矩阵估计在自适应波束形成中需要估计干扰加噪声的协方差矩阵 R^。通常假设在仅包含干扰和噪声的“训练区间”内进行估计。加权向量计算对于静态波束形成w是预先计算好的。对于自适应波束形成根据当前估计的 R^和期望方向的导向矢量 a(θ0​)实时计算 wopt​。波束形成输出计算 y(k)wHx(k)得到合成后的波束输出信号。多波束形成如果需要同时形成 L个指向不同方向的波束只需准备 L个不同的加权向量 wl​(l1,...,L)并行计算 yl​(k)wlH​x(k)即可。3.6 天线与波束形成仿真模块设计在雷达仿真系统中天线与波束形成模块应设计为高度可配置和可复用的组件。图3-2天线与波束形成仿真模块结构框图┌─────────────────────────────────────────────────────────────┐ │ 天线与波束形成仿真模块接口定义 │ ├─────────────────────────────────────────────────────────────┤ │ 输入接口 │ │ 1. 目标/干扰源信息方位角θ, 俯仰角φ, 距离R, 速度v, RCS σ │ │ 2. 发射信号参数波形s(t), 载频f0, 脉冲重复间隔PRI │ │ 3. 控制命令波束指向(θ_beam, φ_beam), 扫描模式, 加权方式 │ │ 4. 环境参数传播损耗L_prop │ ├─────────────────────────────────────────────────────────────┤ │ 内部核心模型 │ │ │ │ ┌────────────────────────────────────────────────────┐ │ │ │ 天线本体模型 │ │ │ │ ┌──────────────┐ ┌──────────────┐ │ │ │ │ │ 阵元位置布局 │ │ 单元方向图 │ │ │ │ │ │ (几何结构) │ │ (查找表/函数)│ │ │ │ │ └──────┬───────┘ └──────┬───────┘ │ │ │ │ │ │ │ │ │ │ └─────────┬───────┘ │ │ │ │ ▼ │ │ │ │ ┌──────────────────┐ │ │ │ │ │ 方向图计算引擎 │ │ │ │ │ │ F_total(θ,φ) │ │ │ │ │ │ Fe(θ,φ) * AF(θ,φ)│ │ │ │ │ └─────────┬────────┘ │ │ │ └───────────────────┼──────────────────────────────┘ │ │ │ │ │ ┌───────────────────▼──────────────────────────────┐ │ │ │ 波束形成与扫描控制模块 │ │ │ │ ┌──────────┐ ┌──────────┐ ┌──────────┐ │ │ │ │ │ 波束控制 │ │ 加权计算 │ │ 扫描调度 │ │ │ │ │ │ 算法 │ │ (静态/ │ │ (机械/ │ │ │ │ │ │ (LCMV等) │ │ 自适应) │ │ 电扫) │ │ │ │ │ └─────┬────┘ └────┬─────┘ └────┬─────┘ │ │ │ │ │ │ │ │ │ │ │ └─────┬──────┴──────┬──────┘ │ │ │ │ ▼ ▼ │ │ │ │ ┌─────────────────┐┌─────────────────┐ │ │ │ │ │ 发射波束赋形 ││ 接收波束形成 │ │ │ │ │ │ 加权向量W_tx ││ 加权向量W_rx │ │ │ │ │ └─────────────────┘└─────────────────┘ │ │ │ └───────────────────┬───────────────────────────┘ │ │ │ │ │ ┌───────────────────▼──────────────────────────────┐ │ │ │ 信号层面仿真 │ │ │ │ 对于每个阵元m (或每个发射通道): │ │ │ │ 1. 计算目标/干扰相对于阵元m的时延τ_m和多普勒f_dm│ │ │ │ 2. 计算阵元m的单元方向图增益 G_m(θ,φ) │ │ │ │ 3. 生成阵元m接收信号: x_m(t) A_m * s(t-τ_m) * │ │ │ │ exp(j2πf_dm t) * G_m(θ,φ) noise_m(t) │ │ │ │ 4. (接收时)应用接收加权w_rx,m合成波束输出y(t) │ │ │ │ 5. (发射时)应用发射加权w_tx,m生成各阵元发射信号│ │ │ └──────────────────────────────────────────────────┘ │ ├─────────────────────────────────────────────────────────────┤ │ 输出接口 │ │ 1. 发射模式各阵元的激励信号用于空间辐射计算 │ │ 2. 接收模式合成后的波束输出信号y(t)送往后级处理 │ │ 3. 实时波束指向与方向图数据用于态势显示 │ └─────────────────────────────────────────────────────────────┘关键实现考虑计算效率对于大规模阵列成千上万个阵元逐个阵元计算信号计算量巨大。需采用并行计算如GPU加速和优化算法如利用FFT计算阵列因子。互耦效应在高精度仿真中需要考虑阵元间的电磁互耦它会改变单元方向图和阵列的有效激励分布。可以通过引入互阻抗矩阵或使用有源单元方向图数据来建模。宽带信号处理对于宽带信号简单的相位偏移会导致波束色散波束指向随频率变化。此时需要使用真时延线TDL代替移相器在仿真中需建模频变导向矢量。校准误差建模实际阵列存在通道幅度/相位不一致、阵元位置误差等需要在加权向量中引入相应的误差模型。3.7 小结本部分深入探讨了雷达天线与波束形成的建模与仿真。我们从天线的基本参数和方向图函数入手建立了描述天线空间辐射特性的数学模型。重点剖析了相控阵天线的核心——阵列因子理论以及实现波束电扫描的原理。进一步我们深入到数字波束形成领域阐述了静态加权和自适应加权MVDR的算法原理与实现方法。最后给出了一个集成化的天线与波束形成仿真模块设计方案。