数学建模D题天然肠衣搭配优化问题答案.docx
《数学建模D题天然肠衣搭配优化问题答案.docx》由会员分享,可在线阅读,更多相关《数学建模D题天然肠衣搭配优化问题答案.docx(17页珍藏版)》请在冰豆网上搜索。
![数学建模D题天然肠衣搭配优化问题答案.docx](https://file1.bdocx.com/fileroot1/2023-2/11/82a4ab95-ea09-4d39-9a97-3ffbc41dfbd8/82a4ab95-ea09-4d39-9a97-3ffbc41dfbd81.gif)
数学建模D题天然肠衣搭配优化问题答案
2011高教社杯全国大学生数学建模竞赛题目
D题天然肠衣搭配问题
摘要
该题主要研究生产天然肠衣及其搭配问题,并且要求在一定的原料情况下,生产的成品捆数越多越好,该问题属于线性规划并且为取整线性规划来求最优解问题。
根据每种规格的规定,在解题的过程中,我们建立线性方程组作为第一层优化,然后将建立的模型带入到lingo软件中,得到第一层优化最优方案,之后又根据实际进行了第二层优化,得到规格一成品捆数的上限为15捆;规格二成品的捆数的上限为37捆;规格三成品的捆数的上限为137捆;总捆数为188捆。
在一定的误差允许范围内,该方案较符合题目所属要求和实际生产情况。
并且生产后的剩余废弃原料少,做到了在限定原料内创造最大利润的好处。
问题简述:
原料按长度分档,通常以0.5米为一档,如:
3-3.4米按3米计算,3.5米-3.9米按3.5米计算,其余的依此类推。
成品规格和原料描述如图所示:
表1成品规格表
最短长度
最大长度
根数
总长度
3
6.5
20
89
7
13.5
8
89
14
∞
5
89
表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.9
10-10.4
10.5-10.9
根数
24
24
20
25
21
23
21
18
长度
11-11.4
11.5-11.9
12-12.4
12.5-12.9
13-13.4
13.5-13.9
14-14.4
14.5-14.9
根数
31
23
22
59
18
25
35
29
长度
15-15.4
15.5-15.9
16-16.4
16.5-16.9
17-17.4
17.5-17.9
18-18.4
18.5-18.9
根数
30
42
28
42
45
49
50
64
长度
19-19.4
19.5-19.9
20-20.4
20.5-20.9
21-21.4
21.5-21.9
22-22.4
22.5-22.9
根数
52
63
49
35
27
16
12
2
长度
23-23.4
23.5-23.9
24-24.4
24.5-24.9
25-25.4
25.5-25.9
根数
0
6
0
0
0
1
本题要求建立数学模型设计一个原料搭配方案,按题中所给规格完成原料搭配方案,并符合如下要求:
(1)对于给定的一批原料,装出的成品捆数越多越好;
(2)对于成品捆数相同的方案,最短长度最长的成品越多,方案越好;
(3)为提高原料使用率,总长度允许有±0.5米的误差,总根数允许比标准少1根;
(4)某种规格对应原料如果出现剩余,可以降级使用。
如长度为14米的原料可以和长度介于7-13.5米的进行捆扎,成品属于7-13.5米的规格;
(5)为了食品保鲜,要求在30分钟内产生方案。
模型的假设:
1、肠衣经过清洗整理后被分割成长度不等的小段(原料),原料在组装过程中长度不发生变化;
2、原料按长度分档,分档后原料不可再被分割;
3、将原料长度视为离散变量;
4、为提高原料使用率,每捆总长度允许有±0.5米的误差,每规格的成品总根数允许比标准少一根。
问题分析:
天然肠衣由于规定的档次(长度)不同,规格也不一样,所以每个规格的每捆肠衣成品长度不同,考虑到要在相同的成品捆数方案里找出最短长度最长的方案,我们想到了整数规划问题[1]的解决办法。
我们首先把肠衣成品的分配问题分开考虑,按下表中的成品规格表的规格将原料分成三类,即:
长度分布在3~6.5米的原料为规格一;长度分布在7~13.5米的原料为规格二;长度分布在14~25.5米的原料为规格三。
每种规格需要满足表中的根数约束,总长度约束,各区间总根数约束及整数约束。
表3成品规格表
规格
最短长度
最大长度
根数
总长度
1
3
6.5
20
89
2
7
13.5
8
89
3
14
∞
5
89
模型建立与求解:
第一层优化
符号声明:
代表三种成品的捆数(取整);
代表从第
个区间取得条数;
代表从第
个区间肠衣的长度,如3-3.4米按3米计算,3.5米-3.9米按3.5米计算,其余的依此类推;
为第
个区间总条数。
输入Lingo求得:
理论上,根据原料总根数和总长度以及每捆成品的根数和总长度,可求得规格一成品捆数的上限为14捆;规格二成品的捆数的上限为37捆;规格三成品的捆数的上限为137捆;总捆数为188捆。
结果分析:
第二层优化
表4原料剩余表
长度
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
剩余根数
0
0
16.625
0
0
0
0
0
长度
7-7.4
7.5-7.9
8-8.4
8.5-8.9
9-9.4
9.5-9.9
10-10.4
10.5-10.9
剩余根数
24
24
10
0
0
0
0
0
长度
11-11.4
11.5-11.9
12-12.4
12.5-12.9
13-13.4
13.5-13.9
14-14.4
14.5-14.9
剩余根数
0
0
0
0
0
0
0
0
长度
15-15.4
15.5-15.9
16-16.4
16.5-16.9
17-17.4
17.5-17.9
18-18.4
18.5-18.9
剩余根数
0
0
0
0
0
0
0
0
长度
19-19.4
19.5-19.9
20-20.4
20.5-20.9
21-21.4
21.5-21.9
22-22.4
22.5-22.9
剩余根数
1.84211
0
0
0
0
0
0
0
长度
23-23.4
23.5-23.9
24-24.4
24.5-24.9
25-25.4
25.5-25.9
剩余根数
0
0
0
0
0
0
根据某种规格对应原料如果出现剩余,可以降级使用的原则。
将4-4.4剩余的15根与7-7.4的4根组成一捆规格一,所以经过第二层优化后,规格一15捆,规格二37捆,规格三137捆,共189捆。
模型稳定性分析
我们所建立的模型通过对目标的最优化问题,使得多目标的规划问题转化为单目标线性规划问题,所以能比较好的反映出各个目标函数的重要程度。
而且模型在计算中作了一些舍入和取整,不可避免的产生了一些误差,但是这些误差的是可以容忍的。
结论
此模型在一定的误差允许范围内,较符合题目所属要求和实际生产情况。
并且生产后的剩余废弃原料少,做到了在限定原料内创造最大利润的好处。
工人可以工人根据这个方案“照方抓药”进行生产,在一定程度上可提高生产效率。
并且此模型易于推广,只需稍加改动就可以推广到解决其他分类封装的问题上。
参考文献
[1]姜启源,《数学模型(第三版)》,北京:
高等教育出版社,2003
附:
1、lingo程序代码
model:
sets:
c/c1..c8/:
a1;
d/d1..d14/:
a2;
e/e1..e24/:
a3;
r/r1..r8/:
b1;
s/s1..s14/:
b2;
t/t1..t24/:
b3;
allowed(r):
q1;
allowed1(s):
q2;
allowed2(t):
q3;
endsets
max=x+y+z;
@for(r(i):
q1(i)<=b1(i));
@for(s(i):
q2(i)<=b2(i));
@for(t(i):
q3(i)<=b3(i));
@sum(r(i):
q1(i)*a1(i))>=88.5*x;
@sum(s(i):
q2(i)*a2(i))>=88.5*y;
@sum(t(i):
q3(i)*a3(i))>=88.5*z;
@sum(r(i):
q1(i)*a1(i))<=89.5*x;
@sum(s(i):
q2(i)*a2(i))<=89.5*y;
@sum(t(i):
q3(i)*a3(i))<=89.5*z;
@sum(r(i):
q1(i))>=19*x;
@sum(s(i):
q2(i))>=7*y;
@sum(t(i):
q3(i))>=4*z;
@sum(r(i):
q1(i))<=20*x;
@sum(s(i):
q2(i))<=8*y;
@sum(t(i):
q3(i))<=5*z;
@gin(x);
@gin(y);
@gin(z);
data:
a1=33.544.555.566.5;
a2=77.588.599.51010.51111.51212.51313.5;
a3=1414.51515.51616.51717.51818.51919.52020.52121.52222.52323.52424.52525.5;
b1=4359394127283421;
b2=2424202521232118312322591825;
b3=35293042284245495064526349352716122060001;
enddata
2、程序运算结果
Globaloptimalsolutionfound.
Objectivevalue:
188.0000
Objectivebound:
188.0000
Infeasibilities:
0.000000
Extendedsolversteps:
0
Totalsolveriterations:
48
VariableValueReducedCost
X14.00000-1.000000
Y37.00000-1.000000
Z137.0000-1.000000
A1(C1)3.0000000.000000
A1(C2)3.5000000.000000
A1(C3)4.0000000.000000
A1(C4)4.5000000.000000
A1(C5)5.0000000.000000
A1(C6)5.5000000.000000
A1(C7)6.0000000.000000
A1(C8)6.5000000.000000
A2(D1)7.0000000.000000
A2(D2)7.5000000.000000
A2(D3)8.0000000.000000
A2(D4)8.5000000.000000
A2(D5)9.0000000.000000
A2(D6)9.5000000.000000
A2(D7)10.000000.000000
A2(D8)10.500000.000000
A2(D9)11.000000.000000
A2(D10)11.500000.000000
A2(D11)12.000000.000000
A2(D12)12.500000.000000
A2(D13)13.000000.000000
A2(D14)13.500000.000000
A3(E1)14.000000.000000
A3(E2)14.500000.000000
A3(E3)15.000000.000000
A3(E4)15.500000.000000
A3(E5)16.000000.000000
A3(E6)16.500000.000000
A3(E7)17.000000.000000
A3(E8)17.500000.000000
A3(E9)18.000000.000000
A3(E10)18.500000.000000
A3(E11)19.000000.000000
A3(E12)19.500000.000000
A3(E13)20.000000.000000
A3(E14)20.500000.000000
A3(E15)21.000000.000000
A3(E16)21.500000.000000
A3(E17)22.000000.000000
A3(E18)22.500000.000000
A3(E19)23.000000.000000
A3(E20)23.500000.000000
A3(E21)24.000000.000000
A3(E22)24.500000.000000
A3(E23)25.000000.000000
A3(E24)25.500000.000000
B1(R1)43.000000.000000
B1(R2)59.000000.000000
B1(R3)39.000000.000000
B1(R4)41.000000.000000
B1(R5)27.000000.000000
B1(R6)28.000000.000000
B1(R7)34.000000.000000
B1(R8)21.000000.000000
B2(S1)24.000000.000000
B2(S2)24.000000.000000
B2(S3)20.000000.000000
B2(S4)25.000000.000000
B2(S5)21.000000.000000
B2(S6)23.000000.000000
B2(S7)21.000000.000000
B2(S8)18.000000.000000
B2(S9)31.000000.000000
B2(S10)23.000000.000000
B2(S11)22.000000.000000
B2(S12)59.000000.000000
B2(S13)18.000000.000000
B2(S14)25.000000.000000
B3(T1)35.000000.000000
B3(T2)29.000000.000000
B3(T3)30.000000.000000
B3(T4)42.000000.000000
B3(T5)28.000000.000000
B3(T6)42.000000.000000
B3(T7)45.000000.000000
B3(T8)49.000000.000000
B3(T9)50.000000.000000
B3(T10)64.000000.000000
B3(T11)52.000000.000000
B3(T12)63.000000.000000
B3(T13)49.000000.000000
B3(T14)35.000000.000000
B3(T15)27.000000.000000
B3(T16)16.000000.000000
B3(T17)12.000000.000000
B3(T18)2.0000000.000000
B3(T19)0.0000000.000000
B3(T20)6.0000000.000000
B3(T21)0.0000000.000000
B3(T22)0.0000000.000000
B3(T23)0.0000000.000000
B3(T24)1.0000000.000000
Q1(R1)43.000000.000000
Q1(R2)59.000000.000000
Q1(R3)22.375000.000000
Q1(R4)41.000000.000000
Q1(R5)27.000000.000000
Q1(R6)28.000000.000000
Q1(R7)34.000000.000000
Q1(R8)21.000000.000000
Q2(S1)0.0000000.000000
Q2(S2)0.0000000.000000
Q2(S3)10.000000.000000
Q2(S4)25.000000.000000
Q2(S5)21.000000.000000
Q2(S6)23.000000.000000
Q2(S7)21.000000.000000
Q2(S8)18.000000.000000
Q2(S9)31.000000.000000
Q2(S10)23.000000.000000
Q2(S11)22.000000.000000
Q2(S12)59.000000.000000
Q2(S13)18.000000.000000
Q2(S14)25.000000.000000
Q3(T1)35.000000.000000
Q3(T2)29.000000.000000
Q3(T3)30.000000.000000
Q3(T4)42.000000.000000
Q3(T5)28.000000.000000
Q3(T6)42.000000.000000
Q3(T7)45.000000.000000
Q3(T8)49.000000.000000
Q3(T9)50.000000.000000
Q3(T10)64.000000.000000
Q3(T11)50.157890.000000
Q3(T12)63.000000.000000
Q3(T13)49.000000.000000
Q3(T14)35.000000.000000
Q3(T15)27.000000.000000
Q3(T16)16.000000.000000
Q3(T17)12.000000.000000
Q3(T18)2.0000000.000000
Q3(T19)0.0000000.000000
Q3(T20)6.0000000.000000
Q3(T21)0.0000000.000000
Q3(T22)0.0000000.000000
Q3(T23)0.0000000.000000
Q3(T24)1.0000000.000000
RowSlackorSurplusDualPrice
1188.00001.000000
20.0000000.000000
30.0000000.000000
416.625000.000000
50.0000000.000000
60.0000000.000000
70.0000000.000000
80.0000000.000000
90.0000000.000000
1024.000000.000000
1124.000000.000000
1210.000000.000000
130.0000000.000000
140.0000000.000000
150.0000000.000000
160.0000000.000000
170.0000000.000000
180.0000000.000000
190.0000000.000000
200.0000000.000000
210.0000000.000000
220.0000000.000000
230.0000000.000000
240.0000000.000000
250.0000000.000000
260.0000000.000000
270.0000000.000000
280.0000000.000000
290.0000000.000000
300.0000000.000000
310.0000000.000000
320.0000000.000000
330.0000000.000000
341.8421050.000000
350.0000000.000000
360.0000000.000000
370.0000000.000000
380.0000000.000000
390.0000000.000000
400.0000000.000000
410.0000000.000000
420.0000000.000000
430.0000000.000000
440.0000000.000000
450.0000000.000000
460.0000000.000000
470.0000000.000000
480.0000000.000000
493.0000000.000000
500.0000000.000000
5114.000000.000000
5234.000000.000000
53137.00000.000000
549.3750000.000000
5537.000000.000000
56127.15790.000000
574.6250000.000000
580.0000000.000000
599.8421050.000000