CN112363106B - 基于量子粒子群的信号子空间波达方向检测方法及系统 - Google Patents
基于量子粒子群的信号子空间波达方向检测方法及系统 Download PDFInfo
- Publication number
- CN112363106B CN112363106B CN202011171199.6A CN202011171199A CN112363106B CN 112363106 B CN112363106 B CN 112363106B CN 202011171199 A CN202011171199 A CN 202011171199A CN 112363106 B CN112363106 B CN 112363106B
- Authority
- CN
- China
- Prior art keywords
- signal
- subspace
- target space
- target
- matrix
- 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
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/04—Details
- G01S3/043—Receivers
-
- 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
- G01S3/46—Systems for determining direction or deviation from predetermined direction using antennas spaced apart and measuring phase or time difference between signals therefrom, i.e. path-difference systems
- G01S3/48—Systems for determining direction or deviation from predetermined direction using antennas spaced apart and measuring phase or time difference between signals therefrom, i.e. path-difference systems the waves arriving at the antennas being continuous or intermittent and the phase difference of signals derived therefrom being measured
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
本发明公开了基于量子粒子群的信号子空间波达方向检测方法和系统,包括:部署信号接收传感器阵列;按预设采样时刻通过所述信号接收传感器阵列对目标空间的信号源进行信号数据接收,接收到的信号形成信号数据矩阵;对信号数据矩阵进行特征值分解,得到噪声子空间;在目标空间中随机选取多个方向角,建立一个随机信号子空间;基于随机信号子空间和噪声子空间之间的正交特性,构造构建代价函数;采用量子粒子群算法对进行求解,得到目标信号空间的各信号源的波达方向角。本发明方法有效克服了利用传统方案进行多目标波达方向检测时,在低信噪比、小样本环境中,容易出现漏检或性能下降问题。
Description
技术领域
本发明属于信号处理技术领域,涉及一种基于量子粒子群的信号子空间波达方向检测方法和检测系统。
背景技术
波达方向(Direction Of Arrival,DOA)估计是信号处理领域的一个重要研究分支,被广泛应用于雷达、声纳、通信以及地震探测等领域中。DOA估计的主要目的是在噪声环境下,分辨两个在方位向上非常接近的目标信号,DOA估计中,对目标信号的估计精度是一项重要技术指标,随着应用场景需求的不断加强,对这项指标的要求越来越高。
在现有的技术中,常用的高分辨子空间算法,如MUSIC算法。MUSIC进行DOA估计的手段是:通过分离噪声子空间,构造基于噪声子空间的空间谱函数,进行谱峰搜索,获取DOA信息,理论上通过减小算法搜索步长可以达到DOA估计的超高精度。然而,在实际应用环境中,由于来自各方面的噪声干扰,信噪比降低,样本数量受限,MUSIC算法的运用并不能达到理想效果,经常会出现目标信息丢失、漏检的情况,而且目标检测的精度不理想。
发明内容
为解决现有技术中存在的不足,本发明提供了一种基于量子PSO的搜索信号子空间波达方向检测方法和检测系统,解决传统方案进行多目标检测时,在低信噪比、小样本环境中,容易出现漏检或检测精度下降问题。
为了解决上述技术问题,本发明采用如下技术方案予以实现:
基于量子粒子群的信号子空间波达方向检测方法,包括:
步骤1,针对待检测的目标空间部署信号接收传感器阵列;
步骤2,按预设采样时刻通过所述信号接收传感器阵列对目标空间的信号源进行信号数据接收,接收到的所有信号形成信号数据矩阵;
步骤3,对信号数据矩阵进行特征值分解,得到噪声子空间UN;
式中,为随机方向角和阵列参数构成的阵列流型矩阵;表示以传感器阵列中第一个传感器为参考,随机方向角为的信号源辐射到其他各传感器相对参考传感器的波程差形成的相位差组成的向量;表示随机选取的P个方向角组成的向量;H表示对矩阵进行共轭转置;M表示传感器阵列中传感器的数量,M≥P+1;d表示相邻传感器之间的间距;[]T表示转置操作;
式中,|| ||2表示2范数;
具体的,所述噪声子空间UN的计算过程为:对阵列接收数据矩阵的协方差矩阵R进行特征值分解,
R=ARsAH+RN (6)
A=[a(θ1),a(θ2),…,a(θp),…,a(θP)] (2)
式中,H表示对矩阵进行共轭转置;A表示目标信号源对应的阵列流型矩阵;a(θp)表示以第一个传感器为参考,第p个信号源辐射到各个传感器相对参考传感器的波程差形成的相位差组成的向量;θp表示第p个信号源的方向角;λ表示信号源的波长;j表示复数单位;Rs表示信号源的自相关矩阵,RN噪声的自相关矩阵。
得到多个特征值,对这些特征值排序:η1≥…≥ηP>ηP+1=…=ηM,其中ηP+1、ηP+2、ηP+3、…、ηM对应的特征向量张成噪声子空间UN。
本发明还公开了基于量子粒子群的信号子空间波达方向检测系统,该系统包括目标空间信号接收模块和数据处理模块;
所述目标空间信号接收模块用于对目标空间的信号源进行信号数据接收;所述目标空间信号接收模块为传感器阵列,且传感器阵列中每两个相邻传感器之间的间距相同;
所述数据处理模块用于对接收到的信号数据矩阵进行处理,具体包括:
噪声子空间估计模块,用于对信号数据矩阵进行特征值分解以获得噪声子空间;
随机信号子空间构建模块,用于在目标空间中随机选取P个方向角建立一个随机信号子空间,P为目标空间中信号源的数量;
代价函数构建模块,用于根据随机信号子空间和噪声子空间之间的正交特性来构建代价函数;
波达方向求解模块,用于采用量子粒子群算法对构建代价函数进行求解,得到目标信号空间的各信号源的波达方向角。
更具体的,所述噪声子空间估计模块包括:
特征值分解模块,用于对信号数据矩阵的协方差矩阵进行特征值分解;
特征值选择模块,用于将数值较小的M个特征值对应的特征向量张成噪声子空间。
与现有技术相比,本发明的有益效果是:
本发明通过构件代价函数,代价函数的任何一个峰值都包含着多目标DOA的全部信息,因此,相当于一个单目标优化,因此使得量子PSO的迭代次数显著降低了。这样,把多目标波达方向角估计的多目标优化问题转化为单目标优化,可以快速搜索到所有目标信号的波达方向角。有效克服了利用传统方案进行多目标波达方向检测时,在低信噪比、小样本环境中,容易出现漏检或性能下降问题。
附图说明
图1是本发明方法中信号接收传感器阵列的模型图。
图2是本发明方法与高分辨子空间MUSIC算法的检测结果。
具体实施方式
以下给出本发明的具体实施方案,需要说明的是,本发明并不局限于以下具体实施方案,凡在本申请技术方案基础上做的等同变换均落入本发明的保护范围。
本发明的具体实施例中,假设目标空间有P个信号源,P个信号源的方向角信息为Θ=[θ1,θ2,…,θp,…θP],其中,θp表示第p个信号源的方向角,P≥1。本发明的目的就是计算这些信号源的角度信息。
步骤1,部署信号接收传感器阵列:
针对上述目标空间,布置信号接收传感器阵列。本发明的传感器阵列由M个阵元(即传感器)组成,如图1所示。M≥P+1,相邻传感器之间的间距为d,且该间距不大于信号源的半波长。
步骤2,按预设采样时刻通过信号接收传感器阵列对目标空间进行信号接收,接收到的信号形成信号数据矩阵,其中t时刻的阵列接收信号数据组成的向量被表示为:
x(t)=[x1(t),x2(t),…,xM(t)]T=As(t)+n(t) (1)
A=[a(θ1),a(θ2),…,a(θp),…,a(θP)] (2)
s(t)=[s1(t),s2(t),…,sp(t),…,sP(t)]T (4)
n(t)=[n1(t),n2(t),…,nm(t),…,nM(t)]T (5)
式中,[]T表示转置操作;A表示目标信号源对应的阵列流型矩阵;a(θp)表示以第一个传感器为参考,第p个信号源到达其他各个传感器相对参考传感器的波程差形成的相位差组成的向量;λ表示信号源的波长;j表示复数单位;s(t)表示t时刻信号源信号值形成的向量;sp(t)表示t时刻第p个信号源的信号值;n(t)表示t时刻传感器接收的噪声值形成的向量;nm(t)第t时刻第m个传感器接收的噪声值。
由T个时刻阵列接收信号数据向量组成矩阵X,X=[x(1),x(2),…,x(t),…,x(T)]。
步骤三,估计噪声子空间:
对阵列接收数据矩阵X的协方差矩阵R进行特征值分解,确定信号子空间US和噪声子空间UN,
R=E[XXH]=ARsAH+RN (6)
式中,H表示对矩阵进行共轭转置,Rs表示信号源的自相关矩阵,RN噪声的自相关矩阵。
对R进行特征值分解后,得到多个特征值,这些特征值存在如下关系:
η1≥…≥ηP>ηP+1=…=ηM
在上述特征值中,由η1、η2、…、ηP对应的特征向量张成的子空间为信号子空间US,由ηP+1、ηP+2、ηP+3、…、ηM对应的特征向量张成的子空间为噪声子空间UN。
式中,表示随机选取的方向角对应的阵列流型矩阵;表示以传感器阵列中第一个传感器为参考,随机方向角为的信号源辐射到其他各传感器相对参考传感器的波程差形成的相位差组成的向量;表示随机选取的P个方向角组成的向量;H表示对矩阵进行共轭转置,CM×P表示维度为M×P的复矩阵。
式中,|| ||2表示2范数;
理论上,上述代价函数的每一个峰值都包含目标的DOA信息。因此,相当于一个单目标优化,这样,把多目标波达方向角估计的多目标优化问题转化为单目标优化,可以快速搜索到所有目标信号的波达方向角。
步骤六,代价函数的优化求解:
该代价函数是一个多变量函数,传统的数学方法难以求解。在对上述代价函数优化求解过程中,本发明选择基于量子行为的粒子群(PSO)算法算法,快速获取目标的DOA信息。因为,量子空间的粒子会被吸引至以局部吸引子为核心的量子势阱中,量子PSO算法涉及到的参数少、兼容性更强、收敛速度更快。并且在量子空间,粒子运动不会受轨道束缚,随机性程度被提高,更适合全局搜索,使得算法具有很好的收敛性。
另外,本发明还公开了一种基于量子粒子群的信号子空间波达方向检测系统,该系统主要包括目标空间信号接收模块和数据处理模块。
其中,目标空间信号接收模块用于对目标空间的信号源进行信号数据接收;具体的,目标空间信号接收模块为如图1所示的传感器阵列,且传感器阵列中每两个相邻传感器之间的间距相同;
数据处理模块用于对接收到的信号数据矩阵进行处理,该模块一般设置在计算机上,数据处理模块具体包括:
噪声子空间估计模块,用于对信号数据矩阵进行特征值分解以获得噪声子空间;具体包括特征值分解模块和特征值选择模块,其中特征值分解模块,用于对信号数据矩阵的协方差矩阵进行特征值分解;特征值选择模块,用于将数值较小的M个特征值对应的特征向量张成噪声子空间。
随机信号子空间构建模块,用于在目标空间中随机选取P个方向角建立一个随机信号子空间,P为目标空间中信号源的数量。本发明具体实施例中构建的随机信号子空间模型见公式(7)所示。
代价函数构建模块,用于根据随机信号子空间和噪声子空间之间的正交特性来构建代价函数。本发明实施例中构建的代价函数模型具体见公式(10)所示。
波达方向求解模块,用于法对构建代价函数进行求解,得到目标信号空间的各信号源的波达方向角。本发明优选的采用量子粒子群算对代价函数进行求解。
以下给出本发明的具体实验,以验证本发明方法的有效性:
实验部分
为验证本发明方法的有效性,将本发明方法与常用的高分辨子空间MUSIC算法进行对比。实验中,传感器阵列由10个(即M=10)阵元组成的均匀线阵,目标信号源为到达方向角为15°、30°的2个远场窄带信号(即P=2,θ1=15°,θ2=30°),快拍次数为30(即T=30),信噪比为-25dB~-10dB;MUSIC算法搜索步长取0.2度,噪声为零均值白复高斯噪声。采用均方根误差(Rmse)作为评价指标进行方案评估。
图1所示为50次独立实验平均结果。从实验结果可以看出,在低信噪比、小样本情况下,传统的高分辨子空间MUSIC方法在对多个目标进行检测时,检测结果偏差较大,而本发明方案可有效降低偏差。
Claims (4)
1.基于量子粒子群的信号子空间波达方向检测方法,其特征在于,包括:
步骤1,针对待检测的目标空间部署信号接收传感器阵列;传感器阵列中每两个相邻传感器之间的间距相同;
步骤2,按预设采样时刻通过所述信号接收传感器阵列对目标空间的信号源进行信号数据接收,接收到的所有信号形成信号数据矩阵;
步骤3,对信号数据矩阵进行特征值分解,得到噪声子空间UN;
式中,为随机方向角和阵列参数构成的阵列流型矩阵;表示以传感器阵列中第一个传感器为参考,随机方向角为的信号源辐射到其他各传感器相对参考传感器的波程差形成的相位差组成的向量;表示随机选取的P个方向角组成的向量;H表示对矩阵进行共轭转置;M表示传感器阵列中传感器的数量,M≥P+1;d表示相邻传感器之间的间距;[]T表示转置操作;λ表示信号源的波长;
式中,||||2表示2范数;
2.如权利要求1所述的基于量子粒子群的信号子空间波达方向检测方法,其特征在于,所述噪声子空间UN的计算过程为:对阵列接收数据矩阵的协方差矩阵R进行特征值分解,
R=ARsAH+RN (6)
A=[a(θ1),a(θ2),…,a(θp),…,a(θP)] (2)
式中,H表示对矩阵进行共轭转置;A表示目标信号源对应的阵列流型矩阵;a(θp)表示以第一个传感器为参考,第p个信号源辐射到其他各个传感器相对参考传感器的波程差形成的相位差组成的向量;θp表示第p个信号源的方向角;λ表示信号源的波长;j表示复数单位;Rs表示信号源的自相关矩阵,RN噪声的自相关矩阵;
得到多个特征值,对这些特征值排序:η1≥…≥ηP>ηP+1=…=ηM,其中ηP+1、ηP+2、ηP+3、…、ηM对应的特征向量张成噪声子空间UN。
3.基于量子粒子群的信号子空间波达方向检测系统,其特征在于,包括目标空间信号接收模块和数据处理模块;
所述目标空间信号接收模块用于对目标空间的信号源进行信号数据接收;所述目标空间信号接收模块为传感器阵列,且传感器阵列中每两个相邻传感器之间的间距相同;
所述数据处理模块用于对接收到的信号数据矩阵进行处理,具体包括:
噪声子空间估计模块,用于对信号数据矩阵进行特征值分解以获得噪声子空间;
随机信号子空间构建模块,用于在目标空间中随机选取P个方向角建立一个随机信号子空间,P为目标空间中信号源的数量;
代价函数构建模块,用于根据随机信号子空间和噪声子空间之间的正交特性来构建代价函数;
波达方向求解模块,用于采用量子粒子群算法对构建代价函数进行求解,得到目标信号空间的各信号源的波达方向角。
4.如权利要求3所述的基于量子粒子群的信号子空间波达方向检测系统,其特征在于,所述噪声子空间估计模块包括:
特征值分解模块,用于对信号数据矩阵的协方差矩阵进行特征值分解;
特征值选择模块,用于将数值较小的M个特征值对应的特征向量张成噪声子空间。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011171199.6A CN112363106B (zh) | 2020-10-28 | 2020-10-28 | 基于量子粒子群的信号子空间波达方向检测方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011171199.6A CN112363106B (zh) | 2020-10-28 | 2020-10-28 | 基于量子粒子群的信号子空间波达方向检测方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112363106A CN112363106A (zh) | 2021-02-12 |
CN112363106B true CN112363106B (zh) | 2022-11-01 |
Family
ID=74511115
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011171199.6A Active CN112363106B (zh) | 2020-10-28 | 2020-10-28 | 基于量子粒子群的信号子空间波达方向检测方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112363106B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104616059A (zh) * | 2015-01-29 | 2015-05-13 | 无锡职业技术学院 | 一种基于量子粒子群的波达方向估计方法 |
CN106093871A (zh) * | 2016-06-01 | 2016-11-09 | 河海大学 | 基于经验模态分解的智能天线波达方向估计系统及方法 |
CN109521393A (zh) * | 2018-11-05 | 2019-03-26 | 昆明理工大学 | 一种基于信号子空间旋转特性的波达方向估计算法 |
CN109946643A (zh) * | 2019-03-18 | 2019-06-28 | 西安电子科技大学 | 基于music求解的非圆信号波达方向角估计方法 |
CN111521968A (zh) * | 2020-05-22 | 2020-08-11 | 南京理工大学 | 基于目标空间分集的欠定doa估计方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11340345B2 (en) * | 2015-07-17 | 2022-05-24 | Origin Wireless, Inc. | Method, apparatus, and system for wireless object tracking |
-
2020
- 2020-10-28 CN CN202011171199.6A patent/CN112363106B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104616059A (zh) * | 2015-01-29 | 2015-05-13 | 无锡职业技术学院 | 一种基于量子粒子群的波达方向估计方法 |
CN106093871A (zh) * | 2016-06-01 | 2016-11-09 | 河海大学 | 基于经验模态分解的智能天线波达方向估计系统及方法 |
CN109521393A (zh) * | 2018-11-05 | 2019-03-26 | 昆明理工大学 | 一种基于信号子空间旋转特性的波达方向估计算法 |
CN109946643A (zh) * | 2019-03-18 | 2019-06-28 | 西安电子科技大学 | 基于music求解的非圆信号波达方向角估计方法 |
CN111521968A (zh) * | 2020-05-22 | 2020-08-11 | 南京理工大学 | 基于目标空间分集的欠定doa估计方法 |
Non-Patent Citations (3)
Title |
---|
一种基于无线传感器阵列的到达方向检测方法;聂卫科等;《传感器与微系统》;20141231;第33卷(第12期);第124-126、129页 * |
基于PSO-MUSIC算法的分布式目标DOA估计;李海等;《雷达科学与技术》;20190430;第17卷(第2期);第119-123、130页 * |
基于子空间拟合波达方向估计方法研究;苏淑靖等;《计算机仿真》;20090115(第01期);第312-316页 * |
Also Published As
Publication number | Publication date |
---|---|
CN112363106A (zh) | 2021-02-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Oispuu et al. | Direct detection and position determination of multiple sources with intermittent emission | |
CN110031794B (zh) | 一种基于差分共性阵重构的相干信源doa估计方法 | |
Shi et al. | On parameter identifiability of diversity-smoothing-based MIMO radar | |
CN111123192B (zh) | 一种基于圆形阵列和虚拟扩展的二维doa定位方法 | |
CN109254272B (zh) | 一种共点式极化mimo雷达的两维角度估计方法 | |
CN111239678A (zh) | 一种基于l型阵列的二维doa估计方法 | |
CN107290732B (zh) | 一种量子大爆炸的单基地mimo雷达测向方法 | |
CN107037398B (zh) | 一种二维music算法估计波达方向的并行计算方法 | |
CN108828504B (zh) | 基于部分相关波形的mimo雷达目标方向快速估计方法 | |
Qi et al. | Time-frequency DOA estimation of chirp signals based on multi-subarray | |
Elbir et al. | Deep-sparse array cognitive radar | |
CN111693947A (zh) | 基于互质阵列doa估计的改进music方法 | |
CN110196417B (zh) | 基于发射能量集中的双基地mimo雷达角度估计方法 | |
CN112363108B (zh) | 信号子空间加权超分辨的波达方向检测方法及系统 | |
CN112881975B (zh) | 基于子阵特征矩阵联合对角化的单脉冲和差波束测角方法 | |
CN112363106B (zh) | 基于量子粒子群的信号子空间波达方向检测方法及系统 | |
CN108594165B (zh) | 一种基于期望最大化算法的窄带信号波达方向估计方法 | |
CN114265004B (zh) | 一种基于子空间对消的干扰下的目标角度估计方法 | |
CN112415469B (zh) | 一种两维数字阵列雷达快速干扰测向方法 | |
CN113109760B (zh) | 一种基于组稀疏的多线谱联合doa估计和聚类方法及系统 | |
Ni et al. | Information-theoretic target localization with compressed measurement using FDA radar | |
CN114325560A (zh) | 波束扫描雷达超分辨目标测向方法 | |
CN113341371A (zh) | 一种基于l阵和二维esprit算法的doa估计方法 | |
CN113093097B (zh) | 一种使用互质阵列的概率假设密度doa跟踪的方法 | |
CN113111304B (zh) | 强冲击噪声下基于量子射线机理的相干分布源测向方法 |
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 |