天然肠衣搭配问题经典

1、池州学院天然肠衣搭配问题组员：陈强赵晋彪 赵海龙目录一、问题重述 错误!未定义书签。问题背景 错误!未定义书签。.问题条件 错误!未定义书签。.问题要求 错误!未定义书签。需要解决的问题 错误!未定义书签。二、问题分析 错误!未定义书签。三、模型假设 错误!未定义书签。四、符号说明 错误!未定义书签。五 模型的建立 错误!未定义书签。、模型建立 错误!未定义书签。、根据要求模型建立 错误!未定义书签。六、模型求解 错误!未定义书签。、问题要求（1）模型求解 错误!未定义书签。、问题要求（2）模型求解 错误!未定义书签。、问题要求（3）模型求解 错误!未定义书签。七、模型的评价与推广 错误!未定

2、义书签。.模型的评价 错误!未定义书签。模型的优点 错误!未定义书签。模型的缺点 错误!未定义书签。模型的推广 错误!未定义书签。八、参考文献 错误!未定义书签。附录 错误!未定义书签。附录A 错误!未定义书签。附录B 错误!未定义书签。附录C 错误!未定义书签。附录D 错误!未定义书签。天然肠衣搭配问题摘要天然肠衣制作加工是我国的一个传统产业，出口量占世界首位，而天然肠衣传统的人工生产方式已不能满足出口量日益增长的需要。因此， 我们从节约生产成本、 提高企业生产效率的角度出发，我们结合原料的供给量、长度及成品规格等约束条件进行了模型设计。根据题目中的表1 中的成品的规格和表2 中的原料，我们

3、所需要解决的问题有： 如何搭配才能使得成品的捆数最多对于针对这一个问题我们采用线性规划建立模型并利用MATLAB以捆数相同，最短长度越长越好的原则，求得模型的最优解。另外，由于所有的原料按长度分档，通常以0.5米为一档，如：3-3.4米按3米计算，3.5米 -3.9米按3.5米计算，其余的依此类推。表1是几种常见成品的规格，长度单位为米，表示没有上限，但实际长度小于26米。再把不同档次的原料按照不同的规格进行搭配，分别搭配成三种规格的成品，依次是成品一（3米 ,20根，总长度89米），成品二（7一米，8根，总长度89米），成品三（148米，5根，总长度89米） 。 运用线性规划分别对成品一、成

4、品二、 成品三建立模型，利用 LINGO编程进行1步，2步，3步优化筛选，得出方案。并且，对各步筛选所剩余原料再进行优化得出方案，另外， 为了提高原材料的使用率，每成品的总长度允许有米的误差，总根数允许比标准少1根，某种成品对应得原材料有剩余，可以降一级使用，这样就会出现每捆总长度米和米，有19根一捆，7根一捆，4根一捆，在满足条件时，计算出最大捆数。关键词：天然肠衣； 线性规划； MATLAB ;LINGO一、问题重述问题背景天然肠衣制作加工是我国的一个传统产业，出口量占世界首位，而天然肠衣传统的生产依靠工人，边丈量长度，边心算，将原材料按指定根数和总长度组装成捆方式已不能满足出 口量日益增

5、长的需要。因此，我们从节约生产成本、提高企业生产效率的角度出发，结合原 料的供给量、长度及成品规格等约束条件进行生产模型设计。.问题条件已知某天然肠衣组装的规格表（见表1）和某批次原料的描述表（见表2）. 表1成品规格表成品最短长度/m取kk度/m根数总长度/m成品一32089成品二7889成品三14OO589表2原料描述表长 度根 数4359394127283421长 度根 数2424202521232118长 度根 数3123225918253529长 度根 数3042284245495064长 度根 数526349352716122长 度根 数060001.问题要求根据生产规格和原料描述

6、表，设计满足以下要求及允许条件下的组装方案 （即对原料进行打捆的搭配方案）。1）对于给定的一批原料，装出的成品时，捆数越多越好；2）对于成品捆数相同的方案，最短长度最长的成品越多，方案越好；3）允许成品总长度有土m的误差，即总长度米，米；4）每捆总根数允许比标准少1根即每捆可以为19根，7根，4根；5）某种成品对应原料如果出现剩余，可以降级使；。6）为了食品安全要求在30分钟内完成。需要解决的问题首先我们要考虑原料应该怎样处理，然后考虑各类成品原料的搭配问题以及 剩余材料的得搭配问题。、问题分析天然肠衣经过人工清理和截取成不等长度，在进行搭配成成品，传统工艺 里，工人们边劳作用心丈量，其操作不

