天然肠衣搭配问题最新版本

2010高教社杯全国大学生数学建模竞赛承诺书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。

我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。

我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。

如有违反竞赛规则的行为，我们将受到严肃处理。

我们参赛选择的题号是（从A/B/C/D中选择一项填写）：我们的参赛报名号为（如果赛区设置报名号的话）：所属学校（请填写完整的全名）：参赛队员(打印并签名) ：1.2.3.指导教师或指导教师组负责人(打印并签名)：日期：年月日天然肠衣搭配问题摘要本文主要是对天然肠衣搭配问题进行探讨与研究，建立线性规划数学模型，利用LINGO编程，得到符合实际的最优方案。

文章以成品规格表和原料描述表为参考依据，采用整数规划，分别从“最大捆数、最优方案和时间限制”三个方面建立优化模型，利用LINGO编程最优求解，最终制作出了一套科学、合理和实用的搭配方案。

本文分三步解决问题，具体如下：首先，通过分析题设，按照要求（1）、（3）和（4），建立“最大捆数”的优化模型。

根据文中给出的三种成品规格，我们建立了三个求最大捆数的整数规划模型。

考虑到剩余原料可以降级规格使用，我们采用“倒序”原则，利用LINGO 编程，先算出第三种规格的最大捆数，接着把剩余原料23.5-23.9米的6根和25.5-25.9米的1根降级到第二种规格搭配使用，以此类推，LINGO运行得到三种规格的最大捆数，分别为134捆,41捆和16捆。

其次，根据最大捆数，本文得到两个具体的搭配方案。

方案一是根据材料使用情况建立最大损失函数模型，通过LINGO编程得到搭配方案（表45-）。

方案二是按照原料的最大利用原则建立优化模型，利用LINGO编程得到搭配方案（表5-），档次低的原料越55-），按照要求（2），比较两个方案的剩余原料（表6多，搭配方案越好。

我们最终选择方案二。

对剩余原料再次替代，得到优化方案（表75-）。

最后，考虑到食品保鲜，要求30分钟内产生方案，而第三规格原料数据太多，给程序运行带来困难，我们对模型改进。

把第三规格原料分成5个批次（表5-、1095-），对每个批次分别建模，求出的最大捆数分别32,32,32,18,19。

然后通过LINGO编程分别求出各个批次的搭配方案（表115-）。

5-、135-、12最后剩余的5根，利用替代原则，再优化1捆，得到最终的搭配方案（表145-）。

关键词：整数规划LINGO搭配方案最优求解天然肠衣1 问题重述天然肠衣的制作加工首先要经过清洗整理后被分割成长度不等的小段（原料），然后进入组装工序。

传统的生产方式就是人工边丈量原料长度边心算，将原材料按指定根数和总长度组装出成品（捆）。

原料按长度分档，以0.5米为一档，附表1是几种常见成品的规格，长度单位为米，∞表示没有上限，但实际长度小于26米。

为了提高生产效率，公司计划改变组装工艺，先丈量所有原料，建立一个原料表。

附表2为某批次原料描述。

根据以上成品和原料描述，设计一个原料搭配方案，工人根据这个方案“照方抓药”进行生产。

公司对搭配方案有以下具体要求：（1）对于给定的一批原料，装出的成品捆数越多越好；（2）对于成品捆数相同的方案，最短长度最长的成品越多，方案越好；（3）为提高原料使用率，总长度允许有0.5±米的误差，总根数允许比标准少1根；（4）某种规格对应原料如果出现剩余，可以降级规格使用。

如长度为14米的原料可以和长度介于13.5-7米的规格；-7米的进行捆扎，成品属于13.5（5）为了食品保鲜,要求在30分钟内产生方案。

根据上述问题建立数学模型，给出求解方法，并对附表1、附表2给出的实际数据进行求解，给出搭配方案。

2 问题分析问题中要求根据成品规格和原料描述，以及公司对搭配方案的具体要求，设计一个原料的最优搭配方案。

2.1对最大捆数的分析制作原料的搭配方案，首先要确定最大捆数。

考虑到文中要求（1）装出的成品捆数越多越好，因此把三个规格的最大捆数作为目标函数。

文中的要求（3）指出总长度可以允许5.0±米的误差，鉴于最大捆数的要求，把要求（3）作为约束条件对目标函数进行优化。

文中（4）要求，如果原料有剩余，可以考虑降级规格使用。

所以，先求第三种规格的最大捆数。

此规格中剩余的原料，降级使用，安排到第二种规格里面，作为新的长度档，参与第二种规格的分配。

同理，第二种规格中剩余的原料降级到第一种规格使用。

对建立的目标函数，根据题目的约束条件，进行LINGO编程，运行得到最大捆数为191捆。

2.2对最优搭配方案的分析在求出最大捆数的前提下，通过建立不同的目标函数得到两种优化模型：一种是按照原料的最大损失函数建立的优化模型，另一种是根据原料的最大利用原则建立优化模型，通过LINGO分别编程，得出两种搭配方案。

鉴于文中要求（2）提出：最短长度最长的成品越多，方案越好，我们选择第二种方案作为最终搭配方案。

要求（3）中提出：总长度允许有5.0±米误差、总根数允许比标准20根少一根。

在满足要求（3）的前提下，又把剩余的原料中档次较高的去替代档次较低的，使得搭配方案更加合理优化，达到了最优的搭配效果。

