运筹学课程设计报告(附代码)

《运筹学》课程设计报告

姓名：

班级：

学号：

一、问题描述

1、机型指派问题

机型指派优化设计是航空公司制定航班计划的重要内容，它要求在满足航班频率和时刻安排以及各机型飞机总数约束的条件下，将各机型飞机指派给相应的航班，使运行成本最小化。本课程设计要求建立机型指派问题的数学模型，应用优化软件Lindo/Lingo进行建模求解，给出决策建议，包括各机型执行的航班子集和相应的运行成本。

2、问题描述

已知某航空公司航班频率和时刻安排如《运筹学课程设计指导书》中表1所示，航班需求数据和运输距离如表2所示，其中，OrignA/P表示起飞机场，Dep.T.表示起飞时间，Dest.A/P表示目标机场，Dist表示轮挡距离，Demand表示航班需求量，Std Dev.表示需求的标准差。该航空公司的机队有两种机型：9架B737-800，座位数162；6架B757-200，座位数200。飞八个机场：A,B,I,J,L,M,O,S。

B737-800的CASM(座英里成本)是0.34元，B757-200是0.36元。两种机型的 RASM（座英里收益）都是 1.2元。以成本最小为目标进行机型指派，在成本方面不仅考虑运行成本，还必须考虑旅客溢出成本，否则将偏向于选取小飞机，使航空公司损失许多旅客。

旅客溢出成本是指旅客需求大于航班可提供座位数时，旅客流失到其他航空公司造成的损失。旅客需求服从N(μ,σ)的正态分布。如果机票推销工作做得好，溢出旅客并不全部损失，有部分溢出旅客将该成本航空公司其他航班，这种现象叫做“再获得”（Recapture）。设有15%的溢出旅客被再获得。

将飞机指派到航班上去，并使飞机总成本最小。

二、分析建模

1．确定决策变量

经过对问题描述的分析得出，要解决飞机机型指派问题，我设定了两类变量：

（1）针对各条航线的机型，令B737-800和B757-200分别为机型1和机型2，设变量Xi,j.其中101≤i≤142,j=1或2。且对于变量Xi,j=0或1，当Xi,j=1，表示第i条航线由第j 种飞机运营。例如，X101，1=1，则第101号航班由第1种机型飞行，且X101，2=0

（2）针对机场时间节点飞机流的变量，设变量Gm,j.表示对于第m个节点上第j种机型的数量，例如，G A1，1表示A机场第1个节点上第1种机型的数量。

2.目标函数

以飞机总成本最小为指派目标，而单个航班的飞机总成本包括两个部分：1.运输成本；2. 旅

客溢出成本；其中运输成本的表达式为：B737-800的架数*162*0.34*该航班的轮挡距离+ B757-200的架数*200*0.36*该航班的轮挡距离；旅客溢出成本的表达式为：航班旅客溢出的期望值*1.2*该航班的轮挡距离*0.85。详细计算公式如下：

（1）营运成本 B737-800 ：C1=

∑i

Xi i Dist S CASM 1,*)(** 101≤i ≤142

B757-200 ：C2=

∑i

Xi i Dist S CASM 2,*)(** 101≤i ≤142

CASM 表示飞机座英里成本，S 表示飞机座位数，Dist （i ）表示第i 条航线的轮挡距离

（2）旅客溢出成本

B737-800 ：C1’=旅客溢出数期望值*机票价格 =

∑-i

Xi i Dist RASM d E 1,*)(**)(%)151( 101≤i ≤142

B757-200 : C2’=旅客溢出数期望值*机票价格 =

∑-i

Xi i Dist RASM d E 2,*)(**)(%)151( 101≤i ≤142

RASM 表示飞机座英里收益

其中，对于两种机型的旅客溢出期望值 E （d ）=

⎰

∞

-c

d d d f c d )()(

=⎰

∞

--+σ

σμπ

/)()(2/1d c c x e -(x2/2)dx

μ表示航班需求量的期望，σ表示需求的标准差， c 表示飞机的座位数

（3）建立目标函数

Min C=C1+C2+C1’+C2’

2.时空网络建模及其约束条件 （1）节点飞机平衡条件

对于每种机型，在时空网络中各节点的飞机流必须保持平衡。如某机型有一定数量航班到达，一定数量航班出发，因此该节点后该机型留下飞机数=原有飞机数+到达飞机数-离开飞机数。下面会对各个机场的具体节点飞机流量状况进行解释说明。（已设定B737-800为机型1， B757-200为机型2）

如分析A机场的各机型飞机流量状况。根据节点平衡条件，节点A1的约束条件：GA1,1=GA6,1-X110,1 （或者为GA1,2=GA6,2-X110,2）。其中GA1,1代表该机场节点现存飞机数目，其中A代表机场，A1,1中前一个1代表机场A的第一个节点，第二个1代表机型1。X110,1中代表飞入或飞出飞机架数，只能为0或1，110代表航班代号，1代表第一种机型。以下约束条件具有相似的意义，将不作详细阐述。并且只详尽列出节点A的约束条件，其他节点的情况可以同理写出。

对于机型1而言有如下约束条件机场A

节点A2的约束条件：GA2,1=GA1,1+x131,1 Array节点A3的约束条件：GA3,1=GA2,1-x111,1

节点A4的约束条件：GA4,1=GA3,1+x132,1

节点A5的约束条件：GA5,1=GA4,1-x112,1

节点A6的约束条件：GA6,1=GA5,1+x133,1

对于机型2而言有如下约束条件

节点A2的约束条件：GA2,2=GA1,2+x131,2）

节点A3的约束条件：GA3,2=GA2,2+x111,2）

节点A4的约束条件：GA4,2=GA3,2+x132,2）

节点A5的约束条件：GA5,2=GA4,2+x112,2）

节点A6的约束条件：GA6,2=GA5,2+x133,2）

机场B机场I