CN113050029A - 一种适用于阵元失效条件下的相位干涉仪测向方法 - Google Patents

Abstract

本发明提出了一种在阵元失效条件下干涉仪测向方法，旨在实现相位干涉仪在阵元失效情况下测向。实现步骤为：构建相位干涉仪测向系统；每个天线阵元接收入射信号并传输；信号处理机判断是否存在失效天线阵元；信号处理机根据阵元失效情况获取基线相位差和虚拟基线相位差；信号处理机对基线相位解模糊；信号处理机获取相位干涉仪的测向结果。本发明设计的非均匀L型天线阵，可以利用基线长度关系构造虚拟基线进行解模糊，实现在阵元失效的情况下测向，同时在阵元完好的情况也能保证测向精度，提高了相位干涉仪测向系统的鲁棒性。

Description

一种适用于阵元失效条件下的相位干涉仪测向方法

技术领域

本发明属于无源测向技术领域，特别涉及一种在阵元失效条件下的相位干涉仪测向方法。

背景技术

传统的有源雷达测向技术过发射大功率信号来实现，具有作用时间长、距离远、精度高等优点，并且受外界环境因素变化的影响小。但是随着电子以及通信信息技术的发展，有源测向技术在被广泛使用的同时，其系统自身隐藏与抗干扰能力差的缺点凸显，无法隐藏容易遭到干扰与攻击。为了弥补有源测向技术的不足与可能产生的严重后果，人们开始研究无源测向技术。

相位干涉仪测向技术即属于无源测向技术的一种，它本身不会辐射信源，只接收处理来自目标的位置及运动状态的相关数据来实现对目标的测向。当入射信号到达干涉仪测向天线阵时，由于天线的空间位置不同，每个天线单元接收到的信号相位也不同，干涉仪就是利用天线单元之间的相位差来获取来波的方向信息。干涉仪测向系统具有灵敏度高、测向精准、算法简单等优点，广泛地用于电磁环境监测、电子对抗等领域。

相位干涉仪主要利用来波在基线间产生的相位差来计算方位，影响其测向精度和效率主要取决于布阵的形式和测向算法。如果出现阵元失效的情况，则会破坏原有阵列的几何分布特征，丢失该通道的数据，导致整个系统的性能下降，甚至无法工作。特别是在机载等天线布置空间有限的平台上，允许天线阵元的尺寸和数量都非常有限，往往出现某个阵元失效就会导致整个设备无法使用。

目前相位干涉仪的相关研究都是基于阵元完好的情况，例如申请公布号CN111323742A，名称为“一种基于曲线阵列的相位干涉仪及其测向方法”的专利申请中，公开了一种基于曲线阵列的相位干涉仪及其测向方法，该方法通过保证一组3个天线的阵列的纵向间隔相等，在横坐标方向和纵坐标方向分别形成虚拟基线，抵消天线位置的纵坐标带来的影响，将曲线阵列模型等效为一维阵列模型进行测向。而干涉仪通常都在室外工作，由于恶劣的自然环境、人为干扰以及阵元老化、阵元的物理损坏等因素的影响，阵元失效问题不可避免，对应阵元失效情况，申请公布号CN111323742A所述的方法就无法保证选取的一组3个天线的阵列的纵向间隔相等，从而导致算法失效。

发明内容

本发明的目的在于克服上述现有技术存在的缺陷，提出一种适用于阵元失效条件下的相位干涉仪测向方法，用于解决现有技术中存在的出现阵元失效情况相位干涉仪无法测向的问题。

为实现上述目的，本发明采取的技术方案包括如下步骤：

(1)构建相位干涉仪测向系统：

