均匀双圆环阵虚拟成阵方法研究.pdf

均匀双圆环阵虚拟成阵方法研究任笑莹1,王英民1,张立琛2,王奇1,廉杰1(1.西北工业大学航海学院,陕西西安,710072;2.中国船舶集团有限公司系统工程研究院,北京,100094)摘要:针对在航空声呐系统有限空间内双圆环阵因孔径受限而导致的目标分辨力和空间分辨率不足的问题,提出了一种可以实现双圆环阵孔径扩展的虚拟成阵方法。该方法利用双圆环阵模态域信号模型和双直线阵阵元域信号模型之间的关系,设计了采样矩阵和相位旋转矩阵,将均匀双圆环阵转换为均匀双直线阵。考虑了实际双圆环阵中的阵元位置误差、阵元间耦合、观察角度误差等的影响,利用最小二乘准则提出了一种将实际双圆环阵近似转换为均匀双直线阵的优化模型,并给出了迭代求解方法。与合成孔径方法相比,该方法是对基阵阵列流形的一种预处理,对使用过程中基阵的运动特性没有特殊要求;与模态域转换方法相比,该方法可以体现目标双直线阵的阵元间隔和两线阵间距,充分利用了双直线阵相对于单直线阵的优势。仿真实验和水池实验结果表明,双圆环阵的孔径得到了扩展,虚拟阵波束图比双圆环阵波束图的主瓣更窄且旁瓣更低,有效提高了目标分辨力和空间分辨率。关键词:双圆环阵;虚拟成阵;孔径扩展;空间分辨率;目标分辨力中图分类号:TJ630.34;U674文献标识码:A文章编号:2096-3920(2023)05-0735-11DOI:10.11993/j.issn.2096-3920.202206012ResearchonVirtualArrayTransformMethodofUniformBicircularArrayRENXiaoying1,WANGYingmin1,ZHANGLichen2,WANGQi1,LIANJie1(1.SchoolofMarineScienceandTechnology,NorthwesternPolytechnicalUniversity,Xian710072,China;2.CSSCSystemsEngineeringResearchInstitute,Beijing100094,China)Abstract:Inordertosolvetheproblemofinsufficienttargetdiscriminationandspatialresolutioncausedbythelimitedapertureofabicirculararrayintheaerialsonarsystem,avirtualarraytransformmethodthatcouldexpandtheapertureofthebicirculararraywasproposed.Intheproposedmethod,therelationshipbetweenthesignalmodelofthebicirculararrayinthemodaldomainandthatofthedoublelineararrayinthearrayelementdomainwasanalyzed,andthesamplingmatrixandphaserotationmatrixweredesignedtotransformtheuniformbicirculararrayintoauniformdoublelineararray.Then,theeffectsofarrayelementpositionerror,mutualcouplingbetweenarrayelements,andobservationangleerrorintheactualbicirculararraywereconsidered.Byusingtheleastsquarecriterion,anoptimizationmodelwasproposedtoapproximatetheactualbicirculararraytoauniformdoublelineararray.Theiterativesolutionmethodwasalsogiven.Comparedwiththesyntheticaperturemethod,the proposed method is a preprocessing of the array manifold and has no special requirements for the motioncharacteristicsofthearrayintheactualapplicationscene.Comparedwiththetransformmethodinthemodaldomain,theproposedmethodcanreflectthearrayelementspacingandthespacingbetweentwolineararraysofthedoublelineararrayandmakefulluseoftheadvantagesofthedoublelineararraycomparedwiththesinglelineararray.Theresultsofthesimulationexperimentandwatertankexperimentshowthattheapertureofthebicirculararrayisexpanded;themainlobeofthevirtualarraybeampatternisnarrower,andthesidelobeislowerthanthatofthebicirculararraybeampattern,whicheffectively收稿日期:2022-06-12;修回日期:2022-07-12.基金项目:国家自然科学基金项目资助(51879221,62001392).作者简介:任笑莹(1989-),女,在读博士,主要研究方向为阵列信号处理.第31卷第5期水下无人系统学报Vol.31 No.52023年10月JOURNALOFUNMANNEDUNDERSEASYSTEMSOct.2023引用格式 任笑莹,王英民,张立琛,等.均匀双圆环阵虚拟成阵方法研究 J.水下无人系统学报,2023,31(5):735-745.水下无人系统学报sxwrxtxb.xml-735improvesspatialresolutionandtargetdiscrimination.Keywords:bicirculararray;virtualarray;apertureexpansion;spatialresolution;targetdiscrimination0引言在海洋环境中,声波因其衰减速度明显小于光和电磁波,而成为海洋探测的重要手段。声呐是利用声波在水中的传播和反射特性对水下目标进行探测和通信的电子设备。针对水听器基阵的阵形和阵列信号处理的研究使声呐的探测能力不断提高。常见的水听器基阵阵形包括均匀直线阵、均匀双直线阵、均匀单圆环阵和均匀双圆环阵等不同的几何结构,这些基阵阵形不仅在声呐系统中得到了广泛应用,经典的均匀间隔分布单直线阵还被应用于导航1、通信系统2、语音处理3和目标方位估计4等研究领域。围绕单直线阵展开的研究主要包括信号处理5-6和阵形结构优化7-8两方面。针对单直线阵自身的几何结构特点导致的左右舷模糊问题,在均匀单直线的基础上增加另一条均匀单直线阵的双线阵结构不仅能提高探测距离还可以有效改善左右舷模糊的问题9。李智忠等10结合舷侧艏艉双线基阵的特点,根据互谱测向原理提出了基于互谱相关的空间谱估计方法。袁静珍11对左右线列阵进行相移处理后,利用相减技术对合成数据进行处理,降低了分辨目标左右舷对两单线阵的间距要求,扩大了同一双线阵对处理信号频带的适用范围。郭乔鹤12在浅海波导条件下,利用双线阵对舰船辐射噪声进行波束形成和空域相关处理,有效提高了对舰船目标的识别能力。双线阵虽然可以有效改善左右舷模糊问题,提高频带适用范围和目标识别能力,但是和单线阵一样无法实现全平面等精度均匀观测,因此具有圆对称结构的单圆环阵引起了学者们的注意。刘振等13针对均匀圆阵基于相位差方法的定位精度较低的问题,首先估计出远场源的方位角和俯仰角后,利用协方差矩阵差分方法提取出近场源差分矩阵,再通过改进的类旋转不变估计信号参数方法计算出近场源的方位角和俯仰角,在远场源和近场源角度相同的情况下能够有效识别混合源,提高了混合源方位估计精度。Wax 等14提出了一种均匀圆环阵利用空间平滑技术实现相干源方位估计的模式变换方法。宋昕等15利用模式转换技术将导向向量限定于确定的椭圆集合中,根据最差性能最优化的设计思想,给出了阵列加权向量的闭式解表达式。Zhao 等16利用蚁群算法和遗传算法对嵌套稀疏圆阵的阵元位置进行优化以降低波束图的旁瓣级。Huang 等17研究了均匀圆环阵的频率不变响应,并分析了同心双圆环阵中内外圆半径变化对模态域波束形成的影响,指出合理设置双圆环阵中的内外半径可以有效避免贝塞尔函数零点的影响18,然后进一步研究了同心双圆环阵的频率不变响应波束形成方法19。相对于单直线阵,双直线阵可以改善左右舷模糊的问题并提高频带适用范围和目标识别能力,单圆环阵可以实现 360全平面的均匀观察,双圆环阵可以提高单圆环阵的方位估计精度。但是在航空声呐系统中,声呐基阵搭载平台的空间往往十分有限,导致基阵的孔径受到限制。合成孔径技术能利用基阵和目标的相对匀速直线运动,实现较小尺寸基阵的孔径扩展,但是该技术对平台运动轨迹和速度要求较高,增加了其在航空声呐平台上的应用难度20。模态域转换方法可以将单圆环阵的阵列流形矩阵转换为具有范德蒙结构的矩阵21,有利于对相干信号的处理,但是模态域转换得到的阵列流形与单直线阵的阵列流形构成并不相同,无法体现单直线阵的阵元间距等物理参数的意义。受单圆环阵相位模式转换方法的启发,文中结合双直线阵相对于单直线阵的优势,考虑将均匀双圆环阵虚拟为首阵元对齐的目标均匀双直线阵,以期扩展双圆环阵的孔径。通过对均匀双圆环阵和均匀双直线阵接收信号模型的分析发现,双圆环阵的模态域信号模型和双直线阵波束域信号模型之间可以通过适当的矩阵实现近似转换,因此文中设计了采样矩阵和相位旋转矩阵,将均匀双圆环阵转换为均匀双直线阵。相对于合成孔径技术和模态域虚拟成阵算法,文中所提出的虚拟成阵算法不仅在降低了应用条件的同时实现了孔径扩展,而且充分利用了双直线阵相对于单直线阵的优势。在仿真实验和水池实验中,分别2023年10月水下无人系统学报第31卷736JournalofUnmannedUnderseaSystemssxwrxtxb.xml-采用常规波束形成(conventionalbeamforming,CBF)方法和旁瓣控制波束形成(sidelobecontrollbeam-forming,SLC)方法22计算阵列加权向量并得到波束图,结果表明虚拟阵波束图比双圆环阵波束图的主瓣更窄且旁瓣更低,有效提高了空间分辨率和目标分辨力。1虚拟成阵算法1.1理想双圆环阵r1r2121假设均匀双圆环阵中各阵元之间都不存在互耦影响,内外圆的阵元数均为 M,内外圆半径分别是 和,参考点是双圆环阵的几何中心,内外均匀单圆环中第 1 个阵元到参考点的连线与 X 轴的夹角分别是和,那么方向上的远场平面波入射到该基阵上的示意图如图 1 所示。远场平面波Y2#2#1#1#121M#M#X图1空间中均匀双圆环阵示意图Fig.1Uniformbicirculararrayinspace此时,该均匀双圆环阵的单位阵列响应向量可表示为ac(1)=exp(jkr1cos(11)exp(jkr1cos(112/M).exp(jkr1cos(112(M1)/M)exp(jkr2cos(12)exp(jkr2cos(122/M).exp(jkr2cos(122(M1)/M)(1)1=22i/M当成立时,双圆环阵中内外2 个单圆环阵相同序号的阵元到几何中心的连线相互重叠。2NN=2n+1假设目标均匀双直线阵的阵元数是,各阵元间不存在互耦,组成双直线阵的 2 个均匀单直线阵的阵元数都是,2 个单直线阵之间的距1离是 D,每个单直线阵中相邻两两阵元间的间距是 d,与两直线等间距的平行线作为 Y 轴,两单直线几何中心连线所在的直线为 X 轴,那么入射到该基阵上的远场平面波来波方向是时,它在空间中的示意图如图 2 所示。远场平面波(2n+1)#(2n+1)#2#1#d12#1#YXD2图2空间中均匀双直线阵示意图Fig.2Uniformdoublelineararrayinspace1此时,目标均匀双直线阵在方向上的单位阵列响应向量是al(1)=exp(jk(D/2)cos1)exp(jk(n)dsin1)exp(jk(D/2)cos1)exp(jk(n1)dsin1).exp(jk(D/2)cos1)exp(jkndsin1)exp(jk(D/2)cos1)exp(jk(n)dsin1)exp(jk(D/2)cos1)exp(jk(n1)dsin1).exp(jk(D/2)cos1)exp(jkndsin1)(2)D/2exp(jk(D/2)cos1)bD式中,为 2 个单直线阵和 Y 轴之间的距离。为便于书写,将用表示,则式(2)简写成al(1)=(1/bD)exp(jk(n)dsin1)(1/bD)exp(jk(n1)dsin1).(1/bD)exp(jkndsin1)bDexp(jk(n)dsin1)bDexp(jk(n1)dsin1).bDexp(jkndsin1)=Q(1)alt(1)(3)其中2023年10月任笑莹,等:均匀双圆环阵虚拟成阵方法研究第5期水下无人系统学报sxwrxtxb.xml-737Q(1)=1/bD.1/bDbD.bD(4)alt(1)=exp(jk(n)dsin1)exp(jk(n1)dsin1).exp(jkndsin1)exp(jk(n)dsin1)exp(jk(n1)dsin1).exp(jkndsin1)(5)那么利用空域离散傅里叶变换(discreteFouriertransform,DFT)原理21可以构造如下矩阵B=G/M00G/M(6)G0(2h+1)M其中,矩阵和 的维度都是,且G=1 exp(j2(h)M)exp(j2(h)2M)exp(j2(h)(M1)M).1 exp(j2(1)M)exp(j2(2)M)exp(j2(1)(M1)M)11111 exp(j2(1)M)exp(j2(2)M)exp(j2(M1)M).1 exp(j2(h)M)exp(j2(2h)M)exp(j2(M1)h)M)(7)利用矩阵 B 可得at(1)=Bac(1)(8)at(1)2(2h+1)1 at(1)其中的维度是,中的第 m个元素是atm(1)=1MM1i=0exp(j2(m(h+1)iM)exp(jkr1cos(112iM),1 m (2h+1)1MM1i=0exp(j2(m(3h+2)iM)exp(jkr2cos(122iM),2h+2 m 2(2h+1)(9)atm(1)=1M+n=jnJn(kr1)exp(jn(11)M1i=0exp(j2(m(h+1)iM)exp(j2niM),1 m 2h+11M+n=jnJn(kr2)exp(jn(12)M1i=0exp(j2(m(3h+2)iM)exp(j2niM),2h+2 m 2(2h+1)=M+=jM+(m(h+1)JM+(m(h+1)(kr1)exp(j(M+(m(h+1)(11),1 m 2h+1M+=jM+(m(3h+2)JM+(m(3h+2)(kr2)exp(j(M+(m(3h+2)(12),2h+2 m 2(2h+1)(10)atm(1)exp(jkr1cos(112iM)exp(jkr2cos(122iM)Jn式中:由式(9)中的和进行傅里叶级数展开获得;为 n 阶第一类贝塞尔函数;为整数。n krJn(kr)(ekr/2n)n2nM kr1M kr2JM+(m(h+1)(kr1)Jm(h+1)(kr1)当时,第一类贝塞尔函数可以用23近似表示。因此,当和成立时,式(10)中 取 0 时的贡献远远大于取其他值时的贡献,此时、JM+(m(3h+2)(kr1)Jm(3h+2)(kr1)成立,式(10)可以表示为atm(1)Mjm(h+1)Jm(h+1)(kr1)exp(j(m(h+1)(11),1 m 2h+1Mjm(3h+2)Jm(3h+2)(kr2)exp(j(m(3h+2)(12),2h+2 m 2(2h+1)(11)构造对角矩阵2023年10月水下无人系统学报第31卷738JournalofUnmannedUnderseaSystemssxwrxtxb.xml-J=diagexp(j(m(h+1)1)Mjm(h+1)Jm(h+1)(kr1),1 m 2h+1diagexp(j(m(3h+2)2)Mjm(3h+2)Jm(3h+2)(kr2),2h+2 m 2(2h+1)(12)at(1)将转换为av(1)=Jat(1)(13)其中av(1)=exp(j(h)1)exp(j(h1)1).exp(j(h1)1)exp(jh1)exp(j(h)1)exp(j(h1)1).exp(j(h1)1)exp(jh1)(14)将式(8)代入式(13)可以得到av(1)=JBac(1)(15)alt(1)av(1)kdnsin11alt(1)av(1)altd=exp(jkndsin1)avd=exp(jh1)向量和中不涉及角度的成分分别是和 h,由角度决定的成分分别是和。为了便于描述,分别记向量和中的最后一个元素为和。altdavd当仅考虑角度成分时,将和分别记为altd=exp(jsin1)(16)avd=exp(j1)(17)三角函数和指数函数之间存在如下关系23sin1=exp(j1)exp(j1)2j(18)因此altd=exp(exp(j1)exp(j1)2)(19)又由于exp(j1)=limm(1+j1m)m=limm1+(m1)(j1m)1+(m2)(j1m)2+(mm1)(j1m)m1+(mm)(j1m)m(20)exp(j1)=limm(1+j1m)m=limm1+(m1)(j1m)1+(m2)(j1m)2+(mm1)(j1m)m1+(mm)(j1m)m(21)mx=m!x!(mx)!x=1,2,m式中,。将式(20)、(21)代入式(19),则altd=explimm(m1)(j1m)1+(m3)(j1m)3+(m2+1)(j1m)2+1=exp(j1)explimm(m3)(j1m)3+(m5)(j1m)5+(m2+1)(j1m)2+1=avdexplimm(m3)(j1m)3+(m5)(j1m)5+(m2+1)(j1m)2+1(22)altdavd式(22)表明和之间可以通过相位旋转实现相互转换。1定义对应于角度的相位旋转量为1=explimm(m3)(j1m)3+(m5)(j1m)5+(m2+1)(j1m)2+1(23)则式(22)可以写为altd=avd1(24)altav因此,仅考虑角度成分的向量和之间的关系可以写为alt=av1(25)altdavd当不考虑角度成分时,记元素和为2023年10月任笑莹,等:均匀双圆环阵虚拟成阵方法研究第5期水下无人系统学报sxwrxtxb.xml-739 altd=exp(jknd)(26)avd=exp(jh)(27)altd avdh/n=kd要想使和之间实现相等关系,需要满足条件,n 和 h 都是整数,因此h=round(kdn)(28)round()式中,表示四舍五入。alt(1)av(1)P由式(28)可知是由中的某些元素组成的子向量,此时必然存在一个采样矩阵 可以实现 alt(1)P av(1)(29)P2(2n+1)2(2h+1)式中,矩阵 的维度是,且P=100|zround(kd)100|zround(kd)1.1(30)alt(1)av(1)结合式(25)和式(29),发现向量和存在如下关系alt(1)Pav(1)1(31)将式(15)代入式(31),则alt(1)PJBac(1)1(32)联立式(3)和式(32),则al(1)Q(1)PJBac(1)1(33)Q1Yal(1)式中,是一个随着声源信号来波方向的不同而发生变化的对角矩阵。期待可以找到一个与方向无关的固定矩阵来近似,那么式(33)变为al(1)YPJBac(1)1(34)当有多个观察角度时,不同方向上的单位阵列响应向量可以组成阵列流形矩阵。双圆环阵的阵列流形矩阵为Ac=ac(1),ac(2),ac(s)(35)均匀双直线阵的阵列流形矩阵为Al=al(1),al(2),al(s)(36)AcAl结合式(34)可以发现和之间存在如下关系Al TAc(37)T2N2M=diag(i)(i)i=1,2,s式中:矩阵 的维度是;相位旋转矩阵的对角元素的绝对值都是1,。1.2实际双圆环阵在航空声呐系统中,双圆环阵中各个阵元间的互耦往往难以避免,尤其是内外圆之间两两阵元间的距离相对于波长一般都很小20,因此 2 个单圆环中两两对应阵元之间的互耦不容忽视;由于制造工艺、测量误差等不可避免的因素,双圆环阵列中各个阵元的实际位置通常也会和理想位置之间存在偏差;在测量实际阵列流形时,实际的来波方向和预设方向之间经常会因为操作不够精确等原因导致实测阵列流形中单位阵列响应的方向存在误差。针对以上各种情况,接下来分析双圆环阵实际阵列流形和目标均匀双直线阵的阵列流形之间存在的关系。当双圆环阵中各阵元间存在互耦时Al TC1Ac(38)Ac=CAcC式中:为存在互耦时的双圆环阵列流形矩阵;为互耦矩阵。当双圆环阵中的阵元位置存在误差时Al TF1Ac(39)Ac=FAcF2M2M式中:为存在位置误差的双圆环阵列流形矩阵,为一个维度的矩阵。当观察方向存在误差时Al TAcU(40)Ac=AcU1Uu(i)式中:为观察方向存在误差时的双圆环阵列流形;矩阵是一个对角阵,其对角元素的绝对值都是 1,作用是实现阵列流形中各个方向导向向量的相位调整。U已知矩阵也是一个仅用来实现相位调整的对角阵,因此矩阵和可以进行合并,式(40)被写为Al TAcL(41)Lss式中,矩阵 是一个维的对角阵,对角元素的绝对值都是 1。AcAl观察式(38)(41)发现,综合考虑阵元间互耦、阵元位置误差和观察角度误差等各方面的影响时,实际双圆环阵的阵列流形和目标均匀双直线阵的阵列流形之间存在如下关系Al DAcL(42)根据最小二乘准则14,构造一个优化模型minD,L?Al DAcL?(43)2023年10月水下无人系统学报第31卷740JournalofUnmannedUnderseaSystemssxwrxtxb.xml-D2N2M式中,矩阵的维度是。DLLL优化问题(43)中需要求解的矩阵有和,由于矩阵 是一个元素绝对值均为 1 的对角矩阵,因此可以初始化 为单位矩阵L0=I(44)LzLD当固定当前矩阵作为 的最优值时,可以得到矩阵此次迭代的最优解Dz+1=AlL1zA1c(45)DzDL当固定当前矩阵作为的最优值时,可以得到矩阵 此次迭代的最优解中各个元素的取值lz+1(i)=exp(jarg(Dzac(i)1al(i)(46)arg()DL式中,表示复数辐角。通过不断地交叉迭代优化24,当满足迭代终止条件时,即可得到优化问题(43)的最优解矩阵和。2仿真实验r1r2假设航空声呐系统中均匀双圆环阵的内圆半径 和外圆半径分别是 0.45m 和 0.50m,单个圆环阵包含的阵元数都是 12,远场窄带声源信号频率为 3000Hz,水下声信号的传播速度是 1500m/s,双圆环水听器基阵的阵元间不存在互耦影响,阵列流形中各导向向量中不存在方向误差,各个阵元均处于理想位置。利用文中所提的虚拟成阵算法可以将双圆环阵转换为目标均匀双直线阵,但是在转换过程中引入了各种近似,因此双圆环阵转换得到的虚拟阵和目标均匀双直线阵之间必然存在一定的差异。为了直观展示双圆环阵利用所提虚拟成阵算法得到的虚拟阵和目标均匀双直线阵之间存在一定的误差以及虚拟阵对双圆环阵的改善效果,仿真实验中针对均匀双圆环阵、目标均匀双直线阵和虚拟阵分别利用 CBF 方法和 SLC 方法进行波束形成,并对比了利用这 2 种波束形成方法所得到的波束图之间的差异。CBF 方法的阵列加权向量为wCBF=a(0)a(0)Ha(0)(47)a(0)式中,为期望方向上的导向向量。SLC 方法的阵列加权向量可以通过求解如下优化模型得到min wHRws.t.wHa(0)=1?wHa(i)?(48)a(i)式中,为旁瓣区域各个方向上的导向向量；为期望响应。wA在已知阵列加权向量 和阵列流形 的情况下,可以得到窄带波束形成器的波束响应p=wHA(49)2N=2(2n+1)n (M1)/2双圆环阵在进行虚拟阵转换的过程中需要选择合适的目标均匀双直线阵,如果目标均匀双直线阵的阵元数过大,那么双圆环阵的虚拟阵和目标均匀双直线阵之间的单位阵列响应差异就会过大。为了将如图 1 所示的均匀双圆环阵虚拟为目标均匀双直线阵,在选择决定均匀双直线阵阵元数的参数 n 时,根据经验可知若目标双直线阵中每个均匀单直线阵的阵元间距都等于半波长,参数 n 的取值需要满足不等式,具体的取值由双圆环阵的阵元数、内外圆半径、目标均匀双直线中 2 个单直线阵之间的距离及 2个单直线阵中各阵元间的距离等参数共同决定。2.1理想双圆环阵1,22i/M1=102i/M2=202i/M 10=18020=170为了更具一般代表性,仿真实验中假设均匀双圆阵的(其中,),以双圆环阵的几何中心为参考点,该均匀双圆环阵中各个阵元在空间中的分布情况如图 3 所示。0.500.50.50X/mY/m0.5图3仿真实验中均匀双圆环阵阵元位置Fig.3Element position of uniform bicircular array insimulationexperiment在文中的仿真实验中,目标均匀双直线阵中2 个单直线阵之间的距离等于双圆环阵内外圆半径之差 0.05m,单个均匀直线阵中两两阵元的间距2023年10月任笑莹,等:均匀双圆环阵虚拟成阵方法研究第5期水下无人系统学报sxwrxtxb.xml-741都是半波长 0.25m,参数 n 取值为 3,目标均匀双直线阵的阵元在空间中的分布情况如图 4 所示。0.500.50.200.2Y/mX/m图4目标均匀双直线阵阵元位置Fig.4Elementpositionofuniformdoublelineararray图 4 中,目标均匀双直线阵中的阵元数是 14,孔径是 1.5m。已知图 3 中均匀双圆环阵的孔径尺寸是 1m,双圆环阵的虚拟阵是目标双直线阵的近似,因此利用文中所提虚拟成阵算法扩展了双圆环阵的孔径。均匀双圆环阵、目标均匀双直线阵和利用文中虚拟成阵算法得到的虚拟阵分别利用 CBF 方法和 SLC 方法进行波束形成,得到的波束结果如图 5和图 6 所示。020406080100500方位角/()双圆环阵双直线阵虚拟阵波束/dB50100图5仿真实验中理想双圆环阵 CBF 波束图Fig.5BeampatternofidealbicirculararrayinsimulationexperimentusingCBFmethod图 5 中,利用 CBF 方法得到的虚拟阵波束图和目标均匀双直线阵的波束图几乎完全重合。图 6中,利用 SLC 方法得到波束图中,目标均匀双直线和虚拟阵波束图仅存在很少的不重叠部分。这说明虚拟阵阵列流形和目标均匀双直线阵的阵列流形之间的差异非常小,虚拟阵可以看成是带有微小误差的目标均匀双直线阵。图 5 中,与理想双圆环阵相比,虚拟阵的波束图主瓣更窄且旁瓣级更低。图 6 中,在获得相同旁瓣级的情况下,虚拟阵的主瓣宽度比理想双圆环阵的明显更窄。2.2带有误差的双圆环阵5%d,5%d0.5,0.5在如图 3 结构的双圆环阵上对阵元位置添加的随机误差,同时在观察角度方向上添加的随机误差,目标双直线阵如图 4所示,考察同时存在阵元误差和观察角度误差时文中所提算法的性能。利用 CBF 方法和 SLC 方法进行波束形成的结果如图 7 和图 8 所示。0204060100500方位角/()双圆环阵双直线阵虚拟阵波束/dB50100图7仿真实验中有误差双圆环阵的 CBF 波束图Fig.7Beam pattern of bicircular array with errors insimulationexperimentusingCBFmethod图 7 和图 8 中,虚拟阵波束图和目标双直线阵波束图的主瓣和旁瓣级都相同。对比虚拟阵和有误差双圆环阵的波束图,得到了和图 5、图 6 一样的结论,证明了所提算法对带有误差的双圆环阵是有效的。综合考虑图 5图 8 可知,不管双圆环阵中是否存在误差,利用所提算法得到的虚拟阵都可以020406080100500方位角/()双圆环阵双直线阵虚拟阵波束/dB50100图6仿真实验中理想双圆环阵 SLC 波束图Fig.6BeampatternofidealbicirculararrayinsimulationexperimentusingSLCmethod2023年10月水下无人系统学报第31卷742JournalofUnmannedUnderseaSystemssxwrxtxb.xml-有效提高双圆环阵的波束性能。3水池实验为了验证实际应用中所提虚拟算法的有效性,在消声水池开展了相关的实验。该消声水池长20m、宽 8m、深 7m,利用圆锥形橡胶尖劈板实现 6 面消声。水池实验中采用由 24 个水听器构成的双圆环阵列,内外圆环的半径分别是 0.45m 和 0.50m,其阵元在空间中的位置如图 9 所示。0.500.50.500.5Y/mX/m图9水池实验中实际双圆环阵阵元位置Fig.9Elementpositionofactualbicirculararrayinwatertankexperiment该双圆环阵水平布放于水下 3m 处,发射声源和双圆环阵位于同一深度。距离双圆环阵几何中心约 5.5m 处的声源信号是频率为3000Hz 的正弦波脉冲信号,脉冲宽度为 50ms,周期为 2s,信号采集电路的采样频率为 20kHz,水池实验系统如图 10 所示。信号源功率放大器发射换能器同心双圆环接收阵列滤波放大电路采集存储电路水空气转台图10水池实验系统Fig.10Watertankexperimentalsystem在水池实验中,通过人为操作转台带动双圆环阵以 2的角度间隔旋转 180,利用采集存储电路对每个角度上各个阵元的接收信号进行采集和存储,之后进行离线处理。在采集的脉冲信号中截取 400 个有效采样快拍,利用鄢社锋25给出的计算方法获取实验中同心双圆环阵的实测阵列流形,之后利用文中所提的虚拟成阵方法得到虚拟阵的阵列流形。与仿真实验中一样,水池实验中的目标均匀双直线阵如图 4 所示,针对实际双圆环阵、目标均匀双直线阵和虚拟阵分别用 CBF 方法和SLC 方法进行波束形成,得到的结果如图 11 和图 12 所示。0204060100500方位角/()双圆环阵双直线阵虚拟阵波束/dB50100图8仿真实验中有误差双圆环阵的 SLC 波束图Fig.8Beam pattern of bicircular array with errors insimulationexperimentusingSLCmethod020408060100500方位角/()双圆环阵双直线阵虚拟阵波束/dB50100图11水池实验中实际双圆环阵 CBF 波束图Fig.11BeampatternofactualbicirculararrayinwatertankexperimentusingCBFmethod2023年10月任笑莹,等:均匀双圆环阵虚拟成阵方法研究第5期水下无人系统学报sxwrxtxb.xml-743观察图 11 和图 12 发现,虚拟阵和目标均匀双直线阵的波束图主瓣宽度和旁瓣级相当。图 11中,虚拟阵比实际双圆环阵的主瓣更窄,同时旁瓣更低。图 12 中,在虚拟阵和实际双圆环阵达到相同旁瓣级的情况下,虚拟阵波束图的主瓣宽度明显更窄。这与仿真实验中的结果一致,证明了在消声水池实验环境中,文中所提虚拟成阵算法的有效性。在仿真实验和水池实验中,虽然虚拟阵的阵元数少于双圆环阵的阵元数,但是虚拟阵的阵列孔径大于双圆环阵,保证了虚拟阵拥有更高的空间分辨率。4结束语文中针对航空声呐系统中双圆环阵因孔径受限而导致的目标分辨力和空间分辨率不足的问题,提出了一种可以实现双圆环阵孔径扩展的虚拟成阵方法。分析了双圆环阵模态域信号模型与均匀双直线阵阵元域信号模型之间存在的固有关系,通过设计可以实现维度变换的采样矩阵和角度转换的相位旋转矩阵,将均匀双圆环阵转换为了均匀双直线阵。利用最小二乘准则提出了一种可以将受到各种误差影响的实际双圆环阵近似为均匀双直线阵的优化模型,并给出了迭代求解方法。与合成孔径技术相比,文中方法是一种对阵列流形的预处理,在实际使用中对基阵的运动轨迹和速度没有特殊要求;与目前已有的模式转换方法相比,文中方法得到的虚拟阵不仅可以体现双直线阵的空间特性,有效利用双直线阵相对于单直线阵的优势,还可以有效扩展双圆环阵的孔径。仿真实验中,分别将内外圆中各阵元不在相同径向上分布的理想双圆环和带有误差的双圆环阵转换为目标均匀双直线阵,证明了所提方法对内外圆中阵元位置的对应关系没有特殊要求。水池实验中将实际双圆环阵近似转换为目标均匀双直线阵,证明了所提算法在实际应用中的有效性。仿真实验和水池实验中,双圆环阵的孔径都得到了扩展,有效改善了波束效果,提高了双圆环阵的空间分辨率和目标分辨力,可为双圆环阵在实际应用中的孔径扩展提供理论依据和实现方法。未来将基于文中工作对以下问题做进一步研究:在实际海洋环境中所提算法的有效性;针对宽带信号,目标均匀双直线阵最优几何参数的合理设置;文中所提算法对多圆环阵(圈数2)孔径扩展的有效性。参考文献:QinL,LiCL,DuYX,etal.DOAestimationandmutu-alcouplingcalibrationalgorithmforarrayinplasmaen-vironmentJ.IEEE Transactions on Plasma Science,2020,48(6):2075-2082.1JounyI.BeamspaceprocessingofcompressivelysampledsignalsreceivedviaauniformlineararrayC/IEEEIn-ternationalSymposiumonAntennasandPropagationandNorthAmericanRadioScienceMeeting.Montreal:IEEE,2020:573-574.2HuangGP,WangYZ,BenestyJ,etal.Combineddiffer-entialbeamformingwithuniformlinearmicrophonear-raysC/IEEE International Conference on Acoustics,SpeechandSignalProcessing(ICASSP).Toronto:IEEE,2021:781-785.3WeiJ,ChenSW,ChangAC.EfficienthybridDOAes-timationformassiveuniformlineararrayJ.IEICETran-sactionsonFundamentalsofElectronicsCommunicationsandComputerSciences,2019,102(5):721-724.4ChenPC,VaidyanathanPP.Convolutionalbeamspacefor l