企业最优生产的数学模型

1、企业最优生产的数学模型第九组：张乐 康倩妮 罗少梅 (西安航空学院，西安 )摘要本文针对企业及工厂应该怎样合理安排生产计划而获得最大利润做了简单分析，主要用于解决企业及工厂在各种互相矛盾，互相排斥的约束条件下如何安排生产获取最大利润，建立了生产量对利润影响的线性规划模型。对于问题一，根据对影响利润的因素的初步分析，综合得出其主要因素有：每种产品的单件利润、生产单位各种产品所需的有关设备台时、生产量、最大需求量、库存量、每月的工作时间、设备维修。综合考虑多种因素，利用线性规划来建立模型解决问题，即将每月各种产品的最大需求量、一月初无库存、任何时候每种产品的存储量均不能超过100件、六月末各种产品

2、各储存50件作为约束条件，最大利润作为目标函数，利用lingo11.0软件求解，得出最大利润为：93.71518万元。对于问题二，要求重新安排维修，并以最大利润作为前提，类比于问题一，并在问题一模型的基础上，添加，分别为第种设备在第个月工作的台数和第种设备在第个月维修的台数。并定义为在不进行维修的情况下工作的台数，则=+；表示第种设备在第月维修的台数等于每种设备可以维修的台数。，利用lingo11.0软件求解，得出结论最大利润为108.8550万元。1月2月3月4月5月6月磨床2立钻11水平钻12镗床1刨床1关键词：线性规划、lingo软件、最大利润问题的提出每个企业都希望在成本最低，工作时间

3、最短的条件下获得最大利润，但各种约束条件总是互相制约，这就需要我们在考虑到实际情况时，酌情筛减。已知某企业要生产7种产品，以，来表示，并给出了每种产品的单件利润，生产单件每种产品的设备所耗费的时间及每种产品在各个月的最大需求量。产品当月销售不了的每件每月存储费为5元，且任何时候每种产品的存储量均不能超过100件。一月初无库存，要求六月末各种产品各存储50件，并且每月均有设备参与维修：一月维修1台磨床，二月维修2台水平钻，三月维修1台镗床，四月维修1台立钻，五月维修1台磨床和1台立钻，六月维修1台刨床和1台水平钻。考虑问题如下：问题一：欲使总利润最大，应如何计划生产？问题二：要在上半年进行设备维

4、修（4台磨床只需安排2台在上半年维修），在使得利润最大的基础上设计一个可行的维修方案。模型的假设1. 除去在维修的设备外，其它设备均能正常运转；2. 每种产品的单件利润随市场波动不大；3. 题目中所给的数据均具有代表性，能够说明问题；4. 该企业按正常工作时间工作，无加班情况；5. 生产单位各种产品所需的有关设备台时不变。符号的说明：第个月第种产品的生产量；：第种产品的单品利润；：第个月第种产品的库存量； ：第个月第种产品的最大需求量；：第种设备生产第种单位产品所需台时；：第种设备在第个月工作的台数；：第种设备在第个月维修的台数；：第种设备在第个月工作的台数（不进行维修的情况下） ；： 第种设

5、备参与维修的总台数；问题的分析问题一的分析： 基于题目的要求，1月初每件产品的库存量为0，因此假设总共有7个月且1月库存量、需求量和生产量均为0，这样就使得原来的1月对应于假设中的2月，以此类推。 根据假设，每种产品在第7个月时的库存量均为50；第个月第种产品的库存量小于等于100；第()个月第种产品的生产量与上一个月的库存量之和减去这个月的库存量小于等于第月第种产品的最大需求量。第种设备在第个月生产第种产品数与第台设备生产第种单位产品所用时间的乘积小于等于第种设备在第个月工作的台数与每月工作的总时间的乘积。问题二的分析： 要求重新安排维修，使得利润最大，类比于问题一，在问题一的基础上设定未知

6、量，分别为第种设备在第个月工作的台数和第种设备在第个月维修的台数。并定义为在不进行维修的情况下工作的台数，则=+；且第种设备在第月维修的台数等于每种设备可以维修的台数。模型的建立问题一：建立最大利润规划模型（模型一） 最大利润与生产量，单件利润，生产时间，库存量等都有联系，基于这几种约束条件，我们由题目中所给条件建立了线性规划模型： 利用lingo11.0软件即可得出获得最大利润的生产方案。（程序见附录1）0问题二：建立获得最大利润的维修方案（模型二） 利用模型一，定义其已知变量为所求变量，添加合理约束条件，建立模型如下： 利用lingo11.0得出获得最大利润的维修方案。（程序见附录2）模型

