天然肠衣搭配问题经典

池州学院
天然肠衣搭配问题
组员：陈强
赵晋彪
赵海龙
目录
一、问题重述 (3)
1.1问题背景 (3)
1.2.问题条件 (4)
1.3.问题要求 (4)
1.4需要解决的问题 (5)
二、问题分析 (5)
三、模型假设 (6)
四、符号说明 (6)
五模型的建立 (6)
5.1、模型建立 (6)
5.2、根据要求模型建立 (9)
六、模型求解 (10)
6.1、问题要求（1）模型求解 (10)
6.2、问题要求（2）模型求解 (12)
6.3、问题要求（3）模型求解 (15)
七、模型的评价与推广 (16)
7.1.模型的评价 (16)
7.1.1模型的优点 (16)
7.1.2模型的缺点 (17)
7.2模型的推广 (17)
八、参考文献 (17)
附录 (17)
附录A (17)
附录B (19)
附录C (23)
附录D (25)
天然肠衣搭配问题
摘要
天然肠衣制作加工是我国的一个传统产业，出口量占世界首位，而天然肠衣传统的人工生产方式已不能满足出口量日益增长的需要。

因此，我们从节约生产成本、提高企业生产效率的角度出发，我们结合原料的供给量、长度及成品规格等约束条件进行了模型设计。

根据题目中的表1中的成品的规格和表2中的原料，我们所需要解决的问题有：如何搭配才能使得成品的捆数最多？对于针对这一个问题我们采用线性规划建立模型并利用MATLAB以捆数相同，最短长度越长越好的原则，求得模型的最优解。

另外，由于所有的原料按长度分档，通常以0.5米为一档，如：3-3.4米按3米计算，3.5米-3.9米按3.5米计算，其余的依此类推。

表1是几种常见成品的规格，长度单位为米，∞表示没有上限，但实际长度小于26米。

再把不同档次的原料按照不同的规格进行搭配，分别搭配成三种规格的成品，依次是成品一(3—6.5米,20根，总长度89米)，成品二(7—13.5米，8根，总长度89米)，成品三（14—∞米，5根，总长度89米）。

运用线性规划分别对成品一、成品二、成品三建立模型，利用LINGO编程进行1步，2步，3步……优化筛选，得出方案。

并且，对各步筛选所剩余原料再进行优化得出方案，另外，为了提高原材料的使用率，每成品的总长度允许有0.5米的误差，总根数允许比标准少1根，某种成品对应得原材料有剩余，可以降一级使用，这样就会出现每捆总长度88.5米和89.5米，有19根一捆，7根一捆，4根一捆，在满足条件时，计算出最大捆数。

关键词：天然肠衣；线性规划；MATLAB ;LINGO
一、问题重述
1.1问题背景
天然肠衣制作加工是我国的一个传统产业，出口量占世界首位，而天然肠衣传统的生产依靠工人，边丈量长度，边心算，将原材料按指定根数和总长度组装成捆方式已不能满足出口量日益增长的需要。

因此，我们从节约生产成本、提高企业生产效率的角度出发，结合原料的供给量、长度及成品规格等约束条件进行生产模型设计。

1.2.问题条件
已知某天然肠衣组装的规格表( 见表1) 和某批次原料的描述表( 见表2) ．
表2 原料描述表
1.3.问题要求
根据生产规格和原料描述表，设计满足以下要求及允许条件下的组装方案( 即对原料进行打捆的搭配方案) 。

1) 对于给定的一批原料，装出的成品时，捆数越多越好；
2) 对于成品捆数相同的方案，最短长度最长的成品越多，方案越好； 3) 允许成品总长度有± 0.5 m 的误差，即总长度88.5米，,89.5米； 4) 每捆总根数允许比标准少1 根即每捆可以为19根，7根，4根； 5) 某种成品对应原料如果出现剩余，可以降级使；。

6）为了食品安全要求在30分钟内完成。

1.4需要解决的问题
首先我们要考虑原料应该怎样处理，然后考虑各类成品原料的搭配问题以及剩余材料的得搭配问题。

二、问题分析
天然肠衣经过人工清理和截取成不等长度，在进行搭配成成品，传统工艺
里，工人们边劳作用心丈量，其操作不仅效率差，有时会造成严重资源浪费，为了解决这一问题即做到同等材料中，加工成成品后，剩余材料最少，捆数最多，我们依据表一，表二的生产规格和所拥有的原料进行生产搭配成成品，作出最有模型。

