CN105824002B - 基于嵌套式子阵阵列的波达方向估计方法 - Google Patents

基于嵌套式子阵阵列的波达方向估计方法 Download PDF

Info

Publication number
CN105824002B
CN105824002B CN201610236372.3A CN201610236372A CN105824002B CN 105824002 B CN105824002 B CN 105824002B CN 201610236372 A CN201610236372 A CN 201610236372A CN 105824002 B CN105824002 B CN 105824002B
Authority
CN
China
Prior art keywords
array
submatrix
msub
mrow
nested type
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
Application number
CN201610236372.3A
Other languages
English (en)
Other versions
CN105824002A (zh
Inventor
杨明磊
陈伯孝
孙磊
曾小路
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xidian University
Original Assignee
Xidian University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Xidian University filed Critical Xidian University
Priority to CN201610236372.3A priority Critical patent/CN105824002B/zh
Publication of CN105824002A publication Critical patent/CN105824002A/zh
Application granted granted Critical
Publication of CN105824002B publication Critical patent/CN105824002B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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

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)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种基于嵌套式子阵阵列的波达方向估计方法，主要解决现有技术中存在自由度、阵列孔径和阵列密度不高的问题。其实现过程是：1.给定总阵元数确定子阵数目和子阵内阵元数目；2.根据子阵内阵元数目选取均匀线阵或最小冗余阵列或嵌套式线阵结构；3.根据子阵数目选取均匀线阵或最小冗余阵列或嵌套式线阵结构；4.根据选取的子阵内阵元结构和子阵结构构造嵌套式子阵阵列；5.根据嵌套式子阵阵列得到接收数据X(t)；6.根据X(t)得到差分合成阵列接收数据zc，进而得到秩恢复数据协方差矩阵RSS；7.对RSS特征值分解，得到波达方向估计角度。本发明在同等条件下具有阵列配置灵活，波达方向估计测角性能好的优点，可用于雷达目标信号检测或功率估值。

Description

