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