-
现代战争是依赖高技术的战争,战场形态正发生着重要变化,主要体现为敌对双方在电磁频谱作战空间的争夺越来越激烈。失去了制电磁权,也必将失去制空和制海权。电子战作为夺取制电磁权的重要手段,面临着来自更先进作战对象、更复杂电磁环境带来的严峻挑战。为适应高强度的对抗环境,电子对抗系统必须具备多功能、多目标对抗、宽带工作等能力[1]。
数字阵列技术[2-3]是数字技术与阵列技术相结合的产物,与传统的阵列系统相比,可以快速改变波束指向和形状,容易形成多个波束,很好地满足了现代信息战争对电子对抗设备多功能、多模式以及多目标宽带对抗的应用需求。目前关于数字阵列技术的研究主要集中在数字阵列接收技术。如文献[4]提出一种基于可变分数时延滤波器的实时数字时延波束形成方法,该方法采用多相结构和Farrow结构相结合的方法实现实时抽取和可变分数时延。文献[5]提出针对有载波宽带雷达信号的接收波束形成实现结构,实现宽带数字阵列各阵元传输时延的精确补偿。文献[6]针对多级抽取结合可变分数时延滤波器的典型接收通道结构,提出了一种新的接收通道优化设计方法。目前关于数字阵列发射技术的研究不足,文献[7]研制开发了数字相控阵雷达发射同频多波束的测试系统,并对发射多波束工作状态下天线的基本特性进行了研究和测试。此外,文献[8-10]也进行了数字阵列发射多波束的研究,但这些研究都是围绕雷达应用而开展的,带宽较窄,无法满足电子对抗系统宽带对抗的要求。针对宽带电子对抗需求,文献[11]研究了频域和时域形成宽带发射波束的原理,进行了理论仿真分析,但并未集成样机和实测验证。
本文针对当前数字阵列发射多波束技术研究的不足,阐述了在数字域实现宽带同时发射多波束的机理,并针对数字延时中关键的分数延时滤波器进行了设计和仿真;在此基础上,构建了原型样机和测试系统,对数字宽带同时发射多波束技术的可行性和有效性进行了评估和验证。
-
考虑均匀线阵构型下单波束发射情形(如图1所示),假设阵元数目为
$N$ ,各阵元间的距离为$d$ ,干扰信号发射方向等波位线与法线方向的夹角为$\theta $ ,输入到各通道的干扰激励信号为$s\left( t \right)$ 。根据相控阵波束形成原理,若要使发射波束指向特定的角度$\theta $ ,对于窄带信号,就必须使“空间相位差”和“阵内相位差”相等;对于宽带信号,就必须使“空间延时差”和“阵内延时差”相等。记波束发射方向为$\theta $ ,那么阵元n(n取值范围为[1,N])与阵元1的的空间波程差为$\Delta {L_{n - 1}} = \left( {n - 1} \right) d \sin \theta $ ,对应的传输延时差为$\Delta {\tau _{n - 1}} = {{\Delta {L_{n - 1}}} \mathord{\left/ {\vphantom {{\Delta {L_{n - 1}}} c}} \right. } c}$ ,即来自阵元$n$ 的信号在空间传输时比来自阵元1的信号要滞后$\Delta {\tau _{n - 1}}$ 。在这种情况下,要保证发射信号等相位面对准角度$\theta $ ,在阵列内部来自通道$n$ 的干扰激励信号就需要比通道1的干扰信号超前$\Delta {\tau _{n - 1}}$ 。假定馈入通道1的中频信号为
$s\left( t \right)$ ,经上变频后形成发射信号:$$ s'\left( t \right) = s\left( t \right) {{\rm{e}}^{{\rm{j}}2{\text π}{f_0}t}} $$ (1) 式中,
${f_0}$ 为本振频率;按照延时差相等原则,那么阵元$ n $ 的发射信号为:$$ s'\left( {t + \Delta {\tau _{n - 1}}} \right) = s\left( {t + \Delta {\tau _{n - 1}}} \right) {{\rm{e}}^{{\rm{j}}2{\text π} {f_0}\left( {t + \Delta {\tau _{n - 1}}} \right)}} $$ (2) 当发射信号为窄带信号,信号带宽
$ B $ 远小于${f_0}$ ,复包络$s\left( t \right)$ 缓慢变化。因此,可以忽略信号包络在各阵元上的差异,即$ s\left( {t + \Delta {\tau _{n - 1}}} \right) \approx s\left( t \right) $ ,式(2)可以简化为:$$ s'\left( {t + \Delta {\tau _{n - 1}}} \right) \approx s\left( t \right) {{\rm{e}}^{{\rm{j}}2{\text π}{f_0}t}} {{\rm{e}}^{{\rm{j}}2{\text π} {f_0}\Delta {\tau _{n - 1}}}} = s'\left( t \right) {{\rm{e}}^{{\rm{j}}2{\text π} {f_0}\Delta {\tau _{n - 1}}}} $$ (3) 比较式(1)和式(3)可以发现,阵元1和n之间存在相位差:
$$ \Delta {\varphi _{n - 1}} = 2{\text π} {f_0}\Delta {\tau _{n - 1}} = 2{\text π} {{\Delta {L_{n - 1}}} \mathord{\left/ {\vphantom {{\Delta {L_{n - 1}}} {{\lambda _0}}}} \right. } {{\lambda _0}}} $$ (4) 通过式(1)可知,窄带情况下改变各阵元辐射信号的相位差值(一般通过功率放大组件中的移相器进行),可实现波束的定向辐射;而这正是当前窄带相控阵的基本原理:通过天线阵元后的移相器的改变,实现波束的电扫。在这种情况下,由于移相器的取值瞬时只能有一个状态,因此每一个瞬时波束只能指向固定的方向,而不能实现同时发射多波束。另外这种通过移相器实现波束扫描的阵列天线,由于存在“孔径渡越”,导致在宽带宽角扫描情况下会出现波束指向偏斜。
从式(4)可以看出,移相器的取值只能取为中心频点对应的相移值
$2{\text π}{f_0}\Delta {\tau _{n - 1}}$ ;当信号具有大宽带时,信号频率${f_0} + {B \mathord{\left/ {\vphantom {B 2}} \right. } 2}$ 处所需移相值为$2{\text π} \left( {{f_0} + {B \mathord{\left/ {\vphantom {B 2}} \right. } 2}} \right) \times \Delta {\tau _{n - 1}}$ ,进而导致“阵内相位差”与“空间相位差”不等,出现色散现象[12]。以${f_0}$ 为基准,窄带移相方式的波束指向为:$$ \theta = \arcsin \left( {\frac{{{f_0}}}{f}\sin {\theta _{{\text{main}}}}} \right) $$ (5) 式中,
$ {\theta _{{\text{main}}}} $ 表示主波束方向。可以发现,当$f$ 逐渐增大,波束指向会逐渐偏移预定指向,即色散现象逐渐严重。鉴于此,在发射信号为宽带的情况下,复包络
$s\left( t \right)$ 在各阵元上的差异不可忽略。根据式(2),对于n通道首先要对输入信号$s\left( t \right)$ 进行延时控制,超前时间值为$\Delta {\tau _{n - 1}}$ ;再进行移相补偿,相位补偿值为$ {{\rm{e}}^{{\rm{j}}2{\text π} {f_0}\Delta {\tau _{n - 1}}}} $ 。理想情况下,超前时间$\Delta {\tau _{n - 1}}$ 可采用线性相位的全通滤波器$ H({{\rm{e}}}^{{\rm{j}}\omega })={{\rm{e}}}^{{\rm{j}}\omega \Delta {\tau }_{n-1}} $ 实现,对应的时域时延滤波器响应函数为$h\left( t \right) = \delta \left( {t + \Delta {\tau _{N - 1}}} \right)$ 。将发射阵列信号写成向量形式进行统一处理,表达式为:
$$ \begin{split} &{\boldsymbol{x}}'\left( t \right) = \left[ {\begin{array}{*{20}{c}} {{{x'}_0}\left( t \right)} \\ {{{x'}_1}\left( t \right)} \\ \vdots \\ {{{x'}_{N - 1}}\left( t \right)} \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} {\delta \left( t \right)} \\ {\delta \left( {t + \Delta \tau } \right) {{\rm{e}}^{{\rm{j}}2{\text π} {f_0} \Delta \tau }}} \\ \vdots \\ {\delta \left( {t + \left( {N - 1} \right) \Delta \tau } \right) {{\rm{e}}^{{\rm{j}}2{\text π} {f_0} \left( {N - 1} \right) \Delta \tau }}} \end{array}} \right] \times \\ &\qquad\qquad\qquad s'\left( t \right) = {{\boldsymbol{B}}_s}\left( \theta \right) s\left( t \right){{\rm{e}}^{{\rm{j}}2{\text π}{f_0}t}} \\[-8pt] \end{split} $$ (6) 其中
${{\boldsymbol{B}}_s}\left( \theta \right) = \big[ {\delta \left( t \right)} ,\delta \left( {t + \Delta \tau } \right) {{\rm{e}}^{{\rm{j}}2{\text π} {f_0} \Delta \tau }},\cdots ,\big.$ $\big.\delta \left( t + \left( {N - 1} \right) \Delta \tau \right) \times {{{\rm{e}}^{{\rm{j}}2{\text π} {f_0} \left( {N - 1} \right) \Delta \tau }} \big]^{\rm{T}}}$ 为信号$s'\left( t \right)$ 的导向矢量,与发射角度、阵元数目、中心频率等参数有关。根据叠加原理,可以推出宽带同时发射数字多波束的表达式。设
$K$ 个信号分别向${\theta _1},{\theta _2},\cdots,{\theta _K}$ 方向发射,那么发射信号为各个信号之和。记和信号为${x'_{{\rm{sum}}}}\left( t \right)$ ,具体表达为:$$ \begin{gathered} {{{\boldsymbol{x}}'}_{{\rm{sum}}}}\left( t \right) = \displaystyle\sum\limits_{i = 1}^K {{{\boldsymbol{B}}_s}\left( {{\theta _i}} \right)} {{s'}_i}\left( t \right)=\end{gathered} $$ $$ \begin{array}{l} \left[ {\begin{array}{*{20}{c}} {\displaystyle\sum\limits_{i = 1}^K {\delta \left( t \right){s_i}\left( t \right)} }\\ {\displaystyle\sum\limits_{i = 1}^K {\delta \left( {t + \Delta {\tau _i}} \right) {s_i}\left( t \right){{\rm{e}}^{{\rm{j}}2{\text π} {f_0} \Delta \tau }}} }\\ \vdots \\ {\displaystyle\sum\limits_{i = 1}^K {\delta \left( {t + \left( {N - 1} \right) \Delta {\tau _i}} \right) {s_i}\left( t \right){{\rm{e}}^{{\rm{j}}2{\text π} {f_0} \left( {N - 1} \right) \Delta \tau }}} } \end{array}} \right]{{\rm{e}}^{{\rm{j}}2{\text π}{f_0}t}}\\[-8pt] \end{array} $$ (7) 式中,延时值
$\Delta {\tau _i}$ 可表示为:$$ \Delta {\tau _i} = \left( {L + \Delta } \right) T $$ (8) 式中,
$T$ 表示采样周期;$L$ 表示整数倍采样周期延时,且$L = {\rm{ceil}}\left( {{{\Delta {L_{n - 1}}} \mathord{\left/ {\vphantom {{\Delta {L_{n - 1}}} T}} \right. } T}} \right)$ ;$ {\rm{ceil}}(·) $ 表示向下取整;$\varDelta $ 表示分数倍采样周期延时,且$\varDelta = {{\Delta {\tau _{N - 1}}} \mathord{\left/ {\vphantom {{\Delta {\tau _{N - 1}}} T}} \right. } T} - L$ 。基于分数延时滤波器的宽带发射波束形成结构如图2所示。其中数字移相的作用是补偿相移值,数字延时线表示整数倍采样周期延时,分数倍时延通过Farrow滤波器等实现。
宽带同时发射多波束技术的优点是各波束都能够获得全孔径增益。相控阵列等效辐射功率计算公式为:
$$ {P_t} = {P_e} N N {G_e} $$ (9) 式中,
${P_t}$ 为等效辐射功率;$N$ 为阵元数;${G_e}$ 为单元天线增益。若整个孔径划分为4个规模相同的子阵,那么子阵的单元数目为整个阵列的1/4,那么等效辐射功率为$P_t^1 = {P_e} {(N \mathord{\left/ {\vphantom {N 4}} \right. } 4)} {(N \mathord{\left/ {\vphantom {N 4}} \right. } 4)} {G_e}$ ,因此子阵发射功率理论上仅为原来全孔径的${1 \mathord{\left/ {\vphantom {1 {16}}} \right. } {16}}$ ;如果基于全孔径产生4个不同指向的干扰波束,每个干扰波束都会利用全阵面进行发射,理论上每个干扰波束的等效辐射功率将下降到原来的1/4。显然,后者等效辐射功率是前者的4倍。推而广之,相对于$M$ 个子阵的情况,采用全孔径$M$ 个数字发射波束,功率会提高$M$ 倍。 -
数字延时分为整数延时和分数延时两个模块。将干扰源的输出数据进行时钟级延时,即可实现整数延时。分数延时模块采用数字滤波器实现,用于产生分数个采样间隔的延时。早期分数延时滤波器的设计主要有窗函数法、拉格朗日插值法等,后来又发展出以多项式为基础的Farrow结构分数延时滤波器[13]。
设分数延时滤波器的频率响应为:
$$ H( {{{\rm{e}}^{{\rm{j}}\omega }}} ) = \sum\limits_{n = 0}^{N - 1} {{C_n}\left( D \right){{\rm{e}}^{ - {\rm{j}}\omega n}}} $$ (10) 式中,
$D$ 表示分数延时参数;${C_n}\left( D \right)$ 为滤波器系数;$\omega $ 表示数字频率上分数延时的带宽。在Farrow滤波器中${C_n}\left( D \right)$ 可表示为关于$D$ 的多项式:${C_n}\left( D \right) = \displaystyle\sum\limits_{m = 0}^{M - 1} {{C_{n,m}}{D^m}} $ ,代入到式(10)中即可得到Farrow滤波器的表达式:$$ H( {{{\rm{e}}^{{\rm{j}}\omega }}} ) = \displaystyle\sum\limits_{n = 0}^{N - 1} {\displaystyle\sum\limits_{m = 0}^{M - 1} {{C_{n,m}}{D^m}{{\rm{e}}^{ - {\rm{j}}\omega n}}} } = \displaystyle\sum\limits_{l = 0}^{NM - 1} {{a_l}{\phi _l}\left( {D,\omega } \right)} $$ (11) 式中,
$l = n + Nm$ ;${a_l} = {C_{n,m}}$ 表示滤波器系数;${\phi _l}\left( {D,\omega } \right) = {D^m}{{\rm{e}}^{ - {\rm{j}}\omega n}}$ 。记系统角频率范围为$\left[ {{\omega _0},{\omega _1}} \right]$ ,分数时延的变化范围为$\left[ {{D_0},{D_1}} \right]$ ,则与理想分数时延滤波器相比,式(11)设计的滤波器频率响应误差为:$$ \chi = \displaystyle\sum\limits_{l = 0}^{NM - 1} {{a_l}{\phi _l}\left( {D,\omega } \right)} - {{\rm{e}}^{ - {\rm{j}}\omega D}} $$ (12) Farrow滤波器多项式系数
${C_{n,m}}$ 的设计准则为最小化设计延时与实际延时之差,即:$$ \begin{split} &\qquad\qquad\min \displaystyle\int\nolimits_{{\omega _0}}^{{\omega _1}} {\displaystyle\int\nolimits_{{D_0}}^{{D_1}} {{{\left| \chi \right|}^2}{\rm{d}}D{\rm{d}}\omega } } \Leftrightarrow \\ & \min \displaystyle\int\nolimits_{{\omega _0}}^{{\omega _1}} {\int\nolimits_{{D_0}}^{{D_1}} {{{\left| {\displaystyle\sum\limits_{l = 0}^{NM - 1} {{a_l}{\phi _l}\left( {D,\omega } \right)} - {{\rm{e}}^{ - {\rm{j}}\omega D}}} \right|}^2}{\rm{d}}D{\rm{d}}\omega } } \Leftrightarrow \\ &\qquad\qquad\quad\min {{\boldsymbol{a}}^{\rm{T}}}{\boldsymbol{\varOmega}} {\boldsymbol{a}} - 2{{\boldsymbol{b}}^{\rm{T}}}{\boldsymbol{a}} + c \\ \end{split} $$ (13) 式中,
${{\boldsymbol{a}}^{\rm{T}}} = {\left[ {{a_0}}\quad{{a_1}}\;\; \cdots \;\;{{a_{NM - 1}}} \right]^{\rm{T}}}$ ,${{\boldsymbol{\varOmega }}_{i,j}} = \displaystyle\int\nolimits_{{\omega _0}}^{{\omega _1}} \displaystyle\int\nolimits_{{D_0}}^{{D_1}} {\phi _l}\left( {D,\omega } \right) {{\overline {{\phi _l}\left( {D,\omega } \right)} {\rm{d}}D{\rm{d}}\omega } }$ ,${b_l} = \displaystyle\int\nolimits_{{\omega _0}}^{{\omega _1}} {\displaystyle\int\nolimits_{{D_0}}^{{D_1}} {{{\rm{Re}}} \left\{ {{{\rm{e}}^{ - {\rm{j}}\omega D}}\overline {{\phi _l}\left( {D,\omega } \right)} } \right\}{\rm{d}}D{\rm{d}}\omega } }$ ,$c =\left( {{\omega _1} - {\omega _0}} \right)\left( {{D_1} - {D_0}} \right)$ 。求取式(13)关于$a$ 的偏导数,可得:$$ \frac{{\partial \left( {{{\boldsymbol{a}}^{\rm{T}}}{\boldsymbol{\varOmega }}{\boldsymbol{a}} - 2{{\boldsymbol{b}}^{\rm{T}}}{\boldsymbol{a}} + c} \right)}}{{\partial {\boldsymbol{a}}}} = 2{\boldsymbol{\varOmega }}{\boldsymbol{a}} - 2{\boldsymbol{b}} $$ (14) 将式(14)置零,利用最小二乘法即可得到式(13)的最优解,即Farrow滤波器的最优设计系数为:
$$ {{\boldsymbol{a}}_{{\rm{LS}}}} = {{\boldsymbol{\varOmega}} ^{ - 1}}{\boldsymbol{b}} $$ (15) 由上述分析可知,Farrow结构是由
$M$ 组$N$ 阶FIR滤波器构成(如图3所示),拟合的组数$M$ 决定了Farrow结构的拟合效果,而相频特性和幅频特性由滤波器长度$N$ 决定,同时也决定了分辨率带宽和时延精度。在对资源的需求上,Farrow结构由$M \times N$ 个乘法器构成,其中$M$ 表示拟合组数,$N$ 表示滤波器长度。该结构实现了输出与滤波器系数的解耦和实时控制,既能满足系统动态调整的需要,又能简化系统架构设计,但其代价是带来了更多乘法器资源的消耗。为减小算法对资源的开销,将输入信号两倍内插后通过分数延时滤波器,如图4所示。此时分数延时滤波器的通带降为原来的一半,从而大大降低了滤波器的阶数,最后通过两倍抽取恢复信号的采样频率;采用多相滤波的方式进行等效变换,可将延时滤波与抽取的顺序互换,从而内插与抽取相抵消,得到基于多相滤波的Farrow滤波器。该滤波器每条支路均工作在初始采样频率,第1条为奇数系数支路,第2条为偶数系数支路。为减小运算量,将插值滤波器设计为高阶的半带滤波器;通过HBF增加十几阶为代价,大大减少了运算时间和运算量。
构建一个原型滤波器长度为20、归一化带宽为0.5的Farrow滤波器,通过MATLAB仿真,该分数延时滤波器的幅频响应特性和相位响应特性如图5和图6所示。由图5可以看出,所设计的延时滤波器归一化频率
$\omega \in \left[ {0,0.65\;{\text π} } \right]$ 时都有平坦的幅度响应,在$\omega $ >$0.65\;{\text π} $ 时幅度响应曲线急剧下降。同样图6显示在$\omega \in \left[ {0,0.65\;{\text π}} \right]$ 时,分别延时0.05~0.5个采样周期,所获得的延时值均较为稳定;在$\omega $ >$0.65\;{\text π} $ 时,延时出现急剧变化。在采样率为2 GHz情况下,滤波器能够获得的带宽为650 MHz,远大于系统设计所需要的400 MHz带宽。因此,在400 MHz工作带宽内可确保各通道的幅度、延时一致性。采用MATLAB软件对单路信号的数字延时进行仿真,算法包括1个整数倍延时单元和1个分数倍延时单元。仿真时采样率设置为2 GHz,输入测试信号为100 MHz载频的正弦波信号。在对初始信号进行整数倍采样点延时的基础上,另外再设置Farrow滤波器的延时值分别为0.2、0.3和0.5倍的采样间隔延时。延时后各信号局部放大曲线如图6b所示,在第33个采样点时,不同的延时曲线已明显发生偏移,相对于0.2倍采样间隔延时曲线,其他两路信号曲线相对延时约为0.1和0.3个采样间隔,说明分数延时滤波器能够产生精确延时。
Testbed for Digital Arrays Assisted Broadband Multi-Beam Transmission
-
摘要: 现有相控阵干扰设备通常采用分孔径或者分时隙方式实现多目标干扰,但前者无法利用全孔径增益,后者在目标数较多时存在干扰资源时域分配紧张的问题。针对这些问题,提出整数倍和分数倍延时分阶可调的方法,整数倍延时通过采样间隔延时,分数倍延时采用Farrow滤波器,实现了大带宽数字发射的目的。此外还对关键延时特性进行了仿真,并构建了小规模数字阵列样机。实验结果表明,该样机瞬时带宽为400 MHz,能够同时发射4波束,干扰方向和功率任意可调,延时特性与仿真结果一致。Abstract: The existing phased array equipment usually adopts sub-aperture or sub-timeslot to realize multi-target jamming, but the former cannot take advantage of full aperture gain, and the latter cannot allocate the time-domain resource smoothly when the number of targets is large. In order to solve this problem, this paper proposes a method to adjust the integer and fractional delays in order to realize the purpose of large bandwidth digital transmission. The integer delay is realized by sampling interval delay, and the fractional delay is realized by Farrow filter. The key delay characteristics are simulated, and the minimum scale digital array prototype is constructed. The test results show that the simultaneous four beams with the instantaneous bandwidth of 400 MHz are realized, the direction and power of jamming beam can be adjusted arbitrarily, and the delay characteristics are also consistent with the simulation results.
-
[1] 丁国如, 孙佳琛, 王海超, 等. 复杂电磁环境下频谱智能管控技术探讨[J]. 航空学报, 2021, 42(4): 200-212. DING G R, SUN J C, WANG H C, et al. Discussion on technologies for intelligent spectrum management and control under complex electromagnetic environments[J]. Acta Aeronautica ET Astronautica Sinica, 2021, 42(4): 200-212. [2] 吴曼青. 数字阵列雷达及其进展[J]. 中国电子科学研究院学报, 2006(1): 11-16. WU M Q. The development of digital array radar[J]. Journal of China Academy of Electronics and Information Technology, 2006(1): 11-16. [3] 朱庆明. 数字阵列雷达述评[J]. 雷达科学与技术, 2004, 2(3): 136-141. ZHU Q M. An overview of digital array radar development[J]. Radar Science and Technology, 2004, 2(3): 136-141. [4] 邹林, 汪学刚. 宽带数字阵列雷达波束形成的优化实现方法[J]. 电子科技大学学报, 2012, 41(3): 368-372. ZOU L, WANG X G. Optimized beamforming method of wideband digital array radar based on time delay[J]. Journal of University of Electronic Science and Technology of China, 2012, 41(3): 368-372. [5] 范占春, 李会勇, 何子述. 基于分数时延的宽带数字阵列波束形成[J]. 雷达科学与技术, 2008, 6(6): 450-453. FAN Z C, LI H Y, HE Z S. Wideband digital array beamforming based on fractional delay[J]. Radar Science and Technology, 2008, 6(6): 450-453. [6] 钱璐, 邹林, 汪学刚. 宽带数字阵列雷达接收通道的优化设计方法研究[J]. 电子科技大学学报, 2022, 51(1): 52-56. QIAN L, ZOU L, WANG X G. Research on optimal design of receiving channel in wideband digital array radar[J]. Journal of University of Electronic Science and Technology of China, 2022, 51(1): 52-56. [7] 张德平, 王超, 袁乃昌. 基于DDS阵列的发射数字波束形成系统设计[J]. 现代防御技术, 2011, 39(1): 125-128. ZHANG D P, WANG C, YUAN N C. Design of transmit digital beam-forming system based on DDS array[J]. Modern Defence Technology, 2011, 39(1): 125-128. [8] 吴鸿超, 万长宁, 熊慎伟, 等. 数字相控阵雷达发射多波束特性研究[J]. 微波学报, 2014, 30(1): 6-9. WU H C, WAN C N, XIONG S W, et al. Transmitting multi-beam characteristics research of digital phase array radar[J]. Journal of Microwaves, 2014, 30(1): 6-9. [9] 彭志清, 姜兴, 谢跃雷. 球面共形阵数字波束发射阵列的实现研究[J]. 电视技术, 2015, 39(3): 89-92. PENG Z Q, JIANG X, XIE Y L. Research and implementation of digital transmitting beamforming for spherical conformal array antenna[J]. Video Engineering, 2015, 39(3): 89-92. [10] 鲁加国. 数字阵列雷达发射多波束的若干关键技术[J]. 雷达科学与技术, 2021, 19(3): 233-240. LU J G. The development of digital array radar[J]. Journal of China Academy of Electronics and Information Technology, 2021, 19(3): 233-240. [11] 张鸿伟. 基于数字波束形成的同时多目标干扰发射技术研究[D]. 西安: 西安电子科技大学, 2018. ZHANG H W. Study on synchronous multi-target jamming technology based on digital beamforming[D]. Xi’an: Xidian University, 2018. [12] 韦文, 李宁, 汤俊, 等. 基于分数时延的宽带自适应波束形成[J]. 清华大学学报(自然科学版), 2011, 51(7): 988-992. WEI W, LI N, TANG J, et al. Broadband adaptive digital beamforming based on fractional delay[J]. Journal of Tsinghua University (Science and Technology), 2011, 51(7): 988-992. [13] 范建俊, 李强, 李广军. 分数倍延时数字滤波器设计[J]. 微电子学, 2011, 41(2): 215-218. FAN J J, LI Q, LI G J. Design of fragment delay interpolation filter[J]. Microelectronics, 2011, 41(2): 215-218. [14] 张光义, 赵玉洁. 相控阵雷达技术[M]. 北京: 电子工业出版社, 2006. ZHANG G Y, ZHAO Y J. Phased array radar technology[M]. Beijing: Publishing House of Electronics Industry, 2006.