构建包括天线阵和信号处理机的相位干涉仪测向系统，所述天线阵采用由正交的第一一维线阵天线I和第二一维线阵天线II组成的L型结构，所述第一一维线阵天线I包括非均匀排布的N个天线阵元(t₁,t₂,...,t_n,...,t_N)，所述第二一维线阵天线II包括非均匀排布的N个天线阵元(t₁,t₂′,...,t_n′,...,t_N′)，N≥4，第一天线阵元t₁位于L型结构两个臂的交汇位置，为I和II共用，2N-1个天线阵元分别通过接收通道与信号处理机连接；t₁与t₂之间的基线长度d₁₂以及t₁与t₂′之间的基线长度d₁₂′、t₂与t₃之间的基线长度d₂₃以及t₂′与t₃′之间的基线长度d₂₃′、t₃与t₄之间的基线长度d₃₄以及t₃′与t₄′之间的基线长度d₃₄′满足d₁₂＝d₃₄，d₁₂′＝d₃₄′，且d₂₃-d₁₂≤λ/2，d₂₃′-d₁₂′≤λ/2，当N＞4时，t₄以及t₄′以后的各天线阵元之间的基线长度满足d_n,1＝kd_n-1,1，d_n,1′＝kd_n-1,1′(n＞4)，其中t_n表示I中的第n个天线阵元，t_n′表示II中的第n个天线阵元，λ表示入射信号的波长，k为正整数；

(2)每个天线阵元接收入射信号并传输：

天线阵中的每个天线阵元通过接收通道将接收的入射信号传输至信号处理机，得到第一一维线阵天线I对应的入射信号集合V＝(v₁,v₂,…,v_n,…,v_N)和第二一维线阵天线II对应的入射信号集合V′＝(v₁,v₂′,…,v_n′,…,v_N′)，其中：

其中，d_n,1和d_n,1′表示t_n和t_n′与t₁的基线长度，A_n和A_n′分别表示t_n和t_n′接收的射频信号，θ和β分别表示入射信号的方位角和俯仰角，F(θ,β)表示天线阵元的方向图函数，e为自然常数，j为虚数符号；

(3)信号处理机判断是否存在失效天线阵元：

信号处理机判断入射信号集合V＝(v₁,v₂,…,v_n,…,v_N)是否存在一个天线阵元的信号量v_a＝0，a∈(1,2,...,N)，同时判断入射信号集合V′＝(v₁,v₂′,…,v_n′,…,v_N′)是否存在一个天线阵元的信号量v_a′＝0，a∈(1,2,...,N)，若是，则存在失效天线阵元，并执行步骤(4)，否则，则不存在失效天线阵元，并执行步骤(5)；

(4)信号处理机获取基线相位差和虚拟基线相位差：

信号处理机获取I的基线相位差(φ₁,φ₂,...,φ_m,...,φ_M)和虚拟基线相位差φ₀以及II的基线相位差(φ₁′,φ₂′,...,φ_m′,...,φ_M′)和虚拟基线相位差φ₀′，并执行步骤(6)，具体步骤为：

(4a)信号处理机获取I的基线相位差(φ₁,φ₂,...,φ_m,...,φ_M)和虚拟基线相位差φ₀：

当n＝1时，v_a＝v₁，信号处理机对V的第二天线阵元t₂与其余天线阵元(t₃,t₄,...,t_n,...t_N)中的每一个之间的基线相位差进行测量，得到M条基线相位差(φ₁,φ₂,...,φ_m,...,φ_M)，同时对第三天线阵元t₃与第四天线阵元t₄之间的基线相位差φ₃₄进行测量，并利用t₂和t₃构成的基线与t₃和t₄构成的基线构造长度为d₀＝d₁₂-d₃₄虚拟基线，再计算虚拟基线的相位差φ₀＝φ₁-φ₃₄，其中，φ_m表示第二天线阵元t₂与第m+2天线阵元之间的基线相位差，M＝N-2；

当n＝2时，v_a＝v₂，信号处理机对V的第一天线阵元t₁与其余天线阵元(t₃,t₄,...,t_n,...t_N)中的每一个之间的基线相位差进行测量，得到M条基线相位差(φ₁,φ₂,...,φ_m,...,φ_M)，同时对第三天线阵元t₃与第四天线阵元t₄之间的基线相位差φ₃₄进行测量，并利用t₁和t₃构成的基线与t₃和t₄构成的基线构造长度为d₀＝d₁₃-2d₃₄的虚拟基线，再计算虚拟基线的相位差φ₀＝φ₁-2φ₃₄，其中，d₁₃表示第一天线阵元t₁与第三天线阵元t₃之间的距离，φ_m表示第一天线阵元t₁与第m+2天线阵元之间的基线相位差；

