CN109696651B - 一种基于m估计的低快拍数下波达方向估计方法 - Google Patents
一种基于m估计的低快拍数下波达方向估计方法 Download PDFInfo
- Publication number
- CN109696651B CN109696651B CN201910086425.1A CN201910086425A CN109696651B CN 109696651 B CN109696651 B CN 109696651B CN 201910086425 A CN201910086425 A CN 201910086425A CN 109696651 B CN109696651 B CN 109696651B
- Authority
- CN
- China
- Prior art keywords
- matrix
- estimation
- covariance matrix
- trace
- arrival
- 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.)
- Expired - Fee Related
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)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明提出一种基于M估计的低快拍数下波达方向估计方法,包括:步骤1.设置天线阵列并获得所述天线阵列的观测数据矢量;步骤2.计算观测数据矢量在快拍数为T的协方差矩阵;步骤3.对所述样本协方差矩阵进行迹正则化得到迹正则化后的矩阵;步骤4.对迹正则化后的矩阵进行托普利兹校正,得到目标矩阵;步骤5.根据迹正则化后的目标矩阵,计算矩阵收缩估计的收缩系数;步骤6.根据所述收缩目标矩阵和所述收缩系数,构造新的协方差矩阵;步骤7.将新的协方差矩阵重新置为样本协方差矩阵,重复步骤2~步骤6,直到得到最终的估计;步骤8.利用得到的协方差矩阵采用子空间分解的方法计算信号的波达方向。
Description
技术领域
本发明属于无线移动通信技术领域,具体涉及一种基于M估计的低快拍数下波达方向估计方法。
背景技术
阵列信号处理经过近几十年的不断发展,已广泛应用于通信、雷达、声呐、电子对抗等许多领域。波达方向估计技术主要研究如何估计空间目标信号的波达方向(Directionof Arrival,简称DOA),是阵列信号处理领域重要的研究内容之一。无论是在军事还是在民用领域,空间目标信号的波达方向都是极为重要的估计参数,是后续定位、探测、导航和成像等诸多技术的先决基础,因此研究波达方向估计技术有着十分重要的现实意义。
随着信息技术的飞速发展,天线阵列在接收目标发射的信号是受到的干扰日趋多样化和复杂化。在这种背景下,信号接收时是面临的噪声常常不再符合高斯白噪声的特性,而经常呈现为存在“拖尾现象”的复椭圆对称分布。另外,随着现在阵列中天线数目的不断增加,保持大量的快拍数会导致较高的计算复杂度和硬件系统实现的复杂度,因此,如何在小快拍数下保持良好的DOA估计精度也是目前阵列信号处理中的一项重要研究。
传统的DOA估计算法,如MUSIC算法等算法,需要对接收信号的协方差矩阵进行子空间分解来实现。但是在实际应用中,该协方差矩阵是难以直接获得的,因此一般的做法是利用接收信号的样本协方差矩阵代替真实协方差矩阵来进行子空间分解,但是在非高斯噪声的条件下,样本协方差矩阵并不能良好地近似真实的接收信号协方差矩阵,因此基于子空间分解的算法在非高斯噪声的背景下性能会严重降低。最近,一种基于M估计的协方差矩阵估计理论被提出,这种估计方法被证明可以有效对抗存在离群值的“拖尾噪声”。另外,在快拍数不足的情况下,对样本协方差矩阵进行子空间分解得到的特征值和特征向量与真实值偏差很大,因此基于子空间分解的DOA估计方法也不能很好地工作。为了解决这个问题,学者们陆续提出了G-MUSIC算法、矩阵收缩估计方法和矩阵托普利兹化的方法,但是它们各有缺陷:G-MUSIC算法只对空间邻近信号有效、矩阵收缩估计在DOA估计中没有合适的收缩目标、矩阵托普利兹化在高信噪比时性能不佳。
另外,在非高斯噪声和快拍数较少两种情况同时存在时,学者们也进行了一定的研究,代表的成果有稳健G-MUSIC算法和将矩阵收缩估计与M估计结合的方法,但是前者对空间分离信号没有性能提升,适用范围很窄,后者选取的目标矩阵为单位矩阵,不适用于DOA估计领域。因此,本问题目前仍没有得到很好的解决。
发明内容
鉴于以上所述现有技术的缺点,本发明的目的在于提供一种基于M估计的低快拍数下波达方向估计方法,以达到准确实现DOA估计的目的。
为实现上述目的及其他相关目的,本发明提供一种基于M估计的低快拍数下波达方向估计方法,该估计方法包括:
步骤1.设置天线阵列并获得所述天线阵列的观测数据矢量X(t);
步骤4.对迹正则化后的矩阵R1进行托普利兹校正,得到目标矩阵RT;
步骤5.根据迹正则化后的目标矩阵RT,计算矩阵收缩估计的收缩系数ρ;
步骤6.根据所述收缩目标矩阵RT和所述收缩系数ρ,构造新的协方差矩阵R3;
步骤8.利用得到的协方差矩阵R3采用子空间分解的方法计算信号的波达方向;
所述收缩系数ρ通过以下方法计算得到:
所述新的协方差矩阵R3通过以下方法计算得到:
其中,xi代表观测数据矢量X(t)的第i列,其代表阵列在第i次快拍中的全部观测数据。
可选地,所述观测数据矢量X(t)的计算方法为:
X(t)=AS(t)+N(t)
其中,X(t)=[x0(t),…,xM-1(t)]T为观测数据矢量,S(t)=[s1(t),…,sK(t)]T为信号源数据矢量,A=[a(θ1),…,a(θK)]为阵列流型矩阵,a(θk)=[a0,k,…,aM-1,k]T为导向矢量,N(t)=[n0(t),…,nM-1(t)]T为噪声矢量。
其中,X(tn)表示在第tn个快拍时的阵列观测数据矢量。
可选地,所述迹正则化后的矩阵R1通过以下方法计算得到:
其中,trace表示计算矩阵的迹,即矩阵对角线元素的和。
可选地,所述目标矩阵RT通过以下方法计算得到:
其中,Jm代表M维的移位矩阵,Jm中只有第m条对角线元素为1,其余元素为0。
可选地,所述步骤8包括以下子步骤:
由MUSIC算法,利用所述信号子空间Us和所述噪声子空间Un的正交性得到空间谱函数为:
其中,θ为谱峰搜索时的角度,a(θ)代表阵列的导向矢量,Un代表对协方差矩阵进行特征值分解后,由较小的M-K个特征值对应的特征向量构成的矩阵,
如上所述,本发明的一种基于M估计的低快拍数下波达方向估计方法,具有以下有益效果:
本发明使用托普利兹校正的方法克服了快拍数不足导致DOA估计精度下降的问题;在构造协方差矩阵时使用了M估计的方法减轻了非高斯噪声对DOA估计精度的影响;同时将托普利兹校正后的矩阵作为目标矩阵,M估计得到的矩阵作为原始矩阵,利用协方差矩阵收缩估计的方法构造新的协方差矩阵,克服了托普利兹校正带来的在DOA估计在高信噪比下性能不佳的问题;综上,本发明可以在小快拍数和非高斯噪声同时存在的条件下实现准确的DOA估计。
附图说明
为了进一步阐述本发明所描述的内容,下面结合附图对本发明的具体实施方式作进一步详细的说明。应当理解,这些附图仅作为典型示例,而不应看作是对本发明的范围的限定。
图1为本发明阵列设置示意图;
图2为本发明实施方式仿真实验DOA的均方根误差在高斯白噪声下随SNR变化关系示意图;
图3为本发明实施方式仿真实验DOA的均方根误差在Student-t分布的噪声下随SNR变化关系示意图;
图4为本发明实施方式仿真实验DOA的均方根误差在威布尔分布的噪声下随SNR变化关系示意图;
图5为本发明实施方式仿真实验DOA的均方根误差在K分布的噪声下随SNR变化关系示意图;
图6为本发明实施例的一种基于M估计的低快拍数下波达方向估计方法的流程图。
具体实施方式
以下通过特定的具体实例说明本发明的实施方式,本领域技术人员可由本说明书所揭露的内容轻易地了解本发明的其他优点与功效。本发明还可以通过另外不同的具体实施方式加以实施或应用,本说明书中的各项细节也可以基于不同观点与应用,在没有背离本发明的精神下进行各种修饰或改变。需说明的是,在不冲突的情况下,以下实施例及实施例中的特征可以相互组合。
需要说明的是,以下实施例中所提供的图示仅以示意方式说明本发明的基本构想,遂图式中仅显示与本发明中有关的组件而非按照实际实施时的组件数目、形状及尺寸绘制,其实际实施时各组件的型态、数量及比例可为一种随意的改变,且其组件布局型态也可能更为复杂。
如图6所示,本发明提供一种基于M估计的低快拍数下波达方向估计方法,包括以下步骤:
步骤1:设置天线阵列并获得所述天线阵列的观测数据矢量X(t)。
具体地,设置一个如图1所示的一维均匀线阵,共有M个阵元,阵列间距d=λ/2。该阵列接收K个波长为λ的窄带平稳独立且非零峰度信号,K≤M-1,阵列采样的快拍数T略大于阵元数M。在本实施方式中阵元数M设置为20,阵列的快拍数T设置为25。本实施方式有K=3两个信号源,其入射角度分别为θ1=-3°,θ2=16°,θ3=55°。阵列阵元上的噪声为复椭圆对称分布的非高斯噪声,且噪声独立于信号。
第m个阵元在t时刻接收到的数据可以表示为:
其中,sk(t)表示第k个信号源,nm(t)表示第m个阵元上的噪声,而am,k表示第m个阵元对第k个信号的阵列流型元素,可以进一步表示为:
其中,ωk有下列形式:
其中,λ是信号波长,θk是第k个信号的DOA。
将阵列观测数据表示成矩阵形式为:
X(t)=AS(t)+N(t)
其中,X(t)=[x0(t),…,xM-1(t)]T为观测数据矢量,S(t)=[s1(t),…,sK(t)]T为信号源数据矢量,A=[a(θ1),…,a(θK)]为阵列流型矩阵,a(θk)=[a0,k,…,aM-1,k]T为导向矢量,N(t)=[n0(t),…,nM-1(t)]T为噪声矢量。
步骤2:设置天线阵列并获得所述天线阵列的观测数据矢量X(t);
其中,T为快拍数,表示在第tn个快拍时的阵列观测数据矢量。
其中,trace代表计算矩阵的迹,即矩阵对角线元素的和。
步骤4:对迹正则化后的矩阵R1进行托普利兹校正,得到目标矩阵RT;
阵列接收数据的真实协方差矩阵是托普利兹矩阵,但当实际接收信号的快拍数较少时,观测数据的样本协方差矩阵不再满足矩阵的托普利兹特性,因此对步骤3中得到的迹正则化后的矩阵R1进行托普利兹校正,得到目标矩阵RT,以减轻低快拍数带来的影响,计算公式为:
其中,Jm代表M维的移位矩阵,Jm中只有第m条对角线元素为1,其余元素为0。
步骤5:根据迹正则化后的矩阵,计算矩阵收缩估计的收缩系数ρ;
由于托普利兹化使阵列协方差矩阵的真实元素发生了变化,因此信噪比较高时,利用托普利兹化后的矩阵计算的DOA的精度会比直接利用样本协方差矩阵计算的DOA的精度低。为了解决这个问题,本发明采用矩阵收缩估计的方法,将托普利兹校正后的目标矩阵RT作为收缩目标矩阵,以提升在高信噪比下的性能,采用收缩估计的方法需要计算矩阵收缩系数,本发明采用“障碍估计”的原则计算收缩系数ρ,其计算公式为:
步骤6:根据所述收缩目标矩阵RT和所述收缩系数ρ,构造新的协方差矩阵R3;
本发明采用M估计的思想重构样本协方差矩阵以对抗非高斯噪声,当阵列的观测矩阵X(t)已知时,M估计的表达式为:
其中xi代表观测数据X(t)的第i列,其代表阵列在第i次快拍中的全部观测数据。u要求是非负、非增、连续的函数。本发明采用Tyler提出的估计法,即取u(x)=1/x。另外,在表达式的左右两边都出现了协方差矩阵R,因此该方法需要一个初始的协方差矩阵R作为初值,然后采取迭代的方式进行求解。另外,为了同时对抗小快拍数带来的影响,本发明同时采用协方差矩阵收缩估计的思路,将M估计得到的协方差矩阵作为收缩估计的主矩阵,将步骤4中托普利兹化后得到的矩阵RT作为收缩目标矩阵,构造新的协方差矩阵,其计算公式为:
其中,收缩系数ρ由步骤5中的计算得到。
如步骤6所述,M估计是一个迭代的过程,因此在构造出新的协方差矩阵后,我们仍需要对新构造的协方差矩阵R3进行迭代更新以取得更精确的协方差矩阵估计结果。具体而言,将步骤5中得到的协方差矩阵R3重新置为样本协方差重复步骤2~步骤6,直到得到最终的估计R3,在本发明中,迭代次数设置为100。
步骤8:利用得到的协方差矩阵R3采用子空间分解的方法计算信号的波达方向。
对步骤7所得新构造的协方差矩阵R3进行特征分解,得到其噪声子空间Un和信号子空间Us。并利用最大的3(K)个特征值所对应的特征向量uz1,uz2,uz3作为列建立特征向量矩阵Uzs={uz1,uz2,uz3}作为信号子空间,利用其余17(M-K)个小特征值对应的特征向量uz4,...,uz20为列建立特征向量矩阵Uzn={uz4,...,uz20}作为噪声子空间,3个大特征值组成对角矩阵Ds,其余17个小特征值组成对角矩阵Dn,如下式所示:
再利用信号子空间和噪声子空间的正交性得到MUSIC算法的空间谱函数为:
其中,θ为谱峰搜索时的角度。
设置第一次谱峰搜索的步长为l1=1,由MUSIC谱函数得到3个峰值,设为θ1',θ’2,θ’3;设置第二次搜索范围为θ1'±2,θ’2±2,θ’3±2,搜索步长为l2=0.1,得到最终的DOA的估计值
为验证算法的性能,本实施方式设计四组仿真实验。仿真选取了四种典型的噪声分布场景,一种是高斯白噪声,其余3种均属于在阵列信号处理和雷达领域常见的非高斯噪声,分布为Student-t分布的噪声、威布尔分布的噪声和K分布的噪声。其中Student-t分布的自由度设置为2.5;威布尔分布的尺度参数设置为1,形状参数设置为0.5;K分布的尺度参数设置为2,形状参数设置为0.25。在仿真实验中,将建议方法与5种其它的方法的性能进行了对比,其余的方法分别为:经典MUSIC算法、仅使用了M估计的方法、稳健G-MUSIC算法、同时结合托普利兹化与M估计的方法和将M估计与矩阵收缩估计结合的方法。
两组实验的随机实验次数均为500,所用均匀线阵相同,阵元数为20,阵列快拍数为25,阵列间距为d=λ/2。入射的信号源数K=3。信号参数分别为θ1=-3°,θ2=16°,θ3=55°。四组实验的结果分别如图2、3、4、5所示。由图3、4、5可知,在非高斯分布的噪声存在的条件下,本发明所提方案在信噪比较低的条件下比其它的方法估计精度更高,同时在信噪比较高的条件下不会出现托普利兹化带来的精度损失问题。
上述实施例仅例示性说明本发明的原理及其功效,而非用于限制本发明。任何熟悉此技术的人士皆可在不违背本发明的精神及范畴下,对上述实施例进行修饰或改变。因此,举凡所属技术领域中具有通常知识者在未脱离本发明所揭示的精神与技术思想下所完成的一切等效修饰或改变,仍应由本发明的权利要求所涵盖。
Claims (6)
1.一种基于M估计的低快拍数下波达方向估计方法,其特征在于,该估计方法包括:
步骤1.设置天线阵列并获得所述天线阵列的观测数据矢量X(t);
步骤4.对迹正则化后的矩阵R1进行托普利兹校正,得到目标矩阵RT;
步骤5.根据迹正则化后的目标矩阵RT,计算矩阵收缩估计的收缩系数ρ;
步骤6.根据所述目标矩阵RT和所述收缩系数ρ,构造新的协方差矩阵R3;
步骤8.利用得到的协方差矩阵R3采用子空间分解的方法计算信号的波达方向;
所述收缩系数ρ通过以下方法计算得到:
所述新的协方差矩阵R3通过以下方法计算得到:
其中,xi代表观测数据矢量X(t)的第i列,其代表阵列在第i次快拍中的全部观测数据。
2.根据权利要求1所述的一种基于M估计的低快拍数下波达方向估计方法,其特征在于,所述观测数据矢量X(t)的计算方法为:
X(t)=AS(t)+N(t)
其中,X(t)=[x0(t),…,xM-1(t)]T为观测数据矢量,S(t)=[s1(t),…,sK(t)]T为信号源数据矢量,A=[a(θ1),…,a(θK)]为阵列流型矩阵,a(θk)=[a0,k,…,aM-1,k]T为导向矢量,N(t)=[n0(t),…,nM-1(t)]T为噪声矢量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910086425.1A CN109696651B (zh) | 2019-01-29 | 2019-01-29 | 一种基于m估计的低快拍数下波达方向估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910086425.1A CN109696651B (zh) | 2019-01-29 | 2019-01-29 | 一种基于m估计的低快拍数下波达方向估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109696651A CN109696651A (zh) | 2019-04-30 |
CN109696651B true CN109696651B (zh) | 2021-03-26 |
Family
ID=66234556
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910086425.1A Expired - Fee Related CN109696651B (zh) | 2019-01-29 | 2019-01-29 | 一种基于m估计的低快拍数下波达方向估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109696651B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111208471B (zh) * | 2020-03-02 | 2023-01-13 | 重庆大学 | 少快拍非线极化电磁波波达方向估计方法 |
CN115792792B (zh) * | 2022-11-22 | 2023-10-27 | 中国人民解放军空军预警学院 | 极低样本下基于协方差矩阵迭代估计的测向方法与系统 |
CN116500625B (zh) * | 2023-06-29 | 2023-10-20 | 天津知海科技有限公司 | 恢复成像方法、装置、系统、电子设备及可读存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104166136A (zh) * | 2014-07-11 | 2014-11-26 | 河海大学 | 一种基于干扰子空间跟踪的高效自适应单脉冲测角方法 |
CN106054123A (zh) * | 2016-06-06 | 2016-10-26 | 电子科技大学 | 一种稀疏l阵及其二维doa估计方法 |
CN106569181A (zh) * | 2016-11-15 | 2017-04-19 | 大连大学 | 基于协方差矩阵重构稳健Capon波束形成的算法 |
CN107329108B (zh) * | 2017-05-03 | 2019-07-23 | 浙江大学 | 基于内插虚拟阵列协方差矩阵Toeplitz化重建的互质阵列波达方向估计方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105022026B (zh) * | 2015-07-08 | 2017-07-07 | 陕西理工学院 | L型阵列的二维到达角估计方法 |
-
2019
- 2019-01-29 CN CN201910086425.1A patent/CN109696651B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104166136A (zh) * | 2014-07-11 | 2014-11-26 | 河海大学 | 一种基于干扰子空间跟踪的高效自适应单脉冲测角方法 |
CN106054123A (zh) * | 2016-06-06 | 2016-10-26 | 电子科技大学 | 一种稀疏l阵及其二维doa估计方法 |
CN106569181A (zh) * | 2016-11-15 | 2017-04-19 | 大连大学 | 基于协方差矩阵重构稳健Capon波束形成的算法 |
CN107329108B (zh) * | 2017-05-03 | 2019-07-23 | 浙江大学 | 基于内插虚拟阵列协方差矩阵Toeplitz化重建的互质阵列波达方向估计方法 |
Non-Patent Citations (4)
Title |
---|
A covariance matrix shrinkage method with Toeplitz rectified target for DOA estimation under the uniform linear array;Yanyan Liu等;《International Journal of Electronics and Communications》;20171231;第50-55页 * |
Regularized -Estimators of Scatter Matrix;Esa Ollila等;《IEEE TRANSACTIONS ON SIGNAL PROCESSING》;20141115;第62卷(第22期);第6059-6070页 * |
Shrinkage Algorithms for MMSE Covariance Estimation;Yilun Chen等;《IEEE TRANSACTIONS ON SIGNAL PROCESSING》;20101031;第58卷(第10期);第5016-5029页 * |
大维随机矩阵理论在阵列信号参数估计中的应用;刘妍妍;《中国博士学位论文全文数据库 信息科技辑》;20180315;第2-6、73-77页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109696651A (zh) | 2019-04-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110109050B (zh) | 嵌套阵列下基于稀疏贝叶斯的未知互耦的doa估计方法 | |
CN109061554B (zh) | 一种基于空间离散网格动态更新的目标到达角度估计方法 | |
CN106054123A (zh) | 一种稀疏l阵及其二维doa估计方法 | |
CN109696651B (zh) | 一种基于m估计的低快拍数下波达方向估计方法 | |
CN109946643B (zh) | 基于music求解的非圆信号波达方向角估计方法 | |
CN109959891A (zh) | 电磁矢量l阵中空间角与极化参数的降维谱峰搜索方法 | |
CN112130111B (zh) | 一种大规模均匀十字阵列中单快拍二维doa估计方法 | |
CN107450047A (zh) | 嵌套阵下基于未知互耦信息的压缩感知doa估计方法 | |
CN113835063B (zh) | 一种无人机阵列幅相误差与信号doa联合估计方法 | |
CN108919178A (zh) | 一种基于对称嵌套阵列的混合场信源定位方法 | |
CN103353588B (zh) | 基于天线均匀平面阵的二维波达方向角估计方法 | |
CN111046591B (zh) | 传感器幅相误差与目标到达角度的联合估计方法 | |
CN104020440B (zh) | 基于l型干涉式线性阵列的二维波达角估计方法 | |
CN110018438B (zh) | 一种波达方向估计方法和装置 | |
CN112379327A (zh) | 一种基于秩损估计的二维doa估计与互耦校正方法 | |
CN113567913B (zh) | 基于迭代重加权可降维的二维平面doa估计方法 | |
CN114814830B (zh) | 一种基于鲁棒主成分分析降噪的米波雷达低仰角测高方法 | |
CN110244272A (zh) | 基于秩一去噪模型的波达方向估计方法 | |
CN110286350A (zh) | 一种l型稀疏阵doa估计的精确配对方法及装置 | |
CN110286351A (zh) | 一种基于l型嵌套阵的二维doa估计方法及装置 | |
CN110275131A (zh) | 一种基于虚拟差分阵列的doa跟踪方法及装置 | |
CN115356678A (zh) | 基于dpnalm算法的稀疏阵列doa估计方法 | |
CN111366893A (zh) | 一种均匀圆阵未知互耦条件下的非圆信号方位角估计方法 | |
CN105572629B (zh) | 一种适用于任意阵列结构的低运算复杂度的二维测向方法 | |
CN112763972B (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210326 Termination date: 20220129 |
|
CF01 | Termination of patent right due to non-payment of annual fee |