- •ТЕРМИНЫ И ОПРЕДЕЛЕНИЯ
- •перечень сокращений и обозначений
- •Введение
- •1. ПРОВЕРКА ГИПОТЕЗЫ О ПРИНАДЛЕЖНОСТИ ОПЫТНЫХ ДАННЫХ ВЫБРАННОМУ ЗАКОНУ РАСПРЕДЕЛЕНИЯ
- •1.1 Расчёт параметров статистического ряда распределения
- •1.2 Определение числовых характеристик наработок
- •1.3. Построение гистограммы
- •1.4. Нормирование случайной велечины
- •1.6. Расчёт критерия согласия χ² пирсона
- •1.7. Расчёт интегральных функций распределения P(t) и F(t)
- •ЗАКЛЮЧЕНИЕ
- •Контрольные вопросы
- •СПИСОК ИСПОЛЬЗОВАННЫХ ИСТОЧНИКОВ
1. ПРОВЕРКА ГИПОТЕЗЫ О ПРИНАДЛЕЖНОСТИ ОПЫТНЫХ ДАННЫХ ВЫБРАННОМУ ЗАКОНУ РАСПРЕДЕЛЕНИЯ
1.1 Расчёт параметров статистического ряда распределения
Рассчитываем параметры статистического ряда распределения:
–размах выборки R = tmax−tmin= 276 − 61 = 215 тыс. км;
–число интервалов k = 1 + 3,32·lg(44) = 6,44 ≈ 7;
–величину интервала h =
≈ 35 тыс. км (для удобства расчётов
принимаем h = 35 тыс. км и строим интервалы, начиная от значения, кратного
5);
–границы интервалов: 60; 95; 130; 165; 200; 235; 270; 305;
–середины интервалов (тыс. км): t1 = 77,5; t2 = 112,5; t3 = 147,5; t4 = 182,5; t5 = 217,5; t6 = 252,5; t7 = 287,5;
–частоту попаданий наработок в эти интервалы: m1 = 1; m2 = 8; m3 = 10; m4 = 18; m5 = 3; m6 = 2; m7 = 2.
Результаты группировки сведены в таблицу 1.
Таблица 1 – Статистический ряд распределения наработок до предельного
состояния |
|
|
|
|
Номер |
Границы интервалов |
Середина интервала |
Частота miоп |
|
интервала |
(ti − ti+1), тыс. км |
ti, тыс. км |
|
|
1-й |
60 |
– 95 |
77,5 |
1 |
2-й |
95 – 130 |
112,5 |
8 |
|
3-й |
130 |
– 165 |
147,5 |
10 |
4-й |
165 |
– 200 |
182,5 |
18 |
5-й |
200 |
– 235 |
217,5 |
3 |
6-й |
235 |
– 270 |
252,5 |
2 |
7-й |
270 |
– 305 |
287,5 |
2 |
Итого |
|
– |
– |
44 |
1.2 Определение числовых характеристик наработок
Определяем числовые характеристики наработок до предельного состояния (по формулам для генеральной совокупности, как в эталоне):
– средняя наработка tcp = 166 тыс. км;
7
–среднеквадратическое отклонение σ(t) = 33,1 тыс. км;
–коэффициент вариации v =
= 0,2.
Вид гистограммы (скос вправо) и значение коэффициента вариации v = 0,2 < 0,33 позволяют предположить, что наработки до предельного состояния объекта распределены по нормальному закону.
1.3. Построение гистограммы
Строим гистограмму распределения частоты наработок объекта до предельного состояния (рисунок 1).
Рисунок 1 – Гистограмма распределения наработок объекта до предельного состояния по интервалам наработки
1.4. Нормирование случайной велечины
Для удобства вычислений проводим нормирование случайной величины наработки t, переходя к величине zi =
, и находим границы новых интервалов. Расчёты сводим в таблицу 2.
8
Таблица 2 – Границы интервалов случайной величины z |
|
|
||||
Интервал |
ti |
ti + 1 |
ti − tcp |
ti + 1 − tcp |
zi = |
zi+1 = |
1-й |
60 |
95 |
-106 |
-71 |
-3,20 |
-2,14 |
2-й |
95 |
130 |
-71 |
-36 |
-2,14 |
-1,09 |
3-й |
130 |
165 |
-36 |
-1 |
-1,09 |
-0,03 |
4-й |
165 |
200 |
-1 |
34 |
-0,03 |
1,03 |
5-й |
200 |
235 |
34 |
69 |
1,03 |
2,08 |
6-й |
235 |
270 |
69 |
104 |
2,08 |
3,14 |
7-й |
270 |
305 |
104 |
139 |
3,14 |
4,20 |
1.5. Расчёт теоретических частот попадания наработок в |
||||||
интервалы |
|
|
|
|
|
|
Рассчитываем теоретические частоты miтеор = N·Pi,где Pi |
= Ф(zi +1)− |
|||||
Ф(zi). Результаты расчета сведены в таблицу 3.
Таблица 3 – Теоретические частоты попадания наработок в интервалы
Интервал |
zi |
zi+1 |
Ф(zi) |
Ф(zi+1) |
Pi = Ф(zi+1) |
miтеор = |
|
|
|
|
|
− Ф(zi) |
N·Pi |
1-й |
-3,20 |
-2,14 |
0,0007 |
0,0162 |
0,0155 |
0,68 |
2-й |
-2,14 |
-1,09 |
0,0162 |
0,1379 |
0,1217 |
5,35 |
3-й |
-1,09 |
-0,03 |
0,1379 |
0,4880 |
0,3501 |
15,40 |
4-й |
-0,03 |
1,03 |
0,4880 |
0,8485 |
0,3605 |
15,86 |
5-й |
1,03 |
2,08 |
0,8485 |
0,9812 |
0,1327 |
5,84 |
6-й |
2,08 |
3,14 |
0,9812 |
0,9992 |
0,0180 |
0,79 |
7-й |
3,14 |
4,20 |
0,9992 |
1,0000 |
0,0008 |
0,04 |
Итого |
– |
– |
– |
– |
1,0000 |
44,0 |
1.6. Расчёт критерия согласия χ² пирсона
Всоответствии с правилами применения критерия χ² Пирсона, опытные
итеоретические частоты в крайних интервалах должны быть не менее 5. В связи с этим производим объединение смежных интервалов:
– Объединяем 1-й и 2-й интервалы: mоп = 1 + 8 = 9; mтеор = 0,68 + 5,35 =
6,03;
– Объединяем 5-й, 6-й и 7-й интервалы: mоп = 3 + 2 + 2 = 7; mтеор = 5,84 +
9
0,79 + 0,04 = 6,67.
После объединения количество интервалов стало равно k = 4.
Расчет критерия согласия χ² для объединенных интервалов сведён в таблицу 4.
Таблица 4 – Результаты расчета критерия согласия χ²
Объединенный |
miоп |
miтеор |
miоп − miтеор |
(miоп − miтеор)2 |
|
интервал |
|
|
|
|
|
1-й и 2-й |
9 |
6,03 |
2,97 |
8,8209 |
1,4630 |
3-й |
10 |
15,40 |
-5,40 |
29,1600 |
1,8935 |
4-й |
18 |
15,86 |
2,14 |
4,5796 |
0,2888 |
5-й, 6-й и 7-й |
7 |
6,67 |
0,33 |
0,1089 |
0,0163 |
СУММА |
44 |
43,96 |
– |
– |
χ²опыт = 3,662 |
Определяем число степеней свободы S = k − r − 1 = 4 − 2 − 1 = 1.
При уровне значимости α = 0,05 и числе степеней свободы S = 1 табличное значение χ2табл = 3,84.
Так как χ2опыт< χ2табл(3,662 < 3,84), расхождение между опытными и теоретическими частотами признаётся незначимым. Принимается гипотеза о принадлежности выборочных данных нормальному закону распределения.
1.7. Расчёт интегральных функций распределения P(t) и F(t)
Рассчитываем вероятности безотказной работы P(t) и вероятности отказов F(t) на наработках, соответствующих серединам исходных
интервалов ti:P(t) = 0,5 − Ф
; F(t) = 1 − P(t), Результаты расчётов
сводим в таблицу 5.
Таблица 5 – Результаты расчета интегральных функций распределения P(t) и
F(t) по интервалам наработки |
|
|
|
|
|
||
Функция |
1-й |
2-й |
3-й |
4-й |
5-й |
6-й |
7-й |
(t = 77,5) (t = 112,5) (t = 147,5) (t = 182,5)(t = 217,5) (t = 252,5) (t = 287,5) |
|||||||
P(t) |
0,9993 |
0,9838 |
0,8621 |
0,5115 |
0,1515 |
0,0188 |
0,0008 |
F(t) |
0,0007 |
0,0162 |
0,1379 |
0,4885 |
0,8485 |
0,9812 |
0,9992 |
По |
данным |
таблицы |
5 строим графики интегральных |
функций |
|||
10
распределения (рисунок 2).
Рисунок 2 – Графики функций вероятности безотказной работы P(t) и вероятности отказов F(t)
11