基于嵌套式子阵阵列的波达方向估计方法
本发明属于信号处理技术领域，特别涉及一种非等距线阵配置方法，用于估计目标个数大于阵元数情况下的波达方向，提高阵列的自由度、测角精度和密度。
背景技术
波达方向DOA估计是雷达、声纳信号处理中的一个重要研究方向。对于一个具有S个阵元的均匀线性阵列，采用传统的波达方向估计方法，如MUSIC，ESPRIT等所能分辨的最大目标个数为S-1个。而欠定波达方向估计问题，即目标个数大于阵元数的DOA估计经常出现并引起广泛的研究兴趣。解决该问题的一种有效方法是利用一个等效的虚拟阵列来提高波达方向估计的自由度DOF，该虚拟阵列是通过对一个特殊设计的阵列接收信号的协方差矩阵向量化来构造的。最近提出的嵌套式阵列NA和互质阵列CA就是利用稀疏阵的差分合成阵列来提高自由度的。而通过实际研究表明,最小冗余阵列MRA通过有效的阵元结构配置，可以获得最大的自由度DOF。在给定阵元数S及虚拟阵列为均匀线性阵列的情况下，最小冗余阵列MRA可以获得最大的阵列孔径。然而，最小冗余阵列的阵元位置及能够获得的自由度都没有明确的闭式解。虽然已有文献，如美国H.L.Van Trees教授的Optimum arrayprocessing:part IV of detection,estimation,and modulation通过穷搜的方法获得了阵元数小于等于17的最小冗余阵列，但却没有简单的方法来预测更大的最小冗余阵列MRA。
相对于最小冗余阵列MRA复杂的阵列结构获取过程，通过设计嵌套式阵列NA和互质阵列CA以获取更高的自由度DOF则很容易。
嵌套式阵列NA结构的阵元位置具有闭式解，并且通过利用接收数据的二阶统计信息，使用S个阵元就可以获得O(S2)的自由度。嵌套式阵列是由两个或多个具有不同阵元间隔的均匀线性子阵组成，它的缺点是第一个子阵的阵元间距通常较小，会引起阵元间的互耦问题。为此，P.P.Vaidyanathan等人后来提出的互质阵列CA可以通过增大阵元间距来减小阵元间的互耦问题。互质阵列是由两个阵元数目分别是M、N的均匀线性子阵构成，其中M、N为互质的整数，它能够获得多于阵元数目的自由度，但是由它构造的虚拟阵列并不是一个完全填充的均匀线性阵列，即在某些位置上没有虚拟阵元。通过利用虚拟阵列而不是原始阵列来进行DOA估计，这样由互质阵列获得的自由度就可以用来实现目标个数大于阵元数的DOA估计。但是在通过阵列接收信号的协方差矩阵构造虚拟阵列时，虚拟阵列中等效信源被实际信源的功率所替代，因此这些等效信源就表现为完全相干的信号。P.P.Vaidyanathan等提出采用一种基于空间平滑的MUSIC算法来解相干。然而由于基于空间平滑的估计算法一般只适用于均匀线阵，DOA估计时就只能利用互质阵列构造的虚拟阵列的一部分，这就会带来自由度和阵列孔径的降低。
然而，即便阵列结构拥有较高的自由度和阵列孔径，还要考虑到阵列密度问题，较低的阵列密度会引起发射功率受限，这在实际工程中是不希望出现的。
综上，现有的几种阵列虽然均能获得多于阵元数目的自由度，但都存在一定的局限，因此需要设计新的综合性能更优的阵列。
发明内容
本发明的目的在于针对上述已有阵列的不足，提出一种嵌套式子阵阵列的配置方法以减少阵元间的互耦，形成完全填充的均匀线性虚拟阵列，并获得较高的阵列密度、自由度和阵列孔径，提高波达方向DOA的估计精度。
本发明的技术思路是:根据总阵元数选择子阵数目及子阵内阵元数，利用现有文献结果得到相应的子阵结构和子阵内阵元结构；通过计算子阵间的最优单元间距，构造嵌套式子阵阵列；计算嵌套式子阵阵列的虚拟阵列以获得差分合成阵列，计算嵌套式子阵阵列的自由度；对差分合成阵列进行进行波达方向DOA估计。其实现方案包括如下：
1)构造嵌套式子阵阵列：
1a)给定总阵元数S，对S进行因式分解，得到子阵内阵元数目M及子阵数目N；
1b)根据子阵内阵元数目M，设计子阵内阵元结构，同时计算子阵内阵元的位置矢量uM：
uM＝[m1,m2,…,mi,…,mM]·d，
其中，mi表示子阵内第i个阵元位置系数，1≤i≤M，且m1＝0，d为入射信号的半波长；
1c)根据子阵数目N，设计子阵结构，同时计算N个子阵的阵位置矢量uN：
uN＝[n1,n2,…,nj,…,nN]·D，
其中，nj表示第j个子阵位置系数，1≤j≤N，且n1＝0，D＝L·d为子阵间的最优间隔单元，L表示子阵间最优间隔单元系数；
1d)根据上述设计参数，计算嵌套式子阵阵列位置矢量v：
其中pk表示嵌套式子阵阵列第k个阵元的位置系数，k＝1,2,…,S，符号表示交叉求和；
2)根据嵌套式子阵阵列，得到嵌套式子阵阵列的接收数据X(t)；
3)根据阵列接收数据X(t)估计嵌套式子阵阵列的协方差矩阵RXX，向量化该协方差矩阵得到向量化接收数据zct，去除zct中重复的元素，得到虚拟的差分合成阵列接收数据zc：
4)将差分合成阵列接收数据zc划分为lc个子阵接收数据，其中zcg是第g个子阵接收数据，g＝1,2,…,lc；计算第g个子阵的协方差矩阵Rg，并对该Rg进行前后向平滑平均，得到秩恢复的数据协方差矩阵RSS，其中lc表示划分的子阵个数，fv表示嵌套式子阵阵列的自由度；
5)根据秩恢复的数据协方差矩阵RSS估计出空间谱函数S(θ)，计算S(θ)中Q个极大值S(θq)，该S(θq)对应的角度θq即为波达方向，其中，θq表示第q个入射信号的入射角，q＝1,2,…,Q，Q表示入射信号的总个数。
本发明与现有阵列结构相比具有以下优点：
1)能获得孔径更大的阵列
最小冗余阵列MRA的阵元位置及可用自由度都没有闭式解，其设计往往需要复杂的计算机穷搜，而本发明阵列的阵元位置及可用自由度可利用已知的阵列的参数预测获得，实现简单，易于获得孔径更大的阵列。
2)测角性能好
现有互质阵列CA构造的虚拟差分合成阵列由于不是完全填充的均匀线阵，会降低DOA估计可用的自由度及测角精度；而本发明阵列构造的差分合成阵列为完全填充的均匀线阵，不仅在同样物理阵元数情况下可获得比CA更大的阵列孔径和更高的自由度，而且能估计更多的目标，具有更好的测角性能。
3)易于生产、更换和扩展
本发明的所有子阵结构相同，便于批量生产，对出现问题的子阵更换简单，可方便安装和扩展成更大的阵列。
4)成本低
本发明阵列是稀疏阵，在保持阵列面积不变的情况下，能够有效减少阵元数目，降低成本，由于各个子阵完全相同，批量生产更利于降低成本。
5)机动性高，角度分辨能力强
本发明的所有子阵都可当成一个子雷达系统，单个子雷达体积较小，机动性好，便于变换阵地；同时由于阵列的稀疏配置，可获得更好的角度分辨能力。
附图说明
图1是本发明的实现流程图；
图2是本发明第2实施例中的16阵元的阵列配置几何结构示意图；
图3是本发明第6实施例中的16阵元的阵列配置几何结构示意图；
图4是本发明第8实施例中的16阵元的阵列配置几何结构示意图；
图5是对用本发明九种实施例阵列配置的差分合成阵列自由度与总阵元数关系比较图；
图6是对用本发明九种实施例阵列配置的阵列孔径与总阵元数关系比较图；
图7是对用本发明九种实施例的阵列配置的阵列密度与总阵元数关系比较图；
图8是对用本发明五种实施例的阵列配置的波达方向DOA估计均方根误差与信噪比关系比较图；
图9是对用本发明五种实施例的阵列配置的波达方向DOA估计均方根误差与快拍数关系比较图。
具体实施方式
以下参照附图对本发明的实施和技术效果作进一步说明：
参照图1，本发明给出如下9种实施例：
实施例1，子阵结构和子阵内阵元结构均采用均匀线阵ULA结构的嵌套式子阵阵列波达方向估计。
步骤1，给定总阵元数S构造子阵结构及子阵内阵元结构。
现有技术中的非均匀阵列结构包括三：一种是最小冗余阵列MRA，其构造是根据总阵元数S得到最小冗余序列和最小冗余阵列MRA的阵列结构；第二种是嵌套式阵列NA，其构造是根据总阵元数S得到阵列的嵌套级数，然后根据嵌套级数确定每一级上的阵元数目，从而得到嵌套式阵列NA的阵列结构；第三种是互质阵列CA，其结构是根据总阵元数S得到两个子阵的阵元数目，从而得到互质阵列CA的阵列结构。
本发明的子阵结构及子阵内阵元结构均采用均匀线阵ULA结构，其设计步骤如下：
1a)给定总阵元数S，对S进行因式分解，得到子阵内阵元数M及子阵数目N：
S＝M·N，
其中M和N的取值最接近，且M≥2,N≥2；
1b)根据子阵内阵元数目M，将子阵内阵元结构设计为与现有的均匀线阵ULA结构相同，同时计算子阵内阵元的位置矢量u1M：
u1M＝[m11,m12,…,m1i,…m1M,]·d，
其中，m1i表示子阵内阵元结构为均匀线阵ULA结构的第i个阵元位置系数，1≤i≤M，且m11＝0，d为入射信号的半波长；
1c)根据子阵数目N，将子阵结构设计为与现有的均匀线阵ULA结构相同，同时计算N个子阵的阵位置矢量u1N：
u1N＝[n11,n12,…,n1j,…n1M,]·D，
其中，n1j表示子阵结构为均匀线阵ULA结构的第j个子阵位置系数，1≤j≤N，且n11＝0，D＝L·d为子阵间的最优间隔单元，L表示子阵间最优间隔单元系数；
1d)根据上述设计参数，计算嵌套式子阵阵列位置矢量v1、嵌套式子阵阵列的阵列孔径a1v和嵌套式子阵阵列的自由度f1v：
a1v＝(N-1)(2M-1)+M-1，
f1v＝(2N-1)(2M-1)，
其中符号表示交叉求和；p1k表示嵌套式子阵阵列第k个阵元的位置系数，k＝1,2,…,S；
1f)根据嵌套式子阵阵列位置矢量v1，得到所有阵元距离参考阵元的间距，从而确定虚拟阵列为完全填充型均匀线阵的嵌套式子阵阵列。
步骤2，根据所设计的嵌套式子阵阵列位置矢量，得到嵌套式子阵阵列的接收数据X(t)。
现有技术中对阵列接收数据X(t)的表达形式有窄带信号接收数据，宽带信号接收数据，本发明采用窄带信号接收数据表达形式，其具体实现如下：
2a)根据嵌套式子阵阵列位置矢量v，计算第q个入射信号的导向矢量a(θq)：
其中，表示第q个入射信号在第k个阵元处的旋转因子，pk表示嵌套式子阵阵列第k个阵元的位置系数，q＝1,2,…,Q，Q表示入射信号的总个数；
2b)根据入射信号的导向矢量a(θq)，构造导向矢量矩阵A：
A＝[a(θ1),a(θ2),…,a(θq),…,a(θQ)]，
2c)根据导向矢量矩阵A，入射信号S(t)，复高斯白噪声n(t)，得到嵌套式子阵阵列的接收数据X(t)：
X(t)＝AS(t)+n(t)，
其中，S(t)＝[s1(t),s2(t),…,sq(t),…,sQ(t)]T，sq(t)表示第q个入射信号，()T表示向量的转置，n(t)表示均值为0，方差为的复高斯白噪声，并且与入射信号不相关，t表示采样时间t＝1,2,…,T，T表示快拍数。
步骤3，根据嵌套式子阵阵列的接收数据X(t)估计阵列协方差矩阵RXX，得到虚拟的差分合成阵列接收数据zc。
3a)根据嵌套式子阵阵列的接收数据X(t)估计阵列协方差矩RXX：
其中，其中E表示求数学期望；(·)H表示共轭转置；Is表示S×S的单位矩阵，表示第q个入射信号的功率；
3b)根据阵列协方差矩RXX，计算向量化接收数据zct
其中，vec表示矩阵的向量化操作；B表示虚拟的差分合成阵列的方位矢量矩阵， 表示Kronecker积，()*表示取共轭，ei表示长度为S,第i个元素为1，其他元素全为0的行向量，()T表示转置，表示入射信号的功率矢量；
3c)根据向量化接收数据zct，得到虚拟的差分合成阵列接收数据zc：
zc＝Dis(zct)，
其中，Dis表示去除向量中重复元素的操作。
步骤4，根据虚拟的差分合成阵列接收数据zc，计算秩恢复的数据协方差矩阵RSS。
4a)根据虚拟的差分合成阵列接收数据zc，得到子阵接收数据zcg，即取差分合成阵列接收数据zc中位于(g+n-1)处的lc个数据构成子阵接收数据zcg，n＝1,2,…,lc,g＝1,2,…,lc，
其中，lc表示划分的子阵个数，fv表示嵌套式子阵阵列的自由度；
4b)根据子阵接收数据zcg，计算第g个子阵的协方差矩阵Rg：
4c)根据第g子阵的协方差矩阵Rg，计算前向平滑的协方差矩阵Rf：
4d)根据子阵接收数据zcg计算第g个子阵的后向接收数据ycg；
ycg＝Trans(zcg)，
其中，Trans表示左右翻转并取共轭；
4e)根据第g个子阵的后向接收数据ycg；计算第g个子阵的后向协方差矩阵Rg b：
4f)根据第g个子阵的后向协方差矩阵Rg b，计算后向平滑的协方差矩阵Rb：
4g)根据前向平滑的协方差矩阵Rf和后向平滑的协方差矩阵Rb，计算秩恢复的数据协方差矩阵RSS：
以上步骤3-步骤4内容皆为现有常用方法。
步骤5，根据RSS估计出空间谱函数S(θ)，计算S(θ)中Q个极大值S(θq)，S(θq)对应的角度θq即为波达方向。
现有技术中利用RSS估计出空间谱函数S(θ)有最大似然算法ML，加权子空间拟合法WSF，波束成形法DBF，基于特征空间的多重信号分类法MUSIC，旋转因子不变法ESPRIT，本发明中采用基于特征空间的多重信号分类法MUSIC，其具体实现如下：
5a)对RSS进行特征值分解：
其中，λt表示矩阵RSS的特征值，Vt表示矩阵RSS对应于λt的特征矢量，t表示采样时间，t＝1,2,…，T，T表示快拍数，Q表示入射信号的总个数；
5b)根据Vt，计算空间谱函数S(θ)：
其中，a(θ)表示入射角度为θ的入射信号的导向矢量，
表示第q个入射信号在第k个阵元处的旋转因子，pk表示嵌套式子阵阵列第k个阵元的位置系数；
5c)根据空间谱函数S(θ)，计算S(θ)中Q个极大值S(θq)，S(θq)对应的角度θq即为波达方向。
实施例2，子阵结构采用均匀线阵ULA结构、子阵内阵元结构采用最小冗余阵列MRA结构的嵌套式子阵阵列波达方向估计。
该实施例2与实施例1的区别在于步骤1构造子阵结构及子阵内阵元结构不同，其他步骤均与实施例1相同，其构造子阵结构及子阵内阵元结构步骤如下：
(1a)给定总阵元数S，对S进行因式分解，得到子阵内阵元数M及子阵数目N：
S＝M·N，
其中M和N的取值最接近，且M≥2,N≥2；
(1b)根据子阵内阵元数目M，将子阵内阵元结构设计为与现有的最小冗余阵列MRA结构相同，同时计算子阵内阵元的位置矢量u2M：
u2M＝[m21,m22,…,m2i,…m2M,]·d，
其中，m2i表示子阵内阵元结构为最小冗余阵列MRA结构的第i个阵元位置系数，1≤i≤M，且m21＝0，d为入射信号的半波长；
(1c)根据子阵数目N，将子阵结构设计为与现有的均匀线阵ULA结构相同，同时计算N个子阵的阵位置矢量u2N：
u2N＝[n21,n22,…,n2j,…n2M,]·D，
其中，n2j表示子阵结构为均匀线阵ULA结构的第j个子阵位置系数，1≤j≤N，且n21＝0，D＝L·d为子阵间的最优间隔单元，L表示子阵间最优间隔单元系数；
(1d)根据上述设计参数，计算嵌套式子阵阵列位置矢量v2、嵌套式子阵阵列的阵列孔径a2v和嵌套式子阵阵列的自由度f2v：
a2v＝(N-1)[M(M-1)+1-2fM]+M(M-1)/2-fM，
f2v＝(2N-1)[M(M-1)+1-2fM]，
其中符号表示交叉求和；p2k表示嵌套式子阵阵列第k个阵元的位置系数，k＝1,2,…,S，fM表示子阵内阵元结构为最小冗余阵列MRA结构的差分合成阵列的冗余度；
(1e)根据嵌套式子阵阵列位置矢量v2，得到所有阵元距离参考阵元的间距，从而确定虚拟阵列为完全填充型均匀线阵的嵌套式子阵阵列，其结构如图2所示。
实施例3，子阵结构采用均匀线阵ULA结构、子阵内阵元结构采用嵌套式线阵NA结构的嵌套式子阵阵列波达方向估计。
该实施例3与实施例1的区别在于步骤1构造子阵结构及子阵内阵元结构不同，其他步骤均与实施例1相同，其构造子阵结构及子阵内阵元结构步骤如下：
1-a)给定总阵元数S，对S进行因式分解，得到子阵内阵元数M及子阵数目N：
S＝M·N，
其中M和N的取值最接近，且M≥2,N≥2；
1-b)根据子阵内阵元数目M，将子阵内阵元结构设计为与现有的嵌套式线阵NA结构相同，同时计算子阵内阵元的位置矢量u3M：
u3M＝[m31,m32,…,m3i,…m3M,]·d，
其中，m3i表示子阵内阵元结构为嵌套式线阵NA结构的第i个阵元位置系数，1≤i≤M，且m31＝0，d为入射信号的半波长；
1-c)根据子阵数目N，将子阵结构设计为与现有的均匀线阵ULA结构相同，同时计算N个子阵的阵位置矢量u3N：
u3N＝[n31,n32,…,n3j,…n3M,]·D，
其中，n3j表示子阵结构为均匀线阵ULA结构的第j个子阵位置系数，1≤j≤N，且n31＝0，D＝L·d为子阵间的最优间隔单元，L表示子阵间最优间隔单元系数；
1-d)根据上述设计参数，计算嵌套式子阵阵列位置矢量v3、嵌套式子阵阵列的阵列孔径a3v和嵌套式子阵阵列的自由度f3v：
其中符号表示交叉求和；p3k表示嵌套式子阵阵列第k个阵元的位置系数，k＝1,2,…,S；
1-e)根据嵌套式子阵阵列位置矢量v3，得到所有阵元距离参考阵元的间距，从而确定虚拟阵列为完全填充型均匀线阵的嵌套式子阵阵列。
实施例4，子阵结构和子阵内阵元结构均采用最小冗余阵列MRA结构的嵌套式子阵阵列波达方向估计。
该实施例4与实施例1的区别在于步骤1构造子阵结构及子阵内阵元结构不同，其他步骤均与实施例1相同，其构造子阵结构及子阵内阵元结构步骤如下：
(1-a)给定总阵元数S，对S进行因式分解，得到子阵内阵元数M及子阵数目N：
S＝M·N，
其中M和N的取值最接近，且M≥2,N≥2；
(1-b)根据子阵内阵元数目M，将子阵内阵元结构设计为与现有的最小冗余阵列MRA结构相同，同时计算子阵内阵元的位置矢量u4M：
u4M＝[m41,m42,…,m4i,…m4M,]·d，
其中，m4i表示子阵内阵元结构为最小冗余阵列MRA结构的第i个阵元位置系数，1≤i≤M，且m41＝0，d为入射信号的半波长；
(1-c)根据子阵数目N，将子阵结构设计为与现有的最小冗余阵列MRA结构相同，同时计算N个子阵的阵位置矢量u4N：
u4N＝[n41,n42,…,n4j,…n4M,]·D，
其中，n4j表示子阵结构为均匀线阵ULA结构的第j个子阵位置系数，1≤j≤N，且n41＝0，D＝L·d为子阵间的最优间隔单元，L表示子阵间最优间隔单元系数；
(1-d)根据上述设计参数，计算嵌套式子阵阵列位置矢量v4、嵌套式子阵阵列的阵列孔径a4v和嵌套式子阵阵列的自由度f4v：
a4v＝[N(N-1)/2-fN][M(M-1)+1-2fM]+M(M-1)/2-fM，
f4v＝[N(N-1)+1-2fN][M(M-1)/2-fM]，
其中符号表示交叉求和；p4k表示嵌套式子阵阵列第k个阵元的位置系数，k＝1,2,…,S，fN表示子阵结构为最小冗余阵列MRA结构的差分合成阵列的冗余度；
(1-e)根据嵌套式子阵阵列位置矢量v4，得到所有阵元距离参考阵元的间距，从而确定虚拟阵列为完全填充型均匀线阵的嵌套式子阵阵列。
实施例5，子阵结构采用最小冗余阵列MRA结构、子阵内阵元结构采用均匀线阵ULA结构的嵌套式子阵阵列波达方向估计。
该实施例5与实施例1的区别在于步骤1构造子阵结构及子阵内阵元结构不同，其他步骤均与实施例1相同，其构造子阵结构及子阵内阵元结构步骤如下：
1A)给定总阵元数S，对S进行因式分解，得到子阵内阵元数M及子阵数目N：
S＝M·N，
其中M和N的取值最接近，且M≥2,N≥2；
1B)根据子阵内阵元数目M，将子阵内阵元结构设计为与现有的均匀线阵ULA结构相同，同时计算子阵内阵元的位置矢量u5M：
u5M＝[m51,m52,…,m5i,…m5M,]·d，
其中，m5i表示子阵内阵元结构为均匀线阵ULA结构的第i个阵元位置系数，1≤i≤M，且m51＝0，d为入射信号的半波长；
1C)根据子阵数目N，将子阵结构设计为与现有的最小冗余阵列MRA结构相同，同时计算N个子阵的阵位置矢量u5N：
u5N＝[n51,n52,…,n5j,…n5M,]·D，
其中，n5j表示子阵结构为最小冗余阵列MRA结构的第j个子阵位置系数，1≤j≤N，且n51＝0，D＝L·d为子阵间的最优间隔单元，L表示子阵间最优间隔单元系数；
1D)根据上述设计参数，计算嵌套式子阵阵列位置矢量v5、嵌套式子阵阵列的阵列孔径a5v和嵌套式子阵阵列的自由度f5v：
a5v＝[N(N-1)/2-fN](2M-1)+M-1，
f5v＝[N(N-1)+1-2fN](2M-1)，
其中符号表示交叉求和；p5k表示嵌套式子阵阵列第k个阵元的位置系数，k＝1,2,…,S；
1E)根据嵌套式子阵阵列位置矢量v5，得到所有阵元距离参考阵元的间距，从而确定虚拟阵列为完全填充型均匀线阵的嵌套式子阵阵列。
实施例6，子阵结构采用最小冗余阵列MRA结构、子阵内阵元结构采用嵌套式线阵NA结构的嵌套式子阵阵列波达方向估计。
该实施例6与实施例1的区别在于步骤1构造子阵结构及子阵内阵元结构不同，其他步骤均与实施例1相同，其构造子阵结构及子阵内阵元结构步骤如下：
(1A)给定总阵元数S，对S进行因式分解，得到子阵内阵元数M及子阵数目N：
S＝M·N，
其中M和N的取值最接近，且M≥2,N≥2；
(1B)根据子阵内阵元数目M，将子阵内阵元结构设计为与现有的嵌套式线阵NA结构相同，同时计算子阵内阵元的位置矢量u6M：
u6M＝[m61,m62,…,m6i,…m6M,]·d，
其中，m6i表示子阵内阵元结构为嵌套式线阵NA结构的第i个阵元位置系数，1≤i≤M，且m61＝0，d为入射信号的半波长；
(1C)根据子阵数目N，将子阵结构设计为与现有的最小冗余阵列MRA结构相同，同时计算N个子阵的阵位置矢量u6N：
u6N＝[n61,n62,…,n6j,…n6M,]·D，
其中，n6j表示子阵结构为最小冗余阵列MRA结构的第j个子阵位置系数，1≤j≤N，且n61＝0，D＝L·d为子阵间的最优间隔单元，L表示子阵间最优间隔单元系数；
(1D)根据上述设计参数，计算嵌套式子阵阵列位置矢量v6、嵌套式子阵阵列的阵列孔径a6v和嵌套式子阵阵列的自由度f6v：
其中符号表示交叉求和；p6k表示嵌套式子阵阵列第k个阵元的位置系数，k＝1,2,…,S；
(1E)根据嵌套式子阵阵列位置矢量v6，得到所有阵元距离参考阵元的间距，从而确定虚拟阵列为完全填充型均匀线阵的嵌套式子阵阵列，其结构如图3所示。
实施例7，子阵结构和子阵内阵元结构均采用嵌套式线阵NA结构的嵌套式子阵阵列波达方向估计。
该实施例7与实施例1的区别在于步骤1构造子阵结构及子阵内阵元结构不同，其他步骤均与实施例1相同，其构造子阵结构及子阵内阵元结构步骤如下：
1-1)给定总阵元数S，对S进行因式分解，得到子阵内阵元数M及子阵数目N：
S＝M·N，
其中M和N的取值最接近，且M≥2,N≥2；
1-2)根据子阵内阵元数目M，将子阵内阵元结构设计为与现有的嵌套式线阵NA结构相同，同时计算子阵内阵元的位置矢量u7M：
u7M＝[m71,m72,…,m7i,…m7M,]·d，
其中，m7i表示子阵内阵元结构为嵌套式线阵NA结构的第i个阵元位置系数，1≤i≤M，且m71＝0，d为入射信号的半波长；
1-3)根据子阵数目N，将子阵结构设计为与现有的嵌套式线阵NA结构相同，同时计算N个子阵的阵位置矢量u7N：
u7N＝[n71,n72,…,n7j,…n7M,]·D，
其中，n7j表示子阵结构为嵌套式线阵NA结构的第j个子阵位置系数，1≤j≤N，且n71＝0，D＝L·d为子阵间的最优间隔单元，L表示子阵间最优间隔单元系数；
1-4)根据上述设计参数，计算嵌套式子阵阵列位置矢量v7、嵌套式子阵阵列的阵列孔径a7v和嵌套式子阵阵列的自由度f7v：
其中符号表示交叉求和；p7k表示嵌套式子阵阵列第k个阵元的位置系数，k＝1,2,…,S；
1-5)根据嵌套式子阵阵列位置矢量v7，得到所有阵元距离参考阵元的间距，从而确定虚拟阵列为完全填充型均匀线阵的嵌套式子阵阵列。
实施例8，子阵结构采用嵌套式线阵NA结构、子阵内阵元结构采用均匀线阵ULA结构的嵌套式子阵阵列波达方向估计。
该实施例8与实施例1的区别在于步骤1构造子阵结构及子阵内阵元结构不同，其他步骤均与实施例1相同，其构造子阵结构及子阵内阵元结构步骤如下：
(1-1)给定总阵元数S，对S进行因式分解，得到子阵内阵元数M及子阵数目N：
S＝M·N，
其中M和N的取值最接近，且M≥2,N≥2；
(1-2)根据子阵内阵元数目M，将子阵内阵元结构设计为与现有的均匀线阵ULA结构相同，同时计算子阵内阵元的位置矢量u8M：
u8M＝[m81,m82,…,m8i,…m8M,]·d，
其中，m8i表示子阵内阵元结构为均匀线阵ULA结构的第i个阵元位置系数，1≤i≤M，且m81＝0，d为入射信号的半波长；
(1-3)根据子阵数目N，将子阵结构设计为与现有的嵌套式线阵NA结构相同，同时计算N个子阵的阵位置矢量u8N：
u8N＝[n81,n82,…,n8j,…n8M,]·D，
其中，n8j表示子阵结构为嵌套式线阵NA结构的第j个子阵位置系数，1≤j≤N，且n81＝0，D＝L·d为子阵间的最优间隔单元，L表示子阵间最优间隔单元系数；
(1-4)根据上述设计参数，计算嵌套式子阵阵列位置矢量v8、嵌套式子阵阵列的阵列孔径a8v和嵌套式子阵阵列的自由度f8v：
其中符号表示交叉求和；p8k表示嵌套式子阵阵列第k个阵元的位置系数，k＝1,2,…,S；
(1-5)根据嵌套式子阵阵列位置矢量v8，得到所有阵元距离参考阵元的间距，从而确定虚拟阵列为完全填充型均匀线阵的嵌套式子阵阵列，其结构如图4所示。
实施例9，子阵结均采用嵌套式线阵NA结构、子阵内阵元结构采用最小冗余阵列MRA结构的嵌套式子阵阵列波达方向估计。
该实施例9与实施例1的区别在于步骤1构造子阵结构及子阵内阵元结构不同，其他步骤均与实施例1相同，其构造子阵结构及子阵内阵元结构步骤如下：
(1-A)给定总阵元数S，对S进行因式分解，得到子阵内阵元数M及子阵数目N：
S＝M·N，
其中M和N的取值最接近，且M≥2,N≥2；
(1-B)根据子阵内阵元数目M，将子阵内阵元结构设计为与现有的最小冗余阵列MRA结构相同，同时计算子阵内阵元的位置矢量u9M：
u9M＝[m91,m92,…,m9i,…m9M,]·d，
其中，m9i表示子阵内阵元结构为最小冗余阵列MRA结构的第i个阵元位置系数，1≤i≤M，且m91＝0，d为入射信号的半波长；
(1-C)根据子阵数目N，将子阵结构设计为与现有的嵌套式线阵NA结构相同，同时计算N个子阵的阵位置矢量u9N：
u9N＝[n91,n92,…,n9j,…n9M,]·D，
其中，n9j表示子阵结构为嵌套式线阵NA结构的第j个子阵位置系数，1≤j≤N，且n91＝0，D＝L·d为子阵间的最优间隔单元，L表示子阵间最优间隔单元系数；
(1-D)根据上述设计参数，计算嵌套式子阵阵列位置矢量v9、嵌套式子阵阵列的阵列孔径a9v和嵌套式子阵阵列的自由度f9v：
其中符号表示交叉求和；p9k表示嵌套式子阵阵列第k个阵元的位置系数，k＝1,2,…,S；
(1-E)根据嵌套式子阵阵列位置矢量v9，得到所有阵元距离参考阵元的间距，从而确定虚拟阵列为完全填充型均匀线阵的嵌套式子阵阵列。
本发明的效果通过以下计算仿真进一步说明：
仿真1：对本发明的九种实施例的阵列配置的差分合成阵列的自由度、有效孔径和阵列密度与总阵元数的关系进行仿真。
1.1)仿真条件：总阵元数为9～36间的一些整数。
1.2)仿真内容与结果
1.2a)在上述1.1)仿真条件下，对本发明九种实施例阵列配置的差分合成阵列自由度与总阵元数的关系进行仿真，结果如图5所示；
1.2b)在上述1.1)仿真条件下，对本发明九种实施例阵列配置的有效孔径与总阵元数关系进行仿真，结果如图6所示；
1.2c)在上述1.1)仿真条件下，对本发明九种实施例阵列配置的阵列密度与总阵元数关系进行仿真，结果如图7所示。
从图5可以看出，差分合成阵列的自由度随总阵元数的增加而增大；当阵元数目相同时，实施例4具有最高的自由度，实施例1的自由度最低。
从图6可以看出，实施例4具有最高的阵列孔径。随着总阵元数的增大，各实施例的阵列孔径差距也愈加明显。
从图7可以看出，阵列密度随总阵元数的增加而减小；实施例1的阵列密度最高。
仿真2：对本发明五种实施例的阵列配置方法的波达方向DOA估计角度均方根误差与信噪比关系进行仿真。
2.1)仿真条件：阵元数S＝12，信源数Q＝16，探测角度θ，-70°≤θ≤70°快拍数T＝500,500次Monte-Carlo实验。
2.2)仿真内容与结果：
在上述2.1)仿真条件下对本发明五种实施例的阵列配置的波达方向DOA估计角度均方根误差与信噪比关系进行仿真，结果如图8所示。
从图8可以看出，DOA估计性能随着信噪比的提高而提高，实施例4具有最小的DOA估计均方根误差。这里只选择五种实施例作比较的原因是总阵元数S＝12时，子阵内阵元数为3的最小冗余阵列MRA和嵌套式阵列NA阵列结构相同，因此实施例2和实施例3结构相同，实施例4和实施例6结构相同，实施例7和实施例9结构相同。
仿真3：对本发明五种实施例的阵列配置的波达方向DOA估计角度均方根误差与快拍数关系进行仿真。
3.1)仿真条件：设阵元数S＝12，信源数Q＝16，探测角度θ，-70°≤θ≤70°，信噪比SNR＝0dB，500次Monte-Carlo实验。
3.2)仿真内容与结果：
在上述3.1)仿真条件下，对本发明五种实施例的阵列配置的波达方向DOA估计角度均方根误差与快拍数关系进行仿真，结果如图9所示。
从图9可以看出，DOA估计性能随着快拍数的增加而提高，实施例4具有最小的DOA估计均方根误差。这里只选择五种实施例作比较的原因是总阵元数S＝12时，子阵内阵元数为3的最小冗余阵列MRA和嵌套式阵列NA阵列结构相同，因此实施例2和实施例3结构相同，实施例4和实施例6结构相同，实施例7和实施例9结构相同。
综上，本发明阵列配置灵活，具有较高的自由度、阵列孔径和阵列密度，较好的波达方向DOA估计测角性能。

