CN109946643A - 基于music求解的非圆信号波达方向角估计方法 - Google Patents
基于music求解的非圆信号波达方向角估计方法 Download PDFInfo
- Publication number
- CN109946643A CN109946643A CN201910203570.3A CN201910203570A CN109946643A CN 109946643 A CN109946643 A CN 109946643A CN 201910203570 A CN201910203570 A CN 201910203570A CN 109946643 A CN109946643 A CN 109946643A
- Authority
- CN
- China
- Prior art keywords
- array
- vector
- covariance
- matrix
- nested
- 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
Classifications
-
- 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
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明提出了一种基于MUSIC求解的非圆信号波达方向角估计方法,用于解决现有技术中存在的非圆信号波达方向角估计精度较低的技术问题,实现步骤为:1)构建嵌套阵列c;2)获取嵌套阵列c的输出信号Y(t);3)计算嵌套阵列c的虚拟均匀阵列协方差向量和椭圆协方差向量4)构造扩展协方差矩阵G;5)计算扩展协方差矩阵G的噪声子空间Un;6)构造目标波达方向角θ的阵列流型矩阵Q(θ);7)基于MUSIC求解的方法,计算目标波达方向角θ的空间谱P(θ);8)绘制目标波达方向角θ和其空间谱P(θ)的幅度谱图。本发明可用于对飞机、舰船运动目标无源定位。
Description
技术领域
本发明属于信号处理技术领域,涉及一种非圆信号波达方向角估计方法,特别涉及一种基于MUSIC求解的非圆信号波达方向角估计的方法,可用于对飞机、舰船运动目标无源定位。
背景技术
在现代通信、雷达和卫星系统中,常用的二相相移键控、M进制幅移键控以及幅度调制等调制信号因其非圆信号特性而受到广大专家学者的关注。非圆信号波达方向角是非圆信号的重要特征参数之一,非圆信号波达方向角估计是指利用天线阵列对电磁空间中具有非圆特性的声学信号、电磁信号进行感应接收,再运用现代信号处理方法快速准确的估计出信号源的方向,在雷达、声纳、无线通信等领域具有重要应用价值。非圆信号波达方向角估计精度是评价估计方法的一项重要指标。针对非圆信号波达方向角估计构建的模型中,应用较为广泛的是基于线性均匀阵列的模型,之前的大部分非圆信号波达方向角的估计都是利用该模型生成的。
目前关于非圆信号波达方向角的估计方法是利用多项式运算的方式估计非圆信号波达方向角,这种方法旨在构建波达角选择矩阵的多项式方程,通过求解多项式方程估计非圆信号波达方向角,与传统的非圆信号波达方向角的估计方法相比,大大提高了非圆信号波达方向角的估计精度。例如申请公布号为CN107300686A,名称为“基于多项式求解的非圆信号波达方向角的估计方法”的专利申请中,公开了一种基于多项式求解的非圆信号波达方向角的估计方法。该方法利用嵌套阵列的输出信号构造出两个波达角选择矩阵,经过计算得到四个子噪声矩阵,利用子噪声矩阵构造多项式方程并求解。该方法的估计精度取决于其构造的目标函数方程,利用方程的求解得到估计的角度,在一定程度上提高了估计精度,但是由于其在构造目标函数方程时,将目标函数的最小值直接设为零,而在实际应用中目标函数的最小值不一定为零,因此估计精度往往达不到预期。
发明内容
本发明的目的在于针对上述现有技术存在的不足,提出一种基于MUSIC求解的非圆信号波达方向角估计方法,用于解决现有技术中非圆信号波达方向角的估计精度较低的技术问题。
为实现上述目的,本发明采取的技术方案包括如下步骤:
(1)构建嵌套阵列c:
(1a)以天线接收机作为阵元,构建包括M个阵元的第一均匀线性阵列a和包括N个阵元且相邻阵元的间距为(M+1)d的第二均匀线性阵列b,M≥1,N≥1,a中相邻阵元的间距为d,b中相邻阵元的间距为(M+1)d,0<d<λ/2,λ为入射到阵元的窄带信号波长,并定义a中各阵元的位置依次为D(1)=1,D(2)=2,…,D(M)=M,定义b中各阵元的位置依次为D(M+1)=M+1,D(M+2)=2(M+1),…,D(M+N)=N(M+1);
(1b)根据第一均匀线性阵列a和第二均匀线性阵列b,构造嵌套阵列c;
(2)获取嵌套阵列c的输出信号Y(t):
将K个非圆目标信号入射到嵌套阵列c中的每一个阵元,K≥1,并对c的输出信号进行采样,然后对所有采样信号滤波,得到嵌套阵列c的输出信号Y(t):
Y(t)=[y1(t),…,yi(t),…,yM+N(t)]T
其中,yi(t)表示c中第i个阵元的输出信号,1≤i≤M+N,t表示离散时间,1≤t≤L,L表示信号在时域上采样点数,[·]T表示矩阵[·]的转置运算;
(3)计算嵌套阵列c的虚拟均匀阵列协方差向量和椭圆协方差向量
(3a)计算等效协方差向量rd和等效椭圆协方差向量rs;
(3b)计算虚拟均匀阵列协方差向量和虚拟均匀阵列椭圆协方差向量
(4)构造扩展协方差矩阵G;
(5)计算扩展协方差矩阵G的噪声子空间Un:
对扩展协方差矩阵G特征分解,得到G的特征值矩阵Λ和G的特征值所对应的特征向量矩阵U,并对Λ中的特征值按照从大到小的顺序进行排列,然后将(L1+L2-K)个较小特征值在U中对应的特征向量作为G的噪声子空间Un;
(6)构造目标波达方向角θ的阵列流型矩阵Q(θ):
(6a)将嵌套阵列c[0°,180°]的观测方向范围按0°到180°的递增方向均匀划分成E+1个子观测方向区域,E≥M+N,并将第e个子观测方向区域和第e+1个子观测方向区域之间的边界射线与0°射线之间的夹角作为第e个目标波达方向角θe,其中e=1,2,…,E,定义E个目标波达方向角θ:
θ=[θ1,θ2,…,θe,…,θE]
(6b)构造目标波达方向角θ的阵列流型向量α1(θ)和阵列流型向量α2(θ),其中θe的阵列流型向量α1(θe)和阵列流型向量α2(θe)的表达式分别为:
其中,j为虚数单位;
(6c)构造目标波达方向角θ的阵列流型矩阵Q(θ),其中θe的阵列流型矩阵Q(θe)的表达式为:
(7)基于MUSIC求解的方法,计算目标波达方向角θ的空间谱P(θ),其中θe的空间谱P(θe)计算公式为:
其中,det[·]表示矩阵[·]的求行列式运算;
(8)绘制目标波达方向角θ与空间谱P(θ)的幅度谱图:
以目标波达方向角θ为x轴坐标,以空间谱P(θ)的幅度值为y轴坐标,绘制幅度谱图,从该幅度谱图中按照幅度值从高到低的顺序寻找幅度值较大的前K个谱峰,K个谱峰的峰值点所对应的x轴坐标即为非圆信号波达方向角的估计值。
本发明与现有技术相比,具有以下优点:
本发明对嵌套阵列输出的扩展协方差矩阵进行特征值分解,由于分解得到的噪声子空间和阵列流型矩阵正交的性质,信号入射方向会出现极小值,这样有效地利用矩阵计算将噪声干扰和目标信号分离出来。通过构造空间谱函数确定方向角的估计精度,并绘制空间谱幅度谱图,然后对空间谱幅度谱图进行谱峰搜索,确认波达方向角的估计值。避免了现有技术中,构造目标函数方程时最小值选取对精度的影响,有效地提高了非圆信号波达方向角的估计精度。
附图说明
图1是本发明的实现流程图;
图2是本发明中嵌套阵列的结构示意图;
图3是本发明中E个波达方向角的角度定义图。
具体实施方式
下面结合附图和具体实施例,对本发明作进一步详细描述:
参照图1、一种基于MUSIC求解的非圆信号波达方向角估计方法,包括如下步骤:
步骤1)构建嵌套阵列c:
步骤1a)以天线接收机作为阵元,构建包括M个阵元的第一均匀线性阵列a和包括N个阵元且相邻阵元的间距为(M+1)d的第二均匀线性阵列b,M≥1,N≥1,a中相邻阵元的间距为d,b中相邻阵元的间距为(M+1)d,0<d<λ/2,λ为入射到阵元的窄带信号波长,并定义a中各阵元的位置依次为D(1)=1,D(2)=2,…,D(M)=M,定义b中各阵元的位置依次为D(M+1)=M+1,D(M+2)=2(M+1),…,D(M+N)=N(M+1);
步骤1b)将b的头部与α的尾部进行拼接,且b中位置为D(M+1)=M+1的阵元与a中位置为D(M)=M的阵元之间的距离为d,得到嵌套阵列c,c的结构如图2所示,使用这样的阵列结构模型,解决了一般的线性均匀阵列模型中估计的信号数目低于阵元数目的问题,增加了阵列可识别信源数目,具体实施例中M=10,N=10,a中相邻阵元的间距d=2米。
步骤2)获取嵌套阵列c的输出信号Y(t):
将K个非圆目标信号入射到嵌套阵列c中的每一个阵元,K≥1,并对c的输出信号进行采样,然后对所有采样信号滤波,得到嵌套阵列c的输出信号Y(t):
Y(t)=[y1(t),…,yi(t),…,yM+N(t)]T
其中,yi(t)表示c中第i个阵元的输出信号,1≤i≤M+N,t表示离散时间,1≤t≤L,L表示信号在时域上采样点数,[·]T表示矩阵[·]的转置运算;
信号在时域上采样点数L在阵列信号处理中称之为快拍数,理论上来说快拍数越大,目标波达方向角的估计精度就会越高。但考虑到大快拍数会增加阵列系统的软硬件复杂度,提高了系统成本,因此合理选取适应系统的快拍数才能有效地提高目标波达方向角的估计精度,具体实施例中上=500。在利用MUSIC求解估计目标波达方向角时,需要对入射到嵌套阵列c中的非圆目标信号数量已知,具体实施例中,K=3;
步骤3)计算嵌套阵列c的虚拟均匀阵列协方差向量和椭圆协方差向量
步骤3a)根据嵌套阵列c的输出信号Y(t),计算等效协方差矩阵Rd和等效椭圆协方差矩阵Rs并将Rd中的各列元素按照由小到大的列号进行拼接,得到等效协方差向量rd,同时将Rs中的各列元素按照由小到大的列号进行拼接,得到椭圆协方差向量rs,
Rd和Rs计算公式分别为:
其中,[·]H表示矩阵[·]的共轭转置运算,[·]T表示矩阵[·]的共轭转置运算;
rd和rs表达式分别为:
rd=[Rd(1,1),Rd(2,1),…,Rd(M+N,1),Rd(1,2),…,Rd(M+N,2),…,Rd(1,M+N),…,Rd(M+N,M+N)]T
rs=[Rs(1,1),Rs(2,1),…,Rs(M+N,1),Rs(1,2),…,Rs(M+N,2),…,Rs(1,M+N),…,Rs(M+N,M+N)]T
其中,Rd(i,j)表示协方差矩阵Rd中位于第i行,第j列的元素,1≤i≤M+N,1≤j≤M+N,M+N为嵌套阵列c的阵元数,Rs(i,j)表示椭圆协方差矩阵Rs中位于第i行,第j列的元素;
步骤3b)计算等效协方差向量rd中所有元素的维数Ei,j和等效椭圆协方差向量rs中所有元素的维数Fi,j,并删除rd中维数相同的元素和维数不连续的元素后,将剩余元素按维数从小到大进行排列,得到元素个数为Cd的虚拟均匀阵列协方差向量同时删除rs中维数相同的元素和维数不连续的元素后,将剩余元素按维数从小到大进行排列,得到元素个数为Cs的虚拟均匀阵列椭圆协方差向量其中,Cd=2MN+2N-1,Cs=MN+M+N;
Ei,j和Fi,j的计算公式分别为:
Ei,j=D(j)-D(i)
Fi,j=D(j)+D(i)
其中,D(i)表示嵌套阵列c中第i个阵元的位置,D(j)表示嵌套阵列c中第j个阵元的位置,1≤i≤M+N,1≤j≤M+N;
步骤4)构造扩展协方差矩阵G:
根据虚拟均匀阵列协方差向量和椭圆协方差向量构造扩展协方差矩阵G;
其中,L1和L2表示中间变量,L1=(Cd+1)/2,L2=Cs+1-(Cd+1)/2,且有L1>L2,(·)*表示向量的共轭运算;
步骤5)计算扩展协方差矩阵G的噪声子空间Un:
对扩展协方差矩阵G特征分解,得到G的特征值矩阵Λ和G的特征值所对应的特征向量矩阵U,并对Λ中的特征值按照从大到小的顺序进行排列,排列后的特征值矩阵中前K个较大特征值对应的特征向量被称为信号子空间,Λ中剩余的(L1+L2-K)个较小特征值对应的特征向量被称为噪声子空间Un;
扩展协方差矩阵G特征分解,其公式为:
G=U.Λ.UH
其中,人为扩展协方差矩阵G的特征值矩阵,U为矩阵G的特征值所对应的特征向量矩阵;
步骤6)构造目标波达方向角θ的阵列流型矩阵Q(θ):
步骤6a)将嵌套阵列c[0°,180°]的观测方向范围按0°到180°的递增方向均匀划分成E+1个子观测方向区域,E≥M+N,并将第e个子观测方向区域和第e+1个子观测方向区域之间的边界射线与0°射线之间的夹角作为第e个目标波达方向角θe,其中e=1,2,…,E,定义E个目标波达方向角θ:
θ=[θ1,θ2,…,θe,…,θE]
E个目标波达方向角的角度定义方法如图3所示,在对嵌套阵列c的观测方向范围的进行均匀划分时,子观测方向区域的数量E+1的选择会影响对目标方向角估计的精度,在确定E+1时通常会根据期望达到的角度估计精度进行设定,E+1取值越大对目标方向角估计的精度就越高,具体实施例中E=180;
步骤6b)构造目标波达方向角θ的阵列流型向量α1(θ)和阵列流型向量α2(θ),其中θe的阵列流型向量α1(θe)和阵列流型向量α2(θe)的表达式分别为:
其中,j为虚数单位;
步骤6c)构造目标波达方向角θ的阵列流型矩阵Q(θ),其中θe的阵列流型矩阵Q(θe)的表达式为:
步骤7)基于MUSIC求解的方法,计算目标波达方向角θ的空间谱P(θ),其中θe的空间谱P(θe)计算公式为:
其中,det[·]表示矩阵[·]的求行列式运算;
MUSIC求解的基本思想是对任意阵列输出信号的协方差矩阵进行特征分解,得到与信号分类相对应的信号子空间和信号分量相正交的噪声子空间后,利用信号子空间和噪声子空间的正交性构造空间谱函数,通过谱峰搜索,对目标信号波达角方向进行估计。本步骤利用步骤5)中得到的扩展协方差矩阵G的噪声子空间Un和步骤6)中的得到的目标波达方向角θ的阵列流型矩阵Q(θ),计算出目标波达方向角θ的空间谱P(θ)。为后续的谱峰搜索以及目标信号波达角方向估计值的确定提供了数据基础;
步骤8)绘制目标波达方向角θ与空间谱P(θ)的幅度谱图:
以目标波达方向角θ为横轴坐标,以空间谱P(θ)的幅度值为纵轴坐标,绘制幅度谱图,从该幅度谱图中按照幅度值从高到低的顺序寻找幅度值较大的前K个谱峰,K个谱峰的峰值点所对应的横轴坐标即为非圆信号波达方向角的估计值;
在绘制θ与P(θ)的幅度谱图的过程中,将每个目标波达方向角θe的角度值作为横坐标并将它对应的空间谱P(θe)的幅度值作为纵坐标,在幅度谱图中以点的形式标出。按这样的方式,绘制E个目标波达方向角θ与空间谱P(θ)的幅度谱图。绘制好的幅度谱图中,总共应含有E个幅度谱点。然后对所有点的进行搜索,寻找出幅度值较大的前K个谱峰,这些谱峰的峰值点所对应的横轴坐标就是本方法MUSIC求解对非圆信号波达方向角的估计值。
Claims (6)
1.一种基于MUSIC求解的非圆信号波达方向角估计方法,其特征在于,包括如下步骤:
(1)构建嵌套阵列c:
(1a)以天线接收机作为阵元,构建包括M个阵元的第一均匀线性阵列a和包括N个阵元且相邻阵元的间距为(M+1)d的第二均匀线性阵列b,M≥1,N≥1,a中相邻阵元的间距为d,b中相邻阵元的间距为(M+1)d,0<d<λ/2,λ为入射到阵元的窄带信号波长,并定义a中各阵元的位置依次为D(1)=1,D(2)=2,…,D(M)=M,定义b中各阵元的位置依次为D(M+1)=M+1,D(M+2)=2(M+1),…,D(M+N)=N(M+1);
(1b)根据第一均匀线性阵列a和第二均匀线性阵列b,构造嵌套阵列c;
(2)获取嵌套阵列c的输出信号Y(t):
将K个非圆目标信号入射到嵌套阵列c中的每一个阵元,K≥1,并对c的输出信号进行采样,然后对所有采样信号滤波,得到嵌套阵列c的输出信号Y(t):
Y(t)=[y1(t),…,yi(t),…,yM+N(t)]T
其中,yi(t)表示c中第i个阵元的输出信号,1≤i≤M+N,t表示离散时间,1≤t≤L,L表示信号在时域上采样点数,[·]T表示矩阵[·]的转置运算;
(3)计算嵌套阵列c的虚拟均匀阵列协方差向量和椭圆协方差向量
(3a)计算等效协方差向量rd和等效椭圆协方差向量rs;
(3b)计算虚拟均匀阵列协方差向量和虚拟均匀阵列椭圆协方差向量
(4)构造扩展协方差矩阵G;
(5)计算扩展协方差矩阵G的噪声子空间Un:
对扩展协方差矩阵G特征分解,得到G的特征值矩阵Λ和G的特征值所对应的特征向量矩阵U,并对Λ中的特征值按照从大到小的顺序进行排列,然后将(L1+L2-K)个较小特征值在U中对应的特征向量作为G的噪声子空间Un;
(6)构造目标波达方向角θ的阵列流型矩阵Q(θ):
(6a)将嵌套阵列c[0°,180°]的观测方向范围按0°到180°的递增方向均匀划分成E+1个子观测方向区域,E≥M+N,并将第e个子观测方向区域和第e+1个子观测方向区域之间的边界射线与0°射线之间的夹角作为第e个目标波达方向角θe,其中e=1,2,…,E,定义E个目标波达方向角θ:
θ=[θ1,θ2,…,θe,…,θE]
(6b)构造目标波达方向角θ的阵列流型向量α1(θ)和阵列流型向量α2(θ),其中θe的阵列流型向量α1(θe)和阵列流型向量α2(θe)的表达式分别为:
其中,j为虚数单位;
(6c)构造目标波达方向角θ的阵列流型矩阵Q(θ),其中θe的阵列流型矩阵Q(θe)的表达式为:
(7)基于MUSIC求解的方法,计算目标波达方向角θ的空间谱P(θ),其中θe的空间谱P(θe)计算公式为:
其中,det[·]表示矩阵[·]的求行列式运算;
(8)绘制目标波达方向角θ与空间谱P(θ)的幅度谱图:
以目标波达方向角θ为x轴坐标,以空间谱P(θ)的幅度值为y轴坐标,绘制幅度谱图,从该幅度谱图中按照幅度值从高到低的顺序寻找幅度值较大的前K个谱峰,K个谱峰的峰值点所对应的x轴坐标即为非圆信号波达方向角的估计值。
2.根据权利要求1所述的基于MUSIC求解的非圆信号波达方向角估计方法,其特征在于,步骤(1b)中所述的构造嵌套阵列c,构造方法为:
将b的头部与a的尾部进行拼接,且b中位置为D(M+1)=M+1的阵元与a中位置为D(M)=M的阵元之间的距离为d,得到嵌套阵列c;
其中a为第一均匀线性阵列,b为第二均匀线性阵列,M为a中阵元数个阵元,M≥1,D(M+1)=M+1为b中第一个阵元的位置,D(M)=M为a中第M个阵元的位置。
3.根据权利要求1所述的基于MUSIC求解的非圆信号波达方向角估计方法,其特征在于,步骤(3a)中所述的计算等效协方差向量rd和等效椭圆协方差向量rs,计算方法如下:
根据嵌套阵列c的输出信号Y(t),计算等效协方差矩阵Rd和等效椭圆协方差矩阵Rs,计算公式分别为:
其中,Y(t)表示嵌套阵列c的输出信号,t表示离散时间,1≤t≤L,L表示信号在时域上采样点数,[·]H表示矩阵[·]的共轭转置运算,[·]T表示矩阵[·]的共轭转置运算;
将Rd中的各列元素按照由小到大的列号进行拼接,得到等效协方差向量rd,同时将Rs中的各列元素按照由小到大的列号进行拼接,得到等效椭圆协方差向量rs,rd和rs的表达式分别为:
rd=[Rd(1,1),Rd(2,1),…,Rd(M+N,1),Rd(1,2),…,Rd(M+N,2),…,Rd(1,M+N),…,Rd(M+N,M+N)]T,
rs=[Rs(1,1),Rs(2,1),…,Rs(M+N,1),Rs(1,2),…,Rs(M+N,2),…,Rs(1,M+N),…,Rs(M+N,M+N)]T,
其中,Rd(i,j)表示协方差矩阵Rd中位于第i行,第j列的元素,1≤i≤M+N,1≤j≤M+N,M+N为嵌套阵列c的阵元数,Rs(i,j)表示椭圆协方差矩阵Rs中位于第i行,第j列的元素,[·]T表示矩阵[·]的共轭转置运算。
4.根据权利要求1所述的基于MUSIC求解的非圆信号波达方向角估计方法,其特征在于,步骤(3b)中所述的计算虚拟均匀阵列协方差向量和虚拟均匀阵列椭圆协方差向量其计算方法如下:
计算等效协方差向量rd中所有元素的维数Ei,j和等效椭圆协方差向量rs中所有元素的维数Fi,j,计算公式分别为:
Ei,j=D(j)-D(i)
Fi,j=D(j)+D(i)
其中,D(i)表示嵌套阵列c中第i个阵元的位置,D(j)表示嵌套阵列c中第j个阵元的位置,1≤i≤M+N,1≤j≤M+N,M+N为嵌套阵列c的阵元数;
删除rd中维数相同的元素和维数不连续的元素后,将剩余元素按维数从小到大进行排列,得到元素个数为Cd的虚拟均匀阵列协方差向量同时删除rs中维数相同的元素和维数不连续的元素后,将剩余元素按维数从小到大进行排列,得到元素个数为Cs的虚拟均匀阵列椭圆协方差向量其中,Cd=2MN+2N-1,Cs=MN+M+N。
5.根据权利要求1所述的基于MUSIC求解的非圆信号波达方向角估计方法,其特征在于,步骤(4)中所述的构造扩展协方差矩阵G,G的表达式如下:
根据虚拟均匀阵列协方差向量和椭圆协方差向量构造扩展协方差矩阵G;
其中,L1和L2表示中间变量,L1=(Cd+1)/2,L2=Cs+1-(Cd+1)/2,且有L1>L2,(·)*表示向量的共轭运算。
6.根据权利要求1所述的基于MUSIC求解的非圆信号波达方向角估计方法,其特征在于,步骤(5)中所述的扩展协方差矩阵G特征分解,其公式为:
G=U·Λ·UH
其中,Λ为扩展协方差矩阵G的特征值矩阵,U为矩阵G的特征值所对应的特征向量矩阵,[·]H表示矩阵[·]的共轭转置运算。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910203570.3A CN109946643B (zh) | 2019-03-18 | 2019-03-18 | 基于music求解的非圆信号波达方向角估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910203570.3A CN109946643B (zh) | 2019-03-18 | 2019-03-18 | 基于music求解的非圆信号波达方向角估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109946643A true CN109946643A (zh) | 2019-06-28 |
CN109946643B CN109946643B (zh) | 2022-08-26 |
Family
ID=67008865
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910203570.3A Active CN109946643B (zh) | 2019-03-18 | 2019-03-18 | 基于music求解的非圆信号波达方向角估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109946643B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111175693A (zh) * | 2020-01-19 | 2020-05-19 | 河北科技大学 | 一种波达方向估计方法及波达方向估计装置 |
CN112363106A (zh) * | 2020-10-28 | 2021-02-12 | 西安电子科技大学 | 基于量子粒子群的信号子空间波达方向检测方法及系统 |
CN112731273A (zh) * | 2020-12-09 | 2021-04-30 | 南京邮电大学 | 一种基于稀疏贝叶斯的低复杂度信号波达方向估计方法 |
CN113391266A (zh) * | 2021-05-28 | 2021-09-14 | 南京航空航天大学 | 基于非圆多嵌套阵降维子空间数据融合的直接定位方法 |
CN114879136A (zh) * | 2022-06-10 | 2022-08-09 | 西安电子科技大学 | 基于非均匀排列的超表面阵列信号波达方向估计方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105403874A (zh) * | 2015-12-25 | 2016-03-16 | 西安电子科技大学 | 非均匀阵列欠定波达方向估计方法 |
CN106019215A (zh) * | 2016-07-15 | 2016-10-12 | 西安电子科技大学 | 基于四阶累量的嵌套阵列波达方向角估计方法 |
CN107037393A (zh) * | 2017-05-19 | 2017-08-11 | 西安电子科技大学 | 基于嵌套阵列的非圆信号波达方向角估计方法 |
JP2017150903A (ja) * | 2016-02-23 | 2017-08-31 | 日本電信電話株式会社 | 音源定位装置、方法、及びプログラム |
CN108957391A (zh) * | 2018-07-24 | 2018-12-07 | 北京理工大学 | 一种基于嵌套阵列的l型天线阵的二维波达方向估计方法 |
-
2019
- 2019-03-18 CN CN201910203570.3A patent/CN109946643B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105403874A (zh) * | 2015-12-25 | 2016-03-16 | 西安电子科技大学 | 非均匀阵列欠定波达方向估计方法 |
JP2017150903A (ja) * | 2016-02-23 | 2017-08-31 | 日本電信電話株式会社 | 音源定位装置、方法、及びプログラム |
CN106019215A (zh) * | 2016-07-15 | 2016-10-12 | 西安电子科技大学 | 基于四阶累量的嵌套阵列波达方向角估计方法 |
CN107037393A (zh) * | 2017-05-19 | 2017-08-11 | 西安电子科技大学 | 基于嵌套阵列的非圆信号波达方向角估计方法 |
CN108957391A (zh) * | 2018-07-24 | 2018-12-07 | 北京理工大学 | 一种基于嵌套阵列的l型天线阵的二维波达方向估计方法 |
Non-Patent Citations (2)
Title |
---|
沈志博等: "信源数未知情况下基于嵌套阵列的DOA估计算法", 《航天电子对抗》 * |
陈浩等: "基于内插阵列变换的非圆信号MUSIC算法", 《电视技术》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111175693A (zh) * | 2020-01-19 | 2020-05-19 | 河北科技大学 | 一种波达方向估计方法及波达方向估计装置 |
CN112363106A (zh) * | 2020-10-28 | 2021-02-12 | 西安电子科技大学 | 基于量子粒子群的信号子空间波达方向检测方法及系统 |
CN112363106B (zh) * | 2020-10-28 | 2022-11-01 | 西安电子科技大学 | 基于量子粒子群的信号子空间波达方向检测方法及系统 |
CN112731273A (zh) * | 2020-12-09 | 2021-04-30 | 南京邮电大学 | 一种基于稀疏贝叶斯的低复杂度信号波达方向估计方法 |
CN112731273B (zh) * | 2020-12-09 | 2023-06-23 | 南京邮电大学 | 一种基于稀疏贝叶斯的低复杂度信号波达方向估计方法 |
CN113391266A (zh) * | 2021-05-28 | 2021-09-14 | 南京航空航天大学 | 基于非圆多嵌套阵降维子空间数据融合的直接定位方法 |
CN114879136A (zh) * | 2022-06-10 | 2022-08-09 | 西安电子科技大学 | 基于非均匀排列的超表面阵列信号波达方向估计方法 |
CN114879136B (zh) * | 2022-06-10 | 2024-08-16 | 西安电子科技大学 | 基于非均匀排列的超表面阵列信号波达方向估计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109946643B (zh) | 2022-08-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109946643A (zh) | 基于music求解的非圆信号波达方向角估计方法 | |
CN103235292B (zh) | 平面相控阵调零保形校正的全维和差测角方法 | |
CN109298383B (zh) | 一种基于变分贝叶斯推断的互质阵波达方向角估计方法 | |
CN109450499B (zh) | 一种基于导向矢量和空间功率估计的鲁棒波束形成方法 | |
CN105335336B (zh) | 一种传感器阵列的稳健自适应波束形成方法 | |
CN103983952A (zh) | 一种非圆信号双基地mimo雷达低复杂度收发角度联合估计方法 | |
CN106021637A (zh) | 互质阵列中基于迭代稀疏重构的doa估计方法 | |
CN110398732B (zh) | 低计算量自适应步长迭代搜索的目标方向检测方法 | |
CN106855622A (zh) | 一种子阵级相控阵雷达的测角方法 | |
CN108828586B (zh) | 一种基于波束域的双基地mimo雷达测角优化方法 | |
CN111736126A (zh) | 一种基于低成本安防雷达的信号处理方法 | |
CN113189592B (zh) | 考虑幅相互耦误差的车载毫米波mimo雷达测角方法 | |
CN109471063B (zh) | 基于延迟快拍的均匀线列阵高分辨波达方向估计方法 | |
CN106226729B (zh) | 基于四阶累量的互质阵列波达方向角估计方法 | |
CN107092005A (zh) | 一种时空联合的阵列处理方法及装置 | |
CN115856767B (zh) | 一种可重构智能超表面辅助的波到达方向估计方法 | |
CN101344582B (zh) | 半盲的最小方差无失真响应波束形成方法 | |
CN110413939B (zh) | 一种基于原子范数的到达角估计方法 | |
CN110196417B (zh) | 基于发射能量集中的双基地mimo雷达角度估计方法 | |
CN106125039B (zh) | 基于局域联合处理的改进空时自适应单脉冲测角方法 | |
CN113671485B (zh) | 基于admm的米波面阵雷达二维doa估计方法 | |
CN107894581A (zh) | 一种宽带阵列波达方向估计方法 | |
CN117527043A (zh) | 一种基于副瓣抑制的卫星通信抗干扰方法及装置 | |
CN108845298B (zh) | 基于杂波映射的自适应波束形成方法 | |
CN102062851A (zh) | 基于改进l阵的星载宽带多目标的测向方法 |
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 |