CN106228234A - 一种基于梯度下降法的多目标粒子群优化方法 - Google Patents
一种基于梯度下降法的多目标粒子群优化方法 Download PDFInfo
- Publication number
- CN106228234A CN106228234A CN201610573200.5A CN201610573200A CN106228234A CN 106228234 A CN106228234 A CN 106228234A CN 201610573200 A CN201610573200 A CN 201610573200A CN 106228234 A CN106228234 A CN 106228234A
- Authority
- CN
- China
- Prior art keywords
- particle
- population
- gradient descent
- swarm optimization
- gradient
- 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
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)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及人工智能研究领域,尤其涉及一种基于梯度下降法的多目标粒子群优化方法,本方法将粒子群优化算法扩展到多目标优化,同时结合梯度下降法,改善多目标粒子群算法的优化性能,解决了非连续函数梯度的计算,加快了算法在计算梯度方面的速度,降低了算法的复杂度。
Description
技术领域
本发明涉及人工智能研究领域,尤其涉及一种基于梯度下降法的多目标粒子群优化方法。
背景技术
优化问题存在于各个方面,是指利用优化技术得到待求解问题的最优解。在理论研究与实际优化问题中,往往存在多个需要同时考虑的目标,即多目标优化问题,此类问题由于各个目标间存在相互制约的关系,很难用常规的数学方法(如线性规划)求解。
目前,对于多目标问题的处理方式大致可以分为两类:第一类是利用加权的思想,将多个目标转化为单个目标,从而利用单目标求解技术进行求解;第二类方法是对多个目标同时优化,利用Pareto最优的思想,构建多目标优化问题的Pareto最优解集,从而给决策者提供参考信息。
粒子群优化算法(particle swarm optimization,PSO)是一种模拟鸟群觅食的群智能优化算法,由Kennedy等提出,通过追随解空间中两个最优值(个体历史最优pBest,全局最优gBest)进行寻优。通过如下两个公式进行更新:
vi(k+1)=ω·vi(k)+c1·γ1·[Pi-xi(k)]+c2·γ2·[Pg-xi(k)];
xi(k+1)=vi(k+1)+xi(k);
i为粒子索引;k为离散时间索引;vi(k)为第i个粒子在第k代的速度;xi(k)为第i个粒子在第k代的位置;ω为惯性权重;c1为认知加 速常数;c2为社会加速常数;γ1γ2为[0,1]区间内的随机数;Pi为第i个粒子的个体最位置(个体最优);Pg为粒子群群体历史最优位置(全局最优)。
粒子群优化算法是目前处理多目标优化问题的一个有效方向,利用随机产生的初始种群,通过更新迭代求解,已被成功应用于多目标优化领域,在工程优化、调度等方向有广泛应用。但是从单目标扩展到多目标问题时,由于进化机制产生了变化,得到的不是单个最优解,而是一组解集,对于解集的更新与维护,pBest,gBest等选取与单目标时有许多差异;目前对多目标优化的研究主要针对其他进化算法,而对多目标粒子群优化算法研究较少。
通常对于单目标优化问题,梯度信息对于寻优的改善作用是有目共睹的,产生了最速下降法、牛顿法、共轭梯度法等求解方法。梯度下降法是利用负梯度方向来决定每次迭代的新的搜索方向,使得每次迭代都能朝着函数最小(对于最小优化问题)的方向前进。梯度下降法是二范数下的最速下降法。梯度下降法的一种简单形式是:
其中α称为学习速率,是一个较小的常数,是x(k)的梯度。单纯的梯度下降法容易导致算法收敛于最近的极值点,而粒子群算法由于它的随机的特性从而容易跳出局部极值点,将两者的特性结合起来对粒子群的寻优求解会有很大帮助。
发明内容
本发明为克服上述的不足之处,目的在于提供一种基于梯度下降 法的多目标粒子群优化方法,本方法将粒子群优化算法扩展到多目标优化,同时结合梯度下降法,改善多目标粒子群算法的优化性能,提升优化效率。
本发明是通过以下技术方案达到上述目的:一种基于梯度下降法的多目标粒子群优化方法,包括如下步骤:
(1)对粒子群种群进行初始化;在目标空间中随机化粒子群,设置每个粒子的个体极值与全局极值为初始位置;设置外部档案集为空;初始化迭代次数t=0;
(2)计算粒子群内各粒子的目标函数值,并判断目标函数值是否满足迭代终止准则,若满足最大迭代次数,则完成优化并输出结果;否则更新迭代次数t=t+1,执行步骤(3);
(3)判断粒子的支配关系,将非支配解放入外部档案集;
(4)在粒子群内选取最优粒子,更新该粒子的个体极值与全局极值;
(5)计算粒子群内所有粒子在各目标函数上的适应值,将新的非支配解更新放入外部档案集;
(6)利用梯度策略更新每个粒子的位置和速度后跳转执行步骤(2)。
作为优选,所述种群初始化包括设置种群大小N、目标个数M、最大迭代次数T。
作为优选,所述步骤(4)选取最优粒子,方法如下:
(i)对于历史最优pBest的选取,利用粒子之间的支配关系确定,从当前粒子和历史最优中选择非支配解为pBest;若粒子互不支配,则随机选取一个为pBest,利用sigma方法选取gBest;
(ii)对于全局最优gBest的选取,利用sigma方法进行选取。
作为优选,所述利用sigma方法选取gBest的方法步骤如下:
1)分别计算外部档案集中粒子的sigma(σ)值,σ表达式如下:
2)从种群中随机选取一个粒子,计算其σi值;
3)计算该粒子σi值与外部档案集的粒子σ值之间的欧氏距离,选取欧式距离最小的外部档案集的粒子作为该粒子的gBest。
作为优选,所述利用梯度策略更新粒子的位置和速度具体为利用梯度下降法结合粒子群算法更新粒子的位置与速度,方法如下:利用个体最优位置pBest与全局最优位置gBest指导粒子的飞行方向,在每次迭代更新寻优方向时,通过利用每个目标函数的梯度信息指导粒子的飞行方向,结合粒子本身的pBest、gBest得到新的粒子群速度。
作为优选,所述粒子速度的更新公式如下:
vi(k+1)=ω·vi(k)+c1·γ1·[Pi-xi(k)]+c2·γ2·[Pg-xi(k)]+c3·gi(f[x(k)])
其中,i为粒子索引;k为离散时间索引;vi(k)为第i个粒子在第k代的速度;xi(k)为第i个粒子在第k代的位置;ω为惯性权重;c1为认知加速常数;c2为社会加速常数;c3是加速常数;γ1γ2为[0,1]区间内的随机数;Pi为第i个粒子的个体最位置;Pg为粒子群群体历史最优位置;f[xi(k)]为目标函数的偏导数。
作为优选,所述f[xi(k)]由如下公式得到:
通过导数定义对上式进行变换得到f[xi(k)]的计算公式,如下所示:
作为优选,所述粒子位置的更新公式如下:
xi(k+1)=vi(k+1)+xi(k)
其中,i为粒子索引;k为离散时间索引;vi(k)为第i个粒子在第k代的速度;xi(k)为第i个粒子在第k代的位置。
本发明的有益效果在于:本方法将梯度下降法与粒子群算法相契合,解决了非连续函数梯度的计算,加快了算法在计算梯度方面的速度,降低了算法的复杂度。
附图说明
图1是本发明实施例优化方法的流程示意图。
具体实施方式
下面结合具体实施例对本发明进行进一步描述,但本发明的保护范围并不仅限于此:
实施例:如图1所示,本发明提出了基于梯度下降法的多目标粒子群优化算法,利用粒子群算法更新迭代的原理,结合梯度下降法,巧用导数定义解决了梯度下降法在非连续函数方面的计算,包括如下步骤:
步骤S1、种群初始化,设置种群大小N,目标个数M,最大迭代次数T;在目标空间中随机化粒子群,设置每个粒子的个体极值与全局极值为初始位置;设置外部档案集为空;初始化迭代次数t=0;
步骤S2、计算各粒子的目标函数值;
步骤S3、判断是否满足迭代终止准则,达到最大迭代次数时停止优化并输出结果;否则更新迭代次数t=t+1,转步骤S4;
步骤S4、判断个体的支配关系,将非支配解更新至外部档案集;
步骤S5、最优粒子选取:
i)对于pBest选取,利用粒子之间的支配关系确定,从当前粒子和历史最优中选择非支配解为pBest;若两者互不支配,则随机选取一个为pBest;利用sigma方法选取gBest。
ii)对于gBest选取,利用sigma方法选取。
其中,所述sigma方法选取gBest,包括如下步骤:
1)分别计算外部档案集中的粒子的sigma(σ)值,对于两个目标的多目标优化来说,σ可表示为:
2)从种群中选取一个粒子,计算其σi值;
3)计算该粒子σi值与外部档案集所有成员σ值之间的欧氏距离,选取欧式距离最小的外部档案集中的成员作为该粒子的gBest。
步骤S6、更新粒子的个体极值与全局极值;
步骤S7、计算所有粒子在各目标函数上的适应值,将新的非支配解更新到外部档案集中;
步骤S8、利用梯度策略更新每个粒子的位置和速度,转步骤S2。具体地,利用梯度下降法结合粒子群算法更新粒子的位置与速度的详细过程如下:
粒子群算法利用一组随机的粒子在解空间中通过迭代来搜索最 优解,利用自身最好位置pBest与群体最好位置gBest指导飞行方向,每次迭代更新寻优方向,我们希望在每次更新迭代中,通过利用每个目标函数的梯度信息来指导粒子的飞行方向,结合粒子本身的pBest、gBest来加快算法的收敛速度,得到新的粒子群速度更新公式如下所示:
vi(k+1)=ω·vi(k)+c1·γ1·[Pi-xi(k)]+c2·γ2·[Pg-xi(k)]+c3·gi(f[x(k)])
其中,c3是加速常数,由于纯梯度下降法可能会导致算法收敛到最近的局部解,所以添加了一个随机系数。f[xi(k)]为目标函数的偏导数,由下公式得到:
但是上述公式只能适用于目标函数是连续可微,对于非连续可微函数,采用导数定义,如下公式所示:
一般来说,对于Δxi的确定比较困难。本文采取比较巧妙的策略,对于粒子群算法来说,每次的更新迭代都会产生新的个体,利用xi-xi-1=Δxi的方法刚好与梯度下降法的思想相契合,解决了非连续函数梯度的计算,也加快了算法在计算梯度方面的速度,降低了算法的复杂度。
以上的所述乃是本发明的具体实施例及所运用的技术原理,若依本发明的构想所作的改变,其所产生的功能作用仍未超出说明书及附图所涵盖的精神时,仍应属本发明的保护范围。
Claims (8)
1.一种基于梯度下降法的多目标粒子群优化方法,其特征在于包括如下步骤:
(1)对粒子群种群进行初始化;在目标空间中随机化粒子群,设置每个粒子的个体极值与全局极值为初始位置;设置外部档案集为空;初始化迭代次数t=0;
(2)计算粒子群内各粒子的目标函数值,并判断目标函数值是否满足迭代终止准则,若满足最大迭代次数,则完成优化并输出结果;否则更新迭代次数t=t+1,执行步骤(3);
(3)判断粒子的支配关系,将非支配解放入外部档案集;
(4)在粒子群内选取最优粒子,更新该粒子的个体极值与全局极值;
(5)计算粒子群内所有粒子在各目标函数上的适应值,将新的非支配解更新放入外部档案集;
(6)利用梯度策略更新每个粒子的位置和速度后跳转执行步骤(2)。
2.根据权利要求1所述的一种基于梯度下降法的多目标粒子群优化方法,其特征在于:所述种群初始化包括设置种群大小N、目标个数M、最大迭代次数T。
3.根据权利要求1所述的一种基于梯度下降法的多目标粒子群优化方法,其特征在于:所述步骤(4)选取最优粒子,方法如下:
(i)对于历史最优pBest的选取,利用粒子之间的支配关系确定,从当前粒子和历史最优中选择非支配解为pBest;若粒子互不支配,则随机选取一个为pBest,利用sigma方法选取gBest;
(ii)对于全局最优gBest的选取,利用sigma方法进行选取。
4.根据权利要求3所述的一种基于梯度下降法的多目标粒子群优化方法,其特征在于:所述利用sigma方法选取gBest的方法步骤如下:
1)分别计算外部档案集中粒子的sigma(σ)值,σ表达式如下:
2)从种群中随机选取一个粒子,计算其σi值;
3)计算该粒子σi值与外部档案集的粒子σ值之间的欧氏距离,选取欧式距离最小的外部档案集的粒子作为该粒子的gBest。
5.根据权利要求1所述的一种基于梯度下降法的多目标粒子群优化方法,其特征在于:所述利用梯度策略更新粒子的位置和速度具体为利用梯度下降法结合粒子群算法更新粒子的位置与速度,方法如下:利用个体最优位置pBest与全局最优位置gBest指导粒子的飞行方向,在每次迭代更新寻优方向时,通过利用每个目标函数的梯度信息指导粒子的飞行方向,结合粒子本身的pBest、gBest得到新的粒子群速度。
6.根据权利要求1所述的一种基于梯度下降法的多目标粒子群优化方法,其特征在于:所述粒子速度的更新公式如下:
vi(k+1)=ω·vi(k)+c1·γ1·[Pi-xi(k)]+c2·γ2·[Pg-xi(k)]+c3·gi(f[x(k)])
其中,i为粒子索引;k为离散时间索引;vi(k)为第i个粒子在第k代的速度;xi(k)为第i个粒子在第k代的位置;ω为惯性权重;c1为认知加速常数;c2为社会加速常数;c3是加速常数;γ1γ2为[0,1]区间内的随机数;Pi为第i个粒子的个体最位置;Pg为粒子群群体历史最优位置;f[xi(k)]为目标函数的偏导数。
7.根据权利要求1所述的一种基于梯度下降法的多目标粒子群优化方法,其特征在于:所述f[xi(k)]由如下公式得到:
通过导数定义对上式进行变换得到f[xi(k)]的计算公式,如下所示:
8.根据权利要求1所述的一种基于梯度下降法的多目标粒子群优化方法,其特征在于:所述粒子位置的更新公式如下:
xi(k+1)=vi(k+1)+xi(k)
其中,i为粒子索引;k为离散时间索引;vi(k)为第i个粒子在第k代的速度;xi(k)为第i个粒子在第k代的位置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610573200.5A CN106228234A (zh) | 2016-07-20 | 2016-07-20 | 一种基于梯度下降法的多目标粒子群优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610573200.5A CN106228234A (zh) | 2016-07-20 | 2016-07-20 | 一种基于梯度下降法的多目标粒子群优化方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106228234A true CN106228234A (zh) | 2016-12-14 |
Family
ID=57531708
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610573200.5A Pending CN106228234A (zh) | 2016-07-20 | 2016-07-20 | 一种基于梯度下降法的多目标粒子群优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106228234A (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107341495A (zh) * | 2017-07-17 | 2017-11-10 | 深圳森阳环保材料科技有限公司 | 一种有效的井下煤尘检测系统 |
CN107589665A (zh) * | 2017-08-25 | 2018-01-16 | 中国人民解放军火箭军装备研究院 | 一种基于梯度粒子群算法的弹道快速优化方法及装置 |
CN107992936A (zh) * | 2017-12-26 | 2018-05-04 | 河北科技大学 | 基于粒子群的避险方法及飞行物的避险方法 |
CN109584099A (zh) * | 2018-10-24 | 2019-04-05 | 华中科技大学 | 电网调峰的水电站短期多目标发电计划编制方法及系统 |
CN109670655A (zh) * | 2019-01-30 | 2019-04-23 | 华北电力大学(保定) | 一种电力系统多目标粒子群优化调度方法 |
CN110298375A (zh) * | 2019-05-16 | 2019-10-01 | 同济大学 | 求解条件非线性最优扰动的并行梯度定义数据处理方法 |
CN110889531A (zh) * | 2019-08-22 | 2020-03-17 | 合肥工业大学 | 一种基于改进gsa-bp神经网络的风电预测方法及预测系统 |
CN113449433A (zh) * | 2021-07-16 | 2021-09-28 | 中存大数据科技有限公司 | 水泥生产工艺模型对应目标函数的约束优化方法及装置 |
US20220108186A1 (en) * | 2020-10-02 | 2022-04-07 | Francisco Daniel Filip Duarte | Niche Ranking Method |
-
2016
- 2016-07-20 CN CN201610573200.5A patent/CN106228234A/zh active Pending
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107341495A (zh) * | 2017-07-17 | 2017-11-10 | 深圳森阳环保材料科技有限公司 | 一种有效的井下煤尘检测系统 |
CN107589665A (zh) * | 2017-08-25 | 2018-01-16 | 中国人民解放军火箭军装备研究院 | 一种基于梯度粒子群算法的弹道快速优化方法及装置 |
CN107992936A (zh) * | 2017-12-26 | 2018-05-04 | 河北科技大学 | 基于粒子群的避险方法及飞行物的避险方法 |
CN109584099A (zh) * | 2018-10-24 | 2019-04-05 | 华中科技大学 | 电网调峰的水电站短期多目标发电计划编制方法及系统 |
CN109584099B (zh) * | 2018-10-24 | 2021-08-31 | 华中科技大学 | 电网调峰的水电站短期多目标发电计划编制方法及系统 |
CN109670655A (zh) * | 2019-01-30 | 2019-04-23 | 华北电力大学(保定) | 一种电力系统多目标粒子群优化调度方法 |
CN109670655B (zh) * | 2019-01-30 | 2023-11-07 | 华北电力大学(保定) | 一种电力系统多目标粒子群优化调度方法 |
CN110298375A (zh) * | 2019-05-16 | 2019-10-01 | 同济大学 | 求解条件非线性最优扰动的并行梯度定义数据处理方法 |
CN110889531A (zh) * | 2019-08-22 | 2020-03-17 | 合肥工业大学 | 一种基于改进gsa-bp神经网络的风电预测方法及预测系统 |
US20220108186A1 (en) * | 2020-10-02 | 2022-04-07 | Francisco Daniel Filip Duarte | Niche Ranking Method |
CN113449433A (zh) * | 2021-07-16 | 2021-09-28 | 中存大数据科技有限公司 | 水泥生产工艺模型对应目标函数的约束优化方法及装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106228234A (zh) | 一种基于梯度下降法的多目标粒子群优化方法 | |
CN105426954A (zh) | 一种基于多策略协同作用的粒子群优化的方法 | |
CN108471143A (zh) | 基于正负反馈粒子群算法的微电网多能源调度优化方法 | |
CN111090941A (zh) | 一种基于多目标优化算法的航天器最优Lambert轨道交会方法 | |
CN112926139A (zh) | 一种基于混沌映射和黄金正弦策略的改进麻雀智能优化方法 | |
CN114217524B (zh) | 一种基于深度强化学习的电网实时自适应决策方法 | |
CN104331083A (zh) | 一种航天器大角度姿态控制参数优化方法 | |
Li et al. | A hybrid assembly sequence planning approach based on discrete particle swarm optimization and evolutionary direction operation | |
CN109002948A (zh) | 基于cda-bp的微电网短期光伏发电功率预测方法 | |
CN104252132A (zh) | 基于自适应遗传算法的行星际轨道控制优化方法 | |
CN110363343A (zh) | 一种混合自适应的水电站群智能优化调度方法及系统 | |
Zheng et al. | A study of aero-engine control method based on deep reinforcement learning | |
CN111199105B (zh) | 一种扑翼运动参数优化方法 | |
Pushpalatha et al. | A succinct summary of the solar MPPT utilizing a diverse optimizing compiler | |
CN104616072A (zh) | 一种基于区间优化的提高谷氨酸发酵产物浓度的方法 | |
CN109189089A (zh) | 用于卫星地面飞控任务的智能辅助决策方法 | |
CN113471989B (zh) | 基于灰狼优化法的微电网智能二次频率控制方法 | |
CN110399697A (zh) | 基于改进遗传学习粒子群算法的飞行器的控制分配方法 | |
Chen et al. | Fuzzy PID controller optimized by improved gravitational search algorithm for load frequency control in multi-area power system | |
CN109034478A (zh) | 一种高速铁路沿线大风迭代竞争高精度预测方法 | |
CN106251031A (zh) | 一种基于生物启发的改进型粒子群优化算法 | |
CN116859728A (zh) | 一种基于改进多目标粒子群算法的船舶动力定位系统推力分配方法 | |
Ding | Research of improved particle swarm optimization algorithm | |
CN106611380A (zh) | 一种改进的帝国主义竞争算法求解作业车间调度问题 | |
CN111552317B (zh) | 一种多航天器四维协同轨迹确定方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20161214 |
|
WD01 | Invention patent application deemed withdrawn after publication |