工厂生产计划模型介绍(正式版)

1、理工大学暑期数学建模强化训练专题三工厂生产计划模型学员： 曹阳许佳利倪迪杭学院：通信工程学院时间：2010.08.18工厂生产计划模型摘要本文根据问题的条件和要求，综合考虑工厂获得最大收益所需的各种条件， 建立了单目标线性规划模型，并通过 LINGO软件的求解，给工厂生产提供了可 行有效的生产计划。对于问题一，针对每月都有不同的机床需要停工维修，且不同的产品所需的 机床工作时间、市场的容量以及市场价格都有所不同的条件， 采用线性规划结合 矩阵的方法，对工厂生产进行了模拟，并得到了获得最大收益时的可行方案。对于问题二，将市场价格的某种变化和引入新机床分开考虑。研究市场价格 的变化时，将新的产品收

2、益矩阵P'替代产品的原始收益矩阵P即可得到获得最大 收益的可行方案；研究引入新机床时，本文采用逆向思维分析确定引入何种新机 床使得收益最大，并给出了收益最大时的生产方案。对于问题三，建立了一个维修矩阵 Y，表示各月所要维修的机床及数量，通 过综合考虑各月获得最大收益的条件，运用线性规划的思想在问题一的模型的基 础上求出了最佳维修方案，进而确定了最佳生产计划，使得总收益最大。关键词：单目标线性规划，矩阵，生产计划，额外奖励，影响力一、问题的提出某厂拥有4台磨床、2台立式钻床、3台卧式钻床、1台镗床和1台刨床， 用以生产7种产品，记作P1至P7。工厂收益规定为产品售价减去原材料费用之 剩余

3、。每种产品单件的收益及所需各机床的加工工时（以小时计）列于下表：产品PlP2巳P4P5P6P7收益106841193磨床0.50.7000.30.20.5立式钻床0.10.200.300.6 :0 H卧式钻床0.200.80000.6镗床0.050.0300.070.10 :0.08 1刨床000.100.0500.05本月（一月）和随后的5个月中，下列机床停工维修:一月二月三月四月五月六月磨床1台卧式钻床2台镗床1台立式钻床1台磨床1台 立式钻床1台卧式钻床1台刨床1台各种产品各月份的市场容量如下表:产品RP2P3P4P5P6P7一月5001000300300800200100二月60050

4、02000400300150三月30060000500P 400100四月2003004005002000100五月01005001001000r 3000六月500500100300110050060每种产品存货最多可到100件，存费为每件每月0.5元。现无存货，要求到6月 底每种产品有存货50件。工厂每周工作6天，每天2班，每班8小时。不需要 考虑排队等待加工的问题。1、为使收益最大，工厂应如何安排各月份各种产品的生产？2、研究市场价格的某种变化及引入新机床对计划和收益的影响。3、若各机床的停工维修时间不作预先规定，而是选择最合适的月份维修。 除磨床外，每台机床在这6个月中的一个月必须停工

5、维修；6个月中4台磨床有 2台需要维修。扩展工厂计划模型，使得可以对灵活安排机床维修时间作出决策。 停工时间的这种灵活性价值如何？二、问题的分析对于问题一，根据题目要求，通过对各月每样产品的生产量进行合理的安排， 使得销售总收益能够达到最大。而收益受到机床的数目、工作时间以及产品的市 场容量和存货量等因素限制，我们可以运用线性规划的方法建立单目标线性规划 模型，结合LINGO软件求出最优生产计划，从而得到最大收益。对于问题二，题目要求模拟出市场价格变化下或者增加新机床时，应如何调 整生产计划以获得最大收益。我们可以对问题一所得的数据进行分析， 分别研究 市场价格变化和机床数量变化对收益的影响。

6、市场价格的变动，实际是改变了单件产品的收益，产品的价格在一个月内的 波动一般不大，假设每个月单件产品收益固定，那么就可以建立一个新的时间一产品收益矩阵，将这个矩阵替换问题一的收益矩阵，即可运用问题一的方法算 出最优的生产方案。增加新机床，一方面会使得工厂产品的生产量增加，但另一方面又使得工厂 的生产成本（新机床的成本）增加，从而影响收益。我们可以逆向思维，在问题 一的模型的基础上，逐一增加各种机床，从而可以得到逐一增加各种机床的最大 收益，通过判断增加的收益与机床成本的关系来确定增加何种机床使得工厂的总 收益最大。对于问题三，问题三较问题一来说机床停工维修时间的自由度变大了，这使得工厂能够更好