当n＝3时，v_a＝v₃，信号处理机对V的第一天线阵元t₁与其余天线阵元(t₂,t₄,t₅,...,t_n,...t_N)中的每一个之间的基线相位差进行测量，得到M条基线相位差(φ₁,φ₂,...,φ_m,...,φ_M)，同时对第二天线阵元t₂与第四天线阵元t₄之间的基线相位差φ₂₄进行测量，并利用t₁和t₂构成的基线与t₂和t₄构成的基线构造长度为d₀＝d₂₄-2d₁₂虚拟基线，再计算虚拟基线的相位差φ₀＝φ₂₄-2φ₁，其中，d₂₄表示第二天线阵元t₂与第四天线阵元t₄之间的距离，当m＝1时，φ_m表示第一天线阵元t₁与第二天线阵元t₂之间的基线相位差，当m＞1时，φ_m表示第一天线阵元t₁与第m+2天线阵元之间的基线相位差；

当n＞3时，v_a＝v_n，信号处理机对V的第一天线阵元t₁与其余天线阵元(t₂,t₃,...,t_n-1,t_n+1,...,t_N)中的每一个之间的基线相位差进行测量，得到M条基线相位差(φ₁,φ₂,...,φ_m,...,φ_M)，同时对第二天线阵元t₂与第三天线阵元t₃之间的基线相位差φ₂₃进行测量，并利用t₁和t₂构成的基线与t₂和t₃构成的基线构造长度为d₀＝d₂₃-d₁₂虚拟基线，再计算虚拟基线的相位差φ₀＝φ₂₃-φ₁，其中，当m＜n时，φ_m表示第一天线阵元t₁与第m+1天线阵元之间的基线相位差，当m≥n时，φ_m表示第一天线阵元t₁与第m+2天线阵元之间的基线相位差；

(4b)信号处理机获取II的基线相位差(φ₁′,φ₂′,...,φ_m′,...,φ_M′)和虚拟基线相位差φ₀′：

信号处理机按照步骤(4a)的方法获取II的基线相位差(φ₁′,φ₂′,...,φ_m′,...,φ_M′)和虚拟基线相位差φ₀′；

(5)信号处理机获取基线相位差和虚拟基线相位差：

信号处理机获取I的基线相位差(φ₁,φ₂,...,φ_m,...,φ_M)和虚拟基线相位差φ₀以及II的基线相位差(φ₁′,φ₂′,...,φ_m′,...,φ_M′)和虚拟基线相位差φ₀′，并执行步骤(6)，具体步骤为：

(5a)信号处理机获取I的基线相位差(φ₁,φ₂,...,φ_m,...,φ_M)和虚拟基线相位差φ₀：

信号处理机对V的第一天线阵元t₁与其余天线阵元(t₂,t₃,...,t_n,...t_N)中的每一个之间的基线相位差进行测量，得到M条基线相位差(φ₁,φ₂,...,φ_m,...,φ_M)，同时对第二天线阵元t₂与第三天线阵元t₃之间的基线相位差φ₂₃进行测量，并利用t₁和t₂构成的基线与t₂和t₃构成的基线构造长度为d₀＝d₂₃-d₁₂虚拟基线，再计算虚拟基线的相位差φ₀＝φ₂₃-φ₁，其中，φ_m表示第一天线阵元t₁与第m+1天线阵元之间的基线相位差，M＝N-1；

(5b)信号处理机获取II的基线相位差(φ₁′,φ₂′,...,φ_m′,...,φ_M′)和虚拟基线相位差φ₀′：

信号处理机对V′的第一天线阵元t₁与其余天线阵元(t₂′,t₃′,...,t_n′,...,t_N′)中的每一个之间的基线相位差进行测量，得到M′条基线相位差(φ₁′,φ₂′,...,φ_m′,...,φ_M′)，同时对第二天线阵元t₂′与第三天线阵元t₃′之间的基线相位差φ₂₃′进行测量，并利用t₁′和t₂′构成的基线与t₂′和t₃′构成的基线构造长度为d₀′＝d₂₃′-d₁₂′虚拟基线，再计算虚拟基线的相位差φ₀′＝φ₂₃′-φ₁′，其中，φ_m′表示第一天线阵元t₁与第m+1天线阵元之间的基线相位差，M′＝N-1；

