CN108520195B - 一种基于引力搜索算法的music谱峰搜索方法 - Google Patents

一种基于引力搜索算法的music谱峰搜索方法 Download PDF

Info

Publication number
CN108520195B
CN108520195B CN201810094288.1A CN201810094288A CN108520195B CN 108520195 B CN108520195 B CN 108520195B CN 201810094288 A CN201810094288 A CN 201810094288A CN 108520195 B CN108520195 B CN 108520195B
Authority
CN
China
Prior art keywords
algorithm
particle
music
array
gravity
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
Application number
CN201810094288.1A
Other languages
English (en)
Other versions
CN108520195A (zh
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.)
Hubei University of Technology
Original Assignee
Hubei University of Technology
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 Hubei University of Technology filed Critical Hubei University of Technology
Priority to CN201810094288.1A priority Critical patent/CN108520195B/zh
Publication of CN108520195A publication Critical patent/CN108520195A/zh
Application granted granted Critical
Publication of CN108520195B publication Critical patent/CN108520195B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/12Classification; Matching
    • G06F2218/14Classification; Matching by matching peak patterns

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明提出一种基于引力搜索算法的MUSIC谱峰搜索方法,该方法将GSA算法和MUSIC算法相结合,具体方法是将MUSIC算法中的空间谱函数设为GSA算法的适应度函数,利用GSA算法优秀的全局寻优能力,迅速地找到空间谱函数的极值点,完成了对谱函数的空间搜索。仿真表明,该方法取得了较好的估计效果,在实验条件相同的情况下,与MUSIC算法中常用的均匀遍历法相比,该方法使得计算量大大降低。

Description