2.3改进的最优搭配方案根据第二种方案，利用LINGO编程，可以很快的求出原料第一、二规格的具体搭配方案。

但是在求第三种规格搭配方案时，由于数据较多，程序运行会占用很多时间，从而不满足要求（5）中的食品保鲜，即要求30分钟内产生方案。

本文对搭配方案模型进行改进，采用“化整为零”的思想优化模型：把第三规格原料分成5个批次，对每个批次分别建模，再用LINGO进行编程，求出最大捆数分别32,32,32,18,19。

然后通过对每个批次建立最优搭配模型，LINGO编程分别求出各个批次的搭配方案。

最后剩余的5根，利用替代原则，再优化1捆，得到最终的搭配方案（表145-）。

3 模型假设（1）假设分割后的肠衣没有损坏；（2）假设原料的丈量误差忽略不计；（3）假设原料的分档无误；（4）假设肠衣在捆装的过程中不出现断裂；（5）假设电脑运行正常；（6）假设在测量长度时不出现重叠；（7）假设工人可以看懂搭配方案；（8）假设肠衣在制作过程中不收缩变形；（9）假设不考虑温度、湿度等外界因素对肠衣质量的影响；4 符号说明1x ：第一种规格中的总捆数3x ：第二种规格中的总捆数5x ：第三种规格中的总捆数i x 1：第一种规格中的第i 档次 ()8,,3,2,1Λ=ii x 3：第二种规格中的第91-档次 ()9,,3,2,1Λ=ii x 4：第二种规格中的第1410-档次 ()4,3,2,1,0=ii x 5：第三种规格中得第91-档次 ()9,,3,2,1Λ=ii x 6：第三种规格中的第1810-档次 ()8,,2,1,0Λ=i ij a ：第一种规格的第i 捆中第j 档的根数 ()11,,3,2,116,,3,2,1ΛΛ==j iij b ：第二种规格的第i 捆中第j 档的根数 ()16,,3,2,141,,3,2,1ΛΛ==j i ij c ：第三种规格的第i 捆中第j 档的根数 ()20,,3,2,1134,,3,2,1ΛΛ==j i5 模型建立与求解1.5最大捆数的模型建立与求解根据附表1,2和公司要求（1），首先建立求最大捆数的优化模型]1[。

考虑原材料分为5.63-米、5.137-米、5.2514-米三个规格，根据每个规格求最大捆数，求和得到原料的总的最大捆数。

要求（4）中提到，剩余原料可以降级使用。

因此先求第三规格的最大捆数，建立目标函数]1[：5max x =要求（3）为提高原料的使用率，总长度允许有5.0±米的误差，以此约束条件：5958575655545352515185.17175.16165.15155.14145.88x x x x x x x x x x ++++++++≤6867666564636261605.22225.21215.20205.19195.18x x x x x x x x x +++++++++570695.895.255.23x x x ≤+ （5.1.1） 第三种规格中每捆5根，得到另一约束条件：560515x x i i =∑= （5.1.2）每档的根数不能超过现有根数，约束为：3551≤x 2952≤x 3053≤x 4254≤x 2855≤x 4256≤x4557≤x 4958≤x 5059≤x 6460≤x 5261≤x 6362≤x4963≤x 3564≤x 2765≤x 1666≤x 1267≤x 268≤x669≤x 170≤x （5.1.3） 利用LINGO 软件进行编程（附录一）计算]2[，得出第三种规格中各档使用根数（表5-1）和总的最大捆数1345=x 捆。

表15- 由表15-知，第三种规格中69x 和70x 有剩余。

根据要求()4中降级使用，把剩余的根数669=x 和170=x 降级到第二种规格进行计算（根据LINGO 内置算法编号不能有断层，所以将69x 和70x 重新编号为45x 和46x ），作为两个新档。

求第二种规格最大捆数，建立目标函数：3max x =约束条件：3938373635343332313115.10105.995.885.775.88x x x x x x x x x x ++++++++≤3464544434241405.895.255.235.13135.12125.11x x x x x x x x ≤+++++++（5.1.4）346318x x i i =∑= （5.1.5）2431≤x 2432≤x 2033≤x 2534≤x 2135≤x 2336≤x 2137≤x 1838≤x 3139≤x 2340≤x 2241≤x 5942≤x 1843≤x 2544≤x 645≤x 146≤x （5.1.6） 利用LINGO 软件进行编程（附录二）计算得出第二种规格各档使用根数（表5-2）和总的最大捆数413=x 捆。

表25-由表25-得出第二种规格中31x 、32x 和33x 有剩余，根据要求()4可降级使用，把剩余的根数2431=x 、732=x 和233=x 降级到第一种规格中进行计算。

得第一种规格最大捆数的目标函数：1m ax x =约束条件：1212019181716151413121115.8985.775.665.555.445.335.88x x x x x x x x x x x x x ≤++++++++++≤ （5.1.7） ∑==2111120i i x x（5.1.8） ∑==2111119i i x x（5.1.9）4311≤x 5912≤x 3913≤x 4114≤x 2715≤x 2816≤x3417≤x 2118≤x 2419≤x 720≤x 221≤x （5.1.10）利用LINGO 软件进行计算（附录三）得出第一种规格中各档使用根数（表5-3）和总的最大捆数为161=x 捆。