(6)信号处理机对I和II的基线相位差分别进行解相位模糊：

信号处理机通过虚拟基线无模糊相位差φ₀对I的基线相位差(φ₁,φ₂,...,φ_m,...,φ_M)逐级进行解相位模糊，得到I的基线的无模糊相位差

同时通过虚拟基线无模糊相位差φ₀′对II的基线相位差(φ₁′,φ₂′,...,φ_m′,...,φ_M′)逐级进行解相位模糊，得到II的基线的无模糊相位差

的计算公式为：

的计算公式为：

其中(φ₁,φ₂,...,φ_m,...,φ_M)，(φ₁′,φ₂′,...,φ_m′,...,φ_M′)表示执行步骤(4)或步骤(5)得到的I和II的基线相位差，d_m表示

对应的基线长度，d_m′表示φ_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型结构天线，基线长度满足：d₁₂＝d₃₄，d₁₂′＝d₃₄′，且d₂₃-d₁₂≤λ/2，d₂₃′-d₁₂′≤λ/2，当N＞4时，t₄以及t₄′以后的各天线阵元之间的基线长度满足d_n,1＝kd_n-1,1，d_n,1′＝kd_n-1,1′(n＞4)，在此基线长度关系下，可通过构造虚拟基线进行解相位模糊，实现在阵元失效的情况下测向，同时在阵元完好的情况下也能保证测向精度，提高了干涉仪测向系统的鲁棒性。

附图说明

图1是本发明的实现流程图；

图2是本发明中构建的相位干涉仪测向系统中天线阵的结构示意图；

具体实施方式

下面结合附图和具体实施例，对本发明作进一步详细描述。

参照图1，本发明包括如下步骤：

步骤1)构建相位干涉仪测向系统：

构建相位干涉仪测向系统，包括天线阵和信号处理机的相位干涉仪测向系统，所述的天线阵结构如图2所示，采用由正交的第一一维线阵天线I和第二一维线阵天线II组成的L型结构，所述第一一维线阵天线I包括非均匀排布的N个天线阵元(t₁,t₂,...,t_n,...,t_N)，所述第二一维线阵天线II包括非均匀排布的N个天线阵元(t₁,t₂′,...,t_n′,...,t_N′)，N≥4，第一天线阵元t₁位于L型结构两个臂的交汇位置，为I和II共用，2N-1个天线阵元分别通过接收通道与信号处理机连接；t₁与t₂之间的基线长度d₁₂以及t₁与t₂′之间的基线长度d₁₂′、t₂与t₃之间的基线长度d₂₃以及t₂′与t₃′之间的基线长度d₂₃′、t₃与t₄之间的基线长度d₃₄以及t₃′与t₄′之间的基线长度d₃₄′满足d₁₂＝d₃₄，d₁₂′＝d₃₄′，且d₂₃-d₁₂≤λ/2，d₂₃′-d₁₂′≤λ/2，当N＞4时，t₄以及t₄′以后的各天线阵元之间的基线长度满足d_n,1＝kd_n-1,1，d_n,1′＝kd_n-1,1′(n＞4)，其中t_n表示I中的第n个天线阵元，t_n′表示II中的第n个天线阵元，λ表示入射信号的波长，k为正整数；

本实施例中，设有如下要求设计天线阵：覆盖频段0.8GHz～6GHz；天线圆盘直径≤300mm；采用平面螺旋天线，直径r＝70mm，取N＝4。系统的最高测向频率为6GHz，对应的λ＝50mm，要获得无模糊相位，显然物理基线长度不可能做到50mm，需要使用虚拟基线的方式，则d₂₃-d₁₂≤25，d₂₃′-d₁₂′≤25，天线圆盘直径≤300mm，则d₁₂+d₂₃+d₃₄≤300，d₁₂′+d₂₃′+d₃₄′≤300，故可设d₁₂＝d₁₂′＝90mm，d₂₃＝d₂₃′＝115mm，d₃₄＝d₃₄′＝90mm。

步骤2)每个天线阵元接收入射信号并传输：

如图2所示，辐射源发射信号，天线阵接收到的入射信号方位角和俯仰角为θ和β，θ表示辐射源的阵面投影和t₁构成直线与II之间的夹角，β表示辐射源和t₁构成直线与阵面之间的夹角。

天线阵中的每个天线阵元通过接收通道将接收的入射信号传输至信号处理机，得到第一一维线阵天线I对应的入射信号集合V＝(v₁,v₂,v₃,v₄)和第二一维线阵天线II对应的入射信号集合V′＝(v₁,v₂′,v₃′,v₄′)，其中：

其中，n∈(1,2,3,4)，d_n,1和d_n,1′表示t_n和t_n′与t₁的基线长度，A_n和A_n′分别表示t_n和t_n′接收的射频信号，F(θ,β)表示天线阵元的方向图函数，e为自然常数，j为虚数符号；

步骤3)信号处理机判断是否存在失效天线阵元：

信号处理机判断入射信号集合V＝(v₁,v₂,v₃,v₄)是否存在一个天线阵元的信号量v_a＝0，a∈(1,2,3,4)，同时判断入射信号集合V′＝(v₁,v₂′,v₃′,v₄′)是否存在一个天线阵元的信号量v_a′＝0，a∈(1,2,3,4)，若是，则存在失效天线阵元，并执行步骤(4)，否则，则不存在失效天线阵元，并执行步骤(5)；

步骤4)信号处理机获取基线相位差和虚拟基线相位差：

信号处理机获取I的基线相位差(φ₁,φ₂)和虚拟基线相位差φ₀以及II的基线相位差(φ₁′,φ₂′)和虚拟基线相位差φ₀′，并执行步骤(6)，具体步骤为：

(4a)信号处理机获取I的基线相位差(φ₁,φ₂)和虚拟基线相位差φ₀：

当v_a＝v₁，信号处理机对V的第二天线阵元t₂与第三天线阵元t₃之间的基线相位φ₁、第二天线阵元t₂与第三天线阵元t₄之间的基线相位φ₂、第三天线阵元t₃与第四天线阵元t₄之间的基线相位差φ₃₄，并利用t₂和t₃构成的基线与t₃和t₄构成的基线构造长度为d₀＝d₁₂-d₃₄＝25mm虚拟基线，再计算虚拟基线的相位差φ₀＝φ₁-φ₃₄；

当v_a＝v₂，信号处理机对V的第二天线阵元t₁与第三天线阵元t₃之间的基线相位φ₁、第二天线阵元t₁与第三天线阵元t₄之间的基线相位φ₂、第三天线阵元t₃与第四天线阵元t₄之间的基线相位差φ₃₄进行测量，并利用t₁和t₃构成的基线与t₃和t₄构成的基线构造长度为d₀＝d₁₃-2d₃₄＝25mm的虚拟基线，再计算虚拟基线的相位差φ₀＝φ₁-2φ₃₄，其中，d₁₃表示第一天线阵元t₁与第三天线阵元t₃之间的距离，d₁₃＝d₁₂+d₂₃＝205mm；

当v_a＝v₃，信号处理机对V的第二天线阵元t₁与第三天线阵元t₂之间的基线相位φ₁、第二天线阵元t₁与第三天线阵元t₄之间的基线相位φ₂、第二天线阵元t₂与第四天线阵元t₄之间的基线相位差φ₂₄进行测量，并利用t₁和t₂构成的基线与t₂和t₄构成的基线构造长度为d₀＝d₂₄-2d₁₂＝25mm虚拟基线，再计算虚拟基线的相位差φ₀＝φ₂₄-2φ₁，其中，d₂₄表示第二天线阵元t₂与第四天线阵元t₄之间的距离，d₂₄＝d₂₃+d₃₄＝205mm；

当v_a＝v₄，信号处理机对V的第二天线阵元t₁与第三天线阵元t₂之间的基线相位φ₁、第二天线阵元t₁与第三天线阵元t₃之间的基线相位φ₂、第二天线阵元t₂与第三天线阵元t₃之间的基线相位差φ₂₃进行测量，并利用t₁和t₂构成的基线与t₂和t₃构成的基线构造长度为d₀＝d₂₃-d₁₂＝25mm虚拟基线，再计算虚拟基线的相位差φ₀＝φ₂₃-φ₁；

(4b)信号处理机获取II的基线相位差(φ₁′,φ₂′)和虚拟基线相位差φ₀′：

