雷达仿真（3）：雷达天线与波束形成的建模与仿真
3.1 引言：天线——雷达系统的“眼睛”与“嘴巴”
在雷达系统中，天线扮演着无可替代的双重角色：它既是将发射机产生的高频电能转换为定向辐射电磁波的“嘴巴”，又是收集目标反射回波的微弱电磁能量的“眼睛”。天线的性能直接决定了雷达的探测距离、角度分辨率、抗干扰能力和空域覆盖范围。因此，在雷达仿真中，对天线及其波束形成网络进行高保真建模，是准确评估雷达整体性能的前提。
随着相控阵技术和数字波束形成（DBF）技术的成熟与发展，现代雷达天线已从传统的机械扫描抛物面天线，演进为能够实现快速电扫描、多波束、自适应调零等先进功能的复杂系统。这给仿真建模带来了新的挑战和机遇。本部分将系统阐述雷达天线的基本原理、关键参数、各类天线模型的构建方法，并重点深入相控阵天线与数字波束形成的仿真实现，为构建高保真雷达仿真系统奠定坚实的天线建模基础。
3.2 天线基础参数与方向图函数建模
天线模型的核心是其辐射方向图（Antenna Pattern），它描述了天线辐射能量（或接收灵敏度）在空间不同方向上的相对强度分布。
3.2.1 关键天线参数
方向图函数 F(θ,ϕ)：在球坐标系中，以天线相位中心为原点，(θ,ϕ)分别表示俯仰角和方位角。方向图函数通常归一化，最大值设为1（0 dB）。
主瓣（Main Lobe）与波束宽度（Beamwidth）：方向图中辐射最强的区域称为主瓣。主瓣的宽度通常用半功率波束宽度（HPBW，即3 dB波束宽度）θ3dB和 ϕ3dB来衡量，它决定了雷达的角度分辨率。
旁瓣（Sidelobes）与栅瓣（Grating Lobes）：主瓣之外的其他辐射瓣。旁瓣电平（SLL）是衡量天线抗干扰和低截获概率性能的关键指标。对于阵列天线，当阵元间距过大时，会在空间其他方向形成与主瓣幅度相当的栅瓣，这是必须避免的。
天线增益（Gain）G：在最大辐射方向上的辐射强度与各向同性天线辐射强度的比值。它与方向性系数 D和辐射效率 η有关：G=ηD。在雷达方程中，增益是核心参数之一。
有效孔径面积（Effective Aperture Area）Ae：表征天线接收电磁波的能力，与增益的关系为 G=λ24πAe。
3.2.2 解析方向图函数模型
在仿真中，我们常用解析函数来近似真实天线的方向图，以平衡计算效率和精度。
余弦模型：适用于许多实际天线的主瓣近似。
F(θ)=cosn(θ)(对于∣θ∣≤π/2)
指数 n决定了波束的宽窄。该模型简单，但无法描述旁瓣。
sinc函数模型：适用于均匀照射的矩形或圆形孔径天线的主瓣及近旁瓣。
矩形孔径（尺寸 Lx×Ly）：
F(θ,ϕ)=sinc(λπLxsinθcosϕ)⋅sinc(λπLysinθsinϕ)
其中 sinc(x)=sin(x)/x。
圆形孔径（直径 D）：
F(θ)=π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(θ)=n=0∑N−1wnexp(jnψ)
其中：
ψ=kdsinθ+β，是相邻阵元间的空间-相位差。
k=2π/λ是波数。
β是相邻阵元间人为引入的相位偏移（移相器设置）。
wn是第 n个阵元的复加权（幅度和相位），用于控制波束形状和旁瓣。
波束指向 θ0由 β决定：为了使所有阵元在 θ0方向同相叠加，需要满足 ψ=0，即：
β=−kdsinθ0
这就是相控阵波束扫描的基本公式。通过电子方式快速改变 β，即可实现波束的无惯性扫描。
3.4.2 平面阵列与二维扫描
对于 M×N个阵元组成的矩形平面阵，阵元位于 (xmn,ymn)。其阵列因子为：
AF(θ,ϕ)=m=0∑M−1n=0∑N−1wmnexp[jk(xmnu+ymnv)]
其中 u=sinθcosϕ, v=sinθsinϕ是方向余弦。为了将波束指向 (θ0,ϕ0)，需要给阵元 (m,n)设置相位补偿：
βmn=−k(xmnu0+ymnv0),u0=sinθ0cosϕ0,v0=sinθ0sinϕ0
3.4.3 单元方向图与总方向图
相控阵天线的总方向图 Ftotal(θ,ϕ)是阵列因子 AF(θ,ϕ)与单个阵元的方向图 Felement(θ,ϕ)（称为“单元方向图”或“有源单元方向图”）的乘积：
Ftotal(θ,ϕ)=Felement(θ,ϕ)⋅AF(θ,ϕ)
单元方向图考虑了阵元间的互耦效应，通常比孤立阵元方向图更宽。在仿真中，单元方向图可以通过电磁仿真软件（如HFSS、CST）计算得到，或以简化模型（如余弦模型）近似。
3.4.4 量化误差与栅瓣抑制
移相器量化误差：数字移相器的相位调整是离散的（如4位移相器有16种相位状态）。这会导致波束指向误差和旁瓣电平升高，需要在仿真中建模。
栅瓣抑制：为了避免出现栅瓣，阵元间距 d必须满足：
d<1+∣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)=m=1∑Mwm∗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的约束下，使阵列输出功率最小（即抑制干扰和噪声）。
wminwHRws.t.wHa(θ0)=1
其解为：
wopt=aH(θ0)R−1a(θ0)R−1a(θ0)
其中：
R=E[x(k)xH(k)]是接收数据的协方差矩阵，需要通过采样数据估计：R^=K1∑k=1Kx(k)xH(k)。
a(θ0)是期望信号方向的导向矢量（Steering Vector）。对于位于 (θ0,ϕ0)的远场窄带信号，第 m个阵元的导向矢量分量为：
am(θ0,ϕ0)=Felement,m(θ0,ϕ0)⋅exp[jk(xmu0+ymv0)]
其中 Felement,m是第 m个阵元在目标方向的单元方向图值。
3.5.4 DBF仿真实现框架
在雷达系统仿真中，DBF模块通常位于接收通道数字化之后，信号处理（如脉冲压缩、CFAR）之前。其仿真流程如下：
阵元信号生成：根据目标、杂波、干扰模型，计算每个阵元接收到的基带信号 xm(k)，其中包含了空间相位差。
协方差矩阵估计：在自适应波束形成中，需要估计干扰加噪声的协方差矩阵 R^。通常假设在仅包含干扰和噪声的“训练区间”内进行估计。
加权向量计算：
对于静态波束形成，w是预先计算好的。
对于自适应波束形成，根据当前估计的 R^和期望方向的导向矢量 a(θ0)，实时计算 wopt。
波束形成输出：计算 y(k)=wHx(k)，得到合成后的波束输出信号。
多波束形成：如果需要同时形成 L个指向不同方向的波束，只需准备 L个不同的加权向量 wl(l=1,...,L)，并行计算 yl(k)=wlHx(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）的算法原理与实现方法。最后，给出了一个集成化的天线与波束形成仿真模块设计方案。