数字阵列宽带多波束发射设计与验证

朱元江，杨 康，鲍秋香，徐 宇，蒋卫祥，宋世千

(1.东南大学信息科学与工程学院 南京 210000；2.中国船舶集团有限公司第七二三研究所 江苏 扬州 225001)

现代战争是依赖高技术的战争，战场形态正发生着重要变化，主要体现为敌对双方在电磁频谱作战空间的争夺越来越激烈。失去了制电磁权，也必将失去制空和制海权。电子战作为夺取制电磁权的重要手段，面临着来自更先进作战对象、更复杂电磁环境带来的严峻挑战。为适应高强度的对抗环境，电子对抗系统必须具备多功能、多目标对抗、宽带工作等能力[1]。

数字阵列技术[2-3]是数字技术与阵列技术相结合的产物，与传统的阵列系统相比，可以快速改变波束指向和形状，容易形成多个波束，很好地满足了现代信息战争对电子对抗设备多功能、多模式以及多目标宽带对抗的应用需求。目前关于数字阵列技术的研究主要集中在数字阵列接收技术。如文献[4]提出一种基于可变分数时延滤波器的实时数字时延波束形成方法，该方法采用多相结构和Farrow 结构相结合的方法实现实时抽取和可变分数时延。文献[5]提出针对有载波宽带雷达信号的接收波束形成实现结构，实现宽带数字阵列各阵元传输时延的精确补偿。文献[6]针对多级抽取结合可变分数时延滤波器的典型接收通道结构，提出了一种新的接收通道优化设计方法。目前关于数字阵列发射技术的研究不足，文献[7]研制开发了数字相控阵雷达发射同频多波束的测试系统，并对发射多波束工作状态下天线的基本特性进行了研究和测试。此外，文献[8-10]也进行了数字阵列发射多波束的研究，但这些研究都是围绕雷达应用而开展的，带宽较窄，无法满足电子对抗系统宽带对抗的要求。针对宽带电子对抗需求，文献[11]研究了频域和时域形成宽带发射波束的原理，进行了理论仿真分析，但并未集成样机和实测验证。

本文针对当前数字阵列发射多波束技术研究的不足，阐述了在数字域实现宽带同时发射多波束的机理，并针对数字延时中关键的分数延时滤波器进行了设计和仿真；在此基础上，构建了原型样机和测试系统，对数字宽带同时发射多波束技术的可行性和有效性进行了评估和验证。

1 宽带发射同时多波束原理

考虑均匀线阵构型下单波束发射情形（如图1所示），假设阵元数目为N，各阵元间的距离为d，干扰信号发射方向等波位线与法线方向的夹角为θ，输入到各通道的干扰激励信号为s(t)。根据相控阵波束形成原理，若要使发射波束指向特定的角度θ，对于窄带信号，就必须使“空间相位差”和“阵内相位差”相等；对于宽带信号，就必须使“空间延时差”和“阵内延时差”相等。记波束发射方向为θ，那么阵元n(n取值范围为[1,N])与阵元1 的的空间波程差为 ΔLn-1=(n-1)dsinθ，对应的传输延时差为 Δ τn-1=ΔLn-1/c， 即来自阵元n的信号在空间传输时比来自阵元1 的信号要滞后Δτn-1。在这种情况下，要保证发射信号等相位面对准角度θ，在阵列内部来自通道n的干扰激励信号就需要比通道1 的干扰信号超前Δ τn-1。

图1 均匀线阵结构模型

假定馈入通道1 的中频信号为s(t)，经上变频后形成发射信号：

式中，f0为本振频率；按照延时差相等原则，那么阵元n的发射信号为：

当发射信号为窄带信号，信号带宽B远小于f0， 复包络s(t)缓慢变化。因此，可以忽略信号包络在各阵元上的差异，即s(t+Δτn-1)≈s(t)，式(2)可以简化为：

