CN112346004A - 一种基于score函数的广义协方差的多信号分类算法 - Google Patents

一种基于score函数的广义协方差的多信号分类算法 Download PDF

Info

Publication number
CN112346004A
CN112346004A CN202011307326.0A CN202011307326A CN112346004A CN 112346004 A CN112346004 A CN 112346004A CN 202011307326 A CN202011307326 A CN 202011307326A CN 112346004 A CN112346004 A CN 112346004A
Authority
CN
China
Prior art keywords
signal
score
formula
function
estimate
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.)
Pending
Application number
CN202011307326.0A
Other languages
English (en)
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.)
Jiangsu Normal University
Original Assignee
Jiangsu Normal 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 Jiangsu Normal University filed Critical Jiangsu Normal University
Priority to CN202011307326.0A priority Critical patent/CN112346004A/zh
Publication of CN112346004A publication Critical patent/CN112346004A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Direction-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/02Direction-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/14Systems for determining direction or deviation from predetermined direction
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization

Abstract

本发明提供了一种基于score函数的广义协方差的多信号分类算法,属于通信技术领域。该方法主要包括以下几点:1.选取适当参数,计算并获得score函数。2.获得天线阵列所输出的混合着噪声的信号,并构建该信号的向量。3.计算并获得基于score函数的广义协方差矩阵。4.计算并获得基于score函数的广义协方差矩阵的估计。5.对基于score函数的广义协方差矩阵进行特征值分解,计算并获得噪声子空间的估计。6.根据噪声子空间的估计,计算并获得空间谱的估计。7.通过谱峰位置进行无线电信号的波达方向估计。实验结果表明该算法能够在脉冲性噪声条件下准确地实现无线电信号的波达方向估计。

Description

