CN111352063B - 一种均匀面阵中基于多项式求根的二维测向估计方法 - Google Patents
一种均匀面阵中基于多项式求根的二维测向估计方法 Download PDFInfo
- Publication number
- CN111352063B CN111352063B CN201911321920.2A CN201911321920A CN111352063B CN 111352063 B CN111352063 B CN 111352063B CN 201911321920 A CN201911321920 A CN 201911321920A CN 111352063 B CN111352063 B CN 111352063B
- Authority
- CN
- China
- Prior art keywords
- finding
- root
- matrix
- polynomial
- follows
- 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
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 59
- 239000011159 matrix material Substances 0.000 claims abstract description 53
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 8
- 230000014509 gene expression Effects 0.000 claims description 17
- 238000013178 mathematical model Methods 0.000 claims description 8
- 239000000654 additive Substances 0.000 claims description 4
- 230000000996 additive effect Effects 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 4
- 239000013598 vector Substances 0.000 claims description 4
- 238000003491 array Methods 0.000 claims description 3
- 238000001228 spectrum Methods 0.000 description 8
- 230000003595 spectral effect Effects 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000004088 simulation Methods 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
-
- 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/74—Multi-channel systems specially adapted for direction-finding, i.e. having a single antenna system capable of giving simultaneous indications of the directions of different signals
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- 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
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种均匀面阵中基于多项式求根的二维测向估计方法,该方法首先由均匀面阵的接收信号得到协方差矩阵;然后对协方差矩阵进行特征值分解并得到信号子空间和噪声子空间,根据方向矩阵和噪声子空间的正交关系确定求根多项式;最后求出多项式的根,完成参数配对,完成二维角度参数估计。本发明的可以充分平衡复杂度和角度估计性能,突破传统二维角度估计方法中,角度估计性能良好但复杂度较高或复杂度较低但角度估计性能一般的局限;并且本发明能够实现较高分辨率的二维DOA估计,角度估计性能优于2D‑PM算法和2D‑ESPRIT算法,同2D‑MUSIC算法基本一致,算法复杂度远低于2D‑MUSIC算法。
Description
技术领域
本发明属于阵列信号处理技术,具体涉及一种均匀面阵中基于多项式求根的二维测向估计方法。
背景技术
早期的空间谱估计研究均基于均匀线阵,但是现实的信源都处在三维空间中,只有通过仰角和方位角才能确定信源的入射位置。虽然均匀线阵具有结构简单的优点,但是仅仅只研究一维的DOA估计显然是具有很大局限性的,所以对多维参数估计的研究具有非常现实的意义。近年来,学者们提出了多种对二维空间谱估计及其多普勒频率的联合估计的方法。二维空间谱估计的常用方法主要包括2D-MUSIC算法、2D-ESPRIT算法、2D-PM算法和PARAFAC技术等。二维空间谱估计中,采用的阵型主要包括L型阵列、面阵、圆阵和平行阵列等,平面阵列利用子空间法实现空间谱估计可以获得较高分辨率。同时,相对于一维线阵,平面阵列中阵元的排布方式也极大缩小了阵列尺寸。
MUSIC算法的提出是在空间谱估计发展的历史上具有里程碑的意义。该算法主要通过噪声子空间和阵列流形之间的正交关系构造谱函数,然后用其极值来实现信源参数估计。其中2D-MUSIC算法是二维DOA估计的典型算法,此方法可以产生渐进无偏估计,但是要在二维参数空间搜索谱峰,所以其计算复杂度非常高。本发明算法通过求根的方式,无需谱峰搜索,获得求根结果后,配对完成信源角度参数估计。在有效降低复杂度的同时保证了角度估计性能。
发明内容
发明目的:为实现均匀面阵中不相干多信源入射条件下进行信号角度参数估计,本发明提供一种基于多项式求根的二维测向估计方法。
技术方案:一种均匀面阵中基于多项式求根的二维测向估计方法,包括如下步骤:
(1)建立均匀面阵接收信号的数学模型;
(2)构建接收信号的协方差矩阵,对协方差矩阵进行特征值分解获得噪声子空间;
(3)由噪声矩阵和方向向量的正交关系确定求根多项式;
(4)基于求根多项式进行求根计算并完成参数配对,确定信源角度。
进一步的,所述均匀面阵接收信号X表达式如下:
所述均匀面阵接收信号的数学模型表达式为:
X=[Ay⊙Ax]S+N=AS+N
式中,均匀面阵共有M×N个阵元,xN(t)表示子阵N的接收信号,AxΦN-1表示子阵N的方向矩阵,nN(t)为子阵N的加性高斯白噪声,Ay表示按Y轴纵向扫描接收信号矩阵,S表示信源矩阵,⊙表示Khatri-Rao积。
进一步的,步骤(2)包括如下过程:
获取L个快拍得到均匀面阵的协方差矩阵的估计如下:
对信号协方差矩阵进行特征值分解,其表达式如下:
其中,Es表示信号子空间,En表示噪声子空间,Ds和Dn均表示对角矩阵。
进一步的,步骤(3)根据方向矩阵和噪声子空间的关系构造关于u和v的求根多项式包括如下计算过程:
求根MUSIC多项式表达式如下:
其中z1=ej2πdu/λ,z2=ej2πdv/λ;
其中根据MUSIC算法构造求根多项式表达式如下:
或者表示如下:
式中,u=sin θ sin φv=cosθ sin φ且θ和φ分别对应信源角度的仰角和方位角,V(u,v)表示二维求根MUSIC多项式,ay(u)和ax(v)表示均匀面阵的阵列流形,(·)H表示该矩阵的共轭转置矩阵。
进一步的,步骤(4)包括如下过程:
(a)确定信源角度对应的二项式根
取Q(u)=0在单位圆内具有最大幅值的K个根的相位获得u的估计,取Q(v)=0在单位圆内具有最大幅值的K个根/>的相位获得v的估计,其计算表达式如下:
(b)完成参数配对,完成二维DOA估计:
通过对和/>进行配对来完成对应角度估计,所构造代价函数表达式如下:
其中表示与/>张成空间相互正交的投影矩阵,/>是根据x轴上M个阵元对应的方向矩阵为Ax和y轴上N个阵元对应的方向矩阵为Ay的估计值经过Khatri-Rao积构造的方向矩阵的估计值;
配对完成后得到二维波达角的估计表达式如下:
其中,和/>分别表示/>和/>中第i和第j个元素。
与现有技术相比,本发明具有四个方面的显著效果:
(1)本发明所提供的估计方法需对信源角度进行谱峰搜索,与2D-MUSIC算法对比,复杂度较低;
(2)该算法将二维求根问题转化为两次一维求根,有效降低了求根复杂度和求根难度;
(3)该算法的角度估计性能优于2D-PM算法和2D-ESPRIT算法,同2D-MUSIC算法的角度估计性能基本一致;
(4)能有效地用于二维DOA估计,同时获得较高精度的角度估计结果。
附图说明
图1是均匀面阵示意图;
图2是本发明所述方法在SNR=5dB时的估计点阵图;
图3是本发明所述方法在SNR=20dB时的估计点阵图;
图4是本发明所述方法的角度估计性能在不同快拍数条件下的对比图;
图5是本发明算法角度估计性能在不同阵元下的曲线图;
图6是本发明所述方法和2D-MUSIC算法、2D-ESPRIT算法、2D-PM算法、RD-MUSIC算法的角度估计性能在相同阵列结构和相同快拍数条件下的对比图。
具体实施方式
为了详细的说明本发明所公开的技术方案,下面结合附图以及具体实施例对本发明的技术方案做进一步的表述。
在本发明中所述的方法中,对于矩阵及其计算表达式中,(·)T,(·)H和(·)-1分别表示为转置,共轭转置和求逆。X表示矩阵,x(·)表示矢量,代表Kronecker积,⊙表示Khatri-Rao积,Rank(·)表示求矩阵的秩,det(·)表示求行列式的值,angle(·)表示求复数的相角。
为实现在均匀面阵中不相干多信源入射条件下进行信号角度参数估计,本发明所述的一种均匀面阵中基于多项式求根的二维测向估计方法的基本思路是:由阵列信号的数学模型确定接收信号的协方差矩阵,通过对协方差矩阵进行特征值分解获得信号子空间和噪声子空间,然后由方向矩阵和噪声子空间的关系构建求根多项式,最后对多项式进行求根并完成参数配对,完成信源角度估计。
本发明所述方法的具体的实施步骤如下:
如图1所示。假设空间有K个窄带远场不相干信源入射到此均匀面阵上,其二维波达方向为(θk,φk),k=1,2,…,K,其中θk,φk分别代表第k个信源的仰角和方位角。首先根据阵列信号数学模型得到噪声子空间,然后根据方向矩阵和噪声子空间的关系构造关于u和v的求根多项式,最后求根并完成参数配对,得出信源信号角度参数估计值。本实施例提供的一种均匀面阵中基于多项式求根的测向估计方法的具体实现如下:
步骤1:建立均匀面阵接收信号的数学模型:
图1所示的均匀面阵,该面阵共有M×N个阵元,均匀分布,两个相邻元素的间距是d,d≤λ/2(λ是波长)。x轴和y轴上信源的方向矢量分别如下:
x轴上M个阵元对应的方向矩阵为Ax=[ax(θ1,φ1),ax(θ2,φ2),…,ax(θK,φK)],具体表示为:
y轴上N个阵元对应的方向矩阵为Ay=[ay(θ1,φ1),ay(θ2,φ2),…,ay(θK,φK)],具体表示为:
面阵中子阵1的接收信号为:
x1(t)=AxS+n1(t) (5)
式中:Ax=[ax(θ1,φ1),ax(θ2,φ2),…,ax(θK,φK)]是子阵1的方向矩阵,n1(t)是子阵1的加性高斯白噪声。是信源矩阵。
第N个子阵的接收信号为:
xN(t)=AxΦN-1S+nN(t) (6)
式中:nN(t)是第n个子阵的加性高斯白噪声。可得整个面阵的接收信号为:
式(7)中的信号也可以用如下数学模型表示:
X=[Ay⊙Ax]S+N=AS+N (8)
其中A=Ay⊙Ax。
根据Khatri-Rao积的定义,方向矩阵A可以表示为:
步骤2:求出协方差矩阵并对/>进行特征值分解,确定噪声子空间:
获取L个快拍得到均匀面阵的协方差矩阵的估计为:
对信号协方差矩阵进行特征值分解,可以表示为:
其中,Es表示信号子空间,En表示噪声子空间,Ds和Dn均表示对角矩阵。
步骤3:根据方向矩阵和噪声子空间的关系构造关于u和v的求根多项式:
2D-MUSIC算法中空间谱函数可表示为:
定义
则,阵列流形可以表示为:
ay(θk,φk)=ay(u)=[1,ej2πdu/λ,…,,ej2π(N-1)du/λ]T (13)
ax(θk,φk)=ax(v)=[1,ej2πdv/λ,…,,ej2π(M-1)dv/λ]T (14)
根据MUSIC算法构造求根多项式:
V(u,v)也可表示为:
或
对其中的u和v的集合进行估计,使方向矩阵在噪声子空间方向上投影最小。等价于
V(u,v)=0 (18)
由于
则
上式即转换为
令,
阵列流形可以改写为:
ay(u)=[1,ej2πdu/λ,…,,ej2π(N-1)du/λ]T=[1,z1,…,z1 N-1]T=ay(z1) (24)
ax(v)=[1,ej2πdv/λ,…,,ej2π(M-1)dv/λ]T=[1,z2,…,z2 M-1]T=ax(z2) (25)
为了消除u*和v*的幂次项,且只对单位圆上的z值感兴趣,所以可以用 代替/>用/>代替/>这就给出了求根MUSIC多项式,即
步骤4:求根,并确定信源角度对应的根:
由于Q(u)和Q(v)均为偶次多项式,它的根相对于单位圆为镜像对。其中,取Q(u)=0在单位圆内具有最大幅值的K个根的相位获得u的估计,取Q(v)=0在单位圆内具有最大幅值的K个根/>的相位获得v的估计。
步骤5:完成参数配对,完成二维DOA估计:
由于对和/>的估计是分开的,所以需要对/>和/>进行配对来完成对应角度估计。构造代价函数,代价函数为:
其中表示与/>张成空间相互正交的投影矩阵。而/>是根据式(3)和式(4)得到的估计值进行Khatri-Rao积构造的方向矩阵的估计值。任取/>和/>其中一个为固定值,不妨令/>固定,则估计值有K!种组合,带入式(30)得到最小值组合就是正确的配对。配对完成后,我们可以得到二维波达角的估计分别为:
其中,和/>分别表示/>和/>中第i和第j个元素,且它们已配对完成。
本发明的方法对于运算实现的复杂度分析如下:
对本发明算法的运算复杂度进行分析,具体如下:由于均匀面阵共有M×N个阵元,信源数为K,快拍数为L,所以本算法的主要复杂度包括:计算接收信号的协方差矩阵需要O{(MN)2L},特征值分解需要O{(MN)3},高阶多项式求根需要O{(2N(M-1))3+(2M(N-1))3+2(MN+1)(K-K)},配对过程需要O{2K3}所以本发明算法的总复杂度为O{(MN)2L+(MN)3+(2N(M-1))3+(2M(N-1))3+2(MN+1)(MN-K)+2K3},由于DOA估计算法的复杂度主要来自于谱峰搜索,由于无需对信源角度进行角度搜索,所以在相同阵列结构下,本算法的复杂度要远低于MUSIC算法。
图2-图3是本发明算法角度估计性能在不同信噪比下的估计点阵图。其中,图2的SNR=5dB,图3的SNR=20dB。由图2-图3表明,该算法可以精确的估计出仰角和方位角,且算法的角度估计性能随着信噪比增加变得更好。其中,入射信号的角度参数为(θ1,φ1)=(20°,30°),(θ2,φ2)=(40°,50°),均匀面阵大小为M=N=6,快拍数L=100。
图4是本发明算法角度估计性能在不同快拍下的曲线图。快拍数增加,即采样数据增多。由图可以得出,算法的角度估计性能随着快拍数增加变得更好。其中,入射信号的角度参数为(θ1,φ1)=(20°,30°),(θ2,φ2)=(40°,50°),均匀面阵大小为M=N=6。
图5是本发明算法角度估计性能在不同阵元下的曲线图。阵元数增加,即分集增益增加。由图可以得出,算法的角度估计性能随着阵元数增加变得更好。其中,入射信号的角度参数(θ1,φ1)=(20°,30°),(θ2,φ2)=(40°,50°),快拍数L=200。
图6是本发明算法和2D-MUSIC算法、2D-PM算法、2D-ESPRIT算法仿真对比结果。由图6表明,本发明算法的角度估计性能优于2D-PM算法、2D-ESPRIT算法且和2D-MUSIC算法的角度估计性能基本一致。其中,入射信号的角度参数(θ1,φ1)=(20°,30°),(θ2,φ2)=(40°,50°),均匀面阵大小为M=N=6,快拍数L=200。
本发明可以充分平衡复杂度和角度估计性能,突破传统二维角度估计方法中,角度估计性能良好但复杂度较高或复杂度较低但角度估计性能一般的局限;本发明所包含的算法能够实现较高分辨率的二维DOA估计,角度估计性能优于2D-PM算法、2D-ESPRIT算法和RD-MUSIC算法,同2D-MUSIC算法基本一致,算法复杂度远低于2D-MUSIC算法。
Claims (3)
1.一种均匀面阵中基于多项式求根的二维测向估计方法,其特征在于,包括如下步骤:
(1)建立均匀面阵接收信号的数学模型;
(2)构建接收信号的协方差矩阵,对协方差矩阵进行特征值分解获得噪声子空间;
(3)由噪声矩阵和方向向量的正交关系确定求根多项式;
包括如下计算过程:
求根MUSIC多项式表达式如下:
其中z1=ej2πdu/λ,z2=ej2πdv/λ;
对于上式,其中根据MUSIC算法构造求根多项式表达式如下:
或者表示如下:
式中,u=sinθsinφ v=cosθsinφ且θ和φ分别对应信源角度的仰角和方位角,V(u,v)表示二维求根MUSIC多项式,ay(u)和ax(v)表示均匀面阵的阵列流形,(·)H表示该矩阵的共轭转置矩阵;
(4)基于求根多项式进行求根计算并完成参数配对,确定信源角度;
包括如下过程:
(a)确定信源角度对应的二项式根
取Q(u)=0在单位圆内具有最大幅值的K个根的相位获得u的估计,取Q(v)=0在单位圆内具有最大幅值的K个根/>的相位获得v的估计,其计算表达式如下:
(b)完成参数配对,完成二维DOA估计:
通过对和/>进行配对来完成对应角度估计,所构造代价函数表达式如下:
其中表示与/>张成空间相互正交的投影矩阵,/>是根据x轴上M个阵元对应的方向矩阵为Ax和y轴上N个阵元对应的方向矩阵为Ay的估计值经过Khatri-Rao积构造的方向矩阵的估计值;
配对完成后得到二维波达角的估计表达式如下:
其中,和/>分别表示/>和/>中第i和第j个元素。
2.根据权利要求1所述的均匀面阵中基于多项式求根的二维测向估计方法,其特征在于,
所述均匀面阵接收信号X表达式如下:
所述均匀面阵接收信号的数学模型表达式为:
X=[Ay ⊙Ax]S+N=AS+N
式中,均匀面阵共有M×N个阵元,xN(t)表示子阵N的接收信号,AxΦN-1表示子阵N的方向矩阵,nN(t)为子阵N的加性高斯白噪声,Ay表示按Y轴纵向扫描接收信号矩阵,S表示信源矩阵,⊙表示Khatri-Rao积。
3.根据权利要求1所述的均匀面阵中基于多项式求根的二维测向估计方法,其特征在于,步骤(2)包括如下过程:
获取L个快拍得到均匀面阵的协方差矩阵的估计如下:
对信号协方差矩阵进行特征值分解,其表达式如下:
其中,Es表示信号子空间,En表示噪声子空间,Ds和Dn均表示对角矩阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911321920.2A CN111352063B (zh) | 2019-12-20 | 2019-12-20 | 一种均匀面阵中基于多项式求根的二维测向估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911321920.2A CN111352063B (zh) | 2019-12-20 | 2019-12-20 | 一种均匀面阵中基于多项式求根的二维测向估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111352063A CN111352063A (zh) | 2020-06-30 |
CN111352063B true CN111352063B (zh) | 2023-08-01 |
Family
ID=71192198
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911321920.2A Active CN111352063B (zh) | 2019-12-20 | 2019-12-20 | 一种均匀面阵中基于多项式求根的二维测向估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111352063B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113219399B (zh) * | 2020-08-05 | 2023-03-10 | 哈尔滨工业大学(威海) | 基于全实值计算的远场窄带无线电信号波达方向估计方法 |
CN113283055B (zh) * | 2021-04-23 | 2024-07-23 | 南京航空航天大学 | 一种基于平行因子模型的多信号分离与测向联合处理方法 |
CN113777554A (zh) * | 2021-08-26 | 2021-12-10 | 南京航空航天大学 | 一种基于求根Capon的二维DOA估计方法 |
CN113960525B (zh) * | 2021-10-15 | 2024-04-12 | 南京航空航天大学 | 一种基于频域toeplitz矩阵重构的跳频信号快速测向方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108872928A (zh) * | 2018-03-19 | 2018-11-23 | 长江大学 | 基于四线性分解的平面阵列角度的估算方法 |
CN108896954B (zh) * | 2018-06-07 | 2022-01-11 | 河海大学 | 互质阵中一种基于联合实值子空间的波达角估计方法 |
CN109613504A (zh) * | 2018-12-07 | 2019-04-12 | 中国人民解放军空军工程大学 | 一种稀疏线性阵列的快速角度估计方法 |
CN109683151A (zh) * | 2019-02-01 | 2019-04-26 | 哈尔滨工程大学 | 非均匀噪声环境下基于矩阵补全的酉求根music角度估计方法 |
CN110031794B (zh) * | 2019-04-16 | 2021-11-09 | 中国人民解放军国防科技大学 | 一种基于差分共性阵重构的相干信源doa估计方法 |
-
2019
- 2019-12-20 CN CN201911321920.2A patent/CN111352063B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN111352063A (zh) | 2020-06-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111352063B (zh) | 一种均匀面阵中基于多项式求根的二维测向估计方法 | |
CN107037392B (zh) | 一种基于压缩感知的自由度增加型互质阵列波达方向估计方法 | |
CN110749858B (zh) | 一种基于多项式求根的展开互质阵测向估计方法 | |
CN106054123A (zh) | 一种稀疏l阵及其二维doa估计方法 | |
CN109597020A (zh) | 一种使用互质线阵进行低复杂度角度估计的方法 | |
CN109375152B (zh) | 电磁矢量嵌套l阵下低复杂度的doa与极化联合估计方法 | |
Zhang et al. | Two-dimensional direction of arrival estimation for coprime planar arrays via polynomial root finding technique | |
CN107092007A (zh) | 一种虚拟二阶阵列扩展的波达方向估计方法 | |
CN109946643B (zh) | 基于music求解的非圆信号波达方向角估计方法 | |
CN109581276A (zh) | 一种基于求和求差嵌套阵的doa估计方法 | |
CN103353588B (zh) | 基于天线均匀平面阵的二维波达方向角估计方法 | |
CN110286350A (zh) | 一种l型稀疏阵doa估计的精确配对方法及装置 | |
CN110954859A (zh) | 基于l型阵列的二维非相干分布式非圆信号参数估计方法 | |
CN106997037A (zh) | 声矢量传感器阵列空间旋转解相干到达角估计方法 | |
CN112255629B (zh) | 基于联合uca的序贯esprit二维不相干分布源参数估计方法 | |
CN110837076A (zh) | 一种基于张量分解的矢量水听器阵列方位估计方法 | |
CN112269172A (zh) | 一种基于张量结构的嵌套mimo雷达角度估计方法和装置 | |
CN111983554A (zh) | 非均匀l阵下的高精度二维doa估计 | |
CN111965591A (zh) | 一种基于四阶累积量矢量化dft的测向估计方法 | |
CN110286351A (zh) | 一种基于l型嵌套阵的二维doa估计方法及装置 | |
CN106980105B (zh) | 电磁矢量传感器阵列空间旋转解相干测向方法 | |
CN112711000B (zh) | 基于最小化准则的电磁矢量互质面阵张量功率谱估计方法 | |
CN112733333A (zh) | 互质面阵中一种基于多项式求根的二维测向估计方法 | |
CN110824416B (zh) | 基于秩损原理的一维非相干分布式非圆信号参数估计方法 | |
CN113721184B (zh) | 基于改进music算法的近场信号源定位方法 |
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 |