比较式(1)和式(3)可以发现，阵元1 和n之间存在相位差：

通过式(1)可知，窄带情况下改变各阵元辐射信号的相位差值（一般通过功率放大组件中的移相器进行），可实现波束的定向辐射；而这正是当前窄带相控阵的基本原理：通过天线阵元后的移相器的改变，实现波束的电扫。在这种情况下，由于移相器的取值瞬时只能有一个状态，因此每一个瞬时波束只能指向固定的方向，而不能实现同时发射多波束。另外这种通过移相器实现波束扫描的阵列天线，由于存在“孔径渡越”，导致在宽带宽角扫描情况下会出现波束指向偏斜。

从式(4)可以看出，移相器的取值只能取为中心频点对应的相移值2 πf0Δτn-1；当信号具有大宽带时，信号频率f0+B/2处 所需移相值为2π(f0+B/2)×Δτn-1，进而导致“阵内相位差”与“空间相位差”不等，出现色散现象[12]。以f0为基准，窄带移相方式的波束指向为：

式中， θmain表示主波束方向。可以发现，当f逐渐增大，波束指向会逐渐偏移预定指向，即色散现象逐渐严重。

鉴于此，在发射信号为宽带的情况下，复包络s(t)在各阵元上的差异不可忽略。根据式(2)，对于n通道首先要对输入信号s(t)进行延时控制，超前时间值为 Δτn-1；再进行移相补偿，相位补偿值为ej2πf0Δτn-1。 理想情况下，超前时间 Δ τn-1可采用线性相位的全通滤波器H(ejω)=ejωΔτn-1实现，对应的时域时延滤波器响应函数为h(t)=δ(t+ΔτN-1)。

将发射阵列信号写成向量形式进行统一处理，表达式为：

其中Bs(θ)=[δ(t),δ(t+Δτ)ej2πf0Δτ,···,δ(t+(N-1)Δτ)×ej2πf0(N-1)Δτ]T为 信号s′(t)的导向矢量，与发射角度、阵元数目、中心频率等参数有关。

根据叠加原理，可以推出宽带同时发射数字多波束的表达式。设K个信号分别向 θ1,θ2,···,θK方向发射，那么发射信号为各个信号之和。记和信号为，具体表达为：

式中，延时值 Δτi可表示为：

式中，T表示采样周期；L表示整数倍采样周期延时，且L=ceil(ΔLn-1/T)； c eil(·)表 示向下取整；Δ表示分数倍采样周期延时，且Δ =ΔτN-1/T-L。

基于分数延时滤波器的宽带发射波束形成结构如图2 所示。其中数字移相的作用是补偿相移值，数字延时线表示整数倍采样周期延时，分数倍时延通过Farrow 滤波器等实现。

图2 基于时延的宽带同时发射多波束形成网络模型

宽带同时发射多波束技术的优点是各波束都能够获得全孔径增益。相控阵列等效辐射功率计算公式为：

式中，Pt为 等效辐射功率；N为 阵元数；Ge为单元天线增益。若整个孔径划分为4 个规模相同的子阵，那么子阵的单元数目为整个阵列的1/4，那么等效辐射功率为=Pe(N/4)(N/4)Ge，因此子阵发射功率理论上仅为原来全孔径的1 /16；如果基于全孔径产生4 个不同指向的干扰波束，每个干扰波束都会利用全阵面进行发射，理论上每个干扰波束的等效辐射功率将下降到原来的1/4。显然，后者等效辐射功率是前者的4 倍。推而广之，相对于M个子阵的情况，采用全孔径M个数字发射波束，功率会提高M倍。

2 数字延时设计

数字延时分为整数延时和分数延时两个模块。将干扰源的输出数据进行时钟级延时，即可实现整数延时。分数延时模块采用数字滤波器实现，用于产生分数个采样间隔的延时。早期分数延时滤波器的设计主要有窗函数法、拉格朗日插值法等，后来又发展出以多项式为基础的Farrow 结构分数延时滤波器[13]。

