CN110082708A - 非均匀阵列设计和波达方向估计方法 - Google Patents
非均匀阵列设计和波达方向估计方法 Download PDFInfo
- Publication number
- CN110082708A CN110082708A CN201910139417.9A CN201910139417A CN110082708A CN 110082708 A CN110082708 A CN 110082708A CN 201910139417 A CN201910139417 A CN 201910139417A CN 110082708 A CN110082708 A CN 110082708A
- Authority
- CN
- China
- Prior art keywords
- array
- matrix
- virtual
- uniform
- obtaining
- 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
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000013461 design Methods 0.000 title claims abstract description 12
- 239000011159 matrix material Substances 0.000 claims abstract description 102
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 11
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 7
- 230000008569 process Effects 0.000 claims abstract description 4
- 239000013598 vector Substances 0.000 claims description 34
- 238000003491 array Methods 0.000 claims description 11
- 238000003786 synthesis reaction Methods 0.000 claims description 10
- 239000000203 mixture Substances 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 5
- 238000011084 recovery Methods 0.000 claims description 3
- 238000000638 solvent extraction Methods 0.000 claims description 2
- 230000008901 benefit Effects 0.000 abstract description 5
- 238000001228 spectrum Methods 0.000 abstract description 2
- 102000002274 Matrix Metalloproteinases Human genes 0.000 abstract 1
- 108010000684 Matrix Metalloproteinases Proteins 0.000 abstract 1
- 238000010276 construction Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 238000000342 Monte Carlo simulation Methods 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000012163 sequencing technique Methods 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
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)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种非均匀阵列设计和波达方向估计方法,主要解决现有技术布阵不灵活,计算复杂度较高的问题;其实现过程是:根据嵌套阵列最大自由度确定非均匀阵列首、尾位置系数;计算满足非均匀阵列差合矩阵中包含所有虚拟阵元位置的位置系数;根据最终得到非均匀阵列阵元位置;根据接收数据计算数据协方差矩阵并矢量化得到虚拟差分合成阵列接收数据r;对r去冗余、排序得到虚拟阵列接收数据进而得到其满秩矩阵v;构造线性算子,通过估计线性算子得到信号子空间;构造选择矩阵得到旋转矩阵,最后通过旋转矩阵估计波达方向;本发明在同等条件下具有阵列配置灵活,无需对数据协方差矩阵进行特征分解,也不需对整个空域角度进行谱峰搜索的优点。
Description
技术领域
本发明属于阵列信号处理技术领域,尤其涉及一种非均匀阵列设计和波达方向估计方法,通过灵活布阵,对多个空域目标源进行波达方向估计。
背景技术
在现代阵列信号处理领域中,波达方向估计(以下简称DOA估计)占有一个重要地位,传统DOA估计主要通过MUSIC方法、ESPRIT方法及其衍生方法进行。对于一个总阵元数为N的均匀线性阵列,采用传统DOA估计方法所能识别的最大空域目标个数为N-1个。在实际中,空域目标个数大于阵元数情况时有发生,采用少于目标个数的均匀线性阵列的传统DOA估计方法将失效。因此,空域目标个数大于阵元数的DOA估计被许多学者广泛研究,最近提出的嵌套阵列和互质阵列通过合理布阵构造一等效的虚拟阵列,提高阵列自由度(Degreeof freedom,DOF),通过利用虚拟阵列而不是原始阵列来进行DOA 估计。嵌套阵列是由两个或多个具有不同阵元间隔的线性子阵组成,由它构造的虚拟阵列是一个完全填充的均匀线性阵列,但其第一个子阵阵元间距通常较小(阵元排布密度较大),会引起阵元间的互耦问题。互质阵列是由两个阵元数分别为M、N的均匀线性子阵组成,其中M、N为互质整数,虽然互质阵列能获得多余阵元数的自由度DOF,但其构造的虚拟阵列不是一个完全填充的均匀线性阵列,导致在DOA估计时,只能利用互质阵列构造的虚拟阵列的一部分,这就带来了自由度和阵列孔径的降低。
由于嵌套阵列、互质阵列每个子阵均由均匀线性子阵组成,这就导致在机载雷达或物理空间受限的情况下,时常出现无法为其子阵阵元安装找到合理位置的问题,同时关于非均匀阵列提出的新算法,发明名称为基于嵌套式子阵阵列的波达方向估计方法,公开(公告)号:CN105824002A,提出的对数据协方差矩阵特征分解后利用MUSIC算法谱峰搜索进行DOA估计;以及发明名称为嵌套阵列基于K-R子空间的波达方向估计方法,公开(公告)号:CN107544051A,提出的稀疏信号重构法利用了信号谱的稀疏性来进行DOA估计;这两种方法在进行DOA估计时,通常需要巨大的计算量。本发明的目的在于克服上述已有非均匀阵列以及DOA估计方法不足,设计一种新的非均匀阵列和波达方向估计方法,该非均匀阵列无需由均匀子阵列组成,且其生成的虚拟阵列是一个完全填充的均匀线性阵列,同时为降低计算量,通过构造一线性算子得到信号子空间进行DOA估计。本发明在保证可识别的空域目标个数大于阵元数前提下,提高了布阵的灵活性,降低DOA估计运算量。
发明内容
为实现上述目的,本发明的技术思路是:根据嵌套阵列经过最优布阵获取阵列最大自由度(最多虚拟阵元个数)思想,来确定非均匀阵列第一个阵元和最后一个阵元位置系数,通过计算非均匀阵列其余阵元位置系数,构造非均匀阵列;对非均匀阵列构造的虚拟差分合成阵列去冗余、排序得到虚拟阵列;构造线性算子,根据虚拟阵列进行DOA估计。其实现方案包括如下:
1)构造非均匀阵列:
1a)根据嵌套阵列最大自由度DOF确定非均匀阵列位置系数x1和xN;
假定非均匀阵列阵元位置为:
D=[d1,d2,...,dn,...,dN]=d×[x1,x2,...,xn,...,xN]
根据嵌套阵列构造的虚拟阵列:
DULA=[-(N2(N1+1)-1)d,...,0,...,(N2(N1+1)-1)d]
获取最大自由度DOF=2N2(N1+1)-1,得到非均匀阵列位置系数则本发明阵列阵元位置为
其中d为阵元间隔,取值为入射信号最小半波长;N1与N2为嵌套阵列子阵阵元个数,N1=N2=N/2;xn表示非均匀阵列第n个阵元的位置系数,n=1,2,...,N,N为总阵元数; 元素为{1~(xN-1)}的随机递增互异整数;
1b)计算满足非均匀阵列差合矩阵中包含所有虚拟阵元位置的位置系数
为保证本发明阵列在差合处理之后包含所有的虚拟阵元位置,且与嵌套阵列自由度一致,构造一向量P=[D,D,...,D],向量P由N个D 按行排列组成,再构造另一矩阵即矩阵为非均匀阵列产生的差合矩阵:
矩阵为反对称矩阵,只研究其上三角形元素即可,将其按行排布得到一个行向量公式如下:
为使非均匀阵列位置所产生的虚拟差合阵列是完备的,即包含所有虚拟阵元位置,构造的差合矩阵所对应的中应包含个不同元素;
1c)根据位置系数获取非均匀阵列阵元位置;
在1a)中,随机从{1~(xN-1)}中选取N-2个互异递增数构成得到一组均匀阵列阵元位置D,进行迭代,当D中阵元位置满足1b)中的差合矩阵所对应中包含个不同值即可停止迭代,否则继续从1a)选取一组阵元位置,进行迭代,直至满足迭代停止条件,退出迭代,最终得到阵元位置
2)根据所设计非均匀阵列得到接收数据X,进而得到数据协方差矩阵 Rx,向量化该协方差矩阵,得到虚拟差分合成阵列接收数据r;
3)根据虚拟差分合成阵列接收数据r,对r进行去冗余、排序操作,最终得到虚拟阵列接收数据
4)根据虚拟阵列接收数据构造选择矩阵Jz对进行秩恢复操作得到满秩矩阵V;
5)根据虚拟阵列的满秩矩阵V,构造一线性算子Q,通过估计线性算子Q 得到信号子空间ES;
6)基于旋转因子不变法的思想定义两个选择矩阵Jg1和Jg2,获得ES1和ES2,通过ES1和ES2的相位关系获得旋转矩阵Φ,最后通过旋转矩阵估计波达方向
本发明与现有阵列结构相比较有如下优点:
1.本发明阵列虽与嵌套阵列自由度相同,但比嵌套布阵、互质布阵方式更灵活,不需要由两个或多个均匀线性子阵构成,且可产生完全填充型均匀线性虚拟阵列进行DOA估计,更适于工程实现;
2.本发明可满足在机载或者物理空间受限的环境下,从本发明中选取一种性能相同但却适合现场位置的非均匀阵列,进行灵活稀疏布阵,得到多余阵元数的自由度,扩大了阵列孔径,使得可估计的空域目标数远超过实际阵元数,具有更好的测角性能;
3.本发明无需对数据协方差矩阵进行特征分解,也不需通过对整个空域角度进行谱峰搜索,通过构造线性算子Q得到信号子空间ES,利用旋转因子不变法思想便可得到信号的DOA估计值,解决了在采样点数增多、空域网格步长划分较细,使得特征分解以及普峰搜索运算量大、耗时长、实时性差的问题;
4.本发明阵列为非均匀阵列,与传统均匀阵列相比,在阵列面积不变的情况下,能够有效减少阵元数目,降低成本,且阵元完全相同,更易于生产和降低成本。
附图说明
图1为本发明的三种六阵元非均匀阵列及虚拟阵列排布方式图;
图2为六阵元二阶嵌套阵列及产生虚拟阵元阵列排布方式图;
图3为本发明方法实现流程图;
图4为本发明方法空域目标DOA估计散点图;
图5为本发明和二阶嵌套阵、均匀线阵对DOA估计均方根误差与信噪比关系比较图;
具体实施方式
为了让本发明的上述和其它目的、特征及优点能更明显,下文特举本发明实施例,并配合所附图示,做详细说明如下:
参照图1为本发明三种六阵元非均匀阵列及产生虚拟阵列排布方式图,从图1可看到存在多种满足迭代条件的阵元位置系数即存在多种阵元排布方式,要说明的是,嵌套阵列排布只是本发明阵列的一个子集;
参照图2为六阵元二阶嵌套阵列及产生虚拟阵元阵列排布方式图,本发明的非均匀阵列与嵌套阵列构造的虚拟阵列是一个完全填充的均匀线性阵列,即全部位置上虚拟阵元都是连续的,二者虚拟阵列的最大自由度DOF相同。因此,本发明阵列具有嵌套阵列的优点,不仅可识别多于阵元数的目标个数而且提高了目标方向识别的空间分辨力,同时具有布阵灵活的优点。
参照图3,由于嵌套阵列构造的虚拟阵列是一个完全填充的均匀线性阵列,可提高DOA估计的自由度DOF,本发明根据嵌套阵列经过最优布阵获取阵列最大自由度即最多虚拟阵元个数思想来设计一种新的非均匀阵列,构造一线性算子,利用连续的虚拟阵列进行DOA估计。为便于叙述,假定本发明阵列与二阶嵌套阵列总阵元数均为N,N为偶数,N为奇数的情况可类比得到,本发明方法步骤如下:
步骤一、根据嵌套阵列构造的虚拟阵列获取最大自由度DOF,计算非均匀阵列阵元位置系数,得到非均匀阵列;
1a)根据嵌套阵列最大自由度DOF确定非均匀阵列位置系数x1和xN;
二阶嵌套阵为获取最大的阵列自由度,则每阶阵元数为N1=N2=N/2,其最大自由度为DOF=2N2(N1+1)-1,由其构造的虚拟阵列可等效为一个阵元位置为:
DULA=[-(N2(N1+1)-1)d,...,0,...,(N2(N1+1)-1)d] 的均匀线阵,
在总阵元数相同情况下为得到与二阶嵌套阵相同的阵列自由度,非均匀阵列位置矢量可表示为:
D=[d1,d2,...,dn,...,dN]=d×[x1,x2,...,xn,...,xN]
第一个阵元为参考阵元其位置系数x1=0,最后一个阵元位置系数则
其中d为阵元间隔,取值为入射信号的最小半波长;xn表示非均匀阵列第 n个阵元的位置系数,n=1,2,...,N,N为总阵元数; 元素为{1~(xN-1)}的随机递增互异整数;
1b)计算满足非均匀阵列差合矩阵中包含所有虚拟阵元位置的位置系数
为保证本发明阵列与嵌套阵列构造的虚拟阵列一致,则本发明阵列必须满足在差合之后包含所有的虚拟阵元位置,即是完备的,其具体实现如下:
首先构造一向量P=[D,D,...,D],构造另一个向量即矩阵为非均匀阵列产生的差合矩阵:
由于矩阵为反对称矩阵,故只研究其上三角形元素即可,将其按行排布得到一个行向量为:
为使非均匀阵列位置所产生的虚拟差合阵列是完备的,即包含所有虚拟阵元位置,构造的差合矩阵所对应的中应包含个不同元素;
1c)根据位置系数获取非均匀阵列阵元位置
在1a)中,随机从{1~(xN-1)}中选取N-2个互异递增数构成得到一组均匀阵列阵元位置D,进行迭代,当D中阵元位置满足1b)中的差合矩阵所对应中包含个不同值即可停止迭代,否则继续从1a)选取一组阵元位置,进行迭代,直至满足迭代停止条件,退出迭代,最终得到阵元位置
步骤二、根据所设计的非均匀阵列,得到阵列接收数据X,进而得到数据协方差矩阵Rx,向量化该协方差矩阵,得到虚拟差分合成阵列接收数据r;
2a)本发明阵列接收数据信号为远场窄带信号,通过L次同步采样获取数据X;
2b)根据非均匀阵列接收数据X,估计得到数据协方差矩阵Rx:
其中L表示快拍数,X∈CN×L,Rx∈CN×N,(·)H表示共轭转置;
2c)根据阵列协方差矩阵Rx,计算虚拟差分合成阵列接收数据r:
r=vec(Rx)
其中,vec(·)表示对矩阵进行向量化操作;
步骤三、根据虚拟差分合成阵列接收数据r,对r进行去冗余、排序操作,最终得到虚拟阵列接收数据
由于接收数据r中包含噪声,故无法直接对r进行去冗余、排序操作,只有通过合理设计才能得到虚拟阵列接收数据其原理是由于接收数据r元素位置与阵列产生的差合阵元位置一一对应,通过对差合阵元位置进行处理得到索引集,由索引集对r内部元素进行选取便可达到对r去冗余、排序目的,进而得到具体实现如下:
3a)根据1)得到的非均匀阵列,构造一列向量列向量由1)中差合矩阵的元素按行叠加组成一个列向量:
再构造另一向量U,其内部个元素均匀连续变化,与非均匀阵列产生的虚拟阵元排布相同:
对列向量中每一个元素与向量U的Uj进行比较,在时,记录在中的索引值i,并令Uj=N2,要说明的是,在进行比较时,只要满足便更新索引集Γ=Γ∪{i},最终得到索引集Γ,其内部包含中个虚拟阵元的索引值,所对应的元素构成一个新的向量其中
i=1,2,...,N2,
式中N为总阵元数;
3b)根据对排序得到新的索引集获取虚拟阵列接收数据
由于索引集Γ所对应中元素并不是按从小到大顺序排列的,因此通过对排序即向量与所对应F中的个元素一一对应;
至此,得到排序后的索引值通过接收数据r元素位置与阵列产生的差合阵元位置一一对应关系,用索引值选取接收数据r便可达到去冗余、排序目的,最终得到虚拟阵列接收数据
其中,sort(·)表示排序运算,F表示内部元素按顺序排列后的向量,表示F中每个元素对应索引值的集合;
步骤四、根据虚拟阵列接收数据构造选择矩阵Jz对进行秩恢复操作得到满秩矩阵V,其具体过程如下:
4a)根据虚拟阵列接收数据构造选择矩阵Jz将分割为个子矩阵:
每个子矩阵用表示,由接收数据的到行组成,其中Jz表示为:
4b)根据叠加得到满秩矩阵V,V可表示为:
其中,
步骤五、非均匀阵列接收数据快拍数L、维度增加时,为了降低由特征分解获得信号子空间的运算量,根据满秩矩阵V,构造一线性算子,通过估计线性算子Q得到信号子空间ES,其具体过程如下:
5a)列出线性算子Q与信号子空间ES关系式:
其中,IK表示K×K的单位矩阵,K为目标个数,(·)H表示共轭转置,因此只需估计出线性算子即可获得信号子空间ES;
5b)根据满秩矩阵V,估计线性算子Q,通过线性算子Q与信号子空间 ES关系式获得信号子空间ES;
为得到线性算子如下:
分割满秩矩阵
其中由V的前K行组成,由V的行组成,则线性算子能被估计如下:
获得线性算子Q后,通过便可得到信号子空间;
步骤六、基于旋转因子不变法的思想定义两个选择矩阵Jg1和Jg2,获得 ES1和ES2,通过ES1和ES2的相位关系获得旋转矩阵Φ,最后利用旋转矩阵估计波达方向其具体实现如下:
6a)定义两个选择矩阵:
其中 表示的单位矩阵,表示的全1向量,和分别表示和零矩阵;
6b)通过选择矩阵Jg1和Jg2获取信号子空间ES1和 ES2,ES1=Jg1ES,ES2=Jg2ES,基于旋转因子不变法的思想,利用信号子空间ES1和ES2关系获得对Ψ∈CK×K进行特征分解 [v,q]=eig(Ψ)得到旋转矩阵Φ∈CK×K对角线元素:
其中eig(·)表示特征分解,表示对矩阵求伪逆,diag(·)表示对角阵,v和q=[q1,q2,...,qk,...,qK]分别表示矩阵Ψ的特征向量矩阵和特征值矩阵,θk表示第k个信号的波达方向,λ为入射信号波长;
6c)通过旋转矩阵Φ估计波达方向
其中表示k个波达方向估计值,Φk,k表示矩阵Φ的第k行第k列的元素,d为入射信号半波长,arg(·)表示取相位运算。
本发明的效果可以通过以下的仿真结果进一步说明:
仿真1:利用本发明非均匀阵列能估计数量大于接收阵列阵元数的目标数,图4仿真为6阵元非均匀阵列,阵元位置矢量D=d×[0,1,3,4,9,11],阵元个数为N=6,d/λ=1/2,对K=11个远场窄带目标进行方向识别,这11个目标的波达方向为:
(-50°,-40°,-30°,-20°,-10°,0°,10°,20°,30°40°,50°)
在信噪比SNR=-10dB,快拍数为L=1024,进行50次Monte-Carlo实验;
从图4中可以看到,在低信噪比下,本发明非均匀阵列在阵元数为6的情况下能正确无误的识别出11个空域目标源,而传统六阵元均匀线阵在此情况下将不能进行目标方向识别;
仿真2:本发明非均匀阵列与二阶嵌套阵、均匀线阵在不同信噪比条件下的性能比较如图5所示,本发明非均匀阵列阵元个数为N=6,阵元位置矢量D=d×[0,1,3,4,9,11];二阶嵌套阵总阵元数为N1+N2=6,N1和N2分别为每阶子阵阵元个数,N1=N2=3,其二阶嵌套阵阵元位置矢量为 Z=d×[0,1,2,3,7,11];均匀线阵总阵元数为6,,对空间K=2个远场窄带目标进行方向识别,两个目标的波达方向为(θ1,θ2)=(20°,40°),将目标 DOA估计的结果和实际值比较,最后得到误差统计。实验信噪比SNR从0dB 到40dB,快拍数为L=1024,进行500次Monte-Carlo实验;
从图5中可以看到,当信噪比在0dB下,本发明的非均匀阵列和二阶嵌套阵DOA估计精度在0.05°附近,均匀线阵的估计精度在0.34°附近,可以得到本发明方法优于均匀线阵;当信噪比在[0,20]dB时,本发明方法的估计精度优于二阶嵌套阵方法,随着信噪比的升高,均匀线阵估计性能上升,曲线下降加快,在40dB时三者性能趋近相同;
以上所述,仅是本发明的较佳实施例而已,并非对本发明做任何形式上的限制,虽然本发明已以较佳实施例揭露如上,然而并非用以限定本发明,任何熟悉本专业的技术人员,在不脱离本发明技术方案范围内,当可利用上述揭示的技术内容做出些许更动或修饰为等同变化的等效实施例,但凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,均仍属于本发明技术方案的范围内。
Claims (6)
1.一种非均匀阵列设计和波达方向估计方法,包含以下步骤:
1)构造非均匀阵列:
1a)根据嵌套阵列最大自由度DOF,确定非均匀阵列位置系数x1和xN;
假定非均匀阵列阵元位置为:
D=[d1,d2,...,dn,...,dN]=d×[x1,x2,...,xn,...,xN]
根据嵌套阵列构造的虚拟阵列
DULA=[-(N2(N1+1)-1)d,...,0,...,(N2(N1+1)-1)d]
获取最大自由度DOF=2N2(N1+1)-1,得到非均匀阵列位置系数x1=0,则本发明阵列阵元位置为
其中d为阵元间隔,取值为入射信号的最小半波长;N1与N2为嵌套阵列子阵阵元个数,N1=N2=N/2;xn表示非均匀阵列第n个阵元的位置系数,n=1,2,...,N,N为总阵元数; 由{1~(xN-1)}的随机递增互异整数组成;
1b)计算非均匀阵列差合矩阵中包含所有虚拟阵元位置的位置系数
为保证本发明阵列在差合处理之后包含所有的虚拟阵元位置,且与嵌套阵列自由度一致,构造一向量P=[D,D,...,D],向量P由N个D按行排列组成,再构造另一矩阵 即矩阵为非均匀阵列产生的差合矩阵:
由于矩阵为反对称矩阵,故只研究其上三角元素即可,将其按行排布得到一个行向量如下:
为使非均匀阵列位置所产生的虚拟差合阵列是完备的,即包含所有虚拟阵元位置,构造的差合矩阵所对应的中应包含个不同元素;
1c)根据位置系数获取非均匀阵列阵元位置;
在1a)中,随机从{1~(xN-1)}中选取N-2个互异递增数构成得到一组均匀阵列阵元位置D,进行迭代,当D中阵元位置满足1b)中的差合矩阵所对应中包含个不同值即可停止迭代,否则继续从1a)选取一组阵元位置,进行迭代,直至满足迭代停止条件,退出迭代,最终得到阵元位置
2)根据所设计非均匀阵列得到接收数据X,进而得到数据协方差矩阵Rx,向量化该协方差矩阵,得到虚拟差分合成阵列接收数据r;
3)根据虚拟差分合成阵列接收数据r,对r进行去冗余、排序操作,最终得到虚拟阵列接收数据
4)根据虚拟阵列接收数据构造选择矩阵Jz对进行秩恢复操作得到满秩矩阵V;
5)根据虚拟阵列的满秩矩阵V,构造一线性算子Q,通过估计线性算子Q得到信号子空间ES;
6)基于旋转因子不变法的思想定义两个选择矩阵Jg1和Jg2,获得ES1和ES2,通过ES1和ES2的相位关系获得旋转矩阵Φ,最后通过旋转矩阵估计波达方向
2.根据权利要求1所述的非均匀阵列设计和波达方向估计方法,其中步骤2)中的虚拟差分合成阵列接收数据r,按如下步骤计算:
2.1)根据非均匀阵列接收数据X,估计得到数据协方差矩阵Rx:
其中,L表示快拍数,X∈CN×L,Rx∈CN×N,(·)H表示共轭转置;
2.2)根据数据协方差矩阵Rx,计算虚拟差分合成阵列接收数据r∈N2×1:
r=vec(Rx)
其中vec(·)表示对矩阵进行向量化操作。
3.根据权利要求1所述的非均匀阵列设计和波达方向估计方法,其中步骤3)中的虚拟阵列数据接收矩阵按如下步骤计算:
接收数据r中包含噪声,无法直接对r进行去冗余、排序操作,通过对差合阵元位置进行处理得到索引集,由索引集对r内部元素进行选取便可达到对r去冗余、排序目的,进而得到具体实现如下:
3.1)根据1)得到的非均匀阵列,构造一列向量列向量由1)中差合矩阵的元素按行叠加组成列向量
再构造另一向量U,其内部个元素均匀连续变化,与非均匀阵列产生的虚拟阵元排布相同:
对列向量中每一个元素与向量U的Uj进行比较,在时,记录在中的索引值i,并令Uj=N2,在进行比较时,只要满足便更新索引集Γ=Γ∪{i},最终得到索引集Γ,其内部包含中个虚拟阵元的索引值,所对应的元素构成一个新的向量
其中,i=1,2,...,N2,
N为总阵元数;
3.2)根据对排序得到新的索引集获取虚拟阵列接收数据
由于索引集Γ所对应中元素并不是按从小到大顺序排列的,因此通过对排序即向量与所对应F中的个元素一一对应;
至此,得到排序后的索引值通过接收数据r元素位置与阵列产生的差合阵元位置一一对应关系,用索引值选取接收数据r便可达到去冗余、排序目的,最终得到虚拟阵列接收数据
其中,sort(·)表示排序运算,F表示内部元素按顺序排列后的向量,表示F中每个元素对应索引值的集合。
4.根据权利要求1所述的非均匀阵列设计和波达方向估计方法,其中步骤4)得到满秩矩阵V,按如下步骤计算:
4.1)根据虚拟阵列接收数据构造选择矩阵Jz将分割为个子矩阵:
每个子矩阵用表示,由接收数据的到行组成,其中Jz表示为:
4.2)根据叠加得到满秩矩阵V,V可表示为:
其中,子阵序号
5.根据权利要求1所述的非均匀阵列设计和波达方向估计方法,其中步骤5)得到线性算子Q以及信号子空间ES,按如下步骤计算:
5.1)根据满秩矩阵V,列出线性算子Q与信号子空间ES关系式:
其中,IK表示K×K的单位矩阵,K为目标个数,(·)H表示共轭转置,因此只需估计出线性算子即可获得信号子空间ES;
5.2)根据线性算子Q与信号子空间ES关系式,通过估计线性算子Q获得信号子空间ES;
为得到线性算子分割满秩矩阵V:
其中由V的前K行组成,由V的行组成,则线性算子能被估计如下:
获得线性算子Q后,通过便可得到信号子空间。
6.根据权利要求1所述的非均匀阵列设计和波达方向估计方法,其中步骤6)得到旋转矩阵Φ,最后利用旋转矩阵Φ估计波达方向按如下步骤计算:
6.1)定义两个选择矩阵:
其中表示单位矩阵,表示的全1向量;表示和零矩阵;
6.2)通过选择矩阵Jg1和Jg2获取信号子空间ES1和ES2,得到对Ψ∈CK×K进行特征分解[v,q]=eig(Ψ)得到旋转矩阵Φ∈CK×K:
其中,信号序号k=1,2,…,K,eig(·)表示特征分解,表示对矩阵求伪逆,diag(·)表示矩阵对角化,v和q=[q1,q2,…,qk,…,qK]分别表示矩阵Ψ的特征向量矩阵和特征值矩阵,θk表示第k个信号的波达方向,λ为入射信号波长;
6.3)通过旋转矩阵Φ估计波达方向
其中表示k个波达方向估计值,Φk,k表示矩阵Φ的第k行第k列元素,d为入射信号半波长,arg(·)表示取相位运算。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910139417.9A CN110082708B (zh) | 2019-02-25 | 2019-02-25 | 非均匀阵列设计和波达方向估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910139417.9A CN110082708B (zh) | 2019-02-25 | 2019-02-25 | 非均匀阵列设计和波达方向估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110082708A true CN110082708A (zh) | 2019-08-02 |
CN110082708B CN110082708B (zh) | 2023-05-02 |
Family
ID=67413068
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910139417.9A Active CN110082708B (zh) | 2019-02-25 | 2019-02-25 | 非均匀阵列设计和波达方向估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110082708B (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110726967A (zh) * | 2019-10-25 | 2020-01-24 | 北京理工大学 | 用于一维阵列测向的单边稀疏嵌套阵设计方法 |
CN110736959A (zh) * | 2019-10-25 | 2020-01-31 | 北京理工大学 | 一种基于和差协同阵构建的平面互质阵列设计方法 |
CN111175691A (zh) * | 2019-11-29 | 2020-05-19 | 北京理工大学 | 一种用于波达方向估计的双边稀疏嵌套阵设计方法 |
CN111751798A (zh) * | 2020-07-22 | 2020-10-09 | 上海英恒电子有限公司 | 一种雷达测角方法 |
CN111896929A (zh) * | 2020-08-28 | 2020-11-06 | 西安电子科技大学 | 非均匀mimo雷达的dod/doa估计算法 |
CN111983554A (zh) * | 2020-08-28 | 2020-11-24 | 西安电子科技大学 | 非均匀l阵下的高精度二维doa估计 |
CN112327303A (zh) * | 2020-10-22 | 2021-02-05 | 四川长虹电器股份有限公司 | 天线虚拟均匀线阵数据获取方法 |
CN112632753A (zh) * | 2020-12-01 | 2021-04-09 | 厦门大学 | 宽带非频变的非均匀间隔阵列天线方向图综合方法及装置 |
CN112630724A (zh) * | 2020-10-30 | 2021-04-09 | 哈尔滨工程大学 | 一种适用于uuv平台的高分辨目标方位估计方法 |
CN113435027A (zh) * | 2021-06-23 | 2021-09-24 | 哈尔滨工程大学 | 一种高自由度低耦合稀疏线性阵列排布方法 |
CN113589224A (zh) * | 2021-08-03 | 2021-11-02 | 宜宾电子科技大学研究院 | 一种基于增强嵌套阵的doa估计方法 |
CN114397619A (zh) * | 2022-01-04 | 2022-04-26 | 西安电子科技大学 | 基于非均匀稀疏阵列二维定位算法 |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002048853A (ja) * | 2000-08-02 | 2002-02-15 | Matsushita Electric Ind Co Ltd | 電波到来方向推定装置及び指向性可変送受信装置 |
US20050195103A1 (en) * | 2004-01-13 | 2005-09-08 | Davis Dennis W. | Phased arrays exploiting geometry phase and methods of creating such arrays |
US20080231505A1 (en) * | 2007-03-23 | 2008-09-25 | Weiqing Zhu | Method of Source Number Estimation and Its Application in Method of Direction of Arrival Estimation |
CN101799535A (zh) * | 2009-11-27 | 2010-08-11 | 西安电子科技大学 | Mimo雷达目标方向的估计方法 |
CN102608565A (zh) * | 2012-03-23 | 2012-07-25 | 哈尔滨工程大学 | 一种基于均匀圆阵列的波达方向估计方法 |
US20130249731A1 (en) * | 2012-03-26 | 2013-09-26 | Fujitsu Ten Limited | Radar apparatus and target detecting method |
CN103344939A (zh) * | 2013-06-14 | 2013-10-09 | 西安交通大学 | 一种非相干及相干混合信号的二维波达方向估计方法 |
US20130300596A1 (en) * | 2012-05-11 | 2013-11-14 | Fujitsu Limited | Detection and ranging apparatus and ranging method |
CN103399291A (zh) * | 2013-07-22 | 2013-11-20 | 西安电子科技大学 | 基于快速稀疏恢复的超分辨波达方向估计方法 |
CN105403856A (zh) * | 2015-10-30 | 2016-03-16 | 西安电子科技大学 | 基于嵌套式最小冗余阵列的波达方向估计方法 |
CN105824002A (zh) * | 2016-04-15 | 2016-08-03 | 西安电子科技大学 | 基于嵌套式子阵阵列的波达方向估计方法 |
CN106019213A (zh) * | 2016-05-09 | 2016-10-12 | 电子科技大学 | 一种部分稀疏l阵及其二维doa估计方法 |
CN108663653A (zh) * | 2018-05-17 | 2018-10-16 | 西安电子科技大学 | 基于l形电磁矢量传感器阵列的波达方向估计方法 |
-
2019
- 2019-02-25 CN CN201910139417.9A patent/CN110082708B/zh active Active
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002048853A (ja) * | 2000-08-02 | 2002-02-15 | Matsushita Electric Ind Co Ltd | 電波到来方向推定装置及び指向性可変送受信装置 |
US20050195103A1 (en) * | 2004-01-13 | 2005-09-08 | Davis Dennis W. | Phased arrays exploiting geometry phase and methods of creating such arrays |
US20080231505A1 (en) * | 2007-03-23 | 2008-09-25 | Weiqing Zhu | Method of Source Number Estimation and Its Application in Method of Direction of Arrival Estimation |
CN101799535A (zh) * | 2009-11-27 | 2010-08-11 | 西安电子科技大学 | Mimo雷达目标方向的估计方法 |
CN102608565A (zh) * | 2012-03-23 | 2012-07-25 | 哈尔滨工程大学 | 一种基于均匀圆阵列的波达方向估计方法 |
US20130249731A1 (en) * | 2012-03-26 | 2013-09-26 | Fujitsu Ten Limited | Radar apparatus and target detecting method |
US20130300596A1 (en) * | 2012-05-11 | 2013-11-14 | Fujitsu Limited | Detection and ranging apparatus and ranging method |
CN103344939A (zh) * | 2013-06-14 | 2013-10-09 | 西安交通大学 | 一种非相干及相干混合信号的二维波达方向估计方法 |
CN103399291A (zh) * | 2013-07-22 | 2013-11-20 | 西安电子科技大学 | 基于快速稀疏恢复的超分辨波达方向估计方法 |
CN105403856A (zh) * | 2015-10-30 | 2016-03-16 | 西安电子科技大学 | 基于嵌套式最小冗余阵列的波达方向估计方法 |
CN105824002A (zh) * | 2016-04-15 | 2016-08-03 | 西安电子科技大学 | 基于嵌套式子阵阵列的波达方向估计方法 |
CN106019213A (zh) * | 2016-05-09 | 2016-10-12 | 电子科技大学 | 一种部分稀疏l阵及其二维doa估计方法 |
CN108663653A (zh) * | 2018-05-17 | 2018-10-16 | 西安电子科技大学 | 基于l形电磁矢量传感器阵列的波达方向估计方法 |
Non-Patent Citations (3)
Title |
---|
CHENGWEI ZHOU ET AL.: "Direction-of-Arrival Estimation for Coprime Array", 《IEEE TRANSACTIONS ON SIGNAL PROCESSING》 * |
张利强: "最小冗余线阵的DOA估计", 《火力与指挥控制》 * |
董健: "一种获取大阵元数低冗余度线阵排列的方法", 《电子学报》 * |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110726967B (zh) * | 2019-10-25 | 2021-08-03 | 北京理工大学 | 用于一维阵列测向的单边稀疏嵌套阵设计方法 |
CN110736959A (zh) * | 2019-10-25 | 2020-01-31 | 北京理工大学 | 一种基于和差协同阵构建的平面互质阵列设计方法 |
CN110726967A (zh) * | 2019-10-25 | 2020-01-24 | 北京理工大学 | 用于一维阵列测向的单边稀疏嵌套阵设计方法 |
CN110736959B (zh) * | 2019-10-25 | 2021-07-09 | 北京理工大学 | 一种基于和差协同阵构建的平面互质阵列设计方法 |
CN111175691A (zh) * | 2019-11-29 | 2020-05-19 | 北京理工大学 | 一种用于波达方向估计的双边稀疏嵌套阵设计方法 |
CN111751798A (zh) * | 2020-07-22 | 2020-10-09 | 上海英恒电子有限公司 | 一种雷达测角方法 |
CN111896929A (zh) * | 2020-08-28 | 2020-11-06 | 西安电子科技大学 | 非均匀mimo雷达的dod/doa估计算法 |
CN111983554A (zh) * | 2020-08-28 | 2020-11-24 | 西安电子科技大学 | 非均匀l阵下的高精度二维doa估计 |
CN111896929B (zh) * | 2020-08-28 | 2023-08-04 | 西安电子科技大学 | 非均匀mimo雷达的dod/doa估计算法 |
CN112327303A (zh) * | 2020-10-22 | 2021-02-05 | 四川长虹电器股份有限公司 | 天线虚拟均匀线阵数据获取方法 |
CN112630724A (zh) * | 2020-10-30 | 2021-04-09 | 哈尔滨工程大学 | 一种适用于uuv平台的高分辨目标方位估计方法 |
CN112632753A (zh) * | 2020-12-01 | 2021-04-09 | 厦门大学 | 宽带非频变的非均匀间隔阵列天线方向图综合方法及装置 |
CN112632753B (zh) * | 2020-12-01 | 2023-08-04 | 厦门大学 | 宽带非频变的非均匀间隔阵列天线方向图综合方法及装置 |
CN113435027A (zh) * | 2021-06-23 | 2021-09-24 | 哈尔滨工程大学 | 一种高自由度低耦合稀疏线性阵列排布方法 |
CN113589224A (zh) * | 2021-08-03 | 2021-11-02 | 宜宾电子科技大学研究院 | 一种基于增强嵌套阵的doa估计方法 |
CN114397619A (zh) * | 2022-01-04 | 2022-04-26 | 西安电子科技大学 | 基于非均匀稀疏阵列二维定位算法 |
Also Published As
Publication number | Publication date |
---|---|
CN110082708B (zh) | 2023-05-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110082708A (zh) | 非均匀阵列设计和波达方向估计方法 | |
CN108872929B (zh) | 基于内插虚拟阵列协方差矩阵子空间旋转不变性的互质阵列波达方向估计方法 | |
CN107092004B (zh) | 基于信号子空间旋转不变性的互质阵列波达方向估计方法 | |
CN108896954B (zh) | 互质阵中一种基于联合实值子空间的波达角估计方法 | |
CN110109050B (zh) | 嵌套阵列下基于稀疏贝叶斯的未知互耦的doa估计方法 | |
CN107544051A (zh) | 嵌套阵列基于k‑r子空间的波达方向估计方法 | |
CN110244272B (zh) | 基于秩一去噪模型的波达方向估计方法 | |
CN109375152B (zh) | 电磁矢量嵌套l阵下低复杂度的doa与极化联合估计方法 | |
CN113032721B (zh) | 一种低计算复杂度的远场和近场混合信号源参数估计方法 | |
CN110837076A (zh) | 一种基于张量分解的矢量水听器阵列方位估计方法 | |
CN111965591A (zh) | 一种基于四阶累积量矢量化dft的测向估计方法 | |
CN111983554A (zh) | 非均匀l阵下的高精度二维doa估计 | |
CN111352063A (zh) | 一种均匀面阵中基于多项式求根的二维测向估计方法 | |
CN113567913A (zh) | 基于迭代重加权可降维的二维平面doa估计方法 | |
CN113673317A (zh) | 基于原子范数最小化可降维的二维离格doa估计方法 | |
CN110174658B (zh) | 基于秩一降维模型和矩阵补全的波达方向估计方法 | |
CN106980105B (zh) | 电磁矢量传感器阵列空间旋转解相干测向方法 | |
CN112255629A (zh) | 基于联合uca阵列的序贯esprit二维不相干分布源参数估计方法 | |
CN113835063B (zh) | 一种无人机阵列幅相误差与信号doa联合估计方法 | |
CN113075610B (zh) | 一种基于互质极化阵列的差分阵列内插的doa估计方法 | |
CN112731280A (zh) | 互质阵列混合噪声环境下的esprit-doa估计方法 | |
CN108872930B (zh) | 扩展孔径二维联合对角化doa估计方法 | |
CN111368256A (zh) | 一种基于均匀圆阵的单快拍测向方法 | |
Yang et al. | A correlation-aware sparse Bayesian perspective for DOA estimation with off-grid sources | |
CN114648041A (zh) | 一种基于平行稀疏阵列的二维欠定doa估计算法 |
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 |