一种基于score函数的广义协方差的多信号分类算法
技术领域
本发明属于通信技术领域,涉及到无线电信号的波达方向估计,特别涉及到基于score函数的广义协方差的多信号分类(SCORE-MUSIC)算法。
背景技术
通信技术的不断发展,在日常的学习、工作和生活中,人们也越来越依赖于无线电通信,与此同时,无线电通信也面临着愈发复杂的电磁环境,此类复杂的电磁环境增加了无线电信号分析和处理的难度,而脉冲性噪声就是此类复杂电磁环境的典型代表。故而,如何提升各类信号处理方法抑制脉冲性噪声的能力就成为亟须解决的问题。
多信号分类(multiple signal classification,MUSIC)算法是无线电信号波达方向估计的典型算法,该方法通常通过对天线阵列所输出的混合着噪声的信号的协方差矩阵进行特征值分解来提取信号子空间和噪声子空间的信息,再利用子空间的正交性来估计无线电信号的波达方向,通常情况下,该方法假设噪声是服从高斯分布的,在这一假设条件下,混合着噪声的信号的协方差是有界的。但是,在某些复杂电磁环境中噪声是服从非高斯分布的,其幅度具有较强的脉冲性,其概率分布具有较重的拖尾,通常可以使用Alpha稳定分布进行表示,此时协方差是不收敛的,故而传统的多信号分类算法的性能会出现下降,甚至完全失效的现象。
为了解决该问题,本发明首先提出了一种基于柯西分布score函数的广义协方差的概念,然后提出了一种基于该广义协方差的多信号分类算法。在本发明所提出的算法中,采用score函数引入非线性映射,并借此抑制脉冲性噪声,实验结果表明该算法能够在脉冲性噪声条件下准确地实现无线电信号的波达方向估计。
发明内容
本发明提供了一种基于score函数的广义协方差的多信号分类算法。
本发明采用的技术方案是:一种基于score函数的广义协方差的多信号分类算法,其特征在于,包括以下步骤:
S1:选取适当参数,计算并获得score函数;
S2:获得天线阵列所输出的混合着噪声的信号,并构建该信号的向量;
S3:计算并获得基于score函数的广义协方差矩阵;
S4:计算并获得基于score函数的广义协方差矩阵的估计;
S5:对基于score函数的广义协方差矩阵进行特征值分解,计算并获得噪声子空间的估计;
S6:根据噪声子空间的估计,计算并获得空间谱的估计。
S7:通过谱峰位置进行无线电信号的波达方向估计。
附图说明
图1是本发明中一种基于score函数的广义协方差的多信号分类算法的总体流程图。
图2是参数取值不同时多条score函数的曲线图。
图3是广义信噪比不同时,基于score函数的广义协方差的多信号分类算法与其他算法性能的对比图。
图4是Alpha稳定分布噪声的特征指数不同时,基于score函数的广义协方差的多信号分类算法与其他算法性能的对比图。
具体实施方案
为了便于理解,下面将结合本发明实施例中的附图对本发明实施例中的技术方案进行详细描述。
如图1所示,一种基于score函数的广义协方差的多信号分类算法,主要包括以下步骤:
S1:选取适当参数,计算并获得score函数:
首先,为score函数的参数α在区间[0.6,2]内取值;
然后,计算并获得score函数,如式(19)所示:
Figure BDA0002788669000000021
式中,fα(x)表示对称Alpha稳定分布的概率密度函数,其中,该分布的中心参数为0,分散系数为1,α表示该分布的特征指数;fα′(x)表示fα(x)的导数。
具体而言,当α≠1时,fα(x)的定义式如式(20)所示:
Figure BDA0002788669000000031
式中,v(θ)可以表示为:
Figure BDA0002788669000000032
当α=1时,fα(x)=f1(x)的定义式如式(22)所示:
Figure BDA0002788669000000033
如图2所示,是特征指数α取值不同时多条score函数的曲线图。
S2:获得的天线阵列所输出的混合着噪声的信号,并构建该信号的向量:
首先,获得线性阵列第m个阵元输出信号xm(t),其表达式如式(23)所示:
Figure BDA0002788669000000034
式中,t表示连续时间变量,K表示无线电信号源的个数,θ12,…θK表示无线电信号的波达方向,sk(t)表示第k个无线电源信号,nm(t)表示第m个阵元上的加性噪声,amk)表示第k个波达方向θk所对应的导向系数,其表达式如式(24)所示:
amk)=exp(-j(2π(m-1)dsin(θk))/λ) (24)
式中,d表示阵元之间的距离,λ表示载波的波长。
然后,构建线性阵列的输出向量x(t),其表达式如式(25)表示:
x(t)=[x1(t),x2(t),…,xM(t)]T (25)
具体而言,线性阵列的输出向量x(t)满足如下关系式:
x(t)=As(t)+n(t) (26)
式中,A=[a(θ1),a(θ2),…,a(θK)]表示导向矩阵,a(θk)=[a1k),a2k),…,aMk)]T表示导向矢量,s(t)=[s1(t),s2(t),…,sK(t)]T表示无线电源信号矢量,n(t)=[n1(t),n2(t),…,nM(t)]T表示噪声矢量。
S3:计算并获得基于score函数的广义协方差矩阵:
首先,计算基于score函数的阵列的输出矢量xs(t),如式(27)所示:
xs(t)=score(x(t))=[score(x1(t)),score(x2(t)),…,score(xM(t))]T (27)
然后,计算基于矢量xs(t)的广义协方差矩阵Rs,如式(28)所示:
Figure BDA0002788669000000041
式中,E表示统计期望操作符,上标T表示转置操作符,
Figure BDA0002788669000000042
表示位于广义协方差矩阵Rs的第i行、第j列的元素,且i,j=1,2,…,M,
Figure BDA0002788669000000043
的表达式如式(29)所示:
Figure BDA0002788669000000044
如果信号xi(t)和xj(t)满足遍历性,则可以利用时间平均代替其统计平均,即利用时间平均计算
Figure BDA0002788669000000045
如式(30)所示:
Figure BDA0002788669000000046
式中,T表示连续时间信号的总时长。
S4:计算并获得基于score函数的广义协方差矩阵的估计:
在实际应用中进行信号处理时,通常情况下,连续时间信号通过采样变为离散时间序列,此外,连续时间信号的总时长和离散时间序列的总长度都是有限的。
首先,计算广义协方差矩阵的估计
Figure BDA0002788669000000047
的第i行、第j列元素
Figure BDA0002788669000000048
其中,
Figure BDA0002788669000000049
的表达式如式(31)所示:
Figure BDA0002788669000000051
式中,l表示离散时间变量,离散时间序列xi[l]和xj[l]分别表示连续时间信号xi(t)和xj(t)的第l个采样点,L表示离散时间序列的总长度。
然后,计算并获得广义协方差矩阵的估计
Figure BDA0002788669000000052
其表达式如式(32)所示:
Figure BDA0002788669000000053
S5:对广义协方差矩阵的估计
Figure BDA0002788669000000054
进行特征值分解,计算并获得噪声子空间的估计
Figure BDA0002788669000000055
首先,对广义协方差矩阵的估计
Figure BDA0002788669000000056
进行特征值分解,如式(33)所示:
Figure BDA0002788669000000057
式中,
Figure BDA0002788669000000058
Figure BDA0002788669000000059
分别表示信号子空间的估计和噪声子空间的估计,
Figure BDA00027886690000000510
Figure BDA00027886690000000511
分别表示信号和噪声所对应特征值组成的对角矩阵的估计,上标H表示共轭转置操作符。
计算并获得噪声子空间的估计
Figure BDA00027886690000000512
如式(34)所示:
Figure BDA00027886690000000513
式中,
Figure BDA00027886690000000514
分别表示与噪声相关的特征值所对应的M-K个特征矢量。
S6:根据噪声子空间的估计,计算并获得空间谱的估计:
首先,将无线电信号的波达方向估计问题转换为优化问题,如式(35)所示:
Figure BDA00027886690000000515
式中,a(θ)表示波达方向θ所对应的导向矢量。
计算并获得空间谱的估计,如式(36)所示:
Figure BDA0002788669000000061
S7:通过谱峰位置进行无线电信号的波达方向估计:
首先,根据无线电信号源的个数K,选取K个空间谱的估计的谱峰;
然后,根据所选取K个谱峰对无线电信号的波达方向进行估计。
采用仿真实验来验证算法性能。在每组实验条件下,进行200次蒙特卡罗实验;性能通过波达方向真实值和估计值的均方根误差来衡量;最终的实验结果是当前实验条件下所有均方根误差的平均值。
参与实验的算法包括:SCORE-MUSIC(本发明的算法),CRCO-MUSIC,FLOM-MUSIC,PFLOM-MUSIC,ROC-MUSIC和GZMNL-MUSIC。
实验的条件为:K=2个无线电信号独立且具有相同功率,调制方式是QAM,他们的波达方向分别为θ1=5°和θ2=15°,上述无线电信号被含有M=5个阵元的均匀线性阵列所接收,相邻阵元的间距d为无线电信号波长λ的1/2,经过采样后,连续时间信号所对应的离散时间序列的总长度为L=1000。
广义信噪比的取值范围为[-10,15]dB,Alpha稳定分布噪声的特征指数为1.5。基于score函数的广义协方差的多信号分类算法与其他算法性能的对比图如图3所示。
由图3可见,当广义信噪比改变时,基于score函数的广义协方差的多信号分类算法始终具有最小的均方根误差,说明相比于其他算法,本发明算法具有最好的性能。
广义信噪比为5dB,Alpha稳定分布噪声的特征指数的取值范围为[0.5,2]。基于score函数的广义协方差的多信号分类算法与其他算法性能的对比图如图4所示。
由图4可见,Alpha稳定分布噪声的特征指数改变时,基于score函数的广义协方差的多信号分类算法始终具有最小的均方根误差,说明相比于其他算法,本发明算法具有很好的性能。