Claims (8)

1.一种基于嵌套式子阵阵列的波达方向估计方法，包括：
1)构造嵌套式子阵阵列：
1a)给定总阵元数S，对S进行因式分解，得到子阵内阵元数目M及子阵数目N；
1b)根据子阵内阵元数目M，设计子阵内阵元结构，同时计算子阵内阵元的位置矢量uM：
uM＝[m1,m2,…,mi,…,mM]·d，
其中，mi表示子阵内第i个阵元位置系数，1≤i≤M，且m1＝0，d为入射信号的半波长；
1c)根据子阵数目N，设计子阵结构，同时计算N个子阵的阵位置矢量uN：
uN＝[n1,n2,…,nj,…,nN]·D，
其中，nj表示第j个子阵位置系数，1≤j≤N，且n1＝0，D＝L·d为子阵间的最优间隔单元，L表示子阵间最优间隔单元系数，L的选取由子阵内阵元结构决定：
当子阵内阵元结构设计为均匀线阵ULA结构时，子阵间最优间隔单元系数为L＝2M-1；
当子阵内阵元结构设计为最小冗余阵列MRA时，子阵间最优间隔单元系数为：
L＝M(M-1)+1-2fM；
当子阵内阵元结构设计为嵌套式线阵NLA时，根据子阵内阵元数目M的奇偶性确定：
若子阵内阵元数目M为偶数，则子阵间最优间隔单元系数L＝M2/2+M-1；
若子阵内阵元数目M为奇数，则子阵间最优间隔单元系数L＝(M+1)2/2-1；
1d)根据上述设计参数，计算嵌套式子阵阵列位置矢量v：
<mrow> <mi>v</mi> <mo>=</mo> <msub> <mi>u</mi> <mi>N</mi> </msub> <mo>&amp;CirclePlus;</mo> <msub> <mi>u</mi> <mi>M</mi> </msub> <mo>=</mo> <mo>{</mo> <msub> <mi>n</mi> <mi>j</mi> </msub> <mo>&amp;CenterDot;</mo> <mi>D</mi> <mo>+</mo> <msub> <mi>m</mi> <mi>i</mi> </msub> <mo>&amp;CenterDot;</mo> <mi>d</mi> <mo>}</mo> <mo>=</mo> <mo>&amp;lsqb;</mo> <msub> <mi>p</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>p</mi> <mn>2</mn> </msub> <mo>,</mo> <mo>...</mo> <mo>,</mo> <msub> <mi>p</mi> <mi>k</mi> </msub> <mo>,</mo> <mo>...</mo> <mo>,</mo> <msub> <mi>p</mi> <mi>S</mi> </msub> <mo>&amp;rsqb;</mo> <mo>&amp;CenterDot;</mo> <mi>d</mi> <mo>,</mo> </mrow>
其中pk表示嵌套式子阵阵列第k个阵元的位置系数，k＝1,2,…,S，符号表示交叉求和；
2)根据嵌套式子阵阵列，得到嵌套式子阵阵列的接收数据X(t)；
3)根据阵列接收数据X(t)估计嵌套式子阵阵列的协方差矩阵RXX，向量化该协方差矩阵得到向量化接收数据zct，去除zct中重复的元素，得到虚拟的差分合成阵列接收数据zc：
4)将差分合成阵列接收数据zc划分为lc个子阵接收数据，其中zcg是第g个子阵接收数据，g＝1,2,…,lc；计算第g个子阵的协方差矩阵Rg，并对该Rg进行前后向平滑平均，得到秩恢复的数据协方差矩阵RSS，其中lc表示划分的子阵个数，fv表示嵌套式子阵阵列的自由度；
5)根据秩恢复的数据协方差矩阵RSS估计出空间谱函数S(θ)，计算S(θ)中Q个极大值S(θq)，该S(θq)对应的角度θq即为波达方向，其中，θq表示第q个入射信号的入射角，q＝1,2,…,Q，Q表示入射信号的总个数。
2.根据权利要求1所述的基于嵌套式子阵阵列的波达方向估计方法，其中步骤1a)中的子阵内阵元数M及子阵数目N，通过按如下公式计算：
S＝M·N，
其中M和N的取值最接近，且M≥2,N≥2。
3.根据权利要求1所述的基于嵌套式子阵阵列的波达方向估计方法，其中步骤1b)中子阵内阵元结构的设计，采用现有均匀线阵ULA结构或最小冗余阵列MRA结构或嵌套式线阵NLA的结构。
4.根据权利要求1所述的基于嵌套式子阵阵列的波达方向估计方法，其中步骤1c)中子阵结构的设计，采用现有均匀线阵ULA结构或最小冗余阵列MRA结构嵌套式线阵NLA的结构。
5.根据权利要求1所述的基于嵌套式子阵阵列的波达方向估计方法，其中步骤2)中的嵌套式子阵阵列的接收数据X(t)，其表示如下：
X(t)＝AS(t)+n(t)，
其中A为导向矢量矩阵，A＝[a(θ1),a(θ2),…,a(θq),…,a(θQ)]，a(θq)表示第q个入射信号的导向矢量， 表示第q个入射信号在第k个阵元处的旋转因子，pk表示嵌套式子阵阵列第k个阵元的位置系数，q＝1,2,…,Q，Q表示入射信号的总个数，S(t)表示入射信号，S(t)＝[s1(t),s2(t),…，sq(t),…，sQ(t)]T，sq(t)表示第q个入射信号，(·)T:表示向量的转置，n(t)表示均值为0，方差为的复高斯白噪声，且与入射信号不相关。
6.根据权利要求1所述的基于嵌套式子阵阵列的波达方向估计方法，其中步骤3)中虚拟的差分合成阵列接收数据zc，按如下步骤计算：
3.1)根据嵌套式子阵阵列的接收数据X(t)估计阵列协方差矩RXX：
RXX＝E[X(t)XH(t)]，
其中，其中E表示求数学期望；(·)H表示共轭转置；
3.2)根据阵列协方差矩RXX，计算向量化接收数据zct：
zct＝vec(RXX)，
其中，vec表示矩阵的向量化操作；
3.3)根据向量化接收数据zct，得到虚拟的差分合成阵列接收数据zc：
zc＝Dis(zct)，
其中，Dis表示去除向量中重复元素的操作。
7.根据权利要求1所述的基于嵌套式子阵阵列的波达方向估计方法，其中步骤4)中秩恢复的数据协方差矩阵RSS，按如下步骤计算：
4.1)根据虚拟的差分合成阵列接收数据zc，得到的子阵接收数据zcg，其表示如下：
<mrow> <msub> <mi>z</mi> <mi>c</mi> </msub> <mo>=</mo> <mo>&amp;lsqb;</mo> <msub> <mi>z</mi> <mrow> <mi>c</mi> <mn>1</mn> </mrow> </msub> <mo>,</mo> <msub> <mi>z</mi> <mrow> <mi>c</mi> <mn>2</mn> </mrow> </msub> <mo>,</mo> <mo>...</mo> <mo>,</mo> <msub> <mi>z</mi> <mrow> <mi>c</mi> <mi>g</mi> </mrow> </msub> <mo>,</mo> <mo>...</mo> <mo>,</mo> <msub> <mi>z</mi> <mrow> <msub> <mi>cl</mi> <mi>c</mi> </msub> </mrow> </msub> <mo>&amp;rsqb;</mo> <mo>,</mo> </mrow>
其中zcg是第g个子阵接收数据，g＝1,2,…,lc；
4.2)根据子阵接收数据zcg，计算第g个子阵的协方差矩阵Rg：
<mrow> <msub> <mi>R</mi> <mi>g</mi> </msub> <mo>=</mo> <msub> <mi>z</mi> <mrow> <mi>c</mi> <mi>g</mi> </mrow> </msub> <msubsup> <mi>z</mi> <mrow> <mi>c</mi> <mi>g</mi> </mrow> <mi>H</mi> </msubsup> <mo>;</mo> </mrow>
4.3)根据第g子阵的协方差矩阵Rg，计算前向平滑的协方差矩阵Rf：
<mrow> <msup> <mi>R</mi> <mi>f</mi> </msup> <mo>=</mo> <mfrac> <mn>1</mn> <msub> <mi>l</mi> <mi>c</mi> </msub> </mfrac> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>g</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>l</mi> <mi>c</mi> </msub> </munderover> <msub> <mi>R</mi> <mi>g</mi> </msub> <mo>;</mo> </mrow>
4.4)根据子阵接收数据zcg计算第g个子阵的后向接收数据ycg；
ycg＝Trans(zcg)，
其中，Trans表示左右翻转并取共轭；
4.5)根据第g个子阵的后向接收数据ycg，计算第g个子阵的后向协方差矩阵Rg b：
<mrow> <msup> <msub> <mi>R</mi> <mi>g</mi> </msub> <mi>b</mi> </msup> <mo>=</mo> <msub> <mi>y</mi> <mrow> <mi>c</mi> <mi>g</mi> </mrow> </msub> <msubsup> <mi>y</mi> <mrow> <mi>c</mi> <mi>g</mi> </mrow> <mi>H</mi> </msubsup> <mo>;</mo> </mrow>
4.6)根据第g个子阵的后向协方差矩阵Rg b，计算后向平滑的协方差矩阵Rb：
<mrow> <msup> <mi>R</mi> <mi>b</mi> </msup> <mo>=</mo> <mfrac> <mn>1</mn> <msub> <mi>l</mi> <mi>c</mi> </msub> </mfrac> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>g</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>l</mi> <mi>c</mi> </msub> </munderover> <msubsup> <mi>R</mi> <mi>g</mi> <mi>b</mi> </msubsup> <mo>;</mo> </mrow>
4.7)根据前向平滑的协方差矩阵Rf和后向平滑的协方差矩阵Rb，计算秩恢复的数据协方差矩阵RSS：
<mrow> <msub> <mi>R</mi> <mrow> <mi>S</mi> <mi>S</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <msup> <mi>R</mi> <mi>f</mi> </msup> <mo>+</mo> <msup> <mi>R</mi> <mi>b</mi> </msup> </mrow> <mn>2</mn> </mfrac> <mo>.</mo> </mrow>
8.根据权利要求1所述的基于嵌套式子阵阵列的波达方向估计方法，其中步骤5)中空间谱函数S(θ)，按如下步骤计算：
5.1)对秩恢复的数据协方差矩阵RSS进行如下特征值分解：
<mrow> <msub> <mi>R</mi> <mrow> <mi>S</mi> <mi>S</mi> </mrow> </msub> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>t</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>T</mi> </munderover> <msub> <mi>&amp;lambda;</mi> <mi>t</mi> </msub> <msub> <mi>V</mi> <mi>t</mi> </msub> <msubsup> <mi>V</mi> <mi>t</mi> <mi>H</mi> </msubsup> <mo>,</mo> </mrow>
其中，λt表示矩阵RSS的特征值，Vt表示矩阵RSS对应于λt的特征矢量，t表示采样时间，t＝1,2,…,T，T表示快拍数，Q表示入射信号的总个数；
5.2)根据特征矢量Vt，计算空间谱函数S(θ)：
<mrow> <mi>S</mi> <mrow> <mo>(</mo> <mi>&amp;theta;</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mo>|</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>t</mi> <mo>=</mo> <mi>Q</mi> <mo>+</mo> <mn>1</mn> </mrow> <mi>T</mi> </munderover> <mi>a</mi> <mrow> <mo>(</mo> <mi>&amp;theta;</mi> <mo>)</mo> </mrow> <msub> <mi>V</mi> <mi>t</mi> </msub> <msubsup> <mi>V</mi> <mi>t</mi> <mi>H</mi> </msubsup> <mi>a</mi> <mrow> <mo>(</mo> <mi>&amp;theta;</mi> <mo>)</mo> </mrow> <mo>|</mo> </mrow> </mfrac> <mo>,</mo> </mrow>
其中，a(θ)表示入射角度为θ的入射信号的导向矢量，pk表示嵌套式子阵阵列第k个阵元的位置系数；
5.3)根据空间谱函数S(θ)，计算S(θ)中Q个极大值S(θq)，S(θq)对应的角度θq即为波达方向。
CN201610236372.3A 2016-04-15 2016-04-15 基于嵌套式子阵阵列的波达方向估计方法 Active CN105824002B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610236372.3A CN105824002B (zh) 2016-04-15 2016-04-15 基于嵌套式子阵阵列的波达方向估计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610236372.3A CN105824002B (zh) 2016-04-15 2016-04-15 基于嵌套式子阵阵列的波达方向估计方法

