CN113050029B - 一种适用于阵元失效条件下的相位干涉仪测向方法 - Google Patents
一种适用于阵元失效条件下的相位干涉仪测向方法 Download PDFInfo
- Publication number
- CN113050029B CN113050029B CN202110326629.5A CN202110326629A CN113050029B CN 113050029 B CN113050029 B CN 113050029B CN 202110326629 A CN202110326629 A CN 202110326629A CN 113050029 B CN113050029 B CN 113050029B
- Authority
- CN
- China
- Prior art keywords
- phase difference
- phi
- antenna
- base line
- baseline
- 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
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
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P10/00—Technologies related to metal processing
- Y02P10/25—Process efficiency
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
本发明提出了一种在阵元失效条件下干涉仪测向方法,旨在实现相位干涉仪在阵元失效情况下测向。实现步骤为:构建相位干涉仪测向系统;每个天线阵元接收入射信号并传输;信号处理机判断是否存在失效天线阵元;信号处理机根据阵元失效情况获取基线相位差和虚拟基线相位差;信号处理机对基线相位解模糊;信号处理机获取相位干涉仪的测向结果。本发明设计的非均匀L型天线阵,可以利用基线长度关系构造虚拟基线进行解模糊,实现在阵元失效的情况下测向,同时在阵元完好的情况也能保证测向精度,提高了相位干涉仪测向系统的鲁棒性。
Description
技术领域
本发明属于无源测向技术领域,特别涉及一种在阵元失效条件下的相位干涉仪测向方法。
背景技术
传统的有源雷达测向技术过发射大功率信号来实现,具有作用时间长、距离远、精度高等优点,并且受外界环境因素变化的影响小。但是随着电子以及通信信息技术的发展,有源测向技术在被广泛使用的同时,其系统自身隐藏与抗干扰能力差的缺点凸显,无法隐藏容易遭到干扰与攻击。为了弥补有源测向技术的不足与可能产生的严重后果,人们开始研究无源测向技术。
相位干涉仪测向技术即属于无源测向技术的一种,它本身不会辐射信源,只接收处理来自目标的位置及运动状态的相关数据来实现对目标的测向。当入射信号到达干涉仪测向天线阵时,由于天线的空间位置不同,每个天线单元接收到的信号相位也不同,干涉仪就是利用天线单元之间的相位差来获取来波的方向信息。干涉仪测向系统具有灵敏度高、测向精准、算法简单等优点,广泛地用于电磁环境监测、电子对抗等领域。
相位干涉仪主要利用来波在基线间产生的相位差来计算方位,影响其测向精度和效率主要取决于布阵的形式和测向算法。如果出现阵元失效的情况,则会破坏原有阵列的几何分布特征,丢失该通道的数据,导致整个系统的性能下降,甚至无法工作。特别是在机载等天线布置空间有限的平台上,允许天线阵元的尺寸和数量都非常有限,往往出现某个阵元失效就会导致整个设备无法使用。
目前相位干涉仪的相关研究都是基于阵元完好的情况,例如申请公布号CN111323742A,名称为“一种基于曲线阵列的相位干涉仪及其测向方法”的专利申请中,公开了一种基于曲线阵列的相位干涉仪及其测向方法,该方法通过保证一组3个天线的阵列的纵向间隔相等,在横坐标方向和纵坐标方向分别形成虚拟基线,抵消天线位置的纵坐标带来的影响,将曲线阵列模型等效为一维阵列模型进行测向。而干涉仪通常都在室外工作,由于恶劣的自然环境、人为干扰以及阵元老化、阵元的物理损坏等因素的影响,阵元失效问题不可避免,对应阵元失效情况,申请公布号CN111323742A所述的方法就无法保证选取的一组3个天线的阵列的纵向间隔相等,从而导致算法失效。
发明内容
本发明的目的在于克服上述现有技术存在的缺陷,提出一种适用于阵元失效条件下的相位干涉仪测向方法,用于解决现有技术中存在的出现阵元失效情况相位干涉仪无法测向的问题。
为实现上述目的,本发明采取的技术方案包括如下步骤:
(1)构建相位干涉仪测向系统:
构建包括天线阵和信号处理机的相位干涉仪测向系统,所述天线阵采用由正交的第一一维线阵天线I和第二一维线阵天线II组成的L型结构,所述第一一维线阵天线I包括非均匀排布的N个天线阵元(t1,t2,...,tn,...,tN),所述第二一维线阵天线II包括非均匀排布的N个天线阵元(t1,t2′,...,tn′,...,tN′),N≥4,第一天线阵元t1位于L型结构两个臂的交汇位置,为I和II共用,2N-1个天线阵元分别通过接收通道与信号处理机连接;t1与t2之间的基线长度d12以及t1与t2′之间的基线长度d12′、t2与t3之间的基线长度d23以及t2′与t3′之间的基线长度d23′、t3与t4之间的基线长度d34以及t3′与t4′之间的基线长度d34′满足d12=d34,d12′=d34′,且d23-d12≤λ/2,d23′-d12′≤λ/2,当N>4时,t4以及t4′以后的各天线阵元之间的基线长度满足dn,1=kdn-1,1,dn,1′=kdn-1,1′(n>4),其中tn表示I中的第n个天线阵元,tn′表示II中的第n个天线阵元,λ表示入射信号的波长,k为正整数;
(2)每个天线阵元接收入射信号并传输:
天线阵中的每个天线阵元通过接收通道将接收的入射信号传输至信号处理机,得到第一一维线阵天线I对应的入射信号集合V=(v1,v2,…,vn,…,vN)和第二一维线阵天线II对应的入射信号集合V′=(v1,v2′,…,vn′,…,vN′),其中:
其中,dn,1和dn,1′表示tn和tn′与t1的基线长度,An和An′分别表示tn和tn′接收的射频信号,θ和β分别表示入射信号的方位角和俯仰角,F(θ,β)表示天线阵元的方向图函数,e为自然常数,j为虚数符号;
(3)信号处理机判断是否存在失效天线阵元:
信号处理机判断入射信号集合V=(v1,v2,…,vn,…,vN)是否存在一个天线阵元的信号量va=0,a∈(1,2,...,N),同时判断入射信号集合V′=(v1,v2′,…,vn′,…,vN′)是否存在一个天线阵元的信号量va′=0,a∈(1,2,...,N),若是,则存在失效天线阵元,并执行步骤(4),否则,则不存在失效天线阵元,并执行步骤(5);
(4)信号处理机获取基线相位差和虚拟基线相位差:
信号处理机获取I的基线相位差(φ1,φ2,...,φm,...,φM)和虚拟基线相位差φ0以及II的基线相位差(φ1′,φ2′,...,φm′,...,φM′)和虚拟基线相位差φ0′,并执行步骤(6),具体步骤为:
(4a)信号处理机获取I的基线相位差(φ1,φ2,...,φm,...,φM)和虚拟基线相位差φ0:
当n=1时,va=v1,信号处理机对V的第二天线阵元t2与其余天线阵元(t3,t4,...,tn,...tN)中的每一个之间的基线相位差进行测量,得到M条基线相位差(φ1,φ2,...,φm,...,φM),同时对第三天线阵元t3与第四天线阵元t4之间的基线相位差φ34进行测量,并利用t2和t3构成的基线与t3和t4构成的基线构造长度为d0=d12-d34虚拟基线,再计算虚拟基线的相位差φ0=φ1-φ34,其中,φm表示第二天线阵元t2与第m+2天线阵元之间的基线相位差,M=N-2;
当n=2时,va=v2,信号处理机对V的第一天线阵元t1与其余天线阵元(t3,t4,...,tn,...tN)中的每一个之间的基线相位差进行测量,得到M条基线相位差(φ1,φ2,...,φm,...,φM),同时对第三天线阵元t3与第四天线阵元t4之间的基线相位差φ34进行测量,并利用t1和t3构成的基线与t3和t4构成的基线构造长度为d0=d13-2d34的虚拟基线,再计算虚拟基线的相位差φ0=φ1-2φ34,其中,d13表示第一天线阵元t1与第三天线阵元t3之间的距离,φm表示第一天线阵元t1与第m+2天线阵元之间的基线相位差;
当n=3时,va=v3,信号处理机对V的第一天线阵元t1与其余天线阵元(t2,t4,t5,...,tn,...tN)中的每一个之间的基线相位差进行测量,得到M条基线相位差(φ1,φ2,...,φm,...,φM),同时对第二天线阵元t2与第四天线阵元t4之间的基线相位差φ24进行测量,并利用t1和t2构成的基线与t2和t4构成的基线构造长度为d0=d24-2d12虚拟基线,再计算虚拟基线的相位差φ0=φ24-2φ1,其中,d24表示第二天线阵元t2与第四天线阵元t4之间的距离,当m=1时,φm表示第一天线阵元t1与第二天线阵元t2之间的基线相位差,当m>1时,φm表示第一天线阵元t1与第m+2天线阵元之间的基线相位差;
当n>3时,va=vn,信号处理机对V的第一天线阵元t1与其余天线阵元(t2,t3,...,tn-1,tn+1,...,tN)中的每一个之间的基线相位差进行测量,得到M条基线相位差(φ1,φ2,...,φm,...,φM),同时对第二天线阵元t2与第三天线阵元t3之间的基线相位差φ23进行测量,并利用t1和t2构成的基线与t2和t3构成的基线构造长度为d0=d23-d12虚拟基线,再计算虚拟基线的相位差φ0=φ23-φ1,其中,当m<n时,φm表示第一天线阵元t1与第m+1天线阵元之间的基线相位差,当m≥n时,φm表示第一天线阵元t1与第m+2天线阵元之间的基线相位差;
(4b)信号处理机获取II的基线相位差(φ1′,φ2′,...,φm′,...,φM′)和虚拟基线相位差φ0′:
信号处理机按照步骤(4a)的方法获取II的基线相位差(φ1′,φ2′,...,φm′,...,φM′)和虚拟基线相位差φ0′;
(5)信号处理机获取基线相位差和虚拟基线相位差:
信号处理机获取I的基线相位差(φ1,φ2,...,φm,...,φM)和虚拟基线相位差φ0以及II的基线相位差(φ1′,φ2′,...,φm′,...,φM′)和虚拟基线相位差φ0′,并执行步骤(6),具体步骤为:
(5a)信号处理机获取I的基线相位差(φ1,φ2,...,φm,...,φM)和虚拟基线相位差φ0:
信号处理机对V的第一天线阵元t1与其余天线阵元(t2,t3,...,tn,...tN)中的每一个之间的基线相位差进行测量,得到M条基线相位差(φ1,φ2,...,φm,...,φM),同时对第二天线阵元t2与第三天线阵元t3之间的基线相位差φ23进行测量,并利用t1和t2构成的基线与t2和t3构成的基线构造长度为d0=d23-d12虚拟基线,再计算虚拟基线的相位差φ0=φ23-φ1,其中,φm表示第一天线阵元t1与第m+1天线阵元之间的基线相位差,M=N-1;
(5b)信号处理机获取II的基线相位差(φ1′,φ2′,...,φm′,...,φM′)和虚拟基线相位差φ0′:
信号处理机对V′的第一天线阵元t1与其余天线阵元(t2′,t3′,...,tn′,...,tN′)中的每一个之间的基线相位差进行测量,得到M′条基线相位差(φ1′,φ2′,...,φm′,...,φM′),同时对第二天线阵元t2′与第三天线阵元t3′之间的基线相位差φ23′进行测量,并利用t1′和t2′构成的基线与t2′和t3′构成的基线构造长度为d0′=d23′-d12′虚拟基线,再计算虚拟基线的相位差φ0′=φ23′-φ1′,其中,φm′表示第一天线阵元t1与第m+1天线阵元之间的基线相位差,M′=N-1;
(6)信号处理机对I和II的基线相位差分别进行解相位模糊:
信号处理机通过虚拟基线无模糊相位差φ0对I的基线相位差(φ1,φ2,...,φm,...,φM)逐级进行解相位模糊,得到I的基线的无模糊相位差同时通过虚拟基线无模糊相位差φ0′对II的基线相位差(φ1′,φ2′,...,φm′,...,φM′)逐级进行解相位模糊,得到II的基线的无模糊相位差的计算公式为:
其中(φ1,φ2,...,φm,...,φM),(φ1′,φ2′,...,φm′,...,φM′)表示执行步骤(4)或步骤(5)得到的I和II的基线相位差,dm表示对应的基线长度,dm′表示φm′对应的基线长度,表示φm对应的无模糊相位差,表示φm′对应的无模糊相位差,M表示I的基线相位差个数,I存在一个阵元失效时,M=N-2,不存阵元失效时M=N-1,M′表示II的基线相位差个数,II存在一个阵元失效时,M′=N-2,不存阵元失效时M′=N-1;
(7)信号处理机获取相位干涉仪的测向结果:
其中S=min(M,M′),min表示取最小值函数。
本发明与现有技术相比,具有以下优点:
本发明采用的非均匀L型结构天线,基线长度满足:d12=d34,d12′=d34′,且d23-d12≤λ/2,d23′-d12′≤λ/2,当N>4时,t4以及t4′以后的各天线阵元之间的基线长度满足dn,1=kdn-1,1,dn,1′=kdn-1,1′(n>4),在此基线长度关系下,可通过构造虚拟基线进行解相位模糊,实现在阵元失效的情况下测向,同时在阵元完好的情况下也能保证测向精度,提高了干涉仪测向系统的鲁棒性。
附图说明
图1是本发明的实现流程图;
图2是本发明中构建的相位干涉仪测向系统中天线阵的结构示意图;
具体实施方式
下面结合附图和具体实施例,对本发明作进一步详细描述。
参照图1,本发明包括如下步骤:
步骤1)构建相位干涉仪测向系统:
构建相位干涉仪测向系统,包括天线阵和信号处理机的相位干涉仪测向系统,所述的天线阵结构如图2所示,采用由正交的第一一维线阵天线I和第二一维线阵天线II组成的L型结构,所述第一一维线阵天线I包括非均匀排布的N个天线阵元(t1,t2,...,tn,...,tN),所述第二一维线阵天线II包括非均匀排布的N个天线阵元(t1,t2′,...,tn′,...,tN′),N≥4,第一天线阵元t1位于L型结构两个臂的交汇位置,为I和II共用,2N-1个天线阵元分别通过接收通道与信号处理机连接;t1与t2之间的基线长度d12以及t1与t2′之间的基线长度d12′、t2与t3之间的基线长度d23以及t2′与t3′之间的基线长度d23′、t3与t4之间的基线长度d34以及t3′与t4′之间的基线长度d34′满足d12=d34,d12′=d34′,且d23-d12≤λ/2,d23′-d12′≤λ/2,当N>4时,t4以及t4′以后的各天线阵元之间的基线长度满足dn,1=kdn-1,1,dn,1′=kdn-1,1′(n>4),其中tn表示I中的第n个天线阵元,tn′表示II中的第n个天线阵元,λ表示入射信号的波长,k为正整数;
本实施例中,设有如下要求设计天线阵:覆盖频段0.8GHz~6GHz;天线圆盘直径≤300mm;采用平面螺旋天线,直径r=70mm,取N=4。系统的最高测向频率为6GHz,对应的λ=50mm,要获得无模糊相位,显然物理基线长度不可能做到50mm,需要使用虚拟基线的方式,则d23-d12≤25,d23′-d12′≤25,天线圆盘直径≤300mm,则d12+d23+d34≤300,d12′+d23′+d34′≤300,故可设d12=d12′=90mm,d23=d23′=115mm,d34=d34′=90mm。
步骤2)每个天线阵元接收入射信号并传输:
如图2所示,辐射源发射信号,天线阵接收到的入射信号方位角和俯仰角为θ和β,θ表示辐射源的阵面投影和t1构成直线与II之间的夹角,β表示辐射源和t1构成直线与阵面之间的夹角。
天线阵中的每个天线阵元通过接收通道将接收的入射信号传输至信号处理机,得到第一一维线阵天线I对应的入射信号集合V=(v1,v2,v3,v4)和第二一维线阵天线II对应的入射信号集合V′=(v1,v2′,v3′,v4′),其中:
其中,n∈(1,2,3,4),dn,1和dn,1′表示tn和tn′与t1的基线长度,An和An′分别表示tn和tn′接收的射频信号,F(θ,β)表示天线阵元的方向图函数,e为自然常数,j为虚数符号;
步骤3)信号处理机判断是否存在失效天线阵元:
信号处理机判断入射信号集合V=(v1,v2,v3,v4)是否存在一个天线阵元的信号量va=0,a∈(1,2,3,4),同时判断入射信号集合V′=(v1,v2′,v3′,v4′)是否存在一个天线阵元的信号量va′=0,a∈(1,2,3,4),若是,则存在失效天线阵元,并执行步骤(4),否则,则不存在失效天线阵元,并执行步骤(5);
步骤4)信号处理机获取基线相位差和虚拟基线相位差:
信号处理机获取I的基线相位差(φ1,φ2)和虚拟基线相位差φ0以及II的基线相位差(φ1′,φ2′)和虚拟基线相位差φ0′,并执行步骤(6),具体步骤为:
(4a)信号处理机获取I的基线相位差(φ1,φ2)和虚拟基线相位差φ0:
当va=v1,信号处理机对V的第二天线阵元t2与第三天线阵元t3之间的基线相位φ1、第二天线阵元t2与第三天线阵元t4之间的基线相位φ2、第三天线阵元t3与第四天线阵元t4之间的基线相位差φ34,并利用t2和t3构成的基线与t3和t4构成的基线构造长度为d0=d12-d34=25mm虚拟基线,再计算虚拟基线的相位差φ0=φ1-φ34;
当va=v2,信号处理机对V的第二天线阵元t1与第三天线阵元t3之间的基线相位φ1、第二天线阵元t1与第三天线阵元t4之间的基线相位φ2、第三天线阵元t3与第四天线阵元t4之间的基线相位差φ34进行测量,并利用t1和t3构成的基线与t3和t4构成的基线构造长度为d0=d13-2d34=25mm的虚拟基线,再计算虚拟基线的相位差φ0=φ1-2φ34,其中,d13表示第一天线阵元t1与第三天线阵元t3之间的距离,d13=d12+d23=205mm;
当va=v3,信号处理机对V的第二天线阵元t1与第三天线阵元t2之间的基线相位φ1、第二天线阵元t1与第三天线阵元t4之间的基线相位φ2、第二天线阵元t2与第四天线阵元t4之间的基线相位差φ24进行测量,并利用t1和t2构成的基线与t2和t4构成的基线构造长度为d0=d24-2d12=25mm虚拟基线,再计算虚拟基线的相位差φ0=φ24-2φ1,其中,d24表示第二天线阵元t2与第四天线阵元t4之间的距离,d24=d23+d34=205mm;
当va=v4,信号处理机对V的第二天线阵元t1与第三天线阵元t2之间的基线相位φ1、第二天线阵元t1与第三天线阵元t3之间的基线相位φ2、第二天线阵元t2与第三天线阵元t3之间的基线相位差φ23进行测量,并利用t1和t2构成的基线与t2和t3构成的基线构造长度为d0=d23-d12=25mm虚拟基线,再计算虚拟基线的相位差φ0=φ23-φ1;
(4b)信号处理机获取II的基线相位差(φ1′,φ2′)和虚拟基线相位差φ0′:
信号处理机按照步骤(4a)的方法获取II的基线相位差(φ1′,φ2′)和虚拟基线相位差φ0′;
步骤5)信号处理机获取基线相位差和虚拟基线相位差:
信号处理机获取I的基线相位差(φ1,φ2,φ3)和虚拟基线相位差φ0以及II的基线相位差(φ1′,φ2′,φ3′)和虚拟基线相位差φ0′,并执行步骤(6),具体步骤为:
(5a)信号处理机获取I的基线相位差(φ1,φ2,φ3)和虚拟基线相位差φ0:
信号处理机对V的第一天线阵元t1与其余天线阵元(t2,t3,t4)中的每一个之间的基线相位差进行测量,得到3条基线相位差(φ1,φ2,φ3),同时对第二天线阵元t2与第三天线阵元t3之间的基线相位差φ23进行测量,并利用t1和t2构成的基线与t2和t3构成的基线构造长度为d0=d23-d12=25mm虚拟基线,再计算虚拟基线的相位差φ0=φ23-φ1;
(5b)信号处理机获取II的基线相位差(φ1′,φ2′,φ3′)和虚拟基线相位差φ0′:
信号处理机对V′的第一天线阵元t1与其余天线阵元(t2′,t3′,t4′)中的每一个之间的基线相位差进行测量,得到3条基线相位差(φ1′,φ2′,φ3′),同时对第二天线阵元t2′与第三天线阵元t3′之间的基线相位差φ23′进行测量,并利用t1′和t2′构成的基线与t2′和t3′构成的基线构造长度为d0′=d23′-d12′=25mm虚拟基线,再计算虚拟基线的相位差φ0′=φ23′-φ1′;
步骤6)信号处理机对I和II的基线相位差分别进行解相位模糊:
信号处理机通过虚拟基线无模糊相位差φ0对I的基线相位差(φ1,φ2,...,φm,...,φM)逐级进行解相位模糊,得到I的基线的无模糊相位差同时通过虚拟基线无模糊相位差φ0′对II的基线相位差(φ1′,φ2′,...,φm′,...,φM′)逐级进行解相位模糊,得到II的基线的无模糊相位差的计算公式为:
其中dm表示对应的基线长度,dm′表示φm′对应的基线长度,表示φm对应的无模糊相位差,表示φm′对应的无模糊相位差,M表示I的基线相位差个数,I存在一个阵元失效时,M=2,不存阵元失效时M=3,M′表示II的基线相位差个数,II存在一个阵元失效时,M′=2,不存阵元失效时M′=3;
步骤7)信号处理机获取相位干涉仪的测向结果:
其中S=min(M,M′),min表示取最小值函数。
Claims (2)
1.一种适用于阵元失效条件下的相位干涉仪测向方法,其特征在于,包括如下步骤:
(1)构建相位干涉仪测向系统:
构建包括天线阵和信号处理机的相位干涉仪测向系统,所述天线阵采用由正交的第一一维线阵天线I和第二一维线阵天线II组成的L型结构,所述第一一维线阵天线I包括非均匀排布的N个天线阵元(t1,t2,...,tn,...,tN),所述第二一维线阵天线II包括非均匀排布的N个天线阵元(t1,t2′,...,tn′,...,tN′),N≥4,第一天线阵元t1位于L型结构两个臂的交汇位置,为I和II共用,2N-1个天线阵元分别通过接收通道与信号处理机连接;t1与t2之间的基线长度d12以及t1与t2′之间的基线长度d12′、t2与t3之间的基线长度d23以及t2′与t3′之间的基线长度d23′、t3与t4之间的基线长度d34以及t3′与t4′之间的基线长度d34′满足d12=d34,d12′=d34′,且d23-d12≤λ/2,d23′-d12′≤λ/2,当N>4时,t4以及t4′以后的各天线阵元之间的基线长度满足dn,1=kdn-1,1,dn,1′=kdn-1,1′,其中tn表示I中的第n个天线阵元,tn′表示II中的第n个天线阵元,λ表示入射信号的波长,k为正整数;
(2)每个天线阵元接收入射信号并传输:
天线阵中的每个天线阵元通过接收通道将接收的入射信号传输至信号处理机,得到第一一维线阵天线I对应的入射信号集合V=(v1,v2,…,vn,…,vN)和第二一维线阵天线II对应的入射信号集合V′=(v1,v2′,…,vn′,…,vN′),其中:
其中,dn,1和dn,1′表示tn和tn′与t1的基线长度,An和An′分别表示tn和tn′接收的射频信号,θ和β分别表示入射信号的方位角和俯仰角,F(θ,β)表示天线阵元的方向图函数,e为自然常数,j为虚数符号;
(3)信号处理机判断是否存在失效天线阵元:
信号处理机判断入射信号集合V=(v1,v2,…,vn,…,vN)是否存在一个天线阵元的信号量va=0,a∈(1,2,...,N),同时判断入射信号集合V′=(v1,v2′,…,vn′,…,vN′)是否存在一个天线阵元的信号量va′=0,a∈(1,2,...,N),若是,则存在失效天线阵元,并执行步骤(4),否则,则不存在失效天线阵元,并执行步骤(5);
(4)信号处理机获取基线相位差和虚拟基线相位差:
信号处理机获取I的基线相位差(φ1,φ2,...,φm,...,φM)和虚拟基线相位差φ0以及II的基线相位差(φ1′,φ2′,...,φm′,...,φM′)和虚拟基线相位差φ0′,并执行步骤(6),具体步骤为:
(4a)信号处理机获取I的基线相位差(φ1,φ2,...,φm,...,φM)和虚拟基线相位差φ0:
当n=1时,va=v1,信号处理机对V的第二天线阵元t2与其余天线阵元(t3,t4,...,tn,...tN)中的每一个之间的基线相位差进行测量,得到M条基线相位差(φ1,φ2,...,φm,...,φM),同时对第三天线阵元t3与第四天线阵元t4之间的基线相位差φ34进行测量,并利用t2和t3构成的基线与t3和t4构成的基线构造长度为d0=d12-d34虚拟基线,再计算虚拟基线的相位差φ0=φ1-φ34,其中,φm表示第二天线阵元t2与第m+2天线阵元之间的基线相位差,M=N-2;
当n=2时,va=v2,信号处理机对V的第一天线阵元t1与其余天线阵元(t3,t4,...,tn,...tN)中的每一个之间的基线相位差进行测量,得到M条基线相位差(φ1,φ2,...,φm,...,φM),同时对第三天线阵元t3与第四天线阵元t4之间的基线相位差φ34进行测量,并利用t1和t3构成的基线与t3和t4构成的基线构造长度为d0=d13-2d34的虚拟基线,再计算虚拟基线的相位差φ0=φ1-2φ34,其中,d13表示第一天线阵元t1与第三天线阵元t3之间的距离,φm表示第一天线阵元t1与第m+2天线阵元之间的基线相位差;
当n=3时,va=v3,信号处理机对V的第一天线阵元t1与其余天线阵元(t2,t4,t5,...,tn,...tN)中的每一个之间的基线相位差进行测量,得到M条基线相位差(φ1,φ2,...,φm,...,φM),同时对第二天线阵元t2与第四天线阵元t4之间的基线相位差φ24进行测量,并利用t1和t2构成的基线与t2和t4构成的基线构造长度为d0=d24-2d12虚拟基线,再计算虚拟基线的相位差φ0=φ24-2φ1,其中,d24表示第二天线阵元t2与第四天线阵元t4之间的距离,当m=1时,φm表示第一天线阵元t1与第二天线阵元t2之间的基线相位差,当m>1时,φm表示第一天线阵元t1与第m+2天线阵元之间的基线相位差;
当n>3时,va=vn,信号处理机对V的第一天线阵元t1与其余天线阵元(t2,t3,...,tn-1,tn+1,...,tN)中的每一个之间的基线相位差进行测量,得到M条基线相位差(φ1,φ2,...,φm,...,φM),同时对第二天线阵元t2与第三天线阵元t3之间的基线相位差φ23进行测量,并利用t1和t2构成的基线与t2和t3构成的基线构造长度为d0=d23-d12虚拟基线,再计算虚拟基线的相位差φ0=φ23-φ1,其中,当m<n时,φm表示第一天线阵元t1与第m+1天线阵元之间的基线相位差,当m≥n时,φm表示第一天线阵元t1与第m+2天线阵元之间的基线相位差;
(4b)信号处理机获取II的基线相位差(φ1′,φ2′,...,φm′,...,φM′)和虚拟基线相位差φ0′:
信号处理机按照步骤(4a)的方法获取II的基线相位差(φ1′,φ2′,...,φm′,...,φM′)和虚拟基线相位差φ0′;
(5)信号处理机获取基线相位差和虚拟基线相位差:
信号处理机获取I的基线相位差(φ1,φ2,...,φm,...,φM)和虚拟基线相位差φ0以及II的基线相位差(φ1′,φ2′,...,φm′,...,φM′)和虚拟基线相位差φ0′,并执行步骤(6),具体步骤为:
(5a)信号处理机获取I的基线相位差(φ1,φ2,...,φm,...,φM)和虚拟基线相位差φ0:
信号处理机对V的第一天线阵元t1与其余天线阵元(t2,t3,...,tn,...tN)中的每一个之间的基线相位差进行测量,得到M条基线相位差(φ1,φ2,...,φm,...,φM),同时对第二天线阵元t2与第三天线阵元t3之间的基线相位差φ23进行测量,并利用t1和t2构成的基线与t2和t3构成的基线构造长度为d0=d23-d12虚拟基线,再计算虚拟基线的相位差φ0=φ23-φ1,其中,φm表示第一天线阵元t1与第m+1天线阵元之间的基线相位差,M=N-1;
(5b)信号处理机获取II的基线相位差(φ1′,φ2′,...,φm′,...,φM′)和虚拟基线相位差φ0′:
信号处理机对V′的第一天线阵元t1与其余天线阵元(t2′,t3′,...,tn′,...,tN′)中的每一个之间的基线相位差进行测量,得到M′条基线相位差(φ1′,φ2′,...,φm′,...,φM′),同时对第二天线阵元t2′与第三天线阵元t3′之间的基线相位差φ23′进行测量,并利用t1′和t2′构成的基线与t2′和t3′构成的基线构造长度为d0′=d23′-d12′虚拟基线,再计算虚拟基线的相位差φ0′=φ23′-φ1′,其中,φm′表示第一天线阵元t1与第m+1天线阵元之间的基线相位差,M′=N-1;
(6)信号处理机对I和II的基线相位差分别进行解相位模糊:
信号处理机通过虚拟基线无模糊相位差φ0对I的基线相位差(φ1,φ2,...,φm,...,φM)逐级进行解相位模糊,得到I的基线的无模糊相位差同时通过虚拟基线无模糊相位差φ0′对II的基线相位差(φ1′,φ2′,...,φm′,...,φM′)逐级进行解相位模糊,得到II的基线的无模糊相位差其中为φm对应的无模糊相位差,为φm′对应的无模糊相位差;
(7)信号处理机获取相位干涉仪的测向结果:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110326629.5A CN113050029B (zh) | 2021-03-26 | 2021-03-26 | 一种适用于阵元失效条件下的相位干涉仪测向方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110326629.5A CN113050029B (zh) | 2021-03-26 | 2021-03-26 | 一种适用于阵元失效条件下的相位干涉仪测向方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113050029A CN113050029A (zh) | 2021-06-29 |
CN113050029B true CN113050029B (zh) | 2022-12-02 |
Family
ID=76515599
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110326629.5A Active CN113050029B (zh) | 2021-03-26 | 2021-03-26 | 一种适用于阵元失效条件下的相位干涉仪测向方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113050029B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6791493B1 (en) * | 2003-08-21 | 2004-09-14 | Northrop Grumman Corporation | Method for using signal frequency change to differentially resolve long baseline interferometer measurements |
CN102419432A (zh) * | 2011-08-25 | 2012-04-18 | 电子科技大学 | 一种基于虚拟基线的圆阵相位干涉仪二维测向方法 |
CN104330767A (zh) * | 2014-11-27 | 2015-02-04 | 中国船舶重工集团公司第七二四研究所 | 基于相位模糊数搜索与最小二乘拟合的干涉仪测向方法 |
CN105259536A (zh) * | 2015-09-17 | 2016-01-20 | 哈尔滨工程大学 | 一种基于l型天线阵列的相位干涉仪测向装置及解算方法 |
CN109270509A (zh) * | 2018-11-09 | 2019-01-25 | 深圳大学 | 数据损失情况下基于矩阵填充的doa估计方法及系统 |
CN109633526A (zh) * | 2019-01-29 | 2019-04-16 | 杭州电子科技大学 | 基于方向函数的非均匀圆阵相位干涉仪测向解模糊方法 |
CN112485756A (zh) * | 2020-11-20 | 2021-03-12 | 中国电子科技集团公司第二十九研究所 | 基于dbf阵列天线的干涉仪信号处理方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4453164A (en) * | 1982-07-26 | 1984-06-05 | Rca Corporation | Method of determining excitation of individual elements of a phase array antenna from near-field data |
KR101357690B1 (ko) * | 2012-10-16 | 2014-02-04 | 국방과학연구소 | 방향탐지용 인터페로미터 배열 안테나 이격비 산출방법 |
CN111044971B (zh) * | 2019-12-23 | 2022-07-01 | 南京长峰航天电子科技有限公司 | 一种二维干涉仪测向系统 |
-
2021
- 2021-03-26 CN CN202110326629.5A patent/CN113050029B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6791493B1 (en) * | 2003-08-21 | 2004-09-14 | Northrop Grumman Corporation | Method for using signal frequency change to differentially resolve long baseline interferometer measurements |
CN102419432A (zh) * | 2011-08-25 | 2012-04-18 | 电子科技大学 | 一种基于虚拟基线的圆阵相位干涉仪二维测向方法 |
CN104330767A (zh) * | 2014-11-27 | 2015-02-04 | 中国船舶重工集团公司第七二四研究所 | 基于相位模糊数搜索与最小二乘拟合的干涉仪测向方法 |
CN105259536A (zh) * | 2015-09-17 | 2016-01-20 | 哈尔滨工程大学 | 一种基于l型天线阵列的相位干涉仪测向装置及解算方法 |
CN109270509A (zh) * | 2018-11-09 | 2019-01-25 | 深圳大学 | 数据损失情况下基于矩阵填充的doa估计方法及系统 |
CN109633526A (zh) * | 2019-01-29 | 2019-04-16 | 杭州电子科技大学 | 基于方向函数的非均匀圆阵相位干涉仪测向解模糊方法 |
CN112485756A (zh) * | 2020-11-20 | 2021-03-12 | 中国电子科技集团公司第二十九研究所 | 基于dbf阵列天线的干涉仪信号处理方法 |
Non-Patent Citations (4)
Title |
---|
A sparse representation based DOA estimation algorithm for a mixture of circular and noncircular signals using sparse arrays;Jingjing Cai et al.;《2017 IEEE International Conference on Signal Processing, Communications and Computing (ICSPCC)》;20180101;第1-5页 * |
DOA Estimation Using Virtual ESPRIT with Successive Baselines and Coprime Baselines;Xiaolin Li et al.;《Circuits, Systems, and Signal Processing》;20201006;第40卷;第2065-2075页 * |
二维干涉仪测向技术研究;刘满朝等;《现代雷达》;20180228;第40卷(第2期);第44-48页 * |
利用基线长度和俯仰角修正改进的二维干涉仪测向算法;肖香彬等;《电子信息对抗技术》;20190131;第34卷(第1期);第5-8、26页 * |
Also Published As
Publication number | Publication date |
---|---|
CN113050029A (zh) | 2021-06-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109669178B (zh) | 一种星载三阵元单脉冲两维测向方法 | |
CN108875099B (zh) | 一种基于长短基线干涉仪测向体制的基线选取方法 | |
CN110187303B (zh) | 一种单脉冲比幅比相联合测向方法 | |
CN101963659A (zh) | 利用三维正交干涉仪实现宽频段测向的方法 | |
Chen et al. | Two-dimensional direction of arrival estimation for improved archimedean spiral array with MUSIC algorithm | |
CN116500543B (zh) | 一种基于参考方向变换的来波角度快速估计方法 | |
CN113050029B (zh) | 一种适用于阵元失效条件下的相位干涉仪测向方法 | |
CN114487993B (zh) | 一种均匀圆阵相关干涉仪测向方法及系统 | |
CN111693935A (zh) | 一种射频测向方法及系统 | |
CN114487992B (zh) | 无须解模糊的多基线相位干涉仪测向方法 | |
Yang et al. | Coprime L‐shaped array connected by a triangular spatially‐spread electromagnetic‐vector‐sensor for two‐dimensional direction of arrival estimation | |
CN105334490A (zh) | 一种宽带信号干涉测向的方法 | |
Mollai et al. | Wideband two dimensional interferometric direction finding algorithm using base-triangles and a proposed minimum planar array | |
CN116070082A (zh) | 一种基于空间谱估计的极化敏感阵列幅相校正方法 | |
CN113504554B (zh) | 基于星载相控阵的非均匀子阵合成干涉仪测向方法 | |
CN112946615B (zh) | 一种相控阵系统幅相误差矫正方法 | |
CN111693936B (zh) | 一种基于射频模拟接收系统的测向方法及系统 | |
CN111431573B (zh) | 一种csi单接入点定位的天线资源配置与阵列设计方法 | |
CN114265005A (zh) | 一种极化相位干涉仪测向方法及装置 | |
Memarian et al. | Multiple signals direction finding of IoT devices through improved correlative interferometer using directional elements | |
CN114460531A (zh) | 一种均匀线阵music空间谱估计方法 | |
CN115906372A (zh) | 一种基于稀疏阵幅相误差校正的无网格参数估计方法 | |
Huang et al. | In-field calibration of passive array receiver using detected target | |
Zhu et al. | Mirrored arrays for direction-of-arrival estimation | |
CN116359835B (zh) | 一种y型基线干涉仪测向装置及测向方法 |
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 | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20230728 Address after: Room A601, National Digital Publishing Base, No. 996 Tiangu 7th Road, Gaoxin, Xi'an City, Shaanxi Province, 710000 Patentee after: XI'AN THUNDER DRIVEN ELECTRONIC TECHNOLOGY CO.,LTD. Address before: 710071 No. 2 Taibai South Road, Shaanxi, Xi'an Patentee before: XIDIAN University |