天然肠衣搭配问题经典.doc

1、池州学院天然肠衣搭配问题组员：陈强 赵晋彪 赵海龙目录一、问题重述 41.1 问题背景 41.2.问题条件 41.3.问题要求 51.4 需要解决的问题 5二、问题分析 5三、模型假设 6四、符号说明 6五 模型的建立 65.1、模型建立 .65.2、根据要求模型建立 .9六、模型求解 106.1、问题要求（1）模型求解 .106.2、问题要求（2）模型求解 .126.3、问题要求（3）模型求解 .15七、模型的评价与推广 177.1.模型的评价 177.1.1 模型的优点 .177.1.2 模型的缺点 .177.2 模型的推广 17八、参考文献 17附录 18附录 A18附录 B20附录 C

2、23附录 D25天然肠衣搭配问题摘要天然肠衣制作加工是我国的一个传统产业，出口量占世界首位，而天然肠衣传统的人工生产方式已不能满足出口量日益增长的需要。因此，我们从节约生产成本、提高企业生产效率的角度出发，我们结合原料的供给量、长度及成品规格等约束条件进行了模型设计。根据题目中的表 1 中的成品的规格和表 2 中的原料，我们所需要解决的问题有：如何搭配才能使得成品的捆数最多？对于针对这一个问题我们采用线性规划建立模型并利用 MATLAB 以捆数相同，最短长度越长越好的原则，求得模型的最优解。另外，由于所有的原料按长度分档，通常以0.5米为一档，如：3-3.4米按3米计算，3.5米-3.9米按3

3、.5米计算，其余的依此类推。表1是几种常见成品的规格，长度单位为米，表示没有上限，但实际长度小于26米。再把不同档次的原料按照不同的规格进行搭配，分别搭配成三种规格的成品，依次是成品一(36.5米,20根，总长度89米)，成品二(713.5米，8根，总长度89米)，成品三（14米，5根，总长度89米） 。运用线性规划分别对成品一、成品二、成品三建立模型，利用LINGO编程进行1步，2步，3步优化筛选，得出方案。并且，对各步筛选所剩余原料再进行优化得出方案，另外，为了提高原材料的使用率，每成品的总长度允许有0.5米的误差，总根数允许比标准少1根，某种成品对应得原材料有剩余，可以降一级使用，这样就

4、会出现每捆总长度88.5米和89.5米，有19根一捆，7根一捆，4根一捆，在满足条件时，计算出最大捆数。关键词：天然肠衣 ； 线性规划 ；MATLAB ;LINGO一、问题重述1.1问题背景天然肠衣制作加工是我国的一个传统产业，出口量占世界首位，而天然肠衣传统的生产依靠工人，边丈量长度，边心算，将原材料按指定根数和总长度组装成捆方式已不能满足出口量日益增长的需要。因此，我们从节约生产成本、提高企业生产效率的角度出发，结合原料的供给量、长度及成品规格等约束条件进行生产模型设计。1.2.问题条件已知某天然肠衣组装的规格表( 见表1) 和某批次原料的描述表( 见表2) 表1 成品规格表成品 最短长度

5、/m 最长长度/m 根数 总长度/m成品一成品二 成品三37146.513.52085898989表 2 原料描述表长度 3-3.4 3.5-3.9 4-4.4 4.5-4.9 5-5.4 5.5-5.9 6-6.4 6.5-6.9根数 43 59 39 41 27 28 34 21长度 7-7.4 7.5-7.9 8-8.4 8.5-8.9 9-9.4 9.5-9.910-10.410.5-10.9根数 24 24 20 25 21 23 21 18长度11-11.411.5-11.912-12.412.5-12.913-13.413.5-13.914-14.414.5-14.9根数 31

6、23 22 59 18 25 35 29长度15-15.415.5-15.916-16.416.5-16.917-17.417.5-17.918-18.418.5-18.9根数 30 42 28 42 45 49 50 64长度19-19.419.5-19.920-20.420.5-20.921-21.421.5-21.922-22.422.5-22.9根数 52 63 49 35 27 16 12 2长度23-23.423.5-23.924-24.424.5-24.925-25.425.5-25.9根数 0 6 0 0 0 11.3.问题要求根据生产规格和原料描述表，设计满足以下要求及允许条

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

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

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

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