7、的求解问题一的求解：根据题中所给数据，利用模型一求解得出结论如下：按下表计划生产可获得最大利润为：93.7151万元每种产品每个月的库存量：库存量产品月份1月008300002月1001000010001003月00000004月00000005月00100010001006月50505050505050每种产品每个月的生产量：生产品产量月份1月50088938330080020002月70060011805003002503月0000040004月20030040050020001005月010060010011003001006月550550035005500问题二的求解利用模型二求解获得

8、最大利润的维修方案如下：按下表安排维修可获得最大利润为：108.8550万元。设备维修安排：维修月台数份设备1月2月3月4月5月6月磨床2立钻11水平钻12镗床1刨床1模型的评价与推广 企业获得利润往往和生产量和销售量有着密切的联系，因本文只给出了生产量和每种产品的大约利润，忽略了单件产品利润随市场波动的影响，并且没有考虑市场需求量的变化及多种因素，这样使得本文的实际应用能力被大大局限。为使本文便于推广，则应该调查单件产品利润变化的大致范围，市场需求量的变化范围，应用于题中，在模型一二种稍做改进，利用lingo编程求解。参考文献1胡红亮 赵芳玲学建模与竞赛辅导第一版 西北大学出版社2010年6

9、月2朱德通 李炳钊最优化模型与实验第一版 同济大学出版社2003年4月3杨桂元 李天胜 徐军 数学模型应用实例第一版 合肥工业大学出版社 2007年6月4 陈根雷 黄继晨 龚佳丽 数学建模竞赛论文 2012/7/20/view/f97ed705cc884.html5 谢金星 lingo教程 2012/7/19 /view/4d2d2c2d915f804d2b16c1f0.html附录1：model:sets:yue/1.7/:t;cp/1.7/:a;sb/1.5/:k;taishi(sb,cp):d;jiqi(

10、sb,yue):b;www(yue,cp):c,v,x;endsetsdata:a=100 60 80 40 110 90 30;b=0 3 4 4 4 3 4 0 2 2 2 1 1 2 0 3 1 3 3 3 2 0 1 1 0 1 1 1 0 1 1 1 1 1 0;d=0.5 0.7 0 0 0.3 0.2 0.5 0.1 0.2 0 0.3 0 0.6 0 0.2 0 0.8 0 0 0 0.6 0.05 0.03 0 0.07 0.1 0 0.03 0 0 0.01 0 0.05 0 0.05;c=0 0 0 0 0 0 0 500 1000 300 300 800 200 100