对于问题要求一，我们运用lingo 软件编程求出优化方案，并作出逐步优化，得出最大捆数；对于问题二，利用matlab 编程，按照最短长度最长的原则选出最优方案，将选出的方案与lingo 结合编程，得出最大捆数，然后在剩余中再按照最短长度最长的原则选择最后所能组成的捆数；对于问题要求三、四、五，对于成品一规格，要求长度1l 满足
m l m 5.63≤≤，[]l l =或[]5.011+=l l 根数1x 满足20=i x 根； 对于成品二规格，要求
长度2l 满足5.1472≤≤l m ，[]22l l =或[]22l l =+0.5,根数2x 满足82=x 根；对于成品三规格，要求长度3l 满足14m ≤3l ≤26m ，[]33l l =或3L =[3L ]+0.5,第三类原料最大长度为∞，但原料的实际长度小于26米,根数3x 满足5根。

在每个档次的原料满足这些约束条件的前提下，建立不定方程组，用lingo 软件求解出每种规格产品的最大捆数；为提高原料使用率，总长度允许有±0.5米的误差，总根数允许比标准少1根。

因此生产条件将变为每捆天然肠衣成品的总长度为89±0.5m ，同时可能会出现成品一每捆十九根，成品二每捆七根，成品三每捆四根，建立不定方程组，用lingo 软件求解出每种规格产品的最大捆数。

三、模型假设
1、假设题目所给的数据真实可靠；
2、所有原料都是新鲜的，没有腐烂；
3、所有原料都是同种原料，不需要对原料进行分类；
4、厂房的环境、车间和工作人员的卫生等不会对原料产生影响；
5、厂房中的温度，光照等不会对原料的保鲜产生影响。

四、符号说明
(1,2,3...46)i x i = 每个档次在组装中所使用的根数 (1,2,3...)i y i = 不同规格成品按某种方案所组装出来的捆数
(1,2,3)i z i = 每种规格所能够组合出来的最大捆数
Z 各规格的最大捆数的总和
(1,2,3...46)s i x i = 各个档次在经过缩小取值范围后的取值
(1,2,3,1,2,3...)j i y i j == 规格i 的成品第j 次按照某个方案组装所的捆数
y i x 各个档次经过组装的剩余量
i n 档次i 的原料开始时的数量
(1,2,3)i L i = 每种规格所能取得的最大理想捆数
(1,2,3)i S i = 规格i 的原料在经过要求（2）的筛选后剩余的总根数
五 模型的建立
5.1、模型建立
成品一
成品一要求是选出20根不同肠衣组装成1捆，原料长度按0.5米为一档，如：3~3.4按3米计算，所选出得肠衣总长度之和为89米，可以得出
()[]89
15.0320
8
1
87654321=-+=+++++++∑=i i x i x x x x x x x x
其中81x x −→−
的限制条件为： 431≤x 592≤x
493≤x 414≤x 275≤x
286≤x
347≤x 288≤x
按照某个方案组装(1,2,3...)i y i =捆之后，18~x x 的根数会减少一部分，减少后限制条件变化为
123443594941i i i
i x y x y x y x y ≤-≤-≤-≤-
567827283421i
i i
i
x y x y x y x y ≤-≤-≤-≤-
经过k 次组装之后，规格1的成品的捆数要求最大，可以得出
11n
i i Maxz y ==∑
成品二
规格2的要求是选出8根不同肠衣组装成1捆，所选出得肠衣总长度之和为89米，可以得出
()[]89
95.078
229
22
9=-+=∑∑==i
i i i
x
i x
922~x x 的根数限制条件为
910111213141524242025212321
x x x x x x x ≤≤≤≤≤≤≤ 16171819202122
18312322591825
x x x x x x x
≤≤≤≤≤≤≤
按照某个方案组装(1,2,3...)i y i =捆之后，922~x x 的根数会减少一部分，减少后限制条件变化为
910111213141524242025212321i i i i i i i
x y x y x y x y x y x y x y ≤-≤-≤-≤-≤-≤-≤- 1617181920
2122
18312322591825i i i i i i
i x y x y x y x y x y x y x
y ≤-≤-≤-≤-≤-≤-≤-
经过k 次组装之后，规格2成品的捆数要求最大，可以得出
21n
i i Maxz y ==∑
成品三
成品3要求是选出5根不同肠衣组装成1捆，，所选出得肠衣总长度之和为89米，可以得出
46
23
232425262728465
1414.51515.51616.5 (2689)
i
i x
x x x x x x x ==++++++=∑
2346~x x 的根数限制条件为
232425262728293031323334352930422842454950645263
x x x x x x x x x x x x ≤≤≤≤≤≤≤≤≤≤≤≤
35363738394041424344454649352716122060001
x x x x x x x x x x x x ≤≤≤≤≤≤≤≤≤≤≤≤
按照某个方案组装(1,2,3...)i y i =捆之后，2346~x x 的根数会减少一部分，减少后限制条件变化为
232425262728293031323334352930422842454950645263i i i i i i i i i i i i
x y x y x y x y x y x y x y x y x y x y x y x y ≤-≤-≤-≤-≤-≤-≤-≤-≤-≤-≤-≤-
35363738394041424344454649352716122060001i i i i i i i i i
x y x y x y x y x y x y x y x y x x x x y ≤-≤-≤-≤-≤-≤-≤-≤-≤≤≤≤-
经过k 次组装之后，成品三成品的产品的捆数要求最大，可以得出
31
n
i i Maxz y ==∑
最后得出成品的最大捆数为
321z z z Z ++=
5.2、根据要求模型建立
成品一
在第一轮筛选后的余量总根数（即满足问题要求2）
115221112S =+++++=
由于题目要求是20跟一捆，即便充分利用问题要求（3）和（4），也达不到题中的要求，所以成品一不能再搭配出多余的捆数。

