测向天线基线选择对测量精度影响的研究_王晨航.pdf

试验与评估航天电子对抗2023年第 3期测向天线基线选择对测量精度影响的研究王晨航1，高加健2，徐璐1，张一嘉1（1.浙江理工大学信息科学与工程学院，浙江 杭州 310018；2.四川九洲电器集团有限责任公司，四川 绵阳 621000）摘要：无线电测向应用广泛，常用的天线阵列类型有线阵、圆阵和椭圆阵等。设计天线阵列时，需要考虑天线尺寸、阵列构成等因素对测向精度的影响，以保证系统性能指标。首先，构建了针对五元线阵和五元均匀圆阵的干涉测向模型，采用了一种基于平行虚拟基线选择完成测向解模糊的方法，将相位测量值与测向理论数据库进行相关，相关最大值所对应的角度即为来波方向。其次，分别分析了测向过程中的频率灵敏度和角度灵敏度所引起的测向误差。最后，通过 MonteCarlo 仿真，分析了相位误差在 15的标准偏差条件下的测向误差，并获取测向精度分布图，为实际应用时天线阵列设计提供参考。关键词：无线电测向；阵列基线；相关法；测量精度中图分类号：TN971+.5 文献标志码:AResearch on the influence of baseline selection of directionfinding antennas on measurement accuracyWang Chenhang1，Gao Jiajian2，Xu Lu1，Zhang Yijia1(1.School of Information Science and Engineering，Zhejiang Sci-Tech University，Hangzhou 310018，Zheji-ang，China；2.Sichuan Jiuzhou Electrical Appliance Group Co.,Ltd，Mianyang 621000，Sichuan，China)Abstract：Radio direction-finding has wide application，and commonly used antenna arrays include linear ar-rays，circular arrays，elliptic arrays，etc.When designing antenna array，it is necessary to consider the influence of antenna size，array composition and other factors on direction-finding accuracy to ensure the system performance index.Firstly，an interference direction-finding model for a five element linear array and a five element uniform cir-cular array is constructed.A method based on parallel virtual baseline selection is used to achieve direction-finding deblurring.The phase measurement values are correlated with the direction-finding theoretical database，and the angle corresponding to the maximum correlation value is the direction of the incoming wave.Secondly，the direc-tion-finding errors caused by frequency sensitivity and angle sensitivity in the direction-finding process are ana-lyzed.Finally，by means of MonteCarlo simulation，the direction-finding error under the condition of 15 standard deviation of phase error is analyzed，and the direction-finding accuracy distribution map is obtained，which provides reference for antenna array design in practical application.Key words：radio direction-finding；array baseline；correlation method；measurement accuracy0 引言近年来，无线通信技术正在快速发展，无线电侦察测向以及精准定位正逐渐成为影响战争胜负的关键因素。无线电测向是在空间中利用接收的信号信息来测定辐射源信号的到达方向，属于无源测向技术1。无源测向是指由测向系统对辐射源进行测向时，测向设备并不需要通过自身先发射一个助于测向的电磁波信号，所以无线电测向的隐蔽性更强2。随着测向设备在海、陆、空等平台的广泛应用，由阵列应用环境不同导致的阵列类型更加多样，本文针对不同的项目需求，通过仿真分析天线阵列分布下的天线基线情况以求得更加良好的测向精度，为工程应用提供支持。收稿日期：20230225；20230402修回。作者简介：王晨航（2000-），男，硕士，主要研究方向为无线电测向。-15DOI:10.16328/j.htdz8511.2023.03.0032023，39（3）航天电子对抗对于测向系统来说，影响测向精度的原因很多，包括外部原因如噪声干扰和障碍物导致的反射波，内部原因如自身通道误差、天线设备互耦等。姜志成分析了干涉仪测向的误差情况，并提出了相位校准方法，但未对不同天线阵列排布对测向系统的影响进行分析3；吴启星等人结合了卡尔曼滤波对误差分析模型进行了优化，使误差分析更为准确真实，但是该研究的测向模型没有考虑一般情况下虚拟基线的使用4。段国文等人对常用天线阵列的测向精度情况进行了分析5，党浩淮等人在其基础上增加了天线相位误差对天线测向的关系6，但是均缺少了阵列排布间距对测向精度误差的影响分析。设计天线阵列时，需要综合考虑天线尺寸、可应用频段等因素，选择合适的天线阵元和分布间距，以满足系统需求7。基于此，本文构建了针对五元线阵和五元均匀圆阵的干涉测向模型，主要工作如下：1）采用基于平行虚拟基线选择完成测向解模糊的方法，将相位测量值与测向理论数据库进行相关查找，相关最大值所对应的角度即为来波方向；2）分析了测向过程中的频率灵敏度和角度灵敏度所引起的测向误差；3）通过 MonteCarlo 仿真，分析了相位误差在 15的标准偏差条件下的测向误差，绘制测向精度分布图，分析系统测向精度。1 阵列测向原理1.1非均匀线阵非均匀线阵是利用电磁波在不同天线阵元之间的波程差与来波方向的关系，进一步求得信号的来波方向，一般来说，辐射源到天线阵列的距离和信号波长相比足够长，因此可以将信号电磁波看作平行波入射。非均匀线阵天线阵元之间距离位置设置十分灵活，在实际应用时更为方便8，因此，在现实中的测向系统中，通常适用于安装平台受限或者应用较窄频段范围的情况。非均匀线阵天线模型如图 1所示。图 1 中 n 为天线阵元个数，d1、d2等分别为其余天线阵元与天线 1之间的间距，为入射波方位角。通过理论模型计算各天线阵元之间相位差，首先可求得天线 2和天线 1之间的相位差为：2，1=2 d1cos（1）进一步可得天线n与天线 1之间的相位差为：n，1=2 dn-1cos（2）由上面的计算方法可以得到一组各天线与天线 1之间的相位差样本库，即顶点基线。由上述公式可灵活变换得到任意 2根天线阵元之间相位差，即：m，n=2(dm-1-dn-1)cos m，n 1（3）对式（3）变换求取得：=arccos(m，n/2(dm-1-dn-1)（4）1.2均匀圆阵均匀圆阵测向天线阵列模型如图 2所示。图 2中n为天线阵元个数，r为天线阵元半径，为相邻阵元夹角，为入射波角度。均匀圆阵模型构建后，可以通过不同阵元之间相位差计算入射波的到达方向，通过计算阵元与参考圆心点之间的相位差，再通过数学计算可得到不同阵元之间的相位差。首先计算出入射波到达第n个阵元和圆心之间的波程差为：d=r cos(-(n-1)（5）进一步可得到天线n相对于圆心的相位差为：dpn，0=2 r cos(-(n-1)（6）式中，为信号的波长。通过数学关系可以求得第n个阵元与第m个阵元之间的相位差为：dpn，m=-2 r sin(-(n+m-2)/2)sin(m-n)/2)（7）均匀圆阵常用的测向方法有相关法和最小二乘法，相关法的基本原理是：构建无测量误差的理论相位差样本库，将实测相位差矢量在样本库中进行相关查找，具有最大相关系数的矢量对应角度即视为入射图 1非均匀线阵模型图 2均匀圆阵模型图-16试验与评估王晨航，等：测向天线基线选择对测量精度影响的研究2023，39（3）波的到达方向9。2 基线选择和相位解模糊2.1基线分布阵列基线是指天线阵元之间组成的连线，一个阵列之中存在多条基线，在使用时根据需要可灵活选择。对于非均匀线性阵列来说，基线情况可以直观看出，只要获取了各天线之间的距离即可计算得出基线长度。以非均匀五元线阵为例，基线示意图如图 3所示。短基线组合由X1X2、X2X3、X3X4、X4X5和X5X1 组成，顶点基线组合有X1X2、X1X3、X1X4、X1X5组成，此外 还 存 在 一 些 基 线 全 部 列 入 全 基 线 组 合 中，包 括X1X2、X2X3、X3X4、X4X5、X5X1、X1X3、X1X4、X2X4、X2X5、X3X5。以均匀五元圆阵为例，基线示意图如图 4所示。短 基 线 组 合 包 括A1A2、A2A3、A3A4、A4A5、A5A1，长 基 线 组 合 包 括A1A3、A2A4、A3A5、A4A1、A5A2以 及 全 基 线 组 合A1A2、A2A3、A3A4、A4A5、A5A1、A1A3、A2A4、A3A5、A4A1、A5A2。2.2解相位模糊由式（1）可知，为求得唯一方位角，需满足和一一对应，则有：2d cos =d 1 2（8）当不满足该条件时，待测相位差不是真实相位差，与之相差了若干个周期2，即相位模糊。则解相位模糊就是计算出相差的周期数将其还原为真实相位差，基本流程如图 5所示。解算时，需要以 1 根长度小于半波长的基线为基准，其对应相位差无相位模糊。线阵中所有基线均平行，以 2根同向基线为例，如X1X2和X2X3，长度分别为d1和d2，对应相位差分别为1和2，假设d1 2，则有：2d1/sin =12d2/sin =2+2k（9）式中，k=0，1，2。对式（9）求解可得：k=|d21-d122d1|Z（10）式中，Z为取整操作。求得k值后代入上式求得相位差?2，以此求得方位角。实际应用时，得到的相位差包含一定的测量误差，根据式（10）计算出的k值实际是：k?=d2?1-d1?22d1=d21-d122d1+d21-d122d1=k+k（11）式中，1、2为测量误差。当|k|小于 0.5 时，测量误差不影响解相位模糊，通过统计分析可得，当d2d1越小时，解模糊成功率越高，因此，可以利用不同基线组合逐次解相位模糊，提高解相位模糊概率。假设所有基线均大于半波长，如基线X3X4、X4X5，对应相位差分别为1和2，长度分别为d1=n1d0和d2=n2d0，且d0 d22（16）进一步有：d3=d2-d3=2d2-d1 0.382d2 d32（18）另有：d4=d3-d4=5d2-3d1 0.382d3 2（19）可得出d3和d4对应相位差为：|2(2d2-d1)/sin =2d3 sin =2(2h2-h1)+(22-1)2(5d2-3d1)/sin =2d4 sin =2(5h2-3h1)+(52-31)（20）为使式（20）有唯一解，对其变换可得：|2h2-h1+(22-1)/(2)+0.5Z=05h2-3h1+(52-31)/(2)+0.5Z=0（21）存在|2-15-3=-1 0，则式（21）中对于h1和h2存在唯一解，因此利用基线d1和d2可以对基线正对方位角90范围内的相位解模糊，进一步利用 5组长基线和短基线可以对全方位相位解模糊，在满足d 2的条件下，可适当增大基线长度，以提高测向精度。3 仿真研究首先，根据输入参数包括尺寸、类型、频率范围等，构建不同的天线模型。其次，在一定的频率范围内，以一定的频率步进和角度步进构建相位差样本库，其中包含每个方位角对应的相位差矢量数据，每个矢量包含阵列中全部基线对应相位差。进行精度仿真测试时，以样本库中的每个相位差矢量为导向矢量，为其增加均方根误差固定的相位误差，并将其归一化到内作为待测相位差，经过解模糊后将其与样本库中理论数据进行相关运算，并求其结果误差。重复以上步骤进行 100 次测试，最后求其均方根差作为最终结果，具体流程如图 6所示。3.1非均匀线阵以非均匀五元线阵为例，相邻阵元间距长度分别为 29 cm、42 cm、71 cm、97 cm，仿真分析不同基线组合对解算结果误差的影响。具体地，分别以短基线组合、长基线组合和全基线组合进行逐次解相位模糊，再进行方位角解算。短基线组合测向精度结果如图 7所示。由图 7 可看出，随着频率增大，测向精度逐渐减小，并且从线阵法线位置入射来波计算误差最小，从线阵水平方向入射来波计算误差较大，进一步得到在测试条件下，最大测向精度为短m=8.794，法线入射方位最大测向精度为法m=2.341，可用方位角范围图 6相关法测试流程图 7线阵短基线测向精度图-18试验与评估王晨航，等：测向天线基线选择对测量精度影响的研究2023，39（3）在法线60范围内。顶点基线组合测向精度结果如图 8所示。图 8中最大测向精度为短m=7.496，法线入射方位最大测向精度为法m=2.094，整体趋势与短基线组合相似，整体误差要小于短基线组合。全基线组合测向精度结果如图 9所示。图 9中最大测向精度为短m=8.051，法线入射方位最大测向精度为法m=2.243，全基线与顶点基线误差结果相近，均小于短基线误差结果。由仿真结果可看出，全基线测量精度随着频率的增大而逐渐降低，对基线相位差进行解模糊时，选取不同的基线组合会影响最后的计算结果，全基线和顶点基线的计算误差相近，短基线误差高于全基线和顶点基线。3.2均匀圆阵构建样本库时，频率步进为 5 MHz，其将产生最大为 2.5 MHz的测频误差，由式（6）对f求导可得：=2f/cr sin(-(n+m-2)2)sin(m-n)/2)（22）当f取 2.5 MHz时，分别计算全方位角范围短基线组合和长基线组合对应误差，通过 Matlab变量空间查询可得最大测频误差分别为短=2.116 0、长=3.423 8。当角度步进为 1时，其将产生最大为 0.5的测角误差，由式（6）对求导可得：=2f/cr cos(-(n+m-2)/2)sin(m-n)/2)（23）以f和为变量，取为 0.5，计算对应测角误差，如图 10所示。由图 10 可看出，在频率一定时，误差随方位角增大呈现正弦性变化；在方位角一定时，误差随频率增大呈现线性变化。给定均匀五元圆阵半径为 50 cm，频率步进为5 MHz，角度步进为 1，相位误差均方根误差为 15，将相位差数据引入相位误差后，依据上述方法进行解模糊，获取每个基线对应的真实相位差数据，再将其与样本库中对应相位差通过相关运算求取方位角，与真实方位角对比得到测量误差值。在本文研究中，短基线组合与长基线组合通过相关运算得到方位角精度，分别对其进行仿真分析，结果如图 11所示。图 11（a）为短基线组合测量精度，最大测量精度为max=3.619，随着频率增大，整体精度随之降低，由于测角误差已经超过引入的相位误差，以最大测角误差作为当时频点的引入误差，频率超过阈值范围后计算误差会逐渐变大。图 11（b）中长基线最大测量精度为max=2.704，整体趋势和短基线解算结果相似，整体误差要小于短基线组合。由仿真结果可以看出，圆阵利用短基线或者虚拟基线进行解模糊，再选择基线组合进行相关运算，长基线测量精度要优于短基线解算，频率逐渐增大时，测量误差逐渐变小，但随着测角误差的增大，测量精度逐渐变差。图 8线阵顶点基线测向精度图图 9线阵全基线测向精度图图 10测角误差示意图-192023，39（3）航天电子对抗4 结束语本文首先介绍了以非均匀线阵和均匀圆阵为阵列天线的方位角解算方法，然后阐述了阵列基线的构成，详细分析了以平行基线解相位模糊的原理，并对不同阵列测向精度进行 Matlab仿真分析，绘制了测向精度分布图。仿真结果表明，非均匀线阵采用全基线组合和顶点基线组合测量精度相近，均高于短基线组合；均匀圆阵测量精度随着频率增大而逐渐降低，使用短基线组合测量精度低于长基线组合测量精度。本文的仿真方法，可为天线设计和工程应用时分析测向误差提供参考。参考文献：1徐立，颜雨.机载无源测向技术研究 J.科学技术创新，2022（19）：188-192.2 崔杨.无线电监测测向与定位技术研究 J.信息记录材料，2022，23（3）：169-171.3姜志成.干涉仪阵列测向误差分析及校准方法研究 J.电子世界，2021（20）：104-106.4吴启星，刘永旭，余银，等.基于阵元相位差模型的任意基线干涉仪定位误差分析算法J.电子信息对抗技术，2022，37（4）：65-69.5段国文，顾俊杰.干涉仪测向天线阵布阵与测向精度关系分析 J.电子信息对抗技术，2019，34（1）：56-60.6党浩淮，陈役涛，周国永.天线相位误差对干涉仪测向影响研究 J.航天电子对抗，2022，38（5）：41-45.7薛昭艺，刘宇，王章，等.一种基于最佳基线组的均匀圆阵测向算法 C 第十七届全国电波传播年会会议论文集，2022：376-379.8谢远举.非均匀阵列测向技术研究与实现 D.重庆：重庆大学，2017.9贾博伦.相关干涉仪测向天线系统研究 D.南京：南京航空航天大学，2018.图 11圆阵测量精度分布图（上接第 4页）4 结束语通过俄乌冲突能看出电磁频谱域的重要地位，必须有效控制和应用电磁频谱资源，以确保作战部队能最大程度地利用电磁域来支援部队的联合作战。虽然俄电子装备在近几次冲突中经受过实战考验，但依旧存在相关理论理念不足、电磁频谱应用管理不精和体系作战能力不强等问题。本文通过研究俄乌冲突中双方用频控频的方法和对抗博弈措施，验证了电磁频谱优势是制胜的关键因素。要从体系顶层设计层面和战术运用层面努力建设好具有自身特色的电磁频谱战装备体系，为未来制胜信息化高端战争打好基础。参考文献：1武坦然，易楷翔，皇甫冬琦，等.俄乌冲突中北约电磁频谱作战介入与俄军教训 J.航天电子对抗，2022，38（3）：1-4.2 许登荣，翁呈祥，杨慧君，等.从近几次局部冲突看俄军电子战能力的发展 J.航天电子对抗，2022，38（2）：59-64.3王全平，刘欣，衣春轮，等.国外新型低轨卫星星座发展及其 潜 在 军 事 应 用 研 究J.战 术 导 弹 技 术，2021（3）：67-74.4姜福涛，赵禄达.美军电磁频谱战发展及现状 J.航天电子对抗，2021，37（4）：60-64.-20