11、大理想捆数(1,23)iL规格 的原料在经过要求（2）的筛选后剩余的总根数iSi五 模型的建立5.1、模型建立成品一成品一要求是选出 20 根不同肠衣组装成 1 捆，原料长度按 0.5 米为一档，如：33.4 按 3米计算，所选出得肠衣总长度之和为 89 米，可以得出8915.032081 876542ii xxx其中 的限制条件为：81x 431592x341x 275x86347x28按照某个方案组装 捆之后， 的根数会减少一部分，减少后限(,2.)iy18x制条件变化为123459iixy56782341iiyx经过 次组装之后，规格 1 的成品的捆数要求最大，可以得出k 1niMaxz

12、y成品二规格 2 的要求是选出 8 根不同肠衣组装成 1 捆，所选出得肠衣总长度之和为 89 米，可以得出 895.07829iii xx的根数限制条件为92x10213454xx1678192035xx按照某个方案组装 捆之后， 的根数会减少一部分，减少后限(,2.)iy92制条件变化为910213454iiiixyxy1678192035iiiixyxy经过 次组装之后，规格 2 成品的捆数要求最大，可以得出k 21niMaxzy成品三成品 3 要求是选出 5 根不同肠衣组装成 1 捆， ，所选出得肠衣总长度之和为89 米，可以得出 462324526728461.1.1.5.9iixxx

13、x的根数限制条件为2346x324562789301234904565xxx3567389401243546492601xxx按照某个方案组装 捆之后， 的根数会减少一部分，减少后(1,.)iy23限制条件变化为2345267829301234904565iiiiiixyxyxy3567389401243546492601iiiiixyxyxy经过 次组装之后，成品三成品的产品的捆数要求最大，可以得出k 31niMaxz最后得出成品的最大捆数为 32Z5.2、根据要求模型建立成品一在第一轮筛选后的余量总根数（即满足问题要求2） 15212S由于题目要求是20跟一捆，即便充分利用问题要求（3）和

14、（4），也达不到题中的要求，所以成品一不能再搭配出多余的捆数。成品二同样经过要求(2)的筛选后剩下的总根数 2403215682S根数可以达到 8 根一捆，再考虑一捆的总长度为 米和每捆为 8.9.5、 、根，8 根不行，再考虑每捆 7 根。成品三也是经过要求(2)筛选后剩下的总根数为 3121627S按题中 5 根一捆，根数能够满足要求，再加上每捆总长度在 89 米，每5.0捆按 5 或求其次 4 根计算。目标函数： maxY约束条件：462310.5(23)89iiix46235iixyxy240yx360yx5 76 8127139280yx40yx9 130 261 432yx0yx3

15、 54 461六、模型求解6.1、问题要求（1）模型求解对于不同规格的成品，会存在一个理想的最大捆数 ，但是由于实际条件S的约束，总捆数不可能超过 ，其中理想的最大捆数等于此种规格的总根数之S和除于此种规格每捆的根数要求成品一的最大理想捆数 14359.34210L余下 12 根成品二的最大理想捆数 2.854余下 2 根成品三的最大理想捆数 3590.135L余下 2 根成品一利用 编程（程序见附录 A） ，对成品一进行最优化筛选，可得方案LINGO表 3 成品一的原料第 1 轮筛选的方案ix12x4x56x78x选取根数2 5 3 3 2 1 3 1方案所得捆数：11经过 1 轮的优化之后

16、，并未达到最优，再进行优化，有方案表 4 筛成品一的原料第 2 筛选的方案ix12x34x56x78x选取根数7 1 2 1 1 5 0 3方案所得捆数：3成品一的实际的总捆数之和 ， 实际的总捆数之和 与最134Maxz1z大理想总捆数 相等，成品一的优化结束。各个档次剩余根数分别为1S表 5 成品一的原料经过筛选所剩余的原料ix12x34x56x78x剩余根数0 1 0 5 0 2 1 1成品二利用同样的方法对成品二进行优化分配处理表 6 成品二的原料第 1 轮筛选的方案ix910x112x1314x15选取根数 0 1 0 1 0 0 0i1671892212选取根数 0 1 1 1 2