7、的根据市场需求来协调各产品每月的生产及存储计划，以达到最大收益。建立一个6 5的维修矩阵，代表各月需要维修的各种机床的数量， 然后 建立一个6 5的表示机床数量的矩阵，在解出最优维修矩阵后将后者减去前者便 可以得到各个月参加生产的机床的数量， 进而运用问题一的方法即可求出此时的 最优生产方案。三、问题的假设1、假设每月只有四周，即每月都只有二十四个工作日；2、不停工维修的机床在工作时间内可以同时工作，即不需要考虑排队等待 加工；3、产品的市场容量不受产品价格和机床数量等因素变化影响；4、市场价格的改变不影响产品的成本，即产品成本是固定不变的；5、认为新增加的机床的性能较好，在这 6个月内不需要

8、维修。四、模型的建立及求解1、符号说明X（i =12川6 j =1,2,川,7）第i个月生产Pj产品的件数aki（k =1,2,川,5）第i个月正常使用第k种机床的台数（k =1,2,“|,5分别表示磨床、立式钻床、卧式钻床、镗床和刨床）tkj单件Pj产品需由第k种机床加工的时间tPj工厂每个月的工作时间单件Pj产品的收益（profit）CjSj第i个月Pj产品的市场容量（capacity）第i个月Pj产品的存货量（stock）n第i个月Pj产品的销售量Zs6个月产品的总存费Zp6个月产品的总毛收益Z6个月产品的总纯收益erk第r次增加第k种机床得到的额外奖励lrk第r次增加第k种机床增加的成

9、本rk第r次增加第k种机床的影响力yki第i个月维修第k种机床的台数bki第i个月该工厂拥有第k种机床的台数2、模型建立2. 1工厂合理安排各月份各种产品的生产使收益最大的模型2.1.1模型的分析为了使收益最大，工厂在总的工作时间内一要尽可能多的生产产品，而且应该尽可能的生产效益高的产品，但是产品的生产数量还需考虑产品的市场容量； 二要尽可能少的存放货物，来减少存费，但是又必须要有一定量的货存量来满足 题意要求。我们可以用单位时间内生产的各种产品所获得的收益的多少来衡量其 效益的高低。通过计算，得到7种产品的效益如表一所示：表一各种产品的效益表产品RP2P3P4P5P6P7效益 （元/小时）1

10、1. 766.458.8910.8124.4411.252.44按照常理，我们预测按工厂获得最大收益的计划生产产品，效益较高的产品的收益应该也是较大的。2.1.2模型的建立与求解根据题目的已知条件，我们可以得到第i个月正常使用第k种机床的台数的 矩阵A为：i'z344434、222112A =(aki k =313332110111J11110单件Pj产品需由第k种机床加工的时间tkj的矩阵T为:0.50.7000.30.20.5 '0.10.200.300.60T =厲爲=0.200.80000.6 。0.050.0300.070.100.08<000.100.0500