7、仅效率差，有时会造成严重资源浪费，为 了解决这一问题即做到同等材料中，加工成成品后，剩余材料最少，捆数最多， 我们依据表一，表二的生产规格和所拥有的原料进行生产搭配成成品，作出最有模型。对于问题要求一，我们运用 lingo软件编程求出优化方案，并作出逐步优化，得出最 大捆数；对于问题二，利用 matlab编程，按照最短长度最长的原则选出最优方案，将选出 的方案与lingo结合编程，得出最大捆数，然后在剩余中再按照最短长度最长的原则选择最后所能组成的捆数；对于问题要求三、四、五，对于成品一规格，要求长度11满足3m l 6.5m, l l或11 11 0.5根数x1满足Xi 20根；对于成品二规

8、格，要求长度l2满足7m l2 14.5, l2l2或l2l2 +根数X2满足X2 8根；对于成品三规格，要求长度l3满足14m< Lw26m, l3 L或L3T L3+,第三类原料最大长度为8,但原料的实际长度小于26米根数x3满足5根。在每个档次的原料满足这些约束条件的前提下，建立不定方程组，用 lingo软件求解出每种规格产品的最大捆数；为提高原料使用率，总长度允许有± 0.5米的误差，总根数允许比标准少1 根。 因此生产条件将变为每捆天然肠衣成品的总长度为89 ± 0.5m,同时可能会出现成品一每捆十九根，成品二每捆七根，成品三每捆lingo 软件求解出每种规

9、格产品的最大捆数。三、模型假设1、 假设题目所给的数据真实可靠；2、 所有原料都是新鲜的，没有腐烂；3、 所有原料都是同种原料，不需要对原料进行分类；4、 厂房的环境、车间和工作人员的卫生等不会对原料产生影响；5、 厂房中的温度，光照等不会对原料的保鲜产生影响。四、符号说明xii (i1,2,3.46)每个档次在组装中所使用的根数yi (i1,2,3.)不同规格成品按某种方案所组装出来的捆数zi (i1,2,3) 每种规格所能够组合出来的最大捆数Z 各规格的最大捆数的总和xis(i1,2,3.46)各个档次在经过缩小取值范围后的取值yij(i1,2,3, j 1,2,3.) 规格i 的成品第

10、j次按照某个方案组装所的捆数xiy各个档次经过组装的剩余量ni 档次 i 的原料开始时的数量Li (i 1,2,3) 每种规格所能取得的最大理想捆数Si (i1,2,3) 规格 i 的原料在经过要求(2)的筛选后剩余的总根数 五 模型的建立、模型建立成品一成品一要求是选出20 根不同肠衣组装成1 捆，原料长度按米为一档，如：3按 3 米计算，所选出得肠衣总长度之和为89 米，可以得出x1 x2 x3 x4 x5 x6 x7 x8208其中x1X8的限制条件为:x143X259X349X441X1 X8 的根数会减少一部分，减少后限X5 27X6 28X734X828按照某个方案组装yi(i 1

11、,2,3.) 捆之后， 制条件变化为X143yiX527yiX259yiX628yiX349yiX734yiX441yiX821yi经过 k 次组装之后，规格1 的成品的捆数要求最大，可以得出nMaXz1yii1成品二规格 2 的要求是选出8 根不同肠衣组装成1 捆， 所选出得肠衣总长度之和为89 米，可以得出22Xi8i922i9X9X22的根数限制条件为X 9 24X 10 24X 11 20X 12 25X 13 21X 14 23X 15 21按照某个方案组装yi(i1,2,3.) 捆之后，X16 18X 17 31X 18 23X 19 22X 20 59X 21 18X 22 25

12、X9 X22的根数会减少一部分，减少后限制条件变化为7 0.5 i 9 Xi89x924yix1618yix 1024yix1731yix1120yix1823yix 1225yix1922yix1321yix 2059yix1423yix 2118yix1521yix 2225yi经过 k 次组装之后，规格2 成品的捆数要求最大，可以得出nMaxz 2yii1成品三成品 3 要求是选出5 根不同肠衣组装成1 捆， ，所选出得肠衣总长度之和为89 米，可以得出46xi 5i 2314x23 14.5x24 15x25 15.5x2616 x2716.5x28 .26x46 89x23 x46