17、 0 1方案所得捆数为：22表 7 成品二的原料第 2 轮筛选的方案 ix910x112x1314x15选取根数 0 0 0 0 0 2 1i16171819212选取根数 1 1 0 0 1 2 0方案所得捆数为：9表 8 成品二的原料第 3 轮筛选的方案ix910x112x1314x15选取根数 0 0 0 0 0 1 2ix1617x1819x2021x2选取根数 2 0 0 0 2 0 1方案所得捆数为：3经 3 轮的筛选之后，各个档次的根数剩余量为表 9 成品二的原料 3 轮筛选的剩余量 ix910x112x1314x15剩余根数 24 2 20 3 21 2 6i161718192

18、012剩余根数 3 0 1 0 0 0 0剩余的原料不能再组装出符合要求的成品，优化结束。最后得出成品二的原料所能够组装的最大实际捆数为 234Maxz成品三成品三的原料采用同样方法，经过九轮筛选之后剩余原料跟数为：表 10 成品三的原料 9 轮筛选的剩余量ix234x256x278x230x132x34x剩余根数0 0 0 0 1 0 1 0 0 0 0 1ix356x378x3940x142x34x546x剩余根数2 0 0 1 12 2 0 6 0 0 0 1经过 9 轮的优化筛选之后，剩余的原料无法再组装出符合要求的成品，最后求得成品三的实际最大捆数为 34+6194+1=29Maxz

19、所以，利用 对题目中给定的原料求解出最大能够组装出的捆数LINGO1237Zzax6.2、问题要求（2）模型求解成品一利用 编程（程序见附录 B） ，按照最短长度最长的原则选出如下方MATLB案，将选出的方案结合 编程（程序见附录 C） ，得出捆数为 11INGO表 11 成品一的原料按照最短长度最长优化的方案ix12x34x56x78x选取根 2 5 3 3 2 1 3 1数在剩余中再按照最短长度最长的方案选择 3 捆，方案如下表 12 成品一的原料按照最短长度最长优化的方案ix12x34x56x78x选取根数7 1 2 1 1 5 0 3成品二在符合成品二的材料中按照最短长度最长的方案选择

20、 22 捆，方案如下表 13 成品二的原料按照最短长度最长优化的方案ix910x112x1314x15选取根数 0 1 0 1 0 0 0i1671892212选取根数 0 1 1 1 2 0 1在剩余中再按照最短长度最长的方案选择 9 捆，方案如下表 14 成品二的原料按照最短长度最长优化的方案ix910x112x1314x15选取根数 0 0 0 0 0 2 1i16171819212选取根数 1 1 0 0 1 2 0再在剩余中再按照最短长度最长的方案选择 3 捆，方案如下表 15 成品二的原料按照最短长度最长优化的方案ix910x112x1314x15选取根数 0 0 0 0 0 1

21、2i1617181922选取根数 2 0 0 0 2 0 1成品三在符合成品三的材料中按照最短长度最长的方案进行选择 42 捆，方案如下表 16 成品三的原料按照最短长度最长优化的方案ix234x256x278x2930x132x34x选取根数0 0 0 0 0 1 1 0 1 1 1 0ix356x378x3940x142x34x546x选取根数0 0 0 0 0 1 0 0 0 0 0 0在剩余中再按照最短长度最长的方案选择 31 捆，方案如下表 17 成品三的原料按照最短长度最长优化的方案ix234x256x278x2930x132x34x选取根数1 0 0 1 0 0 0 0 0 0