成品二
同样经过要求(2)的筛选后剩下的总根数
22422032156182S =+++++++=
根数可以达到8根一捆，再考虑一捆的总长度为88.58889.5、、米和每捆为8根，8根不行，再考虑每捆7根。

成品三
也是经过要求(2)筛选后剩下的总根数为
3111211226127S =++++++++= 按题中5根一捆，根数能够满足要求，再加上每捆总长度在895.0±米，每捆按5或求其次4根计算。

目标函数：max Y 约束条件：
46
23
[140.5(23)]89i
i i x
=+-=∑
46
23
5i
i x
==∑
230yx ≤ 352yx ≤ 240yx ≤ 360yx ≤ 250yx ≤ 370yx ≤ 260yx ≤ 381yx ≤
271yx ≤ 3912yx ≤ 280yx ≤ 402yx ≤
291yx ≤ 410yx ≤ 300yx ≤ 426yx ≤ 310yx ≤ 430yx ≤ 320yx ≤ 440yx ≤ 330yx ≤ 450yx ≤
341yx ≤ 461yx ≤
六、模型求解
6.1、问题要求（1）模型求解
对于不同规格的成品，会存在一个理想的最大捆数S，但是由于实际条件的约束，总捆数不可能超过S，其中理想的最大捆数等于此种规格的总根数之和除于此种规格每捆的根数要求
成品一的最大理想捆数
1435939 (3421)
14
20
L
++++
==
余下12根
成品二的最大理想捆数
2242420 (1825)
44
8
L
++++
==
余下2根
成品三的最大理想捆数
3352930 (01)
135
5
L
++++
==
余下2根
成品一
利用LINGO编程（程序见附录A），对成品一进行最优化筛选，可得方案
经过1轮的优化之后，并未达到最优，再进行优化，有方案
表4 筛成品一的原料第2筛选的方案
成品一的实际的总捆数之和
111314
Maxz=+=，实际的总捆数之和
1
z与最
大理想总捆数
1
S相等，成品一的优化结束。

各个档次剩余根数分别为
表5 成品一的原料经过筛选所剩余的原料
成品二
利用同样的方法对成品二进行优化分配处理
表6 成品二的原料第1轮筛选的方案
表7 成品二的原料第2轮筛选的方案
表8 成品二的原料第3轮筛选的方案
经3轮的筛选之后，各个档次的根数剩余量为
表9成品二的原料3轮筛选的剩余量
料所能够组装的最大实际捆数为
2229334
Maxz=++=
成品三
成品三的原料采用同样方法，经过九轮筛选之后剩余原料跟数为：
求得成品三的实际最大捆数为
342+30+26+14+9+4+2+1+1=129
Maxz=
所以，利用LINGO对题目中给定的原料求解出最大能够组装出的捆数
1231434129177
MaxZ Maxz Maxz Maxz
=++=++=
6.2、问题要求（2）模型求解
成品一
利用MATLAB编程（程序见附录B），按照最短长度最长的原则选出如下方案，将选出的方案结合LINGO编程（程序见附录C），得出捆数为11
表11 成品一的原料按照最短长度最长优化的方案
表12 成品一的原料按照最短长度最长优化的方案
在符合成品二的材料中按照最短长度最长的方案选择22捆，方案如下
表13 成品二的原料按照最短长度最长优化的方案
表14 成品二的原料按照最短长度最长优化的方案
表15成品二的原料按照最短长度最长优化的方案
在符合成品三的材料中按照最短长度最长的方案进行选择42捆，方案如下
表18 成品三的原料按照最短长度最长优化的方案
表24 成品三的原料按照最短长度最长优化的方案
成品最短长度最长，此时各个规格的捆数之和为178，相比要求1所得的捆数多出了1捆，这是因为成品三在选择最短长度最长的方案组装的时候使得材料能够更加充分利用
6.3、问题要求（3）模型求解
通过LINGO编程求解（程序见附录D），有
成品一
没有符合要求(3)的方案。

