重型机械制造企业论文.doc

杭州电子科技大学第十二届大学生数学建模竞赛(承 诺 书)AB选择题号： （在方格内打）我们承诺：我们仔细阅读了校大学生数学建模竞赛的竞赛规则。我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们参赛的题目是： 重型机械制造企业的发展规划问题 我们的参赛报名序号为： 参赛队员 (打印并签名) ：1. 姓名（打印） 潘益坚 学号 09083220 学院 通信工程学院 专业 信息对抗技术 签名 2. 姓名（打印） 蔡钧宇 学号 09061532 学院 自动化 专业 自动化专业 签名 3. 姓名（打印） 刘忠秋 学号 09061538 学院 自动化 专业 自动化专业 签名 日期： 2011 年 5 月22 日 校评阅编号（由校数模组评阅前进行编号）： 27杭州电子科技大学第十二届大学生数学建模竞赛编 号 专 用 页校评阅编号（由校数模组评阅前进行编号）：重型机械制造企业的发展规划问题摘要：本问题要求我们寻找一种最优的资源配置方法以解决企业当前生产经营和下一步发展规划时的资源配置问题，关于资源配置的方式影响因素有很多，主要因素有：各基地产能、各地需求量、运输费用、材料费用、库存能力等。题目要求我们对各种因素进行综合分析，寻找出最优的资源配置方案，另外，题目还要求根据当前钢材市场的价格走势和产品需求，确定钢材的储备量及储备策略。对问题一，我们首先考虑的是运输费用最低的问题，这样才能保证企业利润的最大化。但是考虑到企业根据十二五规划的深入和西部大开发进程以及海外业务的拓展等等，企业考虑长久发展以及适应国家的政策，将会可能在短时间内以放弃最大利润为代价来换取企业长期的利益。由于企业生产的普通机械的市场属于销大于产的情况，所以我们在分配产品是要根据具体的企业发展战略来确定各个销售地的货物分配数量权重。我们在建立模型的时候，首先用权重表示企业根据具体情况对j销售地设置的权重，建立了一个通用模型，只要企业能够给我们足够的信息，我们就能够根据企业的规划来求出最优分配方案。当我们以无权重的时候算出来的最优结果为1474600元。对问题二，根据第一问的计算结果以及具体情况分析，总体上仍然是处于销大于产的情形，没有能够满足市场的需求，所以我们建立数学模型增加两个点进行运算求解，使其产销平衡。运用0-1规划来确定选址。在Lingo运行结果中，我们选用的两个建厂的地点为点8和点13，其运行结果总费用最少为8834000元。第三问是货物储备问题，对于企业的货物储备的决策问题，需要考虑的因素很多。首先就是储存费用，储存费与储存量得多少及储存时间的长短有关，这些储存费又包括自己的机会成本或利息，保藏费，维护费，安全费，场地费等等；同时由于钢材有可能受到腐蚀，所以在实际中也要考虑材料的报废和变质引起的损失；其二就是订购费用，每次订购货物产生的手续费，通信费，采购人员的差旅费以及货物的运输费用等。其特点是与订货次数有关，次数越多订购费用越大。其三是缺货引起的损失：缺货会引起顾客的流失，不能按期交货引起的企业信任度下降，如若签订合同而不能按期交货而产生的罚款等等。而对于本问题，我们所能够考虑的因素受到了限制，由题意可知，我们在此优先考虑当前钢材市场的价格走势和产品的需求对我们确定钢材的储备及储备策略的影响。要达到最优的储备策略，就是要所有的费用达到最少化。在第三问中，我们建立一个价格随时间变化率的约束条件，此条件能够隐形的帮助我们把储存费用和每次进货的费用考虑在内。关键字：0-1规划，线性规划，成本最小，权重，建址问题，库存问题，产销平衡，Dijistra最短路径算法 一、 问题综述： 某重型机械制造企业有8个生产基地分布在全国各地，其中v1,v6是已经建成并正常运作的生产基地，v7和v8正在建设中，预计一年以后可以投入常规产品的正常生产。v1,v18所有顶点都是该企业的主要销售地。该企业产品销售量在过去三年中每年增长速度都超过50%。我们要研究的第一个问题是在给定单位运费、各地间交通状况和各地具体产销数据的情况下，决定各产地产量及产品调运方案。第二个问题是两个在建生产基地建设规模现在还可以根据需要进行调整，设计规模是否合理。如果再增加两个生产基地，如果选择2个地点，多大规模？第三个个问题是该企业产品成本中钢材占50%，根据当前钢材市场的价格走势和产品需求，确定钢材的储备量及储备策略。二、 模型假设【1】生产企业在分配产品给销售地时，我们假设该企业要大力拓展海外业务，所以将会先满足含有出口的销售地对产品的需求。【2】在建厂选址时建厂的面积大小是一定的，我们假定为都是一亩，且已经足够我们建厂所需要的土地。【3】在研究钢材价格时我们假定以H型钢的价格为钢材的价格。【4】在库存问题中我们以价格因素和需求因素为主，其他因素影响相对较小。三、 变量说明：从i地到j地所需的最短距离：从i地到j地运输的机械数量：销售地j所需要的普通机械数量：产地i所能生产的最大普通机械数量：销售地j所需要的高端机械数量：产地i运往有出口地销售地后剩余的普通机械数量：产地i所能生产的最大高端机械数量：销售地j根据企业具体情况所设立的权重：第i地的建设费用：从i地到j地的运输量：第i次购进钢材所需要的价格：第i次购进钢材的数量：第i次购进钢材后能够使用的天数：p厂j季度所要生产的机械台数：第j厂全年的总产量四、 问题的分析及模型的建立问题的简化：由于运输费用是与线路的长度成线性关系，所以我们用Dijistra最短路径算法算出六个生产地分别运往十八个销售地最短路径，再乘以运输费用的单位价格就可将最短距离转换为最小费用。由于我们要在十八个地点选址建厂，所以我们为方便接直接求出每个点到其他各个点的最短距离，转换后的数据见附录的表一。编程时由于这是一个很典型的求两点间最短路的问题，所以我们将Dijistra最短路径算法改变后便成了我们所要的程序，编程时我们运用C语言将数据直接导入导出，具体程序见附表二。经过这一转化，问题大大简化，下面将原问题用数学语言描述.问题一的数据分析和模型建立：将产品合理的分配给各销售地，生产企业首先考虑的是运输费用最低的问题，这样才能保证企业利润的最大化。但是在实际中，我们需要考虑的因素还有很多，比如企业根据十二五规划的深入和西部大开发进程以及海外业务的拓展等等，企业考虑长久发展以及适应国家的政策，将会可能在短时间内以放弃最大利润为代价来换取企业长期的利益。由于企业生产的普通机械的市场属于销大于产的情况，所以我们在分配产品是要根据具体的企业发展战略来确定各个销售地的货物分配数量权重。当各个销售地有权重时，我们把权重分配到各个销售地的运输费用上，具体算法是用权重的倒数去乘以运输费用，这样当权重大的时候运输费用将会减少，这就能够体现出各个销售地对我们的重要程度。带权重的普通机械建模如下： 目标函数： 此为通用模型，在任何场合只要定好权重就可以进行合理的分配。由于我们对该企业的发展规划不清楚，我们知道没有权重时企业所得的利益最大，运费最小，所以我们先以没有权重的算法来分配货物给各销售点，然后分析分配方案，来决定是否应该调整分配方案，使出口所需以及国内所需的货物都能够适当的平衡，使他们都能够满意，我们通过计算得到国内所需与含出口所需所占的比重都接近总所需的百分之五十。下面我们对普通机械建立没有权重时的数学模型如下：目标函数： 运行的程序见附表三。其运行结果为将普通机械全部分配出去的最少运输费用是1376200元。此时的分配方案见下表：从i往j运输运输量销售总数出口总数国内总数X( 1, 17)740484029201920X( 2, 2)400X( 2, 15)280X( 2, 16)120X( 3, 3)800X( 4, 3)380X( 4, 4)420X( 5, 5)520X( 5, 14)480X( 6, 5)80X( 6, 6)300X( 6, 10)320从表格中我们可以看出，供出口的占总销售量的百分之六十，供国内的占百分之四十，也基本符合他们原来各自的相对的需求量，所以我们认为此时能够满足他们的需求。对于高端机械，从数据表我们可以看出总体上是产大于销的，所以我们建立的模型如下： 目标函数： 运行的程序见附表六。其运行结果为将高端机械全部分配出去的最少运输费用是98400.00元。此时的分配方案见下表：第一季度第二季度第三季度第四季度从i地到j地运输量从i地到j地运输量从i地到j地运输量从i地到j地运输量X( 1, 1)2X( 1, 1)2X( 1, 1)3X( 1, 1)1X( 1, 17)12X( 1, 14)3X( 1, 14)10X( 1, 17)25X( 1, 18)6X( 1, 17)16X( 1, 17)20X( 1, 18)12X( 3, 3)10X( 1, 18)8X( 1, 18)10X( 3, 3)12X( 6, 6)6X( 3, 3)15X( 3, 3)20X( 6, 6)10X( 6, 8)4X( 3, 14)5X( 6, 6)12X( 6, 14)8X( 6, 14)8X( 6, 6)12X( 6, 8)6X( 6, 8)6X( 6, 14)2X( 6, 14)2由此我们可以得出此时最少的运输费用为1376200+98400.00=1474600元。此时的总运费是最少的费用，因为我们没有参加任何的限制因素，是最理想化的分配方案。问题二的数据分析和模型建立： 根据我们第一问的计算结果以及具体情况分析， v7和v8的设计规模不合理。虽然v7和v8的产能已经很大，但是v7到v8的每季度的总产量只比第四季度的需求量多，总体上仍然是出于销大于产的情形，没有能够满足市场的需求。而且v7和v8建立的厂址比较偏僻，就算是建立的规模足够的大，也会有很大的运输费用，所以不够合理。按照现在的产能和需求情况，我们需要再选地建立新的生产基地。如果选择建立两个地点，我们应用0-1规划，根据最小费用选择厂址，使增加两个厂址以后能够满足市场的需求，下面我们建立数学模型来综合分析，首先我们将厂址建在各地的工厂面积都设为单位面积，并把建厂用地费用作为建厂费用。我们建立的数学模型如下： 其Lingo的运行程序见附表七，我们将程序运行结果导入到EXCEL表格中，具体见如下表格为生产基地运往各个销售地的运输量，在Lingo运行结果中，我们选用的两个建厂的地点为点8和点13，其运行结果总费用最少为8834000元。从i地运到j地运输量U( 1, 1)248U( 1, 3)63U( 1, 18)689U( 2, 2)400U( 2, 15)280U( 2, 16)120U( 3, 3)800U( 4, 3)380U( 4, 4)420U( 5, 5)386U( 5, 13)290U( 5, 14)322U( 5, 16)2U( 6, 5)200U( 6, 6)8U( 6, 12)208U( 6, 14)235U( 6, 16)49U( 7, 1)12U( 7, 3)57U( 7, 7)220U( 7, 17)1200U( 7, 18)11U( 8, 5)13U( 8, 6)137U( 8, 8)300U( 8, 11)3U( 8, 12)86U( 8, 13)9U( 8, 14)60U( 8, 16)10U( 9, 9)300U( 9, 11)456U( 9, 12)6U( 9, 14)7U( 9, 16)87U( 10, 5)1U( 10, 6)155U( 10, 10)320U( 10, 11)21U( 10, 13)1U( 10, 14)176U( 10, 16)182U( 11, 11)114U( 12, 12)129U( 13, 13)1094问题三的数据分析和模型建立：第三问是货物储备问题，对于企业的货物储备的决策问题，需要考虑的因素很多。首先就是储存费用，储存费与储存量得多少及储存时间的长短有关，这些储存费又包括自己的机会成本或利息，保藏费，维护费，安全费，场地费等等；同时由于钢材有可能受到腐蚀，所以在实际中也要考虑材料的报废和变质引起的损失；其二就是订购费用，每次订购货物产生的手续费，通信费，采购人员的差旅费以及货物的运输费用等。其特点是与订货次数有关，次数越多订购费用越大。其三是缺货引起的损失：缺货会引起顾客的流失，不能按期交货引起的企业信任度下降，如若签订合同而不能按期交货而产生的罚款等等。而对于本问题，我们所能够考虑的因素收到了限制，由题意可知，我们在此优先考虑当前钢材市场的价格走势和产品的需求对我们确定钢材的储备及储备策略的影响。要达到最优的储备策略，就是要所有的费用达到最少化。我们建立数学模型的思想是先通过在中国联合钢铁网数据库网站上查找三种钢铁在上海2010年1月1日到2010年12月32日的价格走势图，见附表九，通过分析这个价格走势图可知这三种刚才虽然都有自己的价格走势，可是在全年中大体的走势情况都是相同的，因此我们在买进钢材的过程中是同时买进三种钢材，使我们在每次使用钢材时能够正好按比例来储存与生产，而且分析图可知H型钢的价格平均会处在另外两种价格的中间，所以我们的钢材价格按照H型钢的价格来确定。并且在价格图中每一个最小极点都可以是一个购买点，则具体的数学模型如下：目标函数：其中最后一个约束条件是用来考虑库存费用以及每次购买的购买成本，这里先规定在图上的每一个极点可以是一个购买点，因为当下一个购买点的价格比现在的购买点价格高不了多少，而他们之间的时间间隔却很长，这样我们就可以在后面的购买点购买，这样会减少储存的费用。对于到底要在哪个购买点购买，就要我们通过对具体数据的分析来确定了，我们这里取的是通过计算所有购买点之间的价格增长率，然后计算这些价格增长率的平均值，并把它设为标准。 可以得到我们的储存方案见下表：时间/天 价格/元 购买量/吨 价格随时间的变化4 3980 524 1.0970 4050 8646 8.14230 3880 26593 2.1273 4230 无330 4350 7467 五、 模型的优缺点及改进方向模型一改进：对于第一个模型，为了能够看出对于每一个销售地点当我们分配不同权重的时候的分配方案以及我们的最少费用不同，我们以下又讨论了当我们先满足有出口的销售地，其他的销售地没有其他权重，所以将剩下的产品选用路费最少的最优分配方案给其他的销售地。所以对于有出口的销售地来说是产量大于销量，其数学模型如下： 目标函数为： 其Lingo应用程序见附表四；其运行结果为运费是3376000元，具体分配方案见下表：第一季度第二季度第三季度第四季度从i地到j地运输量从i地到j地运输量从i地到j地运输量从i地到j地运输量X( 1, 4)1000X( 1, 4)1000X( 1, 4)1000X( 1, 4)350X( 2, 5)600X( 2, 5)710X( 2, 5)700X( 1, 5)650X( 3, 1)600X( 3, 1)800X( 3, 1)600X( 3, 1)350X( 3, 4)200X( 4, 1)550X( 3, 4)200X( 3, 4)450X( 4, 1)700X( 4, 3)180X( 4, 1)700X( 4, 1)450X( 4, 5)100X( 4, 5)70X( 4, 5)100X( 5, 2)200X( 5, 2)200X( 5, 3)1000X( 5, 2)200X( 5, 3)450X( 5, 3)800X( 6, 2)680X( 5, 3)800X( 6, 2)400X( 6, 3)20X( 6, 2)300当我们把机械先分配给有出口的产地后，各地生产地剩余的机械数量附表所示。而此时是销大于产，我们继续将剩余的机械分配给各产地，其模型如下： 其Lingo应用程序见附表五：其运行结果为运费是3376000元，具体分配方案见下表：第一季度第二季度第三季度第四季度从i地到j地运输量从i地到j地运输量从i地到j地运输量从i地到j地运输量X( 2, 2)200X( 2, 2)90X( 2, 2) 100X( 2, 2) 450X( 6, 4)300X( 3, 1) 150 X( 6, 4) 380X( 2, 12)250 X( 3, 3) 650X( 6, 8)20X( 2, 13) 100X( 4, 3)300X( 4, 13)50X( 5, 13)350X( 6, 4)200X( 6, 8)500而此时的高端机械的分配方案和我们的第一问题的一样。由此，我们可以清晰地看出当我们把分配完全倾向于出口时，企业的运费会有很大增加。如若要实际的安排分配方案，就要综合分析公司的发展规划，长久打算以及国家的政策。有这些来分配各地的权重。模型二改进：在我们建立第二个模型时，由于我们不清楚生产规模与建厂面积之间的关系，我们就直接假设成单位面积来计算的，我们可以找到该企业生产规模与建厂面积之间的关系，来更加细化的分析合理建厂问题。模型三改进：在处理第三问题的时候，我们也分析了很多影响储存策略的因素，所以在改进模型时，如果有条件，我们可以尽可能多的影响因素考虑进去，是模型能够更合理，具有普遍意义。参考文献：【1】 百度百科，Dijkstra算法，/view/7839.htm，2011/5/20【2】陈光亭，裘哲勇，数学建模，北京：高等教育出版社，2010年2月。【3】中国联合钢铁网数据库，钢材价格参考数据来源， /data/steel_image.jsp?flag=top&iterms=2&product=20*2000*8000，2011/5/20附录附表表一：Dijistra求出来的每个点到其他各点的最短路径18到18最短距离从点V1出发到各点最短权值分别为：0,1520,1100,1400,1400,1600,1630,5100,2800,2100,3380,2680,1850,1750,1150,1000,730,940从点V2出发到各点最短权值分别为：1520,0,2620,2330,1630,1830,3150,5330,3030,2330,3610,2910,1530,1980,780,1330,2250,2460从点V3出发到各点最短权值分别为：1100,2620,0,300,1600,1800,2730,5300,3000,2300,3580,2880,2050,1250,1850,1300,1830,2040从点V4出发到各点最短权值分别为：1400,2330,300,0,1300,1500,3030,5000,2700,2000,3280,2580,1750,950,1550,1000,2130,2340从点V5出发到各点最短权值分别为：1400,1630,1600,1300,0,200,3030,3700,1400,700,1980,1280,450,350,850,400,2130,2340从点V6出发到各点最短权值分别为：1600,1830,1800,1500,200,0,3230,3500,1200,500,1780,1080,650,550,1050,600,2330,2540从点V7出发到各点最短权值分别为：1630,3150,2730,3030,3030,3230,0,6730,4430,3730,5010,4310,3480,3380,2780,2630,2360,2570从点V8出发到各点最短权值分别为：5100,5330,5300,5000,3700,3500,6730,0,4700,4000,5280,4580,4150,4050,4550,4100,5830,6040从点V9出发到各点最短权值分别为：2800,3030,3000,2700,1400,1200,4430,4700,0,700,2980,2280,1850,1750,2250,1800,3530,3740从点V10出发到各点最短权值分别为：2100,2330,2300,2000,700,500,3730,4000,700,0,2280,1580,1150,1050,1550,1100,2830,3040从点V11出发到各点最短权值分别为：3380,3610,3580,3280,1980,1780,5010,5280,2980,2280,0,700,2430,2330,2830,2380,4110,4320从点V12出发到各点最短权值分别为：2680,2910,2880,2580,1280,1080,4310,4580,2280,1580,700,0,1730,1630,2130,1680,3410,3620从点V13出发到各点最短权值分别为：1850,1530,2050,1750,450,650,3480,4150,1850,1150,2430,1730,0,800,750,850,2580,2790从点V14出发到各点最短权值分别为：1750,1980,1250,950,350,550,3380,4050,1750,1050,2330,1630,800,0,1200,750,2480,2690从点V15出发到各点最短权值分别为：1150,780,1850,1550,850,1050,2780,4550,2250,1550,2830,2130,750,1200,0,550,1880,2090从点V16出发到各点最短权值分别为：1000,1330,1300,1000,400,600,2630,4100,1800,1100,2380,1680,850,750,550,0,1730,1940从点V17出发到各点最短权值分别为：730,2250,1830,2130,2130,2330,2360,5830,3530,2830,4110,3410,2580,2480,1880,1730,0,1670从点V18出发到各点最短权值分别为：940,2460,2040,2340,2340,2540,2570,6040,3740,3040,4320,3620,2790,2690,2090,1940,1670,0附表二：Dijistra最短路径算法程序：注：此程序要根据程序说明使用，我们运用数据导入导出来处理#include#includeusing namespace std;const int MaxNum=1000000; /边权最大值int n; /节点数目int dist501; /到节点1的最短路径值bool state501; /节点被搜索过状态指示int data501501; /邻接矩阵/查找权值最小的节点int findmin()int minnode=0, min=MaxNum;for(int i=1; i=n; i+)if (distimin) & (!statei)min=disti;minnode=i;return minnode;int main()ifstream in(dijkstrainput.txt);ofstream out(18点到18点距离.txt);for(int s=1;s n;for(int p=1; p=n; p+)for(int q=1; q datapq;if (datapq=0) datapq=MaxNum;/初始化for(int i=1; i=n; i+)disti=datasi;states=true;int done=1;while (donen)int node=findmin();if (node!=0)done+; /找到的点的数目加1statenode=true; /标记已经找到了从节点1到节点node的最短路径for(int i=1; idistnode+datanodei) & (!statei)disti=distnode+datanodei;else break;out从点Vs出发到各点最短权值分别为：endl;for(int k=1; k=n; k+)if (distk=MaxNum)out0;elseoutdistk;if (k=n)outendl;elseout,;in.close();out.close();return 0;附表三：不考虑权重的普通机械分配方案不考虑出口 model:sets:chandi/1.6/:a;xiaodi/1.18/:b,c,d,e;route(chandi,xiaodi):s,x;endsetsdata:a=1000,800,800,800,1000,700;b=260,400,1300,420,600,300,220,300,300,320,480,300,300,800,280,450,1200,700;c=300,480,1350,650,680,350,400,600,450,260,650,500,550,1200,400,500,1000,780;d=350,550,1300,400,500,380,500,460,400,600,400,320,400,800,350,780,1200,800;e=280,450,800,300,200,200,650,200,180,650,200,250,320,450,250,700,800,650;s=0,1520,1100,1400,1400,1600,1630,5100,2800,2100,3380,2680,1850,1750,1150,1000,730,940,1520,0,2620,2330,1630,1830,3150,5330,3030,2330,3610,2910,1530,1980,780,1330,2250,2460,1100,2620,0,300,1600,1800,2730,5300,3000,2300,3580,2880,2050,1250,1850,1300,1830,2040,1400,2330,300,0,1300,1500,3030,5000,2700,2000,3280,2580,1750,950,1550,1000,2130,2340,1400,1630,1600,1300,0,200,3030,3700,1400,700,1980,1280,450,350,850,400,2130,2340,1600,1830,1800,1500,200,0,3230,3500,1200,500,1780,1080,650,550,1050,600,2330,2540;enddatamin=sum(route:s*x);for(chandi(i):sum(xiaodi(j):X(i,j)=a(i);for(xiaodi(j):sum(chandi(i):X(i,j)B(j); !for(xiaodi(j):sum(chandi(i):X(i,j)C(j); !for(xiaodi(j):sum(chandi(i):X(i,j)D(j); !for(xiaodi(j):sum(chandi(i):X(i,j)E(j); end 附表四：出口的先考虑：产大于销model:sets:chandi/1.6/:a;chukou/1.5/:B,C,D,E;route(chandi,chukou):S,X;endsetsdata:a=1000,800,800,800,1000,700;B=1300,600,800,1200,700;C=1350,680,1200,1000,780;D=1300,500,800,1200,800;E=800,200,450,800,650;S=1100,1400,1750,730,940, 2620,1630,1980,2250,2460, 0,1600,1250,1830,2730, 300,1300,950,2130,2340, 1600,0,350,2130,2340,1800,200,550,2330,2540;enddatamin=sum(route:S*X);for(chandi(i):sum(chukou(j):X(i,j)a(i);for(chukou(j):sum(chandi(i):X(i,j)=B(j); !for(chukou(j):sum(chandi(i):X(i,j)=C(j); !for(chukou(j):sum(chandi(i):X(i,j)=D(j); !for(chukou(j):sum(chandi(i):X(i,j)=E(j); end附表五:剩余的13个销地 此时销大于产第一季度model:sets:chandi/1.6/:a;xiaodi/1.13/:b,c,d,e;route(chandi,xiaodi):s,x;endsetsdata:a=0,200,0,0,0,300;b=260,400,420,300,220,300,300,320,480,300,300,280,450;c=300,480,650,350,400,600,450,260,650,500,550,400,500;d=350,550,400,380,500,460,400,600,400,320,400,350,780;e=280,450,300,200,650,200,180,650,200,250,320,250,700;s=0,1520,1400,1600,1630,5100,2800,2100,3380,2680,1850,1150,1000 1520,0,2330,1830,3150,5330,3030,2330,3610,2910,1530,780,1330,1100,2620,300,1800,2730,5150,2850,2150,3430,2730,1900,1850,1300,1400,2330,0,1500,3030,4850,2550,1850,3130,2430,1600,1550,1000,1400,1630,1300,200,3030,3700,1400,700,1980,1280,450,850,400,1600,1830,1500,0,3230,3500,1200,500,1780,1080,650,1050,600;enddatamin=sum(route:s*x);for(chandi(i):sum(xiaodi(j):X(i,j)=a(i);for(xiaodi(j):sum(chandi(i):X(i,j)B(j); !for(xiaodi(j):sum(chandi(i):X(i,j)C(j); !for(xiaodi(j):sum(chandi(i):X(i,j)D(j); !for(xiaodi(j):sum(chandi(i):X(i,j)E(j); end 第二季度model:sets:chandi/1.6/:a;xiaodi/1.13/:b,c,d,e;route(chandi,xiaodi):s,x;endsetsdata:a=0,90,800,0,0,0;b=260,400,420,300,220,300,300,320,480,300,300,280,450;c=300,480,650,350,400,600,450,260,650,500,550,400,500;d=350,550,400,380,500,460,400,600,400,320,400,350,780;e=280,450,300,200,650,200,180,650,200,250,320,250,700;s=0,1520,1400,1600,1630,5100,2800,2100,3380,2680,1850,1150,1000 1520,0,2330,1830,3150,5330,3030,2330,3610,2910,1530,780,1330,1100,2620,300,1800,2730,5150,2850,2150,3430,2730,1900,1850,1300,1400,2330,0,1500,3030,4850,2550,1850,3130,2430,1600,1550,1000,1400,1630,1300,200,3030,3700,1400,700,1980,1280,450,850,400,1600,1830,1500,0,3230,3500,1200,500,1780,1080,650,1050,600;enddatamin=sum(route:s*x);for(chandi(i):sum(xiaodi(j):X(i,j)=a(i);!for(xiaodi(j):sum(chandi(i):X(i,j)B(j); for(xiaodi(j):sum(chandi(i):X(i,j)C(j); !for(xiaodi(j):sum(chandi(i):X(i,j)D(j); !for(xiaodi(j):sum(chandi(i):X(i,j)E(j); end 第三季度model:sets:chandi/1.6/:a;xiaodi/1.13/:b,c,d,e;route(chandi,xiaodi):s,x;endsetsdata:a=0,100,0,0,0,400;b=260,400,420,300,220,300,300,320,480,300,300,280,450;c=300,480,650,350,400,600,450,260,650,500,550,400,500;d=350,550,400,380,500,460,400,600,400,320,400,350,780;e=280,450,300,200,650,200,180,650,200,250,320,250,700;s=0,1520,1400,1600,1630,5100,2800,2100,3380,2680,1850,1150,1000 1520,0,2330,1830,3150,5330,3030,2330,3610,2910,1530,780,1330,1100,2620,300,1800,2730,5150,2850,2150,3430,2730,1900,1850,1300,1400,2330,0,1500,3030,4850,2550,1850,3130,2430,1600,1550,1000,1400,1630,1300,200,3030,3700,1400,700,1980,1280,450,850,400,1600,1830,1500,0,3230,3500,1200,500,1780,1080,650,1050,600;enddatamin=sum(route:s*x);for(chandi(i):sum(xiaodi(j):X(i,j)=a(i);!for(xiaodi(j):sum(chandi(i):X(i,j)B(j); !for(xiaodi(j):sum(chandi(i):X(i,j)C(j); for(xiaodi(j):sum(chandi(i):X(i,j)D(j); !for(xiaodi(j):sum(chandi(i):X(i,j)E(j); end 第四季度model:sets:chandi/1.6/:a;xiaodi/1.13/:b,c,d,e;route(chandi,xiaodi):s,x;endsetsdata:a=0,800,0,350,350,700;b=260,400,420,300,220,300,300,320,480,300,300,280,450;c=300,480,650,350,400,600,