油库人员的优化配制(共26页).doc

《油库人员的优化配制(共26页).doc》由会员分享，可在线阅读，更多相关《油库人员的优化配制(共26页).doc（26页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、精选优质文档-倾情为你奉上2011年河南科技大学数学建模竞赛选拔承 诺 书我们仔细阅读了数学建模竞赛选拔的规则.我们完全明白，在做题期间不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人研究、讨论与选拔题有关的问题。我们知道，抄袭别人的成果是违反选拔规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守选拔规则，以保证选拔的公正、公平性。如有违反选拔规则的行为，我们将受到严肃处理。我们选择的题号是（从A/B/C中选择一项填写）： C 队员签名 ：1. 方健 2. 孙亚娜 3. 张腾

2、月 日期： 2011 年 8 月 22 日2011年河南科技大学数学建模竞赛选拔编 号 专 用 页评阅编号（评阅前进行编号）：评阅记录（评阅时使用）：评阅人评分备注专心-专注-专业油库人员的优化配制摘要这篇文章是为了解决某油库人员配置问题的。由于油库工作性质的特殊性，工作人员需要经过合理地调配，油库才能高效的运转。我们讨论的油库中，油库管理者设置了计量与质量检测管理、收发油料管理、装备维护与维修管理、安全保障管理和服务保障管理等大类岗位，工作项目复杂，类目繁多。因此要进行合理的人员配置优化，才能使油库正常运转。第一题要求各大类均设专职，同类目可兼职。因此我们从人员工作变动入手，建立时间节点序列

3、，从而找出了各约束条件，接着我们采用整数分析法，借助lingo软件来求得最优解。第二题我们用与第一题类似的方法进行求解，只是约束条件更加复杂。我们仍借助lingo软件求最优解。第三题由于约束条件过于繁琐，我们采用模拟退火法进行优化，用matlab软件求得最优解。关键字：整数规划，最优解，目标函数，工作量，优化配置。一、 问题重述油库是储存、供应油料及油料装备的重要机构。由于油料是一种易燃、易爆、易挥发、易渗漏，并有一定腐蚀作用的物质。因此，一个油库管理工作要保证正常的运行，必须要设置计量与质量检测管理、收发油料管理、装备维护与维修管理、安全保障管理、服务保障管理等相关的岗位和人员配制。某油库现

4、有不同功能、不同规格的大小储油罐80个，储油量达16万立方米以上，年收发油量达7000多立方米，工作任务十分繁重。根据实际需要按工种分类，油库的工作岗位可以分为大类：（A）计量与质量检测管理；（B）收发油料管理；（C）设备维护与维修管理；（D）安全保障管理；（E）服务保障管理。由于油库工作的性质要求，每一大类都包括若干个具体的工作岗位，每个岗位都需要数量不等的人员和工作量，附表给出了各类工作岗位、所需要的人员数量和全年的工作量。通常油库的所有人员在保障油库正常运行的条件下，还要参加必要业务学习等活动，所以实际要求每个人全年累积从事油库管理相关工作的总工作量不超过175天（每天按小时计算），除节

5、假日外，其他时间用于业务学习等活动。请你们研究解决下列问题：（）根据油库正常的工作任务需要，如果要求（A）、（B）、（C）和（D）类人员都配专职，同类中的各工作岗位人员可以兼职。那么各类工作岗位最少需要配制多少人员？平均年工作量是多少？（）考虑到人员和编制的问题，有关部门提出除了（A）、（B）两大类专业性较强的工作专职专人外，（C）、（D）两大类在时间允许的情况下可以相互兼职。那么这个油库最少需要多少人员才能保证油库的正常运行？并说明各类人员的年总工作量为多少？（）如果油库的所有人员都经过了专业的培训，每个人都能从事（A）、（B）、（C）和（D）类中的任何一项工作，每一个岗位都不设专职人员，那

