CN108133258B - 一种混合全局优化方法 - Google Patents
一种混合全局优化方法 Download PDFInfo
- Publication number
- CN108133258B CN108133258B CN201611094736.5A CN201611094736A CN108133258B CN 108133258 B CN108133258 B CN 108133258B CN 201611094736 A CN201611094736 A CN 201611094736A CN 108133258 B CN108133258 B CN 108133258B
- Authority
- CN
- China
- Prior art keywords
- particle
- search
- particles
- algorithm
- chaotic
- 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.)
- Expired - Fee Related
Links
Images
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]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/0418—Architecture, e.g. interconnection topology using chaos or fractal principles
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Computing Systems (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- Evolutionary Computation (AREA)
- Life Sciences & Earth Sciences (AREA)
- Molecular Biology (AREA)
- Artificial Intelligence (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Health & Medical Sciences (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Feedback Control In General (AREA)
Abstract
本发明涉及一种混合全局优化方法,该方法首先对优化问题利用粒子群算法得到一组当前最优解;然后利用混沌搜索算法使得粒子跳出局部极值;最后,在粒子群算法每一代的迭代过程中引入序列二次规划算法加速局部最优点搜索,得到优化问题的全局最优解。本发明引入粒子群体适应度方差概念,结合混沌搜索和序列二次规划方法,当粒子群体适应度方差小于给定临界值时,粒子易陷入局部最优,对最优粒子进行混沌搜索,以使粒子跳出局部最优,并且根据粒子的进化速度和聚集度自适应的改变惯性权重,进而改变粒子运动状态进一步确保粒子不陷入局部最优;在粒子每一次迭代过程中,引入序列二次规划优化算法,加速粒子局部最优点的搜索,提高算法整体搜索效率。
Description
技术领域
本发明涉及一种优化方法,具体涉及一种混合全局优化方法。
背景技术
优化方法自20世纪60年代以来已取得很大进展和许多理论研究与应用成果。现有的优化方法主要可分为两大类:传统的确定性优化方法和智能优化方法。
以SQP方法为代表的确定性优化方法的实测效果非常好,但也有局限性。首先,只能求解光滑非线性优化问题;其次,获得的最优解只是局部最优解;最后只适合求解中小规模问题。受自然界自适应优化现象的启发,并且伴随着人工智能的发展,模拟生物进化过程、人类思维、生物智能行为等并以此开发求解复杂优化问题的智能优化算法越来越受到人们的关注。
智能优化算法可以分为两类:一类是模拟生物进化的算法,如遗传算法、进化规划等;另一类是基于群体智能的算法,如蚁群算法、粒子群算法等。尽管它们的表现形式和原理各不相同,但它们都有一些共同的特征,即都是群体搜索、随机搜索,具有并行性和全局性。
遗传算法(Genetic Algorithm,GA)是一类借鉴生物界自然选择和遗传机制的随机搜索算法,由美国Holland教授于1975年提出,用数码串来类比生物中的染色体,通过选择、交叉、变异等遗传算子来仿真生物基本进化过程,通过种群不断地“更新换代”获得问题最优解。粒子群算法(Partical Swarm Optimization,PSO)是由Eberhart博士和Kennedy博士于1995年提出的,它是通过模拟鸟群觅食行为而发展起来的一种基于群体协作的群智能随机搜索算法。与进化算法相比,粒子群算法保留了基于种群的全局搜索策略,但它又采用了相对简单的速度-位移模型,避免了复杂遗传算子操作。同时它特有的记忆功能使其可以动态的跟踪当前的搜索情况而适当调整其搜索策略。因此,粒子群算法是一种非常高效的并行搜索算法,自提出以来就受到学术界的广泛重视。
上述智能随机优化算法虽然不受问题种类约束,克服了传统优化方法的局限性,但是普遍存在搜索精度不高、易早熟、对参数依赖性高、搜索效率较低的缺陷。为克服上述缺陷,将几种智能优化方法相混合,对随机优化算法与确定性优化算法进行混合,对某种方法进行改进,从而可以产生有效的全局优化方法,这是优化理论方法中出现的新的研究方向。
基于遗传算法和线性反演方法的混合优化算法思想是:首先在传统遗传算法基础上,对交叉、变异后所得结果的处理以及进化过程的整体分布方面进行改进,并将改进的遗传算法与广义线性反演方法相结合,以改进的遗传算法为基础,在进化的每一代种群中选择目标函数最小的个体,进行一定次数的线性化迭代。混合优化算法克服了线性化方法依赖于初始值和遗传算法局部搜索能力差的缺陷。
基于遗传算法和序列二次规划算法的混合优化算法思想是:首先利用遗传算法进行初始化种群,通过选择、交叉和变异操作进行多次迭代,求得问题初始最优解,接着运用序列二次规划方法以初始最优解为初始搜索点进行迭代搜索,求出最终最优解。该算法能够发挥遗传算法全局寻优和序列二次规划算法局部寻优能力强的优点,有效提高算法搜索精度和搜索效率。
基于粒子群和共轭梯度的混合优化方法思想是:首先利用粒子群算法进行初始化种群,通过粒子群速度和位置多次迭代求出粒子群算法的初始最优解,此时,以初始最优解的最优位置为初始搜索点,以最优位置处的负梯度方向为搜索方向,利用共轭梯度法确定下一个最优位置并继续进行基于共轭梯度法的迭代搜索,直到满足共轭梯度法的终止条件。
上述三种混合优化方法虽然提高了搜索精度和搜索效率,但是仍存在难以跳出局部最优的问题。
发明内容
本发明的目的在于提供一种将粒子群算法、混沌搜索算法和序列二次规划算法相结合的混合全局优化方法,该混合算法具有较高的计算精度和收敛速度,可解决机械优化设计领域的相关问题。
实现本发明目的的技术方案为:一种混合全局优化方法,包括以下步骤:
步骤1,对优化问题利用粒子群算法得到一组当前最优解;
步骤2,利用混沌搜索算法使得粒子跳出局部极值;
步骤3,在粒子群算法每一代的迭代过程中引入序列二次规划算法加速局部最优点搜索,得到优化问题的全局最优解。
进一步的,步骤1具体为:
步骤1-1,初始化下列参数:第一学习因子c1,第二学习因子c2,粒子群算法最大迭代次数MaxDT,粒子数目N,搜索空间维数D,粒子最大飞行速度vmax,最小飞行速度vmin,变量上下限[lb ub],适应度方差临界值[σ2],混沌搜索调节参数β,混沌搜索最大迭代次数M,进化速度因子pspeed,聚集度因子ptogether,进化速度因子的惯性权重调节参数wh和聚集度因子的惯性权重调节参数ws;
步骤1-2,令t=0,t为粒子群算法迭代次数,0≤t≤MaxDT,初始化粒子群位置xi=(xi,1,xi,2,…,xi,D)和速度vi=(vi,1,vi,2,…,vi,D),1≤i≤N;
步骤1-3,计算粒子适应度值fi=f(xi),并初始化粒子群个体最优粒子pi和全局最优粒子pg,计算个体最优值fi,best=f(pi)和全局最优值fbest=f(pg);
步骤1-4,令t=1,按式(1)计算惯性权重w;
w=1-pspeed×wh+ptogether×ws (1)
式中,1≤j≤D;rand1和rand2均为(0,1)的随机数;
步骤1-6,计算粒子适应度值fi=f(xi),并更新pi和pg。
进一步的,步骤2具体为:
步骤2-2,如果σ2<[σ2],则按式(6)计算混沌搜索概率p,否则转步骤2-1;
rand3为(0,1)的随机数,如果rand3<p,则对粒子群中最优粒子pg进行混沌优化搜索,否则转步骤3。
进一步的,混沌优化搜索的具体过程为:
第二步,将混沌变量按式(7)线性映射到优化变量取值区间[lb ub]
rxgj为混沌搜索扰动变量;
第三步,令j=1,对最优解的每一维变量按式(8)进行混沌搜索
第五步,如果k<M,k=k+1,更新混沌变量,并返回第二步;
进一步的,步骤3具体为:
步骤3-1,令最优粒子pg为初始点,运行序列二次规划算法,更新pg和fbest;
步骤3-2,如果t<MaxDT,则t=t+1,否则返回步骤1-5;
步骤3-3,输出pg和fbest。
进一步的,所述步骤1-3中的适应度值fi=f(xi)为待求问题的目标函数值。
进一步的,所述步骤1-3中的个体最优粒子pi和全局最优粒子pg为当前迭代的个体最优解和全局最优解。
与现有技术相比,本发明的显著优点为:本发明引入粒子群体适应度方差概念,并结合混沌搜索和序列二次规划方法,当粒子群体适应度方差小于给定临界值时,说明粒子聚集度较高,此时粒子易陷入局部最优,应该对最优粒子进行混沌搜索,以使粒子跳出局部最优,并且根据粒子的进化速度和聚集度自适应的改变惯性权重,进而改变粒子运动状态进一步确保粒子不陷入局部最优;并且在粒子每一次迭代过程中,引入序列二次规划优化算法,加速粒子局部最优点的搜索,提高算法整体搜索效率。
附图说明
图1为本发明混合全局优化方法的计算流程框图。
图2为基本粒子群算法计算流程框图。
图3为混沌优化算法计算流程框图。
图4为序列二次规划算法计算流程框图。
图5为悬臂梁结构示意图。
具体实施方式
粒子群算法作为经典的群智能优化算法,其基本流程如图2所示,单一的粒子群算法存在前期容易陷入局部最优,而后期搜索缓慢的不足。混沌搜索具有随机性、遍历性和规律性的特点,加入时间足够长,就能够不重复的遍历解空间,其基本流程如图3所示。序列二次规划算法具有局部寻优能力强局部收敛速度快的优点,其计算流程图如图4所示。
在基本粒子群优化算法中,所有粒子都有一个被目标函数决定的适应度值和决定它们位置变化和飞行方向的速度。基本粒子群进化算法的进化方程可以看成是由认知和社会两部分组成,认知代表各粒子本身的思考,粒子跟随粒子本身的最优解飞行,社会代表粒子间的信息共享,粒子跟随全局最优粒子飞行。
假设搜索空间为D维,粒子数目为N,xi=(xi,1,xi,2,…,xi,D)是粒子i的当前位置,vi=(vi,1,vi,2,…,vi,D)是粒子当前飞行速度,则粒子群中的粒子按式(1)和式(2)改变速度和位置:
式中,t为迭代次数;w为惯性权重;为粒子i历史经过的最好位置;为粒子群的全局最好位置;c1、c2分别为认知学习因子和社会学习因子,统称为学习因子;rand1和rand2是在(0,1)上的相互独立的随机数。
混沌现象是一种非周期的具有渐进的自相似有序性的现象。因其具有随机性、遍历性、规律性的独特性质,常作为优化设计领域搜索过程中避免陷入局部极值的一种全局优化机制。Logistic映射一般形式如式(3):
xgj为混沌变量的第j维分量;
分析和研究表明,当控制参数μ=4时,Logistic映射处于完全混沌状态,具有很好遍历性。许多关于混沌的试验研究大都是基于这种完全混沌状态进行的。本发明将在粒子群适应度方差小于给定值,即粒子聚集度较高,可能陷入局部极值的时候,以一定概率p对当前最优粒子中的每一维变量进行混沌搜索。
Logistic混沌搜索算法实现步骤为:
第二步,将混沌变量按式(4)线性映射到优化变量取值区间[lb ub]
rxgj为混沌搜索扰动变量;
第三步,令j=1,对最优解的每一维变量按式(5)进行混沌搜索
第五步,如果k<M,k=k+1,按式(3)更新混沌变量,并返回第二步;
序列二次规划方法由Wilson于1963年提出,后经韩世平以及Powell改进后,在非线性局部寻优的问题上显现出强大的求解能力。其核心思想是利用一系列二次规划子问题逼近原问题的最优解。原问题描述如式(6):
二次规划子问题描述如式(7):
式中,W(x,λ,μ)为原目标函数f(x)的Lagrange函数在迭代点处的Hesse矩阵;g(x)为原函数的梯度向量;JE(x)、JI(x)分别为等式约束cE(x)和不等式约束cI(x)的Jacobi矩阵。式(7)的解δx可以作为原问题下一步迭代的搜索方向;l为序列二次规划方法当前迭代次数。
二次规划算法具体实施步骤为:
步骤1,已知待求问题的目标函数f(x)以及约束函数cE(x)和cI(x),选择初始点z(0)=(x(0),λ(0),μ(0)),并设定精度要求tol,令l=0;
步骤2,计算点x(l)处的目标函数的梯度值g(x(l)),约束函数值cE(x(l))和cI(x(l)),并根据式(8)和式(9)计算约束函数的Jacobi矩阵值JE(x(l))和JI(x(l)),根据式(10)计算拉格朗日函数的Hesse矩阵值W(x(l),λ(l),μ(l)):
m1为等式约束个数,m2为不等式约束个数;
本发明引入粒子群体适应度方差概念,并结合混沌搜索和序列二次规划方法,当粒子群体适应度方差小于给定临界值时,说明粒子聚集度较高,此时粒子易陷入局部最优,应该对最优粒子进行混沌搜索,以使粒子跳出局部最优,并且根据粒子的进化速度和聚集度自适应的改变惯性权重,进而改变粒子运动状态进一步确保粒子不陷入局部最优。并且在粒子每一次迭代过程中,引入序列二次规划优化算法,加速粒子局部最优点的搜索,提高算法整体搜索效率。
本发明在上述技术方案的基础上,从群体适应度方差、自适应惯性权重调整等方面进行改进,具体如下:
一、群体适应度方差
当粒子处于早熟收敛状态时,粒子会出现“靠拢”现象,而粒子聚集位置取决于其适应度值。因此,整个粒子群的状态就可以用粒子群体适应度方差来表示,适应度方差越小,表明粒子越处于收敛状态。假定fi为第i个粒子的适应度,favg为粒子平均适应度,σ2为群体适应度方差,则
式中:f为归一化因子,取值采用如下公式:
二、自适应惯性权重
基本PSO算法中,惯性权重w一般为线性递减,针对非线性度较高的问题,算法易陷入局部最优。为了克服上述弊端,应该根据粒子群的进化速度和聚集度的变化自适应的改变惯性权重以动态调整粒子运动状态,避免粒子陷入局部最优。
在优化过程中,粒子的进化速度可以由全局最优粒子的变化来表示,而粒子的聚集度可以由当前粒子群体中最优粒子适应度和平均适应度的相关表达式来表示。假设粒子群全局最优值为fbest,则进化速度因子pspeed和聚集度因子ptogether可以表示为:
当粒子进化速度减小时,应该使w的值减小,以便粒子在小空间搜索;当粒子聚集度较高时,说明粒子易陷入局部最优,应该使w的值增大,以便粒子增加搜索空间,避免陷入局部最优。综合分析,惯性权重w应该随着进化速度降低而减小,随着聚集度增加而增大,因此,w可以表示为:
w=1-pspeed×wh+ptogether×ws (14)
式中,wh取值区间为[0.4 0.6],ws取值区间为[0.05 0.20]。
如图1所示,本发明的混合全局优化方法具体步骤为:
步骤1,对优化问题利用粒子群算法得到一组当前最优解;
步骤1-1,参数初始化:第一学习因子c1,第二学习因子c2,粒子群算法最大迭代次数MaxDT,粒子数目N,搜索空间维数D,粒子最大飞行速度vmax,最小飞行速度vmin,变量上下限[lb ub],适应度方差临界值[σ2],混沌搜索调节参数β,混沌搜索最大迭代次数M,进化速度因子pspeed,聚集度因子ptogether,进化速度因子的惯性权重调节参数wh和聚集度因子的惯性权重调节参数ws;
步骤1-2,令t=0,初始化粒子群位置xi=(xi,1,xi,2,…,xi,D)和速度vi=(vi,1,vi,2,…,vi,D);
步骤1-3,计算粒子适应度值fi=f(xi),并初始化粒子群个体最优粒子pi和全局最优粒子pg,个体最优值fi,best和全局最优值fbest;
步骤1-4,令t=1,按式(14)计算惯性权重w;
步骤1-6,计算粒子适应度值fi=f(xi),并更新pi和pg;
步骤2,利用混沌搜索算法使得粒子跳出局部极值;
步骤2-1,计算粒子群中粒子平均适应度favg,并按式(12)计算归一化因子f,按式(11)计算粒子群体适应度方差σ2;
步骤2-2,如果σ2<[σ2],则按式(15)计算混沌搜索概率p,否则转步骤7;如果rand<p,则对粒子群中最优粒子pg进行混沌优化搜索,否则转步骤3:
步骤3,在粒子群算法每一代的迭代过程中引入序列二次规划算法加速局部最优点搜索,得到优化问题的全局最优解;具体为:
步骤3-1,对最优粒子pg为初始点,运行序列二次规划算法,更新pg和fbest;
步骤3-2,如果t<MaxDT,则t=t+1,否则返回步骤1-5;
步骤3-3,输出pg和fbest;
下面结合具体实施例对本发明作进一步说明。
实施例
为了验证混合全局优化算法的有效性,采用混沌粒子群、序列二次规划和混合全局优化算法对三个经典函数进行30次寻优,并对求解出的函数最优值进行比较,分析算法计算精度的大小。
算法参数设置为:学习因子c1=c2=1.49,粒子数目N=40,粒子群优化最大迭代次数MaxDT=100,粒子所允许的飞行速度范围为[vmin vmax]=[-10 10],混沌搜索最大迭代次数M=100,混沌搜索调节参数β=0.1,其中适应度方差临界值[σ2]视具体情况再设置。三个基准测试函数信息如表1所示:
表1
使用混沌粒子群、序列二次规划和混合全局优化三种算法分别对三个经典函数寻优30次,其中序列二次规划算法的初始搜索点为在变量范围内随机生成,得到最好全局最优值和平局全局最优值,结果如表2所示:
表2
从表2中可以看出:混沌粒子群算法基本能够求解出全局最优解,但是计算精度不高,当变量维数较多时,算法稳定性不好;序列二次规划算法由于对初始搜索点较敏感,很难搜索到全局最优解;而混合全局优化算法在变量维数较少时,求解结果的计算精度与混沌粒子群算法差异不明显,当变量维数较多时,混合全局优化算法在计算精度上相比混沌粒子群算法有了很大提高,并且求解结果中的最好全局最优值和平均全局最优值差异非常小,说明算法稳定性好。
通过一个机械优化实例来说明混合全局优化方法的具体实现过程:
悬臂梁优化问题:悬臂梁是指梁的一端为不产生轴向、垂直位移和转动的固定支座,另一端为自由端。在实际工程分析中,大部分实际工程受力部件都可以简化为悬臂梁,如图5所示。其设计变量为h(x1)、l(x2)、t(x3)、b(x4),目标函数为最小化总费用f(x)。
其数学模型为如式(16):
本实施例的混合全局优化方法具体步骤为:
步骤1,参数初始化:学习因子c1=c2=1.49,粒子群算法最大迭代次数MaxDT=100,粒子数目N=40,搜索空间维数D=4,粒子最大飞行速度vmax=10和最小飞行速度vmin=-10,适应度方差临界值[σ2]=3,混沌搜索调节参数β=0.1,混沌搜索最大迭代次数M=100,初始进化速度因子pspeed=0.5和聚集度因子ptogether=5,惯性权重调节参数wh和ws;
步骤2,令t=0,初始化粒子群位置xi=(xi,1,xi,2,…,xi,D)和速度vi=(vi,1,vi,2,…,vi,D);
步骤3,计算粒子适应度值fi=f(xi),并初始化粒子群个体最优粒子pi和全局最优粒子pg,个体最优值fi,best和全局最优值fbest;
步骤4,令t=1,计算惯性权重w=1-pspeed×wh+ptogether×ws;
步骤5,更新粒子的速度
步骤6,计算粒子适应度值fi=f(xi),并更新pi和pg;
步骤9,对最优粒子pg为初始点,运行SQP算法,更新pg和fbest;
步骤10,如果t<MaxDT,则t=t+1,否则返回步骤5;
步骤11,输出pg和fbest;
优化结果如表3:
表3
本发明以基本粒子群算法为主体,并结合混沌搜索和序列二次规划方法,当粒子群体适应度方差小于给定临界值时,说明粒子聚集度较高,此时粒子易陷入局部最优,应该对最优粒子进行混沌搜索,以使粒子跳出局部最优,并且根据粒子的进化速度和聚集度自适应的改变惯性权重,进而改变粒子运动状态进一步确保粒子不陷入局部最优。并且在粒子每一次迭代过程中,引入序列二次规划优化算法,加速粒子局部最优点的搜索,提高算法整体搜索效率。
Claims (4)
1.一种混合全局优化方法,其特征在于,包括以下步骤:
步骤1,对悬臂梁优化问题利用粒子群算法得到一组当前最优解;悬臂梁优化问题的目标函数为最小化总费用f(x),其数学模型为如式(16):
步骤1-1,初始化下列参数:第一学习因子c1,第二学习因子c2,粒子群算法最大迭代次数MaxDT,粒子数目N,搜索空间维数D,粒子最大飞行速度vmax,最小飞行速度vmin,变量上下限[lb ub],适应度方差临界值[σ2],混沌搜索调节参数β,混沌搜索最大迭代次数M,进化速度因子pspeed,聚集度因子ptogether,进化速度因子的惯性权重调节参数wh和聚集度因子的惯性权重调节参数ws;
步骤1-2,令t=0,t为粒子群算法迭代次数,0≤t≤MaxDT,初始化粒子群位置xi=(xi,1,xi,2,…,xi,D)和速度vi=(vi,1,vi,2,…,vi,D),1≤i≤N;
步骤1-3,计算粒子适应度值fi=f(xi),并初始化粒子群个体最优粒子pi和全局最优粒子pg,计算个体最优值fi,best=f(pi)和全局最优值fbest=f(pg);
步骤1-4,令t=1,按式(1)计算惯性权重w;
w=1-pspeed×wh+ptogether×ws (1)
进化速度因子pspeed和聚集度因子ptogether表示为:
式中,1≤j≤D;rand1和rand2均为(0,1)的随机数;
步骤1-6,计算粒子适应度值fi=f(xi),并更新pi和pg;
步骤2,利用混沌搜索算法使得粒子跳出局部极值;具体为:
步骤2-2,如果σ2<[σ2],则按式(6)计算混沌搜索概率p,否则转步骤3;
rand3为(0,1)的随机数,如果rand3<p,则对粒子群中最优粒子pg进行混沌优化搜索,否则转步骤3;混沌优化搜索的具体过程为:
第二步,将混沌变量按式(7)线性映射到优化变量取值区间[lb ub]
rxgj为混沌搜索扰动变量;
第三步,令j=1,对最优解的每一维变量按式(8)进行混沌搜索
第五步,如果k<M,k=k+1,更新混沌变量,并返回第二步;
步骤3,在粒子群算法每一代的迭代过程中引入序列二次规划算法加速局部最优点搜索,得到悬臂梁优化问题的全局最优解。
2.根据权利要求1所述的混合全局优化方法,其特征在于,步骤3具体为:
步骤3-1,令最优粒子pg为初始点,运行序列二次规划算法,更新pg和fbest;
步骤3-2,如果t<MaxDT,则t=t+1,返回步骤1-4,否则执行步骤3-3;
步骤3-3,输出pg和fbest。
3.根据权利要求1所述的混合全局优化方法,其特征在于,所述步骤1-3中的适应度值fi=f(xi)为待求问题的目标函数值。
4.根据权利要求1所述的混合全局优化方法,其特征在于,所述步骤1-3中的个体最优粒子pi和全局最优粒子pg为当前迭代的个体最优解和全局最优解。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611094736.5A CN108133258B (zh) | 2016-12-01 | 2016-12-01 | 一种混合全局优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611094736.5A CN108133258B (zh) | 2016-12-01 | 2016-12-01 | 一种混合全局优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108133258A CN108133258A (zh) | 2018-06-08 |
CN108133258B true CN108133258B (zh) | 2020-11-13 |
Family
ID=62387741
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611094736.5A Expired - Fee Related CN108133258B (zh) | 2016-12-01 | 2016-12-01 | 一种混合全局优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108133258B (zh) |
Families Citing this family (28)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109165775A (zh) * | 2018-08-07 | 2019-01-08 | 上海电机学院 | 一种基于二进制粒子群的能源互联网的优化调度方法 |
CN109325267A (zh) * | 2018-08-30 | 2019-02-12 | 广西民族大学 | 一种高维化工数据共轭粒子群算法处理方法 |
CN109657274B (zh) * | 2018-11-16 | 2023-08-29 | 广东省建筑设计研究院 | 建筑结构中基于粒子群优化算法的弦支穹顶索力优化方法 |
CN110145432B (zh) * | 2019-03-29 | 2021-07-09 | 广东工业大学 | 一种基于傅氏分析和改进灰狼算法的波浪发电装置功率控制方法 |
CN110110380B (zh) * | 2019-04-11 | 2023-07-04 | 上海电力学院 | 一种压电执行器迟滞非线性建模方法及应用 |
CN110322075A (zh) * | 2019-07-10 | 2019-10-11 | 上饶市中科院云计算中心大数据研究院 | 一种基于混合优化rbf神经网络的景区客流量预测方法及系统 |
CN110516958B (zh) * | 2019-08-23 | 2022-03-18 | 河海大学常州校区 | 一种面对制造过程的资源调度方法 |
CN110674915B (zh) * | 2019-09-18 | 2022-11-25 | 东北大学 | 基于改进粒子群算法的不规则管道缺陷的反演方法 |
CN111144623B (zh) * | 2019-12-09 | 2023-03-31 | 云南电网有限责任公司临沧供电局 | 一种基于自适应学习因子粒子群的定值调优方法 |
CN111162886A (zh) * | 2019-12-12 | 2020-05-15 | 聂阳 | 数字调幅广播信道估计中导频图案的分配优化方法 |
CN111178599B (zh) * | 2019-12-16 | 2022-05-31 | 浙江工业大学 | 一种基于改进型粒子群算法的订单排产方法 |
CN111080035A (zh) * | 2019-12-31 | 2020-04-28 | 芜湖哈特机器人产业技术研究院有限公司 | 基于改进量子粒子群优化算法的全局路径规划方法 |
CN111222799B (zh) * | 2020-01-13 | 2023-10-27 | 北京工业大学 | 一种基于改进粒子群算法的装配序列规划方法 |
CN111400785A (zh) * | 2020-01-17 | 2020-07-10 | 河海大学 | 一种边坡稳定和安全级别的综合评判方法 |
CN111539117B (zh) * | 2020-04-28 | 2023-06-13 | 南阳师范学院 | 一种基于化学平衡模型的非线性方程系统求解方法及装置 |
CN111695275B (zh) * | 2020-05-11 | 2022-11-11 | 东南大学 | 一种提高超导磁储能线圈优化计算速度的方法 |
CN111709068B (zh) * | 2020-05-14 | 2024-01-12 | 西安理工大学 | 一种用于建筑结构优化设计的一维搜索算法 |
CN111767671B (zh) * | 2020-06-19 | 2022-12-09 | 西安交通大学 | 一种适用于多面转子的不平衡参数辨识方法 |
CN111694396A (zh) * | 2020-07-04 | 2020-09-22 | 湘潭大学 | 一种基于分子运动轨迹搜索算法的mppt控制 |
CN112365115A (zh) * | 2020-08-26 | 2021-02-12 | 天津大学 | 一种电网信息能源系统稳定性评估方法 |
CN112000115B (zh) * | 2020-09-15 | 2022-11-01 | 湖南文理学院 | 一种自主式水下机器人的路径规划方法及相关设备 |
CN112861426A (zh) * | 2021-01-14 | 2021-05-28 | 西北工业大学 | 基于改进粒子群算法的航空发动机加速过程最优控制方法 |
CN113222096B (zh) * | 2021-04-30 | 2022-11-11 | 桂林理工大学 | 一种面向云计算任务调度的改进粒子群算法 |
CN113305834B (zh) * | 2021-05-18 | 2022-06-21 | 广东工业大学 | 一种二自由度机械臂最优控制问题的全局模式搜索方法 |
CN113435722B (zh) * | 2021-06-17 | 2024-01-23 | 上海海事大学 | U型自动化码头多设备混合调度方法及电子设备 |
CN116152316B (zh) * | 2023-04-17 | 2023-07-07 | 山东省工业技术研究院 | 一种基于自适应参数粒子群算法的图像配准方法 |
CN116341390B (zh) * | 2023-05-11 | 2023-11-17 | 西安现代控制技术研究所 | 一种全局搜索快速收敛多约束弹道优化方法 |
CN117149662A (zh) * | 2023-10-27 | 2023-12-01 | 深圳云塔信息技术有限公司 | 一种基于惯性权重混沌粒子群优化技术的测试方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102208046A (zh) * | 2011-05-19 | 2011-10-05 | 山东大学 | 基于二进制粒子群优化算法的混沌时间序列预测方法 |
CN103217900A (zh) * | 2013-02-06 | 2013-07-24 | 浙江工业大学 | 基于实时电价的中压微网混沌pso最优潮流实现方法 |
CN104376360A (zh) * | 2014-11-11 | 2015-02-25 | 浪潮电子信息产业股份有限公司 | 一种新的混沌粒子群优化算法 |
-
2016
- 2016-12-01 CN CN201611094736.5A patent/CN108133258B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102208046A (zh) * | 2011-05-19 | 2011-10-05 | 山东大学 | 基于二进制粒子群优化算法的混沌时间序列预测方法 |
CN103217900A (zh) * | 2013-02-06 | 2013-07-24 | 浙江工业大学 | 基于实时电价的中压微网混沌pso最优潮流实现方法 |
CN104376360A (zh) * | 2014-11-11 | 2015-02-25 | 浪潮电子信息产业股份有限公司 | 一种新的混沌粒子群优化算法 |
Non-Patent Citations (4)
Title |
---|
Constrained optimization using the quasi-chaotic optimization method with the exact penalty function and the sequential quadratic programming;Takashi Okamoto etal.;《2011 IEEE International Conference on Systems, Man, and Cybernetics》;20111121;第1765-1770页 * |
基于SQP 局部搜索的混沌粒子群优化算法;徐文星 等;《控制与决策》;20120430;第27卷(第4期);第557-561页 * |
基于优化神经网络和DGA的变压器故障诊断;苗长新;《高压电学》;20161116;第52卷(第11期);第163-168页 * |
基于混沌搜索的粒子群优化算法;杨俊杰 等;《计算机工程与应用》;20051231;第69-71页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108133258A (zh) | 2018-06-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108133258B (zh) | 一种混合全局优化方法 | |
Valdez et al. | Comparative study of the use of fuzzy logic in improving particle swarm optimization variants for mathematical functions using co-evolution | |
Cheng et al. | An innovative hybrid multi-objective particle swarm optimization with or without constraints handling | |
Grimaldi et al. | PSO as an effective learning algorithm for neural network applications | |
Xiaojing et al. | Center particle swarm optimization algorithm | |
Tang et al. | A hybrid improved whale optimization algorithm | |
Liu et al. | A unified framework for population-based metaheuristics | |
CN110110380B (zh) | 一种压电执行器迟滞非线性建模方法及应用 | |
CN113255873A (zh) | 一种聚类天牛群优化方法、系统、计算机设备和存储介质 | |
Panda et al. | Reduction of linear time-invariant systems using Routh-approximation and PSO | |
Li et al. | An improved binary quantum-behaved particle swarm optimization algorithm for knapsack problems | |
Wang et al. | Electricity load forecasting based on support vector machines and simulated annealing particle swarm optimization algorithm | |
Xia et al. | Research on improved chaotic particle optimization algorithm based on complex function | |
Song et al. | Prediction for chaotic time series of optimized BP neural network based on modified PSO | |
Li et al. | Prediction for short-term traffic flow based on optimized wavelet neural network model | |
Wei et al. | Improved grey wolf optimization based on the Quantum-behaved mechanism | |
Guoqiang et al. | Study of RBF neural network based on PSO algorithm in nonlinear system identification | |
Li et al. | A modified particle swarm optimization with elite archive for typical multi-objective problems | |
CN109657778B (zh) | 一种基于改进的多种群全局最优的自适应鸽群优化方法 | |
Hao et al. | An Improved Glowworm Swarm Optimization Algorithm | |
Yin et al. | An improved artificial bee colony algorithm for staged search | |
Amini et al. | Interval search with quadratic interpolation and stable deviation quantum-behaved particle swarm optimization (IQS-QPSO) | |
Jiao | Simulation Research on the Palm Mechanism of Volleyball Robot Based on Artificial Intelligence and Ant Colony Optimization Algorithm | |
Niu et al. | Quantum-behaved particle swarm optimization algorithm based on dynamic dual-population joint-search mechanism | |
Singh | Test Suite Minimization using Evolutionary Optimization Algorithms |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20201113 Termination date: 20211201 |