13、的根数限制条件为x2335x3549x2429x3635x2530x3727x2642x3816x2728x3912x2842x402x2945x410x3049x426x3150x430x3264x440x3352x450x3463x461按照某个方案组装yi(i1,2,3.) 捆之后，x23 x46 的根数会减少一部分，减少后限制条件变化为x2335yix3549yix2429yix3635yix2530yix3727yix2642yix3816yix2728yix3912yix2842yix402yix2945yix410yix3049yix426yix3150yix430x3264yi

14、x440x3352yix450x3463yix461yi经过 k 次组装之后，成品三成品的产品的捆数要求最大，可以得出nMaxz3yii1最后得出成品的最大捆数为Z z1 z2 z3、根据要求模型建立成品一在第一轮筛选后的余量总根数(即满足问题要求2)S1 1 5 2 2 1 1 12由于题目要求是20跟一捆，即便充分利用问题要求(3)和(4)，也达不到题中的要求，所以成品一不能再搭配出多余的捆数。成品二同样经过要求(2)的筛选后剩下的总根数S2 24 2 20 3 21 5 6 1 82根数可以达到8 根一捆，再考虑一捆的总长度为88.5、 88、 89.5米和每捆为8根， 8 根不行，再考

15、虑每捆7 根。成品三也是经过要求(2)筛选后剩下的总根数为S3 1 1 1 2 1 12 2 6 1 27按题中 5 根一捆，根数能够满足要求，再加上每捆总长度在89 0.5米，每捆按 5 或求其次4 根计算。目标函数：maxY约束条件：464614 0.5(i 23) xi 89xi 5yx230i 23i 23yx35 2yx240yx360yx250yx370yx260yx381yx271yx3912yx280yx402yx291yx410yx300yx426yx310yx430yx320yx440yx330yx，50yx341yx461六、模型求解、问题要求（1）模型求解对于不同规格的

16、成品，会存在一个理想的最大捆数S,但是由于实际条件的 约束，总捆数不可能超过S,其中理想的最大捆数等于此种规格的总根数之和除 于此种规格每捆的根数要求成品一的最大理想捆数43 59 39 34 21 1420成品二的最大理想捆数L2成品三的最大理想捆数L3余下12根24 24 20 18 25 448余下2根35 29 30 .0 1 1355余下2根成品一利用LINGO编程（程序见附录A）,对成品一进行最优化筛选，可得方案表3成品一的原料第1轮筛选的方案xix1x2x3x4x5x6x7x8选取根 数25332131方案所得捆数：11经过1轮的优化之后，并未达到最优，再进行优化，有方案 表4筛

17、成品一的原料第 2筛选的方案XiX1X2X3X4X5X6X7X8选取根 数71211503方案所得捆数：3成品一的实际的总捆数之和 Maxz1 11 3 14,实际的总捆数之和4与最 大理想总捆数G相等，成品一的优化结束。各个档次剩余根数分别为表5成品一的原料经过筛选所剩余的原料XiX1X2X3X4X5X6X7X8剩余根 数01050211成品二利用同样的方法对成品二进行优化分配处理表6成品二的原料第1轮筛选的方案XiX9X10X11X12X13X14X15选取根数0101000XiX16X17X18X19X20X21X22选取根数0111201方案所得捆数为：22表7成品二的原料第2轮筛选的

18、方案XiX9X10X11X12X13X14X15选取根数0000021XiX16X17X18X19X20X21X22选取根数1100120方案所得捆数为:表8成品二的原料第3轮筛选的方案XiX9X10X11X12X13X14X15选取根数0000012xix16x17x18x19x20x21*22选取根数2000201方案所得捆数为：3经3轮的筛选之后，各个档次的根数剩余量为表9成品二的原料3轮筛选的剩余量xix9x10X11X12X13X14X15剩余根数2422032126xiXi6X17X18X19X20X21X22剩余根数3010000剩余的原料不能再组装出符合要求的成品，优化结束。最