6、么在时间允许的情况下，最少需要多少人员能够保证油库的正常运行？并说明各类人员的年总工作量为多少？（）你对该油库的岗位和人员的配置、提高油库的管理水平和运行效率等方面有什么合理化建议。附表：油库大类岗位、需要人员及年工作量二、 问题分析问题一的任务要求是通过形成A、B、C、D类各自配有专职，各类工作中配有兼职的工作模式以保证油库各项工作正常进行，同时满足这年总工作量不超过175天，来完成各类工作的人员优化配置求解。因为全为收油和零发油相关的工作，根据问题要求，所有收油与零发油的相关任务依赖于设备操作，其岗位所需的人数固定不变，且必须同时在岗，故本题可把所有收油任务归为一类，所有零发油任务归为一类

7、，且零发油的时间不确定， 且收油和法油的相关的所有任务可以不再同一时刻进行，故可以让收油的工人兼职与零发油的相关工作，所以对于B类就比较容易求得人员配置的最少人数。题中的问题都有涉及有最少人员的求解。由分析可知，A、C、D几项工作中，能明显的判断出C类与D类工作12月份参加工作的人数最多。又知题目要求求解最少配置人数，所以在最多工作人员参与工作的前提下，求解的最少配置人员，才可以保证工作的正常运行。所以对于C类与D类工作选择12月份来进行变动节点的分析更加合理和准确。然而由分析知道，对于A类，由于单双月工作的人数差不多，所以可以考虑单双月来进行模型的建立。对于问题一和二，我们选择整数规划模型来

8、作为求解的最优模型，并利用lingo软件来求解结果。对于问题三我们选择遗传算法来建立模型并用matlab来求解各项工作人员的最小值。三、 模型的假设1零发油操作相关人员限制 2. 收油相关工作合并3. 工作持续执行时间，员工保持不变4. 员工能正常工作，不考虑请假问题。5. 附表中的信息是2011年的信息。6. 员工的工作时间用一天来进行计算，不足一天按一天来计算。四、 符号说明-A类工作单月优化配置的最少人员数-A类工作双月优化配置的最少人员数-C类工作优化配置的最少人员数-D类工作优化配置的最少人员数-CD类工作优化配置的最少人员数-A类工作中第i项工作在时刻的在岗人数-C类工作中第i项工

9、作在时刻的在岗人数-D类工作中第i项工作在时刻的在岗人数-CD类工作中第i项工作在时刻的在岗人数-第j个人员的变动节点-A项工作岗位上单月需配置的工作人员数-A项工作岗位上双月需配置的工作人员数-C项工作岗位上需配置的工作人员数-D项工作岗位上需配置的工作人员数-A类工作岗位上年总工作量-C类工作岗位上年总工作量-D类工作岗位上年总工作量-CD类工作岗位上年总工作量五、 模型建立与求解问题一前面问题分析中已经说明选择C类12份作为基准来作为人员的变动节点来求解最少人员的数值解。所以首先分析C类工作，来求解各小项任务所需人员的最少人数。由题意可知，为保证油库各项工作的正常进行，同类的各项工作岗位

10、上的人员可以相互兼职，即某工作人员在与自己工作时间不相冲突，且其年总工作量未超过175 天时可帮助完成同类工作中的其他工作，进而可使各类工作配置最优人员。C 类工作的人员优化配置求解对于C类分析可知，可以明显的看出12月份的工作人员最多，C2的工作时间不确定，C6是单月工作，所以选择C1、C3、C4、C5、C7这五项的工作起始时间作为人员的变动节点，其中将C2、C6的人数闲置处理。查阅2011年的日历，结合上面分析的五项工作起始时间可得人员变动节点序列为1、2、5、9、10、12、14、15、16、19、23、26、28、30,根据时间变动序列以工作配置的总工作人员最少为目标函数，建立如下的整

11、数规划模型：我们规定C1、C3、C4、C5、C7在j某时刻的在岗人数为、。 Subject to:，其中：(j=2,3,14) 上述模型中，表示C类工作中所有的员工的年平均总工作量不能超过175天。表示时刻C3项工作已经完成，另外一种情况是没有完成。对于道理一样。表示对于C612月没有人员工作要加上C2项工作中的两人的人数。按照上述约束条件，我们通过lingo软件（附录1）求得最优解14人。 D类工作的人员优化配置对于D类工作，有D1、D2、D3、D4、D5、D6六项基本工作，经分析，12月份上班人数最多，所以采用和C类相同的办法，以12月份来结合附表上的工作起始时间来确定人员变动节点序列。其

