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