《2022年数学建模——工厂计划模定义 .pdf》由会员分享，可在线阅读，更多相关《2022年数学建模——工厂计划模定义 .pdf（14页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、第 1 页 共 14 页YangQiuyan 第 1 页2013/4/17 数学建模工厂生产计划模型学院：数学与统计学院专业：信息与计算科学教师：郑小 洋姓名：杨秋 燕学号： 11001010128 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 1 页，共 14 页 - - - - - - - - - 第 2 页 共 14 页YangQiuyan 第 2 页2013/4/17 摘要本文以工厂所获得的总收益为研究对象，采用了线性规划的分析方法，通过求解不同产品的生产计划以及按计划生产
2、所获得的利润，解决了工厂为达到最大总收益的产品生产计划问题。在问题一的求解过程中，以每月每种产品的销售量和生产量为自变量，以工厂所获得的收益为目标函数，结合各种约束条件，建立了一个动态规划方程组，将各月份各种产品生产的最佳配置转化为动态规划方程组的求解问题，得到了最大收益为6.9256 万元。问题二在问题一的基础上考虑了市场价格的变化及引入新机床两个因素，为使模型简化，首先考虑市场价格的变化对计划和收益的影响。然后假定市场价格不变，利用Lingo软件， 模拟出引入新机床对计划和收益的影响。它是问题一的拓展，通过更改约束方程，利用模型一的计算程序，从而得到拓展模型的最优解。关键字：总收益销售量生
3、产量动态规划一、问题重述某厂拥有4 台磨床、 2 台立式钻床、 3 台卧式钻床、一台镗床和一台刨床，用以生产7种产品，记作P1 至 P7。工厂收益规定为产品售价减去原材料费用之剩余。每种产品单件的收益及所需各机床的加工工时（以小时计）列于下表：产品P1 P2 P3 P4 P5 P6 P7 收益10 6 8 4 11 9 3 磨0.5 0.70 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.08 刨0 0 0.01 0 0.05 0 0.05 本月（一月）和
4、随后的5 个月中，下列机床停工维修：一月磨床一台二月卧式钻床2 台三月镗床一台四月立式钻床一台五月磨床一台，立式钻床一台, 上台下六月刨床一台，卧式钻床一台各种产品各月份的市场容量如下表：产品P1 P2 P3 P4 P5 P6 P7 一月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 六月500 500 100 300 1100 500 60 每种产品存货最多可到100
5、 件。存费每件每月为0.5 。现在无存货。要求到6 月底每种产品有存货 50 件。工厂每周工作6 天，每天 2 班，每班 8 小时。不需要考虑排队等待加工的问题。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 2 页，共 14 页 - - - - - - - - - 第 3 页 共 14 页YangQiuyan 第 3 页2013/4/17 1、为使收益最大，工厂应如何安排各月份各种产品的产量？（考虑价格的某种变化及引入新机床对计划和收益的影响。注意，可假设每月仅有24 个工作日。）
6、2、在工厂计划问题中，各台机床的停工维修不是如问题( I ) 那样规定的月份，而是选择最合适的月份维修。除了磨床外，每台机床在这个月6 中的一个月必须停工维修；6 个月中 4台磨床只有2 台需要维修。 扩展工厂计划模型，以使可作上述灵活安排维修时间的决策。停工时间的这种灵活性价值如何？二、模型假设1.机器除维修时期外，其他时间均能正常工作，2.不需要考虑排队等待加工的问题，3.产品在储存过程中不发生质量问题，4. 该工厂生产的产品除存货外，在不超过市场容量的情况下均能售出，5. 市场价格的变化对产品原料和销售价格均有影响. 三、模型分析及建立模型1、符号说明2、模型的分析21 需解决的问题按照
7、题意，需要解决的问题是给出一个优化的产品加工方案，方案应包括每月每种机床的加工各种产品的数量，以使得工厂所获得的收益最大。2.2 达到的目标工厂的目标是获得最大利润，对于本题，产品生产的利润等于产品的收益乘以产品的符号意义kA k种设备的总数量kjA第 k 种设备第j 个月的可用数量iB第 i 种产品的单件收益ijB第 j 个月第 i 种产品的单件收益ijP第 j 个月第 i 种产品的生产量ijL第 j 个月第 i 种产品的存储量ikT第 k 种设备生产每个i 产品的用时ijV第 j 个月第 i 种产品的市场容量ijS第 j 个月第 i 种产品的销售量kjN第k个设备第j 月能用的数量名师资料
8、总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 3 页，共 14 页 - - - - - - - - - 第 4 页 共 14 页YangQiuyan 第 4 页2013/4/17 数量减去库存产品的存储费用。本题要达到的目标就是通过建立数学模型寻求收益的最大化。 2.3约束条件由于利润等于产品的总收益减去库存产品的存储费用，就是说影响利润有产品的数量、库存产品的存储费用量两个因素，利用其约束条件将转化为求解动态规划方程组的问题，现分别对它的约束条件作如下分析。（ 1）机床每月可工作时间的约
9、束. 工厂每周工作6 天，每天 2 班，每班 8 小时，所以每台机器每月（按30 计算）最多可以工作441 小时，即满足71411ikijkjiTPN；（2）各种产品销售量的约束. 每种产品每月的销售量由市场容量及该种产品的存货量决定，题中已给出各种产品各月份的市场容量表及最大存货量100 件，故销售量应小于等于市场容量与存储量之和即100ijL；(1)()ijijijSLP；（3）六月份每种产品的存货应为50 件 . 650iL；（4）机床停工维修时间的限制. 由于每个月份都有机床停工维修（见题表），经分析，当该种机器处于维修状态时，其可用台数kN会减少从而影响各种产品的生产数量，进而影响总
10、的销售额。3、模型的建立与求解3.1 工厂合理安排各月份各种产品的生产使收益最大的模型3.1.1模型的建立（1）目标函数题目的目标是寻求利润的最大化，而利润Y为六个月的收益的和, 即61jjYY；每月总收益jY由各品种产品的收益之和, 即ij71ijSBYi；库存产品的存储费用为j0.5 L可表示为7j10.50.5ijiLL；名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 4 页，共 14 页 - - - - - - - - - 第 5 页 共 14 页YangQiuyan 第 5
11、 页2013/4/17 目标函数可表示为max6717ji11i0.5()iijjiBSL；化简得max7iji611(0.5)iijjLBS；(2) 约束条件由以上分析可得到如下约束条件71411ikijkjiTPN. ijijSV. (1)()ijijijSLP. 100ijL. ；650iL. (1).ijijijijLLPS00iL目标函数与约束条件的优化模型即为问题一的模型。3.1.2模型的求解我们根据建立的线性规划模型，通过lingo软件编程（程序见附录1） ，得到了工产品加工计划的最大值为6.9256 万元。具体的生产安排计划见表1 产品月份1P2P3P4P5P6P7P一月生产量
12、 (个) 50 50 50 50 50 50 50 销售量 (个) 50 50 50 50 50 50 50 二月生产量 (个) 650 550 293 93 450 393 200 销售量 (个) 600 500 200 0 400 300 150 三生产量 (个) 0 0 0 0 0 394 0 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 5 页，共 14 页 - - - - - - - - - 第 6 页 共 14 页YangQiuyan 第 6 页2013/4/17 月销
13、售量 (个) 50 50 0 0 50 400 50 四月生产量 (个) 275 387 394 494 287 0 175 销售量 (个) 200 300 400 500 200 0 100 五月生产量 (个) 0 88 488 88 988 288 0 销售量 (个) 0 100 500 100 1000 300 0 六月生产量 (个) 475 475 0 275 0 475 0 销售量 (个) 500 500 25 300 25 500 25 表 1 问题一的产品加工计划3.1.3结果分析：加工计划与当月能用的机床、市场容量、 加工产品的时间及加工产品的单项收益有关，尽量在前期避免生产量
14、大于销售量，不然会增加存储费用从而增加成本，尽量增大单项收益高、加工时间短的产品的生产量3.2 问题二的模型建立与求解3.2.1 、 该问研究市场价格的某种变化及引入新机床对计划和收益的影响。为简化模型，我们先分析市场价格对计划和收益的影响，然后在此基础上，由小到大逐个增加五种机床的引进数量并分析总收益的变化趋势来反映其对引入新机床对计划和收益的影响。市场价格对单个产品的售价和原材料成本均有影响，不妨假设， 单个产品的收益月减少率为 x，则第 j 个月单个产品的利润ijB=1iB1(1)jx；其中1iB为第 1 个月第 i 个产品单件的收益。故可以建立总收益的目标7iji 1610.5)(ij
15、ijYmLBaSx；其它约束条件均不变。当 x 从 0.01 每次逐渐增加0.02 时，利用 lingo软件（程序见附录2） ，得到对应每个x 值时的最大总收益，如下表2 所示x 0.01 0.03 0.05 0.07 0.09 0.11 0.13 0.15 0.17 0.19 总收益 ( 元) 68553 67148 65743 64338 62933 61528 60123 58718 57313 55908 表 2 不同 x 值时的最大总收益名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - -
16、- - 第 6 页，共 14 页 - - - - - - - - - 第 7 页 共 14 页YangQiuyan 第 7 页2013/4/17 用 EXCEL软件绘图，得到图1：图 1：最大收益率与市场价格与增长率的关系从表中可以发现，最大总收益随x 的增大而基本呈线性的减小，且总收益为正。接着考虑市场价格不变，即单项产品收益不变时引入新机床对计划和收益的影响，因为磨床加工耗时比较长，所以考虑磨床的引进对总收益的影响。如首先讨论磨床的数量由4逐次增加至8 的过程中，总收益Y 的变化趋势来衡量引入磨床对总收益的影响。用lingo软件进行求解（程序见附录3） ，分析结果发现引进磨床没有对总收益产
17、生影响，总收益仍为 69256 元。再进行研究， 发现引进立式钻床、卧式钻床对总收益也无影响。但当单独引进一台镗床时总收益变为80425 元，比原来增加了11169 元。当单独引进一台刨床时，总收益变为 81786 元，较原来增加了12530 元。当引进一台镗床和一台刨床时，总收益为92955元，较原来增加了23699 元。总收益共这四种结果，在引进机床无影响。从引进机床的结果可以看出，原题中的机床维修对总收益并无很大影响。原因可能是总收益还受到市场容量以及存储量的限制，所以多引进机床并不能多生产。而引进镗床和刨床则会使生产方案发生变化，因为原有的这两台机床就很少，而且在某一月中还要进行维修。
18、3.2.2 （1）模型的改进由于机床的停工维修时间不作预先规定，而是选择最合适的月份维修，对各种产品每月的加工数量的限制改变了。针对这个情况，对于问题三的模型，仍采用问题一的模型，只不过对约束条件进行了修改。目标方程max7iji611(0.5)iijjLBS；（2）约束条件名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 7 页，共 14 页 - - - - - - - - - 第 8 页 共 14 页YangQiuyan 第 8 页2013/4/17 由以上分析，可得到如下约束条件
19、671111411(62)iijjiTPA. 6711411(66),ikijkjiTPA25.kijijSV. (1)()ijijijSLP. 。100ijL. 650iL. (1).ijijijijLLPS00iL四、模型优缺点分析1、模型的优点1. 在模型求解的时候，利用了专门求解规划问题的lingo数学软件， 求解速度很快， 而且结果准确； 2. 在模型的改进中， 采用了动态规划模型，考虑不同变量之间的关系，求解全局最优解。3. 对工厂的最大总收益进行了研究分析，认为增加一台镗床和一台刨床更能提高最大总收益。2、模型的缺点1. 在题目中，我们没有考虑产品每天的存储费用，只是简单地把存储
20、费用看作是月底的存储量与存储价格的积，这样不太全面，有失偏差；2. 模型较为单一，并且没有用很好的检验方法来检验最大收益值。五、模型推广本模型是一个典型的线性规划模型，用来求解最大或最小目标函数值问题。此类问题很多， 也有很多的推广应用价值。优化问题可以说是人们应用科学、工程设计、 商业贸易等领域中最常遇到的一类问题。这种用数学建模的方法来处理优化问题，即建立和求解所谓优化模型。虽然由于建模时要做适当的简化， 可能使得结果不一定完全可行或达到实际上的最优，但是它基于客观规律和数据，模型的建立与求解并不需要耗费太多时间。如果在建模的基础上再赋予其现实意义，就可以期望得到实际问题的一个比较圆满的回
21、答。六、参考文献名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 8 页，共 14 页 - - - - - - - - - 第 9 页 共 14 页YangQiuyan 第 9 页2013/4/17 1.袁新生邵大宏等LINGO和 excel在数学建模中的应用科学出版社 2007年 1月出版2 姜启源，谢金星，叶俊，数学模型（第三版），高等教育出版社，2003 年 8 月. 3 叶其孝，大学生数学建模竞赛辅导教材（三），湖南教育出版社，1998 年 5 月 . 七、附录附录 1.问题一
22、中的程序代码model: sets: m/1.6/:； ! 定义月份下标集；k/1.7/:g； ! 定义产品种类；b/1.5/:； ! 定义机床种类；vol(m,k):v,p,l,s； ! 定义容量矩阵 v, 生产方案矩阵 p, 储存矩阵 l 和销售矩阵 s；time(b,k):t； ! 定义时间矩阵 t ；endsets data: v=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
23、 500 500 100 300 1100 500 60；g=10 6 8 4 11 9 3；t=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.08 0 0 0.1 0 0.05 0 0.05；enddata max=sum(m(i):sum(k(j):g(j)*s(i,j)-0.5*sum(k(j):l(i,j)； ! 目标函数；for(k(j):t(1,j)*p(1,j)=411*3)；for(k(j):t(1,j)*p(2,j)=411*4)；for(k
24、(j):t(1,j)*p(3,j)=411*4)；for(k(j):t(1,j)*p(4,j)=411*4)；for(k(j):t(1,j)*p(5,j)=411*3)；for(k(j):t(1,j)*p(6,j)=411*4)；for(k(j):t(2,j)*p(1,j)=411*2)；for(k(j):t(2,j)*p(2,j)=411*2)；for(k(j):t(2,j)*p(3,j)=411*2)；for(k(j):t(2,j)*p(4,j)=411*1)；for(k(j):t(2,j)*p(5,j)=411*1)；名师资料总结 - - -精品资料欢迎下载 - - - - - - -
25、- - - - - - - - - - - 名师精心整理 - - - - - - - 第 9 页，共 14 页 - - - - - - - - - 第 10 页 共 14 页YangQiuyan 第 10 页2013/4/17 for(k(j):t(2,j)*p(6,j)=411*2)；for(k(j):t(3,j)*p(1,j)=411*3)；for(k(j):t(3,j)*p(2,j)=411*1)；for(k(j):t(3,j)*p(3,j)=411*3)；for(k(j):t(3,j)*p(4,j)=411*3)；for(k(j):t(3,j)*p(5,j)=411*3)；for(k(
26、j):t(3,j)*p(6,j)=411*2)；for(k(j):t(4,j)*p(1,j)=411*1)；for(k(j):t(4,j)*p(2,j)=411*1)；for(k(j):t(4,j)*p(3,j)=411*0)；for(k(j):t(4,j)*p(4,j)=411*1)；for(k(j):t(4,j)*p(5,j)=411*1)；for(k(j):t(4,j)*p(6,j)=411*1)；for(k(j):t(5,j)*p(1,j)=411*1)；for(k(j):t(5,j)*p(2,j)=411*1)；for(k(j):t(5,j)*p(3,j)=411*1)；for(k(
27、j):t(5,j)*p(4,j)=411*1)；for(k(j):t(5,j)*p(5,j)=411*1)；for(k(j):t(5,j)*p(6,j)=411*0)；for(k(j):s(1,j)-v(1,j)=0)；for(k(j):s(2,j)-v(2,j)=0)；for(k(j):s(3,j)-v(3,j)=0)；for(k(j):s(4,j)-v(4,j)=0)；for(k(j):s(5,j)-v(5,j)=0)；for(k(j):s(6,j)-v(6,j)=0)；for(k(j):s(1,j)-p(1,j)=0)；for(k(j):s(2,j)-l(1,j)-p(2,j)=0)；f
28、or(k(j):s(3,j)-l(2,j)-p(3,j)=0)；for(k(j):s(4,j)-l(3,j)-p(4,j)=0)；for(k(j):s(5,j)-l(4,j)-p(5,j)=0)；for(k(j):s(6,j)-l(5,j)-p(6,j)=0)；for(k(j):l(1,j)=100)；for(k(j):l(2,j)=100)；for(k(j):l(3,j)=100)；for(k(j):l(4,j)=100)；for(k(j):l(5,j)=100)；for(k(j):l(6,j)=50)；for(k(j):l(1,j)-p(1,j)-s(1,j)=0)；for(k(j):l(
29、2,j)-l(1,j)+p(2,j)-s(2,j)=0)；for(k(j):l(3,j)-l(2,j)+p(3,j)-s(3,j)=0)；for(k(j):l(4,j)-l(3,j)+p(4,j)-s(4,j)=0)；for(k(j):l(5,j)-l(4,j)+p(5,j)-s(5,j)=0)；for(k(j):l(6,j)-l(5,j)+p(6,j)-s(6,j)=0)；for(vol(i,j):gin(p(i,j)；名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 10 页，共
30、14 页 - - - - - - - - - 第 11 页 共 14 页YangQiuyan 第 11 页2013/4/17 for(vol(i,j):gin(s(i,j)；for(vol(i,j):gin(l(i,j)；for(vol(i,j):gin(v(i,j)；end 附录 2. 问题二中当 x值变化时的 Lingo程序代码model: sets: m/1.6/:； ! 定义月份下标集；k/1.7/:f,g； ! 定义产品种类；b/1.5/:； ! 定义机床种类；vol(m,k):v,p,l,s； ! 定义容量矩阵 v, 生产方案矩阵 p, 储存矩阵 l 和销售矩阵 s；time(b,
31、k):t； ! 定义时间矩阵 t ；endsets data: v=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 500 500 100 300 1100 500 60；f=10 6 8 4 11 9 3；t=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
32、 0.08 0 0 0.1 0 0.05 0 0.05；enddata max=sum(m(i):sum(k(j):g(j)*s(i,j)-0.5*sum(k(j):l(i,j)； ! 目标函数；for(k(j):g(j)-f(j)*(1-x)=0)；x=0.01； !值每次递加 0.02 ，直到 0.19 ；for(k(j):t(1,j)*p(1,j)=411*3)；for(k(j):t(1,j)*p(2,j)=411*4)；for(k(j):t(1,j)*p(3,j)=411*4)；for(k(j):t(1,j)*p(4,j)=411*4)；for(k(j):t(1,j)*p(5,j)=4
33、11*3)；for(k(j):t(1,j)*p(6,j)=411*4)；for(k(j):t(2,j)*p(1,j)=411*2)；for(k(j):t(2,j)*p(2,j)=411*2)；for(k(j):t(2,j)*p(3,j)=411*2)；for(k(j):t(2,j)*p(4,j)=411*1)；for(k(j):t(2,j)*p(5,j)=411*1)；for(k(j):t(2,j)*p(6,j)=411*2)；for(k(j):t(3,j)*p(1,j)=411*3)；for(k(j):t(3,j)*p(2,j)=411*1)；名师资料总结 - - -精品资料欢迎下载 - -
34、 - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 11 页，共 14 页 - - - - - - - - - 第 12 页 共 14 页YangQiuyan 第 12 页2013/4/17 for(k(j):t(3,j)*p(3,j)=411*3)；for(k(j):t(3,j)*p(4,j)=411*3)；for(k(j):t(3,j)*p(5,j)=411*3)；for(k(j):t(3,j)*p(6,j)=411*2)；for(k(j):t(4,j)*p(1,j)=411*1)；for(k(j):t(4,j)*p(2,j)=4
35、11*1)；for(k(j):t(4,j)*p(3,j)=411*0)；for(k(j):t(4,j)*p(4,j)=411*1)；for(k(j):t(4,j)*p(5,j)=411*1)；for(k(j):t(4,j)*p(6,j)=411*1)；for(k(j):t(5,j)*p(1,j)=411*1)；for(k(j):t(5,j)*p(2,j)=411*1)；for(k(j):t(5,j)*p(3,j)=411*1)；for(k(j):t(5,j)*p(4,j)=411*1)；for(k(j):t(5,j)*p(5,j)=411*1)；for(k(j):t(5,j)*p(6,j)=4
36、11*0)；for(k(j):s(1,j)-v(1,j)=0)；for(k(j):s(2,j)-v(2,j)=0)；for(k(j):s(3,j)-v(3,j)=0)；for(k(j):s(4,j)-v(4,j)=0)；for(k(j):s(5,j)-v(5,j)=0)；for(k(j):s(6,j)-v(6,j)=0)；for(k(j):s(1,j)-p(1,j)=0)；for(k(j):s(2,j)-l(1,j)-p(2,j)=0)；for(k(j):s(3,j)-l(2,j)-p(3,j)=0)；for(k(j):s(4,j)-l(3,j)-p(4,j)=0)；for(k(j):s(5,
37、j)-l(4,j)-p(5,j)=0)；for(k(j):s(6,j)-l(5,j)-p(6,j)=0)；for(k(j):l(1,j)=100)；for(k(j):l(2,j)=100)；for(k(j):l(3,j)=100)；for(k(j):l(4,j)=100)；for(k(j):l(5,j)=100)；for(k(j):l(6,j)=50)；for(k(j):l(1,j)-p(1,j)-s(1,j)=0)；for(k(j):l(2,j)-l(1,j)+p(2,j)-s(2,j)=0)；for(k(j):l(3,j)-l(2,j)+p(3,j)-s(3,j)=0)；for(k(j):
38、l(4,j)-l(3,j)+p(4,j)-s(4,j)=0)；for(k(j):l(5,j)-l(4,j)+p(5,j)-s(5,j)=0)；for(k(j):l(6,j)-l(5,j)+p(6,j)-s(6,j)=0)；for(vol(i,j):gin(p(i,j)；for(vol(i,j):gin(s(i,j)；for(vol(i,j):gin(l(i,j)；for(vol(i,j):gin(v(i,j)；名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 12 页，共 14 页 -
39、 - - - - - - - - 第 13 页 共 14 页YangQiuyan 第 13 页2013/4/17 end 附录 3.问题二中引进磨床时的程序代码model: sets: m/1.6/:； ! 定义月份下标集；k/1.7/:g； ! 定义产品种类；b/1.5/:； ! 定义机床种类；vol(m,k):v,p,l,s； ! 定义容量矩阵 v, 生产方案矩阵 p, 储存矩阵 l 和销售矩阵 s；time(b,k):t； ! 定义时间矩阵 t ；endsets data: v=500 1000 300 300 800 200 100 600 500 200 0 400 300 150
40、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；g=10 6 8 4 11 9 3；t=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.08 0 0 0.1 0 0.05 0 0.05；enddata max=sum(m(i):sum(k(j):g(j)*s(i,j)-0.5*sum(k(j):l(i,j)
41、； ! 目标函数；for(k(j):t(1,j)*p(1,j)=411*3)； !每次运行时递加1，直到 7；for(k(j):t(1,j)*p(2,j)=411*4)； !每次运行时递加1，直到 8；for(k(j):t(1,j)*p(3,j)=411*4)； !每次运行时递加1，直到 8；for(k(j):t(1,j)*p(4,j)=411*4)； !每次运行时递加1，直到 8；for(k(j):t(1,j)*p(5,j)=411*3)； !每次运行时递加1，直到 7；for(k(j):t(1,j)*p(6,j)=411*4)； !每次运行时递加1，直到 8；for(k(j):t(2,j)
42、*p(1,j)=411*2)；for(k(j):t(2,j)*p(2,j)=411*2)；for(k(j):t(2,j)*p(3,j)=411*2)；for(k(j):t(2,j)*p(4,j)=411*1)；for(k(j):t(2,j)*p(5,j)=411*1)；for(k(j):t(2,j)*p(6,j)=411*2)；for(k(j):t(3,j)*p(1,j)=411*3)；for(k(j):t(3,j)*p(2,j)=411*1)；for(k(j):t(3,j)*p(3,j)=411*3)；for(k(j):t(3,j)*p(4,j)=411*3)；for(k(j):t(3,j)
43、*p(5,j)=411*3)；for(k(j):t(3,j)*p(6,j)=411*2)；名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 13 页，共 14 页 - - - - - - - - - 第 14 页 共 14 页YangQiuyan 第 14 页2013/4/17 for(k(j):t(4,j)*p(1,j)=411*1)；for(k(j):t(4,j)*p(2,j)=411*1)；for(k(j):t(4,j)*p(3,j)=411*0)；for(k(j):t(4,j)
44、*p(4,j)=411*1)；for(k(j):t(4,j)*p(5,j)=411*1)；for(k(j):t(4,j)*p(6,j)=411*1)；for(k(j):t(5,j)*p(1,j)=411*1)；for(k(j):t(5,j)*p(2,j)=411*1)；for(k(j):t(5,j)*p(3,j)=411*1)；for(k(j):t(5,j)*p(4,j)=411*1)；for(k(j):t(5,j)*p(5,j)=411*1)；for(k(j):t(5,j)*p(6,j)=411*0)；for(k(j):s(1,j)-v(1,j)=0)；for(k(j):s(2,j)-v(2
45、,j)=0)；for(k(j):s(3,j)-v(3,j)=0)；for(k(j):s(4,j)-v(4,j)=0)；for(k(j):s(5,j)-v(5,j)=0)；for(k(j):s(6,j)-v(6,j)=0)；for(k(j):s(1,j)-p(1,j)=0)；for(k(j):s(2,j)-l(1,j)-p(2,j)=0)；for(k(j):s(3,j)-l(2,j)-p(3,j)=0)；for(k(j):s(4,j)-l(3,j)-p(4,j)=0)；for(k(j):s(5,j)-l(4,j)-p(5,j)=0)；for(k(j):s(6,j)-l(5,j)-p(6,j)=0
46、)；for(k(j):l(1,j)=100)；for(k(j):l(2,j)=100)；for(k(j):l(3,j)=100)；for(k(j):l(4,j)=100)；for(k(j):l(5,j)=100)；for(k(j):l(6,j)=50)；for(k(j):l(1,j)-p(1,j)-s(1,j)=0)；for(k(j):l(2,j)-l(1,j)+p(2,j)-s(2,j)=0)；for(k(j):l(3,j)-l(2,j)+p(3,j)-s(3,j)=0)；for(k(j):l(4,j)-l(3,j)+p(4,j)-s(4,j)=0)；for(k(j):l(5,j)-l(4,j)+p(5,j)-s(5,j)=0)；for(k(j):l(6,j)-l(5,j)+p(6,j)-s(6,j)=0)；for(vol(i,j):gin(p(i,j)；for(vol(i,j):gin(s(i,j)；for(vol(i,j):gin(l(i,j)；for(vol(i,j):gin(v(i,j)；end 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 14 页，共 14 页 - - - - - - - - -