一种基于引力搜索算法的MUSIC谱峰搜索方法
技术领域
本发明涉及电子通信技术领域,具体涉及一种基于引力搜索算法的MUSIC谱峰搜索方法。
背景技术
阵列信号处理是电子通信领域的重点方向之一,在雷达、声纳和语音信号处理相关领域得到了广泛应用,最近以来引起了很高的重视,其发展重心主要集中波束形成和DOA估计的问题上。阵列信号处理是指在不同位置的多个传感器按特定的规则组成传感器阵列,利用传感器阵列估计信号的来波方向并且通过波束形成的方法对接收的信号进行处理,随着高分辨率谱估计技术的发展,DOA估计的问题也愈加受到重视。
经典的DOA估计方法中,以MUSIC(Multiple Signal Classification,多重信号分类)算法的应用最为广泛,该算法具有高精度和稳定性好的优点。MUSIC算法是一种典型的子空间结构类算法,其首先将接收数据构造的协方差矩阵分解为噪声子空间和信号子空间,然后利用噪声子空间与信号子空间的方向矢量之间的相互正交性来构造空间谱函数,谱函数对应的极值点即为参数的估计值。由于MUSIC算法具有较高的估计精度和稳定性,且突破了瑞利限制,目前已成为最具研究价值的DOA算法之一。传统的MUSIC算法中通常采用均匀遍历法计算谱函数,即对整个搜索空间进行遍历搜索,估计精度取决于间隔点间的步长大小,虽然该方法能将所有方向都能搜索到,但当精度要求较高时计算量也会相应变大,特别是对于多维、多目标的参数估计情况,MUSIC算法的谱峰搜索更加复杂,所需的时间也更多。由于其谱峰搜索部分的计算复杂度较高,使得该算法在实时性要求较高的场合的应用中受到限制。
MUSIC算法中谱峰搜索部分的计算量往往占该算法总计算量的一半以上,为了提高MUSIC算法的运算速度,使其满足实时性的要求,就要选取合适的方法进行MUSIC谱峰搜索。本发明主要针对MUSIC算法的谱峰搜索方法进行创新,在满足精度要求情况下提高搜索的速度。
下面对信号模型、MUSIC算法原理、引力搜索算法原理进行简单介绍。
1信号模型
假设阵列的信号接收模型为均匀线阵,且阵列天线由M个传感器构成,相邻传感器间隔为d,阵列外有D个波长为λ的窄带平面波,从不同入射角度入射,相邻传感器间隔d=λ/2。则阵列接收数据的矢量形式为:
X(t)=A(θ)S(t)+n(t) (1)
其中,X(t)=[x1(t),x2(t),...,xM(t)]T为M维数据接收向量,S(t)=[s1(t),s2(t),...,sD(t)]T为D维的信号源向量,n(t)是空间和时间均独立的高斯白噪声矢量,均值为0,方差为σ2。阵列流型矩阵A(θ)为M×D阶矩阵,其表达式为:
A(θ)=[a(θ1),a(θ2),...a(θ)D] (2)
其中,第k列a(θk)为第k个信号源sk的导向向量,即:
Figure BDA0001564600990000021
其中,j指虚数单位,
Figure BDA0001564600990000022
θk为第k个信号源的来波方向角。阵列数据的协方差矩阵为:
R=E[X(t)XH(t)]=ARSA+RN (4)
式中,E[·]和[·]H分别指期望值和共轭矩阵,A和X(t)分别指阵列流型矩阵和数据接收向量,RS和RN分别为信源的协方差矩阵和噪源的协方差矩阵。
2MUSIC算法原理
MUSIC算法是基于噪声子空间类的算法,它利用信号子空间和噪声子空间的正交性,构造空间谱函数,通过谱峰搜索,估计信号的来波方向。MUSIC算法是基于以下假设的数学模型上进行的:
假设1:阵列类型为各向同性均匀线阵,且阵元间隔一般小于最高频率信号波长的二分之一;
假设2:空间中的噪声为高斯白噪声,不同阵元间的噪声均为平稳随机过程,且相互独立,各阵元间的噪声均值为零、方差相等;
假设3:空间信号为窄带平面波,它与阵元噪声相互独立;
假设4:信号源数小于阵元数。
在此假设基础上,MUSIC算法对波达方向DOA的估计在理论上可以有任意高的分辨率。由于MUSIC算法的超分辨特点,该算法及其改进算法在近几十年来受到了广泛关注,并已成为超分辨子空间算法中的经典算法之一。
MUSIC算法的实现步骤为:
(1)根据传感器阵列接收的N次快拍数据构造数据协方差矩阵的估计值如下:
Figure BDA0001564600990000031
式中X(i)为第i次接收到的数据。
(2)对数据协方差矩阵进行特征值分解:
R=VΛVH (6)
其中,V是与R同阶的方阵,Λ是由特征值组成的对角矩阵。
(3)对分解得到的特征值按大小进行排序,把与信号个数D相等的特征值所对应的特征向量作为信号子空间,把剩下的(M-D)个特征值对应的特征向量作为噪声子空间:R=VSΛSVS H+VNΛNVN H,VS和VN为别是信号子空间和噪声子空间,并确定R的最小特征值个数nE
(4)构造空间谱函数:
Figure BDA0001564600990000032
其中,a(θ)为信号子空间的方向向量。
对谱函数进行空间搜索,其D个极大值对应的θ就是待估信源方向的估计值。
3引力搜索算法介绍
引力搜索算法(Gravitational Search Algorithm,GSA)是2009年由E.Rashedi根据万有引力现象和牛顿力学定理总结得出的一个仿生智能算法,该算法思想简单,搜索速度快,是全局优化问题的有效解决方案之一。
引力搜索算法(GSA)是一个较为常见的现代元启发式算法,根据算法的描述,有质量的个体之间产生相互作用力,质量较小的个体被吸引并产生加速度移动。每个个体有四个重要的参数:位置、惯性质量、主动惯性质量和被动惯性质量。惯性质量、主动惯性质量和被动惯性质量在算法中被认为是相等的。
与其他的群智能优化算法类似,该算法首先初始化种群的粒子,各个粒子之间存在相互作用力,随着算法的进行,这种作用力会促使所有质量小的粒子朝着相对而言质量更大的粒子移动。各粒子的速度和位置从初代开始改变,粒子所在的位置就是问题的可能解,每个粒子的性能与其惯性质量成一定的比例,当算法终止时,选择性能最好的粒子作为问题的解。
发明内容
本发明的目的是:提出一种基于引力搜索算法的MUSIC谱峰搜索方法,以期能降低运算复杂度,加快算法的执行速度,从而更好地应用于实际场景。
为了达到上述目的,本发明所采用的技术方案是:一种基于引力搜索算法的MUSIC谱峰搜索方法,其特征在于,将引力搜索算法用于MUSIC算法的谱峰搜索问题,则种群中粒子的位置代表着估计方位角在每次迭代过程中的解,所述方法包括如下步骤:
步骤一,初始化所有粒子在给定维度的速度和位置信息,并设置迭代次数和引力搜索算法中的初始参数,将MUSIC空间谱函数PMUSIC设为引力搜索算法的适应度函数,将粒子的位置作为自变量θ,代入MUSIC空间谱函数PMUSIC算出所有粒子的适应度值,并根据适应度值衡量粒子的性能;
步骤二,运用引力搜索算法在每一次迭代过程中更新重力常数,并利用适应度值算出每个粒子的惯性质量,根据重力常数和惯性质量计算出每个粒子的引力之和,再由引力之和和惯性质量得出每个粒子的加速度;
步骤三,利用加速度计算出每次迭代过程中所有粒子的速度和位置,适应度值最大的粒子的位置代表着DOA估计问题的解;
步骤四,检查算法的终止条件,如果满足,终止算法并选择适应度值最大的一组解作为DOA估计问题的解。
进一步地,所述基于引力搜索算法的MUSIC谱峰搜索方法的具体实现过程包括:
步骤1:随机地初始化种群中的粒子,每个粒子都有自己的初始位置,各粒子的位置Xi的表达式为:
Figure BDA0001564600990000051
其中,
Figure BDA0001564600990000052
和D分别表示粒子i在第k维搜索空间的位置和搜索空间的维数,D也等于待估计信源的个数,N为初始种群大小,算法在各搜索空间的搜索范围在[-90°,90°]之间;
步骤2:算出各粒子在不同时刻的惯性质量,计算公式如下:
Mai=Mpi=Mi,i=1,2,..,N (9)
Figure BDA0001564600990000053
Figure BDA0001564600990000054
Mai和Mpi分别代表第i个粒子主动惯性质量和被动惯性质量,Mi(t)代表第i个粒子的惯性质量,fiti(t)代表t时刻第i个粒子的适应度值,是根据MUSIC空间谱函数
Figure BDA0001564600990000055
得出的值,其中,a(θ)为信号子空间的方向向量,VN为噪声子空间,[·]H指共轭矩阵;
依据上述MUSIC空间谱函数PMUSIC公式计算每次迭代中所有粒子的适应度,由于该问题是优化最大值问题,时刻t的最优解best(t)和最差解worst(t)分别为:
Figure BDA0001564600990000056
Figure BDA0001564600990000057
步骤3:从引力定律中可以推出,第k维上一个粒子i在t时刻受到另一个粒子j的作用力的计算公式如下:
Figure BDA0001564600990000061
Rij(t)表示在t时刻两个粒子i和j的欧几里得度量,ε表示一个非常小的常数;
在t时刻的万有引力常数G(t)为:
G(t)=G0×e-α·t/T (18)
G0表示在t0时刻G的取值,一般取值为100,α为下降系数,一般取值为20,t代表当前循环次数,T代表最大循环次数;
则在t时刻,Xi在第k维上受到的引力总和是由种群中其余粒子对其产生的全部引力相加得到,计算公式如下所示:
Figure BDA0001564600990000062
式中,randj表示在[0,1]之间的随机数;
步骤4:从牛顿的力学原理可以得知,粒子产生的加速度等于该粒子受到的引力总和与其惯性质量之比,则粒子i的加速度的计算公式如下:
Figure BDA0001564600990000063
然后在每一次更新过程中,根据计算得到的加速度来更新每个粒子的速度,再根据粒子的速度更新每个粒子的位置,计算公式如下所示:
Figure BDA0001564600990000064
Figure BDA0001564600990000065
与现有技术相比,本发明的有益效果是:本发明提出的基于引力搜索算法的MUSIC谱峰搜索方法,该方法将GSA算法和MUSIC算法相结合,具体方法是将MUSIC算法中的空间谱函数设为GSA算法的适应度函数,利用GSA算法优秀的全局寻优能力,迅速地找到空间谱函数的极值点,完成了对谱函数的空间搜索。仿真表明,该方法取得了较好的估计效果,在实验条件相同的情况下,与MUSIC算法中常用的均匀遍历法相比,该方法使得计算量大大降低。
附图说明
图1为本发明基于引力搜索算法的MUSIC谱峰搜索方法的流程图。
图2为本发明实施例的迭代次数与DOA估计角度的变化曲线。
具体实施方式
为了便于本领域普通技术人员理解和实施本发明,下面结合实施例对本发明作进一步的详细描述,应当理解,此处所描述的实施示例仅用于说明和解释本发明,并不用于限定本发明。
基于引力搜索算法的MUSIC谱峰搜索步骤:
将GSA算法用于MUSIC算法的谱峰搜索问题,则种群中粒子的位置代表着估计方位角在每次迭代过程中的解,该方法的具体描述如下:
步骤1:随机地初始化种群中的粒子,每个粒子都有自己的初始位置,各粒子的位置Xi的表达式为:
Figure BDA0001564600990000071
其中,
Figure BDA0001564600990000072
和D分别表示粒子i在第k维搜索空间的位置和搜索空间的维数,D也等于待估计信源的个数,N为初始种群大小。算法在各搜索空间的搜索范围在[-90°,90°]之间。
步骤2:算出各粒子在不同时刻的惯性质量,计算公式如下:
Mai=Mpi=Mi,i=1,2,..,N (9)
Figure BDA0001564600990000073
Figure BDA0001564600990000074
Mai和Mpi分别代表第i个粒子主动惯性质量和被动惯性质量,Mi(t)代表第i个粒子的惯性质量,fiti(t)代表t时刻第i个粒子的适应度值,是根据MUSIC空间谱函数
Figure BDA0001564600990000081
得出的值。
依据公式(7)计算每次迭代中所有粒子的适应度,由于该问题是优化最大值问题,时刻t的最优解best(t)和最差解worst(t)分别为:
Figure BDA0001564600990000082
Figure BDA0001564600990000083
步骤3:从引力定律中可以推出,第k维上一个粒子i在t时刻受到另一个粒子j的作用力
Figure BDA0001564600990000084
的计算公式如下:
Figure BDA0001564600990000085
Rij(t)表示在t时刻两个粒子i和j的欧几里得度量,ε表示一个非常小的常数,
Figure BDA0001564600990000086
表示粒子i在第k维度t时刻的位置。在t时刻的万有引力常数G(t)为:
G(t)=G0×e-α·t/T (15)
G0表示在t0时刻G的取值,一般取值为100,α为下降系数,一般取值为20,t代表当前循环次数,T代表最大循环次数。
则在t时刻,Xi在第k维上受到的引力总和Fi k(t)是由种群中其余粒子对其产生的全部引力相加得到,计算公式如下所示:
Figure BDA0001564600990000087
式中,
Figure BDA0001564600990000088
为在第k维度中粒子j对粒子i的作用力,randj表示在[0,1]之间的随机数。
步骤4:从牛顿的力学原理可以得知,粒子产生的加速度等于该粒子受到的引力总和与其惯性质量之比,则第k维度的粒子i在t时刻的加速度
Figure BDA0001564600990000089
的计算公式如下:
Figure BDA00015646009900000810
然后在每一次更新过程中,根据计算得到的加速度来更新每个粒子的速度,再根据粒子的速度更新每个粒子的位置,则第k维度中粒子i的速度与位置更新公式如下:
Figure BDA0001564600990000091
Figure BDA0001564600990000092
式中
Figure BDA0001564600990000093
Figure BDA0001564600990000094
分别表示第i个粒子在第k维度的位置和速度。
下面给出该方法的总体思路:
首先初始化所有粒子在各维度的速度和位置信息,并设置迭代次数与算法中的其它参数,MUSIC空间谱函数
Figure BDA0001564600990000095
设为GSA算法的适应度函数。将粒子的位置作为自变量θ,并代入PMUSIC算出所有粒子的适应度值。适应度值是衡量粒子性能的重要指标。
然后在每一次迭代过程中更新重力常数,并利用适应度值算出每个粒子的惯性质量,根据重力常数和惯性质量求出每个粒子的引力之和,再由引力之和和惯性质量得到每个粒子的加速度。
最后利用加速度计算出每次迭代过程中所有粒子的速度和位置,由于本问题是优化最大值问题,所以适应度值最大的粒子的位置代表着DOA估计问题的解。
检查算法的终止条件,如果满足,终止算法并选择公式(19)中适应度值最大的一组解作为DOA估计问题的解。整个方法的流程图如图1所示。
仿真结果
仿真一:假设阵列传感器使用10元等距线阵,相邻传感器的间隔为来波信号波长的一半。并考虑空间中存在一个窄带平面波,其入射方向角是30°。信号的快拍数为100,信噪比(Signal-Noise Ratio,SNR)固定为10dB。在GSA算法中,各参数设置如下:初始种群大小N=10,G0=100,α=20,最大循环次数T=50。
图2给出了迭代次数与DOA估计角度的变化曲线。
由图2可以看到,在参数为一维的情况下,GSA算法迭代到第6次时即停止并成功估计出了单个来波信号的DOA,初始种群大小设置为10,则适应度函数公式(7)计算了10×6=60次,相比较传统谱峰搜索方法的1×180=180次,该方法使得计算量大大降低。
应当理解的是,本说明书未详细阐述的部分均属于现有技术。
应当理解的是,上述针对较佳实施例的描述较为详细,并不能因此而认为是对本发明专利保护范围的限制,本领域的普通技术人员在本发明的启示下,在不脱离本发明权利要求所保护的范围情况下,还可以做出替换或变形,均落入本发明的保护范围之内,本发明的请求保护范围应以所附权利要求为准。