11、.05 丿单件Pj产品的收益Pj的矩阵P为：P =(Pj >17= (1068 411 93 )。第i个月Pj产品的市场容量Cj的矩阵C为:5001000300300800200100"600500200040030015030060000500400100C = (Cj 1 =。' J怡72003004005002000100010050010010003000Q00500100300110050060第i个月P产品的存货量满足以下关系式：0乞Sj乞100 （即每月每种产品存 货最多可到100件），且®=50 （j=1,2,川,7）（即到6月底每种产品都有存

12、货 50 件）。从而我们可以推出，第i个月Pj产品的销售量nj为m =sM+Xj -q，且需 满足第i个月Pj产品的销售量不能超过第i个月Pj产品的市场容量，即nj < Cj。设工厂每个月的工作时间为t，则单种机床在第i个月的工作时间需满足关 系式tkjXj -at岂0，即每个月各种机床生产每一种产品的所花的时间不应超过该 月该机床的总工作时间。6个月产品的总存费Zs为：76Zs 工二 0.5Sij。y j 二6个月产品的总毛收益Zp为：76Z p 二、' Pin j。i =1 jm6个月产品的总纯收益Z为：Z = Zp -Zs。综上所述，得到问题一的单目标线性规划模型：6767

13、maxZ 八、 PjQ 八0.5ji =1 j di =1 j T0<j 兰 100,i =1,2,川,6, j =1,2川,7,© Eq,i=1,2，川,6, j=1,2,IU,7,s.t “tkjXj ant 兰0,i =1,2，川,6, j =1,2，川,7, k=1,2|,5, 山=丹4）% Sj,i =1,211,6, j =1,2,111,7, gj =50,j =1,2,ill,7.使用LINGO软件编程（程序见附件1）对该单目标线性规划模型进行求解， 我们用EXCEL软件将得到的结果（各月份的产量、存货量和销售量）整理如下 列表格所示：表二各种产品各月份的产量表

14、月份'PlP2P3P4P5P6P7一月5008883833008002000二月7006001170500300250三月000004000四月2004004006003000200五月00600003000六月550550035005500总量195024381500125016001750450表三各种产品的各月份销售量表销售产 月份'PP2P3F4P5P6P7一月500888300 :3008002000二月6005002000400300150三月1001000 10100400100四月2003004005002000100五月0100500100503000六月50

15、050050 :300050050总量190023881450120015501700400表四各种产品的各月份存货量表月份、FP2F3P4F5P6F7一月T0r 083 :0000二月100100001000100三月0000000四月0r 10001001000100 :五月001000500100六月50505050505050总量150r 250233 115030050350 1存费（元）75125116.57515025175上述三个表格中的数据，即为收益最大情况下，工厂安排各种机床生产各种 产品的计划。工厂按此计划生产六个月，获得的总毛收益Zp =83278元，总存费为Zs =7

16、41.5元，总收益Z =82536.5元。这就是用我们建立的单目标线性规划模 型计算的最大收益。2.1.3模型的结果分析（1）由表三可以看出，六月份各种产品的存货量都是50，正符合题中的要 求。（2）对比这三个表格，可以看出各种产品各月份的销售量正是该产品该月 份的产量与该产品上个月份的存货量之和减去该产品该月份的存货量。（3）表二中多个产品的产量为0,正是由各个月份中的机床停工维修造成 的，但是产量为0的时候，产品的销售量并不一定为 0,这是由于该产品在上一 个月份中有存货。表三中多个产品的销售量为 0，则是由该产品在该月份的市场 容量为0造成的。（4）对比分析各产品的收益与该产品的效益，两

17、者基本上或是在趋势上是一致的，即效益高的产品收益一般也是比较大的。 这说明模型建立之前我们的预 测是合理的。两者的对比图如图一所示。由于产品的收益数值较大，该图中收益 的曲线是在收益值缩小1000倍的时候作出的，这并不影响两者趋势的分析。2.2.1市场价格的变化对计划和收益的影响市场价格的变动，实际是改变了单件产品的收益。根据最大收益的计算公式:6767maxZ 八 '二二 0.5sij，i咼j咼y j壬市场价格的变化不会导致产品销售量和存货量的改变，若所有产品的市场价 格都增加，或者部分产品的市场价格增加，其他产品的市场价格不变，必然会使 得工厂的最大收益增加；而若所有产品的市场价格

18、都降低，或者部分产品的市场 价格降低，其他产品的市场价格不变，必然导致工厂的最大收益降低。当部分产 品的市场价格增加，其他产品的市场价格降低的时候，工厂的最大收益有可能增 加，也有可能降低。这个时候，我们可以利用问题一的模型， 用单件Pj产品的收 益矩阵p'替代问题一的模型中的收益矩阵 P，计算出产品市场价格改变后工厂的 最大收益，从而来研究市场价格的变化对计划和收益的影响。2.2.2引入新机床对计划和收益的影响分析问题一的数据我们可以得到，有些产品的销售量并没有达到市场容量,也就是说，如果增加新机床来增强生产能力是可以增加收益的。但究竟应该引入 何种机床，应该引入多少数量的机床，才能

19、使得收益最大呢？下面我们来分析研 究。引入新机床，必然带来生产成本的增加(这里增加的成本指的是购买机床的 成本，产品的成本是不变的)。由于购买机床的成本未知，这使得工厂的最大利 益不易求得。因此，我们用逆向思维分析的方法来考虑引入新机床对计划和收益 的影响。在问题一的模型的基础上，增加一台某种机床，然后用LINGO软件计算出其最大收益，更换机床种类便可得到一组最大收益数据。增加一台某种机床后得 到最大收益与问题一的模型计算出的最大收益相减得到的最大利益之差，我们称之为额外奖励。我们在增加机床数目的时候，每次都是增加一台。设第r次增加第k种机床得到的额外奖励为erk。eok即表示不增加机床的额外

20、奖励，显然 eok =0。设第r次增加第k种机床增加的成本为h。定义额外奖励与第k种机床成本的比值为影响力，记为frk，贝Ulrk分析得到，影响力frk具有下面两个性质：(1) 当frk 1时表示购买该机床获得的额外奖励大于成本，能使最大收益增加；当frk <1时表示购买该机床获得的额外奖励小于或等于成本，不值得购买。(2) frk越大，影响力越大，越值得购买。根据这两条性质，分析第r次增加机床的种类。显然，应该增加使得影响力 frk取得最大值的机床。假设各种机床的成本都相同，且较小。利用问题一的模型计算出第r次增加机床的方案如表五所示。表五第r次增加机床的方案类别磨床立式钻床卧式钻床镗

21、床刨床结论r 一1最大收益83508.582536.58257810335395791.5购买镗床额外奖励972041.520816.513255r 一2最大收益104325103353103353103353115483购买刨床额外奖励97200012130r 一3最大收益116455115483115483115483115483购买磨床额外奖励9720000r =4最大收益116455116455116455116455116455无需购买额外奖励00000这表明，在假设的基础上，第一次增加某种机床，镗床的额外奖励大于其它机床，说明镗床的影响力最大，所以第一次应该购买镗床。第二次增加某种

22、机床， 刨床的额外奖励最大，所以第二次购买刨床。第三次增加某种机床，磨床的额外 奖励最大，所以第三次购买磨床。第四次增加某种机床，各种机床的额外奖励都 为0，说明第三次购买机床后该工厂已经已经具有生产满足市场需求的产品的能 力，不需要再增加新机床。即工厂在原有机床的基础上，只需增加一台镗床、一 台刨床和一台磨床就能获得最大利益。 增加一台镗床、一台刨床和一台磨床后工 厂收益最大的生产计划见附录 2。而在实际情况中，机床的成本一般不会很小，而在第三次增加某种机床的时 候，额外奖励已较小（972）,若机床的成本大于972，则其影响力f“加，就不 需要购买机床。2. 3工厂计划的扩展模型2.3.1模

23、型的建立和求解若各机床的停工维修时间不作预先规定， 而是选择最合适的月份维修。在问题一的模型的基础上，我们引入两个新的矩阵 B和Y。设第i个月该工厂拥有第k种机床的台数为bki，其矩阵为：'444444、222222B =(bki k 二33333311111111111设第i个月维修第k种机床的台,数为Yki，其,矩阵为:丫二 yki 56 0则第i个月正常使用第k种机床的台数aki二匕总-丫战。根据问题三的要求，除磨床外，每台机床在这6个月中的一个月必须停工维 修（这里假设只维修一次），所以有：6- yki 二九，i Wk6个月中4台磨床有2台需要维修，有：6二 yki = 2。i

24、 dk d其他约束条件与问题一的模型相同，因此可以得到问题三的单目标线性规划 模型:maxZ 八、pj©'0.5可i 4 j =1i 4 j =1'6送 yki =bki,i 4谆0"j 兰 100,i =1,2,川,6, j =121 山7,% "j,i=1,2,川,6, j=12|(),7,s.t 社必ant 兰0,i=1,2,川6 j=1,2,川,7, k=1,2,山,5,nij = S(i 二 j + Xjj sj i = 1,2,川,6, j = 1,2,1 () ,7,S6j =50, j =1,2,川,7,6送 yki =2.i 二k

25、4同样使用LINGO软件编程（具体程序见附录3）对该单目标线性规划模型 进行求解，我们将计算得到的结果（各月份的产量、存货量、销售量和维修方案） 整理如下列表格所示：表六各种产品各月份的产量表月PP2P3P4P5P6P7一月5001000300300800200100二月6005002000400300150三月400700100100600400200四月0000000五月010050010010003000六月5505501503501150550110总量20502850125085039501750560表七各种产品的各月份销售量表销产品 月份、PP2P3F4P5P6P7一月50010

26、00300300800200100二月6005002000400300150三月30060000500r 400100四月1001001001001000100五月010050010010003000六月5005001003001100:50060总量20002800120080039001700510表八各种产品的各月份存货量表1货产量'品 月鳥PlP2P3P4P5P6P7-月h0r 00 :0000 :二月0000000三月1001001001001000100四月0r 000000 :五月0000000六月50505050505050总量150150150 11501505015

27、0 1存费（元）75757575752575表九各种机床各月份的维修方案月份、12345一月00100二月01000三月00 :000四月00011五月21 1100六月00100上述四个表格中的数据，即为收益最大情况下，工厂安排各种机床生产各种 产品和安排机床进行维修的计划。工厂按此计划生产六个月，获得的总毛收益Zp =109330元，总存费为Zs =475元，总收益Z =108855元。这就是用我们建立 的扩展模型计算的最大收益。停工时间的灵活性价值我们用扩展模型计算的最大 收益与问题一的模型计算的最大收益之差来体现，其值为25305元。2.3.2模型的结果分析（1）比较扩展模型和问题一的

28、模型的结果，我们可以得到扩展模型的各种产品的销售量都比问题一的模型的销售量大，而扩展模型的各种产品的存货量都小于等于问题一的模型的存货量， 这说明在灵活调节各机床的维修时间之后，工厂的收益大大增加了，这正是停工时间的灵活性价值所在。（2）模型扩展后，产品P5的产量和销售量大大增加（均增加了 2350台）， 其他产品的产量和销售量变化不大，这说明调节各机床的维修时间来增加效益高 的产品（艮的效益为24.44元/小时）的产量能使得收益最大化，这也是符合常 理的。五、模型的评价模型的优点：1、本文主要利用矩阵和动态规划的相关知识建立单目标线性规划模型求解， 最终为满足工厂的不同要求给出了最佳生产方案

29、，具有很好的现实指导意义；2、本模型使用矩阵来定义各种影响因素，不仅避免了大量方程的赘写，而 且增加了模型的灵活性，便于对各个参数进行拓展研究，问题二中对市场价格的 研究和问题三的解决就是很好的例子；模型的缺点：本文在解决问题的时候，模型较为单一，并且没有用很好的检验方法来检验 最大收益值。六、模型的进一步探讨在问题二的求解中，在研究引入新机床对计划和收益的影响的时候，我们假设各种机床的成本相同，从而将影响力直接体现在额外奖励上。而在实际情况下， 各种机床的成本是不相同的。因此，在模型的改进中，可以经过市场调查的方式， 得到各种机床的成本，从而精确的计算出其影响力。这样就能更加准确地制定机 床

30、购买方案，更具实际指导意义。参考文献1.袁新生 邵大宏等 LINGO和excel在数学建模中的应用科学出版 社2007年1月出版附录1MODEL:sets:pii/1.7/:p;mon th/1.7/:m;chua n/1.5/:c;time(chua n,pii):t;mon c(mo nth,chua n):mc;ron g(mo nth,pii):r,s, n,x;en dsetsdata:p=1068411 93;t=0.50.7000.30.20.50.1 :0.200.3r 00.6r 0 :0.200.80000.60.050.0300.07r 0.10r 0.08000.100

31、.0500.05mc=00000323p14211142301413114130142210r=000000050010003003008002001006005002000400300150300600005004001002003004005002000100010050010010003000 :500500100300110050060en ddatatz=24*16;for(rong(i,j)|i#gt#1:n(i,j)<=r(i,j);!销售量小于等于市场容量；for(pii(j):s(1,j)=0);for(ro ng(i,j):g in (x(i,j);gi n(s(i,

32、j);for(ro ng(i,j)|i#gt#1: n( i,j)=s(i-1,j)+x(i,j)-s(i,j);!销售量=上月存储量+本月生产量-本月存量；for(rong(i,j):s(i,j)<=100);!存量小于等于 100 件；for(ro ng(i,j)|i#gt#1:s(7,j)=50);zs=0.5*sum(rong(i,j)|i#gt#1:s(i,j);!计算存储费用；zp=sum(ro ng(i,j)|i#gt#1: n( i,j)*p(j);max=zp-zs;!净禾U润;for(mo nth(i)|i#gt#1:for(chua n(k):sum(pii(j):

33、t(k,j)*x(i,j)-mc(i,k)*tz<=0);!对每月机床数量的约束；end附录2表十各种产品各月份的产量表口月份、RP2P3P4P5P6P7一月500 11000300r 300800200100二月6005002000400300150三月300 16000r 0500400100四月2003004005002000100五月010050010010003000六月550 :550150P 3501150550110总量215030501550125040501750560表一各种产品的各月份销售量表销售产月份'、PP2P3P4P5P6P7一月5001000300P 300800200100二月6005002000400300150三月300 16000r 0500400100四月2003004005002000100五月010050010010003000六月500 500100r 300110050060总量210030001500120040001700510表十二各种产品的各月份存货量表1货产月P1P2P3P4P5P6P7一月T000r 0000 :二月0