《运筹学》课程设计报告Word文件下载.docx
- 文档编号:3769597
- 上传时间:2023-05-02
- 格式:DOCX
- 页数:21
- 大小:19.40KB
《运筹学》课程设计报告Word文件下载.docx
《《运筹学》课程设计报告Word文件下载.docx》由会员分享,可在线阅读,更多相关《《运筹学》课程设计报告Word文件下载.docx(21页珍藏版)》请在冰点文库上搜索。
131
111
110
18:
00
17:
05
13:
10
12:
20:
35
8:
A
节点A2的约束条件:
GA2,1=GA1,1+x131,1(或者为GA2,2=GA1,2+x131,2)
节点A3的约束条件:
GA3,1=GA2,1-x111,1(或者为GA3,2=GA2,2+x111,2)
节点A4的约束条件:
GA4,1=GA3,1+x132,1(或者为GA4,2=GA3,2+x132,2)
节点A5的约束条件:
GA5,1=GA4,1-x112,1(或者为GA5,2=GA4,2+x112,2)
节点A6的约束条件:
GA6,1=GA5,1+x133,1(或者为GA6,2=GA5,2+x133,2)
BILM
115
135
114
134
113
124
123
103
122
102
101
21:
30
15:
20
0:
15
14:
9:
45
5:
136
137
138
139
6:
23:
10:
141
142
119
120
121
116
140
117
118
J
130
106
16:
127
109
55
50
19:
105
108
104
107
126
129
128
125
40
25
11:
7:
OS
1:
(2)飞机总数的约束
每基地机场各机型的过夜飞机之和不超过该机型飞机的总数,对于B737-800的数目,有以下约束条件:
GA6,1+GB6,1+GI6,1+GL6,1+GM6,1+GO6,1+GS6,1+GJ39,1<
=9
各项的含义在时空网络建模中已经标明。
同理对于B757-200的数目,有以下约束条件:
GA6,2+GB6,2+GI6,2+GL6,2+GM6,2+GO6,2+GS6,2+GJ39,2<
=6
(3)每条航线飞机数的限制
2
∑Xi,k=1其中i代表航线,如101;
k代表机型,只能是1和2。
具体表达如:
X101,1+X101,2=1,并且X101,1和X101,2只能一个取0,一个取1。
k=1
三.模型求解
model:
sets:
file/@ole('
file.xls'
'
flightNo'
)/:
dist,demand,stddev,type1,type2;
!
引用电子表格,其属性的意义依次是距离,需求,标准差,机型1,机型2;
airport/1..8/;
机场,共有8个;
airline/1..42/;
航班,共有42个;
derive(airline,airport):
M,N;
派生类,将机场和航班组合成一个二维集合;
endsets
data:
dist,demand,stddev=@ole('
);
将电子表格里的数据对应赋给相映的变量;
enddata
min=@sum(airline(i):
dist(i)*0.34*162*type1(i)+dist(i)*0.36*200*type2(i)+0.85*dist(i)*type1(i)*1.2*stddev(i)*@psl
((162-demand(i))/stddev(i))+0.85*dist(i)*type2(i)*1.2*stddev(i)*@psl((200-demand(i))/stddev(i)));
!
目标函数,以成本最小为目标,成本包括运输成本和旅客溢出成本;
@for(airline(i):
@bin(type1(i)));
一个航班机型1的飞机只能取0或1架;
@bin(type2(i)));
一个航班机型2的飞机只能取0或1架;
type1(i)+type2(i)=1);
一个航班两种机型和只能是1;
M(1,1)=M(6,1)-type1(10);
节点约束条件,具体意义见实验报告;
M(2,1)=M(1,1)+type1(31);
M(3,1)=M(2,1)-type1(11);
M(4,1)=M(3,1)+type1(32);
M(5,1)=M(4,1)-type1(12);
M(6,1)=M(5,1)+type1(33);
M(1,2)=M(6,2)-type1(16);
M(2,2)=M(1,2)+type1(37);
M(3,2)=M(2,2)-type1(17);
M(4,2)=M(3,2)+type1(38);
M(5,2)=M(4,2)-type1(18);
M(6,2)=M(5,2)+type1(39);
M(1,3)=M(6,3)+type1(40);
M(2,3)=M(1,3)-type1(19);
M(3,3)=M(2,3)+type1(41);
M(4,3)=M(3,3)-type1(20);
M(5,3)=M(4,3)+type1(42);
M(6,3)=M(5,3)-type1(21);
M(1,4)=M(39,4)-type1(40);
M(2,4)=M(1,4)-type1(25);
M(3,4)=M(2,4)-type1(22);
M(4,4)=M(3,4)-type1(37);
M(5,4)=M(4,4)+type1(16);
M(6,4)=M(5,4)+type1(19);
M(7,4)=M(6,4)+type1(7)-type1(31);
M(8,4)=M(7,4)-type1(28);
M(9,4)=M(8,4)+type1
(1);
M(10,4)=M(9,4)+type1(4)-type1(34);
M(11,4)=M(10,4)+type1(10);
M(12,4)=M(11,4)+type1(17);
M(13,4)=M(12,4)-type1(41);
M(14,4)=M(13,4)+type1(13);
M(15,4)=M(14,4)-type1(38);
M(16,4)=M(15,4)+type1(8);
M(17,4)=M(16,4)-type1(32);
M(18,4)=M(17,4)-type1(29);
M(19,4)=M(18,4)-type1(35);
M(20,4)=M(19,4)+type1
(2)-type1(42);
M(21,4)=M(20,4)+type1(5);
M(22,4)=M(21,4)+type1(20);
M(23,4)=M(22,4)-type1(26);
M(24,4)=M(23,4)+type1(11);
M(25,4)=M(24,4)-type1(23);
M(26,4)=M(25,4)+type1(18);
M(27,4)=M(26,4)+type1(14);
M(28,4)=M(27,4)-type1(33);
M(29,4)=M(28,4)-type1(36);
M(30,4)=M(29,4)-type1(24);
M(31,4)=M(30,4)+type1(9);
M(32,4)=M(31,4)+type1(21);
M(33,4)=M(32,4)-type1(27);
M(34,4)=M(33,4)+type1(12);
M(35,4)=M(34,4)+type1(3);
M(36,4)=M(35,4)+type1(6);
M(37,4)=M(36,4)-type1(30);
M(38,4)=M(37,4)+type1(15);
M(39,4)=M(38,4)-type1(39);
M(1,5)=M(6,5)-type1
(1);
M(2,5)=M(1,5)-type1
(2);
M(3,5)=M(2,5)+type1(22);
M(4,5)=M(3,5)-type1(3);
M(5,5)=M(4,5)+type1(23);
M(6,5)=M(5,5)+type1(24);
M(1,6)=M(6,6)-type1(13);
M(2,6)=M(1,6)+type1(34);
M(3,6)=M(2,6)-type1(14);
M(4,6)=M(3,6)+type1(35);
M(5,6)=M(4,6)-type1(15);
M(6,6)=M(5,6)+type1(36);
M(1,7)=M(6,7)-type1(7);
M(2,7)=M(1,7)+type1(28);
M(3,7)=M(2,7)-type1(8);
M(4,7)=M(3,7)+type1(29);
M(5,7)=M(4,7)-type1(9);
M(6,7)=M(5,7)+type1(30);
M(1,8)=M(6,8)-type1(4);
M(2,8)=M(1,8)-type1(5);
M(3,8)=M(2,8)+type1(25);
M(4,8)=M(3,8)-type1(6);
M(5,8)=M(4,8)+type1(26);
M(6,8)=M(5,8)+type1(27);
M(6,1)+M(6,2)+M(6,3)+M(39,4)+M(6,5)+M(6,6)+M(6,7)+M(6,8)<
=9;
过夜飞机数约束条件,具体意义见实验报告;
N(1,1)=N(6,1)-type2(10);
N(2,1)=N(1,1)+type2(31);
N(3,1)=N(2,1)-type2(11);
N(4,1)=N(3,1)+type2(32);
N(5,1)=N(4,1)-type2(12);
N(6,1)=N(5,1)+type2(33);
N(1,2)=N(6,2)-type2(16);
N(2,2)=N(1,2)+type2(37);
N(3,2)=N(2,2)-type2(17);
N(4,2)=N(3,2)+type2(38);
N(5,2)=N(4,2)-type2(18);
N(6,2)=N(5,2)+type2(39);
N(1,3)=N(6,3)+type2(40);
N(2,3)=N(1,3)-type2(19);
N(3,3)=N(2,3)+type2(41);
N(4,3)=N(3,3)-type2(20);
N(5,3)=N(4,3)+type2(42);
N(6,3)=N(5,3)-type2(21);
N(1,4)=N(39,4)-type2(40);
N(2,4)=N(1,4)-type2(25);
N(3,4)=N(2,4)-type2(22);
N(4,4)=N(3,4)-type2(37);
N(5,4)=N(4,4)+type2(16);
N(6,4)=N(5,4)+type2(19);
N(7,4)=N(6,4)+type2(7)-type2(31);
N(8,4)=N(7,4)-type2(28);
N(9,4)=N(8,4)+type2
(1);
N(10,4)=N(9,4)+type2(4)-type2(34);
N(11,4)=N(10,4)+type2(10);
N(12,4)=N(11,4)+type2(17);
N(13,4)=N(12,4)-type2(41);
N(14,4)=N(13,4)+type2(13);
N(15,4)=N(14,4)-type2(38);
N(16,4)=N(15,4)+type2(8);
N(17,4)=N(16,4)-type2(32);
N(18,4)=N(17,4)-type2(29);
N(19,4)=N(18,4)-type2(35);
N(20,4)=N(19,4)+type2
(2)-type2(42);
N(21,4)=N(20,4)+type2(5);
N(22,4)=N(21,4)+type2(20);
N(23,4)=N(22,4)-type2(26);
N(24,4)=N(23,4)+type2(11);
N(25,4)=N(24,4)-type2(23);
N(26,4)=N(25,4)+type2(18);
N(27,4)=N(26,4)+type2(14);
N(28,4)=N(27,4)-type2(33);
N(29,4)=N(28,4)-type2(36);
N(30,4)=N(29,4)-type2(24);
N(31,4)=N(30,4)+type2(9);
N(32,4)=N(31,4)+type2(21);
N(33,4)=N(32,4)-type2(27);
N(34,4)=N(33,4)+type2(12);
N(35,4)=N(34,4)+type2(3);
N(36,4)=N(35,4)+type2(6);
N(37,4)=N(36,4)-type2(30);
N(38,4)=N(37,4)+type2(15);
N(39,4)=N(38,4)-type2(39);
N(1,5)=N(6,5)-type2
(1);
N(2,5)=N(1,5)-type2
(2);
N(3,5)=N(2,5)+type2(22);
N(4,5)=N(3,5)-type2(3);
N(5,5)=N(4,5)+type2(23);
N(6,5)=N(5,5)+type2(24);
N(1,6)=N(6,6)-type2(13);
N(2,6)=N(1,6)+type2(34);
N(3,6)=N(2,6)-type2(14);
N(4,6)=N(3,6)+type2(35);
N(5,6)=N(4,6)-type2(15);
N(6,6)=N(5,6)+type2(36);
N(1,7)=N(6,7)-type2(7);
N(2,7)=N(1,7)+type2(28);
N(3,7)=N(2,7)-type2(8);
N(4,7)=N(3,7)+type2(29);
N(5,7)=N(4,7)-type2(9);
N(6,7)=N(5,7)+type2(30);
N(1,8)=N(6,8)-type2(4);
N(2,8)=N(1,8)-type2(5);
N(3,8)=N(2,8)+type2(25);
N(4,8)=N(3,8)-type2(6);
N(5,8)=N(4,8)+type2(26);
N(6,8)=N(5,8)+type2(27);
N(6,1)+N(6,2)+N(6,3)+N(39,4)+N(6,5)+N(6,6)+N(6,7)+N(6,8)<
=6;
end
四.结果分析
1.通过约束条件对程序运行结果进行检验。
例如,各个机场的机型1过夜航班之和是否符合约束条件:
M(6,1)+M(6,2)+M(6,3)+M(40,4)+M(6,5)+M(6,6)+M(6,7)+M(6,8)<
过夜飞机数约束条件,具体意义见实验报告
过夜航班(前面的数字代表节点,后面的代表机场)
飞机架数
M(6,1)0.000000
M(6,2)1.000000
M(6,3)0.000000
M(39,4)3.000000
M(6,5)2.000000
M(6,6)0.000000
M(6,7)1.000000
M(6,8)2.000000
各项之和为9,符合该条件。
又有约束条件如下,为各个机场的机型2过夜航班之和
N(6,1)+N(6,2)+N(6,3)+N(40,4)+N(6,5)+N(6,6)+N(6,7)+N(6,8)<
对于该约束条件,有以下测试结果:
N(6,1)1.000000
N(6,2)0.000000
N(6,3)0.000000
N(39,4)2.000000
N(6,5)0.000000
N(6,6)2.000000
N(6,7)0.000000
N(6,8)0.000000
各项之和为5,小于6,符合条件。
2.各个航线的机型指派如下:
航班与对应的机型
TYPE1(101)1.000000
TYPE1(102)1.000000
TYPE1(103)1.000000
TYPE1(104)1.000000
TYPE1(105)1.000000
TYPE1(106)1.000000
TYPE1(107)1.000000
TYPE1(108)1.000000
TYPE1(109)1.000000
TYPE1(110)0.000000
TYPE1(111)0.000000
TYPE1(112)1.000000
TYPE1(113)0.000000
TYPE1(114)0.000000
TYPE1(115)1.000000
TYPE1(116)1.000000
TYPE1(117)0.000000
TYPE1(118)1.000000
TYPE1(119)1.000000
TYPE1(120)1.000000
TYPE1(121)1.000000
TYPE1(122)1.000000
TYPE1(123)1.000000
TYPE1(124)1.000000
TYPE1(125)1.000000
TYPE1(126)1.000000
TYPE1(127)1.000000
TYPE1(128)1.000000
TYPE1(129)1.000000
TYPE1(130)1.000000
TYPE1(131)0.000000
TYPE1(132)1.000000
TYPE1(133)0.000000
TYPE1(134)1.000000
TYPE1(135)0.000000
TYPE1(136)0.000000
TYPE1(137)0.000000
TYPE1(138)1.000000
TYPE1(139)1.000000
TYPE1(140)1.000000
TYPE1(141)1.000000
TYPE1(142)1.000000
TYPE2(101)0.000000
TYPE2(102)0.000000
TYPE2(103)0.000000
TYPE2(104)0.000000
TYPE2(105)0.000000
TYPE2(106)0.000000
TYPE2(107)0.00000
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 运筹学 课程设计 报告