CN111461385A - 一种基于Lévy飞行的粒子群算法的生产计划优化方法 - Google Patents
一种基于Lévy飞行的粒子群算法的生产计划优化方法 Download PDFInfo
- Publication number
- CN111461385A CN111461385A CN201911357022.2A CN201911357022A CN111461385A CN 111461385 A CN111461385 A CN 111461385A CN 201911357022 A CN201911357022 A CN 201911357022A CN 111461385 A CN111461385 A CN 111461385A
- Authority
- CN
- China
- Prior art keywords
- production plan
- formula
- particle
- particle swarm
- swarm algorithm
- 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
Links
- 239000002245 particle Substances 0.000 title claims abstract description 69
- 238000004519 manufacturing process Methods 0.000 title claims abstract description 53
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000005457 optimization Methods 0.000 title claims abstract description 21
- 230000006870 function Effects 0.000 claims description 15
- 238000003860 storage Methods 0.000 claims description 12
- 238000009826 distribution Methods 0.000 claims description 9
- 230000008859 change Effects 0.000 claims description 7
- 150000001875 compounds Chemical class 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 5
- 210000001744 T-lymphocyte Anatomy 0.000 claims description 3
- 230000007423 decrease Effects 0.000 claims description 3
- 240000007594 Oryza sativa Species 0.000 claims 1
- 235000007164 Oryza sativa Nutrition 0.000 claims 1
- 235000009566 rice Nutrition 0.000 claims 1
- 230000008901 benefit Effects 0.000 description 3
- 238000007792 addition Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 239000003960 organic solvent Substances 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- 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
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/04—Manufacturing
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/30—Computing systems specially adapted for manufacturing
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Marketing (AREA)
- Tourism & Hospitality (AREA)
- Health & Medical Sciences (AREA)
- General Business, Economics & Management (AREA)
- Game Theory and Decision Science (AREA)
- Biomedical Technology (AREA)
- Development Economics (AREA)
- Manufacturing & Machinery (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Life Sciences & Earth Sciences (AREA)
- Artificial Intelligence (AREA)
- Primary Health Care (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- Evolutionary Computation (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明为一种基于Lévy飞行的粒子群算法的生产计划优化方法,包括以下步骤:对生产计划进行数学建模;初始化粒子群的各种参数以及粒子的位置;计算各粒子的适应度值,对个体最优和全局最优进行更新;进行粒子位置更新,前期大概率采用本发明提出的基于Lévy飞行的改进粒子群算法,后期逐渐减小使用概率;不断迭代直至达到收敛条件,从而得到生产计划的解。本发明的优点是:将Lévy飞行与粒子群算法相结合并运用于生产计划的研究,有效防止粒子陷入局部最优,提高解的精度和收敛速度,可以适应于非线性问题。
Description
技术领域
本发明属于生产计划优化领域,尤其涉及一种基于Lévy飞行的粒子群算法的生产计划优化方法。
背景技术
随着近些年来流程工业的整体市场情况的波动,间歇式生产模式逐渐兴起,其能够带给流程工业的优势也愈发明显。流程工业的企业面对复杂多变的市场环境,想要提高自身的竞争力,就必须根据市场需求制定合理的生产计划,生产计划的目的是要在充分考虑市场需求与企业生产能力的情况下,调整企业各部门间的协作关系来提高效益。生产计划的模型应由不同部门间的费用模型组建,包括生产费用、运输费用、短缺费用和仓储费用,并进行统筹规划,使得这些费用的总和最小。随着生产计划规模的扩大,生产计划求解的难度急剧加大,传统粒子群算法难以求得高质量的解,对于一些非线性问题传统解法也无能为力。
发明内容
本发明主要解决了基本粒子群算法优化生产计划易陷入局部最优、收敛速度慢、精度低的缺点,提供了一种能够提高生产计划求解的速度和精度,对于非线性问题亦有一定优势的基于Lévy飞行的粒子群算法的生产计划优化方法。
本发明解决其技术问题所采用的技术方案是,一种基于Lévy飞行的粒子群算法的生产计划优化方法,包括以下步骤:
步骤1:建立生产计划模型;
步骤2:初始化参数,该参数包括种群的粒子数N、维度d、学习因子c1和c2、迭代次数ger、初始概率P0、结束概率P1、初始惯性权重ω0和结束惯性权重ω1;
步骤3:初始化粒子的位置,并将粒子的个体最优位置pid设置为当前位置;
步骤4:计算各个粒子的目标函数值,如果fit(xid)优于fit(pid),则更新pid为xid的当前值,否则pid不变;
其中,fit(xid)和fit(pid)分别为目标函数在xid和pid的适应度值,xid为粒子在当前迭代中第d维的值,pid为粒子在当前迭代过程中第d维的个体最优值;
步骤5:如果fit(xid)优于fit(pgd),则更新pgd为xid的当前值,否则pgd不变;
其中,fit(pgd)为目标函数在当前全局最优点的适应度值;
步骤6:进行粒子位置更新;
步骤7:判断是否达到最大迭代次数/达到目标函数的收敛精度,若是,则结束,输出结果;若否,返回步骤4。
以生产计划模型中各个影响生产成本的因素集合为粒子群,通过粒子群算法求解各个因素最优值,以实现生产成本最优。
作为上述方案的一种优选方案,所述的目标函数如下:
式中,n表示生产调度周期,N表示周期集合,s表示产品,S表示产品集合,α,β,γ,δ分别表示所对应种类产品生产费用单价、运输费用单价、短缺费用单价及贮存费用单价,分别表示在调度周期n内产品s的生产量、运输量、短缺量和贮存量。
作为上述方案的一种优选方案,所述生产计划模型如下:
作为上述方案的一种优选方案,所述步骤6中采用加入Lévy飞行的位置更新公式和基本粒子群位置更新公式中的一种位置更新公式对粒子位置进行更新。采用两种位置更新公式更新粒子位置,解决了传统粒子群算法易陷入局部最优的问题。
作为上述方案的一种优选方案,所述采用加入Lévy飞行的位置更新公式进行粒子位置更新的概率为P,采用基本粒子群位置更新公式进行粒子位置更新的概率为1-P。概率P的公式如下:
式中,P0和P1分别为开始和结束时的概率值,i为当前的迭代次数,ger为总迭代次数。根据概率来决定采用何种位置更新公式,且这个概率随着迭代次数的改变而改变,加快了收敛速度,提高了收敛精度。
作为上述方案的一种优选方案,所述加入Lévy飞行的位置更新公式如下:
式中,c1和c2是学习因子,为常数;r1和r2为[0,1]内均匀分布的随机数;和分别是粒子在第k+1次和第k迭代中第d维的位置,两者均有上下限,是粒子pi在第k次迭代中第d维的个体最优位置,pgd是种群的第d维的全局最优位置,ω为惯性权重,α0为系数,常取α0=0.01,为点乘积,Levy(β)为随机搜索路径,表示服从参数为β的莱维分布,Levy(β)~|t|-β,0<β≤3。
作为上述方案的一种优选方案,采用Mantegna对莱维分布进行模拟,表达式如下:
式中,μ,v为正态分布,定义μ~N(0,σ2),v~N(0,1);
式中,Γ(·)为Gamma函数,常取β=1.5。
作为上述方案的一种优选方案,所述惯性权重ω随着迭代次数增加逐渐减小,其公式如下:
式中,ω0和ω1分别为初始和最终的惯性权重值,i为当前的迭代次数,ger为总迭代次数。
本发明的有益效果是:采用两种粒子群位置更新方式,根据概率来决定采用何种更新方式,能够更好的跳出局部最优,收敛到全局最优,在加快了收敛速度的同时提高了收敛精度,能够更快的获得最优生产计划。
附图说明
图1为实施例中基于Lévy飞行的粒子群算法的生产计划优化方法的一种流程示意图。
具体实施方式
下面通过实施例,并结合附图,对本发明的技术方案作进一步的说明。
实施例:
本实施例一种基于Lévy飞行的粒子群算法的生产计划优化方法,如图1所示,包括以下步骤:
步骤1:建立生产计划模型,本实施例中考虑的因素包括生产量、运输量、短缺量和贮存量,最终建立的生产计划模型如下:
步骤2:初始化参数,该参数包括种群的粒子数N、维度d、学习因子c1和c2、迭代次数ger、初始概率P0、结束概率P1、初始惯性权重ω0和结束惯性权重ω1;
步骤3:初始化粒子的位置,并将粒子的个体最优位置pid设置为当前位置;
步骤4:计算各个粒子的目标函数值,目标函数如下:
式中,n表示生产调度周期,N表示周期集合,s表示产品,S表示产品集合,α,β,γ,δ分别表示所对应种类产品生产费用单价、运输费用单价、短缺费用单价及贮存费用单价,分别表示在调度周期n内产品s的生产量、运输量、短缺量和贮存量;如果fit(xid)优于fit(pid),则更新pid为xid的当前值,否则pid不变;其中,fit(xid)和fit(pid)分别为目标函数在xid和pid的适应度值,xid为粒子在当前迭代中第d维的值,pid为粒子在当前迭代过程中第d维的个体最优值,本实施例中fit(xid)小于fit(pid)即为fit(xid)优于fit(pid);
步骤5:如果fit(xid)优于fit(pgd),则更新pgd为xid的当前值,否则pgd不变;fit(pgd)为目标函数在当前全局最优点的适应度值;
步骤6:进行粒子位置更新;粒子位置更新采用采用加入Lévy飞行的位置更新公式和基本粒子群位置更新公式中的一种,采用加入Lévy飞行的位置更新公式进行粒子位置更新的概率为P,采用基本粒子群位置更新公式进行粒子位置更新的概率为1-P。概率P的公式如下:
式中,P0和P1分别为开始和结束时的概率值,i为当前的迭代次数,ger为总迭代次数;本实施例中采用生成随机数r的形式实现位置更新公式的选择,当随机数r>P时采用原基本粒子群位置更新公式,当r≤P时采用加入Lévy飞行的位置更新公式,r为[0,1]内的均匀随机数,本实施例中,P0=0.8,P1=0。
其中,加入Lévy飞行的位置更新公式如下:
式中,c1和c2是学习因子,为常数;r1和r2为[0,1]内均匀分布的随机数;和分别是粒子在第k+1次和第k迭代中第d维的位置,两者均有上下限,是粒子pi在第k次迭代中第d维的个体最优位置,pgd是种群的第d维的全局最优位置,ω为惯性权重,α0为系数,常取α0=0.01,为点乘积,Levy(β)为随机搜索路径,表示服从参数为β的莱维分布,Levy(β)~|t|-β,0<β≤3,莱维分布十分复杂,实现难度高,因此本实施例中采用Mantegna方法模拟进行模拟,Mantegna方法表达式如下:
式中,μ,v为正态分布,定义μ~N(0,σ2),v~N(0,1)。
式中,Γ(·)为Gamma函数,常取β=1.5。
惯性权重ω随着迭代次数增加逐渐减小,其公式如下:
式中,ω0和ω1分别为初始和最终的惯性权重值,i为当前的迭代次数,ger为总迭代次数。
步骤7:判断是否达到最大迭代次数或达到目标函数的收敛精度,若是,则结束,输出结果,输出的结果pgd包括当前调度周期内产品的生产量、运输量、短缺量和贮存量;若否,返回步骤4。
本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代,但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。
Claims (8)
步骤1:建立生产计划模型;
步骤2:初始化参数,该参数包括种群的粒子数N、维度d、学习因子c1和c2、迭代次数ger、初始概率P0、结束概率P1、初始惯性权重ω0和结束惯性权重ω1;
步骤3:初始化粒子的位置,并将粒子的个体最优位置pid设置为当前位置;
步骤4:计算各个粒子的目标函数值,如果fit(xid)优于fit(pid),则更新pid为xid的当前值,否则pid不变;
其中,fit(xid)和fit(pid)分别为目标函数在xid和pid的适应度值,xid为粒子在当前迭代中第d维的值,pid为粒子在当前迭代过程中第d维的个体最优值;
步骤5:如果fit(xid)优于fit(pgd),则更新pgd为xid的当前值,否则pgd不变;
其中,fit(pgd)为目标函数在当前全局最优点的适应度值;
步骤6:进行粒子位置更新;
步骤7:判断是否达到最大迭代次数/达到目标函数的收敛精度,若是,则结束,输出结果;若否,返回步骤4。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911357022.2A CN111461385B (zh) | 2019-12-25 | 2019-12-25 | 一种基于Lévy飞行的粒子群算法的生产计划优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911357022.2A CN111461385B (zh) | 2019-12-25 | 2019-12-25 | 一种基于Lévy飞行的粒子群算法的生产计划优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111461385A true CN111461385A (zh) | 2020-07-28 |
CN111461385B CN111461385B (zh) | 2022-04-22 |
Family
ID=71683291
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911357022.2A Active CN111461385B (zh) | 2019-12-25 | 2019-12-25 | 一种基于Lévy飞行的粒子群算法的生产计划优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111461385B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114819577A (zh) * | 2022-04-18 | 2022-07-29 | 杭州昌邦科技有限公司 | 面向市场需求不确定性的生产计划与调度集成优化方法 |
CN116415480A (zh) * | 2022-09-06 | 2023-07-11 | 中国人民解放军海军航空大学 | 一种基于ipso的飞机海上平台出动离场规划方法 |
CN116757446A (zh) * | 2023-08-14 | 2023-09-15 | 华中科技大学 | 基于改进粒子群算法的梯级水电站调度方法及系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104332995A (zh) * | 2014-11-14 | 2015-02-04 | 南京工程学院 | 基于改进粒子群算法的配电网重构优化方法 |
CN107437813A (zh) * | 2017-07-17 | 2017-12-05 | 安徽大学 | 基于布谷鸟‑粒子群的配电网无功功率优化方法 |
CN109544234A (zh) * | 2018-11-22 | 2019-03-29 | 安徽师范大学 | 一种云环境下基于帕累托最优的资源定价方法 |
KR102009608B1 (ko) * | 2018-09-19 | 2019-08-09 | 인하대학교 산학협력단 | 센서 네트워크에서 uav의 센싱 데이터 정보 최대화를 위한 pso 기반 uav 위치 탐색 방법 및 장치 |
-
2019
- 2019-12-25 CN CN201911357022.2A patent/CN111461385B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104332995A (zh) * | 2014-11-14 | 2015-02-04 | 南京工程学院 | 基于改进粒子群算法的配电网重构优化方法 |
CN107437813A (zh) * | 2017-07-17 | 2017-12-05 | 安徽大学 | 基于布谷鸟‑粒子群的配电网无功功率优化方法 |
KR102009608B1 (ko) * | 2018-09-19 | 2019-08-09 | 인하대학교 산학협력단 | 센서 네트워크에서 uav의 센싱 데이터 정보 최대화를 위한 pso 기반 uav 위치 탐색 방법 및 장치 |
CN109544234A (zh) * | 2018-11-22 | 2019-03-29 | 安徽师范大学 | 一种云环境下基于帕累托最优的资源定价方法 |
Non-Patent Citations (1)
Title |
---|
王荣 等: "基于改进莱维飞行粒子群算法的光伏系统MPPT方法", 《南昌大学学报(工科版)》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114819577A (zh) * | 2022-04-18 | 2022-07-29 | 杭州昌邦科技有限公司 | 面向市场需求不确定性的生产计划与调度集成优化方法 |
CN116415480A (zh) * | 2022-09-06 | 2023-07-11 | 中国人民解放军海军航空大学 | 一种基于ipso的飞机海上平台出动离场规划方法 |
CN116415480B (zh) * | 2022-09-06 | 2024-01-19 | 中国人民解放军海军航空大学 | 一种基于ipso的飞机海上平台出动离场规划方法 |
CN116757446A (zh) * | 2023-08-14 | 2023-09-15 | 华中科技大学 | 基于改进粒子群算法的梯级水电站调度方法及系统 |
CN116757446B (zh) * | 2023-08-14 | 2023-10-31 | 华中科技大学 | 基于改进粒子群算法的梯级水电站调度方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN111461385B (zh) | 2022-04-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111461385B (zh) | 一种基于Lévy飞行的粒子群算法的生产计划优化方法 | |
CN108133258B (zh) | 一种混合全局优化方法 | |
Li et al. | A multi-information fusion “triple variables with iteration” inertia weight PSO algorithm and its application | |
CN104408518B (zh) | 基于粒子群优化算法的神经网络学习优化方法 | |
CN108710905A (zh) | 一种基于多模型联合的备件数量预测方法和系统 | |
CN106022521A (zh) | 基于Hadoop架构的分布式BP神经网络的短期负荷预测方法 | |
CN102831269A (zh) | 一种流程工业过程工艺参数的确定方法 | |
CN102510059A (zh) | 基于bp神经网络的超短期风电功率预测方法 | |
CN108596242A (zh) | 基于小波神经网络和支持向量机的电网气象负荷预测方法 | |
CN108062585A (zh) | 一种基于飞蛾扑火算法计算函数最值的方法 | |
CN107832983B (zh) | 基于动态规划与遗传算法的铸造熔炼批计划与调度方法 | |
CN111553118A (zh) | 基于强化学习的多维连续型优化变量全局优化方法 | |
CN114880806A (zh) | 基于粒子群优化的新能源汽车销量预测模型参数优化方法 | |
CN112884236A (zh) | 一种基于vdm分解与lstm改进的短期负荷预测方法及系统 | |
CN109145399A (zh) | 一种基于改进的粒子滤波算法的疲劳裂纹扩展预测方法 | |
CN112836885B (zh) | 组合负荷预测方法、装置、电子设备和存储介质 | |
CN113901728A (zh) | 基于迁移遗传算法的电脑第二类装配线平衡优化方法 | |
CN109670634A (zh) | 一种基于非等长粒度化特征的煤气流量数据预测方法 | |
CN109902850A (zh) | 确定库存控制策略的方法、装置及存储介质 | |
WO2024113585A1 (zh) | 一种面向离散制造系统的智能交互式决策方法 | |
CN112766609A (zh) | 一种基于云计算的用电量预测方法 | |
CN104090932A (zh) | 一种内容推荐方法及装置 | |
CN105320991A (zh) | 基于因素的钢铁企业工序能耗预测方法 | |
CN114648178B (zh) | 一种基于ddpg算法的电能计量装置运维策略优化方法 | |
CN108288105A (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 |