CN112346004A - 一种基于score函数的广义协方差的多信号分类算法 - Google Patents
一种基于score函数的广义协方差的多信号分类算法 Download PDFInfo
- 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
Links
- 238000007635 classification algorithm Methods 0.000 title claims abstract description 23
- 239000011159 matrix material Substances 0.000 claims abstract description 25
- 230000003595 spectral effect Effects 0.000 claims abstract description 8
- 238000001228 spectrum Methods 0.000 claims abstract description 8
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 7
- 150000001875 compounds Chemical class 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000012545 processing Methods 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 239000000654 additive Substances 0.000 claims description 2
- 230000000996 additive effect Effects 0.000 claims description 2
- 239000006185 dispersion Substances 0.000 claims description 2
- 238000005457 optimization Methods 0.000 claims description 2
- 239000004576 sand Substances 0.000 claims description 2
- 239000000126 substance Substances 0.000 claims description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims 1
- 238000004422 calculation algorithm Methods 0.000 abstract description 17
- 238000004891 communication Methods 0.000 abstract description 5
- 238000000034 method Methods 0.000 abstract 1
- 238000002474 experimental method Methods 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
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
-
- 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
- G06F17/16—Matrix 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-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)所示:
式中,fα(x)表示对称Alpha稳定分布的概率密度函数,其中,该分布的中心参数为0,分散系数为1,α表示该分布的特征指数;fα′(x)表示fα(x)的导数。
具体而言,当α≠1时,fα(x)的定义式如式(20)所示:
式中,v(θ)可以表示为:
当α=1时,fα(x)=f1(x)的定义式如式(22)所示:
如图2所示,是特征指数α取值不同时多条score函数的曲线图。
S2:获得的天线阵列所输出的混合着噪声的信号,并构建该信号的向量:
首先,获得线性阵列第m个阵元输出信号xm(t),其表达式如式(23)所示:
式中,t表示连续时间变量,K表示无线电信号源的个数,θ1,θ2,…θK表示无线电信号的波达方向,sk(t)表示第k个无线电源信号,nm(t)表示第m个阵元上的加性噪声,am(θk)表示第k个波达方向θk所对应的导向系数,其表达式如式(24)所示:
am(θk)=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)=[a1(θk),a2(θk),…,aM(θk)]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)所示:
式中,T表示连续时间信号的总时长。
S4:计算并获得基于score函数的广义协方差矩阵的估计:
在实际应用中进行信号处理时,通常情况下,连续时间信号通过采样变为离散时间序列,此外,连续时间信号的总时长和离散时间序列的总长度都是有限的。
式中,l表示离散时间变量,离散时间序列xi[l]和xj[l]分别表示连续时间信号xi(t)和xj(t)的第l个采样点,L表示离散时间序列的总长度。
S6:根据噪声子空间的估计,计算并获得空间谱的估计:
首先,将无线电信号的波达方向估计问题转换为优化问题,如式(35)所示:
式中,a(θ)表示波达方向θ所对应的导向矢量。
计算并获得空间谱的估计,如式(36)所示:
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:通过谱峰位置进行无线电信号的波达方向估计。
3.如权利1所述的一种基于score函数的广义协方差的多信号分类算法,其特征在于,所述步骤S2具体包括以下步骤:
首先,获得线性阵列第m个阵元输出信号xm(t),其表达式如式(5)所示:
式中,t表示连续时间变量,K表示无线电信号源的个数,θ1,θ2,…θK表示无线电信号的波达方向,sk(t)表示第k个无线电源信号,nm(t)表示第m个阵元上的加性噪声,am(θk)表示第k个波达方向θk所对应的导向系数,其表达式如式(6)所示:
am(θk)=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)=[a1(θk),a2(θk),…,aM(θk)]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)所示
式中,T表示连续时间信号的总时长。
8.如权利1所述的一种基于score函数的广义协方差的多信号分类算法,其特征在于,所述步骤S7具体包括:首先,根据无线电信号源的个数K,选取K个空间谱的估计的谱峰;然后,根据所选取K个谱峰对无线电信号的波达方向进行估计。
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)
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)
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 | 大连理工大学 | 一种无直达信号条件下的无线电波达方向估计方法 |
-
2020
- 2020-11-20 CN CN202011307326.0A patent/CN112346004A/zh active Pending
Patent Citations (3)
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)
Title |
---|
张金凤: ""脉冲噪声环境下波达方向估计方法研究"", 《中国博士学位论文全文数据库信息科技辑, pages 14 * |
栾声扬: ""有界非线性协方差与相关熵及在无线定位中的应用"", 《中国博士学位论文全文数据库工程科技Ⅱ辑》, pages 21 - 33 * |
Cited By (4)
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 |