Publications (2)

Publication Number Publication Date
CN105824002A CN105824002A (zh) 2016-08-03
CN105824002B true CN105824002B (zh) 2018-05-04

Family

ID=56526956

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610236372.3A Active CN105824002B (zh) 2016-04-15 2016-04-15 基于嵌套式子阵阵列的波达方向估计方法

Country Status (1)

Country Link
CN (1) CN105824002B (zh)

Families Citing this family (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106291473B (zh) * 2016-08-19 2018-07-17 西安电子科技大学 嵌套式天线阵列设置方法
CN106707257B (zh) * 2016-12-01 2019-06-21 西安电子科技大学 基于嵌套阵列的mimo雷达波达方向估计方法
CN106646388B (zh) * 2016-12-16 2019-01-25 西安电子科技大学 基于嵌套阵列的mimo雷达抗干扰方法
CN106646344B (zh) * 2016-12-16 2019-02-01 西北工业大学 一种利用互质阵的波达方向估计方法
CN106501773B (zh) * 2016-12-23 2018-12-11 云知声（上海）智能科技有限公司 基于差分阵列的声源方向定位方法
CN106990385B (zh) * 2017-03-22 2019-11-05 西安电子科技大学 基于最小冗余线阵的非均匀面阵设计方法
CN107450047B (zh) * 2017-08-28 2020-06-16 电子科技大学 嵌套阵下基于未知互耦信息的压缩感知doa估计方法
CN108008348B (zh) * 2017-11-16 2023-07-18 华南理工大学 基于可调夹角均匀线阵的水下波达方向估计方法及装置
CN110082708B (zh) * 2019-02-25 2023-05-02 西安电子科技大学 非均匀阵列设计和波达方向估计方法
CN109932681B (zh) * 2019-04-04 2022-08-02 哈尔滨工程大学 一种基于空-时信息的降冗余嵌套阵列设置方法
CN110095749B (zh) * 2019-05-05 2022-08-02 哈尔滨工程大学 基于非圆信号的多重孔径嵌套阵列设置及doa估计方法
CN110275131B (zh) * 2019-05-30 2021-11-30 电子科技大学 一种基于虚拟差分阵列的doa跟踪方法及装置
CN110726967B (zh) * 2019-10-25 2021-08-03 北京理工大学 用于一维阵列测向的单边稀疏嵌套阵设计方法
CN111025236A (zh) * 2019-12-25 2020-04-17 西安电子科技大学昆山创新研究院 一种基于改进型嵌套线阵的非均匀面阵设计方法
CN111736118B (zh) * 2020-07-22 2020-11-17 东南大学 一种线列阵阵列扩展方法
CN112731282B (zh) * 2020-12-09 2023-12-26 浙江工业大学 一种基于虚拟子阵交织的声阵列空间滤波方法及系统
CN116540242B (zh) * 2023-07-03 2023-10-20 天津知海科技有限公司 干涉图像生成方法、装置、电子设备及可读存储介质

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102608565B (zh) * 2012-03-23 2014-01-08 哈尔滨工程大学 一种基于均匀圆阵列的波达方向估计方法
CN103731189B (zh) * 2014-01-08 2016-10-12 桂林电子科技大学 共形天线阵动态子阵划分方法及波达方向估计方法
CN105182285A (zh) * 2015-10-14 2015-12-23 中国电子科技集团公司第二十八研究所 一种基于声矢量二维嵌套阵列的目标测向方法
CN105403856B (zh) * 2015-10-30 2017-10-24 西安电子科技大学 基于嵌套式最小冗余阵列的波达方向估计方法
CN105445696A (zh) * 2015-12-22 2016-03-30 天津理工大学 一种嵌套l型天线阵列结构及其波达方向估计方法

Also Published As

Publication number Publication date
CN105824002A (zh) 2016-08-03

Similar Documents

Publication Publication Date Title
CN105824002B (zh) 基于嵌套式子阵阵列的波达方向估计方法
CN105403856B (zh) 基于嵌套式最小冗余阵列的波达方向估计方法
Wang et al. Nuclear norm minimization framework for DOA estimation in MIMO radar
CN106772226B (zh) 基于压缩感知时间调制阵列的doa估计方法
Yan et al. Real-valued root-MUSIC for DOA estimation with reduced-dimension EVD/SVD computation
CN103971029B (zh) 一种用于网格失配下的doa估计的交替迭代方法
CN104865556B (zh) 基于实域加权最小化l1范数方法的MIMO雷达系统DOA估计方法
CN107544051A (zh) 嵌套阵列基于k‑r子空间的波达方向估计方法
CN107907855A (zh) 一种互素阵列转化为均匀线阵的doa估计方法及装置
CN105403874A (zh) 非均匀阵列欠定波达方向估计方法
CN108802705A (zh) 一种基于稀疏的空时自适应处理方法及系统
CN109239651B (zh) 互质面阵下的二维doa跟踪方法
Wu et al. Fast 3-D imaging algorithm based on unitary transformation and real-valued sparse representation for MIMO array SAR
CN111965591A (zh) 一种基于四阶累积量矢量化dft的测向估计方法
Srinivas et al. Performance comparison of reconstruction algorithms in compressive sensing based single snapshot doa estimation
CN115236584A (zh) 基于深度学习的米波雷达低仰角估计方法
CN113655444B (zh) 一种阵元失效下基于重加权先验的mimo雷达doa估计方法
Zheng et al. DOA estimation via coarray tensor completion with missing slices
CN112733333A (zh) 互质面阵中一种基于多项式求根的二维测向估计方法
Meng et al. Efficient real-valued DOA estimation based on the trigonometry multiple angles transformation in monostatic MIMO radar
Lai et al. Improved DFT method for DOA estimation with extended coprime array: based on large difference coarray
CN112016037A (zh) 一种互质面阵中基于降维Capon求根的二维测向估计方法
Liu et al. Real-valued reweighted l 1 norm minimization method based on data reconstruction in MIMO radar
Yao et al. Degree-of-freedom strengthened cascade array for DOD-DOA estimation in MIMO array systems
Chen et al. A gridless DOA estimation algorithm based on unsupervised deep learning

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant