CN112858995A - 一种基于分布式阵列的联合测角方法和装置 - Google Patents
一种基于分布式阵列的联合测角方法和装置 Download PDFInfo
- Publication number
- CN112858995A CN112858995A CN202110081253.6A CN202110081253A CN112858995A CN 112858995 A CN112858995 A CN 112858995A CN 202110081253 A CN202110081253 A CN 202110081253A CN 112858995 A CN112858995 A CN 112858995A
- Authority
- CN
- China
- Prior art keywords
- array
- matrix
- sub
- phase difference
- arrays
- 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.)
- Granted
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
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
本发明提出一种基于分布式阵列的联合测角方法和装置,属于天线阵列测角技术领域。该方法首先将多个均匀线性子阵列水平摆放在同一直线上构成分布式阵列;然后在任一采样时刻利用该分布式阵列对多个目标辐射的电磁波信号进行采样得到单快拍采样信号,根据盲源分离算法估计各子阵列间的相位差,最后通过匹配滤波得到该采样时刻各目标方位角的估计结果。本发明可以实现在各个子阵列空间位置关系未知的情况下,使得该分布式阵列的角度分辨能力接近各子阵列彼此相参时的角度分辨能力,即等效于单个大孔径阵列,并实现精准测角。
Description
技术领域
本发明属于天线阵列测角技术领域,特别提出一种基于分布式阵列的联合测角方法和装置。
背景技术:
阵列测角是一种经典的目标方位角测量方法,该方法利用目标到阵列各阵元的距离差与各阵元接收信号间相位差的对应关系,将角度测量问题转化为阵列接收信号的频率估计问题,从而通过频率估计方法实现测角。
在对多目标同时进行角度测量时,阵列测角方法的角度分辨力取决于阵列的孔径,即对于单个阵列,阵列的孔径越大,其角度分辨能力越高。实际应用中,为实现更高的角度分辨力,阵列孔径设计已可以达到百米量级。但对于这样的阵列系统,由于规模、尺寸庞大,一方面导致系统设计上相比小规模阵列更为复杂,另一方面导致其难以应用到运动平台,从而给广泛使用带来困难。
为解决单个大孔径阵列存在的问题,一种方式是采用分布式相参阵列作为代替。在角度测量的场景中,分布式相参阵列利用空间分布临近的多个小孔径子阵列,在各子阵列空间位置关系精确已知的情况下,实现等效于各子阵列孔径之和的大孔径阵列的角度分辨能力,同时避免了单个大孔径阵列系统复杂度高的问题。但实际情况中,特别是在运动平台,各子阵列空间位置关系很难精确已知,这给分布式相参阵列系统的应用带来挑战。
因此,如何在各子阵列空间位置关系未知或不精确已知的情况下,实现分布式阵列联合测角,并达到接近相参或单个大孔径阵列的角度分辨力,这成为目前阵列测角的一个关键问题。已有的算法尝试从子空间投影、低秩矩阵恢复、稀疏恢复自校正等角度去解决这个问题,但这些方法都存在相应的问题。如子空间投影和低秩矩阵恢复等方法,其在各子阵列空间位置关系未知的情况下,直接对目标方位角进行估计,这类方法通常不能实现各子阵列空间位置关系的精确估计,从而无法达到接近相参阵列的角度分辨力;而稀疏恢复自校正等方法需假设不同目标方位角非常接近,在实际情况中,这样的假设会带来模型失配,导致估计存在偏差。因此,目前相关方法都不能很好地解决分布式阵列联合测角的问题。
盲源分离技术通过利用不同信号源辐射电磁波波形的独立性,实现在阵列流形未知的情况下,对不同源信号的波形进行分离和估计。盲源分离技术目前在通信信号分离、雷达抗干扰等领域中有广泛的应用。
发明内容:
本发明的目的是为克服已有技术的不足之处,提出一种基于分布式阵列的联合测角方法和装置。本发明可以实现在各个子阵列空间位置关系未知的情况下,使得该分布式阵列的角度分辨能力接近各子阵列彼此相参时的角度分辨能力,即等效于单个大孔径阵列,并实现精准测角。
本发明提出一种基于分布式阵列的联合测角方法,其特征在于,包括以下步骤:
1)设置分布式阵列;
将M≥2个均匀线性子阵列水平摆放在同一直线上,该M个子阵列组成分布式阵列;其中每个子阵列有N个阵元,N≥2;各子阵列内相邻阵元的间距为d;
在该阵列的远场,有U个目标持续辐射波长为ξ的电磁波,U≥2,各目标相对阵列的方位角表示为θu∈(-π/2,π/2),u=1,...,U;
2)数据采集;
将t时刻第m个子阵列中第n个阵元接收到的信号表示为ym,n(t),m=1,...,M,n=1,...,N,则t时刻第m个子阵列接收信号构成矢量ym(t),表示为:
将各子阵列接收的信号在Ts>0时刻进行采样,得到该时刻的单快拍采样信号构成矩阵Y,表示为:
3)利用步骤2)得到的矩阵Y,采用盲源分离算法估计各子阵列间的相位差,得到Ts时刻的子阵列间相位差估计矩阵;
4)利用步骤3)得到的子阵列间相位差估计矩阵,采用匹配滤波的方式得到Ts时刻各目标方位角的估计结果。
本发明的特点及有益效果:
本发明提出一种基于分布式阵列的联合测角方法和装置,通过挖掘所有子阵列的单快拍数据中不同目标对应信号之间的独立性,实现对各个子阵列间相位差的估计和补偿,从而使分布式阵列系统在各子阵列空间位置关系未知的情况下达到接近相参阵列的角度分辨力。本发明创新性地将盲源分离技术应用到基于分布式阵列的联合测角问题中,有效地提升系统角度分辨能力。本发明的分布式阵列的角度分辨能力等效于单个大孔径阵列,测角精准,易于实现。
附图说明
图1是本发明方法的整体流程图;
图2是本发明装置的结构示意图。
具体实施方式:
本发明提出一种基于分布式阵列的联合测角方法和装置,下面结合附图和具体实施例对本发明进行进一步详细说明。
本发明提出一种基于分布式阵列的联合测角方法,整体流程图如图1所示,包括以下步骤:
1)设置分布式阵列;
本发明考虑在一个二维平面上,将M≥2个临近的均匀线性子阵列水平摆放在同一直线上,该M个子阵列组成分布式阵列;其中每个子阵列有N个阵元,N≥2,用d表示子阵列内相邻阵元的间距,用ηm表示第m个子阵中第1个阵元相对第1个子阵列中第1个阵元的相对空间位置,m=1,...,M,其中,η1=0。在阵列的远场,有U个目标持续辐射窄带信号的电磁波,本发明将所有目标辐射的电磁波的波长记为ξ,U≥2,各目标相对阵列的方位角表示为θu∈(-π/2,π/2),u=1,...,U。场景中,已知参数M,N,d,ξ,U,未知参数ηm,θu,待估计参数为θu。
2)数据采集;
将t时刻第m个子阵列中第n个阵元接收到的信号表示为ym,n(t),m=1,...,M,n=1,...,N,则t时刻第m个子阵列接收信号构成矢量ym(t),表示为:
将各子阵列接收的信号在Ts>0时刻进行采样,得到该时刻的单快拍采样信号构成矩阵Y,表示为:
3)相位差估计;
采用盲源分离算法估计各子阵列间的相位差。本发明以盲源分离中JointApproximate Diagonalization of Eigen-matrices(JADE)算法为例,具体步骤如下:
3-1)数据白化:首先,利用单快拍采样信号构成的矩阵Y计算协方差矩阵,即:
构建白化矩阵W,表示为:
将采样信号Y进行白化处理,得到:
[QZ]a,b表示矩阵QZ的第a行第b列元素,a=i+(j-1)U,b=k+(l-1)U,i,j,k,l=1,...,U;
而Cum(·)定义为:
其中,E代表期望;
在本发明考虑的场景中,可以用多个子阵列的单快拍数据估计上述期望,即:
3-3)对Z的四阶积累矩阵QZ进行特征值分解,得到前U个最大的特征值和对应的特征向量,分别记为λu>0和定义一组矩阵其中每个矩阵满足vec(Hu)=hu,其中vec(·)表示对矩阵做列向量化。建立矩阵集合
其中,|diag(·)|2表示求矩阵的对角元素平方和。矩阵近似联合对角化是一个经典的矩阵优化问题,而上述优化问题可以采用Jacobi方法进行求解。针对U=2的情况,具体方法为:
gu=[[λuHu]1,1-[λuHu]2,2,[λuHu]1,2+[λuHu]2,1,j([λuH]2,1-[λuHu]1,2)]T,
4)角度估计;
在实现各个子阵列间的相位差的估计后,采用匹配滤波的方式对Ts时刻各目标方位角进行估计,得到该时刻各目标方位角的估计结果,具体表示为:
本发明还提出一种采用上述方法的基于分布式阵列的联合测角装置,其结构如图2所示,包括:M个子阵列信号采集模块、数据中心存储模块、相位差估计模块和角度估计模块。其中,每个子阵列信号采集模块的输出端分别连接数据中心存储模块的输入端,数据中心存储模块的输出端通过物理链路分别连接到相位差估计模块的输入端与角度估计模块的输入端,相位差估计模块的输出端通过物理链路连接到角度估计模块的输入端。
所述每个子阵列信号采集模块用于分别获取各子阵列接收到的目标辐射电磁波信号,并在同一时刻进行采样,然后将各子阵列采样后的单快拍采样信号数据发送给数据中心存储模块。
所述数据中心存储模块用于存储各子阵列信号采集模块发送的单快拍采样信号数据,然后将存储的数据分别发送给相位差估计模块和角度估计模块。
所述相位差估计模块用于根据数据中心存储模块发送的单快拍采样信号数据,通过盲源分离的算法思想,估计各子阵列之间由于空间分布带来的相位差,然后将子阵列间相位差的估计值发送给角度估计模块。
所述角度估计模块用于根据数据中心存储模块发送的单快拍采样信号数据和相位差估计模块发送的子阵列间相位差的估计值,利用匹配滤波的思想,计算得到对应采样时刻的各个目标的方位角的估计结果并输出。
Claims (4)
1.一种基于分布式阵列的联合测角方法,其特征在于,包括以下步骤:
1)设置分布式阵列;
将M≥2个均匀线性子阵列水平摆放在同一直线上,该M个子阵列组成分布式阵列;其中每个子阵列有N个阵元,N≥2;各子阵列内相邻阵元的间距为d;
在该阵列的远场,有U个目标持续辐射波长为ξ的电磁波,U≥2,各目标相对阵列的方位角表示为θu∈(-π/2,π/2),u=1,...,U;
2)数据采集;
将t时刻第m个子阵列中第n个阵元接收到的信号表示为ym,n(t),m=1,...,M,n=1,...,N,则t时刻第m个子阵列接收信号构成矢量ym(t),表示为:
将各子阵列接收的信号在Ts>0时刻进行采样,得到该时刻的单快拍采样信号构成矩阵Y,表示为:
3)利用步骤2)得到的矩阵Y,采用盲源分离算法估计各子阵列间的相位差,得到Ts时刻的子阵列间相位差估计矩阵;
4)利用步骤3)得到的子阵列间相位差估计矩阵,采用匹配滤波的方式得到Ts时刻各目标方位角的估计结果。
2.如权利要求1所述的方法,其特征在于,所述步骤3)具体步骤如下:
3-1)利用矩阵Y计算协方差矩阵:
构建白化矩阵W,表示为:
将矩阵Y进行白化处理,得到:
[QZ]a,b表示矩阵QZ的第a行第b列元素,a=i+(j-1)U,b=k+(l-1)U,i,j,k,l=1,...,U;
Cum(·)定义为:
其中,E代表期望;
则
其中,|diag(·)|2表示求矩阵的对角元素平方和;
4.一种采用如权利要求3所述方法的基于分布式阵列的联合测角装置,其特征在于,该装置包括:M个子阵列信号采集模块、数据中心存储模块、相位差估计模块和角度估计模块;其中,每个子阵列信号采集模块的输出端分别连接数据中心存储模块的输入端,数据中心存储模块的输出端分别连接相位差估计模块的输入端与角度估计模块的输入端,相位差估计模块的输出端连接角度估计模块的输入端;
所述每个子阵列信号采集模块用于分别获取各子阵列接收到的目标辐射电磁波信号,并在同一时刻进行采样,然后将各子阵列采样后的单快拍采样信号数据发送给数据中心存储模块;
所述数据中心存储模块用于存储各子阵列信号采集模块发送的单快拍采样信号数据,然后将存储的数据分别发送给相位差估计模块和角度估计模块;
所述相位差估计模块用于根据数据中心存储模块发送的单快拍采样信号数据,估计各子阵列之间的相位差,然后将子阵列间相位差的估计值发送给角度估计模块;
所述角度估计模块用于根据数据中心存储模块发送的单快拍采样信号数据和相位差估计模块发送的子阵列间相位差的估计值,计算得到对应采样时刻的各个目标的方位角的估计结果并输出。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110081253.6A CN112858995B (zh) | 2021-01-21 | 2021-01-21 | 一种基于分布式阵列的联合测角方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110081253.6A CN112858995B (zh) | 2021-01-21 | 2021-01-21 | 一种基于分布式阵列的联合测角方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112858995A true CN112858995A (zh) | 2021-05-28 |
CN112858995B CN112858995B (zh) | 2023-01-31 |
Family
ID=76008732
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110081253.6A Active CN112858995B (zh) | 2021-01-21 | 2021-01-21 | 一种基于分布式阵列的联合测角方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112858995B (zh) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5383164A (en) * | 1993-06-10 | 1995-01-17 | The Salk Institute For Biological Studies | Adaptive system for broadband multisignal discrimination in a channel with reverberation |
US20160011297A1 (en) * | 2014-07-11 | 2016-01-14 | Electronics And Telecommunications Research Institute | Positioning method and apparatus using wireless signal |
CN105548952A (zh) * | 2015-12-03 | 2016-05-04 | 中国电子科技集团公司第五十四研究所 | 基于盲源分离的同频多信号测向方法 |
CN108287333A (zh) * | 2018-03-19 | 2018-07-17 | 电子科技大学 | 一种联合jade和clean的主瓣抗干扰方法 |
CN109061597A (zh) * | 2018-08-23 | 2018-12-21 | 哈尔滨工业大学 | 基于盲源分离与时频脊波域滤波的电离层杂波抑制方法 |
CN109471063A (zh) * | 2018-11-06 | 2019-03-15 | 江西师范大学 | 基于延迟快拍的均匀线列阵高分辨波达方向估计方法 |
CN111044979A (zh) * | 2019-12-13 | 2020-04-21 | 电子科技大学 | 一种基于盲源分离的主瓣干扰对消与目标角度估计方法 |
CN111273218A (zh) * | 2020-03-09 | 2020-06-12 | 上海无线电设备研究所 | 一种基于多层互质阵列的相干源波达方向估计方法 |
-
2021
- 2021-01-21 CN CN202110081253.6A patent/CN112858995B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5383164A (en) * | 1993-06-10 | 1995-01-17 | The Salk Institute For Biological Studies | Adaptive system for broadband multisignal discrimination in a channel with reverberation |
US20160011297A1 (en) * | 2014-07-11 | 2016-01-14 | Electronics And Telecommunications Research Institute | Positioning method and apparatus using wireless signal |
CN105548952A (zh) * | 2015-12-03 | 2016-05-04 | 中国电子科技集团公司第五十四研究所 | 基于盲源分离的同频多信号测向方法 |
CN108287333A (zh) * | 2018-03-19 | 2018-07-17 | 电子科技大学 | 一种联合jade和clean的主瓣抗干扰方法 |
CN109061597A (zh) * | 2018-08-23 | 2018-12-21 | 哈尔滨工业大学 | 基于盲源分离与时频脊波域滤波的电离层杂波抑制方法 |
CN109471063A (zh) * | 2018-11-06 | 2019-03-15 | 江西师范大学 | 基于延迟快拍的均匀线列阵高分辨波达方向估计方法 |
CN111044979A (zh) * | 2019-12-13 | 2020-04-21 | 电子科技大学 | 一种基于盲源分离的主瓣干扰对消与目标角度估计方法 |
CN111273218A (zh) * | 2020-03-09 | 2020-06-12 | 上海无线电设备研究所 | 一种基于多层互质阵列的相干源波达方向估计方法 |
Non-Patent Citations (5)
Title |
---|
周必雷等: "基于盲分离的空时联合处理抗复合干扰方法", 《系统工程与电子技术》 * |
朱伟等: "强相干干扰下弱信号波达方向估计", 《雷达科学与技术》 * |
杨雪亚等: "一种新的二维角度估计的高分辨算法", 《电子与信息学报》 * |
汤小为等: "集中式多输入多输出雷达信号盲分离算法研究", 《宇航学报》 * |
王俊等: "单快拍数据的分布式二维阵列测角方法研究", 《电子与信息学报》 * |
Also Published As
Publication number | Publication date |
---|---|
CN112858995B (zh) | 2023-01-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110208735B (zh) | 一种基于稀疏贝叶斯学习的相干信号doa估计方法 | |
WO2018094565A1 (zh) | 脉冲噪声下的波束成形方法及装置 | |
CN109212526B (zh) | 用于高频地波雷达的分布式阵列目标角度测量方法 | |
CN106680815B (zh) | 基于张量稀疏表示的mimo雷达成像方法 | |
CN108683619B (zh) | 一种低复杂度的大规模mimo信道参数估计方法 | |
CN108957387B (zh) | 一种卫星信号二维到达角估计方法及系统 | |
CN112269172B (zh) | 一种基于张量结构的嵌套mimo雷达角度估计方法和装置 | |
CN110515033B (zh) | 一种基于Toeplitz矩阵恢复的欠通道测向系统及方法 | |
CN112379327A (zh) | 一种基于秩损估计的二维doa估计与互耦校正方法 | |
CN113673317B (zh) | 基于原子范数最小化可降维的二维离格doa估计方法 | |
CN110927711A (zh) | 一种双基地emvs-mimo雷达的高精度定位算法及装置 | |
CN112305495A (zh) | 一种基于原子范数最小的互质阵列协方差矩阵重构方法 | |
CN111965591A (zh) | 一种基于四阶累积量矢量化dft的测向估计方法 | |
CN111352063A (zh) | 一种均匀面阵中基于多项式求根的二维测向估计方法 | |
CN106680779B (zh) | 脉冲噪声下的波束成形方法及装置 | |
CN109946663B (zh) | 一种线性复杂度的Massive MIMO目标空间方位估计方法和装置 | |
CN113671485B (zh) | 基于admm的米波面阵雷达二维doa估计方法 | |
CN108614234B (zh) | 基于多采样快拍互质阵列接收信号快速傅里叶逆变换的波达方向估计方法 | |
CN113075649B (zh) | 一种适用于分布式网络化雷达的信号级直接定位方法 | |
CN112711000B (zh) | 基于最小化准则的电磁矢量互质面阵张量功率谱估计方法 | |
CN113567913A (zh) | 基于迭代重加权可降维的二维平面doa估计方法 | |
CN112763972A (zh) | 基于稀疏表示的双平行线阵二维doa估计方法及计算设备 | |
CN112858995B (zh) | 一种基于分布式阵列的联合测角方法和装置 | |
CN109557503B (zh) | 基于相关矩阵重建解相干的mimo互质阵列doa估计方法 | |
WO2022166477A1 (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 |