CN111079615A - 一种基于莱维飞行细菌觅食进化的盲源分离方法 - Google Patents

一种基于莱维飞行细菌觅食进化的盲源分离方法 Download PDF

Info

Publication number
CN111079615A
CN111079615A CN201911258192.5A CN201911258192A CN111079615A CN 111079615 A CN111079615 A CN 111079615A CN 201911258192 A CN201911258192 A CN 201911258192A CN 111079615 A CN111079615 A CN 111079615A
Authority
CN
China
Prior art keywords
population
individual
evolution
maximum
formula
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
Application number
CN201911258192.5A
Other languages
English (en)
Other versions
CN111079615B (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.)
Harbin Engineering University
Original Assignee
Harbin Engineering 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 Harbin Engineering University filed Critical Harbin Engineering University
Priority to CN201911258192.5A priority Critical patent/CN111079615B/zh
Publication of CN111079615A publication Critical patent/CN111079615A/zh
Application granted granted Critical
Publication of CN111079615B publication Critical patent/CN111079615B/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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/004Artificial life, i.e. computing arrangements simulating life
    • G06N3/006Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Artificial Intelligence (AREA)
  • General Engineering & Computer Science (AREA)
  • Biomedical Technology (AREA)
  • Data Mining & Analysis (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Biophysics (AREA)
  • Computational Linguistics (AREA)
  • Signal Processing (AREA)
  • Evolutionary Computation (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Computing Systems (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明属于辐射源信号分离领域,目的在于提供一种基于莱维飞行细菌觅食进化的盲源分离方法:首先,对群智能交互机制参数和种群进行初始化,然后开始迭代和趋向循环,判断是否达到最大趋向次数,若未达到继续进行趋向操作,否则开始复制操作,完成复制操作后开始种群迁徙,完成迁徙后,对种群实施差分进化并进行多样性评估,若子代平均个体间距大于父代,则将子代种群替换父代;最后,判定是否达到最大迭代次数,若未达到则继续迭代,否则将最优参数组合输出至VMD中得到本征模态分量,相关运算后选取相关系数最大且与源信号具有相同数量的分量作为最终源信号。本发明提高了盲源分离的收敛速度及可靠性,增强了盲源分离的实时性及准确性。

Description

一种基于莱维飞行细菌觅食进化的盲源分离方法
技术领域
本发明属于辐射源信号分离领域,具体涉及一种基于莱维飞行细菌觅食进化的盲源分离方法。
背景技术
盲源分离是盲信号处理中极为重要的一个分支,是在已知源信号数目等少量先验信息或完全未知先验信息的前提下实现混合信号分离的信号处理方法,广泛应用于智能家居、生物医学以及无线通信等领域。盲源分离善于处理正定或超定的情况,即观测信号数目不少于源信号数目。然而,由于信号环境趋于复杂和接收设备日益微型化,使得观测信号的数目往往小于源信号的数目,甚至只有一路观测信号,故在多数情况下面临的是单通道欠定盲源分离问题。针对该问题应用最为广泛的是虚拟多通道技术,该技术借助经验模态分解(EmpiricalMode Decomposition,EMD)和变分模态分解(Variational ModeDecomposition,VMD)等算法将单通道转为多通道,再应用多通道盲源分离技术处理信号。
EMD将源信号分解为若干个本征模态分量之和,然后从源信号中提取每个本征模态分量。该方法适用于处理非线性、非平稳的时间序列,但是在处理间断信号和包络拟合时会出现模态混叠与端点效应问题,极大地降低了分离的准确性。
VMD作为一种新的自适应信号处理方法,将信号的分解问题转换成约束模型求取最优解的问题,可以有效的克服EMD的端点效应和模态混叠现象。该方法需要预先设定分解层数和惩罚因子,其中分解层数过大会导致过分解,反之则导致欠分解;惩罚因子过大导致分量信号带宽过窄,反之则导致分量信号过宽。此外,人为预设VMD参数随机性大且耗费资源多,导致分离的实时性与可靠性差。
发明内容
本发明的目的在于提供一种解决单通道欠定盲源分离存在收敛速度慢、准确性低和实时性差的基于莱维飞行细菌觅食进化的盲源分离方法。
为了实现本发明的目的,采用的技术方案为:
一种基于莱维飞行细菌觅食进化的盲源分离方法,包括如下步骤:
步骤一:对群智能交互机制参数和种群进行初始化,并设置最大总迭代次数、最大趋向循环次数和适应度函数,然后开始迭代;
步骤二:开始趋向循环过程,首先计算初始适应度值并保存,然后产生一个随机方向向量,让种群个体在该方向上以莱维飞行产生的随机步长进行游动,并根据适应度值的大小比较结果决定是否更换方向,最后达到最大转向次数后处理下一个种群个体,直到所有种群个体被处理;
步骤三:判断是否达到最大趋向次数,如果未达到则返回到第二步继续进行趋向操作,否则开始复制操作过程,首先计算每个种群个体适应度值的累加和(能量值)并进行排序,然后淘汰能量较差的半数种群,并依据能量较好的半数种群应用分布式算法估计再生得到新的半数种群,最后与存留下的半数种群组成新种群;
步骤四:完成复制操作后开始种群迁徙,每个种群个体在自适应概率下以高斯分布的形式进行迁徙;
步骤五:完成迁徙操作后,对种群实施差分进化并进行多样性评估,如果子代平均个体间距大于父代,则将子代种群替换原有父代,而原有父代舍弃;
步骤六:判定是否达到最大迭代次数,如果未达到则返回继续迭代,否则将寻优的最优参数组合输出至VMD中;
步骤七:VMD对从外界环境中观测到的辐射源信号进行分离操作;
步骤八:将本征模态分量与源信号做相关运算,选取相关系数最大且与源信号数目具有相同数量的分量作为最终分离出的源信号。
所述种群的初始化按照如下公式生成:
θi(1)=θmin+rand(-1,1)·(θmaxmin)
式中,θi(1)表示初始化第i个种群个体的位置,其中i=1,2,···,S;rand(-1,1)表示[-1,1]之间的随机数;θmin和θmax分别为优化区间的最小值与最大值。
所述初始适应度值的计算公式如下:
J(i,j)=fitness(θi(j))
式中,J(i,j)表示种群个体i在第j次趋向操作中的适应度值;θi(j)表示在第j次趋向操作后种群个体i的位置信息。
所述对能量较高的半数种群个体进行分布式估计再生的公式如下所示:
Figure BDA0002310886480000021
式中,θμ,q表示分布式估计再生的半数种群位置信息,η1和η2是[0,1]之间的均匀分布的随机数,μ和q分别是种群较优位置的分维度均值和标准差向量。
所述每个种群的自适应概率计算为:
Figure BDA0002310886480000031
式中,Pself(i)表示第i个种群个体的自适应迁徙概率;Ped为基本迁徙概率;
Figure BDA0002310886480000032
表示种群个体能量的最大值;
Figure BDA0002310886480000033
表示种群个体能量的最小值;rand(0,1)表示[0,1]之间的随机数。
所述以高斯分布的进行迁徙的公式为:
Figure BDA0002310886480000034
其中,
Figure BDA0002310886480000035
表示最优种群个体的第m维分量,N(0,2)表示均值为0,方差为2的高斯分布。
所述对种群进行差分变异的公式如下:
vi(g+1)=θr1(g)+F(θr2(g)-θr3(g))
式中,
Figure BDA0002310886480000036
表示一个新中间种群的位置信息;
Figure BDA0002310886480000037
表示第g代中第i个种群个体的位置信息;θr1(g)、θr2(g)和θr3(g)分别表示第g代中第r1、r2和r3个种群个体的位置信息;i≠r1≠r2≠r3,r1、r2和r3分别是区间[1,S]之间的随机整数;F为变异概率。
所述对种群进行多样性评估的计算公式为:
Figure BDA0002310886480000038
式中,L表示搜索空间的最大对角长度;
Figure BDA0002310886480000039
表示第g代种群中第i个个体位置的第m维分量;
Figure BDA00023108864800000310
表示第g代种群中所有个体位置的第m维平均值;Div(g)为第g代种群的平均个体间距。
所述VMD运算的计算公式为:
Figure BDA00023108864800000311
式中,Pn为An的归一化值;An为VMD在不同位置下分离出的本征模态分量经过Hilbert解调后的包络信号;n=1,2,···,M表示求取包络时将本征模态分量瞬时包络范围划分的等级数目;{Ed|d=1,2,···,K}表示An的包络熵;fitness(o)表示计算种群个体位置为o的适应度函数;min表示取最小值。
与现有技术相比,本发明的有益效果是:
本发明应用群智能交互机制、采用自适应概率迁徙方式、融合莱维飞行的方法、结合差分进化的思想及联合分布式估计的思想,极大的提高了盲源分离的收敛速度及可靠性,增强了盲源分离的实时性及准确性。
附图说明
图1为本发明的方案设计图。
具体实施方式
下面结合附图一对本发明做进一步说明。
本发明属于辐射源信号分离领域,具体涉及一种基于莱维飞行细菌觅食进化的盲源分离方法,目的在于提供一种解决单通道欠定盲源分离存在收敛速度慢、准确性低和实时性差的基于莱维飞行细菌觅食进化的盲源分离方法。
为了实现本发明的目的,采用的技术方案为:
一种基于莱维飞行细菌觅食进化的盲源分离方法,包括如下步骤:
步骤一:对群智能交互机制参数和种群进行初始化,并设置最大总迭代次数、最大趋向循环次数和适应度函数,然后开始迭代;
步骤二:开始趋向循环过程,首先计算初始适应度值并保存,然后产生一个随机方向向量,让种群个体在该方向上以莱维飞行产生的随机步长进行游动,并根据适应度值的大小比较结果决定是否更换方向,最后达到最大转向次数后处理下一个种群个体,直到所有种群个体被处理;
步骤三:判断是否达到最大趋向次数,如果未达到则返回到第二步继续进行趋向操作,否则开始复制操作过程,首先计算每个种群个体适应度值的累加和(能量值)并进行排序,然后淘汰能量较差的半数种群,并依据能量较好的半数种群应用分布式算法估计再生得到新的半数种群,最后与存留下的半数种群组成新种群;
步骤四:完成复制操作后开始种群迁徙,每个种群个体在自适应概率下以高斯分布的形式进行迁徙;
步骤五:完成迁徙操作后,对种群实施差分进化并进行多样性评估,如果子代平均个体间距大于父代,则将子代种群替换原有父代,而原有父代舍弃;
步骤六:判定是否达到最大迭代次数,如果未达到则返回继续迭代,否则将寻优的最优参数组合输出至VMD中;
步骤七:VMD对从外界环境中观测到的辐射源信号进行分离操作;
步骤八:将本征模态分量与源信号做相关运算,选取相关系数最大且与源信号数目具有相同数量的分量作为最终分离出的源信号。
所述种群的初始化按照如下公式生成:
θi(1)=θmin+rand(-1,1)·(θmaxmin)
式中,θi(1)表示初始化第i个种群个体的位置,其中i=1,2,···,S;rand(-1,1)表示[-1,1]之间的随机数;θmin和θmax分别为优化区间的最小值与最大值。
所述初始适应度值的尖酸公式如下:
J(i,j)=fitness(θi(j))
式中,J(i,j)表示种群个体i在第j次趋向操作中的适应度值;θi(j)表示在第j次趋向操作后种群个体i的位置信息。
所述对能量较高的半数种群个体进行分布式估计再生的公式如下所示:
Figure BDA0002310886480000051
式中,θμ,q表示分布式估计再生的半数种群位置信息,η1和η2是[0,1]之间的均匀分布的随机数,μ和q分别是种群较优位置的分维度均值和标准差向量。
所述每个种群的自适应概率计算为:
Figure BDA0002310886480000052
式中,Pself(i)表示第i个种群个体的自适应迁徙概率;Ped为基本迁徙概率;
Figure BDA0002310886480000053
表示种群个体能量的最大值;
Figure BDA0002310886480000054
表示种群个体能量的最小值;rand(0,1)表示[0,1]之间的随机数。
所述以高斯分布的进行迁徙的公式为:
Figure BDA0002310886480000055
其中,
Figure BDA0002310886480000056
表示最优种群个体的第m维分量,N(0,2)表示均值为0,方差为2的高斯分布。
所述对种群进行差分变异的公式如下:
vi(g+1)=θr1(g)+F(θr2(g)-θr3(g))
式中,
Figure BDA0002310886480000061
表示一个新中间种群的位置信息;
Figure BDA0002310886480000062
表示第g代中第i个种群个体的位置信息;θr1(g)、θr2(g)和θr3(g)分别表示第g代中第r1、r2和r3个种群个体的位置信息;i≠r1≠r2≠r3,r1、r2和r3分别是区间[1,S]之间的随机整数;F为变异概率。
所述对种群进行多样性评估的计算公式为:
Figure BDA0002310886480000063
式中,L表示搜索空间的最大对角长度;
Figure BDA0002310886480000064
表示第g代种群中第i个个体位置的第m维分量;
Figure BDA0002310886480000065
表示第g代种群中所有个体位置的第m维平均值;Div(g)为第g代种群的平均个体间距。
所述VMD运算的计算公式为:
Figure BDA0002310886480000066
式中,Pn为An的归一化值;An为VMD在不同位置下分离出的本征模态分量经过Hilbert解调后的包络信号;n=1,2,···,M表示求取包络时将本征模态分量瞬时包络范围划分的等级数目;{Ed|d=1,2,···,K}表示An的包络熵;fitness(o)表示计算种群个体位置为o的适应度函数;min表示取最小值。
进一步详细描述如下:
一种基于莱维飞行细菌觅食进化的盲源分离方法具体过程如下:
设定优化区间的最小值与最大值θmin和θmax、最大迭代次数Gmax、种群大小S、趋向操作最大次数Nc、一次趋向中种群个体的最大转向次数Ns、搜索空间维数D、基本迁徙概率Ped、变异概率F以及交叉概率CR等参数。以搜索范围内种群个体的位置信息(K,α)作为寻优问题的候选解,其中K表示分解层数,K0表示最优分解层数;α表示惩罚因子,α0表示最优惩罚因子。
步骤一:按照下式生成初始种群:
θi(1)=θmin+rand(-1,1)·(θmaxmin)
式中,θi(1)表示初始化第i个种群个体的位置,其中i=1,2,···,S;rand(-1,1)表示[-1,1]之间的随机数;θmin和θmax分别为优化区间的最小值与最大值。
在不同种群个体位置条件下做VMD运算,选取分离出分量的局部极小熵值作为信号的适应度函数,计算公式为:
Figure BDA0002310886480000071
式中,Pn为An的归一化值;An为VMD在不同位置下分离出的本征模态分量经过Hilbert解调后的包络信号;n=1,2,···,M表示求取包络时将本征模态分量瞬时包络范围划分的等级数目;{Ed|d=1,2,···,K}表示An的包络熵;fitness(o)表示计算种群个体位置为o的适应度函数;min表示取最小值。
步骤二:开始总迭代g=1,g表示当前迭代次数即第g代种群。
步骤三:开始趋向循环j=1,j表示当前趋向次数。
(1)计算种群个体初始位置对应的适应度值,计算公式如下:
J(i,j)=fitness(θi(j))
式中,J(i,j)表示种群个体i在第j次趋向操作中的适应度值;θi(j)表示在第j次趋向操作后种群个体i的位置信息,即寻优问题的候选解。
Figure BDA0002310886480000081
式中,Jcc(θ,R(j))表示在第j次趋向中种群个体之间感应值之和;Jcc(θ,θi(j))表示每两个种群个体之间的感应值;R(j)={θi(j)|i=1,2···,S}表示在第j次趋向操作后种群中每个成员的位置信息;θ=[θ12,···,θD]T表示D维搜索空间,(θ12,···,θD)是θ中的分量;D表示搜索空间的维度;S表示种群的大小;datt为引力的深度;watt为引力的宽度;hrep为斥力的高度;wrep为斥力的宽度;
Figure BDA0002310886480000082
为种群个体i位置信息的第m维分量且m=1,2,···,D;θm表示整个种群中其他个体的第m维分量;一般情况下,取datt=hrep
(2)存储种群个体i当前最好的适应度值Jlast=J(i,j),其中Jlast表示种群个体目前最好的适应度值。
(3)令种群个体i以步骤三中第(3)-(5)步所述方式趋向一步。种群个体i生成一个随机方向向量Δ(i),其中的每一个元素都是分布在[-1,1]上的随机数。
(4)设置种群个体方向转向次数l=1。结合莱维飞行的思想,生成第i个种群个体游动的步长C(i)。
Figure BDA0002310886480000083
式中,ε~N(0,σ2),δ~N(0,1),N(·)表示高斯分布,即ε服从均值为0、方差为σ2的高斯分布,δ服从均值为0、方差为1的高斯分布;σ为ε的标准差,β一般取[1,3]之间的随机数。
Figure BDA0002310886480000084
式中,Γ为标准的Gamma函数。
(5)种群个体在选定方向后,以步骤三第(4)步中的步长C(i)进行游动。
θi(j+1)=θi(j)+C(i)Φ(i)
式中,Φ(i)表示种群个体i游动的方向。
Figure BDA0002310886480000091
式中,ΔT(i)是方向向量Δ(i)的转置。
(6)对游动后的种群个体计算其适应度值J(i,j+1),如果J(i,j+1)<Jlast,则令Jlast=J(i,j+1),且该种群个体返回步骤三中第(4)步继续沿此方向进行莱维飞行式游动,否则l=l+1返回步骤三中第(3)步转换方向,直至l=Ns后返回步骤三中第(1)步处理种群个体i+1,直到所有种群个体被处理,此时所有种群个体完成一个趋向操作,趋向次数j=j+1。
步骤四:如果j≤Nc,返回步骤三继续趋向操作,否则进行下一步。
步骤五:计算种群个体的适应度值的累加和为种群个体的能量值。
Figure BDA0002310886480000092
式中,Nc为趋向行为的次数,
Figure BDA0002310886480000093
表示第i个种群个体的能量值。
步骤六:对种群个体能量值排序,淘汰能量较低的半数种群个体,对能量较高的半数种群个体进行分布式估计再生。按照如下公式进行复制再生:
Figure BDA0002310886480000094
式中,θμ,q表示分布式估计再生的半数种群位置信息,η1和η2是[0,1]之间的均匀分布的随机数,μ和q分别是种群较优位置的分维度均值和标准差向量。
步骤七:经过复制操作后,种群个体采用自适应概率迁徙操作,每个种群个体以高斯分布的形式重新分布,即判断随机数rand(0,1)是否大于种群个体i的自适应迁徙概率Pself(i)。其自适应概率计算为:
Figure BDA0002310886480000095
式中,Pself(i)表示第i个种群个体的自适应迁徙概率;Ped为基本迁徙概率;
Figure BDA0002310886480000096
表示种群个体能量的最大值;
Figure BDA0002310886480000101
表示种群个体能量的最小值;rand(0,1)表示[0,1]之间的随机数。
若rand(0,1)>Pself(i)则该种群个体位置信息不变,否则以如下公式进行迁徙。
Figure BDA0002310886480000102
其中,
Figure BDA0002310886480000103
表示最优种群个体的第m维分量,N(0,2)表示均值为0,方差为2的高斯分布。
步骤八:经过迁徙操作后,为减小种群的退化程度,采用差分进化的思想增加种群的多样性。
(1)以如下公式对种群进行差分变异得到一个新的中间种群。
vi(g+1)=θr1(g)+F(θr2(g)-θr3(g))
式中,
Figure BDA0002310886480000104
表示一个新中间种群的位置信息;
Figure BDA0002310886480000105
表示第g代中第i个种群个体的位置信息;θr1(g)、θr2(g)和θr3(g)分别表示第g代中第r1、r2和r3个种群个体的位置信息;i≠r1≠r2≠r3,r1、r2和r3分别是区间[1,S]之间的随机整数;F为变异概率。
(2)对g代种群θi(g)和其变异的中间种群vi(g+1)进行相应个体之间的交叉操作,具体如下:
Figure BDA0002310886480000106
式中,CR是交叉概率;
Figure BDA0002310886480000107
Figure BDA0002310886480000108
分别表示ui(g+1)和vi(g+1)中的第m维分量;
Figure BDA0002310886480000109
表示第g+1代新种群中第i个个体的位置信息;
Figure BDA00023108864800001010
表示
Figure BDA00023108864800001011
中的第m维分量。
(3)对种群进行多样性评估,其计算公式为:
Figure BDA00023108864800001012
式中,L表示搜索空间的最大对角长度;
Figure BDA00023108864800001013
表示第g代种群中第i个个体位置的第m维分量;
Figure BDA00023108864800001014
表示第g代种群中所有个体位置的第m维平均值;Div(g)为第g代种群的平均个体间距,该值反映种群离散程度并与种群多样性成正比。
(4)种群选择机制如下:
Figure BDA0002310886480000111
式中,θi(g+1)为g+1代种群中第i个种群个体的位置信息。即如果子代平均个体间距Div(g+1)大于父代平均个体间距Div(g),则舍弃原有父代种群替换为子代种群,此时种群完成一个总迭代g=g+1。
步骤九:判断是否达到最大总迭代次数Gmax。如果g≤Gmax,则返回步骤二,否则输出适应度值最小的种群个体位置信息,即最优候选解(K00)作为VMD所需参数组合。
步骤十:VMD首先借助群智能交互机制寻优得到最优参数组合(K00)和从外界环境中观测到的辐射源信号进行分离操作,得到K0个本征模态分量,然后将该分量与源信号做相关运算,最后选取相关系数最大的前λ(源信号数目λ≤K0)个分量作为最终分离出的信号。
综上所述,本发明属于辐射源信号分离领域,目的在于提供一种基于莱维飞行细菌觅食进化的盲源分离方法:首先,对群智能交互机制参数和种群进行初始化,然后开始迭代和趋向循环,判断是否达到最大趋向次数,若未达到继续进行趋向操作,否则开始复制操作,完成复制操作后开始种群迁徙,完成迁徙后,对种群实施差分进化并进行多样性评估,若子代平均个体间距大于父代,则将子代种群替换父代;最后,判定是否达到最大迭代次数,若未达到则继续迭代,否则将最优参数组合输出至VMD中得到本征模态分量,相关运算后选取相关系数最大且与源信号具有相同数量的分量作为最终源信号。本发明提高了盲源分离的收敛速度及可靠性,增强了盲源分离的实时性及准确性。

Claims (9)

1.一种基于莱维飞行细菌觅食进化的盲源分离方法,其特征在于,包括如下步骤:
步骤一:对群智能交互机制参数和种群进行初始化,并设置最大总迭代次数、最大趋向循环次数和适应度函数,然后开始迭代;
步骤二:开始趋向循环过程,首先计算初始适应度值并保存,然后产生一个随机方向向量,让种群个体在该方向上以莱维飞行产生的随机步长进行游动,并根据适应度值的大小比较结果决定是否更换方向,最后达到最大转向次数后处理下一个种群个体,直到所有种群个体被处理;
步骤三:判断是否达到最大趋向次数,如果未达到则返回到第二步继续进行趋向操作,否则开始复制操作过程,首先计算每个种群个体适应度值的累加和(能量值)并进行排序,然后淘汰能量较差的半数种群,并依据能量较好的半数种群应用分布式算法估计再生得到新的半数种群,最后与存留下的半数种群组成新种群;
步骤四:完成复制操作后开始种群迁徙,每个种群个体在自适应概率下以高斯分布的形式进行迁徙;
步骤五:完成迁徙操作后,对种群实施差分进化并进行多样性评估,如果子代平均个体间距大于父代,则将子代种群替换原有父代,而原有父代舍弃;
步骤六:判定是否达到最大迭代次数,如果未达到则返回继续迭代,否则将寻优的最优参数组合输出至VMD中;
步骤七:VMD对从外界环境中观测到的辐射源信号进行分离操作;
步骤八:将本征模态分量与源信号做相关运算,选取相关系数最大且与源信号数目具有相同数量的分量作为最终分离出的源信号。
2.根据权利要求1所述的一种基于莱维飞行细菌觅食进化的盲源分离方法,其特征在于,所述的所述种群的初始化按照如下公式生成:
θi(1)=θmin+rand(-1,1)·(θmaxmin)
式中,θi(1)表示初始化第i个种群个体的位置,其中i=1,2,···,S;rand(-1,1)表示[-1,1]之间的随机数;θmin和θmax分别为优化区间的最小值与最大值。
3.根据权利要求1所述的一种基于莱维飞行细菌觅食进化的盲源分离方法,其特征在于,所述初始适应度值的计算公式如下:
J(i,j)=fitness(θi(j))
式中,J(i,j)表示种群个体i在第j次趋向操作中的适应度值;θi(j)表示在第j次趋向操作后种群个体i的位置信息。
4.根据权利要求1所述的一种基于莱维飞行细菌觅食进化的盲源分离方法,其特征在于,所述对能量较高的半数种群个体进行分布式估计再生的公式如下所示:
Figure FDA0002310886470000021
式中,θμ,q表示分布式估计再生的半数种群位置信息,η1和η2是[0,1]之间的均匀分布的随机数,μ和q分别是种群较优位置的分维度均值和标准差向量。
5.根据权利要求1所述的一种基于莱维飞行细菌觅食进化的盲源分离方法,其特征在于,所述每个种群的自适应概率计算为:
Figure FDA0002310886470000022
式中,Pself(i)表示第i个种群个体的自适应迁徙概率;Ped为基本迁徙概率;
Figure FDA0002310886470000023
表示种群个体能量的最大值;
Figure FDA0002310886470000024
表示种群个体能量的最小值;rand(0,1)表示[0,1]之间的随机数。
6.根据权利要求1所述的一种基于莱维飞行细菌觅食进化的盲源分离方法,其特征在于,所述以高斯分布的进行迁徙的公式为:
Figure FDA0002310886470000025
其中,
Figure FDA0002310886470000026
表示最优种群个体的第m维分量,N(0,2)表示均值为0,方差为2的高斯分布。
7.根据权利要求1所述的一种基于莱维飞行细菌觅食进化的盲源分离方法,其特征在于,所述对种群进行差分变异的公式如下:
vi(g+1)=θr1(g)+F(θr2(g)-θr3(g))
式中,
Figure FDA0002310886470000027
表示一个新中间种群的位置信息;
Figure FDA0002310886470000028
表示第g代中第i个种群个体的位置信息;θr1(g)、θr2(g)和θr3(g)分别表示第g代中第r1、r2和r3个种群个体的位置信息;i≠r1≠r2≠r3,r1、r2和r3分别是区间[1,S]之间的随机整数;F为变异概率。
8.根据权利要求1所述的一种基于莱维飞行细菌觅食进化的盲源分离方法,其特征在于,所述对种群进行多样性评估的计算公式为:
Figure FDA0002310886470000029
式中,L表示搜索空间的最大对角长度;
Figure FDA0002310886470000031
表示第g代种群中第i个个体位置的第m维分量;
Figure FDA0002310886470000032
表示第g代种群中所有个体位置的第m维平均值;Div(g)为第g代种群的平均个体间距。
9.根据权利要求1所述的一种基于莱维飞行细菌觅食进化的盲源分离方法,其特征在于,所述VMD运算的计算公式为:
Figure FDA0002310886470000033
式中,Pn为An的归一化值;An为VMD在不同位置下分离出的本征模态分量经过Hilbert解调后的包络信号;n=1,2,···,M表示求取包络时将本征模态分量瞬时包络范围划分的等级数目;{Ed|d=1,2,···,K}表示An的包络熵;fitness(o)表示计算种群个体位置为o的适应度函数;min表示取最小值。
CN201911258192.5A 2019-12-10 2019-12-10 一种基于莱维飞行细菌觅食进化的盲源分离方法 Active CN111079615B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911258192.5A CN111079615B (zh) 2019-12-10 2019-12-10 一种基于莱维飞行细菌觅食进化的盲源分离方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911258192.5A CN111079615B (zh) 2019-12-10 2019-12-10 一种基于莱维飞行细菌觅食进化的盲源分离方法

Publications (2)

Publication Number Publication Date
CN111079615A true CN111079615A (zh) 2020-04-28
CN111079615B CN111079615B (zh) 2023-03-31

Family

ID=70313982

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911258192.5A Active CN111079615B (zh) 2019-12-10 2019-12-10 一种基于莱维飞行细菌觅食进化的盲源分离方法

Country Status (1)

Country Link
CN (1) CN111079615B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111865452A (zh) * 2020-06-29 2020-10-30 哈尔滨工程大学 一种基于自适应细菌觅食的单通道盲源分离方法
CN112199996A (zh) * 2020-09-04 2021-01-08 西安交通大学 基于参数自适应VMD及快速Hoyer谱图指标的滚动轴承诊断方法
CN115017934A (zh) * 2022-03-21 2022-09-06 哈尔滨工程大学 种群熵量子雷莫拉鱼搜索机制的欠定盲源分离方法
CN116165542A (zh) * 2023-03-01 2023-05-26 上海玫克生储能科技有限公司 一种电池参数的辨识方法、装置、设备及存储介质

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103530709A (zh) * 2013-11-04 2014-01-22 上海海事大学 一种基于细菌觅食优化方法的集装箱码头泊位和岸桥分配方法
US20150269314A1 (en) * 2014-03-20 2015-09-24 Rudjer Boskovic Institute Method and apparatus for unsupervised segmentation of microscopic color image of unstained specimen and digital staining of segmented histological structures
CN106898361A (zh) * 2017-03-16 2017-06-27 杭州电子科技大学 基于反馈变分模式分解的单通道盲源分离方法
CN109409194A (zh) * 2018-08-30 2019-03-01 中国航发湖南动力机械研究所 多模态时域信号模态分离、阻尼参数辨识方法及存储介质

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103530709A (zh) * 2013-11-04 2014-01-22 上海海事大学 一种基于细菌觅食优化方法的集装箱码头泊位和岸桥分配方法
US20150269314A1 (en) * 2014-03-20 2015-09-24 Rudjer Boskovic Institute Method and apparatus for unsupervised segmentation of microscopic color image of unstained specimen and digital staining of segmented histological structures
CN106898361A (zh) * 2017-03-16 2017-06-27 杭州电子科技大学 基于反馈变分模式分解的单通道盲源分离方法
CN109409194A (zh) * 2018-08-30 2019-03-01 中国航发湖南动力机械研究所 多模态时域信号模态分离、阻尼参数辨识方法及存储介质

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
PRIYANKA DEY,ET AL: "Single channel blind source separation based on variational mode decomposition and PCA", 《2015 ANNUAL IEEE INDIA CONFERENCE (INDICON)》 *
张前图等: "具有Levy飞行特征的双子群果蝇优化算法", 《计算机应用》 *
汤杰等: "基于IVMD的单通道盲源分离方法及其应用", 《组合机床与自动化加工技术》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111865452A (zh) * 2020-06-29 2020-10-30 哈尔滨工程大学 一种基于自适应细菌觅食的单通道盲源分离方法
CN111865452B (zh) * 2020-06-29 2022-04-05 哈尔滨工程大学 一种基于自适应细菌觅食的单通道盲源分离方法
CN112199996A (zh) * 2020-09-04 2021-01-08 西安交通大学 基于参数自适应VMD及快速Hoyer谱图指标的滚动轴承诊断方法
CN115017934A (zh) * 2022-03-21 2022-09-06 哈尔滨工程大学 种群熵量子雷莫拉鱼搜索机制的欠定盲源分离方法
CN116165542A (zh) * 2023-03-01 2023-05-26 上海玫克生储能科技有限公司 一种电池参数的辨识方法、装置、设备及存储介质
CN116165542B (zh) * 2023-03-01 2023-10-20 上海玫克生储能科技有限公司 一种电池参数的辨识方法、装置、设备及存储介质

Also Published As

Publication number Publication date
CN111079615B (zh) 2023-03-31

Similar Documents

Publication Publication Date Title
CN111079615B (zh) 一种基于莱维飞行细菌觅食进化的盲源分离方法
Yang et al. A robust EM clustering algorithm for Gaussian mixture models
CN110866526A (zh) 图像分割方法、电子设备及计算机可读存储介质
CN106203625B (zh) 一种基于多重预训练的深层神经网络训练方法
CN106886986B (zh) 基于自适应组结构稀疏字典学习的图像融合方法
CN111192211B (zh) 一种基于单个深度神经网络的多噪声类型盲去噪方法
CN110516537B (zh) 一种基于自步学习的人脸年龄估计方法
JP2015052832A (ja) 重み設定装置および方法
CN113298191B (zh) 基于个性化半监督在线联邦学习的用户行为识别方法
CN109697451A (zh) 相似图像聚类方法及装置、存储介质、电子设备
CN111950455A (zh) 一种基于lffcnn-gru算法模型的运动想象脑电信号特征识别方法
CN110263236B (zh) 基于动态多视图学习模型的社交网络用户多标签分类方法
WO2022126448A1 (zh) 一种基于演化学习的神经网络结构搜索方法和系统
CN108921037A (zh) 一种基于BN-inception双流网络的情绪识别方法
TWI512685B (zh) 移動物體偵測方法及其裝置
CN110599443A (zh) 一种使用双向长短期记忆网络的视觉显著性检测方法
Wu et al. Image denoising with rectified linear units
Liu et al. Deep MIMO detection based on belief propagation
CN107240100B (zh) 一种基于遗传算法的图像分割方法和系统
CN115330620A (zh) 一种基于循环生成对抗网络的图像去雾方法
CN116862024A (zh) 基于聚类和知识蒸馏的可信个性化联邦学习方法和装置
CN114189899B (zh) 一种基于随机聚合波束成形的用户设备选择方法
CN114117945A (zh) 基于用户-服务交互图的深度学习云服务QoS预测方法
CN110995375B (zh) 一种无线信道的衰落特征提取方法及提取装置
CN116543259A (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