19、后得出成品二的原料所能够组装的最大实际捆数为 Maxz2 22 9 3 34成品二成品三的原料采用同样方法，经过九轮筛选之后剩余原料跟数为: 表10成品三的原料9轮筛选的剩余量XiX23X24X25X26X27X28X29X30X31X32X33X34剩余 根数000010100001XiX35X36X37X38X39X40X41X42X43x44X45X46剩余 根数2001122060001经过9轮的优化筛选之后，剩余的原料无法再组装出符合要求的成品， 最后 求得成品三的实际最大捆数为Maxz3 42+30+26+14+9+4+2+1+1=129所以，利用LINGO对题目中给定的原料求解出

20、最大能够组装出的捆数MaxZ Maxz1 Maxz2 Maxz3 14 34 129 177、问题要求（2）模型求解利用MATLAB编程（程序见附录B）,按照最短长度最长的原则选出如下方 案，将选出的方案结合LINGO编程（程序见附录C）,得出捆数为11表11成品一的原料按照最短长度最长优化的方案XiXiX2X3X4X5X6X7X8选取根 数25332131在剩余中再按照最短长度最长的方案选择 3捆，方案如下表12成品一的原料按照最短长度最长优化的方案XiX1X2X3X4X5X6X7X8选取根 数71211503成品二在符合成品二的材料中按照最短长度最长的方案选择22捆，方案如下表13成品二的

21、原料按照最短长度最长优化的方案XiX9X10X11X12X13X14X15选取根数0101000XiX16X17X18X19X20X21X22选取根数0111201在剩余中再按照最短长度最长的方案选择 9捆，方案如下 表14成品二的原料按照最短长度最长优化的方案XiX9X10X11X12X13X14X15选取根数0000021XiX16X17X18X19X20X21X22选取根数1100120再在剩余中再按照最短长度最长的方案选择 3捆，方案如下 表15成品二的原料按照最短长度最长优化的方案XiX9X10X11X12X13X14X15选取根数0000012XiX16X17X18X19X20X2

22、1X22选取根数2000201成品二在符合成品三的材料中按照最短长度最长的方案进行选择42捆，方案如下表16成品三的原料按照最短长度最长优化的方案XiX23X24X25X26X27X28X29X30X31X32X33X34选取 根数000001101110XiX35X36X37X38X39X40X41X42X43X44X45X46选取 根数000001000000在剩余中再按照最短长度最长的方案选择 31捆，方案如下 表17成品三的原料按照最短长度最长优化的方案XiX23X24X25X26X27X28X29X30X31X32X33X34选取 根数100100000002XiX35X36X37X

23、38X39X40X41X42X43X44X45X46选取 根数010000000000在剩余中再按照最短长度最长的方案选择 27捆，方案如下 表18成品三的原料按照最短长度最长优化的方案XiX23X24X25X26X27X28X29X30X31X32X33X34选取 根数010010010000XiX35X36X37X38X39X40X41X42X43X44X45X46选取 根数101000000000在剩余中再按照最短长度最长的方案选择 15捆，方案如下 表19成品三的原料按照最短长度最长优化的方案XiX23X24X25X26X27X28X29X30X31X32X33X34选取 根数0020

24、00010000XiX35X36X37X38X39X40X41X42X43X44X45X46选取 根数100100000000在剩余中再按照最短长度最长的方案选择 7捆，方案如下 表20成品三的原料按照最短长度最长优化的方案XiX23X24X25X26X27X28X29X30X31X32X33X34选取 根数000100001300XiX35X36X37X38X39X40X41X42X43X44X45X46选取 根数000000000000在剩余中再按照最短长度最长的方案选择 4捆，方案如下 表21成品三的原料按照最短长度最长优化的方案XiX23X24X25X26X27X28X29X30X31

25、X32X33X34选取 根数100100000010XiX35X36X37X38X39X40X41X42X43X44X45X46选取 根数110000000000在剩余中再按照最短长度最长的方案选择 2捆，方案如下 表22成品三的原料按照最短长度最长优化的方案XiX23X24X25X26X27X28X29X30X31X32X33X34选取 根数010000010030XiX35X36X37X38X39X40X41X42X43X44X45X46选取 根数000000000000在剩余中再按照最短长度最长的方案选择 1捆，方案如下 表23成品三的原料按照最短长度最长优化的方案XiX23X24X25