12、中，D2的工作时间不确定，所以将其人数闲置。同时D6项工作需一人工作241天，又知要求是没人的工作量不得超过175天，所以将闲置一人作为专职。所以由D1、D3、D4、D5这四项的工作起始时间作为人员的变动节点。所以可以得到1，3，5，7，9，10，11，12，13，15，17，19，21，23，25，26，27，29，31，然后建立人员优化的配置模型，其中我们规定，D1、D3、D4、D5在j时刻的在岗人数为、所以可得整数规划模型如下：其中，述模型中，表示D工作中所有的员工的年平均总工作量不能超过175天。表示时刻D3工作已经完成，另外一种情况是没有完成。对于道理一样。表示D2闲置的6人加上D6

13、闲置的1人总共7人，然后加上其他所有的最多在岗人数。同样我们通过lingo软件，求得最优解33人。（附录1）A类工作的人员优化配置对A类工作，有A1、A2、A3、A4、A5、A6共计6 项具体的工作，经分析，这类工作在单双月中上班的人数差不多，故分别以单双月中这6 项工作的起始时间作为人员的变动节点序列，并建立与C 类工作人员优化配制相同的整数规划模型。（1）A类的单月的人员优化配置单月的人员变动节点序列是10,15，25 A 类工作在单月的人员优化配制模型如下：Subject to 其中，求得结果：最优解为3人。（2）A类双月的人员优化配置双月的人员变动节点为1，10，25，所以在双月的人员

14、优化配置如下：通过lingo软件可求得最优解：单月为3人，双月为4人。（附录1）因此，A类应具备最少人数4人。B类的人员优化配置B类工作是依赖于设备的收发油的相关工作，其岗位所需的人数固定不变，且必须同时在岗，所以，把所有收油任务归为一类，所有发油任务归为一类，且收油发油的工作可以不在同一时间进行，所以让收油工人代替发油工人工作，可得知与收油相关的需要的最少工作人员是23人，即B类工作最少需要人员为23人。E类的人员优化配置对于E 类工作，其具体的4 项工作均为全专职工作，为避免员工总的工作量超过175 天，故可根据公式求得y（y表示使E类某一项具体工作所需最少人数）根据上述函数，将题给数据代

15、入求得E2为 5 人， E3 为9 人， E4为17 人。通过上述的求解，得出油库的除伙食服务保障外所需的员工数为125 人，根据题给要求可以得出E1所需人数为20 人。因此E类共需51人。问题二 问题二中，C.D两大类工作在时间允许的情况下可以相互兼职，故本文将C、D归为一类，称为CD类。由问题一可知，CD两大类均在12月份得工作人数最多。且C2、D2项工作的工作人员需闲置，又因为D6项工作要求除节假日外每天一次，为避免其工作人员的工作量超过175天，所以D6项需要让一人闲置出来。综上，应以CD类除C2、D2、D6项工作的其他项工作在12月份得起始时间为变动节点，得到变动节点序列1，2,3,

16、5,7,9,10,11,12,13,14,15,16,17,19,21,23,25,26,27,28,29,30,31并建立整数规划模型如下:特别说明，C1、C3、C4、C5、C7、D1、D3、D4、D5在j时刻的在岗人数分别为、Subject to:因为A、B类工作的工作人员只能在同类工作岗位上兼职，所以与问题一的结果相同，E类工作所需的最少工作的人员的求解和问题一也相同。对于CD类工作，通过lingo软件可以计算出最终结果，其最小人数为69人。（附表1）ABCDE总人数4236951145问题三：如果油库的所有人员都经过了专业的培训，每个人都能从事（A）、（B）、（C）和（D）类中的任何一