22、0 2ix356x378x3940x142x34x546x选取根数0 1 0 0 0 0 0 0 0 0 0 0在剩余中再按照最短长度最长的方案选择 27 捆，方案如下表 18 成品三的原料按照最短长度最长优化的方案ix234x256x278x2930x132x34x选取根数0 1 0 0 1 0 0 1 0 0 0 0ix356x378x3940x142x34x546x选取根数1 0 1 0 0 0 0 0 0 0 0 0在剩余中再按照最短长度最长的方案选择 15 捆，方案如下表 19 成品三的原料按照最短长度最长优化的方案ix234x256x278x2930x132x34x选取根数0 0

23、2 0 0 0 0 1 0 0 0 0ix356x378x3940x142x34x546x选取根数1 0 0 1 0 0 0 0 0 0 0 0在剩余中再按照最短长度最长的方案选择 7 捆，方案如下表 20 成品三的原料按照最短长度最长优化的方案ix234x256x278x2930x132x34x选取根数0 0 0 1 0 0 0 0 1 3 0 0ix356x378x3940x142x34x546x选取根数0 0 0 0 0 0 0 0 0 0 0 0在剩余中再按照最短长度最长的方案选择 4 捆，方案如下表 21 成品三的原料按照最短长度最长优化的方案ix234x256x278x2930x1

24、32x34x选取根数1 0 0 1 0 0 0 0 0 0 1 0ix356x378x3940x142x34x546x选取根数1 1 0 0 0 0 0 0 0 0 0 0在剩余中再按照最短长度最长的方案选择 2 捆，方案如下表 22 成品三的原料按照最短长度最长优化的方案ix234x256x278x2930x132x34x选取根数0 1 0 0 0 0 0 1 0 0 3 0ix356x378x3940x142x34x546x选取根数0 0 0 0 0 0 0 0 0 0 0 0在剩余中再按照最短长度最长的方案选择 1 捆，方案如下表 23 成品三的原料按照最短长度最长优化的方案ix234x

25、256x278x2930x132x34x选取根数0 0 0 0 0 0 0 3 1 1 0 0ix356x378x3940x142x34x546x选取根数0 0 0 0 0 0 0 0 0 0 0 0在剩余中再按照最短长度最长的方案选择 1 捆，方案如下表 24 成品三的原料按照最短长度最长优化的方案ix234x256x278x2930x132x34x选取根数0 0 0 0 0 0 2 2 0 0 0 0ix356x378x3940x142x34x546x选取根数1 0 0 0 0 0 0 0 0 0 0 0按照上述成品三选择原材料进行组装，即可使得这批原材料的所组装出来的成品最短长度最长，此

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

27、要求(3)的优化方案ix234x256x278x2930x132x34x选取根数0 0 0 0 0 0 0 0 0 0 0 0ix356x378x3940x142x34x546x选取根数0 0 0 0 3 0 0 1 0 0 0 0在剩余材料按要求(3)选择 1 捆时，方案为表 26 成品三的原料按照要求(3)的优化方案ix234x256x278x2930x132x34x选取根数0 0 0 0 0 0 0 0 0 0 0 0ix356x378x3940x142x34x546x选取根数1 0 0 0 0 0 0 2 0 0 0 0在剩余材料按要求(3)选择 1 捆时，方案为表 27 成品三的原料

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

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

30、版）北京 科学出版社 20032谢金星、薛毅.优化建模与 LINDO 软件M.北京：清华大学出版社，2005 3傅远德.线性规划和整数规划M.成都：成都科技大学出版社，1989.12 4卢开澄、卢华明.线性规划M.北京：清华大学出版社，2009 附录附录 A、B、C、D 每种程序均附上一个模板，后面的采用递减重复运行即可附录 A成品一max=y;3*x1+3.5*x2+4*x3+4.5*x4+5*x5+5.5*x6+6*x7+6.5*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=28;

31、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+7.5*x10+8*x11+8.5*x12+9*x13+9.5*x14+10*x15+10.5*x16+11*x17+11.5*x18+12*x19+12.5*x20

32、+13*x21+13.5*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(x9);bnd(0,x9,24);gin(x10);bnd(0,x10,24);gin(x11);bnd(0,x11,20);gin(x12);bnd(0,x12,25);g