信号处理机按照步骤(4a)的方法获取II的基线相位差(φ₁′,φ₂′)和虚拟基线相位差φ₀′；

步骤5)信号处理机获取基线相位差和虚拟基线相位差：

信号处理机获取I的基线相位差(φ₁,φ₂,φ₃)和虚拟基线相位差φ₀以及II的基线相位差(φ₁′,φ₂′,φ₃′)和虚拟基线相位差φ₀′，并执行步骤(6)，具体步骤为：

(5a)信号处理机获取I的基线相位差(φ₁,φ₂,φ₃)和虚拟基线相位差φ₀：

信号处理机对V的第一天线阵元t₁与其余天线阵元(t₂,t₃,t₄)中的每一个之间的基线相位差进行测量，得到3条基线相位差(φ₁,φ₂,φ₃)，同时对第二天线阵元t₂与第三天线阵元t₃之间的基线相位差φ₂₃进行测量，并利用t₁和t₂构成的基线与t₂和t₃构成的基线构造长度为d₀＝d₂₃-d₁₂＝25mm虚拟基线，再计算虚拟基线的相位差φ₀＝φ₂₃-φ₁；

(5b)信号处理机获取II的基线相位差(φ₁′,φ₂′,φ₃′)和虚拟基线相位差φ₀′：

信号处理机对V′的第一天线阵元t₁与其余天线阵元(t₂′,t₃′,t₄′)中的每一个之间的基线相位差进行测量，得到3条基线相位差(φ₁′,φ₂′,φ₃′)，同时对第二天线阵元t₂′与第三天线阵元t₃′之间的基线相位差φ₂₃′进行测量，并利用t₁′和t₂′构成的基线与t₂′和t₃′构成的基线构造长度为d₀′＝d₂₃′-d₁₂′＝25mm虚拟基线，再计算虚拟基线的相位差φ₀′＝φ₂₃′-φ₁′；

步骤6)信号处理机对I和II的基线相位差分别进行解相位模糊：

信号处理机通过虚拟基线无模糊相位差φ₀对I的基线相位差(φ₁,φ₂,...,φ_m,...,φ_M)逐级进行解相位模糊，得到I的基线的无模糊相位差

同时通过虚拟基线无模糊相位差φ₀′对II的基线相位差(φ₁′,φ₂′,...,φ_m′,...,φ_M′)逐级进行解相位模糊，得到II的基线的无模糊相位差

的计算公式为：

的计算公式为：

其中d_m表示

对应的基线长度，d_m′表示φ_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个天线阵元(t₁,t₂,...,t_n,...,t_N)，所述第二一维线阵天线II包括非均匀排布的N个天线阵元(t₁,t₂′,...,t_n′,...,t_N′)，N≥4，第一天线阵元t₁位于L型结构两个臂的交汇位置，为I和II共用，2N-1个天线阵元分别通过接收通道与信号处理机连接；t₁与t₂之间的基线长度d₁₂以及t₁与t₂′之间的基线长度d₁₂′、t₂与t₃之间的基线长度d₂₃以及t₂′与t₃′之间的基线长度d₂₃′、t₃与t₄之间的基线长度d₃₄以及t₃′与t₄′之间的基线长度d₃₄′满足d₁₂＝d₃₄，d₁₂′＝d₃₄′，且d₂₃-d₁₂≤λ/2，d₂₃′-d₁₂′≤λ/2，当N＞4时，t₄以及t₄′以后的各天线阵元之间的基线长度满足d_n,1＝kd_n-1,1，d_n,1′＝kd_n-1,1′(n＞4)，其中t_n表示I中的第n个天线阵元，t_n′表示II中的第n个天线阵元，λ表示入射信号的波长，k为正整数；

(2)每个天线阵元接收入射信号并传输：

天线阵中的每个天线阵元通过接收通道将接收的入射信号传输至信号处理机，得到第一一维线阵天线I对应的入射信号集合V＝(v₁,v₂,…,v_n,…,v_N)和第二一维线阵天线II对应的入射信号集合V′＝(v₁,v₂′,…,v_n′,…,v_N′)，其中：