26、X26X27X28X29X30X31X32X33X34选取 根数000000031100XiX35X36X37X38X39X40X41X42X43X44X45X46选取 根数000000000000在剩余中再按照最短长度最长的方案选择 1捆，方案如下表24成品三的原料按照最短长度最长优化的方案XiX23X24X25X26X27X28X29X30X31X32X33X34选取 根数000000220000XiX35X36X37X38X39X40X41X42X43X44X45X46选取 根数100000000000按照上述成品三选择原材料进行组装，即可使得这批原材料的所组装出来的 成品最短长度最长，

27、此时各个规格的捆数之和为 178,相比要求1所得的捆数多 出了 1捆，这是因为成品三在选择最短长度最长的方案组装的时候使得材料能够 更加充分利用、问题要求（3）模型求解通过LINGO编程求解（程序见附录 D）,有成品一没有符合要求（3）的方案。成品二8根一捆，选出剩余的最长的8根总长度为11.5 6 10.5 9.5 84米，没有 达到最低要求88.5米，每捆减少1根就更不能满足要求，所以没有搭配方案。 成品三用lingo软件运行每捆总长度为88.589.5米，5根一捆的搭配，没有符合的 方案，再用4根一捆，结果可以搭配出6捆。在选择4捆时，具体的方案为表25成品三的原料按照要求（3）的优化方

28、案XiX23X24X25X26X27X28X29X30X31X32X33X34选取 根数000000000000XiX35X36X37X38X39X40X41X42X43X44X45X46选取 根数000030010000在剩余材料按要求（3）选才¥ 1捆时，方案为表26成品三的原料按照要求（3）的优化方案XiX23X24X25X26X27X28X29X30X31X32X33X34选取 根数000000000000XiX35X36X37X38X39X40X41X42X43X44X45X46选取 根数100000020000在剩余材料按要求（3）选才¥ 1捆时，方案为表27成

29、品三的原料按照要求（3）的优化方案XiX23X24X25X26X27X28X29X30X31X32X33X34选取 根数000000000000XiX35X36X37X38X39X40X41X42X43X44X45X46选取 根数010010100001七、模型的评价与推广.模型的评价7.1.1模型的优点(1)本模型的建立在相应的理论指导下，理论基础较为成熟，并且有相应 的专用软件支持，因此可信度较高。(2)利用数学工具，通过Lingo , Matlab编程的方法，严格地对模型进行 求解，具有较高的科学性。(3)模型的建立与实际紧密联系，充分考虑实际情况的多样性，从而使得 模型更贴进实际，通用

30、性、推广性较强。7.1.2模型的缺点(1)选取的数据难免在现实生活中存在诸多不协调，带来不可避免误差；(2)在对各档的数据进行处理时，都采用各档的标准值，这样会存在系统 误差。(3)在建模的过程中，虽然考虑了很多因素，但仍然忽略了其他不确定因 素等次要因素的影响，比如关于原料切割损坏问题等等，在模型假设里都忽略了 模型的推广该模型是一个典型的整数线性规划模型，在实际生活中有着广泛的使用空问。该模型不仅可以对天然肠衣搭配， 还可以对类似的问题进行求解。 以本模型 为基础可以进一步对复杂问题，用类似的方法对模型进行改进，因此模型具有较 强的普遍适用性。八、参考文献1万福永戴浩晖潘建瑜数学实验教程(

31、matlab版)北京科学出版社 20032谢金星、薛毅.优化建模与LINDO软件M.北京：清华大学出版社，20053傅远德.线性规划和整数规划M.成都：成都科技大学出版社，4卢开澄、卢华明.线性规划M.北京：清华大学出版社，2009附录附录A、B、C、D每种程序均附上一个模板，后面的采用递减重复运行即可附录A成品一max=y;3*x1+*x2+4*x3+*x4+5*x5+*x6+6*x7+*x8=89;x1+x2+x3+x4+x5+x6+x7+x8=20;y*x1<=43;y*x2<=59;y*x3<=39;y*x4<=41;y*x5<=27;y*x6<=2

32、8;y*x7<=34;y*x8<=21;gin(x1);bnd(0,x1,43);gin(x2);bnd(0,x2,59);gin(x3);bnd(0,x3,39);gin(x4);bnd(0,x4,41);gin(x5);bnd(0,x5,27);gin(x6);bnd(0,x6,28);gin(x7);bnd(0,x7,34);gin(x8);bnd(0,x8,21);gin(y);bnd(0,y,14);成品二max=y;7*x9+*x10+8*x11+*x12+9*x13+*x14+10*x15+*x16+11*x17+*x18+12*x19+*x20+13*x21+*x2

