Автор работы: Пользователь скрыл имя, 15 Декабря 2011 в 13:52, курсовая работа
В данной курсовой работе были обработаны статистические данные наработки до отказа турбобура в часах. На основании чего было выбрано закон распределения Вейбулла и построены графики дифференциальной, интегральной, обратно интегральной функции и функции интенсивности.
Определим
математическое ожидание:
Определим среднеквадратичное
отклонение:
Определим значение
дисперсии:
Определим значение
коэффициента вариации:
Проверим член
124 и 127 по критерию Романовского по формуле:
Проверяем t=127:
Проверяем t=124:
Следовательно
член 124 и 127 исключаем из дальнейшего рассмотрения.
3.3 Проверка резко
выделяющихся величин по критерию Ирвина
по формуле:
Следовательно, анализируемые величины оставляем при дальнейшем рассмотрении.
3.4 Проверка резко
выделяющихся величин по критерию Груббса
по формуле:
Так как для обеих точек при n=189 заведомо (таблица 5 приложения), то оставляем крайние точки в рассматриваемой совокупности.
Таблица №7 – расчетная таблица
№ | Интервал,
ч |
∆t | Середина | n*i | p*i |
1 | 0-10 | 10 | 5 | 23 | 0,121693 |
2 | 10-20 | 10 | 15 | 20 | 0,10582 |
3 | 20-30 | 10 | 25 | 19 | 0,100529 |
4 | 30-40 | 10 | 35 | 20 | 0,10582 |
5 | 40-50 | 10 | 45 | 20 | 0,10582 |
6 | 50-60 | 10 | 55 | 17 | 0,089947 |
7 | 60-70 | 10 | 65 | 15 | 0,079365 |
8 | 70-80 | 10 | 75 | 23 | 0,121693 |
9 | 80-90 | 10 | 85 | 8 | 0,042328 |
10 | 90-100 | 10 | 95 | 10 | 0,05291 |
11 | 100-110 | 10 | 105 | 7 | 0,037037 |
12 | 110-120 | 10 | 115 | 7 | 0,037037 |
Определим частность
в i-ом интервале:
4 Построение графиков статистических функций(гистограмм)
4.1Определим значение
интегральной функции распределения:
4.2Определим эмпирическую
плотность распределения:
4.3Определим значение
обратной интегральной функции распределения:
4.4Определим значение
функции интенсивности:
Сведем полученные данные в таблицу:
Таблица №8
N | ∆t | ticp | n*i | |||||
1 | 10 | 5 | 23 | 0,121693 | 0,012169 | 0,121693 | 0,878307 | 0,013855 |
2 | 10 | 15 | 20 | 0,10582 | 0,010582 | 0,227513 | 0,772487 | 0,013699 |
3 | 10 | 25 | 19 | 0,100529 | 0,010053 | 0,328042 | 0,671958 | 0,014961 |
4 | 10 | 35 | 20 | 0,10582 | 0,010582 | 0,433862 | 0,566138 | 0,018692 |
5 | 10 | 45 | 20 | 0,10582 | 0,010582 | 0,539683 | 0,460317 | 0,022989 |
6 | 10 | 55 | 17 | 0,089947 | 0,008995 | 0,62963 | 0,37037 | 0,024286 |
7 | 10 | 65 | 15 | 0,079365 | 0,007937 | 0,708995 | 0,291005 | 0,027273 |
8 | 10 | 75 | 23 | 0,121693 | 0,012169 | 0,830688 | 0,169312 | 0,071875 |
9 | 10 | 85 | 8 | 0,042328 | 0,004233 | 0,873016 | 0,126984 | 0,033333 |
10 | 10 | 95 | 10 | 0,05291 | 0,005291 | 0,925926 | 0,074074 | 0,071429 |
11 | 10 | 105 | 7 | 0,037037 | 0,003704 | 0,962963 | 0,037037 | 0,1 |
12 | 10 | 115 | 7 | 0,037037 | 0,003704 | 1 | 0 | ∞ |
Рисунок №1 – гистограмма эмпирической плотности распределения
Рисунок №2 – гистограмма интегральной функции распределения
Рисунок №3 – гистограмма
обратной интегральной функции
Рисунок №4 – гистограмма функции интенсивности
5. Выбор теоретического
закона распределения.
При обработке статистического материала важной задачей является подбор теоретического
закона распределения,
наилучшим образом описывающего
статистическое распределение.
Теоретический
закон подбирают, принимая во внимание:
- физическую природу отказов;
- опыт отработки деталей и изделий аналогичного назначения;
- форму кривой плотности распределения;
- совпадение опытных точек с теоретической кривой интегральной функции безотказности;
- коэффициент
вариации.
Значение коэффициента вариации, характеризующего рассеивание показателя надежности,
уже позволяет
судить об условиях эксплуатации машин
и их технологии изготовления.
5.1 Определим
коэффициент вариации:
Так как V=0,64,
то выбираем закон распределения Вейбулла.
5.2 Определим параметры
распределения Вейбулла.
Так как коэффициент
вариации V=0,64, то по таблице 2 приложения
определим находим:
b=1.6
5.3 Для определения
параметра a нужно найти Кb.
По таблице находим
0,640=V;
0,897=Kb;
Следовательно
уравнения закона распределения
Вейбулла примут вид:
6.Рассчитаем все теоретические функции.
6.1 определим значение интегральной функции распределения:
6.2 Определим
дифференциальную функцию распределения:
6.3 Расчет обратной
интегральной функции распределения:
6.4 Расчет
значение функции интенсивности отказов:
Найдем
разность между функциями F*(t)
и F(t):
Определим
вероятность попадания
Сведем полученные данные в таблицу:
Таблица №6
t | f(t) | P(t) | pi | F*(t) | F(t) | λ(t) | ||
0 | 0 | 1 | 0 | 0 | ||||
10 | 23 | 0,001795 | 0,936 | 0,0635 | 0,121693 | 0,064 | 0,001917 | 0,058 |
20 | 20 | 0,004762 | 0,819 | 0,1169 | 0,227513 | 0,181 | 0,005812 | 0,047 |
30 | 19 | 0,007596 | 0,683 | 0,1362 | 0,328042 | 0,317 | 0,011119 | 0,011 |
40 | 20 | 0,009633 | 0,547 | 0,1364 | 0,433862 | 0,453 | 0,017618 | 0,019 |
50 | 20 | 0,010624 | 0,422 | 0,1247 | 0,539683 | 0,578 | 0,025178 | 0,038 |
60 | 17 | 0,010619 | 0,315 | 0,1069 | 0,62963 | 0,685 | 0,033707 | 0,055 |
70 | 15 | 0,009837 | 0,228 | 0,086 | 0,708995 | 0,772 | 0,043135 | 0,063 |
80 | 23 | 0,008565 | 0,160 | 0,0676 | 0,830688 | 0,840 | 0,053409 | 0,009 |
90 | 8 | 0,007075 | 0,110 | 0,05 | 0,873016 | 0,890 | 0,064485 | 0,017 |
100 | 10 | 0,005581 | 0,073 | 0,0365 | 0,925926 | 0,927 | 0,076326 | 0,001 |
110 | 7 | 0,004225 | 0,048 | 0,0255 | 0,962963 | 0,952 | 0,0889 | 0,010 |
120 | 7 | 0,003081 | 0,030 | 0,0173 | 1 | 0,970 | 0,102179 | 0,030 |
f(t) - дифференциальная функция распределения (эмпирическая плотность распределения)
P(t) - обратная интегральной функции распределения
pi - вероятность попадания случайной величины в i-ый интервал
F(t) - интегральная функция распределения
λ(t) - функция интенсивности отказов
|D|
- разность между функциями F*(t)
и F(t)
Рисунок №5 – гистограмма и график теоретической
эмпирической плотности
Рисунок №6 – гистограмма и график теоретической интегральной функции
Рисунок №7 – гистограмма и график теоретической обратной интегральной функции распределения
Рисунок №8 – гистограмма и график теоретической функции интенсивности отказов
7. Проверка
гипотезы о соответствии эмпирического
и теоретического распределений с помощью
критериев согласия.
Назначение критерия
χ2 - критерия Пирсона
Критерий χ2 применяется
в двух целях:
1) для сопоставления
эмпирического распределения
2) для сопоставления
двух, трех или более эмпирических
распределений одного и того
же признака (в скрипте до 10).
Описание критерия
Критерий χ2 отвечает
на вопрос о том, с одинаковой ли
частотой встречаются разные значения
признака в эмпирическом и теоретическом
распределениях или в двух и более эмпирических
распределениях.
Преимущество
метода состоит в том, что он позволяет
сопоставлять распределения признаков,
представленных в любой шкале, начиная
от шкалы наименований. В самом простом
случае альтернативного распределения
"да - нет", "допустил брак - не допустил
брака", "решил задачу - не решил задачу"
и т. п. мы уже можем применить критерий
χ2.