设分数延时滤波器的频率响应为：

式中，D表示分数延时参数；Cn(D)为滤波器系数； ω表示数字频率上分数延时的带宽。在Farrow滤波器中Cn(D)可 表示为关于D的 多项式：Cn(D)=，代入到式(10)中即可得到Farrow 滤波器的表达式：

式中，l=n+Nm;al=Cn,m表 示滤波器系数；ϕl(D,ω)=Dme-jωn。 记系统角频率范围为 [ω0,ω1]，分数时延的变化范围为 [D0,D1]，则与理想分数时延滤波器相比，式(11)设计的滤波器频率响应误差为：

Farrow 滤波器多项式系数Cn,m的设计准则为最小化设计延时与实际延时之差，即：

将式(14)置零，利用最小二乘法即可得到式(13)的最优解，即Farrow 滤波器的最优设计系数为：

由上述分析可知，Farrow 结构是由M组N阶FIR 滤波器构成（如图3 所示），拟合的组数M决定了Farrow 结构的拟合效果，而相频特性和幅频特性由滤波器长度N决定，同时也决定了分辨率带宽和时延精度。在对资源的需求上，Farrow 结构由M×N个 乘法器构成，其中M表示拟合组数，N表示滤波器长度。该结构实现了输出与滤波器系数的解耦和实时控制，既能满足系统动态调整的需要，又能简化系统架构设计，但其代价是带来了更多乘法器资源的消耗。

图3 Farrow 滤波器结构示意图

为减小算法对资源的开销，将输入信号两倍内插后通过分数延时滤波器，如图4 所示。此时分数延时滤波器的通带降为原来的一半，从而大大降低了滤波器的阶数，最后通过两倍抽取恢复信号的采样频率；采用多相滤波的方式进行等效变换，可将延时滤波与抽取的顺序互换，从而内插与抽取相抵消，得到基于多相滤波的Farrow 滤波器。该滤波器每条支路均工作在初始采样频率，第1 条为奇数系数支路，第2 条为偶数系数支路。为减小运算量，将插值滤波器设计为高阶的半带滤波器；通过HBF 增加十几阶为代价，大大减少了运算时间和运算量。

图4 基于多相滤波的Farrow 滤波器结构图

构建一个原型滤波器长度为20、归一化带宽为0.5 的Farrow 滤波器，通过MATLAB 仿真，该分数延时滤波器的幅频响应特性和相位响应特性如图5 和图6 所示。由图5 可以看出，所设计的延时滤波器归一化频率 ω ∈[0,0.65 π]时都有平坦的幅度响应，在 ω＞0.65 π时幅度响应曲线急剧下降。同样图6 显示在 ω ∈[0,0.65 π]时，分别延时0.05～0.5个采样周期，所获得的延时值均较为稳定；在ω＞ 0.65 π时，延时出现急剧变化。在采样率为2 GHz情况下，滤波器能够获得的带宽为650 MHz，远大于系统设计所需要的400 MHz 带宽。因此，在400 MHz 工作带宽内可确保各通道的幅度、延时一致性。

图5 Farrow 滤波器性能仿真图

图6 Farrow 滤波器信号延时图

采用MATLAB 软件对单路信号的数字延时进行仿真，算法包括1 个整数倍延时单元和1 个分数倍延时单元。仿真时采样率设置为2 GHz，输入测试信号为100 MHz 载频的正弦波信号。在对初始信号进行整数倍采样点延时的基础上，另外再设置Farrow 滤波器的延时值分别为0.2、0.3 和0.5倍的采样间隔延时。延时后各信号局部放大曲线如图6b 所示，在第33 个采样点时，不同的延时曲线已明显发生偏移，相对于0.2 倍采样间隔延时曲线，其他两路信号曲线相对延时约为0.1 和0.3 个采样间隔，说明分数延时滤波器能够产生精确延时。