其中，d_n,1和d_n,1′表示t_n和t_n′与t₁的基线长度，A_n和A_n′分别表示t_n和t_n′接收的射频信号，θ和β分别表示入射信号的方位角和俯仰角，F(θ,β)表示天线阵元的方向图函数，e为自然常数，j为虚数符号；

(3)信号处理机判断是否存在失效天线阵元：

信号处理机判断入射信号集合V＝(v₁,v₂,…,v_n,…,v_N)是否存在一个天线阵元的信号量v_a＝0，a∈(1,2,...,N)，同时判断入射信号集合V′＝(v₁,v₂′,…,v_n′,…,v_N′)是否存在一个天线阵元的信号量v_a′＝0，a∈(1,2,...,N)，若是，则存在失效天线阵元，并执行步骤(4)，否则，则不存在失效天线阵元，并执行步骤(5)；

(4)信号处理机获取基线相位差和虚拟基线相位差：

信号处理机获取I的基线相位差(φ₁,φ₂,...,φ_m,...,φ_M)和虚拟基线相位差φ₀以及II的基线相位差(φ₁′,φ₂′,...,φ_m′,...,φ_M′)和虚拟基线相位差φ₀′，并执行步骤(6)，具体步骤为：

(4a)信号处理机获取I的基线相位差(φ₁,φ₂,...,φ_m,...,φ_M)和虚拟基线相位差φ₀：

当n＝1时，v_a＝v₁，信号处理机对V的第二天线阵元t₂与其余天线阵元(t₃,t₄,...,t_n,...t_N)中的每一个之间的基线相位差进行测量，得到M条基线相位差(φ₁,φ₂,...,φ_m,...,φ_M)，同时对第三天线阵元t₃与第四天线阵元t₄之间的基线相位差φ₃₄进行测量，并利用t₂和t₃构成的基线与t₃和t₄构成的基线构造长度为d₀＝d₁₂-d₃₄虚拟基线，再计算虚拟基线的相位差φ₀＝φ₁-φ₃₄，其中，φ_m表示第二天线阵元t₂与第m+2天线阵元之间的基线相位差，M＝N-2；

当n＝2时，v_a＝v₂，信号处理机对V的第一天线阵元t₁与其余天线阵元(t₃,t₄,...,t_n,...t_N)中的每一个之间的基线相位差进行测量，得到M条基线相位差(φ₁,φ₂,...,φ_m,...,φ_M)，同时对第三天线阵元t₃与第四天线阵元t₄之间的基线相位差φ₃₄进行测量，并利用t₁和t₃构成的基线与t₃和t₄构成的基线构造长度为d₀＝d₁₃-2d₃₄的虚拟基线，再计算虚拟基线的相位差φ₀＝φ₁-2φ₃₄，其中，d₁₃表示第一天线阵元t₁与第三天线阵元t₃之间的距离，φ_m表示第一天线阵元t₁与第m+2天线阵元之间的基线相位差；

当n＝3时，v_a＝v₃，信号处理机对V的第一天线阵元t₁与其余天线阵元(t₂,t₄,t₅,...,t_n,...t_N)中的每一个之间的基线相位差进行测量，得到M条基线相位差(φ₁,φ₂,...,φ_m,...,φ_M)，同时对第二天线阵元t₂与第四天线阵元t₄之间的基线相位差φ₂₄进行测量，并利用t₁和t₂构成的基线与t₂和t₄构成的基线构造长度为d₀＝d₂₄-2d₁₂虚拟基线，再计算虚拟基线的相位差φ₀＝φ₂₄-2φ₁，其中，d₂₄表示第二天线阵元t₂与第四天线阵元t₄之间的距离，当m＝1时，φ_m表示第一天线阵元t₁与第二天线阵元t₂之间的基线相位差，当m＞1时，φ_m表示第一天线阵元t₁与第m+2天线阵元之间的基线相位差；

当n＞3时，v_a＝v_n，信号处理机对V的第一天线阵元t₁与其余天线阵元(t₂,t₃,...,t_n-1,t_n+1,...,t_N)中的每一个之间的基线相位差进行测量，得到M条基线相位差(φ₁,φ₂,...,φ_m,...,φ_M)，同时对第二天线阵元t₂与第三天线阵元t₃之间的基线相位差φ₂₃进行测量，并利用t₁和t₂构成的基线与t₂和t₃构成的基线构造长度为d₀＝d₂₃-d₁₂虚拟基线，再计算虚拟基线的相位差φ₀＝φ₂₃-φ₁，其中，当m＜n时，φ_m表示第一天线阵元t₁与第m+1天线阵元之间的基线相位差，当m≥n时，φ_m表示第一天线阵元t₁与第m+2天线阵元之间的基线相位差；

