ГОУ ВПО
ДАЛЬНЕВОСТОЧНЫЙ ГОСУДАРСТВЕННЫЙ
УНИВЕРСИТЕТ ПУТЕЙ СООБЩЕНИЯ
кафедра
«Управление
эксплуатационной работой»
Расчетно-графическая работа
«Обработка статистических данных и
установление закона распределения
случайных величин»
Выполнил: Роднов Е.А.
студент 232 группы
Проверил: Виноградова Л.Л.
Хабаровск, 2007
В табл. 1 приведены моменты фактического прибытия поездов на станцию, полученные в результате натурных наблюдений.
Необходимо:
1. Составить статистический ряд интервалов прибытия грузовых поездов на сортировочную станцию.
2. Установить основные временные параметры входящего на станцию поездопотока (среднее значение, дисперсию, среднее квадратическое отклонение и коэффициент вариации).
3. Построить гистограмму распределения вероятностей появления интервалов между поездами.
4. Определить параметр в эрланговском распределении интервалов прибытия поездов в парк приема.
5. По внешнему виду гистограммы подобрать теоретический закон распределения.
6. Проверить по критерию согласия Пирсона и условию Романовского правдоподобность гипотезы о выбранном теоретическом распределении интервалов прибытия поездов в расформирование.
7. Рассчитать число бригад ПТО в парке приема (время на техническое обслуживание t
то
принять равным 20 мин).
Таблица 1
Моменты фактического прибытия поездов на станцию
№ п/п |
Время прибытия (ч, мин) |
№ п/п |
Время прибытия (ч, мин) |
№ п/п |
Время прибытия (ч, мин) |
№ п/п |
Время прибытия (ч, мин) |
№ п/п |
Время прибытия (ч, мин) |
1 |
0-04 |
19 |
6-01 |
37 |
12-43 |
55 |
20-45 |
73 |
4-02 |
2 |
0-10 |
20 |
6-26 |
38 |
13-00 |
56 |
22-45 |
74 |
4-20 |
3 |
0-25 |
21 |
7-12 |
39 |
14-10 |
57 |
22-49 |
75 |
4-30 |
4 |
0-45 |
22 |
7-22 |
40 |
14-22 |
58 |
22-54 |
76 |
4-55 |
5 |
1-15 |
23 |
7-40 |
41 |
14-24 |
59 |
22-59 |
77 |
5-08 |
6 |
1-30 |
24 |
8-01 |
42 |
14-50 |
60 |
23-15 |
78 |
5-15 |
7 |
2-01 |
25 |
8-15 |
43 |
14-55 |
61 |
23-25 |
79 |
5-22 |
8 |
2-26 |
26 |
9-35 |
44 |
16-25 |
62 |
23-47 |
80 |
6-15 |
9 |
2-40 |
27 |
9-45 |
45 |
16-35 |
63 |
23-51 |
81 |
7-05 |
10 |
2-45 |
28 |
9-53 |
46 |
16-38 |
64 |
0-04 |
82 |
7-25 |
11 |
2-50 |
29 |
10-05 |
47 |
16-50 |
65 |
0-35 |
83 |
7-35 |
12 |
3-48 |
30 |
10-15 |
48 |
17-40 |
66 |
0-48 |
84 |
7-44 |
13 |
3-52 |
31 |
10-26 |
49 |
18-05 |
67 |
1-01 |
85 |
8-10 |
14 |
3-58 |
32 |
10-34 |
50 |
18-10 |
68 |
1-48 |
86 |
8-21 |
15 |
4-15 |
33 |
11-05 |
51 |
18-36 |
69 |
2-38 |
87 |
8-27 |
16 |
4-40 |
34 |
11-37 |
52 |
18-50 |
70 |
3-01 |
88 |
8-38 |
17 |
5-09 |
35 |
12-04 |
53 |
19-58 |
71 |
3-12 |
89 |
10-00 |
18 |
5-43 |
36 |
12-20 |
54 |
20-11 |
72 |
3-50 |
Решение
Интервалы прибытия определил путем вычитания предыдущего времени прибытия поезда из последующего и представил в табл. 2.
Таблица 2
Интервалы прибытия поездов
№ п/п |
Интервал, мин |
№ п/п |
Интервал, мин |
№ п/п |
Интервал, мин |
№ п/п |
Интервал, мин |
№ п/п |
Интервал, мин |
№ п/п |
Интервал, мин |
№ п/п |
Интервал, мин |
№ п/п |
Интервал, мин |
№ п/п |
Интервал, мин |
1 |
6
|
11 |
58
|
21 |
10
|
31 |
8
|
41 |
26
|
51 |
14
|
61 |
22
|
71 |
38
|
81 |
20
|
2 |
15
|
12 |
4
|
22 |
18
|
32 |
31
|
42 |
5
|
52 |
68
|
62 |
4
|
72 |
12
|
82 |
10
|
3 |
20
|
13 |
6
|
23 |
21
|
33 |
32
|
43 |
90
|
53 |
13
|
63 |
13
|
73 |
18
|
83 |
9
|
4 |
30
|
14 |
17
|
24 |
14
|
34 |
27
|
44 |
10
|
54 |
34
|
64 |
31
|
74 |
10
|
84 |
26
|
5 |
15
|
15 |
25
|
25 |
80
|
35 |
16
|
45 |
3
|
55 |
120
|
65 |
13
|
75 |
25
|
85 |
11
|
6 |
31
|
16 |
29
|
26 |
10
|
36 |
23
|
46 |
12
|
56 |
4
|
66 |
13
|
76 |
13
|
86 |
6
|
7 |
25
|
17 |
34
|
27 |
8
|
37 |
17
|
47 |
50
|
57 |
5
|
67 |
47
|
77 |
7
|
87 |
11
|
8 |
14
|
18 |
18
|
28 |
12
|
38 |
70
|
48 |
25
|
58 |
5
|
68 |
50
|
78 |
7
|
88 |
82
|
9 |
5
|
19 |
25
|
29 |
10
|
39 |
12
|
49 |
5
|
59 |
16
|
69 |
23
|
79 |
53
|
10 |
5
|
20 |
46
|
30 |
11
|
40 |
2
|
50 |
26
|
60 |
10
|
70 |
11
|
80 |
50
|
Группировка происходит по классам (разрядам). Количество классов К определил по формуле:
К = (1 + 3,21 ∙ lgn
), (1)
где n
– общее число наблюдений.
К = (1 + 3,21 ∙ lg88) = 7, 24.
Принимаем количество классов К равным 8.
Величину интервала (шаг класса) группирования Iопределил по формуле:
, (2)
где хmax
, хmin
– наибольшее и наименьшее значения случайной величины.
.
Далее произвел группирование интервалов по разрядам. В процессе группирования установил, сколько интервалов mi
попало в разряд ti
– ti+1
. Последующие расчеты основных параметров статистического ряда выполнил в форме табл. 3, в которую свел все промежуточные результаты вычислений.
Таблица 3
Обработка статистического ряда интервалов между моментами
прибытия поездов на станцию
№ п/п |
Границы разрядов, ti
- ti+1
|
Число интервалов в разряде, mi
|
ЧастостьPi
|
Среднее значение в разряде, ti
|
ti
∙ P |
ti
2
∙ P |
1 |
2-16,75 |
46 |
0,523 |
9,375 |
4,901 |
45,943 |
2 |
16,75-31,5 |
25 |
0,284 |
24,125 |
6,854 |
165,345 |
3 |
31,5-46,25 |
5 |
0,057 |
38,875 |
2,209 |
85,867 |
4 |
46,25-61 |
6 |
0,068 |
53,625 |
3,656 |
196,066 |
5 |
61-75,75 |
2 |
0,023 |
68,375 |
1,554 |
106,253 |
6 |
75,75-90,5 |
3 |
0,034 |
83,125 |
2,834 |
235,560 |
7 |
90,5-105,25 |
0 |
0 |
97,875 |
0 |
0 |
8 |
105,25-120 |
1 |
0,011 |
112,625 |
1,280 |
144,141 |
Итого |
88 |
1 |
23,287 |
979,176 |
Для каждого разряда наблюдаемых величин подсчитал их количество и определил частость.
Математическое ожидание М (х)
определил по формуле:
, (3)
где – среднее значение разряда i
.
М (х)
= 9,375 ∙ 0,523 + 24,125 ∙ 0,284 + 38,875 ∙ 0,057 + 53,625 ∙ 0,068 + 68,375 ∙ 0,023 + 83,125 ∙ 0,034 + 97,875 ∙ 0 + 112,625 ∙ 0,011 = 23,287.
Дисперсию D
(
x
)
определил по формуле:
D
(
x
) =
M
2
(
x
) – (
M
[
x
])2
, (4)
где – второй начальный момент случайной величины.
M
2
(
x
)
= 9,3752
∙ 0,523 + 24,1252
∙ 0,284 + 38,8752
∙ 0,057 + 53,6252
∙ 0,068 + 68,3752
∙ 0,023 + 83,1252
∙ 0,034 + 97,8752
∙ 0 + 112,6252
∙ 0,011 = 979,176.
D
(
x
) =
979,176 – 23,2872
= 436,892.
Среднее квадратическое отклонение – это квадратный корень из дисперсии:
. (5)
= 20,902.
Коэффициент вариации определяется как отношение среднего квадратического отклонение к математическому ожиданию:
. (6)
.
На основании расчетных характеристик (табл. 3) строится гистограмма распределения интервалов прибытия поездов (рис. 1).
Ординату гистограммы определил по формуле:
, (7)
где ∆
ti
=
ti
+1
–
ti
шаг конкретного i
-го разряда.
Рис. 1. Гистограмма распределения интервалов прибытия поездов
Параметр Эрланга определил по следующей формуле:
, (8)
.
Исходя из найденного значения параметра Эрланга и внешнего вида гистограммы, сделал предположение, что для данного распределения наиболее близок закон показательного теоретического распределения.
Теоретическая вероятность Р*i
интервалов определенной величины в их общей совокупности равна:
Р*i
= F(ti
) – F(ti
–
1
), (9)
где F(ti
) и F(ti
–
1
) –
функция показательного распределения.
F(ti
) = 1 –
е–
λ
ti
,
F(ti–
1
) = 1 –
е–
λ
ti–
1
,
где λ – интенсивность поступления поездов на станцию
. (10)
.
Рассчитанные данные занесены в табл. 4.
Таблица 4
Характеристики распределения интервалов
между поступающими в переработку поездами
Границы разрядов, ti
- ti+1
|
Число интервалов в разряде, mi
|
λt |
е–
λt
|
F(ti
) |
Р*i
|
nР*i
|
mi
–nP*i
|
(mi
–nP*i
)2
|
(mi
–nP*i
)2
/nР*i
|
2-16,75 |
46 |
0,086 |
0,918 |
0,082 |
0,429 |
37,731 |
8,269 |
68,376 |
1,8122 |
16,75-31,5 |
25 |
0,719 |
0,489 |
0,511 |
0,229 |
20,113 |
4,887 |
23,885 |
1,1875 |
31,5-46,25 |
5 |
1,353 |
0,261 |
0,739 |
0,122 |
10,721 |
-5,721 |
32,733 |
3,0531 |
46,25-61 |
6 |
1,986 |
0,139 |
0,861 |
0,065 |
5,715 |
0,285 |
0,081 |
0,0142 |
61-75,75 |
2 |
2,619 |
0,074 |
0,926 |
0,035 |
3,046 |
-1,046 |
1,095 |
0,3595 |
75,75-90,5 |
3 |
3,253 |
0,040 |
0,960 |
0,018 |
1,624 |
1,376 |
1,894 |
1,1660 |
90,5-105,25 |
0 |
3,886 |
0,021 |
0,979 |
0,010 |
0,866 |
-0,866 |
0,749 |
0,8657 |
105,25-120 |
1 |
4,520 |
0,011 |
0,989 |
0,005 |
0,461 |
0,539 |
0,290 |
0,6286 |
5,153 |
0,006 |
0,994 |
88 |
9,0867 |
Как видно из последней графы табл. 4 критерий согласия Пирсона χ2
= 9,0867.
Число степеней свободы rопределил по формуле:
r = R– S, (11)
где R– число разрядов; S – число наложенных связей.
r = 8 – 2 = 6.
Пользуясь специальной таблицей, определяется вероятность Р(χ2
) = 0,1736. Значит, гипотеза о показательном распределении интервалов поступления поездов на станцию не совсем правдоподобна.
По условию Романовского, гипотеза о принятом теоретическом законе распределения считается правдоподобной, если соблюдается следующее неравенство:
. (12)
, следовательно, расхождение между теоретическим и эмпирическим распределениями не столь существенно и гипотеза о показательном законе распределения интервала прибытия поездов на станцию правдоподобна.
Число бригад ПТО в парке приема рассчитал, исходя из условия:
, (13)
где Ip
– расчетный интервал прибытия поездов; t
ТО
– время на техническое обслуживание поезда одной бригадой; Б
– число бригад.
, (14)
где Imin
– минимальный интервал между поездами, прибывающими на станцию, Imin
= 2 мин.
Icp
= M(t)
= 23, 287.
= 12,64 мин.
Время на техническое обслуживание t
ТО
принимается равным 20 мин.
Число бригад рассчитал из формулы (13) и округлил до целого числа, следовательно, принял 2 бригады ПТО.
|