留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

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

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

朱元江, 杨康, 鲍秋香, 徐宇, 蒋卫祥, 宋世千. 数字阵列宽带多波束发射设计与验证[J]. 电子科技大学学报, 2023, 52(6): 818-825. doi: 10.12178/1001-0548.2023077
引用本文: 朱元江, 杨康, 鲍秋香, 徐宇, 蒋卫祥, 宋世千. 数字阵列宽带多波束发射设计与验证[J]. 电子科技大学学报, 2023, 52(6): 818-825. doi: 10.12178/1001-0548.2023077
ZHU Yuanjiang, YANG Kang, BAO Qiuxiang, XU Yu, JIANG Weixiang, SONG Shiqian. Testbed for Digital Arrays Assisted Broadband Multi-Beam Transmission[J]. Journal of University of Electronic Science and Technology of China, 2023, 52(6): 818-825. doi: 10.12178/1001-0548.2023077
Citation: ZHU Yuanjiang, YANG Kang, BAO Qiuxiang, XU Yu, JIANG Weixiang, SONG Shiqian. Testbed for Digital Arrays Assisted Broadband Multi-Beam Transmission[J]. Journal of University of Electronic Science and Technology of China, 2023, 52(6): 818-825. doi: 10.12178/1001-0548.2023077

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

doi: 10.12178/1001-0548.2023077
详细信息
    作者简介:

    朱元江(1985 − ),男,博士生,研究员,主要从事数字发射多波束技术方面的研究

    通讯作者: 蒋卫祥,E-mail:wxjiang81@seu.edu.cn
  • 中图分类号: TN947

