CN111142064B - 一种提高阵列自由度和虚拟孔径的阵列设置方法 - Google Patents
一种提高阵列自由度和虚拟孔径的阵列设置方法 Download PDFInfo
- Publication number
- CN111142064B CN111142064B CN202010010591.6A CN202010010591A CN111142064B CN 111142064 B CN111142064 B CN 111142064B CN 202010010591 A CN202010010591 A CN 202010010591A CN 111142064 B CN111142064 B CN 111142064B
- Authority
- CN
- China
- Prior art keywords
- array
- elements
- array element
- virtual
- subarrays
- 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
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
-
- 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)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
本发明属于阵列信号处理技术领域,具体涉及一种提高阵列自由度和虚拟孔径的阵列设置方法。首先,根据给定的任意阵元数W确定两个线性子阵的阵元数量m和n;其次,根据m和n确定两个线性子阵的阵列间隔dm和dn;最后根据m、n、dm和dn分别确定两个线性子阵的阵元位置。本发明与传统互质阵列相比具有更大的虚拟阵元个数、连续虚拟阵元个数和阵列物理孔径,可以实现任意数阵元的阵列设置,可实现入射信号源数大于阵元数情况下的DOA估计。在相同阵元数情况下,本发明与传统互质阵列相比具有更多的信源到达角估计能力和更低的DOA估计均方根误差。本发明设计简单、合理、可行性高,能有效解决传统互质阵列物理孔径小的实际问题。
Description
技术领域
本发明属于阵列信号处理领域,涉及一种提高阵列自由度和虚拟孔径的阵列设置方法。
背景技术
波达方向(Direction of arrival,DOA)估计是雷达、声呐、地震和无线电通信处理领域内的一个重要研究内容。对于一个线性阵列采用传统的子空间类方法如多重信号分类(Multiple signal classification,MUSIC)和多尺度旋转不变因子(Estimation ofsignal parameters via rotational invariant techniques,ESPRIT)可以实现最多W-1个信号源的DOA估计。对于信号源数大于阵元数的欠定DOA估计问题,通常是利用等效的虚拟阵列提高阵列的自由度(Degrees of freedom,DOFs),该虚拟阵列是通过对阵列接收信号的协方差矩阵向量化来构造的,通常,增大阵列物理孔径、虚拟阵元个数和连续虚拟阵元个数可以实现DOA估计性能的改善。目前主要有最小冗余阵列、最小空洞阵列、嵌套阵列和互质阵列等四种稀疏阵列可以有效提高阵列孔径、虚拟阵元个数和连续虚拟阵元个数。最小冗余阵列和最小空洞阵列由于没有闭式解,因而无法给出任意阵元数的阵元位置;嵌套阵列的第一级子阵为密集阵列,因而容易产生严重的互耦,限制了其应用的范围;互质阵列相对来说阵元间距较大,产生的互耦较小,但是其物理孔径较小、虚拟阵列个数少、连续虚拟阵元个数少,对提高阵列角度分辨率和信源个数不利,且阵列设计的前提是寻找满足条件的互质数,不易给出任意数阵元的阵元位置闭式解。
发明内容
本发明的目的在于提出一种提高阵列自由度和虚拟孔径的阵列设置方法,解决了现有互质阵列虚拟阵元个数少、连续虚拟阵元个数少和物理孔径不够大的问题。
为解决上述技术问题,本发明基于互质数原理提出了一种提高虚拟阵列的任意数阵元的阵列设置方法,可实现任意数阵元的阵列设置,且该阵列具有更高的物理孔径、虚拟阵元个数和连续虚拟阵列个数。
现将本发明构思及技术解决方案叙述如下:
步骤1:根据给定的阵元数W确定改进型互质阵列两个子阵的阵元数量m和n。
步骤2:根据m和n确定两个子阵的阵列间隔dm和dn。
步骤3:根据m、n、dm和dn确定两个子阵的阵元位置Sm和Sn。
步骤4:构成最终的改进型互质阵列,该阵列阵元位置可表示为S=Sm∪Sn。
步骤1所述两子阵列的阵元数量m和n为:
(1)当总阵元数W为3的倍数时,即表示为W=3k,则令R=k,n=W-m,q=0;
(2)当总阵元数W为3的倍数加1时,即表示为W=3k+1,则令R=k,n=W-m,令,q=1;
(3)当总阵元数W为3的倍数加2时,即表示为W=3k+2,则令R=k+1,n=W-m,q=2。
步骤2所述的两个子阵列的阵列间隔dm和dn为:
(1)子阵列1为含有m个阵元且阵元间距为(R+1)d的均匀线阵,其中d为入射信号半波长;
(2)子阵列2为含有n个阵元且阵元间距为Rd的均匀线阵。
步骤3所述的两个子阵的阵元位置Sm和Sn为:
(1)子阵列2为含有n个阵元,其阵元位置表示为:Sn={0,R,…,(n-1)Rd}。
(2)子阵列1为含有m个阵元,其阵元位置表示为:
当q=0或q=1时,Sm={0,(R+1)d,…,(m-1)(R+1)d}+(R(R+q)+1)d。
当q=2时,Sm={0,(R+1)d,…,(m-1)(R+1)d}+(R2+1)d。
本发明与现有互质阵列设置方法相比具有以下优点:
1、本发明提出的阵列设置方法,基于互质数原理,通过合理的阵元设置,同时提升了阵列的虚拟阵元个数、连续虚拟阵元个数和物理孔径,且可以实现任意数阵元的阵列设置。
2、本发明提出的阵列设置方法,提高了阵列的物理孔径,即整个阵列长度更大,因此可以提高测向精度。
3、本发明提出的阵列设置方法,提高了虚拟阵元个数和连续虚拟阵列个数,即可以实现更多信源的DOA估计。
4、本发明设计简单、合理、可行性高,能有效解决传统互质阵列物理孔径小的实际问题。
附图说明
图1、为本发明阵列设置方法的工艺流程图;
图2、为两种阵列物理孔径随阵元数变化图;
图3、为两种阵列的虚拟阵元个数随阵元数变化图;
图4、为两种阵列的连续虚拟阵元个数随阵元数变化图。
具体实施方式
下面结合附图对本发明做更进一步描述:
本发明针对传统互质阵列中存在的阵元物理孔径较小、虚拟阵元个数和连续虚拟阵元个数少的问题,提出了一种提高阵列自由度和虚拟孔径的阵列设置方法。可用于入射信源个数大于阵元数情况下的波达方向估计、相比于传统的互质阵列具有更好的DOA估计性能。
参照图1,即本发明阵列设置方法的工艺流程图,具体实现步骤如下:
步骤1:根据给定的阵元数W确定改进型互质阵列两个子阵的阵元数量m和n。
1a)当总阵元数W为3的倍数时,即表示为W=3k,则令R=k,n=W-m,q=0;
1b)当总阵元数W为3的倍数加1时,即表示为W=3k+1,则令R=k,n=W-m,令,q=1;
1c)当总阵元数W为3的倍数加2时,即表示为W=3k+2,则令R=k+1,n=W-m,q=2。
步骤2:根据m和n确定两个子阵的阵列间隔dm和dn。
2a)子阵列1为含有m个阵元且阵元间距为(R+1)d的均匀线阵,其中d为入射信号半波长;
2b)子阵列2为含有n个阵元且阵元间距为Rd的均匀线阵。
步骤3:根据m、n、dm和dn确定两个子阵的阵元位置Sm和Sn。
3a)子阵列2为含有n个阵元,其阵元位置表示为:Sn={0,R,…,(n-1)Rd}。
3b)子阵列1为含有m个阵元,其阵元位置表示为:
当q=0或q=1时,Sm={0,(R+1)d,…,(m-1)(R+1)d}+(R(R+q)+1)d。
当q=2时,Sm={0,(R+1)d,…,(m-1)(R+1)d}+(R2+1)d。
步骤4:构成最终的改进型互质阵列,该阵列阵元位置可表示为S=Sm∪Sn。
表1:一种新的阵列阵元位置分布表
参照表1,分别是阵元数为30、31和32的阵元位置分布表。带下划线的为子阵1的阵元,不带下划线的为子阵2的阵元。子阵1的阵元间距比子阵2的阵元间距大半个波长。
参照图2,是两种物理阵列孔径随阵元数变化图。从图2中可以显然看到,随着总阵元数的变化,本发明对应的阵列结构的物理阵列孔径大于传统互质阵列。
参照图3和图4,是两种质阵列的虚拟阵元个数和连续虚拟阵元个数随总阵元数变化图。从图3和图4中可以显然看到,随着阵元数的变化,本发明对应的阵列结构的虚拟阵元个数和连续虚拟阵元个数均大于传统互质阵列。
综上,本发明具有更大的虚拟阵元个数、连续虚拟阵元个数和物理阵列孔径,因而具有更优的DOA估计性能。
Claims (4)
1.一种提高阵列自由度和虚拟孔径的阵列设置方法,该方法包括以下步骤:
步骤1:根据给定的阵元数W确定改进型互质阵列两个子阵的阵元数量m和n;
步骤2:根据m和n确定两个子阵的阵列间隔dm和dn;
步骤3:根据m、n、dm和dn确定两个子阵的阵元位置Sm和Sn;
步骤4:构成最终的改进型互质阵列,该阵列阵元位置可表示为S=Sm∪Sn。
2.根据权利要求1所述的一种提高阵列自由度和虚拟孔径的阵列设置方法,其特征在于,步骤1所述两子阵列的阵元数量m和n为:
(1)当总阵元数W为3的倍数时,即表示为W=3k,则令R=k,n=W-m,q=0;
(2)当总阵元数W为3的倍数加1时,即表示为W=3k+1,则令R=k,n=W-m,令,q=1;
(3)当总阵元数W为3的倍数加2时,即表示为W=3k+2,则令R=k+1,n=W-m,q=2。
3.根据权利要求1所述的一种提高阵列自由度和虚拟孔径的阵列设置方法,其特征在于,步骤2所述的两个子阵列的阵列间隔dm和dn为:
(1)子阵列1为含有m个阵元且阵元间距为(R+1)d的均匀线阵,其中d为入射信号半波长;
(2)子阵列2为含有n个阵元且阵元间距为Rd的均匀线阵。
4.根据权利要求1所述的一种提高阵列自由度和虚拟孔径的阵列设置方法,其特征在于,步骤3所述的两个子阵的阵元位置Sm和Sn为:
(1)子阵列2为含有n个阵元,其阵元位置表示为:Sn={0,R,…,(n-1)Rd};
(2)子阵列1为含有m个阵元,其阵元位置表示为:
当q=0或q=1时,Sm={0,(R+1)d,…,(m-1)(R+1)d}+(R(R+q)+1)d;
当q=2时,Sm={0,(R+1)d,…,(m-1)(R+1)d}+(R2+1)d。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010010591.6A CN111142064B (zh) | 2020-01-06 | 2020-01-06 | 一种提高阵列自由度和虚拟孔径的阵列设置方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010010591.6A CN111142064B (zh) | 2020-01-06 | 2020-01-06 | 一种提高阵列自由度和虚拟孔径的阵列设置方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111142064A CN111142064A (zh) | 2020-05-12 |
CN111142064B true CN111142064B (zh) | 2023-09-15 |
Family
ID=70523682
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010010591.6A Active CN111142064B (zh) | 2020-01-06 | 2020-01-06 | 一种提高阵列自由度和虚拟孔径的阵列设置方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111142064B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115865193B (zh) * | 2023-02-27 | 2023-05-09 | 中国人民解放军火箭军工程大学 | 一种用于反射内存组网性能测试的装置及其方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011128088A (ja) * | 2009-12-18 | 2011-06-30 | Keio Gijuku | 位置推定システム及びプログラム |
CN106415931A (zh) * | 2014-05-29 | 2017-02-15 | 丰田自动车株式会社 | 阵列天线装置 |
CN109752687A (zh) * | 2019-01-15 | 2019-05-14 | 哈尔滨工程大学 | 一种降低阵元互耦效应的改进型互质阵列设置方法 |
CN109932680A (zh) * | 2019-04-04 | 2019-06-25 | 哈尔滨工程大学 | 一种基于平移互质阵列的非圆信号波达方向估计方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9912074B2 (en) * | 2014-12-12 | 2018-03-06 | The Boeing Company | Congruent non-uniform antenna arrays |
-
2020
- 2020-01-06 CN CN202010010591.6A patent/CN111142064B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011128088A (ja) * | 2009-12-18 | 2011-06-30 | Keio Gijuku | 位置推定システム及びプログラム |
CN106415931A (zh) * | 2014-05-29 | 2017-02-15 | 丰田自动车株式会社 | 阵列天线装置 |
CN109752687A (zh) * | 2019-01-15 | 2019-05-14 | 哈尔滨工程大学 | 一种降低阵元互耦效应的改进型互质阵列设置方法 |
CN109932680A (zh) * | 2019-04-04 | 2019-06-25 | 哈尔滨工程大学 | 一种基于平移互质阵列的非圆信号波达方向估计方法 |
Non-Patent Citations (1)
Title |
---|
张伟 ; 刘畅 ; 贾勇 ; 石荣 ; .基于双频虚拟互质阵列的均匀稀疏阵列DOA估计算法.现代雷达.2018,(04),全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN111142064A (zh) | 2020-05-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109932680B (zh) | 一种基于平移互质阵列的非圆信号波达方向估计方法 | |
CN108344967B (zh) | 基于互质面阵的二维波达方向快速估计方法 | |
Zheng et al. | DOA estimation for coprime linear arrays: An ambiguity-free method involving full DOFs | |
CN107315160B (zh) | 基于内插虚拟阵列信号原子范数最小化的互质阵列波达方向估计方法 | |
CN108957391B (zh) | 一种基于嵌套阵列的l型天线阵的二维波达方向估计方法 | |
CN107102291B (zh) | 基于虚拟阵列内插的无网格化互质阵列波达方向估计方法 | |
CN107092004B (zh) | 基于信号子空间旋转不变性的互质阵列波达方向估计方法 | |
CN112698264B (zh) | 增广互质阵列脉冲噪声环境下相干信源的doa估计方法 | |
CN109669156B (zh) | 冲击噪声下基于量子帝王蝶的圆阵模式空间动态测向方法 | |
CN111736118B (zh) | 一种线列阵阵列扩展方法 | |
CN107783078B (zh) | 一种波束-多普勒酉esprit多目标角度估计方法 | |
CN107561484A (zh) | 基于内插互质阵列协方差矩阵重建的波达方向估计方法 | |
CN112115595A (zh) | 一种用于非圆信号波达方向估计的三级扩展变换嵌套阵设计方法 | |
CN104931923A (zh) | Grid Iterative ESPRIT,一种可扩展的用于均匀圆阵二维到达角的快速估计算法 | |
CN111142064B (zh) | 一种提高阵列自由度和虚拟孔径的阵列设置方法 | |
CN113189538A (zh) | 一种基于互质稀疏排列的三元阵列及其空间谱估计方法 | |
CN113075610B (zh) | 一种基于互质极化阵列的差分阵列内插的doa估计方法 | |
CN110531312A (zh) | 一种基于稀疏对称阵列的doa估计方法和系统 | |
CN108614235B (zh) | 一种多鸽群信息交互的单快拍测向方法 | |
CN110703185B (zh) | 一种基于多级扩展嵌套阵列的波达方向估计方法 | |
Xiao et al. | A weighted forward-backward spatial smoothing DOA estimation algorithm based on TLS-ESPRIT | |
CN113391261A (zh) | 一种基于泰勒展开的增广互质阵列快速doa估计方法 | |
CN110095749B (zh) | 基于非圆信号的多重孔径嵌套阵列设置及doa估计方法 | |
Yao et al. | Degree-of-freedom strengthened cascade array for DOD-DOA estimation in MIMO array systems | |
CN110726967A (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 |