Определение оценок параметров закона распределения Вейбулла
Распределения Вейбулла занимает промежуточное положение между нормальным и экспоненциальным распределениями. В частном случае при b =1, с = 0 распределение Вейбулла переходит в экспоненциальное; при этом параметр а =Tср представляет собой среднее значение (математическое ожидание) случайной величины Т.
Анализ формул позволяет утверждать:
– при b > 1 интенсивность событий монотонно возрастает;
– при b = 1 интенсивность событий не изменяется во времени;
– при b < 1 интенсивность событий монотонно убывает.
Скорость изменения интенсивности событий определяется значением b. Эта особенность распределения позволяет использовать его для описания безотказности объектов в течение трех периодов их эксплуатации: приработки, нормальной эксплуатации и старения.
Распределение Вейбулла имеет ресурс для подшипников качения (b = 1,4÷1,5), зубчатых колес редукторов (b = 1,4 ÷ 1,8), зубчатых муфт (b =1,5), тормозных обкладок (b =1,4), тормозных шкивов (b = 1,5), ходовых колес кранов (b = 2) и т. д.
Распределение Вейбулла имеет две разновидности: двухпараметрическое и трехпараметрическое. Двухпараметрическое распределение получается из трехпараметрического при параметре сдвига с = 0
Для распределения Вейбулла [12] справедливо уравнение, задаваемое интегральной функции
где а — параметр масштаба, характеризующий степень растянутости кривой распределения вдоль оси x и связанный со средним значением случайной величины; b — параметр формы; с — параметр сдвига, являющийся минимальным возможным значением случайной величины x.
Дифференциальная функция [12] плотности распределения Вейбулла
Интенсивность событий
.
Математическое ожидание
.
Вычисление этого интеграла приводит к выражению
,
где гамма-функция.
Гамма-функцией от аргумента n называется функция
.
Оценку параметров масштаба а, формы b и сдвига с автотракторной техники удобно определять методом моментов [6,7,8] по выборке независимых наблюдений случайной величины X.
Определение производится по значению коэффициента асимметрии с помощью специальных таблиц, в следующем порядке:
– вычисляют оценку коэффициента асимметрии по формуле:
(3.8)
для группированных данных оценку коэффициента асимметрии св вычисляют по формуле:
(3.9)
– по полученному значению рв из таблицы 1 приложения 2 при помощи линейной интерполяции находят оценку параметра b и значения коэффициентов gb и Kb;
– определяют оценку для параметра (а) по формуле:
(3.10)
– находят значение c по формуле:
(3.11)
в качестве оценки параметра с берут одно из двух значений, взятых по абсолютной величине :
(3.12)
Задача 7
Определить оценки параметров исправного состояния тормозной системы из закона Вейбулла для данных из задачи 2.
Пример решения
Вычислить значение коэффициента асимметрии по формуле (3.9):
По таблице 1 приложения 2 с помощью линейной интерполяции найти значения оценок параметра b и значения коэффициентов gb и Kb.
Интерполяцию произвести по схеме Эйткина, по зависимости значений параметра формы от коэффициента асимметрии в двух точках
взятых так, что .
Найти с помощью схемы Эйткина значение b для ρb = 0,68 по формуле:
3.13
Аналогичным образом найти значения коэффициентов gb и Kb:
Определить значение параметра масштаба а:
Найти значение c :
Так как 54,6 < 70, принимаем с = 54,6.