(4b)信号处理机获取II的基线相位差(φ₁′,φ₂′,...,φ_m′,...,φ_M′)和虚拟基线相位差φ₀′：

信号处理机按照步骤(4a)的方法获取II的基线相位差(φ₁′,φ₂′,...,φ_m′,...,φ_M′)和虚拟基线相位差φ₀′；

(5)信号处理机获取基线相位差和虚拟基线相位差：

信号处理机获取I的基线相位差(φ₁,φ₂,...,φ_m,...,φ_M)和虚拟基线相位差φ₀以及II的基线相位差(φ₁′,φ₂′,...,φ_m′,...,φ_M′)和虚拟基线相位差φ₀′，并执行步骤(6)，具体步骤为：

(5a)信号处理机获取I的基线相位差(φ₁,φ₂,...,φ_m,...,φ_M)和虚拟基线相位差φ₀：

信号处理机对V的第一天线阵元t₁与其余天线阵元(t₂,t₃,...,t_n,...t_N)中的每一个之间的基线相位差进行测量，得到M条基线相位差(φ₁,φ₂,...,φ_m,...,φ_M)，同时对第二天线阵元t₂与第三天线阵元t₃之间的基线相位差φ₂₃进行测量，并利用t₁和t₂构成的基线与t₂和t₃构成的基线构造长度为d₀＝d₂₃-d₁₂虚拟基线，再计算虚拟基线的相位差φ₀＝φ₂₃-φ₁，其中，φ_m表示第一天线阵元t₁与第m+1天线阵元之间的基线相位差，M＝N-1；

(5b)信号处理机获取II的基线相位差(φ₁′,φ₂′,...,φ_m′,...,φ_M′)和虚拟基线相位差φ₀′：

信号处理机对V′的第一天线阵元t₁与其余天线阵元(t₂′,t₃′,...,t_n′,...,t_N′)中的每一个之间的基线相位差进行测量，得到M′条基线相位差(φ₁′,φ₂′,...,φ_m′,...,φ_M′)，同时对第二天线阵元t₂′与第三天线阵元t₃′之间的基线相位差φ₂₃′进行测量，并利用t₁′和t₂′构成的基线与t₂′和t₃′构成的基线构造长度为d₀′＝d₂₃′-d₁₂′虚拟基线，再计算虚拟基线的相位差φ₀′＝φ₂₃′-φ₁′，其中，φ_m′表示第一天线阵元t₁与第m+1天线阵元之间的基线相位差，M′＝N-1；

(6)信号处理机对I和II的基线相位差分别进行解相位模糊：

信号处理机通过虚拟基线无模糊相位差φ₀对I的基线相位差(φ₁,φ₂,...,φ_m,...,φ_M)逐级进行解相位模糊，得到I的基线的无模糊相位差

同时通过虚拟基线无模糊相位差φ₀′对II的基线相位差(φ₁′,φ₂′,...,φ_m′,...,φ_M′)逐级进行解相位模糊，得到II的基线的无模糊相位差

其中

为φ_m对应的无模糊相位差，

为φ_m′对应的无模糊相位差；

(7)信号处理机获取相位干涉仪的测向结果：

信号处理机利用

和

分别对方位角θ和俯仰角β进行最小二乘估计，得到入射信号的方位角θ和俯仰角β：

其中S＝min(M,M′)，min表示取最小值函数，d_m表示

对应的基线长度，d_m′表示φ_m′对应的基线长度。

2.根据权利要求1所述的一种适用于阵元失效条件下的相位干涉仪测向方法，其特征在于，步骤(6)中所述的信号处理机对I的基线相位差(φ₁,φ₂,...,φ_m,...,φ_M)逐级进行解相位模糊，以及对II的基线相位差(φ₁′,φ₂′,...,φ_m′,...,φ_M′)逐级进行解相位模糊，公式为：

其中int表示取整函数。

Images