在线时间 4 小时 最后登录 2017-2-18 注册时间 2009-7-27 听众数 5 收听数 0 能力 0 分 体力 36 点 威望 2 点 阅读权限 20 积分 53 相册 0 日志 1 记录 0 帖子 43 主题 7 精华 0 分享 0 好友 1
升级 50.53%
TA的每日心情 奋斗 2014-2-8 12:28
签到天数: 5 天
[LV.2]偶尔看看I
) {* o( e7 |4 K 本文就巴黎凯旋门环岛的交通问题,建立了控制进入、环绕、流出此环岛车流量的红绿灯模型,目的是使环岛内交通顺畅,并且尽量让堵车时间短,堵车数量少。 / W4 n, }: D8 L
通过分析,发现环岛内的最大车流量为 1000, 波动范围为 +200 ,还可根据车道宽度计算出每个路口的最大车流量。这两个因素对环岛交通有着很大影响。因此,主要考虑车流量和环岛内的车辆数目的影响。并设定,在建立模型时,环内车辆总数最好不超过 1000 辆。
$ q$ i6 R# n, h 根据各时段车流量的多少,本文将车流分布为四种情况:高峰期、次高峰、一般情况、稀疏情况。再根据各时期的车流量,建立了环岛内车辆总数 Q 关于流入量与流出量的方程: 。通过随机模拟,得出环岛 12 个路口的车流量,并根据堵车时间尽量短,堵车数量尽量少的原则,找出所有可能的红绿灯组合(前提是每个路口都有红绿灯),通过比较,得出最优化的组合(具体组合见模型建立与求解部分)。 1 H3 c, k2 I4 R
通过随机模拟,对于不同时期,得到不同最佳方案:
% F3 t5 j3 U* V8 V 1 .对于高峰期,将红绿灯时间分为四个阶段: 1. 编号为 1 3 5 7 9 11 (见图一)的红灯亮,其余的绿灯亮,持续时间 T1=65 秒; 2. 红灯灭,所有绿灯亮,持续时间 T2=27 秒; 3. 编号为 2 4 6 8 10 12 的红灯亮,其余绿灯亮,持续时间 T3=65 秒; 4. 绿灯全亮,持续时间 T4=27 秒。之后重复上述循环。红绿灯总周期为 T=184 秒。 % S7 n2 _7 T! V3 m# G
2 .对于次高峰的方案,红绿灯组合与开关顺序与高峰期完全相同只是各时段持续时间不同: T1=T3=35 秒, T2=T4=23 秒。 8 F X" L+ Z' B# \% H- p* h
3 .对于一般情况和稀疏情况,红绿灯顺序为:在所有路口,先红灯亮,持续时间为 T=30 秒;之后绿灯亮,持续时间为 T=50 秒。之后,重复循环。
# t: t- a( X9 m' t2 b% ] 由以上方案来模拟计算一天内环岛内车流量 Q ,其值超过最大容量的平均概率不超过 5.00% ,较为理想科学,所以此方案可行性较高。
; z8 r g. |- b1 U. M' w2 g6 l 最后,对模型进行了改进与评价。
3 l4 K2 A( M! j a1 L # ~) }# ^# u3 A6 W2 v' l
' N9 l- {. \/ T) r/ l+ P9 D3 T' k 关键词 :环岛车流控制 红绿灯控制 排列组合 随机模拟 等待时间 + Z; g# p1 v9 Q, [
; D5 C1 s2 e1 \% z
3 @1 N1 p1 O* m$ S4 }2 B( R
}( v: W' T& `8 O# o 8 F" a, S0 a% Z# ?( ]
2 }/ `( O! [, t' r& U
3 q1 J1 Y9 x- o# H9 V 巴黎凯旋门环岛有 12 个路口,其中有 2 条主道, 10 条支道。在进入该环岛的道路入口处可以设计有一些信号灯,或其他标志来控制车辆的流通。即为环岛制定车流控制模型,要综合考虑各时刻的车流量,环岛内最大车流量,天气情况,工作日与周末情况等因素。 ( U0 w0 f* m7 b F0 E7 t. z
我们的目标是,根据已知的信息,建立控制环岛车流量的具体模型,并分析该模型的优劣与稳定情况。
$ ?6 T+ |* |9 B, ?$ A 8 X1 S, O3 O0 m9 v+ S; w
! W: s! L1 f0 [8 V& Y4 V" D
. r5 \/ G" G* M% ?
* F+ M) ]& L: ?1 h
9 D: U+ }; H$ l2 A ) o% s. r* ^# \0 }9 }/ g! i
) W Z- N' _+ l' w6 L3 `& ^
1 .假设每个路口的进入车辆服从均匀分布(具体的分布情况见问题分析)。
, N, X; @+ r5 D: `+ k' M$ M8 I 2 .假设每个路口的离开车辆也服从均匀分布。
% d g" y5 L" b7 x8 C3 h9 a 3 .假设每个路口都配置有红绿灯装置。
! M+ W8 q5 z( O0 |' m 4 .假设环岛为单行道,只允许进入车辆沿着俯视逆时针方向行驶。
: E- o6 ^+ h7 Y1 p7 x! ^+ t* G* `# ~ 5 .假设进入环岛的车辆最多只在环内行驶一圈,不能多次在环内循环。 ; J7 Z3 p: B0 b5 c
6 .设环岛内与各路口处的车速为 20km/h ~ 30km/h 即 5.6m /s~8.3m/s 。 & a' g) L9 [& }/ B3 R. ]0 B
7 .假设只考虑正常情况下的交通,不考虑发生车祸和路面维修等意外情况。 ) D5 `! k/ H8 b: D# o7 X
! ]4 M% j3 E k$ ]
% v% H/ `5 U. g7 h+ J
# s; K- |% w2 E3 r: W
:环岛内半径。 " a0 g. u6 Y3 j. R5 _
:环岛外半径。
% _% t* B) M1 ~: {. W, ~2 ~* a :车底面积。
7 ^$ ]1 ]- A" E( n :环岛路面面积。其值应该为两圆面积之差, 。 & Y& e4 z8 }8 J* A f
:环岛可容最大车辆数。 取整。
" o3 i7 D8 p) N% l4 \! p Y :环岛内车辆总数。
0 ~& j/ `7 _6 E5 O1 A :环岛内车辆总数的当前值。
7 d/ Y7 n ^% C% R% }: H :各路口进入的车流量。( 1<i<12 ) ; R' ]3 _! q. R0 r4 _
:各路口离开的车流量。( 1<i<12 )
4 {4 H2 |( c% w" ` P :逻辑控制变量,用于表示各路口的通堵情况。 =1 表示通路,即绿灯亮; =0 表示堵车,即红灯亮。( 1<i<12 ) 0 R% x" K, `0 C- S
% d2 Q. j, C2 A* a) k: E :表示所有路口的流出车流量。
8 s+ e0 {8 y `6 d" E
" a0 `6 Z; L+ e. ^+ ]! w1 O :表示红绿灯持续时间,具体是红灯或绿灯,模型中会具体说明。
+ F" k' H. b; k( f* k & t8 A5 ?# @6 S4 d( C1 N9 g
:为某种情况下的堵车数量,具体模型中会说明。 7 c& W/ L2 e/ Z& m
6 r0 S2 ^ L8 b4 W } :车流密度,作为参考因素,将影响对车流量的模拟。 ; F7 p E7 h; t" D/ m. I3 O4 U
) w# _& c# L, v$ y
( h, s% j6 ?, u* W" D
- }, k5 Z5 i# [2 y) ?$ H& Z% V
此问题属于交通流问题,我们在初步考虑这个问题时,参考了交通流模型的结构方程 。我们认为影响环岛车流量的因素有很多:红绿灯,时刻(高峰期,平时等),天气情况,游客人数(虽然凯旋门游客时从地下进入凯旋门的,但是每个路口还是设置了人行道,所以游客的多少也会影响到车流速度,因而影响到车流量)。正常情况下我们不再考虑路面维修和车祸的影响。 2 h1 t# U6 i0 J8 v8 O
由上分析,我们需要做的是通过对交通流情况的模拟,找到最优化的红绿灯控制情况,从而达到车辆最优化控制的目的。 - Z: H& C B" X6 `0 ?. x5 Z7 p
因此,我们将此问题归类为最优规划类问题。 3 x2 c" B( A1 L4 [' R S
我们查找到了以下参数:
3 P+ _' `0 G- g! s9 h 凯旋门环岛每天平均车流量: 110 万 / 天。 % N ]% T3 @: n4 F/ d& h
环岛外半径: 80m 。
* v% \3 B% _% _) N1 c O 环岛内半径: 53m 。
: \& Y2 z- S. |; b: Z 一般中型车的底座面积:( 7~10 ) m2 0 ~ u. B" j1 }+ M5 C& b1 J
主道可以同时并行 3~4 辆车;支道可以同时并行 1~2 两车。 ; P. A) O8 M% G2 h( Z% k% G' \
: m6 S/ z" I2 ]2 z1 C1 H
3 {9 \ Z4 X8 N; i; l
4.1 环岛最大车容量: ) Q4 G9 w+ H8 x- R
由上面搜集的数据,我们可以计算出环岛最大车容量。 8 O% l9 Q) r) Y
环岛内半径为 ,外半径为 ,车底面积为
K' F, v) k6 \7 u( N+ m! l 则环岛路面面积应该为两圆面积之差:
p2 ~! f$ o; F7 r6 ^7 H, ? 。
, k" u4 ]1 v$ v 则环岛可容最大车辆数为: (取整)
/ x/ |8 {- Y! N* F 可得环岛最大可容车辆数目为: =1327 (辆)。 % W2 A( d0 Y: t' |
考虑到车之间应该有一定的间距,并且应保证环岛有一定的畅通,流畅性,我们设定环岛最大可容车辆数为 N=1000 辆 ( 稍微超过 1000 也行,我们只要保证严格地不超过 1200) 。 6 h; L6 }; J7 b6 Y4 a% y% J
- r* N' U+ M2 R! Y1 u5 K3 Q
; B" ^8 W* l+ E2 ]+ W( ^ 4.2 各时段的车流情况
& A8 {% y" _+ n3 e9 p& b
( w: w7 { i+ Y* T' R" _# k 工作日 2 I1 h' D4 \8 y
时间分布 # g7 `/ t% M3 Q# k
时期分布 / q& R2 R# m2 Q! b0 W) c7 |* A5 z
0 : 00~5 : 00
1 j' n5 I `3 M) Y 稀疏情况 , r/ \' G( o) i8 a6 r8 l2 X
5 : 00~6 : 00 - F0 c+ a' N8 _: K* N- v
一般情况
, [) _8 p$ c6 C 6 : 00~7 : 30
. G: z4 z8 z! ?/ ^ 次高峰 " Q e, Z9 Z. }+ W
7 : 30~9 : 00 " w j5 Y( Y: b8 Q4 X) a- e5 Z- q
高峰期
2 P7 v# P' g+ U 9 : 00~17 : 30 2 _; J8 B7 X' I. W( A5 w* R9 d
次高峰 . n# ~0 ?& u3 O0 \# ?
17 : 30~19 : 30
$ Q( U9 R8 j% G7 w4 n 高峰期
" I+ s5 G0 D; _# | \% k! r 19 : 30~21 : 00 ! ?9 [) b7 C$ x6 l7 {0 `5 `, g
次高峰 g6 N9 I# }! m
21 : 00~23 : 00 , H/ w) T5 a9 p# T5 X* E
一般情况
+ i& ^( z4 i' U: l+ Z 23 : 00~24 : 00
' W8 l' U7 `& F6 K% [9 I, a4 e" `; P 稀疏情况 2 a' }3 G3 P8 o1 K# V9 t% C
% k/ t0 n' a5 F# m4 ?; f
8 b' a. Z: l. J; t
周末 + G+ H E) i# |# S3 Z
时间分布
* f2 `/ g6 d: {" V! ?+ P4 C% j 时期分布 8 ^2 p0 m' X [- F) o# b, g
0 : 00~5 : 00 A7 c4 P# \; o s
稀疏情况
0 q' d Y5 o. \* d 5 : 00~6 : 00
# v4 o% H: I8 Q4 r6 u$ U 一般情况
+ _- F: d4 C& I6 e" T 6 : 00~8 : 00
9 _8 \/ @. N k6 a, l2 n" ] 次高峰 " [7 R, a% C1 i* u/ d: p6 R
8 : 00~17 : 30
& E: G7 b; m2 B; @. K, u( L 高峰期
9 k& r% x4 V% H6 p 17 : 30~23 : 00
6 I* C: _7 v7 w" e 次高峰 ( j$ O- k7 {+ o( o1 ~
23 : 00~0 : 00 0 a7 x# V, t4 }- |) p! x
一般情况
% v- N. p" L1 q$ N1 R* w7 O0 W& |
' ~* O, @7 }4 }6 S* u+ Q 说明:
$ W" H+ A; f3 h; M/ f4 Z5 n 在巴黎和法国其他主要城市,高峰时段的交通最为挤塞。法国每日的交通高峰时段是早上7时30分至9时及下午5时30分至7时30分的上下班时间。在星期五法国人一般都会外出旅游,所以交通高峰期会较平日来得更早,在下午4时起便开始阻塞,其中尤以离开巴黎的各条公路最为繁忙,而非高峰时段的交通一般非常顺畅。 5 B& b) k- u9 o1 _' p
8 w1 o1 I7 T \0 s3 f
4.3 对于交通模型的假设与估计
- ~6 w$ p( o/ O* B 对于交通流模型: ) e) c" {% C* ?4 R( Y( ]
其中: q 为车流量(即单位时间内通过的车辆数);
- [0 e! Q( B# c- [9 F) b" p2 W6 J9 { 0 A9 T3 U: N8 r8 S; w) l
9 e0 }. d, r# Y
为车流密度(单位路长的车辆数);
3 L8 W* g: H, s ! w% r6 q+ q& m3 e C. S5 e
+ P) P5 W; f+ k6 `7 A- C* ^ - L9 V; Z, k* B+ g7 A
为最大车流密度。
: g* W* L2 Y1 Z# ~/ {5 ~5 h
5 S5 V) x) o4 w8 e$ R ` - R; v& a2 N9 q+ E
0 t$ e! D. ?$ w/ P+ g1 q; b 为最大车速(注意:车速时车流密度的函数, )。 " c- g4 n9 w$ c3 T, r) c
根据上面的方程,我们可以估计出每个路口不同阶段的车流量,这包括流入与流出。 6 ^4 H9 @2 k& D2 l% d. A# [" [4 c
为了保证总塞车量最小,同时等车时间最短。我们针对不同时段对车流量进行了不同的划分:
% n' s/ Q p% z* s8 D 环岛内车辆总数 Q + @9 Y1 }9 _( o8 F# {
时 期
/ _. Q, V1 E/ t1 i1 f 有红灯亮 . _6 R+ ?. V# b8 A4 ~
无红灯亮
# j- X- t8 c6 ~# q! ` `9 E+ u, r 进 1 Z' {/ f7 S! n
出
7 K1 j5 s( U/ [, i 进 5 a' |1 A7 V0 X3 C8 b+ R2 Q( {' M- u: X
出 $ q* c0 u$ a& {' l6 s3 ?" S
主道 y
* a" U0 }5 g* T) i 支道 y
$ u2 y) d6 k" Q- `0 E% e" l" [( z5 w" g 主道 w
8 r# _' s1 u9 G$ Q8 u3 K 支道 w # X# T' C% i" b, v2 T
主道 y ) w, x$ W1 V3 Y' I) }( @
支道 y * l4 L& P! p% g( y1 x) u7 j, V. E' x
主道 w
5 _6 h# I8 g8 d3 C 支道 w
1 P- ~4 h* m0 H% ` 800~1000 & l5 b- a1 {& v$ H$ x( O) [
高峰期 1 W8 A9 g- ]7 [! i9 l
3~4 8 y: d' Y) Z9 I2 Y, Y2 S0 p* N
1~2 2 I# [3 R* F. P( z7 ^; K+ `% a
0~4 * \7 | I- \" j( w# I& x$ o
0~2
* ^8 ^! ^, x. S6 H 3~4 # U( U& \' T- x* a7 ?" l/ e
1~2
. |. ^ t) y9 ~ 0~4
9 _7 ~- J% ~) Q% Y+ L 0~2
! H# |9 p$ A* U2 w8 c 500~800
) J6 b+ c0 y0 W Y. } l 次高峰 : `6 \) E. k8 C9 l
2~4 / \8 u8 ?8 y2 K/ l! @) `2 W
0~2
; z! s$ ]5 `- ?* x2 ` 0~4 7 k/ }( A( ^1 i- d
0~2
( T i* H6 ]$ {8 I; h1 ` 2~4 9 A5 H. O' P2 Q" ^- T
0~2
6 o6 e& n8 C# q 0~2
1 k) H3 {( T8 g t( X; ~ 0~1 ' J4 S e! e \( X o; y( o" r
200~500 ) D0 F3 ?% d9 f# R/ @
一般情 况
# |; @! i6 k" G6 Q! D 1~2
* O6 K2 U: h4 K6 N( v 0~2
! X5 M1 Y2 J( M# | 0~4
5 S [& w$ H0 w8 c 0~2
; w4 F8 @* G. M1 n+ `% | 1~2
) [! c4 [* U0 R( ~' P- r M: M 0~2 ' J$ J M- a- j) v A8 n
0~2 ) Z4 [" o( t8 O; R0 T: i
0~1
* ~$ V+ q; R, I* X! w 0~200 1 I' o$ Z0 D# h- z- E
稀疏情 况 : j" c! t( U# H4 o
*
+ w ]) S4 Y8 ~% s5 V* K5 V * / j% V0 A# o) |$ G8 j; u
* 4 S4 y V2 H" y7 A% t' B' `
* $ g; B! _+ `6 [" u3 z
* : l: i9 k3 C8 z: T7 O0 W
*
, G' V+ J: D5 z; J Q *
: T: m' k7 i |& p2 u! k * + S! B) O3 S' `0 w0 U v& q
我们先设立一个逻辑控制变量 ,
' S1 w. e. k) d5 f( M0 v: G) m 对第 i 个 l 路口,当有车进入时, =1 (即认为绿灯亮)。 2 l8 ?0 N$ R3 s$ H# ]# v6 @
当没有车进入时, =0 (即认为绿灯灭)。 4 ^/ J8 S# D; W8 d. E- H9 ]
又设 为第 i 个路口的车流量(辆 / 秒)。
# z# Y; Y7 W( G, S) W* z 则我们可以列出下列等式: 1 h' o* V% A4 p2 g$ Q1 O
根据:单位时间内,环路车流量的增量 = 流入的车流量—流出的车流量。 ; x; S* c. n: R* j. E
0 d d4 b B+ Q
/ \5 q% u# x+ f% u) _! J$ n: K % y7 {! \- X3 y, c. _: Z
dq 表示单位时间内环路车流量的增量。 + S! H5 y$ v! w
0 |+ g1 @; t$ d 对于 以及 我们可以用 rand 模拟。
! A4 n; M8 M5 O$ t2 Y) s J( _0 y8 Q9 `/ A6 ?# { Y
因此,环岛内车辆总数 Q 满足: ! S+ P5 h) [7 o& o
4 y1 e$ ]* p1 x7 Q3 S
注:
8 F; L. K: e+ C" l7 j 由于 的组合有很多种,我们加入限定条件,即要保证等待时间最短(红灯亮的时间最短),以及等车数量最少。
$ j, L) a+ O/ L& a* l
! g5 ~5 Q) ` R. g: g$ @* B& T) h- a 因为等待时间就是红灯亮的时间 T ,等车数量又与车速和等车时间有关。
+ \ R6 P9 f. H
: c* o1 c* W5 S9 Y$ n 为此,我们设立下列函数: ( s5 x3 ]+ D" y" i" a
0 d+ V. f/ S5 ~/ F5 K
7 `2 y& ?( g( N# A! R
: m J& W' o; X7 Z
3 g& z V6 T3 `/ Q; r 说明: 2 Y( U* ]/ c$ s6 B* {; k
为各路口的逻辑值(通为 1 ,不通为 0 ) 0 \ Z/ X1 V8 ~5 a
0 T- S g' H1 K. p 为第 i 个路口的车流量(辆 / 秒) ! P/ m/ }) w' {4 S7 {. Q
为循环中第一次亮红灯时的堵车量, 为第二次亮红灯时( 的对立面)的堵车辆。 . S' Q$ |/ `( P& f; h2 U. p
为总堵车辆。
3 L3 I. Y5 Q- \9 t0 x
0 t8 H/ [3 Y+ R: r" K. b 上面的分析可能需用到下列参数值: ' X& h, [6 @9 B. R$ R
1.8 B% i3 ^5 c/ h" U4 E7 n# C7 s' F
每条路段上的最大车流量。 6 \. p4 I% B; I8 d
2.
' `+ { F; W. F2 B' b8 S+ s6 _3 Q- ? 每天路段上的最大车流密度。 - z, C. Q+ P, O
3.
C; u9 \3 y7 I" a/ E 每条路口进入的车流量(辆 / 秒)。 ! ^8 t0 K# n: Z; a4 W
4.; [; L& G$ V0 I- y3 K
每条路口开出的车流量(辆 / 秒)。 ( C2 P0 A& ]4 z9 g5 @
! d6 ?( U' F+ I, e1 z 通过模拟,我们将在不同的组合中找出最佳的红绿灯方案,并通过多次模拟,确定时间分配。 3 K/ F0 m$ J4 ^: h4 H+ ^4 K
x, u4 P$ I( Y& T6 W7 G
" g0 [/ ?# c3 ]7 A `: I) L 一、对于高峰期时,我们对于下列组合进行了模拟(程序见附件):
# Z/ x9 N4 Z& F0 ? 红灯亮的个数(盏)
# R: p4 M( t2 c9 S0 z/ U5 L9 Z- R 12 & R: S h1 }( |' P3 g
11
; O* h! z$ u1 ?2 X' b 10 , Y0 M( j' J6 S9 V0 ~0 q6 k
9
- D; U) { E( N; p+ c2 g" q' N 8
& r k' W: O1 ^5 M _ 7 ) s% m: B# e+ F
6
A+ L$ n5 q) l4 V6 W2 f 5 - L" H% V& x& v
4
; @5 p! P4 a5 Z7 q/ I4 [+ _ 3 4 {9 N \5 F1 Q: O& _ V
2
0 U* t7 m3 m, A6 M- Q8 u( ~8 {) f 1 . v. g8 \0 A8 z. a' O3 `* m% Y& a
0 & J5 T# }( f0 F: x- K
平均最短等待时间(秒)
! o) ?! P/ x+ t 16 * U$ [3 Y4 U% K' N$ v7 q
20
: K3 S7 {, Y; u! r) M/ E: y 22
9 Z5 |4 g' Q* ]. G3 n8 a# z 27 6 r; b. V- C, p3 g9 G* a( g; ?& [
30 . s7 X! b# }$ E$ D% ~' O
42 1 C4 t: [# \ i
70
0 d" P, I& I! g5 Z. W) | 154 - d* ^4 s& a+ D) F
Inf 7 }/ w6 x! Z! n& f# T
( 无穷大) . p0 C/ _ V, U
Inf % }$ e" ~0 q/ `; A9 C
Inf 0 v5 i% Q7 j x" [1 a* _
Inf " y l% b$ Q( I) L5 c$ d7 `4 c
Inf & \. x4 f% [9 B E) n- P
注释:
h" P% L/ T, S& |& [ 对于红灯亮的盏数,我们可以有很多种组合方式,比如红灯亮 1 盏,可以是 1~12 编号中任何一个亮,但这些组合中总是有一种或多种为最优组合,这从我们程序结果可以看出。
1 E- w, A& d, S u3 K# `- o: b 0 H/ U7 @. @9 @# e. C, ~
分析: ^4 V b" l8 J# L6 z
# a' q# n. ]' h. V 有 0 盏红灯亮: ) r3 F+ V- g4 i$ ~9 n$ m/ h
此情况显然不合理,因为没有红灯就无法控制环岛总量。
& @7 L( P5 V+ C" _7 n1 E7 l' c7 U
4 Y0 q6 z. q2 l 有 1 盏红灯亮: $ T4 L+ g+ O) n: K
对于此种情况,经过模拟发现不可能达到降低环岛车辆的目的,反而,环岛内将更加拥堵。(过程见程序) ' h' ]+ R& y- M, \% t
, v- t+ V t, C' p" S* G# i 有 2 盏红灯亮:
3 c* v* q2 }2 z 此情况结果同上,不可能达到降低车辆的目的。 0 r! V; y0 Z2 m8 m6 V" A# f# l
+ {% B8 i k! X" g/ L/ f 有 3 盏红灯亮: : I' Q/ M- |# c/ L1 c, s
此情况结果同上,不可能达到降低车辆的目的。 ) g3 H. S; j6 m3 m
; |' M& R& P I8 X5 l; Z# I4 {% ?
有 4 盏红灯亮: & |9 \6 j& \! S
此情况结果同上,不可能达到降低车辆的目的。 # k+ q l, \' P$ H' v2 b- \$ ?
- m( s# i" }4 {9 f& N 由上分析说明,红灯至少应该亮五盏以上。 ' y& D1 W* X) z- P
9 r; p6 t: F M, y0 D 为此,我们排出一下组合: : @# @6 n" P1 B4 ]8 S5 @* H
5 —— 7 : $ i7 k# f5 j4 X. ?- P4 c
此种组合方式下,可以分为:
! a$ p, Z. {" l" q1 O |" T a. 开五盏红灯时(包括一个主道)的等待时间为 144 秒,开七盏红灯时(包括另一个主道及其他支道)的等待时间为 48 。
8 j& Y& n9 O. p* c; T) g, j s 此时,总塞车量为:
, G- F+ v8 @) r0 y& \4 N7 s/ D 3 G3 S) E0 d: ^6 h5 x7 x
b. 开五盏红灯时(不包括主道)的等待时间为 inf 秒,故此情况不成立。 # b) b3 O- r0 |' Y" _
0 q; D% d0 k- e6 ]( J; m
6 —— 6 :
; U- r1 ?6 F+ m7 p0 A( W; U 此种组合方式下,开六盏红灯时(包括一个主道)的等待时间为 68 秒,开另外六盏红灯时(包括另一个主道及其他支道)的等待时间也为 67 (由于分布相同,其实两个等待时间应该一样,但由于是模拟,不可避免的造成一定的差异)。 3 }/ {6 `5 F$ g' @
此时,总塞车量为:
2 w3 ~+ Q1 {. V! m& u: Z . k3 R" w( C% ]% e9 P `1 N
8 a9 M7 M; {. M/ U 在保持堵车时间尽量短,堵车数量尽量少的原则下:
d3 G; W, V+ t) c* z 只有选择 6 —— 6 组合是最优的。 " a9 E. {, A& o1 b
根据等概率原理,各条支道应看做概率上相同的路口,而两条大道也是等效的,因此,在模拟时,我们就可以人为地设定组合,只需保证总数按红灯亮的盏数分布即可。比如:对有 6 盏红灯亮,我们选定组合编号为: 1 3 5 7 9 11 ( 1 为主道),此组合方式与 2 4 6 8 10 12 等效。
$ H& r/ E8 t, u% S3 o9 r- t $ U$ H( b) O; @! n! F% M4 t3 N
这时我们可以确定红绿灯的循环模式。
3 M; {* k2 H* @7 U 不妨设定,先使编号为 1 3 5 7 9 11 的红灯亮,在经过 T(T=68 秒 ) 后,打开所有绿灯(包括原来的绿灯),再等环岛内车辆上升至限定值后(经计算 t=25 ),再打开另外路口 6 盏红灯,其编号为 2 4 6 8 10 12 。之后,重复上述循环。
' T+ g8 q5 P+ V" Z1 W
6 Q8 ?2 {7 ?: ^2 e h @; g0 ?5 S* z" L/ z
二、对次高峰,模拟结果如下(程序见附件):
S5 J+ h8 {* Z8 U 亮红灯个数(盏) - X" E- v0 T d7 w) @# K
12 ' h5 U3 J0 I4 s2 a- A/ L
11 # `* r4 V; m( H( C. b* i3 Q6 V
10 P3 P5 e L" B, B# d1 L5 k- y7 @
9
U: V) N+ l# L5 r( h 8 7 @! v. O- l3 A1 b2 M8 l6 b
7 ! ?0 P* G8 \% x8 D# @
6
3 t) J7 V& d1 H4 L1 X8 @ 5 4 H4 ]; A" r! x3 W! |1 ~7 g
4 6 M m( T7 g) l+ n- o
3
9 [) S, Z# Y/ H6 l 2
0 ]: E$ J5 N( C ^! I+ I4 ~ 1 7 ^/ c) U; L6 l. f: I' Q# K
0 ; W) W! [7 C/ v7 G: h1 u+ y& X
平均等待时间(秒)
) v; S, K2 |; k& a 24 6 ^4 @2 L1 X: r% S" G# \
30
% r/ p* h) z! {* j3 |4 T8 Q 31 , Y8 Q: Q( z+ T. Y& o7 z A: |
32 : l; X5 m7 T c: {7 [' a7 O& ^! ]( B
35
0 ^' d B5 g; e1 C1 f 43
/ M: T6 f' H* H- y 57 # y: x/ v( [) q' e
68 8 A4 h' P/ @: R& D2 J& v9 D# j
96
8 j. E( I4 T' O/ {6 q Inf
: q7 n e% d- l0 q5 M1 N ; z0 r/ [0 ~' k) y, a) e' V
Inf ( _* }' W5 b+ R
Inf 6 B: H( B8 C4 ]: h7 j5 C" U
Inf
9 L7 k6 |" o% E: d
4 w# Z$ Q; I* ^3 R
说明: $ y% G x1 C5 j5 W) ^
对于红灯数目小于 4 的情况,实际上有的模拟值满足要求,但由于等待时间太长( 100 秒),并且情况及其不稳定,多次出现 inf ,也就是不能达到降低车辆的效果,我们认为这些情况都不现实,均统一成 inf 类。
7 V+ I( G3 N7 I+ C# L+ K; g: U Z" e; f C6 J% Y; G4 O
由上分析说明:红灯至少应该亮四盏以上。 ; u7 ~ N' M# y) L0 r& X
$ V" E3 K% W: Y1 g3 N( p5 b 为此,我们排出下列组合:
4 @: f# g* P( L% m5 ~6 f- L, c 4 —— 4 —— 4 :
, x/ v8 I' z$ H 此种组合方式下,开 4 盏红灯时(编号为 1 3 4 5 )的等待时间为 80 秒,开另外 4 盏红灯(编号为 2 6 7 8 )的等待时间为 77 秒,开最后剩下的 4 盏红灯( 9 10 11 12 )的等待时间为 147 秒。 & X- D4 w. N/ G
此时,总塞车量为: 1 w1 _- V) _- E1 q
6 i5 }8 p+ D, n6 C 4 —— 8 : ! p9 i5 b( Y; x% B' S4 r/ J8 f% e
此种组合方式下,开 4 盏红灯时(编号为 1 3 4 5 )的等待时间为 80 秒,开另外 8 盏红灯(编号为 2 6 7 8 9 10 11 12 )的等待时间为 39 秒。 % A' w$ h, p9 o% B
此时,总塞车量为: + n5 {* a1 \6 s; h8 j
p& `7 P( z; A. @! H% m( [8 g % J+ ]5 n. z4 M
5 —— 7 : / g: _& I- w' k0 w0 X. R
开五盏红灯时(包括一个主道)的等待时间为 58 秒,开七盏红灯时(包括另一个主道及其他支道)的等待时间为 45 。 + ~ L2 u3 e- F1 A6 k7 r; C0 s
此时,总塞车量为: 1 \- M* U2 k$ |0 Y4 D
. y; Q% E: |/ p: X
0 |4 S4 |0 f& | C. I% \8 B 6 —— 6 :
# I3 Y9 p; T( X, o0 d- b2 E6 Y. u 此种组合方式下,开六盏红灯时(包括一个主道)的等待时间为 50 秒,开另外六盏红灯时(包括另一个主道及其他支道)的等待时间也为 50 (由于分布相同,其实两个等待时间应该一样,但由于时模拟,不可避免的造成一定的差异)。
; x1 q" \+ U' b( |0 r# I1 z 此时,总塞车量为:
: q6 {# U7 V( H5 k : }1 [: y" ?8 R* L; @) J6 f9 t
5 g; i1 ?. M, Q1 ^2 y. \9 E 由上可知:
: ~* m) a; r/ y% P& T 对于高峰期和次高峰期都应该选取 6 —— 6 的组合,并且将两条主道分配到不同的组合中。 ( S0 i& P) L4 J# n
& R* e. \& Z8 @# v9 F: B" Q
5 t8 f9 D. I; G 说明:(为什么选取组合时两条大道不能同时选取?) 6 }; T# V$ H) P* D' ]0 ]" ]$ D
下面只针对高峰期说明: 3 \& Y% D; p. M5 x' {! I
对于高峰期 同时选取两条大道的情况:
! P* K- d7 v2 t. X 有 2 盏红灯亮:如果同时选取两条大道通路时,经模拟,其等待时间 T=inf ,也就是说不可能达到降低车辆的目的。
' j* Z! f3 I) ~/ w! q8 C$ C 有 3 盏红灯亮:如果同时选取两条大道通路时,经模拟,其等待时间 T=inf ,也就是说不可能达到降低车辆的目的。 : ^8 k0 ^0 d( I1 a& G
有 4 盏红灯亮:如果同时选取两条大道通路时,经模拟,其等待时间 T=158 秒,并且多次模拟发现其等待时间出现为 inf 的几率很大,也就是说不可能达到降低车辆的目的。
- k8 W% \3 d1 p( m) n; w 有 5 盏红灯亮:如果同时选取两条大道通路时,经模拟,其等待时间 T=101 秒,并且多次模拟发现其等待时间出现为 inf 的几率很大,也就是说不可能达到降低车辆的目的。 , G0 u2 ^* z2 C$ B
有 6 盏红灯亮:如果同时选取两条大道通路时,经模拟,其等待时间 T=63 秒。
5 ^& \# @# i+ u+ @8 | 有 7 盏红灯亮:如果同时选取两条大道通路时,经模拟,其等待时间 T=50 秒。
9 ?/ e1 ~# F( X 有 8 盏红灯亮:如果同时选取两条大道通路时,经模拟,其等待时间 T=45 秒。 6 F! ^3 c* f+ ^0 ~! M
有 9 盏红灯亮:如果同时选取两条大道通路时,经模拟,其等待时间 T=35 秒。
8 _( b' S0 x: W; G 有 10 盏红灯亮:如果同时选取两条大道通路时,经模拟,其等待时间 T=30 秒。 ! ?+ B6 `4 ?: U% P5 u
5 @1 G/ A6 t# N+ g5 B
同样地,考虑到我们设计的算法,对于高峰期 ,不可能不选取某一条大道,所以我们只需考虑对称选取,即组合时尽可能的将大道分配在不同组合中。
. \$ M7 R1 _9 g- O" W( h 有 2 盏红灯亮:只选取一个大道通路,另一个大道堵塞时,经模拟,其等待时间多次出现 T=inf ,说明此种情况不可能大道降低车辆的目的。
4 L4 S. C4 ^, `' c 有 3 盏红灯亮:只选取一个大道通路,另一个大道堵塞时,经模拟,其等待时间 T=96 秒,但也多次出现 inf 的情况,因此不考虑此种情形。
$ ]! A9 w! Y2 K7 Z$ @ 有 4 盏红灯亮:只选取一个大道通路,另一个大道堵塞时,经模拟,其等待时间 3 l9 g0 G6 ]2 p6 ~9 x
T=65 秒,但也有很大的几率出现 inf 的现象,也不考虑。
2 I1 Y9 Y# A- R0 Y7 l5 d4 f 有 5 盏红灯亮:只选取一个大道通路,另一个大道堵塞时,经模拟,其等待时间
1 w7 X& n/ c( ?: d T=45 秒。 / q9 z- K7 V M+ u
有 6 盏红灯亮:只选取一个大道通路,另一个大道堵塞时,经模拟,其等待时间
U- I) Z# Y% _/ z( }% s T=35 秒。 - z. [6 T4 Y4 H8 n/ j
有 7 盏红灯亮:只选取一个大道通路,另一个大道堵塞时,经模拟,其等待时间 8 N9 z! }& S3 H2 u
T=31 秒。
9 Y. N5 A; c6 Q. P* n4 y 有 8 盏红灯亮:只选取一个大道通路,另一个大道堵塞时,经模拟,其等待时间 8 W" q$ R8 Z D% }3 C0 e
T=27 秒。 6 B j P+ C9 ?- d
有 9 盏红灯亮:只选取一个大道通路,另一个大道堵塞时,经模拟,其等待时间
* U, ~0 ]% a& g4 L5 f/ X( n$ u T=25 秒。 8 G7 l: N* r2 m! V! M: E# K" g
有 10 盏红灯亮:只选取一个大道通路,另一个大道堵塞时,经模拟,其等待时间 =23 秒。
5 Y/ ]* @: \+ U0 v( N * i) ?2 l1 A) N9 [6 v
对比上面的两种组合下的结果,显然第二种情况更为节约时间,对于所有红灯亮的情况,只选取一个大道通畅的情况能保证等车时间。因此,我们认为,选取组合时两条大道不能同时选取。
5 y. L3 q; ~( y/ A% ^ 0 Y5 l2 ~9 g% ^& X) p
由此,我们可以得出高峰期和次高峰期的红绿灯控制方案: 2 N8 r( {! o; `$ q
对于高峰期的方案: 1 N4 A* l- N m+ w
先亮 6 盏红灯(包括主道的那个红灯),其持续时间为 T=65 秒,此后,红灯全灭,绿灯全部打开(包括原来的),持续时间为 T=27 秒(在高峰期,对绿灯全部打开的情况,即所有路口通畅时,经模拟,通畅时间为 T=27 秒,此后若不打开红灯限制车流入,将超过环岛最大车容量,因此,时间不能超过 27 秒,但是,我们为方便设计考虑,将时间定为 27 秒)。之后,又打开另外一组没亮过的红灯,持续时间为 T=65 秒。结束后重复上述过程。 ) _/ i/ S& a6 O9 z- j5 j( `
" S/ {) D. z4 v! u( E' a
对于次高峰的方案:
% h- h s( p" R: M9 U5 ? 先亮 6 盏红灯(包括主道的那个红灯),其持续时间为 T=35 秒,此后,红灯全灭,绿灯全部打开(包括原来的),持续时间为 T=23 秒(在高峰期,对绿灯全部打开的情况,即所有路口通畅时,经模拟,通畅时间为 T=25 秒,此后若不打开红灯限制车流入,将错过环岛最大车容量,因此,时间不能超过 25 秒,为此,我们为方便设计考虑,将时间定为 20 秒)。之后,又打开另外一组没亮过的红灯,持续时间为 T=35 秒。结束后重复上述过程。
1 w3 c# p# y3 h! p
" S1 |4 J- }$ f' z 1 k) K$ m o9 T5 J4 k
三、对于一般情况与稀疏情况的说明: + S( e* U' j. L7 }+ m
# B8 v& J9 V% P6 r/ } A.
6 q: m/ I$ |; f# j$ u# ~) f% \. U/ v5 c 一般情况:
7 [3 v- P$ y9 P" C/ ~- y- g 对于 6 盏灯的组合(每个组合只分配有一个大道),其等待时间 T>150 秒,如果红灯时间仍然按此时间设计的话,肯定是不科学的,因为不可能让汽车等待如此之久。因此,我们从尽可能减少等待时间为标准,经过模拟,发现当所有路口均亮红灯时,其等待时间最少,为 T=23 秒,而这也符合一般城市中非高峰期的等车时间,我们为了方便设计,将此时间定位 30 秒,而 30 秒也是一个可承受的等候时间。对于绿灯全开时的情况,更趋前面的假设,经过模拟,畅通时间为 T=50 秒。因此,我们选定一般情况时,红绿灯亮灭的原则时,所有路口红灯全亮,持续时间为 T=30 秒,此后红灯灭,绿灯开,持续时间为 50 秒。之后,重复循环。
1 A3 F- @ S+ F) c% M B. 稀疏情况: ( D6 W1 o* Q, J) } Y
对于稀疏情况,车流量具有不确定性,我们无法估计具体的车流量,但由于此种情况下车流量很小,我们可以将之归到一般情况,并且以一般情况的红绿灯规则来控制。 9 @7 ^9 l1 V2 D& o) A
( p0 p) _6 W0 T( u9 t- J9 r
( Y6 Z- H( D3 c+ Z- P' |* P2 N
3 j4 ^5 M0 g0 j$ }$ X' M, j" R 根据我们的方案,我们采取随机模拟的方法,分别对高峰期,次高峰,一般情况和稀疏情况进行随机模拟。 * P% f! w4 g! L7 T) i- b* N+ L
为了保证环岛内交通的流畅,我们设定环岛内的车辆总数 Q 不能超过 1000 ,但实际上换岛内最大车容量为 1327 ,因此,我们在考虑交通流畅性的前提下,可以适当地放宽这个限制,严格规定 Q 不能超过 1200 。
4 u% ]3 _& M2 E9 B$ i 我们检验的目的是为了了解模型的稳定性,为此,我们对四种情况分别进行了 24 小时的模拟,其结果如下:
% A# h1 y% p3 t2 j$ D5 ?" ~7 u 1. 高峰期:(程序见附录) ! Q4 c( ~: R% ^
第一阶段红灯持续时间 t=65 秒
; g# r5 g `# L7 q2 b+ q# \ 第二阶段绿灯持续时间 t=27 秒
2 n) p7 z! `- c: F 第三阶段红灯持续时间 t=65 秒 ( X& e1 e4 Q" F" Y! L9 a3 H; \
第四阶段绿灯持续时间 t=27 秒
: i8 r' Q0 V2 \1 E 总周期 T=184 秒 1 E0 A! x9 W! b) n( H9 c" Q6 e9 F
* K8 q1 i) @* c- M$ @) N
对于此方案,我们在模拟时发现,由于每周期都会累积一定的车辆,也就是误差,在很长时间后,其累积的误差将达到非常大并且不合理(超出最大容量)的数值。因此,我们需要增加一个修正时间,并且此时间应该很小,只在车辆超过一定数量时才加入。 5 |( D. a+ Y; s+ I0 m% A
我们的做法是,当环岛内车辆大于 1000 时就对红灯持续时间加 3 秒钟,即此时红灯持续时间 t=65+3=68 秒。在车总量 Q 没有超过 1000 时,我们仍然以 65 秒的规定时间运行红灯。
9 H% s2 Z1 S# l: L. E. y- o 这样,我们模拟 24 小时高峰期后:超过 1200 的车辆次数为 37 ,占一天内车辆总数的比例为 1.97% 。(这只是模拟一次的情况,在模型改进中,我们模拟八次后取平均,得出更加准确的比例: 2.74% )
2 T% T% ^7 E9 W. @( X9 B 对于此比例,我们认为是相当小的,也就是说,发生环岛堵车的概率时非常小的,因为我们是对 1 天进行模拟,累积误差显然会相当大。而一般的高峰期只持续 2 小时左右,累积误差必然很小,其堵车概率也应该低于 1.97% 。 6 G L3 C# f$ _1 q8 U# U4 g1 p
: S5 }* m8 @. y5 I7 V) I
2. 次高峰期:(程序见附录)
' `4 }8 Z3 f5 S" l0 b: z 第一阶段红灯持续时间 t=35 秒
0 h& H9 K; {* _. p7 @ 第二阶段绿灯持续时间 t=23 秒
7 L$ s- r. ]) f' O9 D 第三阶段红灯持续时间 t=35 秒 - {+ X9 A" a; y v5 `: v8 { L" J
第四阶段绿灯持续时间 t=23 秒 # V+ N7 ^! M( i& H
总周期 T=116 秒
4 F! t0 Y8 S, q$ B4 y 对于此方案,我们为了保证环岛被最大利用,同时又能使交通运转顺畅,设定环岛内最大车辆数不超过 800 ,经我们模拟 24 小时次高峰:超过 800 辆的几率为:
" m$ ~( R0 Z; K# U- T4 R ,显然这是非常好的方案,鉴于此,我们不对此方案做修正,即沿用模型建立中确定的红绿灯持续时间。 ( F5 D' w, ~- L' ?
3. 一般情况和稀疏情况: & i ^9 @8 K4 K( G: G
因为车流量的原因,不可能造成交通的拥堵,因此,我们不在对此情况做模型检验。为了说明时间安排的科学性,可参考其他大城市的一般情况的红绿灯时间。
+ J3 h: ?8 P$ Q; B
" }: e6 n. p' a, \# C( { * S8 M5 Z. |/ J/ i' a' m
/ N5 A' q3 e0 w 1 .对于工作日和非工作日,由于车流量的分布不同,我们可以根据表 1 来设计红绿灯时间安排。
9 x, v% B) Q- a7 y 2 .我们只考虑了每个路口流入与流出的关系,并没有考虑到车辆在环岛内的绕行情况。所以可以增加限制条件:环岛内并行车辆不碰撞,这样可以选出更加优化的方案。
* {" ?7 H6 }' ~3 H$ U5 F 3 .不妨考虑车辆在环路中的相位问题,这项可以细化到每辆车的行驶情况,但这样相对来说较为复杂,我们不予考虑。 ; v9 M- `' Z0 o4 |. u2 m5 B( G
4 .对高峰期时间的修正:
; y* o% y& U% J, a' Q/ c k 若不对高峰期的红灯持续时间作修正,则经长时间后,累积误差将使环岛内车总量超过 1200 (我们称之为危险),这是非常可怕和不安全的。为此,我们对红灯持续时间做一点微小的修正。经过我们的模拟:(程序见附录) , h+ `' l- ?0 |/ J
修正时间 t=0 时,出现危险的几率: 89.62% 。
* o# h B; u- p& I 修正时间 t= -1 时,出现危险的几率: 88.56% 。
3 m& ~! T4 r- @3 O+ Z$ \/ V) h% X 修正时间 t= -2 时,出现危险的几率: 98.03% 。 ; V! V& s$ y+ s1 K- H2 f+ b6 i' c
其实,如果减少红灯时间,显然,这时在这段时间内进入的车辆数目就会增加,在不修正时已经危险的情况下当然就会照成危险几率变大。
3 G2 `8 m3 T8 k5 d+ b5 ]: S* v 所以,我们应该将修正时间调为正值。 $ F) r% r9 `4 S( d& |8 v- a. [
修正时间 t=1 时,出现危险的几率: 93.33% 。 2 n7 ~+ x7 N$ L6 p
修正时间 t=2 时,出现危险的几率: 13.74 % 。
1 F% w+ k% e0 C8 O' T; o: B" c# J# L( D 修正时间 t=3 时,出现危险的几率: 2.74% 。
) r; h7 I; Q/ A0 x0 C5 o 因此,我们以 5% 为限定,确定出修正时间为 3 秒。
" k& {, _& ] I% Z3 {' J4 |
, [4 }$ W) R, V$ X: l, o9 b 8.1 优点
3 S, j T; q! ?- D. b. O$ p
W0 C" c0 u* H0 t" @3 @ 1.本文对不同车流时段(高峰期、次高峰、一般情况与稀疏情况)模型分别进行了模拟计算,得到了最优组合下的红绿灯循环时间。由于车流量是基于模拟的,并且环岛内车辆总数也是先设定的,因此,我们的模型可以适用于很多情况。并且,根据我们的模型,对于已知车总量和具体车流分布情况,可以重新确定出最优化的红绿灯控制模型。 0 a" u1 @1 n* n4 p+ Y
1 x/ o: f, Y2 W8 y* v1 a 2.在建模过程中,我们对所有可能出现的红绿灯组合情况进行了模拟,这样最终得到的最优组合的方法是很科学的。
$ Z$ {+ E P3 U7 \7 b
$ d Y- P8 V; [1 d) |5 D6 }- B 8.2 缺点
& Q \3 e; B2 J8 H& W 7 F* x! g) W# C
1.在模拟模型的过程中,我们假定车流量服从均匀分布,这带有一定的主观性,并且我们并没有考虑每一辆车的具体行驶情况,比如车辆在环路中的相位问题,这可能造成某些紧急事件发生时不能及时疏通道路的问题。
zan
总评分: 体力 + 1
查看全部评分