33、in(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+14.5*x24+15*x25+15.5*x26+16*

34、x27+16.5*x28+17*x29+17.5*x30+18*x31+18.5*x32+19*x33+19.5*x34+20*x35+20.5*x36+21*x37+21.5*x38+22*x39+22.5*x40+23*x41+23.5*x42+24*x43+24.5*x44+25*x45+25.5*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

35、*x28=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,x25,30);gin(x26);bnd(0,x26,42);gin(x27);bnd(0,x27,28); gin(

36、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,x38,16);gin(x39);bnd(0,x39,12);gin(x40);bnd(0,x40,2);gin(x

37、41);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/11 for x4=0:41/11 for x5=0:27/11for x6=0:28/11 for x7=0:34/11 for x8=0:21/11 k1=x1*3+x2*3.

38、5+x3*4+x4*4.5+x5*5+x6*5.5+x7*6+x8*6.5;k2=x1+x2+x3+x4+x5+x6+x7+x8;if (k1=89)endendendendendendendendend成品二function yy2 c=1;for x9=0:24/22 for x10=0:24/22 for x11=0:20/22 for x12=0:25/22 for x13=0:21/22 for x14=0:23/22 for x15=0:21/22 for x16=0:18/22 for x17=0:31/22 for x18=0:23/22 for x19=0:22/22 for

39、 x20=0:59/22 for x21=0:18/22 for x22=0:25/22 k1=7*x9+7.5*x10+8*x11+8.5*x12+9*x13+9.5*x14+10*x15+10.5*x16+11*x17+11.5*x18+12*x19+12.5*x20+13*x21+13.5*x22;k2=x9+x10+x11+x12+x13+x14+x15+x16+x17+x18+x19+x20+x21+x22;if (k1=89) endendendendendendendendendendendendendendend成品三function yy3 c=1; for x23=0:35

40、/42 for x24=0:29/42 for x25=0:30/42 for x26=0:42/42 for x27=0:28/42 for x28=0:42/42 for x29=0:45/42 for x30=0:49/42 for x31=0:50/42 for x32=0:64/42 for x33=0:52/42 for x34=0:63/42 for x35=0:49/42 for x36=0:35/42 for x37=0:27/42 for x38=0:16/42 for x39=0:12/42 for x40=0:2/42 for x41=0:0/42 for x42=0:

41、6/42 for x43=0:0/42 for x44=0:0/42 for x45=0:0/42 for x46=0:1/42k1=14*x23+14.5*x24+15*x25+15.5*x26+16*x27+16.5*x28+17*x29+17.5*x30+18*x31+18.5*x32+19*x33+19.5*x34+20*x35+20.5*x36+21*x37+21.5*x38+22*x39+22.5*x40+23*x41+23.5*x42+24*x43+24.5*x44+25*x45+25.5*x46;k2=x23+x24+x25+x26+x27+x28+x29+x30+x31+x3

42、2+x33+x34+x35+x36+x37+x38+x39+x40+x41+x42+x43+x44+x45+x46;if (k1=89)endendendendendendendendendendendendendendendendendendendendendendendendend附录 C产品一max=y;3*x1+3.5*x2+4*x3+4.5*x4+5*x5+5.5*x6+6*x7+6.5*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=28;y*x7=34;y*x8=21;gi

43、n(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+7.5*x10+8*x11+8.5*x12+9*x13+9.5*x14+10*x15+10.5*x16+11*x17+11.5*x18+12*x19+12.5*x20+13*x21+13.5*x22=8

44、9;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(0,x11,20);gin(x12);bnd(0,x12,25);gin(x13);bnd(0,x13,

45、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+14.5*x24+15*x25+15.5*x26+16*x27+16.5*x28+17*x2

46、9+17.5*x30+18*x31+18.5*x32+19*x33+19.5*x34+20*x35+20.5*x36+21*x37+21.5*x38+22*x39+22.5*x40+23*x41+23.5*x42+24*x43+24.5*x44+25*x45+25.5*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*x28=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;