一种基于变参考电平的零点凹口双重可控方向图综合改进方法.pdf

一种基于变参考电平的零点凹口双重可控方向图综合改进方法刘海波1,2，韩宇新1,2，武兴恺1,2，陈广茂1,2（1.北京理工大学信息与电子学院,北京100081；2.北京理工大学 信息与电子学院 国家重点实验室电磁感知研究中心,北京100081）摘 要：为了解决基于变参考电平的自适应波束形成方法无法同时控制方向图零点和凹口的问题，研究了基于最大信噪比准则的任意阵列方向图综合方法，重点对主瓣宽度固定条件下能够实现最低等副瓣电平的方向图综合方法进行了理论分析，并对前人提出的可变参考电平参量进行了重新定义，新定义的变参考电平能够在波束主瓣宽度固定的前提下，同时产生等副瓣电平、可控凹口宽度和深度和可控零点深度的抗干扰波束，并通过计算机仿真验证了算法的有效性.该方法对于提升相控阵雷达对抗杂波和干扰性能具有较高的应用价值.关键词：波束形成；共形阵；抗干扰中图分类号：TN958.92 文献标志码：A 文章编号：1001-0645(2024)04-0421-07DOI：10.15918/j.tbit1001-0645.2023.070An Improved Method for Synthesizing Dual-Control AdaptivePattern with Nulls and Gaps Based on Variable Reference VoltageLIU Haibo1,2，HAN Yuxin1,2，WU Xingkai1,2，CHEN Guangmao1,2（1.School of Information and Electronics,Beijing Institute of Technology,Beijing 100081,China;2.Electromagnetic Sensing Research Center of CEMEE State Key Laboratory,School ofInformation and Electronics,Beijing Institute of Technology,Beijing 100081,China）Abstract：To solve the problem that the pattern synthesis method based on the adaptive beam forming with vari-able reference voltage cannot simultaneously control the nulls and gaps of the pattern,a pattern synthesis meth-od based on the maximum signal-to-noise ratio(MSNR)for arbitrary arrays was proposed in this paper,focusingon the theoretical analysis of pattern synthesis methods to achieve the lowest side lobe level under a fixed mainlobe width.Redefining the variable reference voltage parameters proposed by previous studies,the new variablereference voltage was designed to generate an anti-interference beam with controllable gaps and nulls,as well asequal side lobe levels under the condition of a fixed beam main lobe width.The effectiveness of the algorithmwas verified by computer simulations.This method has provided with a significant practical value to improve theperformance of phased array radar in resistance of clutters and interferences.Key words：beamforming；conformal array；anti-interference 近几十年以来，相控阵雷达由于其先进的性能而被广泛应用 1 3，相控阵天线方向图综合是阵列信号处理的热点研究问题4 6.常见的方向图综合方法主要分为 3 类：基于凸优化理论的方向图综合算法、基于加权最小二乘的方向图综合算法和基于自适应阵理论的方向图综合算法，基于凸优化理论的方向图综合算法能得到的权值向量的最优解，但需要将方向图综合问题转化为凸问题；基于加权最小二乘 收稿日期：2023 03 29基金项目：国家自然科学基金资助项目(31727901)作者简介：刘海波（1980），男，博士，讲师，E-mail：.第 44 卷第 4 期北 京 理 工 大 学 学 报Vol.44No.42024 年 4 月Transactions of Beijing Institute of TechnologyApr.2024的方向图综合算法通过调整最小二乘代价函数实现，调整中需要对不同参数通过试凑法选取，无法保证阵列响应满足期望条件；基于自适应阵理论的方向图综合算法，依据最大信噪比准则迭代求解权值向量，能够实现任意阵型的方向图综合.1990 年，OLEN 等7首次系统地将自适应波束形成原理运用于阵列的方向图综合，该算法假设干扰位于方向图低副瓣区域所对应的空间方向，并基于最大信噪比准则实现了方向图的低副瓣，且算法收敛后形成的方 向 图 具 有 等 副 瓣 特 性.1999 年 在 Olen-Compton算法基础上，提出一种能够在实现等副瓣的同时控制波束的主瓣形状的波束方向图综合方法8.2003 年，刁跃龙9对该算法的迭代公式进行了修改，降低了迭代次数缩短了算法收敛所需时间.此外，还引入了变参考电平的概念，基于变参考电平该算法不仅能使方向图在固定主瓣宽度下达到最低的等副瓣电平，还能同时产生宽度和深度可控的凹口.2014 年将该算法应用于 MTD 滤波器的设计10，2018 年基于 Olen-Compton 算法解决了恒定响应宽带波束形成问题11.本文在 Olen-Compton 算法基础上，通过对可变参考电平进行修改，得到了一种新的波束方向图综合算法，新算法能够在共形阵方向图主瓣宽度一定的前提下，同时实现方向图的等低副瓣、可控凹口和可控零点，从而使波束能灵活地对抗杂波和干扰.与最新的等低副瓣方向图综合算法相比12 14，本文所提算法能够保持方向图凹口电平相等，且加入零点和凹口后，对其他区域的副瓣电平影响较小，最后通过计算机仿真验证了算法的有效性.1 算法模型 1.1 基于自适应波束形成的方向图综合方法简述wB()阵列天线的方向图综合问题就是在阵元位置、阵元数量和单元方向图已知的前提下，通过某种方法来找到一组复数权值，该权值能使方向图在期望波束指向处具有最高的增益，且旁瓣电平满足某些特定的要求.阵列天线的方向图函数可以表示为p()=?wHv()?（1）v()wv()式中：和分别代表阵列的导向矢量和复权值；H 代表 Hermite 转置.导向矢量可表示为v(i)=g1(i)ej1(i)g2(i)ej2(i)gN(i)ejN(i)T（2）iNgi()ii()ii()=kxikxii式中：表示来波方向也是期望波束指向；表示阵元个数；表示第 个阵元的单元方向图；T 表示矩阵转置；为电磁波在空间中传播时，由于波程差导致的第 个阵元的相位延迟，即，其中 表示波数，表示阵列中的第 个阵元在空间中的坐标位置.在阵列天线方向图的综合问题中，运用最大信噪比(MSNR)准则进行波束形成，就是寻找一组复权值，该权值能使阵列天线的方向图在波束指向的方向获得的增益尽可能大，而干扰方向阵列增益尽量小，从而使得阵列的输出满足信干噪比(SINR)最大9.阵列的输入可以表示为x=xs+xn+xf（3）xsxnxf式中：表示期望信号；表示天线单元本身包含的噪声；表示干扰信号.将非期望信号合并，有xu=xn+xf（4）xu式(4)中，表示无用信号.根据自适应阵原理，在 MSNR 准则下，最优权矢量表示为wopt=q（5）qR1uRsR1uRsRsRu式中：为一个大于 0 的常数；为的最大特征值所对应的特征向量，而和分别为无用信号和期望信号的自相关矩阵，其数学表达式为Ru=ExuxTu（6）Rs=ExsxTs（7）当期望信号只来自一个方向时，式(5)可表示为wsopt=R1vs（8）vs式中为所期望波束指向的导向矢量.sfAsAfsf12假设存在一个期望信号和一个干扰信号，它们同时入射到阵列上，其来波方向分别记为和，信号振幅记为和，初相记为和，角频率分别为和，则这两个信号可表示为xs=Asej(1t+s)vs（9）xf=Afej(2t+f)vf（10）vsvf式中和分别为期望信号和干扰信号的导向矢量.根据式(6)和(7)，有Rs=ExsxTs=A2svsvTs（11）Rf=ExfxTf=A2fvfvTf（12）xn设噪声矢量可表示为xn=n1(t)n2(t)ni(t)nN(t)（13）ni(t)ii=1,2,.,N式中为第 个阵元的噪声，.设噪声功422北 京 理 工 大 学 学 报第 44 卷2率均为且相互独立，则噪声的自相关矩阵为Rn=ExnxTn=2I（14）假设信号，噪声和干扰均满足两两独立分布7，则有Ru=Rf+Rn=2I+A2fvfvTf（15）s=A2ff=A2f令，则信噪比和干噪比可写为RSN=A2s2=s2（16）RIN=A2f2=f2（17）2=1RSN=s RIN=f为讨论简便，不妨设，则，.因此，期望信号和无用信号的自相关矩阵9可表示为Rs=ExsxTs=svsvTs（18）Rf=ExfxTf=fvfvTf（19）RsRfwsopt此时便可以利用和来确定 MSNR 准则下阵列的最优权矢量.Olen-Compton 算法的核心主要基于以下表达式wsopt=R1uvs（20）Ru=2I+Mi=1fivfivTfi=2I+VfFVTf（21）fiivfiVfF式(20)和(21)所描述的也是 MSNR 准则下的最优权适量求解问题，只不过将干扰信号由一个变为多个的同时，式(19)也就相应地变为式(21).根据式(17)，既能表示第 个干扰信号，也表示其干噪比.表示其导向矢量，和分别可表示为F=diagf1f2 fm（22）Vf=vf1vf2 vfm（23）D()d()现在简述 Olen-Compton 算法的内容.首先根据式(2)来确定期望信号的导向矢量.然后确定一个期望的副瓣参考电平，若要求其比主瓣峰值低dB，则期望副瓣电平(线性值)可换算为d()=Vpeak10D()/20（24）Vpeakp()p()式中，表示方向图的峰值电平，可由式(1)确定.m假设个干扰一次性地入射到方向图的副瓣区fi(i=1,2,.,m)fifik域，这些干扰的入射角分别记为，接着对赋初值，为了方便计算，可将其直接取 1.然后，在处对当前方向图的电平和期望方向图的电平进行做差比较，若差值大于 0，则表明此处干扰功率不足，需要增大当前位置的干扰功率，若差值小于 0，则需要减小当前位置的干扰功率.随着干扰功率在不断调整，阵列天线方向图的副瓣电平也在不断下降，因此主瓣宽度也在不断展宽.因此，在每次迭代时，务必要保证当前方向图的主瓣区域内干扰功率为 0.由于每次迭代后，方向图的峰值电平会发生变化，这里将式(24)进行修正，使其为迭代次数 和角度 的函数，即d(,k)=Vpeak(k)10D()/20（25）Vpeak(k)kp(,k)式中表示第 次迭代时方向图的最大峰值电平.Olen 算法的干扰功率迭代公式为9fi(k+1)=max0,fi(k)+K(fi,k),0,fi 0(fi,k)kfi式中：表示方向图主瓣区，该区域内干扰功率被置 0；表示第 次迭代时方向的干扰功率；为迭代系数，；表示第 次迭代后在方向上当前方向图的电平值与期望方向图的电平值之差，用公式表示为(fi,k)=p(fi,k)d(fi,k)（27）这里再次强调：由于干扰功率在随着迭代不断变化，导致基于 MSNR 准则的波束形成器的最优权值也变化.因此，方向图的主瓣区也在随之改变.因此，每次迭代后务必要重新确定主瓣区域.KKKOlen-Compton 算法存在的问题是：算法收敛速度太慢，一般需要上百次才能收敛，且值的选取不具有通用性，而值的大小又和算法的收敛性和稳定性密切相关，因此每次实验只能通过“累尝法”来寻找合适的9.刁跃龙9在对式(26)所述的迭代公式进行详细分析后，将迭代公式修改为fi(k+1)=max0,fi(k)(1+K(fii,k)/d(fii,k),fi(Left(k),Right(k)0,fi(Left(k),Right(k)（28）(fi,k)由于式(28)考虑了干扰间的相对大小并对电平差进行了归一化处理，因此算法的收敛速度K大幅加快，且值易于选取(一般为 1 或者 0.1).自此，该算法收敛速度问题和系数选取问题在一定程度上第 4 期刘海波等：一种基于变参考电平的零点凹口双重可控方向图综合改进方法423被解决.sLL但 Olen-Compton 算法所解决的问题是在固定波束指向和主副瓣电平差的前提下的任意阵列方向图综合问题，但这种规定可能是无解的，且得到的方向图也不再具有等副瓣特性.例如：希望在波束指向为时得到旁瓣峰值电平比主瓣低 dB 的均匀副瓣方向图，但如果 的值较大时，无论算法迭代多少次，都无法得到期望的方向图，而只能得到较为令人满意的方向图，且得到的方向图其实是根据最小均方误差获得的一个模糊的解，而且在这种情况下方向图的副瓣不一定是均匀的9.事实上，在进行阵列天线的方向图综合时，由于波束宽度与测角精度息息相关，往往更加关心这样一个问题：在波束的指向和主瓣宽度都确定的前提下，如何获得最低的均匀副瓣方向图.为解决该问题，刁跃龙9提出变参考电平的概念，即将式(25)修改为d 旁瓣区域内最小的峰值电压（29）d并且在每次迭代前，均对集合 进行搜寻.这种峰值搜索对于现代计算机仿真软件来说是很容易完成的.dd引入变参考电平后，随着迭代不断进行，高的副瓣总是试图与最低副瓣趋于同一水平，而参考电平又是旁瓣区域内最低的峰值电平，因此当所有的副瓣峰值与 趋于一致时，即可得到均匀等副瓣的方向图.引入变参考电平后的算法步骤可以描述为9：s(1)确定期望波束指向，并利用式(2)求出其导向矢量；(Left,Right)(Left,Right)(2)根据预设的主瓣宽度和波束指向，确定主瓣区域；(3)在主瓣区以外的副瓣区引入大量的虚拟干扰，并对其赋初始功率(为方便操作，可直接令所有干扰的功率初值为 1)；kp(,k)(4)根据式(11)和(12)，求出期望信号和无用信号的自相关矩阵，并根据 MSNR 准则计算在当前迭代次数 下的最优权值，然后利用式(1)求出阵列的方向图；d(,k)(,k)(5)根据式(29)计算当前的参考电平,并确定当前方向图的电平与期望方向图的电平差；(6)根据式(28)所述的迭代公式更新下次迭代时的干扰功率，这里要注意的是，在引入变参考电平后，主瓣区是被固定的，因此不需要重新搜索主瓣区范围；(7)重复步骤(4)(6)，直至迭代次数达到预设上限或满足以下表达式20lgmax(p(fi,k)d(fi,k)（30）wsoptwsopt式中 为所能容忍的最大副瓣峰值电平差.此时，MS-NR 准则下求得的权值 即为满足公式(30)限制的权矢量.根据式(5)中最大信噪比准则的最优权矢量表达式，就是所求的最优权矢量，这里的值可以根据阵列的峰值增益需求来确定.刁跃龙9在提出变参考电平概念后，还对其进行了一定程度的应用扩展.如为了对抗干扰，往往希望阵列天线的方向图能够在某些方向形成凹口区.针对这个问题，刁跃龙将式(29)所述的变参考电平修改为d(,k)=10L0/20d1(,k),在凹口区以外的副瓣区d1(,k),在凹口区（31）d1(,k)而此时的选择与式(27)所述不同，此时修改为d1(,k)凹口区域内最小的副瓣峰值电压（32）L0式中表示凹口深度与等副瓣电平的电平差，dB.通过这样设置变参考电平，收敛后的的方向图不仅能保持等副瓣特性，还能够在指定的区域产生凹口，且凹口深度可控.1.2 算法改进实际情况中，在抗杂波和干扰时会存在以下几种需求：(1)方向图在一个或多个固定的区域形成宽度和深度可控的凹口；(2)方向图只在某个或某几个固定的角度形成深度可控的零点，而不形成具有一定宽度的凹口；(3)方向图在某个或某几个固定区域形成宽度和深度可控的凹口，同时在某几个固定的角度形成深度可控的零点.对于需求(1)，采用式(32)所述的可变参考电平即可实现.而对于需求(2)和需求(3)，由于式(32)的搜索方法为对指定的副瓣区进行最低旁瓣峰值搜索，而零点位置只是一个或几个固定的点，因此不存在副瓣峰值，参考电平在这些位置的峰值搜索结果一定是空的.为实现零点控制，将变参考电平修改为424北 京 理 工 大 学 学 报第 44 卷d(,k)=10L0/20d1(,k),在主瓣区和零点位置以外的区域d1(,k),在零点位置（33）其中d1(,k)方向图在零点位置的电压（34）如此，便可以在指定位置产生深度可控的零点.以上方法也适用于同时产生多个深度不同的零点，只需在式(33)中对多个零点进行分别控制即可.下面以两个零点为例，给出变参考点的定义公式d(,k)=10L1/20d1(,k),在所有零点以外的副瓣区d1(,k),在零点1位置10(L1L2)/20d1(,k),在零点2位置（35）L1L2d1(,k)式中，零点 1 和零点 2 的零点深度比收敛后的正常副瓣电平分别低和，dB.此时修改为d1(,k)零点1位置的方向图电压（36）对于需求(3)，可以考虑将变参考电平修改为d(,k)=10L1/20d1(,k),在凹口和零点位置以外的副瓣区d1(,k),在凹口位置10(L1L2)/20d1(,k),在零点位置（37）d1(,k)此时与式(32)相同，为凹口区的最低峰值电平.同样地，也可以按照式(35)中对参考电平的修改方式来修改式(37)，以此同时控制多个凹口和多个零点.2 仿真验证圆弧阵是极具代表性的共形阵11，本节以均匀圆弧阵为示例对本文所述算法进行仿真验证，完整的 66 元均匀圆环阵(UCA)布阵示意方式如图 1所示.Y333412364656632313536XO图 1 均匀圆环阵列阵元排布示意图Fig.1 Schematic diagram of uniform circular array element layout 取图 1 所示圆环阵相邻的 22 个阵元(5611 号阵元)为均匀圆弧阵.对圆弧阵进行仿真实验，具体参数见表 1，算法收敛条件取式(30)，阵列的单元方向图取f(m,)=cos(m),90 m 900,其他（38）mOmX式中为以为圆心，第个阵元与轴的夹角.2.1 凹口控制选择参考电平如式(31)所述，可以实现固定主瓣宽度下最低等副瓣方向图的同时控制波束的零陷宽度和深度.如要求 1：在波束指向为 0时形成等副瓣波束，波束宽度限制为 20，且在 3050处产生比副瓣电平低 30 dB 的凹口；要求 2：在波束指向为 10时形成等副瓣波束，波束宽度限制为 20，且在4535和 3550处产生比副瓣电平分别低20 dB 和 30 dB 的凹口.采用表 1 所述的圆环阵进行仿真，凹口控制结果如图 2 所示.图 2 的仿真结果表明，方向图在均匀低副瓣的同时能将主瓣宽度保持在 20，而且在指定位置产生 表 1 圆弧阵仿真参数Tab.1 Simulation parameters of arc array参数参数值波长/m0.10阵元数22圆心角/()120半径/m0.55方向图采样间隔/()0.50方位空域范围/()9090K迭代系数1最大迭代次数200第 4 期刘海波等：一种基于变参考电平的零点凹口双重可控方向图综合改进方法425了期望宽度和深度的凹口.2.2 零点控制选择参考电平如式(33)所述，可以实现固定主瓣宽度下最低等副瓣方向图的同时控制波束的零点位置和深度.如要求 3：在波束指向为 0时形成等副瓣波束，波束宽度限制为 20，且在60，40，50和70处产生比副瓣电平低 40 dB 的零点；要求 4：在波束指向为 20时形成等副瓣波束，波束宽度限制为20，且在85、50、30和 60处产生比副瓣电平分别低 40 dB、40 dB、50 dB 和 30 dB 的零点.采用表 1所述的圆环阵进行仿真，凹口控制结果如图 3 所示.根据图 3 的仿真结果可知，方向图在均匀低副瓣的同时能将主瓣宽度保持在 20，而且在指定位置产生了期望深度的零点.80 60 40 20020406080方位角/()(a)要求3仿真结果80 60 40 20020406080方位角/()(b)要求4仿真结果初始方向图多零点方向图等副瓣电平值初始方向图多零点方向图等副瓣电平值8090706050403020100归一化增益/dB80706050403020100归一化增益/dB图 3 零点控制结果Fig.3 The results of nulls controlling 2.3 凹口与零点同时控制选择参考电平如式(37)所述，可以实现固定主瓣宽度下最低等副瓣方向图的同时控制波束的零点和凹口.如要求 5：在波束指向为25时形成等副瓣波束，波束宽度限制为 20，在60处产生比副瓣电平低 40 dB 的零点，且在 3550处产生比副瓣低 30 dB的凹口；要求 6：在波束指向为 0时形成等副瓣波束，波束宽度限制为 20，且在6035和 40处产生比 80 60 40 20020406080方位角/()80 60 40 20020406080方位角/()80706050403020100归一化增益/dB80706050403020100归一化增益/dB初始方向图单凹口方向图等副瓣电平值初始方向图单凹口方向图等副瓣电平值(a)要求1仿真结果(b)要求2仿真结果图 2 凹口控制结果Fig.2 The results of gaps controlling 初始方向图零点凹口同时存在方向图等副瓣电平值初始方向图零点凹口同时存在方向图等副瓣电平值(a)要求5仿真结果(b)要求6仿真结果80 60 40 20020406080方位角/()80 60 40 20020406080方位角/()80706050403020100归一化增益/dB80706050403020100归一化增益/dB图 4 凹口和零点同时控制的结果Fig.4 The results of nulls and gaps controlling simultaneously426北 京 理 工 大 学 学 报第 44 卷副瓣电平分别低 30 dB 和 40 dB 的凹口和零点.采用表 1 所述的圆环阵进行仿真，凹口控制结果如图 4所示.通过图 4 可知，方向图在均匀低副瓣的同时能将主瓣宽度保持在 20，而且在指定位置产生了期望深度的零点与期望宽度和深度的凹口.3 结论本文对已有的基于自适应阵原理的低副瓣方向图综合方法进行了改进.相比于原算法，新算法能够实现波束在固定主瓣宽度下同时具备等低副瓣、可控深度和宽度的凹口与可控深度的零点，这对于相控阵雷达对抗杂波和干扰有一定的实用价值.参考文献：徐峰,孙雨泽,杨小鹏,等.MIMO 雷达波束形成的低副瓣LCMV 算法 J.信号处理,2017,33(6):805 810.XU Feng,SUN Yuze,YANG Xiaopeng,et al.LCMValgorithm with low sidelobe in MIMO radar beamformingJ.Journal of Signal Processing,2017,33(6):805 810.(inChinese)1 蔡炯,王锐,胡程,等.基于新型扫描昆虫雷达的迁飞昆虫目标检测及密度反演 J.信号处理,2022,38(7):1333 1352.CAI Jiong,WANG Rui,HU Cheng,et al.Target detection anddensity inversion of migrating insects based on a novelscanning insect radarJ.Signal Processing,2022,38(7):1333 1352.(in Chinese)2 徐峰,薛珺亓,杨小鹏.相控阵 MIMO 雷达抗干扰波束形成算法 J.信号处理,2017,33(12):1531 1535.XU feng,XUE Junqi,YANG Xiaopeng,et al.Phased-MIMOradar anti-jamming beamforming algorithmJ.SignalProcessing,2017,33(12):1531 1535.(in Chinese)3 曹渊,刘威,崔东华.适用于任意阵列的鲁棒波束形成算法J.北京理工大学学报,2019,39(12):1263 1267,1297.CAO Yuan,LIU Wei,CUI Donghua.A robust beamformingalgorithm for arbitrary arrayJ.Transactions of BeijingInstitute of Technology,2019,39(12):1263 1267,1297.(inChinese)4 吕岩,曹菲,杨剑,等.基于改进广义线性组合算法的极化阵列稳健波束形成 J.北京理工大学学报,2022,42(12):1305 1311.5 L Yan,CAO Fei,YANG Jian,et al.Robust beamforming ofpolarization array based on an improved general linearcombination algorithmJ.Transactions of Beijing Institute ofTechnology,2022,42(12):1305 1311.(in Chinese)罗雪.基于未知互耦模型的波达方向估计 J.北京理工大学学报,2021,41(12):1286 1292.LUO Xue.DOA estimation based on unknown mutualcoupling modelJ.Transactions of Beijing Institute ofTechnology,2021,41(12):1286 1292.(in Chinese)6 OLEN C A,COMPTON R T.A numerical pattern synthesisalgorithm for arraysJ.IEEE Transactions on Antennas andPropagation,1990,38(10):1666 1676.7 ZHOU P Y,INGRAM M A.Pattern synthesis for arbitraryarrays using an adaptive array methodJ.IEEE Transactionson Antennas and Propagation,1999,47(5):862 869.8 刁跃龙.基于自适应阵原理的任意阵方向图综合 D.西安:西安电子科技大学出版社,2003.DIAO Yuelong.Arbitrary array pattern synthesis based onadaptive array principleD.Xian:Xidian University Press,2003.(in Chinese)9 贾可新,柳桃荣.一种基于最大信杂比的 MTD 滤波器设计算法 J.雷达科学与技术,2014,12(5):506 509.JIA Kexin,LIU Taorong.An algorithm of MTD filter basedon maximum signal-to-clutter ratioJ.Radar Science andTechnology,2014,12(5):506 509.(in Chinese)10 毕杨，甘甜.基于 Olen-Compton 的恒定响应宽带波束形成算法 J.无线电工程,2018,48(1):22 24.BI Yang,GAN Tian.Broadband constant response beamforming based on Olen-Compton algorithmJ.RadioEngineering,2018,48(1):22 24.(in Chinese)11 SHEN Zhubin,LI Jianfeng,WU Qihui.Robust low sidelobesynthesis for arbitrary arrays based on gradient ascent with min-oracleJ.IEEE Transactions on Antennas and Propagation,2023,71(1):450 462.12 SHI Xue.Pattern synthesis with closed-form array responsecontrol algorithmsJ.IEEE Antennas and WirelessPropagation Letters,2022,21(8):1507 1511.13 WANG Wenxia,YAN Shefeng,MAO Linlin,et al.Robustminimum variance beamforming with sidelobe-level controlusing the alternating direction method of multipliersJ.IEEETransactions on Aerospace and Electronic Systems,2021,57(5):3506 3519.14（责任编辑：刘芳）第 4 期刘海波等：一种基于变参考电平的零点凹口双重可控方向图综合改进方法427