Testbed for Digital Arrays Assisted Broadband Multi-Beam Transmission

  • 摘要: 现有相控阵干扰设备通常采用分孔径或者分时隙方式实现多目标干扰,但前者无法利用全孔径增益,后者在目标数较多时存在干扰资源时域分配紧张的问题。针对这些问题,提出整数倍和分数倍延时分阶可调的方法,整数倍延时通过采样间隔延时,分数倍延时采用Farrow滤波器,实现了大带宽数字发射的目的。此外还对关键延时特性进行了仿真,并构建了小规模数字阵列样机。实验结果表明,该样机瞬时带宽为400 MHz,能够同时发射4波束,干扰方向和功率任意可调,延时特性与仿真结果一致。
  • 图  1  均匀线阵结构模型

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

    图  3  Farrow滤波器结构示意图

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

    图  5  Farrow滤波器性能仿真图

    图  6  Farrow滤波器信号延时图

    图  7  验证系统框图

    图  8  实验验证平台

    图  9  阵列法向波束方向图

    图  10  阵列扫描方向图

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

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

  • [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.
  • [1] 张先超, 任天时, 赵耀, 樊锐.  移动边缘计算时延与能耗联合优化方法 . 电子科技大学学报, 2022, 51(5): 737-742. doi: 10.12178/1001-0548.2021244
    [2] 钱璐, 邹林, 汪学刚.  宽带数字阵列雷达接收通道的优化设计方法研究 . 电子科技大学学报, 2022, 51(1): 52-56. doi: 10.12178/1001-0548.2021180
    [3] 李志鹏, 李兴和, 黄虎, 饶申宇.  无数字时延滤波器的宽带大规模阵列雷达去斜算法 . 电子科技大学学报, 2022, 51(3): 371-376. doi: 10.12178/1001-0548.2021395
    [4] 赵东成, 谢景昭, 孙罡.  面向时延敏感业务的服务功能链部署 . 电子科技大学学报, 2021, 50(6): 852-860. doi: 10.12178/1001-0548.2021178
    [5] 尹建勇, 汪云超, 杨鹏.  一种双圆极化宽带天线及其阵列应用 . 电子科技大学学报, 2020, 49(5): 695-699. doi: 10.12178/1001-0548.2019245
    [6] 王小菊, 敦涛, 马祥云, 徐如祥, 祁康成, 曹贵川, 林祖伦.  基于ANSYS的场发射阵列阴极焊接应力研究 . 电子科技大学学报, 2018, 47(1): 157-160. doi: 10.3969/j.issn.1001-0548.2018.01.024
    [7] 陈伟, 秦志光, 邓蔚, 聂鹏.  信道竞争感知的多速率多跳无线网络时延路由测量 . 电子科技大学学报, 2012, 41(5): 668-670. doi: 10.3969/j.issn.1001-0548.2012.05.005
    [8] 邹林, 汪学刚.  宽带数字阵列雷达波束形成的优化实现方法 . 电子科技大学学报, 2012, 41(3): 368-372. doi: 10.3969/j.issn.1001-0548.2012.03.008
    [9] 陈庆利, 蒲亦非, 黄果, 周激流.  数字图像的0~1阶Riemann-Liouville分数阶微分增强模板 . 电子科技大学学报, 2011, 40(5): 772-776. doi: 10.3969/j.issn.1001-0548.2011.05.026
    [10] 王洪, 吕幼新, 汪学刚, 刘磊.  宽带数字接收机的高效FPGA设计 . 电子科技大学学报, 2008, 37(3): 364-365,369.
    [11] 覃岭, 邓小炜, 何子述, 段军棋.  宽带数字下变频的高效实现方法研究 . 电子科技大学学报, 2008, 37(5): 698-700,708.
    [12] 谢昭贤, 黄大贵, 明爱国, 李智军, 下条诚.  基于同时发射声纳环移动机器人导航系统 . 电子科技大学学报, 2007, 36(2): 309-311,321.
    [13] 孟庆宇, 刘本永, 姚宏达.  基于分数本征特征的手写数字识别 . 电子科技大学学报, 2006, 35(3): 289-291.
    [14] 龙宁, 张凤荔.  基于FFT的数字多波束测向算法研究 . 电子科技大学学报, 2005, 34(1): 16-18.
    [15] 甘荣兵, 王建国.  改进的对星载SAR的应答式欺骗干扰 . 电子科技大学学报, 2005, 34(5): 614-617.
    [16] 徐成节, 陈祝明, 黄巍.  宽带LFM信号去斜率数字脉压方法性能分析 . 电子科技大学学报, 2004, 33(5): 492-495.
    [17] 滕旭东, 袁晓, 赵元英, 魏永豪.  数字分数微分器系数的快速算法 . 电子科技大学学报, 2004, 33(4): 457-460.
    [18] 林云松, 黄勇, 肖先赐.  宽带数字正交技术及性能分析 . 电子科技大学学报, 1999, 28(1): 20-23.
    [19] 唐瑜, 严梅, 洪福明.  自适应阵列抗CDMA多址干扰的研究 . 电子科技大学学报, 1998, 27(1): 21-24.
    [20] 刘乃琦, 左志宏.  计算机病毒对抗技术CVCM . 电子科技大学学报, 1997, 26(3): 283-288.
  • 加载中
图(12)
计量
  • 文章访问数:  4629
  • HTML全文浏览量:  1256
  • PDF下载量:  58
  • 被引次数: 0
出版历程
  • 收稿日期:  2023-03-21
  • 修回日期:  2023-07-15
  • 录用日期:  2023-05-03
  • 网络出版日期:  2023-11-29
  • 刊出日期:  2023-11-25

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

doi: 10.12178/1001-0548.2023077
    作者简介:

    朱元江(1985 − ),男,博士生,研究员,主要从事数字发射多波束技术方面的研究

    通讯作者: 蒋卫祥,E-mail:wxjiang81@seu.edu.cn
  • 中图分类号: TN947

摘要: 现有相控阵干扰设备通常采用分孔径或者分时隙方式实现多目标干扰,但前者无法利用全孔径增益,后者在目标数较多时存在干扰资源时域分配紧张的问题。针对这些问题,提出整数倍和分数倍延时分阶可调的方法,整数倍延时通过采样间隔延时,分数倍延时采用Farrow滤波器,实现了大带宽数字发射的目的。此外还对关键延时特性进行了仿真,并构建了小规模数字阵列样机。实验结果表明,该样机瞬时带宽为400 MHz,能够同时发射4波束,干扰方向和功率任意可调,延时特性与仿真结果一致。

English Abstract

朱元江, 杨康, 鲍秋香, 徐宇, 蒋卫祥, 宋世千. 数字阵列宽带多波束发射设计与验证[J]. 电子科技大学学报, 2023, 52(6): 818-825. doi: 10.12178/1001-0548.2023077
引用本文: 朱元江, 杨康, 鲍秋香, 徐宇, 蒋卫祥, 宋世千. 数字阵列宽带多波束发射设计与验证[J]. 电子科技大学学报, 2023, 52(6): 818-825. doi: 10.12178/1001-0548.2023077
ZHU Yuanjiang, YANG Kang, BAO Qiuxiang, XU Yu, JIANG Weixiang, SONG Shiqian. Testbed for Digital Arrays Assisted Broadband Multi-Beam Transmission[J]. Journal of University of Electronic Science and Technology of China, 2023, 52(6): 818-825. doi: 10.12178/1001-0548.2023077
Citation: ZHU Yuanjiang, YANG Kang, BAO Qiuxiang, XU Yu, JIANG Weixiang, SONG Shiqian. Testbed for Digital Arrays Assisted Broadband Multi-Beam Transmission[J]. Journal of University of Electronic Science and Technology of China, 2023, 52(6): 818-825. doi: 10.12178/1001-0548.2023077
  • 现代战争是依赖高技术的战争,战场形态正发生着重要变化,主要体现为敌对双方在电磁频谱作战空间的争夺越来越激烈。失去了制电磁权,也必将失去制空和制海权。电子战作为夺取制电磁权的重要手段,面临着来自更先进作战对象、更复杂电磁环境带来的严峻挑战。为适应高强度的对抗环境,电子对抗系统必须具备多功能、多目标对抗、宽带工作等能力[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  均匀线阵结构模型

      假定馈入通道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滤波器等实现。

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

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

      $$ {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$表示滤波器长度。该结构实现了输出与滤波器系数的解耦和实时控制,既能满足系统动态调整的需要,又能简化系统架构设计,但其代价是带来了更多乘法器资源的消耗。

      图  3  Farrow滤波器结构示意图

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

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

      构建一个原型滤波器长度为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工作带宽内可确保各通道的幅度、延时一致性。

      图  5  Farrow滤波器性能仿真图

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

      图  6  Farrow滤波器信号延时图

    • 按照上述原理,搭建了一个实物验证平台进行测试验证。整个阵列发射系统由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  多波束增益控制下发射方向图

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

参考文献 (14)

目录

    /

    返回文章
    返回