经济均衡问题及其应用 (三）：交通流均衡问题

浅夏110

例4  某地有如图2所示的一个公路网，每天上班时间有6千辆小汽车要从居民区 A 前往工作区D。经过长期观察，我们得到了图中5条道路上每辆汽车的平均行驶时间和 汽车流量之间的关系，如表4所示，那么，长期来看，这些汽车将如何在每条道路上分布？



+ p9 \, Q& L+ j
7 Y8 `; _3 R% f4 x/ _* U6 i$ g
（1）问题分析

: t2 Y9 z* m. P这个问题看起来似乎与前面几个例子完全不同，但实际上交通流与市场经济活动类似，也存在着均衡。

我们可以想象有一个协调者，正如前面几个例子中的所谓中间商可以理解为市场规律一样，实际上这里的所谓协调者也可以认为是交通流的规律。交通流的规律就是每辆 汽车都将选择使自己从 A到D运行时间少的路线，其必然的结果是无论走哪条路线 从 A到D，终花费的时间应该是一样的（否则,花费时间较长的那条线路上的部分 汽车就会改变自己的路线，以缩短自己的行驶时间）。
/ N! P( Y) y4 _. O- ^
也就是说，长期来看，这些汽车在每条道路上的分布将达到均衡状态（所谓均衡， 这里的含义就是每辆汽车都不能仅仅通过自身独自改变道路，节省其行驶时间）。在这种想法下，我们来建立线性规划模型。
0 _- C, p& S* L( q5 c
8 q; K+ T( j* N3 x/ @. `（2）优化模型

交通流的规律要求所有道路上的流量达到均衡，我们仍然类似例1和例2来考虑问 题。如果车流量是一辆车一辆车增加的，那么在每条道路上车流量小于2时，车流量会有一个分布规律；当某条道路上流量正好超过2时，新加入的一辆车需要选择使自己堵 塞时间短的道路。这就提示我们把同一条道路上的流量分布分解成不同性质的三个部分。也就是说，我们用 Y(AB) 表示道路 AB 上的总的流量，并进一步把它分解成三部分：

i）道路 AB 上的流量不超过2时的流量，用X (2, AB)表示；

: K% R5 F+ e% Z( g0 P4 sii）道路 AB 上的流量超过2但不超过3时，超过2的流量部分用 X (3, AB) 表示；
, ]% K, a8 o- [3 x( A+ G+ j6 c
" t* \0 o1 \/ f/ s9 x. T7 ]iii）道路 AB 上的流量超过3但不超过4时，超过3的流量部分用 X (4, AB) 表示。

' S6 z2 _# j2 s3 v依次类推，对道路AC,BC,BD,CD , ,, 上同理可以定义类似的决策变量。因此，问题中总共有20个决策变量  和   

+ x. l) ?6 T( a9 {1 b问题的目标应当是使总的堵塞时间小。用  表示流量 对应的堵塞时间（ 即表3中的数据，是对每辆车而言的），我们看看用     [ j为道路] 作为总堵塞时间是否合适。很容易理解：后面加入道路的车辆可能又会造成前面进入道路的车辆的进一步堵塞，如流量为3时，原先流量为2的车辆实际上也只能按 的时间通过，而不是 。也就是说，  并不是总堵塞时间。但是我们也可以发现 关于i是单调增加的，即不断增加的车流只会使以前的堵塞加剧而不可能使以前的堵塞减缓。所以，关于决策变量 而言，     [ j为道路]  与我们希望优化的目标的单调性是一致的。因此，可以用     [ j为道路] 作为目标函数进行优化。
6 Q0 S- ?, n; i/ |4 a$ @2 M
约束条件有三类：
3 z5 R) N* M: p8 _
i）每条道路上的总流量Y 等于该道路上的分流量 X 的和；
' y! _, s9 ~* E6 {
( A5 ?! n5 i+ sii）道路交汇处A,B,C,D（一般称为节点）的流量守恒（即进入量等于流出量）；

iii）决策变量的上限限制，如 等。

8 }3 r: [7 E( L 于是对应的优化模型很容易直接写出（略）。
. @5 p0 c0 X; g5 \
（3）模型求解

编写LINGO程序如下：

; ~+ Z7 U) {. j: lMODEL:
TITLE 交通流均衡;
; U# D4 t. P$ A. P/ Z7 \SETS:  
    ROAD/AB,AC,BC,BD,CD/:Y;  
8 t, N# _, Y' t) x6 a  Q    CAR/2,3,4/;  
: {* O( l% D6 t: ^+ U8 [( Q( m+ m    LINK(CAR,ROAD): T, X;
ENDSETS
DATA:  
! 行驶时间(分钟) ;  
; I+ b1 p, _/ @4 |6 `% d9 a+ {5 ~    T=20,52,12,52,20        
      30,53,13,53,30   
; Z: P! u  U) C      40,54,14,54,40;
ENDDATA
[OBJ] MIN=@SUM(LINK: T*X); ! 目标函数;
) G9 e5 E/ `, ^& Y1 Y! 四个节点的流量守恒条件;
3 f6 B) p& J2 Z8 r[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));
) {: _( c+ C, W$ d" j6 b9 f2 p[NODE_D] Y(@INDEX(BD))+Y(@INDEX(CD))=6;
+ N2 ^- T% y! T4 Y2 I( \! 每条道路上的总流量Y等于该道路上的分流量X的和;
@FOR( ROAD(I): [ROAD_LIM] @SUM(CAR(J): X(J,I)) = Y(I));
/ Q& T- P1 E1 s) |/ F/ O3 c  t! 每条道路的分流量X的上下界设定;
@FOR(LINK(I,J)|I#EQ#1: @BND(0,X(I,J),2) );
@FOR(LINK(I,J)|I#GT#1: @BND(0,X(I,J),1) );
END
: h2 i# ?1 Y" Q  L- R. r! B: z' V可以指出的是，上面4个节点的流量守恒条件中，其实只有3个是独立的（也就是说，第4个条件总可以从其它3个方程推导出来），因此从中去掉任何一个都不会影响到计算结果。
0 P6 y1 ?* E- b9 d+ @/ L/ V2 X
. w' M# Z. r5 H% D) r3 h/ X（4）结果解释

LINGO的运行结果表明，均衡时道路 AB,AC,BC,BD,CD的流量分别是4，2，2，2，4（千辆）车。但是要注意，正如我们建立目标函数时所讨论过的，这时得到的目 标函数值452并不是真正的总运行和堵塞时间，而是一个用来表示目标函数趋势的虚拟 的量，没有太多实际物理意义。事实上，可以求出这时的真正运行时间是：每辆车通过AB,AC,BC,BD,CD 道路分别需要40，52，12，52，40（min），也就是在图中三条路线 ABD, ACD,ABCD 上都需要92min，所以这也说明交通流确实达到了均衡。
于是，均衡时真正的总运行时间应该是 6  × 92 =552（千辆车·min）。

（5）模型讨论
7 q( m3 X. W( M8 t: V: Y
; l& }4 H. L# T仔细想想就会发现，上面的解并不是最优解，即均衡解并不一定是最优的流量分配方案。为了求出使所有汽车的总运行时间小的交通流，应该如何做呢？也就是说，这相当于假设有一个权威的机构来统筹安排，最优地分配这些交通流，而不是像求均衡解 时那样认为各个个体（每辆车）都可以自己选择道路，自然达到平衡状态。
* s3 K) F- F' ~# q
为了进行统筹规划，我们需要把新增的流量   造成的实际堵塞时间计算出来（仍按每辆车计算），而不是像上面那样不考虑对原有车流造成的堵塞效应。以道路 AB 为例。

* v, }4 b0 ?* T6 Hi）当流量为2千辆时，每辆车的通过时间为20min，所以总通过时间是40（千辆 车·min）；

ii）当流量增加一个单位（本题中一个单位就是1千辆）达到3千辆时，每辆车的 通过时间为30min，所以总通过时间是90（千辆车·min）；

iii）当流量再增加一个单位达到4千辆时，每辆车的通过时间为40min，所以总 通过时间是160（千辆车·min）。
1 |  y7 Y0 ]* C. }: a6 m
) Q' H8 B* D  F# u由此可见，流量超过2而不超过3时，单位流量的增加导致的总通过时间的变化为 90－40＝50（千辆车·min）；流量超过3而不超过4时，单位流量的增加导致的总通 过时间的变化为160－90＝70（千辆车·min）。 类似地，对所有道路，都可以得到单位流量的增加导致总行驶时间的增量和汽车流 量之间的关系（参加表5）。
; e) ^; X* W2 h" W6 O! w$ b5 y( k
5 P6 R" ~5 t' ?/ Y6 `

用表5中的总行驶时间的增量数据代替前面模型中的每辆车的行驶时间数据 , 模型的其它部分完全不用变。重新求解LINGO模型，LINGO程序如下：

MODEL:
TITLE 交通流均衡;
& i6 n5 D2 p1 JSETS:  
- F" B  v3 U: v% ?    ROAD/AB,AC,BC,BD,CD/:Y;  
    CAR/2,3,4/;
    LINK(CAR,ROAD): T, X;
ENDSETS
DATA:  
" p' T3 R3 S7 ^/ N! 行驶时间(分钟) ;  
+ a5 ]3 d5 t  F5 v: @4 ~. Y! cT= 20 52 12 52 20
- ^5 z# n8 H4 I/ t0 `1 x6 S* T   50 55 15 55 50
   70 57 17 57 70 ;
ENDDATA
[OBJ] MIN=@SUM(LINK: T*X); ! 目标函数;
! 四个节点的流量守恒条件;
[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));
( h+ I8 c4 `) x1 H[NODE_D] Y(@INDEX(BD))+Y(@INDEX(CD))=6;
! 每条道路上的总流量Y等于该道路上的分流量X的和;
@FOR( ROAD(I): [ROAD_LIM] @SUM(CAR(J): X(J,I)) = Y(I));
* |/ Z# M0 G$ A! 每条道路的分流量X的上下界设定;
! b# b( k  u: N$ N) Y$ ^) t@FOR(LINK(I,J)|I#EQ#1: @BND(0,X(I,J),2) );
- Q7 n6 T( ]6 s( u: |& q5 X+ y" j@FOR(LINK(I,J)|I#GT#1: @BND(0,X(I,J),1) );
END
求得的最优车流分配方式是：道路 AB,AC,BD,CD的流量都是3千辆，而道路BC上没有流量；总（加权）运行时间为498（千辆车·min），优于均衡时的结果552（千 辆车·min）。此时，每辆车的运行时间＝498/6=83（min），少于均衡时的92min。 当然，这个最优解必须强制执行，否则 AB 道路上的一些车到底B 点时，发现当前走 BCD的时间只需要 42 3012 = + （min），比走BD的时间（53min）短很多，所以 他们就会改走BCD，导致走BCD的时间（主要是走道路CD的时间）增加；如此下 去，最后终将到达前面我们得到的均衡状态。
* q! P- [/ v0 b: l; |# A$ {0 S8 e
这是一个非常有趣的结果：当一个系统中的每个个体都独自追求个体利益大化 时，整体的利益却没有达到最大化。 更令人惊讶的是：这个例子的道路网中如果没有道路BC ，从 A到D的平均时间 是83min；而新开了一条道路BC 以后，从 A到D的平均时间居然变成92min，不是 加快反而减慢了。由此也可以理解，做出一个科学、合理的交通网的规划是一件相当复杂的工作。
————————————————
3 o0 L1 R; x2 J版权声明：本文为CSDN博主「wamg潇潇」的原创文章，遵循CC 4.0 BY-SA版权协议，转载请附上原文出处链接及本声明。
原文链接：https://blog.csdn.net/qq_29831163/java/article/details/89404182