CN110009088A - 一种基于种群熵、种群方差改进的钢轨裂纹检测算法 - Google Patents
一种基于种群熵、种群方差改进的钢轨裂纹检测算法 Download PDFInfo
- Publication number
- CN110009088A CN110009088A CN201910058482.9A CN201910058482A CN110009088A CN 110009088 A CN110009088 A CN 110009088A CN 201910058482 A CN201910058482 A CN 201910058482A CN 110009088 A CN110009088 A CN 110009088A
- Authority
- CN
- China
- Prior art keywords
- population
- individual
- value
- entropy
- variance
- 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
- 238000001514 detection method Methods 0.000 title claims abstract description 18
- 230000035772 mutation Effects 0.000 claims description 20
- 238000000034 method Methods 0.000 claims description 15
- 238000000354 decomposition reaction Methods 0.000 claims description 8
- 108090000623 proteins and genes Proteins 0.000 claims description 6
- 239000006185 dispersion Substances 0.000 claims description 3
- 230000014759 maintenance of location Effects 0.000 claims description 3
- 230000002068 genetic effect Effects 0.000 abstract description 21
- 238000003745 diagnosis Methods 0.000 abstract description 3
- 230000002028 premature Effects 0.000 abstract description 3
- 230000007423 decrease Effects 0.000 abstract description 2
- 238000012360 testing method Methods 0.000 description 6
- 230000003044 adaptive effect Effects 0.000 description 3
- 238000013528 artificial neural network Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 239000012141 concentrate Substances 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000009659 non-destructive testing Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000005215 recombination Methods 0.000 description 1
- 230000006798 recombination Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000013024 troubleshooting Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial 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)
- Data Mining & Analysis (AREA)
- General Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Evolutionary Computation (AREA)
- Artificial Intelligence (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Health & Medical Sciences (AREA)
- Complex Calculations (AREA)
Abstract
本发明公开了一种基于种群熵、种群方差改进的钢轨裂纹检测算法,提出将Sigmoid曲线与种群熵、种群方差结合起来改进遗传算子,进而推导出基于种群熵以及种群方差的遗传算法,最后用于钢轨裂纹故障诊断。该基于种群熵、种群方差改进的钢轨裂纹检测算法,解决了现有算法中种群多样性的急剧下降导致的早熟问题。
Description
技术领域
本发明属于钢轨裂纹无损检测算法技术领域,特别是涉及一种基于种群熵、种群方差改进的钢轨裂纹检测算法。
背景技术
铁路交通是现代人类社会必不可缺的运输工具。随着科学技术的日新月异,铁路交通呈现高速和重载两大发展趋势,因此在列车高速、货车重载运行的情况下,钢轨裂纹等缺陷更易导致重大的安全事故,甚至造成重大人员伤亡等灾难性的后果。
针对故障识别常用的常规方法主要有:时域特征法、频域特征法、小波包分解法等。时域的能量幅值判据,当信号的幅值大于某一设定值则认为是故障;又如频域的能量判据,当某一特定的频谱带中出现一定能量的峰值时则认为是故障;再如小波包分解法,通过小波包分解将信号分解到特定的小波层,再通过判断该层上能量是否超过设定值,来判别系统是否存在故障。但是,上述故障识别方法存在的主要缺点是依赖于人工经验来设定判据参数(如阈值、频带、小波层数等)。而采用常规故障诊断方法进行钢轨裂纹轮轨噪声诊断的人工经验难以获得,导致其判据参数无法确定,因此需要利用人工智能算法进行改进。目前,人工智能算法主要有遗传算法改进的小波包分解法、遗传算法改进的人工神经网络算法、遗传算法优化的小波神经网络算法等,其能够通过数据自动学习获取到判据的参数。但是,普通遗传算法存在效率低、容易早熟、易陷入局部最优解等缺点,导致了上述三种方法总体效率不高、诊断准确率偏低等缺点。
对于遗传算法而言,无论是早熟问题还是陷入局部最优解的问题,都是因为群体中个体结构多样性的急剧减少而造成的,所以就要在进化过程中增加种群的多样性。
发明内容
本发明的目的是提供一种基于种群熵、种群方差改进的钢轨裂纹检测算法,解决了现有算法中种群多样性的急剧下降导致的早熟问题。
为了达到上述目的,本发明所采用的技术方案是,一种基于种群熵、种群方差改进的钢轨裂纹检测算法,具体按照以下步骤依次实施:
步骤1,设置算法的参数;
步骤2,初始化原种群,得到初始化种群;
步骤3,计算步骤2得到的初始化种群中个体的适应度值;
步骤4,根据步骤3得到的个体的适应度值,求出满足设定值的最优个体并加以保存;
步骤5,根据步骤4得到的最优个体的总数,与目标值作比较,达到要求则进行步骤10,反之依次重复进行步骤6、步骤7、步骤8、步骤9、步骤3、步骤4,直至达到要求;
步骤6,采用轮赌盘方法对步骤4得到的初始化种群进行选择操作;
步骤7,对步骤6得到的初始化种群进行交叉操作;
步骤8,对步骤7得到的初始化种群进行变异操作;
步骤9,计算步骤8得到的初始化种群中最优个体和最差个体,采用精英保留策略更新初始化种群,之后返回步骤3;
步骤10,输出结果,结束程序。
本发明的技术方案,还具有以下特点:
在所述步骤1中,设置算法的参数为:原种群规模popsize=50,最小交叉概率Pcmin=0.6,最大交叉概率Pcmax=0.9,最小变异概率Pmmin=0.01,最大变异概率Pmmax=0.1;
3.根据权利要求2所述的基于种群熵、种群方差改进的钢轨裂纹检测算法,其特征在于,所述步骤2具体按照以下步骤依次实施:
步骤2.1,利用故障信号的时域、频域特征结合起来构建钢轨裂纹的故障判据如公式所示:
其中:xi为待优化系数,Ei为三层小波包分解第三层第i频段的能量函数,为输入信号样本数据,Amax为时域幅值函数,k为权重系数;
根据故障判据条件,得出下列故障判据函数如公式:
将m组样本数据求得的故障判据函数值与数据样本属性值相比较,如果故障判据函数值与数据样本属性值相同,说明判断正确,代价函数值增加;数据样本属性函数和代价函数分别描述为:
4.根据权利要求3所述的基于种群熵、种群方差改进的钢轨裂纹检测算法,其特征在于,在所述步骤3中:个体的适应度值根据公式(4)得到。
5.根据权利要求1所述的基于种群熵、种群方差改进的钢轨裂纹检测算法,其特征在于,所述步骤7和步骤8具体为:
设进化过程中第t代原种群P(t)中第i个体为N为原种群规模,个体的适应度为ft i,第t代种群个体的平均适应度为那么第t代种群的方差为:
式中:方差Vt反映了种群个体的离散程度,Vt的值大就说明个体比较分散,Vt的值小说明个体比较集中;
将信息熵引入种群多样性的判定中,用以反映种群内部进化状态;将实际问题中可行域等分为L个小区间,第i个小区间为Ai(i=1,…,L),在第t代进化中,落在中的个体数目记为定义第t代种群的熵为:
式中:表示每个小区间内的个体占种群规模N的概率大小,其中Ht的取值范围为(0,logN);
结合种群方差、种群熵在进化过程中的特征,种群进化到第t代的交叉变异公式如式(7)、(8),种群方差、种群熵联合方程如公式(9):
式中:a=9.903438,Pcmax是最大交叉概率,Pcmin是最小交叉概率,Pmmax是最大变异概率,Pmmin是最小变异概率,logN是种群熵理论上的最大值,max{Vt}是种群进化到第t代时种群方差的最大值;
根据交叉概率Pc判断是否交叉,交叉的子代采用一致交叉算子,不交叉的子代从父代获取基因;根据变异概率Pm判断是否变异,变异采用一致变异操作,不发生变异的子代从父代获取基因。
附图说明
图1是本发明的一种基于种群熵、种群方差改进的钢轨裂纹检测算法的工作原理图;
图2是本发明的一种基于种群熵、种群方差改进的钢轨裂纹检测算法的自适应调整曲线图;
图3是测试函数f1在不同迭代次数下采用本发明算法、自适应遗传算法以及传统遗传算法的总收敛次数比较图;
图4是测试函数f2在不同迭代次数下采用本发明算法、自适应遗传算法以及传统遗传算法的总收敛次数比较图。
具体实施方式
以下结合附图说明和具体实施例对本发明的技术方案作进一步地详细说明。
如图1所示,一种基于种群熵、种群方差改进的钢轨裂纹检测算法,具体按照以下步骤依次实施:
步骤1,设置算法的参数,该设置算法的参数为:原种群规模popsize=50,最小交叉概率Pcmin=0.6,最大交叉概率Pcmax=0.9,最小变异概率Pmmin=0.01,最大变异概率Pmmax=0.1;
步骤2,初始化原种群,得到初始化种群;步骤2具体按照以下步骤依次实施:
步骤2.1,利用故障信号的时域、频域特征结合起来构建钢轨裂纹的故障判据如公式所示:
其中:xi为待优化系数,Ei为三层小波包分解第三层第i频段的能量函数,为输入信号样本数据,Amax为时域幅值函数,k为权重系数;
根据故障判据条件,得出下列故障判据函数如公式:
将m组样本数据求得的故障判据函数值与数据样本属性值相比较,如果故障判据函数值与数据样本属性值相同,说明判断正确,代价函数值增加;数据样本属性函数和代价函数分别描述为:
步骤3,以公式(4)为依据,计算步骤2得到的初始化种群中个体的适应度值;
步骤4,根据步骤3得到的个体的适应度值,求出满足设定值的最优个体并加以保存;
步骤5,根据步骤4得到的最优个体的总数,与目标值作比较,达到要求则进行步骤10,反之依次重复进行步骤6、步骤7、步骤8、步骤9、步骤3、步骤4,直至达到要求;
步骤6,采用轮赌盘方法对步骤4得到的初始化种群进行选择操作;
步骤7,对步骤6得到的初始化种群进行交叉操作;
步骤8,对步骤7得到的初始化种群进行变异操作;
步骤7和步骤8具体为:
设进化过程中第t代原种群P(t)中第i个体为N为原种群规模,个体的适应度为ft i,第t代种群个体的平均适应度为那么第t代种群的方差为:
式中:方差Vt反映了种群个体的离散程度,Vt的值大就说明个体比较分散,Vt的值小说明个体比较集中;
将信息熵引入种群多样性的判定中,用以反映种群内部进化状态;将实际问题中可行域等分为L个小区间,第i个小区间为Ai(i=1,…,L),在第t代进化中,落在中的个体数目记为定义第t代种群的熵为:
式中:表示每个小区间内的个体占种群规模N的概率大小,其中Ht的取值范围为(0,logN);
结合种群方差、种群熵在进化过程中的特征,种群进化到第t代的交叉变异公式如式(7)、(8),种群方差、种群熵联合方程如公式(9):
式中:a=9.903438,Pcmax是最大交叉概率,Pcmin是最小交叉概率,Pmmax是最大变异概率,Pmmin是最小变异概率,logN是种群熵理论上的最大值,max{Vt}是种群进化到第t代时种群方差的最大值;
根据交叉概率Pc判断是否交叉,交叉的子代采用一致交叉算子,不交叉的子代从父代获取基因;根据变异概率Pm判断是否变异,变异采用一致变异操作,不发生变异的子代从父代获取基因;
步骤9,计算步骤8得到的初始化种群中最优个体和最差个体,采用精英保留策略更新初始化种群,之后返回步骤3;
步骤10,输出结果,结束程序。
交叉、变异概率随Φt的变化如图2所示,在进化初期Φt趋近于Φmax,此时种群多样性高,应该加大交叉概率,减小变异概率,促进个体之间交叉重组,增加全局搜索能力;进化中期,交叉、变异概率随着Φt按照Sigmoid曲线自适应变化;进化后期,Φt趋近于Φmin,种群集中,应该增大变异概率,增加种群的多样性,减小交叉概率,保护搜索到的较优解,使得算法保持稳定,提高算法收敛速度。
图3是测试函数f1在不同迭代次数下采用本发明算法、自适应遗传算法以及传统遗传算法的总收敛次数比较图;图4是测试函数f2在不同迭代次数下采用本发明算法、自适应遗传算法以及传统遗传算法的总收敛次数比较图,其中VHGA表示本发明算法,SGA表示传统遗传算法,IAGA表示自适应遗传算法。从图3和图4可以看出,本发明的算法曲线均在传统遗传算法和自适应遗传算法曲线之上,表明本发明的算法有更高的收敛概率,并且该算法的总收敛次数在每次统计中都具有较强的稳定性。
表1
通过表1可以总结出:对于测试函数f1,本发明的算法(VHGA)的总收敛次数比自适应遗传算法(IAGA)多95次,比传统遗传算法(SGA)多349次;对于测试函数f2,本发明的算法(VHGA)的总收敛次数比IAGA多57次,比SGA多184次,即本算法总收敛次数比SGA最多高45%,比IAGA最多高28%,所以本文设计的该算法在总收敛次数上优于其他两种算法。
Claims (5)
1.一种基于种群熵、种群方差改进的钢轨裂纹检测算法,其特征在于,具体按照以下步骤依次实施:
步骤1,设置算法的参数;
步骤2,初始化原种群,得到初始化种群;
步骤3,计算步骤2得到的初始化种群中个体的适应度值;
步骤4,根据步骤3得到的个体的适应度值,求出满足设定值的最优个体并加以保存;
步骤5,根据步骤4得到的最优个体的总数,与目标值作比较,达到要求则进行步骤10,反之依次重复进行步骤6、步骤7、步骤8、步骤9、步骤3、步骤4,直至达到要求;
步骤6,采用轮赌盘方法对步骤4得到的初始化种群进行选择操作;
步骤7,对步骤6得到的初始化种群进行交叉操作;
步骤8,对步骤7得到的初始化种群进行变异操作;
步骤9,计算步骤8得到的初始化种群中最优个体和最差个体,采用精英保留策略更新初始化种群,之后返回步骤3;
步骤10,输出结果,结束程序。
2.根据权利要求1所述的基于种群熵、种群方差改进的钢轨裂纹检测算法,其特征在于,在所述步骤1中,设置算法的参数为:原种群规模popsize=50,最小交叉概率Pcmin=0.6,最大交叉概率Pcmax=0.9,最小变异概率Pmmin=0.01,最大变异概率Pmmax=0.1;
3.根据权利要求2所述的基于种群熵、种群方差改进的钢轨裂纹检测算法,其特征在于,所述步骤2具体按照以下步骤依次实施:
步骤2.1,利用故障信号的时域、频域特征结合起来构建钢轨裂纹的故障判据如公式所示:
其中:xi为待优化系数,Ei为三层小波包分解第三层第i频段的能量函数,为输入信号样本数据,Amax为时域幅值函数,k为权重系数;
根据故障判据条件,得出下列故障判据函数如公式:
将m组样本数据求得的故障判据函数值与数据样本属性值相比较,如果故障判据函数值与数据样本属性值相同,说明判断正确,代价函数值增加;数据样本属性函数和代价函数分别描述为:
4.根据权利要求3所述的基于种群熵、种群方差改进的钢轨裂纹检测算法,其特征在于,在所述步骤3中:个体的适应度值根据公式(4)得到。
5.根据权利要求1所述的基于种群熵、种群方差改进的钢轨裂纹检测算法,其特征在于,所述步骤7和步骤8具体为:
设进化过程中第t代原种群P(t)中第i个体为N为原种群规模,个体的适应度为ft i,第t代种群个体的平均适应度为那么第t代种群的方差为:
式中:方差Vt反映了种群个体的离散程度,Vt的值大就说明个体比较分散,Vt的值小说明个体比较集中;
将信息熵引入种群多样性的判定中,用以反映种群内部进化状态;将实际问题中可行域等分为L个小区间,第i个小区间为Ai(i=1,…,L),在第t代进化中,落在中的个体数目记为定义第t代种群的熵为:
式中:表示每个小区间内的个体占种群规模N的概率大小,其中Ht的取值范围为(0,logN);
结合种群方差、种群熵在进化过程中的特征,种群进化到第t代的交叉变异公式如式(7)、(8),种群方差、种群熵联合方程如公式(9):
式中:a=9.903438,Pcmax是最大交叉概率,Pcmin是最小交叉概率,Pmmax是最大变异概率,Pmmin是最小变异概率,logN是种群熵理论上的最大值,max{Vt}是种群进化到第t代时种群方差的最大值;
根据交叉概率Pc判断是否交叉,交叉的子代采用一致交叉算子,不交叉的子代从父代获取基因;根据变异概率Pm判断是否变异,变异采用一致变异操作,不发生变异的子代从父代获取基因。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910058482.9A CN110009088A (zh) | 2019-01-22 | 2019-01-22 | 一种基于种群熵、种群方差改进的钢轨裂纹检测算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910058482.9A CN110009088A (zh) | 2019-01-22 | 2019-01-22 | 一种基于种群熵、种群方差改进的钢轨裂纹检测算法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110009088A true CN110009088A (zh) | 2019-07-12 |
Family
ID=67165457
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910058482.9A Pending CN110009088A (zh) | 2019-01-22 | 2019-01-22 | 一种基于种群熵、种群方差改进的钢轨裂纹检测算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110009088A (zh) |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108564592A (zh) * | 2018-03-05 | 2018-09-21 | 华侨大学 | 基于动态多种群集成差分进化算法的图像分割方法 |
-
2019
- 2019-01-22 CN CN201910058482.9A patent/CN110009088A/zh active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108564592A (zh) * | 2018-03-05 | 2018-09-21 | 华侨大学 | 基于动态多种群集成差分进化算法的图像分割方法 |
Non-Patent Citations (2)
Title |
---|
刘芳等: "基于种群多样性的自适应遗传算法优化仿真", 《计算机仿真》 * |
赵姣: "基于轮轨噪声的钢轨裂纹故障诊断算法研究", 《中国优秀硕士论文全文数据库 工程科技Ⅱ辑》 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103926526A (zh) | 一种基于改进的rbf神经网络的模拟电路故障诊断方法 | |
CN104573820A (zh) | 一种用于带约束条件的工程优化问题的遗传算法 | |
CN107277830A (zh) | 一种基于粒子群优化和变异算子的无线传感器网络节点部署方法 | |
CN107392919A (zh) | 基于自适应遗传算法的灰度阈值获取方法、图像分割方法 | |
CN106228979A (zh) | 一种公共场所异常声音特征提取及识别方法 | |
CN110222830A (zh) | 一种基于自适应遗传算法优化的深度前馈网络故障诊断方法 | |
CN115980826A (zh) | 一种基于加权的元启发式组合模型的岩爆烈度预测方法 | |
CN109492816A (zh) | 一种基于混合智能的煤与瓦斯突出动态预测方法 | |
CN117473424A (zh) | 基于随机森林的变压器故障诊断方法、系统、设备及介质 | |
CN109842614B (zh) | 基于数据挖掘的网络入侵检测方法 | |
CN116383630A (zh) | 基于改进灰狼算法的概率神经网络电弧故障检测方法 | |
CN110009088A (zh) | 一种基于种群熵、种群方差改进的钢轨裂纹检测算法 | |
CN106446376B (zh) | 一种考虑风险等级划分的海洋平台溜桩评估方法 | |
CN105512726B (zh) | 基于免疫遗传优化的可靠性分配方法及装置 | |
CN105654498A (zh) | 基于动态局部搜索和免疫克隆自动聚类的图像分割方法 | |
CN110175698A (zh) | 一种基于改进粒子群优化bp神经网络的水泵设备状态预测方法 | |
CN115996135A (zh) | 一种基于特征组合优化的工业互联网恶意行为实时检测方法 | |
CN1312636C (zh) | 形态学滤波器自动目标检测方法 | |
Yue et al. | Power system short-term load forecasting based on neural network with artificial immune algorithm | |
CN108446718A (zh) | 一种动态深度置信网络分析方法 | |
CN113361709A (zh) | 基于变异的深度神经网络模型修复方法 | |
An | Application of genetic algorithm based on f-ratio rule in signal feature selection | |
Chen et al. | Adaptive immune genetic algorithm for tire tread pattern pitch parameters optimization | |
Murata et al. | Many-objective optimization for knapsack problems using correlation-based weighted sum approach | |
Kowatari et al. | Analysis on population size and neighborhood recombination on many-objective optimization |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190712 |