33、2=89;x9+x10+x11+x12+x13+x14+x15+x16+x17+x18+x19+x20+x21+x22=8;y*x9<=24;y*x10<=24;y*x11<=20;y*x12<=25;y*x13<=21;y*x14<=23;y*x15<=21;y*x16<=18;y*x17<=31;y*x18<=23;y*x19<=22;y*x20<=59;y*x21<=18;y*x22<=25;gin(x9);bnd(0,x9,24);gin(x10);bnd(0,x10,24);gin(x11);bnd(

34、0,x11,20);gin(x12);bnd(0,x12,25);gin(x13);bnd(0,x13,21);gin(x14);bnd(0,x14,23);gin(x15);bnd(0,x15,21); gin(x16);bnd(0,x16,18);gin(x17);bnd(0,x17,31); gin(x18);bnd(0,x18,23); gin(x19);bnd(0,x19,22);gin(x20);bnd(0,x20,59); gin(x21);bnd(0,x21,18);gin(x22);bnd(0,x22,25);gin(y);bnd(0,y,44);成品三max=y;14*x2

35、3+*x24+15*x25+*x26+16*x27+*x28+17*x29+*x30+18*x31+*x32+19*x33+*x34+20*x35+*x36+21*x37+*x38+22*x39+*x40+23*x41+*x42+24*x43+*x44+25*x45+*x46=89;x1+x2+x3+x4+x5+x6+x7+x8+x9+x10+x11+x12+x13+x14+x15+x16+x17+x18+x19+x20+ x21+x22+x23+x24=5;y*x23<=35;y*x24<=29;y*x25<=30;y*x26<=42;y*x27<=28;y*x

36、28<=42;y*x29<=45;y*x30<=49;y*x31<=50;y*x32<=64;y*x33<=52;y*x34<=63;y*x35<=49;y*x36<=35;y*x37<=27;y*x38<=16;y*x39<=12;y*x40<=2;y*x41<=0;y*x42<=6;y*x43<=0;y*x44<=0;y*x45<=0;y*x46<=1;gin(x23);bnd(0,x23,35);gin(x24);bnd(0,x24,29);gin(x25);bnd(0,x2

37、5,30);gin(x26);bnd(0,x26,42);gin(x27);bnd(0,x27,28);gin(x28);bnd(0,x28,42);gin(x29);bnd(0,x29,45);gin(x30);bnd(0,x30,49);gin(x31);bnd(0,x31,50);gin(x32);bnd(0,x32,64); gin(x33);bnd(0,x33,52);gin(x34);bnd(0,x34,63);gin(x35);bnd(0,x35,49);gin(x36);bnd(0,x36,35);gin(x37);bnd(0,x37,27);gin(x38);bnd(0,x3

38、8,16);gin(x39);bnd(0,x39,12);gin(x40);bnd(0,x40,2);gin(x41);bnd(0,x41,0);gin(x42);bnd(0,x42,6);gin(x43);bnd(0,x43,0);gin(x44);bnd(0,x44,0);gin(x45);bnd(0,x45,0);gin(x46);bnd(0,x46,1);gin(y);bnd(0,y,136);附录 B成品一function yy1c=1;for x1=0:43/11for x2=0:59/11for x3=0:39/11for x4=0:41/11for x5=0:27/11for

39、x6=0:28/11for x7=0:34/11for x8=0:21/11k1=x1*3+x2*+x3*4+x4*+x5*5+x6*+x7*6+x8*;k2=x1+x2+x3+x4+x5+x6+x7+x8;if (k1=89)&(k2=20)fprintf(' 第%d 种n %d,%d,%d,%d,%d,%d,%d,%d n',c,x1,x2,x3,x4,x5,x6,x7,x8)c=c+1;endendendendendendendendend成品二function yy2c=1;for x9=0:24/22for x10=0:24/22for x11=0:20/22

40、for x12=0:25/22for x13=0:21/22for x14=0:23/22for x15=0:21/22for x16=0:18/22for x17=0:31/22for x18=0:23/22for x19=0:22/22for x20=0:59/22for x21=0:18/22for x22=0:25/22k1=7*x9+*x10+8*x11+*x12+9*x13+*x14+10*x15+*x16+11*x17+*x18+12*x19+*x20+13*x21+*x22;k2=x9+x10+x11+x12+x13+x14+x15+x16+x17+x18+x19+x20+x2