Claims (2)

1.一种基于引力搜索算法的MUSIC谱峰搜索方法,应用于电子通信领域的阵列信号处理的DOA估计的问题上,解决MUSIC算法谱峰搜索部分的计算复杂度高的问题,其特征在于:
首先,定义信号模型:
假设阵列的信号接收模型为均匀线阵,且阵列天线由M个传感器构成,相邻传感器间隔为d,阵列外有D个波长为λ的窄带平面波,从不同入射角度入射,D表示信号个数,相邻传感器间隔d=λ/2,则阵列接收数据的矢量形式为:
X(t)=A(θ)S(t)+n(t) (1)
其中,t为时刻,X(t)=[x1(t),x2(t),...,xM(t)]T为M维数据接收向量,S(t)=[s1(t),s2(t),...,sD(t)]T为D维的信号源向量,n(t)是空间和时间均独立的高斯白噪声矢量,均值为0,方差为σ2,阵列流型矩阵A(θ)为M×D阶矩阵,其表达式为:
A(θ)=[a(θ1),a(θ2),...a(θ)D] (2)
其中,第k列a(θk)为第k个信号源sk的导向向量,即:
Figure FDA0002571769270000011
其中,j指虚数单位,
Figure FDA0002571769270000012
θk为第k个信号源的来波方向角,阵列数据的协方差矩阵为:
R=E[X(t)XH(t)]=ARSA+RN (4)
式中,E[·]和[·]H分别指期望值和共轭矩阵,A和X(t)分别指阵列流型矩阵和数据接收向量,RS和RN分别为信源的协方差矩阵和噪源的协方差矩阵;
其次,关于MUSIC算法:
MUSIC算法是基于噪声子空间类的算法,它利用信号子空间和噪声子空间的正交性,构造空间谱函数,通过谱峰搜索,估计信号的来波方向,MUSIC算法是基于以下假设的数学模型上进行的:
假设1:阵列类型为各向同性均匀线阵,且阵元间隔一般小于最高频率信号波长的二分之一;
假设2:空间中的噪声为高斯白噪声,不同阵元间的噪声均为平稳随机过程,且相互独立,各阵元间的噪声均值为零、方差相等;
假设3:空间信号为窄带平面波,它与阵元噪声相互独立;
假设4:信号源数小于阵元数;
在此假设基础上,MUSIC算法的实现步骤为:
(1)根据传感器阵列接收的N次快拍数据构造数据协方差矩阵的估计值如下:
Figure FDA0002571769270000021
式中X(i)为第i次接收到的数据;
(2)对数据协方差矩阵进行特征值分解:
R=VΛVH (6)
其中,V是与R同阶的方阵,Λ是由特征值组成的对角矩阵;
(3)对分解得到的特征值按大小进行排序,把与信号个数D相等的特征值所对应的特征向量作为信号子空间,把剩下的(M-D)个特征值对应的特征向量作为噪声子空间:R=VSΛSVS H+VNΛNVN H,VS和VN为别是信号子空间和噪声子空间,并确定R的最小特征值个数nE
(4)构造空间谱函数:
Figure FDA0002571769270000022
其中,a(θ)为信号子空间的方向向量;
对谱函数进行空间搜索,其D个极大值对应的θ就是待估信源方向的估计值;
再次,关于引力搜索算法:
引力搜索算法是一个现代元启发式算法,根据算法的描述,有质量的个体之间产生相互作用力,质量较小的个体被吸引并产生加速度移动;每个个体有四个重要的参数:位置、惯性质量、主动惯性质量和被动惯性质量,惯性质量、主动惯性质量和被动惯性质量在算法中被认为是相等的;
该算法首先初始化种群的粒子,各个粒子之间存在相互作用力,随着算法的进行,这种作用力会促使所有质量小的粒子朝着相对而言质量更大的粒子移动;各粒子的速度和位置从初代开始改变,粒子所在的位置就是问题的可能解,每个粒子的性能与其惯性质量成一定的比例,当算法终止时,选择性能最好的粒子作为问题的解;
所述基于引力搜索算法的MUSIC谱峰搜索方法,将引力搜索算法用于MUSIC算法的谱峰搜索问题,则种群中粒子的位置代表着估计方位角在每次迭代过程中的解,所述方法包括如下步骤:
步骤一,初始化所有粒子在给定维度的速度和位置信息,并设置迭代次数和引力搜索算法中的初始参数,将MUSIC空间谱函数PMUSIC设为引力搜索算法的适应度函数,将粒子的位置作为自变量θ,代入MUSIC空间谱函数PMUSIC算出所有粒子的适应度值,并根据适应度值衡量粒子的性能;
步骤二,运用引力搜索算法在每一次迭代过程中更新重力常数,并利用适应度值算出每个粒子的惯性质量,根据重力常数和惯性质量计算出每个粒子的引力之和,再由引力之和和惯性质量得出每个粒子的加速度;
步骤三,利用加速度计算出每次迭代过程中所有粒子的速度和位置,适应度值最大的粒子的位置代表着DOA估计问题的解;
步骤四,检查算法的终止条件,如果满足,终止算法并选择适应度值最大的一组解作为DOA估计问题的解。
2.如权利要求1所述的基于引力搜索算法的MUSIC谱峰搜索方法,其特征在于,所述方法的具体实现过程包括:
步骤1:随机地初始化种群中的粒子,每个粒子都有自己的初始位置,各粒子的位置Xi的表达式为:
Figure FDA0002571769270000031
其中,
Figure FDA0002571769270000032
和D分别表示粒子i在第k维搜索空间的位置和搜索空间的维数,D也等于待估计信源的个数,N为初始种群大小,算法在各搜索空间的搜索范围在[-90°,90°]之间;
步骤2:算出各粒子在不同时刻的惯性质量,计算公式如下:
Mai=Mpi=Mi,i=1,2,..,N (9)
Figure FDA0002571769270000041
Figure FDA0002571769270000042
Mai和Mpi分别代表第i个粒子主动惯性质量和被动惯性质量,Mi(t)代表第i个粒子的惯性质量,t为时刻,fiti(t)代表t时刻第i个粒子的适应度值,是根据MUSIC空间谱函数
Figure FDA0002571769270000043
得出的值,其中,a(θ)为信号子空间的方向向量,VN为噪声子空间,[·]H指共轭矩阵;
依据上述MUSIC空间谱函数PMUSIC公式计算每次迭代中所有粒子的适应度,由于该问题是优化最大值问题,时刻t的最优解best(t)和最差解worst(t)分别为:
Figure FDA0002571769270000044
Figure FDA0002571769270000045
步骤3:从引力定律中可以推出,第k维上一个粒子i在t时刻受到另一个粒子j的作用力的计算公式如下:
Figure FDA0002571769270000046
Rij(t)表示在t时刻两个粒子i和j的欧几里得度量,ε表示一个非常小的常数;
在t时刻的万有引力常数G(t)为:
G(t)=G0×e-α·t/T (18)
G0表示在t0时刻G的取值,取值为100,α为下降系数,取值为20,t对应着当前循环次数,T代表最大循环次数;
则在t时刻,Xi在第k维上受到的引力总和是由种群中其余粒子对其产生的全部引力相加得到,计算公式如下所示:
Figure FDA0002571769270000051
式中,randj表示在[0,1]之间的随机数;
步骤4:从牛顿的力学原理可以得知,粒子产生的加速度等于该粒子受到的引力总和与其惯性质量之比,则粒子i的加速度的计算公式如下:
Figure FDA0002571769270000052
然后在每一次更新过程中,根据计算得到的加速度来更新每个粒子的速度,再根据粒子的速度更新每个粒子的位置,计算公式如下所示:
Figure FDA0002571769270000053
Figure FDA0002571769270000054
CN201810094288.1A 2018-01-31 2018-01-31 一种基于引力搜索算法的music谱峰搜索方法 Active CN108520195B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810094288.1A CN108520195B (zh) 2018-01-31 2018-01-31 一种基于引力搜索算法的music谱峰搜索方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810094288.1A CN108520195B (zh) 2018-01-31 2018-01-31 一种基于引力搜索算法的music谱峰搜索方法

