Добавил:
МТУСИ Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Лабораторная работа 6 / ЛР_6_Мягков_БАП2201.pdf
Скачиваний:
0
Добавлен:
21.04.2026
Размер:
640.1 Кб
Скачать

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