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