Publications (2)

Publication Number Publication Date
CN108520195A CN108520195A (zh) 2018-09-11
CN108520195B true CN108520195B (zh) 2020-12-29

Family

ID=63432732

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810094288.1A Active CN108520195B (zh) 2018-01-31 2018-01-31 一种基于引力搜索算法的music谱峰搜索方法

Country Status (1)

Country Link
CN (1) CN108520195B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113050037B (zh) * 2021-03-23 2022-10-04 上海交通大学 一种变电站设备异常声源定位方法及系统
CN116400028B (zh) * 2023-05-29 2023-08-22 湖南汇湘轩生物科技股份有限公司 一种基于气味传感器的香精质量检测方法、系统及介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2007101451A1 (en) * 2006-03-09 2007-09-13 Fundacio Privada Centre Tecnologic De Telecomunicacions De Catalunya Method and system for estimating directions-of-arrival in low power or low sample size scenarios
CN103901394A (zh) * 2014-03-21 2014-07-02 哈尔滨工程大学 一种冲击噪声环境下的量子万有引力搜索动态doa估计方法
CN105246005A (zh) * 2015-09-15 2016-01-13 国家电网公司 基于混合引力搜索算法的立体传声器阵列优化设计方法
CN107168050A (zh) * 2017-06-01 2017-09-15 武汉科技大学 一种引力搜索算法及基于该算法的污水处理控制方法
CN107451562A (zh) * 2017-07-31 2017-12-08 湖北工业大学 一种基于混沌二进制引力搜索算法的波段选择方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2007101451A1 (en) * 2006-03-09 2007-09-13 Fundacio Privada Centre Tecnologic De Telecomunicacions De Catalunya Method and system for estimating directions-of-arrival in low power or low sample size scenarios
CN103901394A (zh) * 2014-03-21 2014-07-02 哈尔滨工程大学 一种冲击噪声环境下的量子万有引力搜索动态doa估计方法
CN105246005A (zh) * 2015-09-15 2016-01-13 国家电网公司 基于混合引力搜索算法的立体传声器阵列优化设计方法
CN107168050A (zh) * 2017-06-01 2017-09-15 武汉科技大学 一种引力搜索算法及基于该算法的污水处理控制方法
CN107451562A (zh) * 2017-07-31 2017-12-08 湖北工业大学 一种基于混沌二进制引力搜索算法的波段选择方法

