- •Реферат
- •Введение и постановка задач исследования
- •1. Теоретические основы методов моментов и максимального правдоподобия
- •1.1. Оценка первых четырех выборочных начальных и центральных моментов и построение эмпирического распределения
- •1.2. Метод моментов
- •1.3 Метод максимального правдоподобия (ммп)
- •.4 О методе «вероятностной бумаги»
- •2. Алгоритмы расчета параметров распределений методом моментов
- •2.1. Алгоритмы расчета параметров гамма-распределения методом моментов и его развитие
- •.2. Алгоритмы расчета параметров распределения Вейбулла
- •. Алгоритмы расчета параметров распределений методом максимального правдоподобия
- •3.1. Алгоритмы расчета параметров гамма-распределения
- •3.2. Алгоритмы расчета параметров распределения Вейбулла для полной выборки
- •3.3. Алгоритмы расчета параметров распределения Вейбулла для цензурированной выборки
- •4 Статистическая проверка гипотез. Критерий «стареющих» распределений
- •4.1 Критерий согласия а.Н. Колмогорова
- •4.2 Критерии «стареющего» распределения
- •5. Программы и результаты решения конкретной задачи обработки статистических данных
- •5.1. Постановка задачи
- •5.2. Расчет первых четырех выборочных моментов и построение эмпирической функции распределения
- •{Расчет центральных моментов по формулам теоретической связи}
- •5.3. Расчет параметров предполагаемых теоретических распределений методом моментов
- •5.4. Расчет параметров предполагаемых теоретических распределениях методом максимального правдоподобия
- •5.5 Развитие метода максимального правдоподобия
- •5.5 А. Метод вероятностной бумаги График функции распределения можно представить в виде совокупности точек (t, p) на плоскости, где
- •5.6 Сводка результатов, их анализ и выводы
- •1 Экспоненциальное распределение
- •2 Гамма-распределение
- •3 Распределение Эрланга
- •4 Распределение Вейбулла-Гнеденко
- •6 Распределение экстремальных (максимальных) значений
- •Метод вероятностной бумаги График функции распределения можно представить в виде совокупности точек (t, p) на плоскости, где
1 Экспоненциальное распределение
Плотность
вероятности
,
где
,
.
Функция
распределения
Вероятность
безотказной работы
Интенсивность
отказов
.
Соотношения между моментами и параметром распределения "λ":
Для
среднего времени наработки до отказа
.
Для
дисперсии и среднеквадратичного
отклонения
;
.
Для
коэффициентов асимметрии и эксцесса
;
.
Для
медианы
.
Коэффициент вариации
.
2 Гамма-распределение
Плотность
вероятности наработки до отказа
где
— параметр масштаба (
),
— параметр формы (
),
— гамма-функция или эйлеров интеграл
второго рода
или
.
Аналитического выражения для функции распределения наработки на отказ не существует (аналитическое выражения для нее существует только для целых положительных значений параметра ; см. ниже распределение Эрланга).
Известны формулы связи моментов с параметрами "a" и "λ" гамма-распределения:
; ; ; ; .
Коэффициент
вариации при этом
.
Мода:
для значений
.
Квантиль
находится из уравнения
для
Точка перегиба
Начальные
моменты таковы
В
задачах обработки статистических
экспериментальных данных потребуются
выражения и оценки численных значений
производной
и функции
или в развернутом виде
.
Нетрудно
показать, что
=
,
а
функция примет вид
.
При
—
целом (
)
функция
принимает
значения
,
то есть,
...
и т.д. Таким образом, гамма-функция Эйлера
- это распространение функции или
операции «факториал» на случай нецелых
чисел, в том числе и отрицательных.
Табличные
значения функции
в
справочниках ограничены значениями
аргумента "a",
принадлежащими интервалу
.
Значения
гамма-функции для
,
но при этом
... ,
При больших
значениях a
(
)
по формуле
,
которые следуют из функционального
уравнения Эйлера
.
Примеры:
;
;
Полезные соотношения и значения гамма-функции:
;
.
Функция распределения времени наработки до отказа, как отмечалось выше, не имеет аналитического выражения. В общем виде она может быть представлена таким образом
где
— неполная гамма-функция.
при
:
; при
:
;
Примечание. Из гамма-распределения «вытекают»:
при
— экспоненциальное распределение;
при
и a,
кратном
, будем иметь χ2-распределение
(при этом
—
число степеней свободы);
при a — целом: a = 1; 2; ... ; k; ... — распределение Эрланга.
Существенное уменьшение вычислительных трудностей может быть достигнуто применением асимптотических разложений (аппроксимационных формул) Стирлинга:
.
Ряд
Стирлинга полезен для больших
.
Для действительных положительных a
абсолютные
величины ошибки меньше, чем абсолютная
величина последнего из взятых членов.
Другие
полезные аппроксимации Стирлинга:
при
.
.
Известно важное неравенство:
<
<
.
3 Распределение Эрланга
Плотность распределения наработки до отказа:
для
;
;
— целое.
Функция распределения времени наработки до отказа:
Вероятность безотказной работы:
.
Интенсивность
отказов системы в целом
Соотношения между моментами и параметрами распределения определяются как и у гамма-распреденения, но с заменой параметра на .
Распределение
Эрланга порядка k
описывает распределение случайной
величины
как суммы k
штук независимых случайных величин,
каждая из которых распределена по
показательному (экспоненциальному)
закону с параметром λ.
Распределению
Эрланга удовлетворяет время наработки
до отказа резервированной системы с
включением «холодного» резерва по
способу замещения при условии, что
наработка до отказа включенного элемента
подчинена экспоненциальному закону.
При этом
,
где m
— число резервных элементов. Из
соотношения
вытекает свойство структур с «холодным»
резервом – средняя наработка системы
до отказа линейно возрастает от числа
резервных элементов.