17、项工作，每一个岗位都不设专职人员，那么在时间允许的情况下，最少需要多少人员能够保证油库的正常运行？并说明各类人员的年总工作量为多少？一：模型的建立该问题要求A、B、C、D四类工作合为一类，因此我们将这四类工作作为一个大类进行优化。而这四项工作中的类目繁多，约束条件冗杂，不适合进行简单的整数规划或者枚举法。因此我们采用现代优化算法中的模拟退火算法（SA）来优化这个较大规模的模型。运用模拟退火算法，首先要确定研究范围的解空间，确定目标函数，经多次降温后求得最优解。这个问题我们分别从单月和双月进行分析，下面我们以双月为例，优化过程如下：（一）、建立解空间：本题的解空间由一个阶矩阵构成：，其中为第项工

18、作在时刻工作的人数。由附表可知，双月中12月的工作量最大，因此我们只分析12月的最小人数，即可适用于所有双月的最小人数。由于零发油及其相关工作时间不确定，为了避免该类工作发生时无人上岗，我们将这些人数预置，不参与优化；同样由于“油库环境保养与卫生工作”只在节假日休息，我们也将1人作为专职（该人不参与优化），以避免个人工作量超出175 天。此外，由于双月10 号的收油及其相关工作依赖于设备操作，这些岗位所需人员数固定不变，必须同时在岗，因此，我们将这一类工作整合为一项工作进行分析。最后，由模型假设，我们将不到一天的工作时间也记为一天。我们以题目附表中要求的最小人数作为初始解,=（二）、确定目标函

19、数依题意，优化的最终目标是求得总人数的最小值。由此建立目标函数：（三）、产生新解和接受新解1、运用matlab软件，可随机的对某一项工作的人数进行调整，调整后可得新解；2、将新解与原解进行比差，得到代价函数差：3、新解接收准则：当新解对应的函数值较原解小时（说明函数值减小），此时新解被无条件接受；当新解对应的函数值较原解大时（说明函数值增大），此时新解被以某一概率接受。即对于某一个温度 和该优化问题的一个解，可以生成。接受作为下一个新解的概率为：其中，为对应初始解的一个参数，为对应解的一个参数。4、降温过程：在温度 下，经过很多次的转移之后，降低温度得到 。在 下重复上述过程。因此整个优化过程

20、就是不断寻找新解和缓慢降温的交替过程。最终的解是对该问题寻优的结果。在matlab软件中，我们设降温系数=0.999，这样缓慢的降温过程会使结果更加精确。（四）、得到最优解由上述过程，我们看到在每个下，所得到的一个新状态完全依赖于前一个状态，可以和前面的状态, 无关，因此这是一个马尔可夫过程。使用马尔可夫过程对上述模拟退火的步骤进行分析，结果表明：从任何一个状态生成的概率，在的邻域中是均匀分布的，且新状态被接受的概率满足上述公式，那么经过有限次的转换，在温度下的平衡态的分布由下式给出：当温度最终降为0时，的分布变为：D为矩阵S的解空间。通常在matlab软件中，我们选定的终止温度t，判断退火过