Claims (8)

1.一种基于score函数的广义协方差的多信号分类算法,其特征在于,包括以下步骤:
S1:选取适当参数,计算并获得score函数;
S2:获得天线阵列所输出的混合着噪声的信号,并构建该信号的向量;
S3:计算并获得基于score函数的广义协方差矩阵;
S4:计算并获得基于score函数的广义协方差矩阵的估计;
S5:对基于score函数的广义协方差矩阵进行特征值分解,计算并获得噪声子空间的估计;
S6:根据噪声子空间的估计,计算并获得空间谱的估计。
S7:通过谱峰位置进行无线电信号的波达方向估计。
2.如权利1所述的一种基于score函数的广义协方差的多信号分类算法,其特征在于,所述步骤S1具体包括以下步骤:
首先,为score函数的参数α在区间[0.6,2]内取值;然后,计算并获得score函数,如式(1)所示:
Figure FDA0002788668990000011
式中,fα(x)表示对称Alpha稳定分布的概率密度函数,其中,该分布的中心参数为0,分散系数为1,α表示该分布的特征指数;f′α(x)表示fα(x)的导数。
具体而言,当α≠1时,fα(x)的定义式如(2)所示:
Figure FDA0002788668990000012
式中,v(θ)可以表示为:
Figure FDA0002788668990000013
当α=1时,fα(x)=f1(x)的定义式如(4)所示:
Figure FDA0002788668990000014
3.如权利1所述的一种基于score函数的广义协方差的多信号分类算法,其特征在于,所述步骤S2具体包括以下步骤:
首先,获得线性阵列第m个阵元输出信号xm(t),其表达式如式(5)所示:
Figure FDA0002788668990000021
式中,t表示连续时间变量,K表示无线电信号源的个数,θ12,…θK表示无线电信号的波达方向,sk(t)表示第k个无线电源信号,nm(t)表示第m个阵元上的加性噪声,amk)表示第k个波达方向θk所对应的导向系数,其表达式如式(6)所示:
amk)=exp(-j(2π(m-1)dsin(θk))/λ) (6)
式中,d表示阵元之间的距离,λ表示载波的波长。
然后,构建线性阵列的输出向量x(t),其表达式如式(7)所示:
x(t)=[x1(t),x2(t),…,xM(t)]T (7)
具体而言,线性阵列的输出向量x(t)满足如下关系式:
x(t)=As(t)+n(t) (8)
式中,
A=[a(θ1),a(θ2),…,a(θK)]表示导向矩阵,a(θk)=[a1k),a2k),…,aMk)]T表示导向矢量,s(t)=[s1(t),s2(t),…,sK(t)]T表示无线电源信号矢量,n(t)=[n1(t),n2(t),…,nM(t)]T表示噪声矢量。
4.如权利1所述的一种基于score函数的广义协方差的多信号分类算法,其特征在于,所述步骤S3具体包括以下步骤:
首先,计算基于score函数的阵列的输出矢量xs(t),如式(9)所示:
xs(t)=score(x(t))=[score(x1(t)),score(x2(t)),…,score(xM(t))]T (9)
然后,计算基于矢量xs(t)的广义协方差矩阵Rs,如式(10)所示
Figure FDA0002788668990000031
式中,E表示统计期望操作符,上标T表示转置操作符,
Figure FDA0002788668990000032
表示位于广义协方差矩阵Rs的第i行、第j列的元素,且i,j=1,2,…,M,
Figure FDA0002788668990000033
的表达式如式(11)所示:
Figure FDA0002788668990000034
如果信号xi(t)和xj(t)满足遍历性,则可以利用时间平均代替其统计平均,即利用时间平均计算
Figure FDA0002788668990000035
如式(12)所示:
Figure FDA0002788668990000036
式中,T表示连续时间信号的总时长。
5.如权利1所述的一种基于score函数的广义协方差的多信号分类算法,其特征在于,所述步骤S4具体包括以下步骤:
在实际应用中进行信号处理时,通常情况下,连续时间信号通过采样变为离散时间序列,此外,连续时间信号的总时长和离散时间序列的总长度都是有限的。
首先,计算广义协方差矩阵的估计
Figure FDA0002788668990000037
的第i行、第j列元素
Figure FDA0002788668990000038
其中,
Figure FDA0002788668990000039
的表达式如式(13)所示:
Figure FDA00027886689900000310
式中,l表示离散时间变量,离散时间序列xi[l]和xj[l]分别表示连续时间信号xi(t)和xj(t)的第l个采样点,L表示离散时间序列的总长度。
然后,计算并获得广义协方差矩阵的估计
Figure FDA00027886689900000311
其表达式如式(14)所示:
Figure FDA00027886689900000312
6.如权利1所述的一种基于score函数的广义协方差的多信号分类算法,其特征在于,所述步骤S5具体包括以下步骤:
首先,对广义协方差矩阵的估计
Figure FDA0002788668990000041
进行特征值分解,如式(15)所示:
Figure FDA0002788668990000042
式中,
Figure FDA0002788668990000043
Figure FDA0002788668990000044
分别表示信号子空间的估计和噪声子空间的估计,
Figure FDA0002788668990000045
Figure FDA0002788668990000046
分别表示信号和噪声所对应特征值组成的对角矩阵的估计,上标H表示共轭转置操作符。计算并获得噪声子空间的估计
Figure FDA0002788668990000047
如式(16)所示:
Figure FDA0002788668990000048
式中,
Figure FDA0002788668990000049
分别表示与噪声相关的特征值所对应的M-K个特征矢量。
7.如权利1所述的一种基于score函数的广义协方差的多信号分类算法,其特征在于,所述步骤S6具体包括以下步骤:
首先,将无线电信号的波达方向估计问题转换为优化问题,如式(17)所示:
Figure FDA00027886689900000410
式中,a(θ)表示波达方向θ所对应的导向矢量。计算并获得空间谱的估计,如式(18)所示:
Figure FDA00027886689900000411
8.如权利1所述的一种基于score函数的广义协方差的多信号分类算法,其特征在于,所述步骤S7具体包括:首先,根据无线电信号源的个数K,选取K个空间谱的估计的谱峰;然后,根据所选取K个谱峰对无线电信号的波达方向进行估计。
CN202011307326.0A 2020-11-20 2020-11-20 一种基于score函数的广义协方差的多信号分类算法 Pending CN112346004A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011307326.0A CN112346004A (zh) 2020-11-20 2020-11-20 一种基于score函数的广义协方差的多信号分类算法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011307326.0A CN112346004A (zh) 2020-11-20 2020-11-20 一种基于score函数的广义协方差的多信号分类算法