41、1+x22;if(k1=89)&(k2=8)fprintf('第 %d 种n %d,%d,%d,%d,%d,%d,%d,%d,%d,%d,%d,%d,%d,%dn',c,x9,x10,x11,x12,x13,x14,x15,x16,x17,x18,x19,x20,x21,x22)c=c+1;endendendendendendendendendendendendendendend成品三function yy3c=1;for x23=0:35/42for x24=0:29/42for x25=0:30/42for x26=0:4242for x27=0:2842for x

42、28=0:4242for x29=0:4542for x30=0:49/42for x31=0:50/42for x32=0:6442for x33=0:5242for x34=0:6342for x35=0:49/42for x36=0:3542forx37=0:2742for x38=0:1642forx39=0:1242for x40=0:242for x41=0:0/42for x42=0:842for x43=0:0/42for x44=0:0/42for x45=0:0/42for x46=0:“42k1=14*x23+*x24+15*x25+*x26+16*x27+*x28+17

43、*x29+*x30+18*x31+*x32+19*x33+*x34+20*x35+*x36+21*x37+*x38+22*x39+*x40+23*x41+*x42+24*x43+*x44+25*x4 5+*x46;k2=x23+x24+x25+x26+x27+x28+x29+x30+x31+x32+x33+x34+x35+x36+x37+x38+x39+ x40+x41+x42+x43+x44+x45+x46;if(k1=89)&(k2=5)fprintf('第d种n %d,%d,%d,%d,%d,%d,%d,%d,%d,%d,%d,%d,%d,%d,%d,%d,%d,%d,%

44、d,%d,%d,%d, %d,%dn',c,x23,x22,x25,x26,x27,x28,x29,x30,x31,x32,x33,x34,x35,x36,x37,x38,x39,x40,x41, x42,x43,x44,x45,x46)c=c+1;endendendendendendendendendendendendendendendendendendendendendendendendend附录 C产品一max=y;3*x1+*x2+4*x3+*x4+5*x5+*x6+6*x7+*x8=89;x1+x2+x3+x4+x5+x6+x7+x8=20;y*x1<=43;y*x2&

45、lt;=59;y*x3<=39;y*x4<=41;y*x5<=27;y*x6<=28;y*x7<=34;y*x8<=21;gin(x1);bnd(0,x1,43);gin(x2);bnd(0,x2,59);gin(x3);bnd(0,x3,39);gin(x4);bnd(0,x4,41);gin(x5);bnd(0,x5,27);gin(x6);bnd(0,x6,28);gin(x7);bnd(0,x7,34);gin(x8);bnd(0,x8,21);gin(y);bnd(0,y,14);产品二max=y;7*x9+*x10+8*x11+*x12+9*x1

46、3+*x14+10*x15+*x16+11*x17+*x18+12*x19+*x20+13*x21+*x22=89;x9+x10+x11+x12+x13+x14+x15+x16+x17+x18+x19+x20+x21+x22=8;y*x9<=24;y*x10<=24;y*x11<=20;y*x12<=25;y*x13<=21;y*x14<=23;y*x15<=21;y*x16<=18;y*x17<=31;y*x18<=23;y*x19<=22;y*x20<=59;y*x21<=18;y*x22<=25;gin(

47、x9);bnd(0,x9,24);gin(x10);bnd(0,x10,24);gin(x11);bnd(0,x11,20);gin(x12);bnd(0,x12,25);gin(x13);bnd(0,x13,21);gin(x14);bnd(0,x14,23);gin(x15);bnd(0,x15,21); gin(x16);bnd(0,x16,18);gin(x17);bnd(0,x17,31); gin(x18);bnd(0,x18,23); gin(x19);bnd(0,x19,22);gin(x20);bnd(0,x20,59); gin(x21);bnd(0,x21,18);gin(x22);bnd(0,x22,25);gin(y);bnd(0,y,44);产品三max=y;14*x23+*x24+15*x25+*x26+16*x27+*x28+17*x29+*x30+18*x31+*x32+19*x33+*x34+20*x35+*x36+21*x37+*x38+22*x39+*x40+23*x41+*x42+24*x43+*x44+25*x45+* x46=89