Живоглядов В.Г. Теория движения транспортных и пешеходных потоков
Подождите немного. Документ загружается.
2.7. Экспериментальное исследование
Следует заметить, что в данном разделе осуществляется проверка тео-
ретических разработок на достоверность и работоспособность в аспекте оп-
тимальности посредством аппроксимированного моделирования.
2.7.1. Блок – схема алгоритма расчета основных характеристик
светофорного объекта.
Начало
Ввод
n, tr, ν
i:=1÷4
Ввод
λ
i
jk
, τ
i
jk
j
:=1÷4
k
:=1÷n
i:=1÷4
Ввод
T
i
jk
j
:=1÷2
k
:=1÷4
Ввод
S
j
,
μ
пj
, μ
лj
, k
пj
, k
лj
j
:=1÷4
(1)
(2)
(3)
(4)
(5)
j
:=1÷4
k
:=
1÷
4
S
1
:=0; m:=0
i
:=1÷4
S
1
:=S
1
+λ
i
jk
m:=m+τ
i
jk
λ
jk
:=S
1
;
τ
jk
:=m/4
;
(6)
231
j:=1÷4
k:=1÷n
t0l
jk
:=0;
rl
jk
:=0;
(7)
j:=1÷4
k:=1÷n
(12)
i:=1÷4
j:=1÷4
k:=1÷n
h
i
jk
:=1
-
e
-
λi
j
k
;
t0l
jk
:=t0l
jk
+h
i
jk
*T
i
jk
;
rl
jk
:=rl
jk
+h
i
jk
;
rl
jk
<>0
да нет
t0l
jk
:=t0l
jk
/rl
jk
t0l
jk
:=0
c1
j
:=max2(t0l,j,n);
c1[1]>c1[3]
да нет
T01:=c1[1] T01:=c1[3];
c1[2]>c1[4]
да нет
T02:=c1[2] T02:=c1[4];
T0:=T01+T02;
(8)
(9)
(10)
(11)
232
n
(13)
(12)
c_1:=0;
c_2:=0;
c_3:=0; c_4:=0;
i:=1÷4
c_1:=c_1+
λ
i
11
τ
i
11
c_2:=c_2+ λ
i
31
τ
i
31
c_3:=c_3+ λ
i
21
τ
i
21
c_4:=c_4+ λ
i
41
τ
i
41
c
_1:=
c
_1+(λ
3
31
τ
3
31
+
λ
4
31
τ
4
31
);
c_2:=c_2+(λ
3
11
τ
3
11
+ λ
4
11
τ
4
11
);
c_3:=c_3+(λ
3
41
τ
3
41
+ λ
4
41
τ
4
41
);
c_4:=c_4+(λ
3
21
τ
3
21
+ λ
4
21
τ
4
21
);
да нет
c
_1:=
c
_1+((μ
п1
+ μ
л3
)*ν)/(
k
п1
+
k
л3
);
c_2:=c_2+((μ
п3
+ μ
л1
)*ν)/( k
п3
+k
л1
);
c_3:=c_3+((μ
п2
+ μ
л4
)*ν)/( k
п2
+k
л4
);
c_4:=c_4+((μ
п4
+ μ
л2
)*ν)/( k
п4
+k
л2
)
ошибка
(
k
п1
+
k
л3
)<>0
или
(k
п3
+k
л1
)<>0 или (k
п2
+k
л4
)<>0
или (k
п4
+k
л2
)<>0)
c_1>c_2
да нет
q:=c_1 q:=c_2
c_3>c_4
да нет
r:=c_3 r:=c_4
( )
(19)
n=1:
конец
233
c1[1]:=
λ
1
11
τ
1
11
+
λ
2
11
τ
2
11
+
λ
3
32
τ
3
32
+
λ
4
32
τ
4
32
+(
μ
п1
+ μ
л3
)*ν
/(
k
п1
+
k
л3
);
c1[2]:= λ
1
12
τ
1
12
+ λ
3
12
τ
3
12
+ λ
3
32
τ
3
32
+ λ
4
32
τ
4
32
+( μ
п1
+ μ
л3
)*ν /( k
п1
+k
л3
);
c1[3]:= λ
1
31
τ
1
31
+ λ
2
31
τ
2
31
+ λ
3
12
τ
3
12
+ λ
4
12
τ
4
12
+( μ
п3
+ μ
л1
)*ν /( k
п3
+k
л1
);
c1[4]:= λ
1
32
τ
1
32
+ λ
3
32
τ
3
32
+ λ
3
12
τ
3
12
+ λ
4
12
τ
4
12
+( μ
п3
+ μ
л1
)*ν /( k
п3
+k
л1
);
c2[1]:= λ
1
21
τ
1
21
+ λ
2
21
τ
2
21
+ λ
3
42
τ
3
42
+ λ
4
42
τ
4
42
+( μ
п2
+ μ
л4
)*ν /( k
п2
+k
л4
);
c2[2]:= λ
1
22
τ
1
22
+ λ
3
22
τ
3
22
+ λ
3
42
τ
3
42
+ λ
4
42
τ
4
42
+( μ
п2
+ μ
л4
)*ν /( k
п2
+k
л4
);
c2[3]:= λ
1
41
τ
1
41
+ λ
2
41
τ
2
41
+ λ
3
22
τ
3
22
+ λ
4
22
τ
4
22
+( μ
п4
+ μ
л2
)*ν /( k
п4
+k
л2
);
c2[4]:= λ
1
42
τ
1
42
+ λ
3
42
τ
3
42
+ λ
3
22
τ
3
22
+ λ
4
22
τ
4
22
+( μ
п4
+ μ
л2
)*ν /( k
п4
+k
л2
);
q:=c1[1]; r:=c2[1];
i:=2÷4
q<c1[i]
q:=c1[i]
да
нет
r<c2[i]
r:=c2[i]
да
нет
(19)
(15)
n=3:
c
1[1]:=
λ
2
11
τ
2
11
+(
μ
п1
+ μ
л3
)*ν
/(
k
п1
+
k
л3
);
c1[2]:= λ
1
12
τ
1
12
+ λ
3
33
τ
3
33
+( μ
п1
+ μ
л3
)*ν /( k
п1
+k
л3
);
c1[3]:= λ
1
13
τ
1
13
+ λ
3
13
τ
3
13
+ λ
3
33
τ
3
33
+ λ
4
33
τ
4
33
+( μ
п1
+ μ
л3
)*ν /( k
п1
+k
л3
);
c1[4]:= λ
2
31
τ
2
31
+( μ
п3
+ μ
л1
)*ν /( k
п3
+k
л1
);
c1[5]:= λ
1
32
τ
1
32
+ λ
3
13
τ
3
13
+( μ
п3
+ μ
л1
)*ν /( k
п3
+k
л1
);
c1[6]:= λ
1
33
τ
1
33
+ λ
3
33
τ
3
33
+ λ
3
13
τ
3
13
+ λ
4
13
τ
4
13
+( μ
п3
+ μ
л1
)*ν /( k
п3
+k
л1
);
c2[1]:= λ
2
21
τ
2
21
+( μ
п2
+ μ
л4
)*ν /( k
п2
+k
л4
);
c2[2]:= λ
1
22
τ
1
22
+ λ
3
43
τ
3
43
+( μ
п2
+ μ
л4
)*ν /( k
п2
+k
л4
);
c2[3]:= λ
1
23
τ
1
23
+ λ
3
23
τ
3
23
+ λ
3
43
τ
3
43
+ λ
4
43
τ
4
43
+( μ
п2
+ μ
л4
)*ν /( k
п2
+k
л4
);
c2[4]:= λ
2
41
τ
2
41
+ ( μ
п4
+ μ
л2
)*ν /( k
п4
+k
л2
);
c2[5]:= λ
1
42
τ
1
42
+ λ
3
23
τ
3
23
+( μ
п4
+ μ
л2
)*ν /( k
п4
+k
л2
);
c2[6]:= λ
1
43
τ
1
43
+ λ
3
43
τ
3
43
+ λ
3
23
τ
3
23
+ λ
4
23
τ
4
23
+( μ
п4
+ μ
л2
)*ν /( k
п4
+k
л2
);
q:=c1[1]; r:=c2[1];
(14)
n=2:
234
i:=2÷6
q<c1[i]
q:=c1[i]
да
нет
r<c2[i]
r:=c2[i]
да
нет
(19)
c
1[1]:=
λ
2
11
τ
2
11
+(
μ
п1
+ μ
л3
)*ν
/(
k
п1
+
k
л3
);
c1[2]:= λ
1
12
τ
1
12
+ λ
3
34
τ
3
34
+ λ
4
34
τ
4
34
+( μ
п1
+ μ
л3
)*ν /( k
п1
+k
л3
);
c1[3]:= λ
1
13
τ
1
13
+ λ
3
34
τ
3
34
+ λ
4
34
τ
4
34
+( μ
п1
+ μ
л3
)*ν /( k
п1
+k
л3
);
c1[4]:= λ
3
14
τ
3
14
+ λ
4
14
τ
4
14
+( μ
п1
+ μ
л3
)*ν /( k
п1
+k
л3
);
c1[5]:= λ
2
31
τ
2
31
+( μ
п3
+ μ
л1
)*ν /( k
п3
+k
л1
);
c1[6]:= λ
1
32
τ
1
32
+ λ
3
14
τ
3
14
+ λ
4
14
τ
4
14
+ ( μ
п3
+ μ
л1
)*ν /( k
п3
+k
л1
);
c1[7]:= λ
1
33
τ
1
33
+ λ
3
14
τ
3
14
+ λ
4
14
τ
4
14
+ ( μ
п3
+ μ
л1
)*ν /( k
п3
+k
л1
);
c1[8]:= λ
1
34
τ
1
34
+ λ
4
34
τ
4
34
+( μ
п3
+ μ
л1
)*ν /( k
п3
+k
л1
);
c2[1]:= λ
2
21
τ
2
21
+( μ
п2
+ μ
л4
)*ν /( k
п2
+k
л4
);
c2[2]:= λ
1
22
τ
1
22
+ λ
3
44
τ
3
44
+ λ
4
44
τ
4
44
+( μ
п2
+ μ
л4
)*ν /( k
п2
+k
л4
);
c2[3]:= λ
1
23
τ
1
23
+ λ
3
44
τ
3
44
+ λ
4
44
τ
4
44
+( μ
п2
+ μ
л4
)*ν /( k
п2
+k
л4
);
c2[4]:= λ
3
24
τ
3
24
+ λ
4
24
τ
4
24
+( μ
п2
+ μ
л4
)*ν /( k
п2
+k
л4
);
c2[5]:= λ
2
41
τ
2
41
+ ( μ
п4
+ μ
л2
)*ν /( k
п4
+k
л2
);
c2[6]:= λ
1
42
τ
1
42
+ λ
3
24
τ
3
24
+ λ
4
24
τ
4
24
+( μ
п4
+ μ
л2
)*ν /( k
п4
+k
л2
);
c2[7]:= λ
1
43
τ
1
43
+ λ
3
24
τ
3
24
+ λ
4
24
τ
4
24
+( μ
п4
+ μ
л2
)*ν /( k
п4
+k
л2
);
c2[8]:= λ
3
44
τ
3
44
+ λ
4
44
τ
4
44
+( μ
п4
+ μ
л2
)*ν /( k
п4
+k
л2
);
q:=c1[1]; r:=c2[1];
i:=2÷8
q<c1[i]
q:=c1[i]
да
нет
r<c2[i]
r:=c2[i]
да
нет
(19)
n=4:
(16)
235
c
1[1]:=
λ
2
11
τ
2
11
+(
μ
п1
+ μ
л3
)*ν
/(
k
п1
+
k
л3
);
c1[2]:= λ
1
12
τ
1
12
+ λ
2
12
τ
2
12
+( μ
п1
+ μ
л3
)*ν /( k
п1
+k
л3
);
c1[3]:= λ
1
13
τ
1
13
( μ
п1
+ μ
л3
)*ν /( k
п1
+k
л3
);
c1[4]:= λ
1
14
τ
1
14
+( μ
п1
+ μ
л3
)*ν /( k
п1
+k
л3
);
c1[5]:= λ
1
15
τ
1
15
+( μ
п1
+ μ
л3
)*ν /( k
п1
+k
л3
);
c1[6]:= λ
2
31
τ
2
31
+( μ
п3
+ μ
л1
)*ν /( k
п3
+k
л1
);
c1[7]:= λ
1
32
τ
1
32
+ λ
2
32
τ
2
32
+( μ
п3
+ μ
л1
)*ν /( k
п3
+k
л1
);
c1[8]:= λ
1
33
τ
1
33
+ ( μ
п3
+ μ
л1
)*ν /( k
п3
+k
л1
);
c1[9]:= λ
1
34
τ
1
34
+( μ
п3
+ μ
л1
)*ν /( k
п3
+k
л1
);
c1[10]:= λ
1
35
τ
1
35
+( μ
п3
+ μ
л1
)*ν /( k
п3
+k
л1
);
c2[1]:= λ
2
21
τ
2
21
+( μ
п2
+ μ
л4
)*ν /( k
п2
+k
л4
);
c2[2]:= λ
1
22
τ
1
22
+ λ
2
22
τ
2
22
+( μ
п2
+ μ
л4
)*ν /( k
п2
+k
л4
);
c2[3]:= λ
1
23
τ
1
23
+( μ
п2
+ μ
л4
)*ν /( k
п2
+k
л4
);
c2[4]:= λ
1
24
τ
1
24
+ ( μ
п2
+ μ
л4
)*ν /( k
п2
+k
л4
);
c2[5]:= λ
1
25
τ
1
25
+ ( μ
п2
+ μ
л4
)*ν /( k
п2
+k
л4
);
c2[6]:= λ
2
41
τ
2
41
+ ( μ
п4
+ μ
л2
)*ν /( k
п4
+k
л2
);
c2[7]:= λ
1
42
τ
1
42
+ λ
2
42
τ
2
42
+( μ
п4
+ μ
л2
)*ν /( k
п4
+k
л2
);
c2[8]:= λ
1
43
τ
1
43
+( μ
п4
+ μ
л2
)*ν /( k
п4
+k
л2
);
c2[9]:= λ
1
44
τ
1
44
+( μ
п4
+ μ
л2
)*ν /( k
п4
+k
л2
);
c2[9]:= λ
1
45
τ
1
45
+( μ
п4
+ μ
л2
)*ν /( k
п4
+k
л2
);
q:=c1[1]; r:=c2[1];
i:=2÷10
q<c1[i]
q:=c1[i]
да
нет
r<c2[i]
r:=c2[i]
да
нет
(17)
n=5:
(19)
236
(18)
n=6:
c
1[1]:=
λ
2
11
τ
2
11
+(
μ
п1
+ μ
л3
)*ν
/(
k
п1
+
k
л3
);
c1[2]:= λ
1
12
τ
1
12
+ λ
2
12
τ
2
12
+( μ
п1
+ μ
л3
)*ν /( k
п1
+k
л3
);
c1[3]:= λ
1
13
τ
1
13
( μ
п1
+ μ
л3
)*ν /( k
п1
+k
л3
);
c1[4]:= λ
1
14
τ
1
14
+( μ
п1
+ μ
л3
)*ν /( k
п1
+k
л3
);
c1[5]:= λ
1
15
τ
1
15
+( μ
п1
+ μ
л3
)*ν /( k
п1
+k
л3
);
c1[6]:= λ
1
16
τ
1
16
+( μ
п1
+ μ
л3
)*ν /( k
п1
+k
л3
);
c1[7]:= λ
2
31
τ
2
31
+( μ
п3
+ μ
л1
)*ν /( k
п3
+k
л1
);
c1[8]:= λ
1
32
τ
1
32
+ λ
2
32
τ
2
32
+( μ
п3
+ μ
л1
)*ν /( k
п3
+k
л1
);
c1[9]:= λ
1
33
τ
1
33
+ ( μ
п3
+ μ
л1
)*ν /( k
п3
+k
л1
);
c1[10]:= λ
1
34
τ
1
34
+( μ
п3
+ μ
л1
)*ν /( k
п3
+k
л1
);
c1[11]:= λ
1
35
τ
1
35
+( μ
п3
+ μ
л1
)*ν /( k
п3
+k
л1
);
c1[12]:= λ
1
36
τ
1
36
+( μ
п3
+ μ
л1
)*ν /( k
п3
+k
л1
);
c2[1]:= λ
2
21
τ
2
21
+( μ
п2
+ μ
л4
)*ν /( k
п2
+k
л4
);
c2[2]:= λ
1
22
τ
1
22
+ λ
2
22
τ
2
22
+( μ
п2
+ μ
л4
)*ν /( k
п2
+k
л4
);
c2[3]:= λ
1
23
τ
1
23
+( μ
п2
+ μ
л4
)*ν /( k
п2
+k
л4
);
c2[4]:= λ
1
24
τ
1
24
+ ( μ
п2
+ μ
л4
)*ν /( k
п2
+k
л4
);
c2[5]:= λ
1
25
τ
1
25
+ ( μ
п2
+ μ
л4
)*ν /( k
п2
+k
л4
);
c2[11]:= λ
1
25
τ
1
25
+ ( μ
п2
+ μ
л4
)*ν /( k
п2
+k
л4
);
c2[6]:= λ
2
41
τ
2
41
+ ( μ
п4
+ μ
л2
)*ν /( k
п4
+k
л2
);
c2[7]:= λ
1
42
τ
1
42
+ λ
2
42
τ
2
42
+( μ
п4
+ μ
л2
)*ν /( k
п4
+k
л2
);
c2[8]:= λ
1
43
τ
1
43
+( μ
п4
+ μ
л2
)*ν /( k
п4
+k
л2
);
c2[9]:= λ
1
44
τ
1
44
+( μ
п4
+ μ
л2
)*ν /( k
п4
+k
л2
);
c2[10]:= λ
1
45
τ
1
45
+( μ
п4
+ μ
л2
)*ν /( k
п4
+k
л2
);
c2[12]:= λ
1
46
τ
1
46
+( μ
п4
+ μ
л2
)*ν /( k
п4
+k
л2
);
q:=c1[1]; r:=c2[1];
i:=2÷12
q<c1[i]
q:=c1[i]
да
нет
r<c2[i]
r:=c2[i]
да
нет
(19)
(20)
237
α13
:=q/(q+r);
β13
:=1
-
α13
;
α24:=r/(q+r);β24:=1-α24;
(20)
Uslovie:=0
j:=1÷4
нет
j<=2
да
d:=
λ
2
j1
; w:=
τ
2
j1
;
f:= λ
2
j+21
; m:= τ
2
j+21
;
d< λ
2
jk
d:=λ
2
jk
да
нет
k:=2÷n
w< τ
2
jk
w:=τ
2
jk
да
нет
f< λ
2
j+2k
f:=λ
3
j+2k
да
нет
m< τ
3
j+2k
m:=τ
3
j+2k
да
нет
да
j mod 2 <>0
нет
N
j
:= (((1
-
q
-
r)*S
j
)/
/(1-q-r+α13-((μ
пj
+μ
лj+2
)*ν)/
/(k
пj
+k
лj+2
)+d*w+f*m))-2*tr
N
j
:= (((1
-
q
-
r)*S
j
)/
/(1-q-r+α24-((μ
пj
+μ
лj+2
)*ν)/
/(k
пj
+k
лj+2
)+d*w+f*m))-2*tr
да
j mod 2 <>0
нет
N
j
:= (((1
-
q
-
r)*S
j
)/
/(1-q-r+α13-((μ
пj
+μ
лj-2
)*ν)/
/(k
пj
+k
лj-2
)+d*w+f*m))-2*tr
N
j
:= (((1
-
q
-
r)*S
j
)/
/(1-q-r+α24-((μ
пj
+μ
лj-2
)*ν)/
/(k
пj
+k
лj-2
)+d*w+f*m))-2*tr
d:=
λ
2
j1
; w:=
τ
2
j1
;
f:= λ
2
j-21
; m:= τ
2
j-21
;
d< λ
2
jk
d:=λ
2
jk
да
нет
k:=2÷n
w< τ
2
jk
w:=τ
2
jk
да
нет
f< λ
2
j-2k
f:=λ
3
j-2k
да
нет
m< τ
3
j-2k
m:=τ
3
j-2k
да
нет
Nj <= T
0
Uslovie:=uslovie+1
да
нет
(21)
(22)
(23)
238
А. Алгоритм расчета основных характеристик СО с двухфазахным циклом
регулирования движения
Uslovie=4
да нет
(24)
(35)
q+r<>1
да нет
C:=(2*tr+T
0
)/(1-q-r)
ошибка
конец
t
зел13
:=
α
13
*(
С
-
T
0
);
t
зел24
:=β13*(С-T
0
);
t
кр
13
:=t
зел
24
; t
кр
24
:=t
зел
13
;
j:=1÷4
k:=1÷n
n
0цh
:=n
0цh
+ n
0jk
да
j mod 2 <>0
нет
n
0jk
:=(tr+T
0
+β13*(С-T
0
))/(1-q)*λ
jk
n
0jk
:=(tr+T
0
+β24*(С-T
0
))/(1-r)*λ
jk
n
0h
:=
n
0ц
h
*3600/С
(25)
(26)
(2
7
)
(28)
Z
1TE13
:=(2*tr+T
0
+
β13
*(
С
-
T
0
))/2;
Z
2TE24
:=(2*tr+T
0
+β24*(С-T
0
))/2;
n
13
:=0; n
24
:=0;
(29)
k:=1÷n
n
13
:=
n
13
+(n
01k
+
n
03k
);
n
24
:= n
24
+(n
02k
+n
04k
);
(30)
Z
13
:=
n
13
*
Z
1TE13
;
Z
24
:=n
24
*Z2
TE24
;
Z
цh
:=(Z
13
+Z
24
)/3600;
Z
h
:=(Z
13
+Z
24
)/C;
(31)
(32)
239
j:=1÷4
k:=1÷n
да
j mod 2 <>0
нет
t
0jk
:=((2*tr+T
0
+β13*(С-T
0
))*q)/(1-q) t
0jk
:=((2*tr+T
0
+β24*(С-T
0
))*q)/(1-q)
(34)
(35)
j:=1÷4
k:=1÷n
да
j mod 2 <>0
нет
t
xjk
:=t
зел13
- t
0jk
t
xjk
:=t
зел24
- t
0jk
f:=0
k:=1÷n
Z
jkЖив
:= Z
jkЖив
+ n
0jk
m:=1÷round(n
0jk
)
Z
jkЖив
:=
(2*
Z
jkЖив
-
2* n
0jk
*round(n
0jk
))/(2*λ
jk
)
Z
jЖив
:= Z
jЖив
+ Z
jkЖив
Z
1TEjЖив
:= Z
jЖив
/n
0j
; Z
hЖив
:= Z
hЖив
+ Z
jЖив
j:=1÷4
Z
цhЖив
:=
Z
hЖив
*
3600
/C
;
(32)
(33)
240