CN111693936B - 一种基于射频模拟接收系统的测向方法及系统 - Google Patents

一种基于射频模拟接收系统的测向方法及系统 Download PDF

Info

Publication number
CN111693936B
CN111693936B CN202010566250.7A CN202010566250A CN111693936B CN 111693936 B CN111693936 B CN 111693936B CN 202010566250 A CN202010566250 A CN 202010566250A CN 111693936 B CN111693936 B CN 111693936B
Authority
CN
China
Prior art keywords
phase
antenna
path
phase difference
angle
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202010566250.7A
Other languages
English (en)
Other versions
CN111693936A (zh
Inventor
左乐
聂剑坤
应钱诚
张浩斌
唐勇
王秀君
胡泽华
范保华
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
CETC 29 Research Institute
Original Assignee
CETC 29 Research Institute
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by CETC 29 Research Institute filed Critical CETC 29 Research Institute
Priority to CN202010566250.7A priority Critical patent/CN111693936B/zh
Publication of CN111693936A publication Critical patent/CN111693936A/zh
Application granted granted Critical
Publication of CN111693936B publication Critical patent/CN111693936B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
    • G01S3/02—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
    • G01S3/14—Systems for determining direction or deviation from predetermined direction
    • G01S3/46—Systems for determining direction or deviation from predetermined direction using antennas spaced apart and measuring phase or time difference between signals therefrom, i.e. path-difference systems
    • G01S3/48—Systems for determining direction or deviation from predetermined direction using antennas spaced apart and measuring phase or time difference between signals therefrom, i.e. path-difference systems the waves arriving at the antennas being continuous or intermittent and the phase difference of signals derived therefrom being measured
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
    • G01S3/02—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
    • G01S3/04—Details
    • G01S3/043—Receivers

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明提供了一种基于射频模拟接收系统的测向方法，对于N元均匀圆阵天线，通过鉴相器获得每一路天线与相邻天线的相位差，并记为该路鉴相器的相位差；在来波以仰角为0°入射圆阵天线时，记录每一路鉴相器的相位差为初始相位差；构造每一路鉴相器对应的复数；根据仰角
Figure DDA0002547934970000013
和方位角
Figure DDA0002547934970000011
，k为迭代次数，计算每路鉴相器的理论相位差；引入向量X和方向向量G，结合构造的复数计算迭代次数k+1的二维入射角的仰角
Figure DDA0002547934970000014
和方位角
Figure DDA0002547934970000012
；计算迭代次数为k+1与迭代次数为k的二维入射角之差，并于门限值比较，小于门限值进入输出仰角和方位角；否则，k的值加1，计算理论相位差。本方法优势在于避免了逻辑判断等运算、相位模糊、鉴相器输入信号幅度不均衡的问题，提高了测向精度。

Description

一种基于射频模拟接收系统的测向方法及系统
技术领域
本发明涉及无线电监测技术领域，特别涉及一种基于射频模拟接收系统的测向方法及系统。
背景技术
干涉仪测向体制属于相位测向体制，利用来波到达测向天线阵时，不同天线空间位置不同导致各天线单元接收的信号相位不同来求解来波的方向信息。干涉仪测向技术因其精度高、原理清晰、频带宽等特点在无源测向领域得到了广泛的应用。
相位干涉仪由于具有很高的测向精度，是目前测向系统中普遍采用的测向方法。线性最小二乘法由于运算简单，采用最小二乘对圆阵的相位信息进行处理，可得二维入射角(见文献：Y.Wu,H.C.So.Simple and accurate two-dimensional angle estimationfor a single source with uniform circular array[J].IEEE Antennas Wireless andPropagation Letters,2008,7(1):78–80；B.Liao,Y.T.Wu,S.C.Chan.A generalizedalgorithm for fast two-dimensional angle estimation of a single source withuniform circular arrays[J].IEEE Antennas Wireless Propagation Letters,2012,11(1):984–986)。但基于相位的最小二乘法未考虑相位模糊问题，不适用于大口径圆阵。
对于大口径圆阵，由于鉴相器的鉴相范围只有2π，当天线间距足够大而导致实际相位差超过2π时，鉴相器输出的相位值会以2π为周期翻转，出现多值模糊，这就造成真实相位与实际接收相位之间是多对一映射。当圆阵半径时，相位值域即会超过2π，从而出现相位模糊。要获得高的测角精度，两天线间的距离要越长；要实现相位不超过2π，两天线间的距离就要足够小，这与提高测向精度的条件相矛盾。
目前已报道的文献报道采用了多基线法(见文献：夏韶俊,杨晶.一维搜索与长短基线相结合的干涉仪设计方法[J].火控雷达技术,2013,42(02):33–37；季晓光,高晓光.一种机载无源定位方法——干涉仪定位[J].火力与指挥控制,2008,33(11):158–161；张刚兵,刘渝,刘宗敏.基线比值法相位解模糊算法[J].南京航空航天大学学报,2008,(05):665–669；李兴华,顾尔顺.干涉仪解模糊技术研究[J].现代防御技术,2008,(03):92-96；姜超.基于干涉仪测向的机载单站无源定位系统研究与应用[D].导师：陆起涌；赵波.复旦大学,2008；谢立允,王广松,戴旭初.圆阵相位干涉仪二维测向解模糊新方法[J].遥测遥控,2007,(05):53–59；赵国伟.高精度机载单站无源定位技术研究[D].导师：李勇.西北工业大学,2007；袁孝康.相位干涉仪测向定位研究[J].上海航天,1999,(03):3–9)、平行基线法(见文献：张艺帆.弹载雷达测向技术研究[D].导师：江朝抒.电子科技大学,2018；蒲刚.平行基线解模糊干涉仪测向算法及实现[D].导师：韩春林.电子科技大学,2013；电子科技大学.一种基于平行基线的圆阵相位干涉仪二维测向方法:中国,CN201110235023.7[P].2012-04-18；中国航天科技集团公司第五研究院第五一三研究所.一种组合相位差测向和空间谱测向的测向定位系统:中国,CN201510182536.4[P].2015-07-01)、扩展基线法(见文献：电子科技大学.一种扩展基线解模糊的相位干涉仪测向方法:中国,CN201110226585.5[P].2012-04-11)、相关法(见文献：谈文韬,林明,黎仁刚.奇偶阵元数均匀圆阵测向性能研究[J].现代雷达,2016,38(11):24–29；桂新涛.基于均匀圆阵的二维测向算法研究[D].导师：何子述.电子科技大学,2014；刘满超.无线电测向方法研究[D].导师：杨建红.兰州大学,2013；韩广,王斌,王成.相关运算在相位干涉仪解模糊中的应用[J].声学技术,2010,29(05):538–542；电子科技大学.一种基于相位差增量的相关干涉仪测向方法:中国,CN201310116050.1[P].2013-08-07)、虚拟基线法(见文献：蒋林鸿.基于GPU的宽带干涉仪信号处理及测向算法研究[D].导师：何子述.电子科技大学,2012；吴奉微,程婷,贾可新,何子述.基于虚拟阵列变换的干涉仪测向算法[J].现代雷达,2012,34(3):42–45；Li Peng-Fei.Broadband direction of arrival estimation based on virtual baselinetransformation and RBFNN[J].Journal of Astronautics,2012,33(2):210-216)等。还有最优基线法，根据方位角的不同，选取不同的最优基线组，求解二维入射角(见文献：潘玉剑；张晓发；黄敬健；杨骏；袁乃昌.模拟鉴相圆阵干涉仪测向性能的提高及其验证[J].系统工程与电子技术,37(6),2015)，存在计算复杂以及采用天线应用不充分的问题。
对高精度测向而言，要求天线间的距离足够大，这就会造成相位的模糊区间很大，一个相位差对应很多种可能的入射角，这就要求增加多组天线去除入射角的多值。而对于上述多基线、平行基线、扩展基线、虚拟基线等解模糊方法，精度越高，天线数越多，而角度测量仅使用了部分天线，其余的天线只是为了排除多值，导致硬件资源的不充分利用，无法通过增加的天线提高测角精度。
模拟体制接收机在干涉仪测向中广泛使用，其采用模拟鉴相器测量两路天线接收的相位差(见文献：J.H.Lee,J.M.Woo.Interferometer direction-finding system withimproved DF accuracy using two different array configurations[J].IEEEAntennas and Wireless Propagation Letters,2015,14:719–722)。但其采用所有天线对某一路天线鉴相，在电路实现上要求该路天线后接一分多的功分器，造成该路信号强度小、与其余路接收电平幅度不均衡的问题。
总之，以上方法存在以下问题：
(a)步骤复杂，需搜索、聚类门限设置、逻辑判断等步骤；
(b)计算量大；
(c)天线应用不充分导致精度不高；
(d)输入鉴相器的信号强度不均衡。
发明内容
本发明所要解决的技术问题是：针对上述存在的问题，提供了一种基于射频模拟接收系统的测向方法及系统，达到高精度实时二维测向的目的。
本发明采用的技术方案如下：一种基于射频模拟接收系统的测向方法包括以下过程：
步骤1、对于N元均匀圆阵天线，通过射频模拟接收系统的鉴相器获得每一路天线与相邻一路天线的电压相位差记为该路鉴相器的电压相位差，其中最后一路鉴相器的电压差为第一路天线与最后一路天线的电压相位差；
步骤2、在来波以仰角为0°入射圆阵天线时，记每一路鉴相器的相位差为初始相位差；
步骤3、根据电压相位差和初始相位差构造每一路鉴相器对应的复数；
步骤4、设初始仰角为
Figure BDA0002547934950000031
初始方位角为
Figure BDA0002547934950000032
迭代次数k为0；
步骤5、迭代次数为k时，根据仰角
Figure BDA0002547934950000033
和方位角
Figure BDA0002547934950000034
计算每路鉴相器的理论相位差；
步骤6、引入向量X和方向向量G结合构造的复数计算迭代次数k+1的二维入射角的仰角
Figure BDA0002547934950000035
和方位角
Figure BDA0002547934950000036
步骤7、计算迭代次数为k+1与迭代次数为k的二维入射角之差，并将二维入射角之差结果与门限值比较，若小于门限值，则进入步骤8；否则，迭代次数k的值加1，进入步骤5；
步骤8、输出二维入射角
Figure BDA0002547934950000037
和
Figure BDA0002547934950000038
作为辐射源的仰角和方位角。
进一步的，所述步骤1中，N元均匀圆阵天线中N≥3，第i路鉴相器的电压相位差输出为：
Figure BDA0002547934950000039
其中，Φi表示第i路天线的电压相位。
进一步的，所述步骤3中，构造复数的具体方法为：
Figure BDA00025479349500000310
其中，Vi表示第i路鉴相器对应的复数，ΔΦi表示第i路鉴相器的电压相位差，
Figure BDA00025479349500000311
表示来波以仰角为0°入射圆阵天线时，第i路鉴相器输出的电压相位差；i＝1,2,3,…,N。
进一步的，所述步骤5中，计算第i路鉴相器的理论相位差的具体方法为：
Figure BDA0002547934950000041
其中，
Figure BDA0002547934950000042
表示第i路鉴相器的理论相位差，c为空间中波速，ρ表示圆阵天线的半径，αi表示第i路天线的角度位置，其中，第一路天线的角度位置为0；两天线角度间距
Figure BDA0002547934950000043
进一步的，所述步骤6中，计算二维入射角的仰角
Figure BDA0002547934950000044
和方位角
Figure BDA0002547934950000045
的具体步骤为：
步骤61、设置向量
Figure BDA0002547934950000046
Figure BDA0002547934950000047
其中，T表示矩阵转置；
步骤62、计算方向向量，G＝[g1,g2]T；其中：
Figure BDA0002547934950000048
Figure BDA0002547934950000049
步骤63、更新向量X(k)，得到第k+1次向量X(k+1)的值；
X(k+1)＝X(k)+G
其中，
Figure BDA00025479349500000410
步骤64、根据X(k+1)计算第k+1次的二维入射角的仰角
Figure BDA00025479349500000411
和方位角
Figure BDA00025479349500000412
Figure BDA00025479349500000413
Figure BDA00025479349500000414
其中，||表示复数模值，real表示取复数实部运算，arg表示求复数相位运算，其中，mod(x,2π)为x对2π取模后的余数。
进一步的，所述步骤7中，计算迭代次数为k+1与迭代次数为k的二维入射角之差的具体方法为：
Figure BDA00025479349500000415
进一步的，所述步骤7中，门限值为所需测向精度。
进一步的，所述射频模拟接收系统包括：N路射频前端、N路一分二功分器、N路鉴相器，每一路天线依次连接一路射频前端、一路一分二功分器，天线接收的信号经射频前端放大、滤波后，进入一分二功分器，分为两路信号，所述鉴相器分别接收相邻两路天线对应的一分二功分器输出的信号，输出相邻两路的电压相位差，第一路一分二功分器和最后一路一分二功分器的其中一路信号连接同一路鉴相器。
本发明还提供了一种基于射频模拟接收系统的测向系统，包括：N元均匀圆阵天线、N路一分二功分器、N路鉴相器及信号处理模块，每一路天线依次连接一路射频前端、一路一分二功分器，天线接收的信号经射频前端放大、滤波后，进入一分二功分器，分为两路信号，所述鉴相器分别接收相邻两路天线对应的一分二功分器输出的信号，输出相邻两路天线的电压相位差，第一路一分二功分器和最后一路一分二功分器输出的一路信号连接同一路鉴相器；所述信号处理模块接收鉴相器输出的电压相位差，用于执行上述测向方法，输出辐射源的仰角和方位角。
与现有技术相比，采用上述技术方案的有益效果为：采用复数运算，从而避免了相位模糊问题；利用的圆阵上全部天线的信息，提高了测向精度高；采用迭代算法，避免了逻辑判断等运算；相邻两天线间采集相位差，避免了鉴相器输入信号幅度不均衡的问题。
附图说明
图1是本发明的测向系统架构图。
图2是本发明的均匀圆阵示意图。
图3是本发明的测向坐标系阵示意图。
图4是本发明的仰角均方根误差与入射仰角关系图。
图5是本发明的方位角均方根误差与入射仰角关系图。
图6是本发明的仰角均方根误差与入射方位角关系图。
图7是本发明的方位角均方根误差与入射方位角关系图。
具体实施方式
下面结合附图对本发明做进一步描述。
本发明提供了一种实现简单、精度高的二维入射角测量方法，采用如图1所示射频模拟接收系统，其主要功能是将空间电磁波转换为相位差。天线接收信号后经射频前端放大、滤波后，进入一分二功分器，分路后采用两两鉴相，第2路天线与第1路天线鉴相，记相位差为ΔΦ1；第3路天线与第2路天线鉴相，记相位差为ΔΦ2；以此类推至第N路天线与第N-1路天线鉴相，记相位差为ΔΦN-1。第1路天线与第N路天线鉴相，记相位差为ΔΦN。结合图1所示接线关系，相位差输出为
Figure BDA0002547934950000061
Φi表示第i路天线的电压相位。所采用的测向系统可采用、但不仅限于图1的形式。图中N为天线数量。
对一个N元均匀圆阵，N≥3，半径为ρ，第i个天线的角度位置为αi，如图2所示。两天线角度间距Δα＝2π/N。不失一般性，令第1个天线角度位置为0。不失一般性，天线阵位于xoy面，仰角θs为来波方向与z轴夹角，方位角为
Figure BDA0002547934950000062
来波方向与x轴夹角，如图3所示。具体测向方法如下：
一种基于射频模拟接收系统的测向方法包括以下过程：
步骤1、对于N元均匀圆阵天线，如图2所示，通过射频模拟接收系统的鉴相器获得每一路天线与后一路天线的电压相位差记为该路鉴相器的电压相位差，其中最后一路鉴相器的电压差为第一路天线与最后一路天线的电压相位差；
步骤2、在来波以仰角为0°入射圆阵天线时，记每一路鉴相器的相位差为初始相位差
Figure BDA0002547934950000063
步骤3、根据电压相位差和初始相位差构造每一路鉴相器对应的复数；
步骤4、设初始仰角为
Figure BDA0002547934950000064
初始方位角为
Figure BDA0002547934950000065
迭代次数k为0；
步骤5、根据仰角
Figure BDA0002547934950000066
和方位角
Figure BDA0002547934950000067
计算每路鉴相器的理论相位差；
步骤6、引入向量X和方向向量G结合构造的复数计算迭代次数k+1的二维入射角的仰角
Figure BDA0002547934950000068
和方位角
Figure BDA0002547934950000069
步骤7、计算迭代次数为k+1与迭代次数为k的二维入射角之差，并将二维入射角之差结果与门限值比较，若小于门限值，则进入步骤8；否则，迭代次数k的值加1，进入步骤5；
步骤8、输出二维入射角
Figure BDA00025479349500000610
和
Figure BDA00025479349500000611
作为辐射源的仰角和方位角。
在一个优选实施例中，所述步骤3中，构造复数的具体方法为：
Figure BDA00025479349500000612
其中，Vi表示第i路鉴相器对应的复数，ΔΦi表示第i路鉴相器的电压相位差，
Figure BDA00025479349500000613
表示来波以仰角为0°入射圆阵天线时，第i路鉴相器的相位差；i＝1,2,3,…,N。
在一个优选实施例中，所述步骤5中，计算每路鉴相器的理论相位差的具体方法为：
Figure BDA0002547934950000071
其中，
Figure BDA0002547934950000072
表示第i路鉴相器的理论相位差，c为空间中波速，ρ表示圆阵天线的半径，如图3所示，αi表示第i路天线的角度位置，其中，第一路天线的角度位置为0；两天线角度间距
Figure BDA0002547934950000073
在一个优选实施例中，所述步骤6中，计算二维入射角的仰角
Figure BDA0002547934950000074
和方位角
Figure BDA0002547934950000075
的具体步骤为：
步骤61、设置向量
Figure BDA0002547934950000076
Figure BDA0002547934950000077
其中，T表示矩阵转置；
步骤62、计算方向向量，G＝[g1,g2]T；其中：
Figure BDA0002547934950000078
Figure BDA0002547934950000079
步骤63、更新向量X(k)，得到第k+1次向量X(k+1)的值；
X(k+1)＝X(k)+G
其中，
Figure BDA00025479349500000710
步骤64、根据X(k+1)计算第k+1次的二维入射角的仰角
Figure BDA00025479349500000711
和方位角
Figure BDA00025479349500000712
Figure BDA00025479349500000713
Figure BDA00025479349500000714
其中，||表示复数模值，real表示取复数实部运算，arg表示求复数相位运算，其中，mod(x,2π)为x对2π取模后的余数。
在一个优选实施例中，所述步骤7中，计算迭代次数为k+1与迭代次数为k的二维入射角之差的具体方法为：
Figure BDA0002547934950000081
在一个优选实施例中，所述步骤7中，门限值为所需测向精度。
在一个优选实施例中，所述射频模拟接收系统包括：N路射频前端、N路一分二功分器、N路鉴相器，每一路天线依次连接一路射频前端、一路一分二功分器，天线接收的信号经射频前端放大、滤波后，进入一分二功分器，分为两路信号，所述鉴相器分别接收相邻两路天线对应的一分二功分器输出的信号，输出相邻两路的电压相位差，第一路一分二功分器和最后一路一分二功分器的其中一路信号连接同一路鉴相器。
本发明还提供了一种基于射频模拟接收系统的测向系统，包括：N元均匀圆阵天线、N路一分二功分器、N路鉴相器及信号处理模块，每一路天线依次连接一路射频前端、一路一分二功分器，天线接收的信号经射频前端放大、滤波后，进入一分二功分器，分为两路信号，所述鉴相器分别接收相邻两路天线对应的一分二功分器输出的信号，输出相邻两路天线的电压相位差，第一路一分二功分器和最后一路一分二功分器输出的一路信号连接同一路鉴相器；所述信号处理模块接收鉴相器输出的电压相位差，用于执行上述测向方法，输出辐射源的仰角和方位角。
本发明的方法的优势在于采用复数运算，从而避免了相位模糊问题；利用的圆阵上全部天线的信息，提高了测向精度高；采用迭代算法，避免了逻辑判断等运算；相邻两天线间采集相位差，避免了鉴相器输入信号幅度不均衡的问题。
为了验证本发明的效果，通过1000次蒙特卡洛仿真，计算测向误差的均方根，结果显示于图4至图7中。工作频率3GHz，信号噪声为高斯白噪声。
图4和图5分别为在不同相位噪声条件下、仰角均方根误差和方位角均方根误差与入射仰角的关系。仿真条件为：阵元数N＝11，组阵半径r＝360mm，相位噪声均方根分别取5°、10°、15°、20°，入射方位角60°。
图6和图7分别为在不同相位噪声条件下、仰角均方根误差和方位角均方根误差与入射方位角的关系。仿真条件为：阵元数N＝7，组阵半径r＝400mm，相位噪声均方根分别取5°、10°、15°、20°，入射仰角30°。
本发明并不局限于前述的具体实施方式。本发明扩展到任何在本说明书中披露的新特征或任何新的组合，以及披露的任一新的方法或过程的步骤或任何新的组合。如果本领域技术人员，在不脱离本发明的精神所做的非实质性改变或改进，都应该属于本发明权利要求保护的范围。
本说明书中公开的所有特征，或公开的所有方法或过程中的步骤，除了互相排斥的特征和/或步骤以外，均可以以任何方式组合。
本说明书中公开的任一特征，除非特别叙述，均可被其他等效或具有类似目的的替代特征加以替换。即，除非特别叙述，每个特征只是一系列等效或类似特征中的一个例子而已。

Claims (8)

1.一种基于射频模拟接收系统的测向方法，其特征在于，包括以下过程：
步骤1、对于N元均匀圆阵天线，通过射频模拟接收系统的鉴相器获得每一路天线与相邻天线的电压相位差，并记为该路鉴相器的电压相位差，其中最后一路鉴相器的电压相位差为第一路天线与最后一路天线的电压相位差；
步骤2、在来波以仰角为0°入射圆阵天线时，记每一路鉴相器的电压相位差为初始相位差
Figure FDA0003725800190000011
步骤3、根据电压相位差和初始相位差构造每一路鉴相器对应的复数；
步骤4、设初始仰角为
Figure FDA0003725800190000012
初始方位角为
Figure FDA0003725800190000013
迭代次数k为0；
步骤5、迭代次数为k时，根据仰角
Figure FDA0003725800190000014
和方位角
Figure FDA0003725800190000015
计算每路鉴相器的理论相位差；
步骤6、引入向量X和方向向量G，结合构造的复数计算迭代次数k+1的二维入射角的仰角
Figure FDA0003725800190000016
和方位角
Figure FDA0003725800190000017
步骤7、计算迭代次数为k+1与迭代次数为k的二维入射角之差，并将二维入射角之差结果与门限值比较，若小于门限值，则进入步骤8；否则，迭代次数k的值加1，进入步骤5；
步骤8、输出二维入射角
Figure FDA0003725800190000018
和
Figure FDA0003725800190000019
作为辐射源的仰角和方位角；
计算二维入射角的仰角
Figure FDA00037258001900000110
和方位角
Figure FDA00037258001900000111
的具体步骤为：
步骤61、设置向量
Figure FDA00037258001900000112
Figure FDA00037258001900000113
其中，T表示矩阵转置；
Figure FDA00037258001900000114
表示两天线角度间距；
步骤62、计算方向向量，G＝[g1,g2]T；其中：
Figure FDA00037258001900000115
Figure FDA00037258001900000116
其中，f均表示信号频率，c为空间中波速，ρ表示圆阵天线的半径，αi表示第i路天线的角度位置，Vi表示第i路鉴相器对应的复数，
Figure FDA00037258001900000117
表示第i路鉴相器的理论相位差；
步骤63、更新向量X(k)，得到第k+1次向量X(k+1)的值；
X(k+1)＝X(k)+G
其中，
Figure FDA0003725800190000021
步骤64、根据X(k+1)计算第k+1次的二维入射角的仰角
Figure FDA0003725800190000022
和方位角
Figure FDA0003725800190000023
Figure FDA0003725800190000024
Figure FDA0003725800190000025
其中，||表示复数模值，real表示取复数实部运算，arg表示求复数相位运算，mod(x,2π)为x对2π取模后的余数。
2.根据权利要求1所述的基于射频模拟接收系统的测向方法，其特征在于，所述步骤1中，N元均匀圆阵天线中N≥3，第i路鉴相器的电压相位差输出为：
Figure FDA0003725800190000026
其中，Φi表示第i路天线的电压相位。
3.根据权利要求2所述的基于射频模拟接收系统的测向方法，其特征在于，所述步骤3中，构造复数的具体方法为：
Figure FDA0003725800190000027
其中，Vi表示第i路鉴相器对应的复数，ΔΦi表示第i个鉴相器的电压相位差，
Figure FDA0003725800190000028
表示来波以仰角为0°入射圆阵天线时，第i个鉴相器输出的电压相位差；i＝1,2,3,…,N。
4.根据权利要求3所述的基于射频模拟接收系统的测向方法，其特征在于，所述步骤5中，计算每路鉴相器的理论相位差的具体方法为：
Figure FDA0003725800190000029
其中，
Figure FDA00037258001900000210
表示第i路鉴相器的理论相位差，c为空间中波速，ρ表示圆阵天线的半径，f均表示信号频率，αi表示第i路天线的角度位置，其中，第一路天线的角度位置为0；两天线角度间距
Figure FDA00037258001900000211
5.根据权利要求1所述的基于射频模拟接收系统的测向方法，其特征在于，所述步骤7中，计算迭代次数为k+1与迭代次数为k的二维入射角之差的具体方法为：
Figure FDA00037258001900000212
6.根据权利要求1所述的基于射频模拟接收系统的测向方法，其特征在于，所述步骤7中，门限值为所需测向精度。
7.根据权利要求1所述的基于射频模拟接收系统的测向方法，其特征在于，所述射频模拟接收系统包括：N路射频前端、N路一分二功分器、N路鉴相器，每一路天线依次连接一路射频前端、一路一分二功分器，天线接收的信号经射频前端放大、滤波后，进入一分二功分器，分为两路信号，所述鉴相器分别接收相邻两路天线对应的一分二功分器输出的信号，输出相邻两路的电压相位差，第一路一分二功分器和最后一路一分二功分器的其中一路信号连接同一路鉴相器。
8.一种基于射频模拟接收系统的测向系统，其特征在于，包括：N元均匀圆阵天线、N路一分二功分器、N路鉴相器及信号处理模块，每一路天线依次连接一路射频前端、一路一分二功分器，天线接收的信号经射频前端放大、滤波后，进入一分二功分器，分为两路信号，所述鉴相器分别接收相邻两路天线对应的一分二功分器输出的信号，输出相邻两路天线的电压相位差，第一路一分二功分器和最后一路一分二功分器输出的一路信号连接同一路鉴相器；所述信号处理模块接收鉴相器输出的电压相位差，用于执行权利要求1-6任一项的测向方法，输出辐射源的仰角和方位角。
CN202010566250.7A 2020-06-19 2020-06-19 一种基于射频模拟接收系统的测向方法及系统 Active CN111693936B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010566250.7A CN111693936B (zh) 2020-06-19 2020-06-19 一种基于射频模拟接收系统的测向方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010566250.7A CN111693936B (zh) 2020-06-19 2020-06-19 一种基于射频模拟接收系统的测向方法及系统

Publications (2)

Publication Number Publication Date
CN111693936A CN111693936A (zh) 2020-09-22
CN111693936B true CN111693936B (zh) 2022-08-09

Family

ID=72482213

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010566250.7A Active CN111693936B (zh) 2020-06-19 2020-06-19 一种基于射频模拟接收系统的测向方法及系统

Country Status (1)

Country Link
CN (1) CN111693936B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111693935B (zh) * 2020-06-19 2023-04-21 中国电子科技集团公司第二十九研究所 一种射频测向方法及系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5463397A (en) * 1993-10-25 1995-10-31 Hughes Aircraft Company Hyper-precision SAR interferometry using a dual-antenna multi-pass SAR system
CN102175988A (zh) * 2011-01-21 2011-09-07 电子科技大学 一种基于维度拆分的相关干涉仪测向方法
CN102419432A (zh) * 2011-08-25 2012-04-18 电子科技大学 一种基于虚拟基线的圆阵相位干涉仪二维测向方法

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5030959A (en) * 1986-06-03 1991-07-09 Southwest Research Institute Enhanced automatic bearing indicator system
US6421008B1 (en) * 2001-02-12 2002-07-16 The Aerospace Corporation Method to resolve interferometric ambiguities
US7201054B2 (en) * 2004-10-18 2007-04-10 The Boeing Company System and method for resolving phase ambiguity of a transducer array to determine direction of arrival of received signals
CN102411136B (zh) * 2011-08-09 2013-05-08 电子科技大学 一种扩展基线解模糊的相位干涉仪测向方法
DE102012208153A1 (de) * 2012-05-15 2013-11-21 Plath Gmbh Verfahren und Vorrichtung zum Bestimmen des Elevationswinkels und des Azimuts bei Interferometerpeilern
US10228443B2 (en) * 2012-12-02 2019-03-12 Khalifa University of Science and Technology Method and system for measuring direction of arrival of wireless signal using circular array displacement
CN104122527B (zh) * 2014-07-14 2016-08-17 中国人民解放军国防科学技术大学 一种基于查表法的圆阵相位干涉仪宽带瞬时测向方法
CN109164408B (zh) * 2018-07-17 2022-09-13 中国电子科技集团公司第二十九研究所 一种采用两传感器的频变信号二维测向方法及设备
CN109270486B (zh) * 2018-07-25 2020-09-22 中国电子科技集团公司第二十九研究所 一种采用两传感器旋转的频变信号二维测向方法
CN109633526B (zh) * 2019-01-29 2020-09-01 杭州电子科技大学 基于方向函数的非均匀圆阵相位干涉仪测向解模糊方法
CN110007267B (zh) * 2019-01-29 2020-08-18 杭州电子科技大学 一种基于混合基线的均匀圆阵干涉仪测向解模糊方法
CN109901104A (zh) * 2019-04-18 2019-06-18 中国电子科技集团公司第二十九研究所 一种通过估计时差解干涉仪测向模糊的方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5463397A (en) * 1993-10-25 1995-10-31 Hughes Aircraft Company Hyper-precision SAR interferometry using a dual-antenna multi-pass SAR system
CN102175988A (zh) * 2011-01-21 2011-09-07 电子科技大学 一种基于维度拆分的相关干涉仪测向方法
CN102419432A (zh) * 2011-08-25 2012-04-18 电子科技大学 一种基于虚拟基线的圆阵相位干涉仪二维测向方法

Also Published As

Publication number Publication date
CN111693936A (zh) 2020-09-22

Similar Documents

Publication Publication Date Title
CN111693934B (zh) 一种基于数字干涉系统的测向方法及系统
CN102411136B (zh) 一种扩展基线解模糊的相位干涉仪测向方法
CN109164408B (zh) 一种采用两传感器的频变信号二维测向方法及设备
CN109270486B (zh) 一种采用两传感器旋转的频变信号二维测向方法
CN103344940B (zh) 低复杂度的doa估计方法及系统
CN106227701B (zh) 一种阵列信号的接收通道的幅相误差的自校正方法
CN108594164B (zh) 一种平面阵列doa估计方法及设备
CN107576951B (zh) 基于嵌套式电磁矢量传感器阵列的波达方向估计方法
CN102830386A (zh) 一种强干扰背景下的任意阵型微弱信号源角度估计方法
CN104931923A (zh) Grid Iterative ESPRIT,一种可扩展的用于均匀圆阵二维到达角的快速估计算法
CN111693935B (zh) 一种射频测向方法及系统
CN108363048B (zh) 一种基于块稀疏的极化mimo雷达的角度估计方法
CN111693936B (zh) 一种基于射频模拟接收系统的测向方法及系统
CN111983554A (zh) 非均匀l阵下的高精度二维doa估计
CN103323810B (zh) 一种l阵方位角和俯仰角配对的信号处理方法
CN110907923A (zh) 基于平行因子算法的双基地emvs-mimo雷达角度估计算法及装置
CN113391261A (zh) 一种基于泰勒展开的增广互质阵列快速doa估计方法
CN112363108A (zh) 信号子空间加权超分辨的波达方向检测方法及系统
Zang et al. Research on hybrid algorithm based on TDOA
He et al. Extended aperture 2-D direction finding with a two-parallel-shape-array using propagator method
CN105334490A (zh) 一种宽带信号干涉测向的方法
Ning et al. A velocity independent MUSIC algorithm for DOA estimation
Wang et al. Estimating the DOA and polarization parameters with sparse collocated loop and dipole cross array
Wang et al. The DOA and polarization estimation of coherent signals via UCA-IEVS
CN113050029B (zh) 一种适用于阵元失效条件下的相位干涉仪测向方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant