资源描述:
《运筹学,课程设计,.doc》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、交通流均衡问题一题目分析1.1问题某地有如图1所示的一个公路网,每天上班时间有6千辆小汽车要从居民区A前往居民区D,经过长期观察,我们得到了图中5条道路上每辆汽车的品均行驶时间和汽车流量之间的关系,如表1-1所示,那么,长期来看,这些汽车将如何在每条道路上分布?BDAC居民区工作区图1一个公路网示意图表1-1平均行驶时间和汽车流量之间的关系道路ABACBCBDCD行驶时间/min流量≤220521252202<流量≤330531353303<流量≤440541454401.2问题分析这个问题我们可以想
2、象有一个协调组,实际上所谓的协调者可以认为是交通流规律,交通流规律就是每辆汽车都将选择使自己从A到D运行时间最少的路线,其必然结果是无论走那条路线从A到D,最终花费的时间应该是一样的(否则,花费时间较长的那条线路上的部分汽车就会改变自己的路线,以缩短自己的行驶时间).也就是说,长期来看,这些汽车在每条道路上的分布将达到均衡状态(所谓均衡状态,这里的含义就是每辆汽车都不能仅仅通过自身独自改变道路节省其行驶时间),在这种想法下,我们来建立线性规划模型。二总体设计交通流的规律要求所有道路上的流量达到均衡,如
3、果车流量是一辆车一辆车增加的,那么在每条道路上车流量少于2时,车流量会有一个分布规律;当某路道路上流量正好超过2时,新加入的一辆车需要选择使自己堵塞时间最短的道路,这就提示我们吧同一条道路上的流量分布分解成不同性质的三个部分,也就是说,我们用Y(AB)表示道路AB上的总的流量,并进一步吧他分解成三个部分;(1)道路AB上的流量不超过2时的流量,用X(2,AB)表示;(2)道路AB上的流量超过2但不超过3时,超过2的流量部分用X(3,AB)表示;(3)道路AB上的流量超过3但不超过4时,超过3的流量部分
4、用X(4,AB)表示;依此类推,对道路AC,BC,BD,CD上同理可以定义类似的决策变量。因此,问题中总共有20个决策变量Y(j)和X(i,j)(I=2,3,4;j=AB,AC,BC,BD,CD)。问题的目标应当是使总的堵塞时间最少。用T(i,j)表示流量X(i,j)对应的堵塞时间(即表1-1的数据,是对每辆车而言的),我们看看用T(i,j)X(i,j)作为总堵塞时间是否合适,很容易理解;后面加入道路的车辆可能又会造成前面进入道路的车辆的进一步堵塞,如流量为3时,原先流量为2的车辆实际上也只能按T(3
5、,j)的时间通过,而不是T(2,j),也就是说,T(i,j)X(i,j)并不是总的堵塞时间,但是我们也可以发现,T(i,j)关于I单调增加的,即不断增加的车流只会使以前的堵塞加剧而不可能使以前的堵塞减缓,所以,关于决策变量X(i,j)而言,T(i,j)X(i,j)与我们希望优化的目标的单调性是一致的。因此,可以用T(i,j)X(i,j)作为目标函数进行优化。约束条件有三类:(1)每条道路上的总流量Y等于该道路上的分流量X的和;(2)道路交汇初A,B,,C,D(一般称为节点)的流量守恒(即流入量等于流出
6、量);(3)决策变量的上限限制,如X(2,AB)≤2,X(3,AB)≤1,X(4,AB)≤1等。三程序这个问题可以用LINDO软件求解,不过这里用LINGO软件求解更方便些,LINGO模型如下:SETS:ROAD/AB,AC,BC,BD,CD/:Y;CAR/2,3,4/;LINK(CAR,ROAD):T,X;ENDSETSDATA:!行驶时间;T=20,52,12,52,2030,53,13,53,3040,54,14,54,40;ENDDATA[OBJ]MIN=@SUM(LINK:T*X);!目标函
7、数;!四个节点的流量守恒条件;[NODE_A]Y(@INDEX(AB))+Y(@INDEX(AC))=6;[NODE_B]Y(@INDEX(AB))=Y(@INDEX(BC))+Y(@INDEX(BD));[NODE_C]Y(@INDEX(AC))+Y(@INDEX(BC))=Y(@INDEX(CD));[NODE_D]Y(@INDEX(BD))+Y(@INDEX(CD))=6;!每条道路上的总流量Y等于该道路上的分流量X的和;@FOR(ROAD(I):[ROAD_LIM]@SUM(CAR(J):X(
8、J,I))=Y(I));!每条道路的分流量X的上下界设定;@FOR(LINK(I,J)
9、I#EQ#1:@BND(0,X(I,J),2));@FOR(LINK(I,J)
10、I#GT#1:@BND(0,X(I,J),1));END可以指出的是,上面四个节点的流量守恒条件条件中,其实只有3个是独立的(也就是说,第四个总可以从其他3个方程推导出来),因此从中去掉任何一个都不会影响到计算结果。求解这个模型,得到如下解答(其中略去已知参数T的显示结果):Globalo