Publications (1)

Publication Number Publication Date
CN112346004A true CN112346004A (zh) 2021-02-09

Family

ID=74364303

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011307326.0A Pending CN112346004A (zh) 2020-11-20 2020-11-20 一种基于score函数的广义协方差的多信号分类算法

Country Status (1)

Country Link
CN (1) CN112346004A (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113052081A (zh) * 2021-03-26 2021-06-29 江苏师范大学 一种基于柯西Score星座图的调制方式识别方法
CN113205000A (zh) * 2021-04-07 2021-08-03 江苏师范大学 一种基于柯西Score极坐标图调制方式识别方法
CN114884781A (zh) * 2022-03-30 2022-08-09 西安电子科技大学 一种基于线性规划和多信号分类算法的载波频偏估计方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103901395A (zh) * 2014-03-28 2014-07-02 哈尔滨工程大学 一种冲击噪声环境下相干信号波达方向动态跟踪方法
CN105913044A (zh) * 2016-05-04 2016-08-31 大连理工大学 一种基于Sigmoid协方差矩阵的多重信号分类方法
CN110320490A (zh) * 2019-06-05 2019-10-11 大连理工大学 一种无直达信号条件下的无线电波达方向估计方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103901395A (zh) * 2014-03-28 2014-07-02 哈尔滨工程大学 一种冲击噪声环境下相干信号波达方向动态跟踪方法
CN105913044A (zh) * 2016-05-04 2016-08-31 大连理工大学 一种基于Sigmoid协方差矩阵的多重信号分类方法
CN110320490A (zh) * 2019-06-05 2019-10-11 大连理工大学 一种无直达信号条件下的无线电波达方向估计方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
张金凤: ""脉冲噪声环境下波达方向估计方法研究"", 《中国博士学位论文全文数据库信息科技辑, pages 14 *
栾声扬: ""有界非线性协方差与相关熵及在无线定位中的应用"", 《中国博士学位论文全文数据库工程科技Ⅱ辑》, pages 21 - 33 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113052081A (zh) * 2021-03-26 2021-06-29 江苏师范大学 一种基于柯西Score星座图的调制方式识别方法
CN113052081B (zh) * 2021-03-26 2024-03-15 江苏师范大学 一种基于柯西Score星座图的调制方式识别方法
CN113205000A (zh) * 2021-04-07 2021-08-03 江苏师范大学 一种基于柯西Score极坐标图调制方式识别方法
CN114884781A (zh) * 2022-03-30 2022-08-09 西安电子科技大学 一种基于线性规划和多信号分类算法的载波频偏估计方法

Similar Documents

Publication Publication Date Title
CN112346004A (zh) 一种基于score函数的广义协方差的多信号分类算法
CN110927659B (zh) 互耦条件下任意阵列流形doa估计与互耦校准方法及系统
CN110045321B (zh) 基于稀疏和低秩恢复的稳健doa估计方法
Jung et al. PCA consistency in high dimension, low sample size context
CN110208735A (zh) 一种基于稀疏贝叶斯学习的相干信号doa估计方法
US11567161B2 (en) Method for estimating the direction-of-arrival of a coprime array based on virtual domain statistics reconstruction of single-bit quantized signal
CN106501770B (zh) 基于幅相误差阵列的远近场宽带混合源中近场源定位方法
WO2022165872A1 (zh) 一种针对毫米波3d mimo信道的路径参数提取方法
CN109639333B (zh) 一种基于有效重构协方差矩阵的波束形成方法
CN114726385B (zh) 基于功率估计的卫星导航接收机空域抗干扰方法
CN110212951B (zh) 一种基于Butler矩阵的大规模MIMO信道估计方法
Molaei et al. Two‐dimensional DOA estimation for multi‐path environments by accurate separation of signals using k‐medoids clustering
Karnik et al. Thomson’s multitaper method revisited
CN114884841A (zh) 基于高阶统计和非均匀阵列的欠定参数联合估计方法
CN112668155B (zh) 一种基于二次重构的稳健波束形成方法及系统
CN114460535A (zh) 一种基于加权欧拉esprit的多阵列非圆源直接定位方法
CN114879133A (zh) 一种多径和高斯色噪声环境下的稀疏角度估计方法
Garcia-Ariza et al. Computing the closest positive definite correlation matrix for experimental MIMO channel analysis
Wei et al. Exploring conventional enhancement and separation methods for multi‐speech enhancement in indoor environments
López et al. Estimation of information in parallel Gaussian channels via model order selection
Wang et al. DOA estimation of smart antenna signal based on MUSIC algorithm
CN113805139A (zh) 基于聚焦变换的宽带信号稀疏表示波达方向估计方法
Kumar et al. Performance analysis of digital beamforming with modified ESPRIT-MUSIC direction of arrival estimation algorithm for multi input multi output non orthogonal multiple access system over Rayleigh fading channel
CN113740797B (zh) 一种透镜阵列下的高精度单快拍目标到达角度估计方法
Jiang et al. An efficient null-broadening adaptive beamformer based on magnitude response constraints

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