11、 600 500 200 0 400 300 150 300 600 0 0 500 400 100 200 300 400 500 200 0 100 0 100 500 100 1000 300 0 500 500 100 300 1100 500 60;enddatamax=sum(www:x*a-5*v)-sum(cp:a*50); for(www:v=100); for(cp(j):v(7,j)=50); for(cp(j):v(1,j)=0);for(www(i,j)|i#gt#1: x(i,j)+v(i-1,j)-v(i,j)=c);for(sb(w): for(yue(i):s

12、um(cp(j):x(i,j)*d(w,j)=b*384);end执行结果：Global optimal solution found. Objective value: .8 Infeasibilities: 0. Total solver iterations: 47 V( 1, 1) 0. 0. V( 1, 2) 0. 0. V( 1, 3) 0. 0. V( 1, 4) 0. 0. V( 1, 5) 0. 0. V( 1, 6) 0. 0. V( 1, 7) 0. 0. V( 2, 1) 0. 46.60714 V( 2, 2) 0. 65.00000 V( 2, 3) 82.5000

13、0 0. V( 2, 4) 0. 5. V( 2, 5) 0. 30.71429 V( 2, 6) 0. 22.14286 V( 2, 7) 0. 31.25000 V( 3, 1) 100.0000 0. V( 3, 2) 100.0000 0. V( 3, 3) 0. 10.00000 V( 3, 4) 0. 0. V( 3, 5) 100.0000 0. V( 3, 6) 0. 5. V( 3, 7) 100.0000 0. V( 4, 1) 0. 105.0000 V( 4, 2) 0. 65.00000 V( 4, 3) 0. 5. V( 4, 4) 0. 10.00000 V( 4

14、, 5) 0. 115.0000 V( 4, 6) 0. 5. V( 4, 7) 0. 35.00000 V( 5, 1) 0. 5. V( 5, 2) 0. 5. V( 5, 3) 0. 5. V( 5, 4) 0. 5. V( 5, 5) 0. 5. V( 5, 6) 0. 5. V( 5, 7) 0. 5. V( 6, 1) 0. 5. V( 6, 2) 0. 5. V( 6, 3) 100.0000 0. V( 6, 4) 0. 5. V( 6, 5) 100.0000 0. V( 6, 6) 0. 5. V( 6, 7) 100.0000 0. V( 7, 1) 50.00000 0

15、. V( 7, 2) 50.00000 0. V( 7, 3) 50.00000 0. V( 7, 4) 50.00000 0. V( 7, 5) 50.00000 0. V( 7, 6) 50.00000 0. V( 7, 7) 50.00000 0. X( 1, 1) 0. 153.5714 X( 1, 2) 0. 272.1429 X( 1, 3) 0. 0. X( 1, 4) 0. 0. X( 1, 5) 0. 482.1429 X( 1, 6) 0. 0. X( 1, 7) 0. 612.1429 X( 2, 1) 500.0000 0. X( 2, 2) 888.5714 0. X

16、( 2, 3) 382.5000 0. X( 2, 4) 300.0000 0. X( 2, 5) 800.0000 0. X( 2, 6) 200.0000 0. X( 2, 7) 0. 12.85714 X( 3, 1) 700.0000 0. X( 3, 2) 600.0000 0. X( 3, 3) 117.5000 0. X( 3, 4) 0. 0. X( 3, 5) 500.0000 0. X( 3, 6) 300.0000 0. X( 3, 7) 250.0000 0. X( 4, 1) 0. 0. X( 4, 2) 0. 0. X( 4, 3) 0. 0. X( 4, 4) 0

17、. 135.0000 X( 4, 5) 0. 90.00000 X( 4, 6) 400.0000 0. X( 4, 7) 0. 30.00000 X( 5, 1) 200.0000 0. X( 5, 2) 300.0000 0. X( 5, 3) 400.0000 0. X( 5, 4) 500.0000 0. X( 5, 5) 200.0000 0. X( 5, 6) 0. 0. X( 5, 7) 100.0000 0. X( 6, 1) 0. 0. X( 6, 2) 100.0000 0. X( 6, 3) 600.0000 0. X( 6, 4) 100.0000 0. X( 6, 5

18、) 1100.000 0. X( 6, 6) 300.0000 0. X( 6, 7) 100.0000 0. X( 7, 1) 550.0000 0. X( 7, 2) 550.0000 0. X( 7, 3) 0. 0. X( 7, 4) 350.0000 0. X( 7, 5) 0. 290.0000 X( 7, 6) 550.0000 0. X( 7, 7) 0. 370.0000 附录2：model:sets:yue/1.7/:t;cp/1.7/:a;sb/1.5/:k,s;taishi(sb,cp):d;jiqi(sb,yue):b,z,p;www(yue,cp):c,v,x;en

19、dsetsdata:a=100 60 80 40 110 90 30;d=0.5 0.7 0 0 0.3 0.2 0.5 0.1 0.2 0 0.3 0 0.6 0 0.2 0 0.8 0 0 0 0.6 0.05 0.03 0 0.07 0.1 0 0.03 0 0 0.01 0 0.05 0 0.05;c=0 0 0 0 0 0 0 500 1000 300 300 800 200 100 600 500 200 0 400 300 150 300 600 0 0 500 400 100 200 300 400 500 200 0 100 0 100 500 100 1000 300 0

20、500 500 100 300 1100 500 60;s=2 2 3 1 1;p=0 4 4 4 4 4 4 0 2 2 2 2 2 2 0 3 3 3 3 3 3 0 1 1 1 1 1 1 0 1 1 1 1 1 1;enddatamax=sum(www:x*a-5*v)-sum(cp:a*50); for(www:v=100); for(cp(j):v(7,j)=50); for(cp(j):v(1,j)=0);for(www(i,j)|i#gt#1:x(i,j)+v(i-1,j)-v(i,j)=c);for(sb(w): for(yue(i):sum(cp(j):x(i,j)*d(w,j)=b(w,i)*384);FOR(sb(i): sum(yue(j):z(i,j)=s(i);for(jiqi:gin(b);for(jiqi:z+b=p);end执