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