成品二
8根一捆，选出剩余的最长的8根总长度为11.5610.59.584
+⨯+=米，没有达到最低要求88.5米，每捆减少1根就更不能满足要求，所以没有搭配方案。

成品三
用lingo软件运行每捆总长度为88.5~89.5米，5根一捆的搭配，没有符合的方案，再用4根一捆，结果可以搭配出6捆。

在选择4捆时，具体的方案为
七、模型的评价与推广
7.1.模型的评价
7.1.1模型的优点
（1）本模型的建立在相应的理论指导下，理论基础较为成熟，并且有相应的专用软件支持，因此可信度较高。

（2）利用数学工具，通过Lingo，Matlab编程的方法，严格地对模型进行求解，具有较高的科学性。

（3）模型的建立与实际紧密联系，充分考虑实际情况的多样性，从而使得模型更贴进实际，通用性、推广性较强。

7.1.2模型的缺点
（1）选取的数据难免在现实生活中存在诸多不协调，带来不可避免误差；
（2）在对各档的数据进行处理时，都采用各档的标准值，这样会存在系统误差。

（3）在建模的过程中，虽然考虑了很多因素，但仍然忽略了其他不确定因素等次要因素的影响，比如关于原料切割损坏问题等等，在模型假设里都忽略了
7.2模型的推广
该模型是一个典型的整数线性规划模型，在实际生活中有着广泛的使用空间。

该模型不仅可以对天然肠衣搭配，还可以对类似的问题进行求解。

以本模型为基础可以进一步对复杂问题，用类似的方法对模型进行改进，因此模型具有较强的普遍适用性。

八、参考文献
[1]万福永戴浩晖潘建瑜数学实验教程（matlab版）科学出版社2003
[2]谢金星、薛毅.优化建模与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*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+7.5*x10+8*x11+8.5*x12+9*x13+9.5*x14+10*x15+10.5*x16+11*x17+11.5*x 18+12*x19+12.5*x20+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);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+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=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*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;
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(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(x3 6);
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(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 yy1
c=1;
for x1=0:43/11
for x2=0:59/11
for x3=0:39/11
for x4=0:41/11
for x5=0:27/11
for x6=0:28/11
for x7=0:34/11
for x8=0:21/11
k1=x1*3+x2*3.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)&(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;
end
end
end
end
end
end
end
end
end
成品二
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 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)&(k2==8)
fprintf('第%d种\n %d,%d,%d,%d,%d,%d,%d,%d,%d,%d,%d,%d,%d,%d \n',c,x9,x10,x11,x12,x13,x14,x15,x16,x17,x18,x19,x20,x21,x22)
c=c+1;
end
end
end
end
end
end
end
end
end
end
end
end
end
end
end
成品三
function yy3
c=1;
for x23=0:35/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
.
kszl
.
for x40=0:2/42
for x41=0:0/42
for x42=0:6/42
for x43=0:0/42
for x44=0:0/42
for x45=0:0/42
for x46=0:1/42
k1=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+2 2.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+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,%d,%d,%d,%d ,%d,%d
\n',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;
end
end
end
end
end
end
.
kszl
.
end
end
end
end
end
end
end
end
end
end
end
end
end
end
end
end
end
end
end
附录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;
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*x 18+12*x19+12.5*x20+13*x21+13.5*x22=89;
.
kszl
.
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,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*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*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;
.
kszl
.
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(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(x3 6);
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(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);
附录D
max=y;
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=89;
x23+x24+x25+x26+x27+x28+x29+x30+x31+x32+x33+x34+x35+x36+x37+x38+x3 9+x40+x41+x42+x43+x44+x45+x46=4;
y*x23<=0;
y*x24<=0;
y*x25<=0;
y*x26<=0;
y*x27<=1;
y*x28<=0;
y*x29<=1;
y*x30<=0;
y*x31<=0;
y*x32<=0;
y*x33<=0;
y*x34<=1;
y*x35<=2;
y*x36<=0;
.
kszl
.
y*x37<=0;
y*x38<=1;
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,0);gin(x24);bnd(0,x24,0);gin(x25);bnd(0,x25,0);
gin(x26);bnd(0,x26,0);gin(x27);bnd(0,x27,1);
gin(x28);bnd(0,x28,0);gin(x29);bnd(0,x29,1);
gin(x30);bnd(0,x30,0);gin(x31);bnd(0,x31,0);gin(x32);bnd(0,x32,0);
gin(x33);bnd(0,x33,0);gin(x34);bnd(0,x34,1);gin(x35);bnd(0,x35,2);gin(x36); bnd(0,x36,0);gin(x37);bnd(0,x37,0);gin(x38);bnd(0,x38,1);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);
.
kszl。