CN112115595A - 一种用于非圆信号波达方向估计的三级扩展变换嵌套阵设计方法 - Google Patents
一种用于非圆信号波达方向估计的三级扩展变换嵌套阵设计方法 Download PDFInfo
- Publication number
- CN112115595A CN112115595A CN202010927619.2A CN202010927619A CN112115595A CN 112115595 A CN112115595 A CN 112115595A CN 202010927619 A CN202010927619 A CN 202010927619A CN 112115595 A CN112115595 A CN 112115595A
- Authority
- CN
- China
- Prior art keywords
- array
- virtual
- nested
- arrays
- freedom
- 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.)
- Withdrawn
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- 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
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种用于非圆信号波达方向估计的三级扩展变换嵌套阵设计方法。该方法首先交换了传统嵌套阵两个子阵的位置,使稀疏子阵在前,密集子阵在后。之后通过实施对密集子阵的稀疏化布局,使其虚拟和差阵具有大范围连续段,进而获得自由度的提升,其实现步骤是:计算阵元间距基本单元,确定阵元数和阵元间隔参数,根据解析式确定阵元的物理位置。本发明可以在相同阵元数的条件下,有效增大虚拟和差阵的自由度,提高波达方向的估计精度,提高空间目标探测能力,减少阵列运行成本。
Description
技术领域
本发明属于阵列信号处理领域,涉及到一种适用于非圆信号波达方向估计的稀疏阵列优化布局及设计方法,具体是一种基于嵌套阵的阵列孔径扩展结构,可用于生成高自由度的虚拟和差阵列。
背景技术
波达方向估计是阵列信号处理领域的一个重要课题,在雷达、声呐、无线通信、导航等应用中扮演了重要的角色。传统的波达方向估计采用了均匀排布的传感器阵列,传感器间隔小于半波长。这种均匀排布方式,保证了信号的空间谱形式唯一,空间信息无混叠,并能通过谱函数唯一辨识空间探测目标。然而,均匀线阵的可识别信源数受到了自由度的严重限制,一个N传感器的均匀线阵最多只能识别N-1个信源。
为了解决信源数多于阵元数的欠定波达方向估计问题,使用了较少传感器的非均匀线阵,也被称作稀疏阵,被提出。稀疏阵通过对各阵元接收信号计算协方差,再对协方差矩阵进行矢量化操作,获得等效虚拟相控阵。其虚拟阵元位置位于实际阵元的差分位置处,如果能令虚拟差阵连续排布,就可以利用虚拟阵实现目标方位估计。然而,之前的研究都是基于圆信号发展而来的。事实上,大量的非圆信号,如二进制相移键控信号和脉冲振幅调制信号,都被广泛应用于无线通信系统中。非圆信号的一个重要特性是它的椭圆协方差矩阵不等于零。对椭圆协方差矩阵进行矢量化操作,得到的虚拟阵阵元位置是实际阵元位置之和。因此,结合收到非圆信号的协方差与椭圆协方差矩阵,得到的虚拟阵将同时包含了差阵与和阵,虚拟阵元数相较于传统差阵能够提升一倍。对于非圆信号的欠定波达方向估计问题,一些新的稀疏阵被提出。例如,有和差虚拟阵的改进嵌套阵通过同时移位传统嵌套阵的稀疏子阵与密集子阵扩大了虚拟阵孔径,但是这个阵型需要在零点布置一个额外传感器作为参考传感器。变换嵌套阵基于降低和阵与差阵之间冗余的思想,交换了传统嵌套阵两个子阵的位置,进一步提高了阵列自由度。但是,目前提出的这些稀疏阵列结构,其自由度仍有很大的提升空间。因此,如何设计一种一维稀疏阵,使其虚拟和差阵具有高自由度(即更长连续段),对于实际应用有着重要的意义。
发明内容
本发明的目的在于提供一种适用于非圆信号波达方向估计的稀疏阵列优化布局及设计方法,在变换嵌套阵的基础上,对其密集子阵进行稀疏化处理。在具备相同数目的传感器情况下,获得比现有稀疏阵更大的虚拟和差阵、更高的阵列自由度和更强的空间目标探测能力。
为了解决上述技术问题,本发明是通过以下技术方案实现的,它包括步骤如下:
(1)用d表示阵元间距的基本单元,将其取为半波长,即d=λ/2;λ表示阵列入射信号的波长;
(2)根据阵列给定的总阵元数N,选取两个正整数N1和N2作为阵列参数,满足N=N1+N2,其中N1≥13,N2≥1;
(7)找到虚拟和差阵的最大连续段[-Lu,Lu]d,即可得到阵列自由度Lu;利用这个虚拟连续均匀线阵,可实施各类波达方向估计算法,对空间信号的来波方向进行准确估计。
本发明与现有技术相比具有以下优点:
(1)本发明通过将变换嵌套阵的密集子阵稀疏化,即将密集子阵分成了三部分,并将其中两部分分别向右移动(N1+1)d和2(N1+1)d,降低了密集子阵的自差集与密集和稀疏子阵的互差集之间的冗余以及密集子阵的自和集与密集和稀疏子阵的互和集之间的冗余;本发明可以扩大阵列孔径并有效降低天线之间的互藕效应;
(2)相同阵元数条件下,本发明提出的三级扩展变换嵌套阵结构可以通过虚拟和差阵获得更大的虚拟连续孔径和阵列自由度,有效提高波达方向估计精度,增大空间分辨力,提高空间目标探测数量,减少阵列运行成本;
(3)根据阵元总数N就可以获得具有最大自由度的阵列排布位置,方便理论设计和实际工程应用。
附图说明
图2是本发明提出的三级扩展变换嵌套阵的虚拟和差阵示意图,其中N1=13,N2=3。
图3是本发明设计结构与其他四种稀疏结构的空间目标方位估计结果对比图。
图4(a)是本发明设计结构与其他四种稀疏结构的归一化波达方向估计均方误差随信噪比变化对比图。
图4(b)是本发明设计结构与其他四种稀疏结构的归一化波达方向估计均方误差随框架数量变化对比图。
具体实施方式
下面结合附图及实施例对本发明的内容做进一步的说明。
以入射信号频率20GHz,阵元总数N=16为例,对本发明的三级扩展变换嵌套阵(如图1所示)构建进行详细说明。
(1)计算入射信号波长λ=c/f=3.0×108/20×109=1.5cm;用d表示阵元间距的基本单元,将其取为半波长,即d=λ/2=7.5mm;
(2)根据阵列给定的总阵元数N=16,选择一组N1和N2,满足N=N1+N2,其中N1≥13,N2≥1;在本例中选择了N1=13,N2=3;
实际阵元布放位置为
{0,105,210,225,232.5,240,300,322.5,360,367.5,375,382.5,390,397.5,412.5,457.5}mm;
(6)计算三级扩展变换嵌套阵的虚拟和差阵(如图2所示)。由于虚拟和差阵关于零点对称,图2中只画出了虚拟和差阵的正半部分;
下面结合仿真实例对本发明的效果做进一步描述。
仿真实例:设定阵元总数N=18不变,对以下五种阵型结构进行比较。基本互质阵(Coprime Array,CA)选取N1=5,N2=9,有多周期子阵的和差互质阵(Diff–sum CoprimeArray with Multi-period Subarrays,DsCAMpS)选取P2=4,N1=4,N2=3,第一类变换嵌套阵(Transformed Nested Array-1,TNA-1)选取N1=13,N2=5,第二类变换嵌套阵(Transformed Nested Array-2,TNA-2)选取N1=13,N2=5,本发明的三级扩展变换嵌套阵(Three-level Extended Transformed Nested Array,ThETNA)选取N1=13,N2=5.根据上述嵌套阵的解析表达式,可以获得它们的归一化阵元位置(省略间距单元d)分别为
假设空间有76个不相关非圆信源入射到上述稀疏阵上,入射角均匀分布在-60°和60°之间。入射信号的信噪比为5dB,框架数为F=200,框架长度为L=200。五种稀疏阵对波达方向的估计结果如图3所示。可以看到,本发明提出的三级扩展变换嵌套阵ThETNA能够准确估计出所有76个信源,TNA-1和TNA-2的估计结果在[-60°,-50°]和[50°,60°]之间出现了明显的错误,而两类互质阵无法估计出正确的入射角度。这是由于,CA、DsCAMpS、TNA-1、TNA-2、ThETNA的自由度分别是94、90、138、152、170。自由度越大,估计准确性越高。因此本发明提出的阵型估计结果最准确。
假设51个不相关非圆信号入射到上述五种稀疏结构上,入射角度均匀分布在-60°和60°之间。实施500次蒙特卡洛实验,框架数与框架长度分别为F=200和L=200。归一化波达方向估计的均方误差(Mean Square Error,MSE)随入射信号信噪比(Signal to noiseratio,SNR)变化的曲线绘制在图4(a)中。图4(b)中给出了归一化波达方向估计的均方误差作为框架数的函数,这里信噪比被设为0dB,框架长度为L=100。由图可见,本发明提出的阵型的均方误差比其他四种稀疏阵都低,由于它的自由度最高。因此提出阵型的目标估计精度最高。
综上所述,本发明所提出的阵列结构适用于非圆信号的波达方向估计,它的虚拟和差阵相较于现有稀疏阵有更长的连续范围,因此能提供更高的自由度。该结构在阵列波达方向估计方面和空间目标探测能力方面具有更优越的性能。
Claims (4)
1.用于一维阵列测向的三级扩展变换嵌套阵设计方法,其特征在于:包括以下步骤:
(1)用d表示阵元间距的基本单元,将其取为半波长,即d=λ/2;λ表示阵列入射信号的波长;
(2)根据阵列给定的总阵元数N,选取两个正整数N1和N2作为阵列参数,满足N=N1+N2,其中N1≥13,N2≥1;
(7)找到虚拟和差阵的最大连续段[-Lu,Lu]d,即可得到阵列自由度Lu;利用这个虚拟连续均匀线阵,可实施各类波达方向估计算法,对空间信号的来波方向进行准确估计。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010927619.2A CN112115595A (zh) | 2020-09-07 | 2020-09-07 | 一种用于非圆信号波达方向估计的三级扩展变换嵌套阵设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010927619.2A CN112115595A (zh) | 2020-09-07 | 2020-09-07 | 一种用于非圆信号波达方向估计的三级扩展变换嵌套阵设计方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112115595A true CN112115595A (zh) | 2020-12-22 |
Family
ID=73802709
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010927619.2A Withdrawn CN112115595A (zh) | 2020-09-07 | 2020-09-07 | 一种用于非圆信号波达方向估计的三级扩展变换嵌套阵设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112115595A (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113010992A (zh) * | 2020-12-30 | 2021-06-22 | 中国航天科工集团八五一一研究所 | 一种基于展宽方式的嵌套阵列流形优化方法 |
CN113050026A (zh) * | 2021-03-29 | 2021-06-29 | 西安电子科技大学 | 拓展阵元间距移位子阵的结构性运动稀疏阵列设计方法 |
CN113075611A (zh) * | 2021-04-02 | 2021-07-06 | 西北大学 | 基于扩展和移位嵌套阵列构建四阶差分阵列的方法 |
CN113296049A (zh) * | 2021-04-28 | 2021-08-24 | 南京航空航天大学 | 互质阵列脉冲环境下非圆信号的共轭增广doa估计方法 |
CN114859301A (zh) * | 2022-03-31 | 2022-08-05 | 西安电子科技大学 | 用于波达方向估计的运动平台拓展空间分形阵列设计方法 |
CN114996922A (zh) * | 2022-05-23 | 2022-09-02 | 南京航空航天大学 | 一种稀疏嵌套阵列及其设计方法 |
CN112986900B (zh) * | 2021-02-07 | 2023-09-19 | 中国人民解放军国防科技大学 | 一种用于doa估计的灵活互质阵列 |
-
2020
- 2020-09-07 CN CN202010927619.2A patent/CN112115595A/zh not_active Withdrawn
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113010992A (zh) * | 2020-12-30 | 2021-06-22 | 中国航天科工集团八五一一研究所 | 一种基于展宽方式的嵌套阵列流形优化方法 |
CN112986900B (zh) * | 2021-02-07 | 2023-09-19 | 中国人民解放军国防科技大学 | 一种用于doa估计的灵活互质阵列 |
CN113050026A (zh) * | 2021-03-29 | 2021-06-29 | 西安电子科技大学 | 拓展阵元间距移位子阵的结构性运动稀疏阵列设计方法 |
CN113075611A (zh) * | 2021-04-02 | 2021-07-06 | 西北大学 | 基于扩展和移位嵌套阵列构建四阶差分阵列的方法 |
CN113075611B (zh) * | 2021-04-02 | 2023-09-08 | 西北大学 | 基于扩展和移位嵌套阵列构建四阶差分阵列的方法 |
CN113296049A (zh) * | 2021-04-28 | 2021-08-24 | 南京航空航天大学 | 互质阵列脉冲环境下非圆信号的共轭增广doa估计方法 |
CN113296049B (zh) * | 2021-04-28 | 2024-02-20 | 南京航空航天大学 | 互质阵列脉冲环境下非圆信号的共轭增广doa估计方法 |
CN114859301A (zh) * | 2022-03-31 | 2022-08-05 | 西安电子科技大学 | 用于波达方向估计的运动平台拓展空间分形阵列设计方法 |
CN114996922A (zh) * | 2022-05-23 | 2022-09-02 | 南京航空航天大学 | 一种稀疏嵌套阵列及其设计方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112115595A (zh) | 一种用于非圆信号波达方向估计的三级扩展变换嵌套阵设计方法 | |
CN107102291B (zh) | 基于虚拟阵列内插的无网格化互质阵列波达方向估计方法 | |
CN108344967B (zh) | 基于互质面阵的二维波达方向快速估计方法 | |
CN107315160B (zh) | 基于内插虚拟阵列信号原子范数最小化的互质阵列波达方向估计方法 | |
CN108896954B (zh) | 互质阵中一种基于联合实值子空间的波达角估计方法 | |
CN109490820B (zh) | 一种基于平行嵌套阵的二维doa估计方法 | |
CN110927661A (zh) | 基于music算法的单基地展开互质阵列mimo雷达doa估计方法 | |
CN112014792B (zh) | 一种改进的二维互质阵列doa估计方法 | |
CN107092004A (zh) | 基于信号子空间旋转不变性的互质阵列波达方向估计方法 | |
CN107703478B (zh) | 基于互相关矩阵的扩展孔径二维doa估计方法 | |
CN104931923A (zh) | Grid Iterative ESPRIT,一种可扩展的用于均匀圆阵二维到达角的快速估计算法 | |
CN109613504A (zh) | 一种稀疏线性阵列的快速角度估计方法 | |
CN113189538A (zh) | 一种基于互质稀疏排列的三元阵列及其空间谱估计方法 | |
Patwari | Sparse linear antenna arrays: a review | |
CN110726967B (zh) | 用于一维阵列测向的单边稀疏嵌套阵设计方法 | |
CN116224219A (zh) | 一种阵列误差自校正原子范数最小化doa估计方法 | |
CN112100837A (zh) | 一种用于非圆信号波达方向估计的二级扩展变换嵌套阵设计方法 | |
CN115728722A (zh) | 4d雷达的天线阵列、数据探测方法及4d雷达 | |
CN110095749B (zh) | 基于非圆信号的多重孔径嵌套阵列设置及doa估计方法 | |
CN113050026A (zh) | 拓展阵元间距移位子阵的结构性运动稀疏阵列设计方法 | |
CN114609580A (zh) | 一种基于非圆信号的无孔互质阵列设计方法 | |
CN109507634A (zh) | 一种任意传感器阵列下的基于传感算子的盲远场信号波达方向估计方法 | |
CN111830458B (zh) | 一种平行线阵单快拍二维测向方法 | |
CN111142064B (zh) | 一种提高阵列自由度和虚拟孔径的阵列设置方法 | |
CN111650552A (zh) | 一种基于和差协同阵构建的改进l型阵列设计方法 |
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 | ||
WW01 | Invention patent application withdrawn after publication |
Application publication date: 20201222 |
|
WW01 | Invention patent application withdrawn after publication |