3 原型样机及测试平台构建

按照上述原理，搭建了一个实物验证平台进行测试验证。整个阵列发射系统由4 部分组成：天线、微波模块、干扰源及波束形成板、电源板以及控制计算机，具体如图7 所示。整个测试系统预留在线校正接口，但实测时采用仪表事先校正并形成校正表存放于干扰源及波束形成板。

图7 验证系统框图

天线采用2×7 规模的开口波导天线阵列（如图8a 所示），取其中水平向6 个单元构成一维线阵进行测试，天线阵元水平方向间距为39 mm。微波模块主要为一个放大耦合组件（如图8b 所示），主通道为6 路，预留1 路用于在线校正，工作频率为2.7～3.5 GHz，增益为19～21 dB，各通道间相位一致性≤±4°，耦合度为20±1 dB。干扰源及波束合成板包括3 个双通道数模转换器（digital to analog converter, DAC）、1 个大规模的现场可编程逻辑门阵列（field programmable gate array, FPGA）芯片以及一些对外接口芯片；在FPGA 中主要实现4 路直接数字频率合成（digital frequency synthesizer,DDS）干扰源数据，以及数据的移相、延时、叠加等运算。整个数字处理插箱如图8c 所示。控制计算机作为上位机，完成参数设置、波束发控等操作。

图8 实验验证平台

按照图7 所示构建测试样机，并在暗室环境下进行远场发射方向图的测试。测试所用仪表为罗德施瓦茨FSW 频谱仪。记阵列法向为90°，在此角度下顺时针方向方位角依次增大到180°，逆时针方向方位角依次减小到0°。首先测试一个发射波束的形成能力，形成一个法向波束，2.7 GHz 和3.1 GHz 频点下发射方向图如图9 所示。

图9 阵列法向波束方向图

在发射阵列扫描到+20°（即方位向110°）的情况下，2.7、2.9、3.1 GHz 频点下的发射方向图如图10a 所示。在发射阵列扫描到-45°（即方位向45°）的情况下，各频点下的发射方向图如图10b所示。两种情况下，各频点下发射波束的等效辐射功率有差异，但波束指向较为精准，都准确指向了110°和45°。两种扫描角度下，阵面未进行任何加权。文献[14]中判断宽带阵列的标准之一就是波束指向偏差需要小于波束宽度的四分之一，而本文测试结果远远超出这一标准，实现了宽带信号的发射能力。

图10 阵列扫描方向图

如图11 所示，测试样机能够有效生成4 个载频、指向不同的干扰波束，发射信号形式为连续波，信号频点分别为2.7、2.8、2.9 和3 GHz，波束指向分为90°、60°、10°、80°。为进一步验证该样机辐射功率重分配能力，调整各路干扰源的幅度调制系数，图12a 所示为调整2.7 GHz 和3 GHz干扰源通道幅度系数形成的方向图，从图中可直观看出2.7 GHz 和3 GHz 波束辐射功率明显下降，由于能量守恒，2.8 GHz 和2.9 GHz 波束辐射功率上升。图12b 所示为调整2.9 GHz 干扰源通道幅度系数形成的方向图，从图中可看出2.9 GHz 波束辐射功率明显上升。

图11 多波束同时发射方向图

图12 多波束增益控制下发射方向图

4 结 束 语

本文针对数字阵列干扰机同时多目标对抗的需求，利用小规模低频段阵列构建了原理样机，验证了宽带数字同时发射多波束技术的可行性。在仿真验证关键延时特性的基础上进行了暗室测试，获得了瞬时带宽400 MHz 的单干扰波束天线方向图，也获得了同时4 干扰波束的天线方向图，验证了数字域实现不同指向、不同载频、不同功率的多个独立可控宽带干扰波束的技术途径。该技术可根据不同应用场景实时动态调整发射频率、波束指向和辐射功率，在雷达、通信和电子对抗领域有广阔的应用前景。