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