Also Published As

Publication number Publication date
CN108520195A (zh) 2018-09-11

Similar Documents

Publication Publication Date Title
CN109932680B (zh) 一种基于平移互质阵列的非圆信号波达方向估计方法
CN110031794B (zh) 一种基于差分共性阵重构的相干信源doa估计方法
CN108375751B (zh) 多信源波达方向估计方法
CN109490819B (zh) 一种基于稀疏贝叶斯学习的离格波达方向估计方法
CN107092004B (zh) 基于信号子空间旋转不变性的互质阵列波达方向估计方法
CN110007266B (zh) 一种冲击噪声下的任意阵列相干源测向方法
Bhuiya et al. Analysis of direction of arrival techniques using uniform linear array
CN108520195B (zh) 一种基于引力搜索算法的music谱峰搜索方法
CN107907855A (zh) 一种互素阵列转化为均匀线阵的doa估计方法及装置
WO2023070994A1 (zh) 一种阵列超分辨波达方向估计方法
CN104616059A (zh) 一种基于量子粒子群的波达方向估计方法
Gupta et al. Comparison of conventional and subspace based algorithms to estimate Direction of Arrival (DOA)
CN115119142A (zh) 一种基于传感器网络的分布式直接定位方法
KR101958337B1 (ko) 신호의 도래각을 추정하는 방법 및 장치
Magdy et al. Direction of arrival estimation based on maximum likelihood criteria using gravitational search algorithm
CN106646347B (zh) 基于小生境差分进化的多重信号分类谱峰搜索方法
CN115932714A (zh) 一种基于双曲正切核相关熵的单快拍波达方向估计方法
CN113341371B (zh) 一种基于l阵和二维esprit算法的doa估计方法
Sharma et al. Deterministic maximum likelihood direction of arrival estimation using GSA
CN114185002B (zh) 一种基于波束空间矩阵束的三维参数估计方法
CN106886627B (zh) 一种m-uca估计m-1个信源的建模方法
Kumar et al. Genetic Algorithm based optimization of uniform circular array
CN103792509B (zh) 电磁信号的二维波达方向角估计方法
CN113552530A (zh) 一种基于量子鼠群的近场和远场源混合测向方法
CN113093111A (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