天然肠衣搭配优化模型

1天然肠衣搭配优化模型摘要本文研究了天然肠衣在组装出成品时的搭配问题，在建立数学模型时，主要考虑了采用何种捆绑模式进行搭配才为最优方案，从而使得公司在天然肠衣的组装过程中捆绑效率提高，且在生产中可靠，实用。在模型假设允许的条件下，方案的设计中考虑到了两个标准：一是装出的成品捆数越多越好，二是对于成品捆数相同的方案，最短长度最长的成品越多，方案越好。为了使模型设计得到的方案在实际生产中简单，有效的进行。在解决肠衣搭配模型的问题中，将不同的三种规格CBA,,依次从后往前推理，先取最长规格的（我们假设C规格）肠衣，再将C规格剩余的肠衣降级使用。而后，我们再从较长的规格（我们假设为B规格）中取，再剩余的肠衣降级使用。最后从最短规格（我们假设为A规格）中取，剩余的作为废料。由于C规格中变量相对绝对值较大，所以我将C规格分成两种情况处理，20个变量中，前十个为C1规格，后十个为2C规格，其他限制条件不变。本文在问题的解决中采用数学的非线性规划模型进行求解，并运用lingo软件，使其与之相结合使用。在非线性规划中，为了满足条件中成品捆数越多越好，我们确定目标函数为iXmax，从而在问题的解决中能够得到最优的答案，通过lingo软件对模型中的约束性条件进行处理，使得非线性规划和lingo软件在模型的处理中起到了关键性作用，并使方案得以优化。关键词：肠衣搭配数学模型非线性规划lingo软件最优方案21问题重述肠衣在组装的时候总长度是一定的，但组装前的肠衣长度是不等的小段。为了提高生产效率，公司计划改变组装工艺，丈量了所有库存原料。且为了使肠衣长度均匀，将原料合理的分为三个规格进行搭配。现在要对这批原料进行搭配，为了使原料的使用率提高，总长度允许有±0.5的误差，总根数允许比标准少一根；当某种规格对应原料出现剩余时，可以降级使用，如长度为14米的原料可以和长度介于7-13.5米的进行捆扎，成品属于7-13.5米的规格；当成品捆数相同时，最短长度最长的成品越多越好；最后组装的成品捆数越多越好；为了食品保鲜，要求在30分钟内产生方案。根据以上条件，设计一个最优的原料搭配方案，使工人能“照方抓药”。2问题分析对于肠衣搭配问题首先要确定采用哪些捆绑模式，所谓捆绑模式，是指工人在现有的原料中所要取肠衣根数的方案。于是问题化为在简化工人工作的条件下，按照哪几种合理的模式，每种模式捆绑多少根怎样的原料肠衣最为优化。而所谓优化，可以有两种标准，一是装出的成品捆数越多越好，二是对于成品捆数相同的方案，最短长度最长的成品越多，方案越好。再者，为了问题简单化，题中允许总长度允许有±0.5米的误差，总根数允许比标准少1根；如果某种规格对应原料如果出现剩余，可以降级使用。如长度为14米的原料可以和长度介于7-13.5米的进行捆扎，成品属于7-13.5米的规格。那么为了满足以上要求，首先我们要从后往前推理，先取最长规格的（我们假设为C规格），再将C规格剩余的，降级使用。再者，我们再从较长的规格（我们假设为B规格）中取，剩余的肠衣降级使用。最后从最短规格（我们假设为A规格）中取，剩余的作为废料。由于C规格中变量相对绝对值较大，所以我将C规格变成两种，20个变量中，前十个为C1规格，后十个为C2规格。所以我们选择建立非线性规划模型分析求解，同时确定捆绑具体方案。33模型假设在加工肠衣时工人正常工作，原料被正确的按规格分配。生产的肠衣成品有CBA,,三种规格。不同规格的肠衣分配到不同的车间加工。零售商从肠衣厂进货时所获得的肠衣成品均为合格品。不考虑因操作不当而引起的肠衣损坏，即所生产的肠衣成品都是令人满意的。搭配过程中，肠衣不会变质。取j次正好能取尽。原料按长度分档，以0.5米为一档，如3-3.4米按3米计算，依此类推。`4符号假设ijr:表示按C种规格在第i（i=1,2,…20）种长度段中，第j（j=1,2,3）次取时，取得的根数；ijy:表示按B种规格在第i种长度段中，第j次取时，，取得的根数。ijt:表示按A种规格在第i种长度段中，第j次取时，，取得的根数。jx:表示第j次取得的梱数。iAsiBsiCs:表示表2中，相对应的根数。如：Cs1(即为14-14.4米的总根数)=35,Cs20（即为25.5-25.9米的总根数）=1，Bs1（即为7-7.4米的总根数）=24，Bs14(即为13.5-13.9米的总根数)=25，As1(即为3-3.4米的总根数)=43，As8(即为6.5-6.9米的总根数)=21。5模型建立Ⅰ、决策变量我们假设从14-14.4米段中第一次取11r根，从14.5-14.9米段中第一次取12r根，从15-15.4米段中第一次取31r根，依此类推，从25.5-25.9米段中第一次取201r根。设从7-7.4米段中第一次取11y根，从7.5-7.9米段中第一次取21y根，从8-8.4米段4中第一次取31y根，依此类推，从13.5-13.9米段中第一次取141y根；从7-7.4米段中第二次取12y根，从7.5-7.9米段中第二次取22y根，从8-8.4米段中第二次取32y根，依此类推，从13.5-13.9米段中第二次取142y根。设从3-3.4米段中第一次取11t根，从3.5-3.9米段中第一次取21t根，从4-4.4米段中第一次取31t根，依此类推，从6.5-6.9米段中第一次取81t根；从3-3.4米段中第二次取12t根，从3.5-3.9米段中第二次取22t根，从4-4.4米段中第二次取32t根，依此类推，从6.5-6.9米段中第二次取82t根;从3-3.4米段中第三次取13t根，从3.5-3.9米段中第三次取23t根，从4-4.4米段中第三次取33t根，依此类推，从6.5-6.9米段中第三次取83t根。显然，ijr,ijy,ijt均为正整数。Ⅱ、决策目标为了满足条件中成品捆数越多越好，我们确定目标函数;XiMaxⅢ、约束条件根据表1和表2(如下图)，我们对A,B,C1,C2四种规格进行如下约束：（1）在按C1和C2种规格捆绑时，为了使捆数最多，应有：5.89*5.8854)*(11ijiCijiCirsrsxr5（2）在按B种规格捆绑时，为了使捆数最多，应有：5.89*5.8887)*(ijiBijiBjijysysxy（3）在按A种规格捆绑时，为了使捆数最多，应有：5.89*5.882019)*(ijiAijiAjijtstsxt表1成品规格表最短长度最大长度根数总长度36.52089713.588914∞589表2原料描述表长度3-3.43.5-3.94-4.44.5-4.95-5.45.5-5.96-6.46.5-6.9根数4359394127283421长度7-7.47.5-7.98-8.48.5-8.99-9.49.5-9.910-10.410.5-10.9根数2424202521232118长度11-11.411.5-11.912-12.412.5-12.913-13.413.5-13.914-14.414.5-14.9根数3123225918253529长度15-15.415.5-15.916-16.416.5-16.917-17.417.5-17.918-18.418.5-18.9根数3042284245495064长19-19.19.5-1920-20.20.5-2021-21.21.5-2122-22.22.5-226度4.94.94.94.9根数526349352716122长度23-23.423.5-23.924-24.424.5-24.925-25.425.5-25.9根数060001根据约束条件和目标函数建立整数非线性规划模型。6模型求解将上述建立的整数非线性规划模型运用lingo软件求解，其求解程序见附录1，附录3,附录5,运行结果显示见附录2，附录4，附录6由求解结果可知，C1种规格最多能捆69捆，C2种规格最多能捆8捆，B种规格最多能捆32捆，A种规格最多能捆14捆。对于lingo求得方案，用下表表示：在C1规格中，进行两次搭配，搭配方案如表:表C1C1规格1919．52020．52121．52222．523．525．5W10000000221W20000001040在C2规格中，进行两次搭配，搭配方案如表:表C2C2规格1919．52020．52121．52222．523．525．5R10000003100R200000030107在B规格中，进行两次搭配，搭配方案如表:表BB规格77.588.599.51010.51111.51212.51313.5Y100010111011011Y201001000200400在A规格中，进行三次搭配，搭配方案如表:表AA规格33.544.555.566.5T170021180T234333022T3143305217模型讨论在本次搭配方案中，主要运用lingo软件，根据公司对搭配方案的要求，并且使目标函数XiMax取得最大值，我们对函数进行了四方面约束条件，一是，对取得的根数分三种规格进行要求；二是，对每捆肠衣的总长度进行了限制；三是，对每一小段中所取肠衣根数进行了要求；最后，在运行lingo语句时，要求取整。但是在操作中会有大量剩余，造成严重损失，并且开始假设3-3.4米按3米计算等，也会产生误差和损失。8参考文献[1]赫孝良等.数学建模竞赛赛题简析与论文点评.西安：西安交通大学出版社，2002.[2]陈理荣主编.数学建模导论.北京：北京邮电大学出版社，1999.[3]吴翊，吴孟达，成礼智编著.数学建模的理论与实践.长沙：国防科技大学出版社，1999[4]韩中庚.数学建模竞赛:获奖论文精选与点评.北京：科学出版社:第1版(2008年4月7日).9附录1：C1种规格中所取根数的lingo总程序代码：model:max=(x1+x2);r11*x1+r12*x2=35;r21*x1+r22*x2=29;r31*x1+r32*x2=30;r41*x1+r42*x2=42;r51*x1+r52*x2=28;r61*x1+r62*x2=42;r71*x1+r72*x2=45;r81*x1+r82*x2=49;r91*x1+r92*x23=50;r101*x1+r102*x2=64;r11+r21+r31+r41+r51+r61+r71+r81+r91+r101=4;r11+r21+r31+r41+r51+r61+r71+r81+r91+r101=5;r12+r22+r32+r42+r52+r62+r72+r82+r92+r102=4;r12+r22+r32+r42+r52+r62+r72+r82+r92+r102=5;14*r11+14.5*r21+15*r31+15.5*r41+16*r51+16.5*r61+17*r71+17.5*r81+18*r91+18.5*r101=88.5;14*r11+14.5*r21+15*r31+15.5*r41+16*r51+16.5*r61+17*r71+17.5*r81+18*r91+18.5*r101=89.5;14*r12+14.5*r22+15*r32+15.5*r42+16*r52+16.5*r62+17*r72+17.5*r82+18*r92+18.5*r102=88.5;14*r12+14.5*r22+15*r32+15.5*r42+16*r52+16.5*r62+17*r72+17.5*r82+18*r92+18.5*r102=89.5;@gin(x1);@gin(x2);@gin(x3);@gin(r11);@gin(r21);@gin(r31);@gin(r41);@gin(r51);@gin(r61);@gin(r71);@gin(r81);@gin(r91);@gin(r101);@gin(r12);@gin(r22);@gin(r32);@gin(r42)