21、程是否结束。若T =2; c(2,1)=12;c(4,8)=3;for(time(j)|j#ne#1#and#j#ne#2#and#j#ne#5#and#j#ne#7#and#j#ne#8#and#j#ne#13:c(3,j)=12); for(time(j)|j#eq#7#or#j#eq#13:c(5,j)=2);c(4,9)=if(c(4,8)*(t(9)-t(8)#ge#m(5),0,c(4,8); for(time(j)|j#ge#2:c(2,j)=if(t(j)-t(1)*c(2,1)#ge#m(3),0,c(2,j-1); for(time(j):total(j)=sum(work

22、(i):c(i,j); c1=max(time(j):total(j); min=c1+2; (c1+2)*175=mc; End结果：D类：model: title D组12月最少人数分配; sets: work/1,2,3,4/:m; time/1.19/:t,total; links(work,time):d; endsets data: m=6,9,3,20; t=1,3,5,7,9,10,11,12,13,15,17,19,21,23,25,26,27,29,31; md=4921; md3=3;enddata for(time(j)|j#ne#6:d(1,j)=0); for(ti

23、me(j)|j#ne#3#and#j#ne#8#and#j#ne#12#and#j#ne#16#and#j#ne#17:d(2,j)=0); for(time(j)|j#ne#3#and#j#ne#8#and#j#ne#12#and#j#ne#16:d(3,j)=0); for(time(j)|j#eq#6#or#j#eq#8#or#j#eq#16:d(4,j)=0); d(1,6)=6; for(time(j)|j#ne#6#and#j#ne#8#and#j#ne#16:d(4,j)=20); d(2,3)=6;d(2,8)=6;d(2,12)=6;d(2,16)=6; d(2,9)=if(

24、t(9)-t(8)*d(2,8)#ge#md3,0,d(2,8); d(2,17)=if(t(17)-t(16)*d(2,16)#ge#md3,0,d(2,16); for(time(j):total(j)=sum(work(i):d(i,j); d1=max(time(j):total(j); min=d1+7; (d1+7)*175=md; End A单月：model: title A组单月最少人数分配; sets: work/1,2,3/:m; time/1,2,3/:t,total; links(work,time):a; endsets data:m=4,30,5; t=10,15,

25、25; Ma=258; enddata a(1,1)=0;a(1,2)=0;a(2,1)=0;a(3,2)=0;a(3,3)=0; a(1,3)=2;a(2,2)=2;a(3,1)=2; a(2,3)=if(a(2,2)*(t(3)-t(2)#ge#m(2),0,a(2,2); for(time(j):total(j)=sum(work(i):a(i,j); a1=max(time(j):total(j); 175*a1=Ma; min=a1; end A双月：model: title A组双月最少人数分配; sets: work/1,2,3,4/; time/1,2,3/:t,total;

26、links(work,time):a; endsets data: t=1,10,25; Ma=258; enddata a(1,1)=0;a(1,2)=0; a(2,2)=0;a(2,3)=0;a(3,1)=0; a(3,3)=0; a(4,1)=0; a(4,3)=0; a(1,3)=2;a(2,1)=2;a(3,2)=2;a(4,2)=2; for(time(j):total(j)=sum(work(i):a(i,j); a1=max(time(j):total(j); 175*a1=Ma; min=a1; end Cd:model: title CD 组混合人数最优配置; sets:

27、work /1.9/:m; time /1.24/:t,total; links(work,time):cd; endsets data: m=2,360,12,12,1,6,9,3,20; t=1,2,3,5,7,9,10,11,12,13,14,15,16,17,19,21,23,25,26,27,28,29,30,31; mcd=6879; enddata for(time(J)|J#ne#7:cd(1,J)=0); cd(2,24)=0;for(time(J)|J#ne#2 #and# J#ne#4 #and# J#ne#6#and# J#ne#9 #and# J#ne#13#and#

28、 J#ne#15#and# J#ne#17 #and# J#ne#19 #and# J#ne#23 :cd(3,J)=0); for(time(J)|J#ne#12#and#j#ne#13#and#j#ne#14:cd(4,J)=0); for(time(J)|J#ne#11 #and# J#ne#21 :cd(5,J)=0); for(time(J)|J#ne#7:cd(6,J)=0); for(time(J)|J#ne#4 #and# J#ne#9 #and# J#ne#10 #and# J#ne#15 #and# J#ne#19 #and# J#ne#20:cd(7,J)=0); for

29、(time(J)|J#ne#4 #and# J#ne#9 #and# J#ne#15 #and# J#ne#19 :cd(8,J)=0); cd(9,2)=0;cd(9,7)=0;cd(9,9)=0;cd(9,11)=0;cd(9,13)=0;cd(9,19)=0;cd(9,21)=0;cd(9,23)=0; cd(3,2)=12;cd(3,4)=12;cd(3,7)=12; cd(3,9)=12;cd(3,11)=12;cd(3,13)=12; cd(3,16)=12;cd(3,18)=12;cd(3,20)=12; cd(1,7)=2;cd(2,1)=12;cd(4,10)=3; cd(5

30、,1)=2;cd(5,10)=2;cd(5,19)=2;cd(6,7)=6; cd(7,4)=6;cd(7,9)=6;cd(7,13)=6;cd(7,18)=6; cd(8,4)=6;cd(8,9)=6;cd(8,13)=6;cd(8,18)=6; cd(4,11)=if(cd(4,10)*(t(11)-t(10)#ge#m(4),0,cd(4,10); for(time(J)|J#ne#4 #and# J#ne#7 #and# J#ne#13 #and# J#ne#16 :cd(9,J)=20); cd(7,5)=if(cd(7,4)*(t(5)-t(4)#ge#m(7),0,cd(7,4)

31、; cd(7,14)=if(cd(7,13)*(t(14)-t(13)#ge#m(7),0,cd(7,13); for(time(J)|J#ge#2: cd(2,J)=if(t(J)-t(1)*cd(2,1)#ge#m(2),0,cd(2,1); for(time(J): total(j)=sum(work(I):cd(I,J); cd1=max(time(J):total(J); min=cd1+9; (cd1+9)*175=mcd; end （2）：退火法clear;clc;tic;D_org=D1;D0=D1;e=0.130;L=10000;at=0.999;T=1;for k=1:LD

32、1=D0;rand(state,sum(clock);num1=randperm(9);row=num1(5);R_row,R_col,val=find(D1(row,:) ;rand(state,sum(clock);len=length(R_col);num2=randperm(len);col_org=R_col(num2(ceil(length(R_col)/2);col=col_org;if col_org=1col=col_org;elsewhile D1(row,col_org-1)0col_org=col_org-1;col=col_org;if col_org=1col=1;

33、break;endendendrand(state,sum(clock);num3=randperm(D1(row,col);num4=num3(1);D1(row,col)=D1(row,col)+num4;n=work(row)/D1(row,col);n=ceil(n);tail=col+1;while D1(row,tail)0tail=tail+1;endfor i=0:n-1;D1(row,col+i)=D1(row,col);endfor i=n:tail;D1(row,col+i)=0;endtotal0=max(sum(D0);total1=max(sum(D1);df=to

34、tal1-total0;if dfrand(1)D0=D1;endendT=T*at;if Tebreak;endD1=D0(:,1:31);endtoc;D0=D0(:,1:31);delta=D0-D_org;disp(the number of worker: );disp(max(sum(D0);附录2油库大类岗位、需要人员及年工作量工作岗位类别工作岗位代号工作岗位名称工作岗位的年工作量备注（A）计量与质量检测管理A1月计量检测12次2人2天/次每月25日开始A2半年计量检测2次2人2天/次每年6、12月1日开始A3储油质量检验6次2人15天/次每单月15日开始A4全面检测2次2人2.

35、5天/次每年1、7月10日开始A5收油计量检测6次2人0.5天/次每双月10日A6收油质量检验6次2人0.5天/次每双月10日（B）收发油料管理B1收油操作6次16人0.5天/次每双月10日B2收油操作6次6人0.5天/次每双月10日B3收油操作6次1人1天/次每双月10日B4零发油操作120次5人0.5天/次平均每年120次B5零发油操作120次1人0.5天/次平均每年120次B6零发油操作120次1人0.5天/次平均每年120次（C）设备维护与维修管理C1收油设备维护值班6次2人1天/次每双月10日C2零发油设备维护值班120次2人0.5天/次平均每年120次C3常规维护保养1次12人30

36、天/次每年12月1日开始C4设备日常维护保养104次12人1天/次每周一、周五进行C5常用设备维护保养12次3人4天/次每月15日开始C6安全设备维护保养6次2人4天/次每单月20日开始C7配电设备维护保养26次2人0.5天/次每年双周的周三进行（D）安全保障管理D1收油消防值班6次6人1天/次每双月10日D2零发油消防值班120次6人0.5天/次平均每年120次D3消防设施维护52次6人1.5天/次每周一开始D4消防车辆维护52次6人0.5天/次每周一进行D5日常安全检查与维护183次20人1天/次每两天次D6油库环境保养与卫生241次1人1天/次除节假日外每天次（E）服务保障管理E1伙食服务保障员额每15人编制人全专职E2医疗服务保障2人/天全专职E3车辆服务保障4人/天全专职E4警卫服务保障8人/天全专职说明：（1）年工作量格式：次人天/次，表示该项任务每年要执行次，“至少需要人”同时执行，每次人需要连续工作z天。（）有些任务人数越多