Построение графика статистической функции распределения. Вариационный ряд. Построение гистограммы наработок между отказами

1. Исходные данные:

Наработка до отказа имеет распределение приведенное в таблице 1. Определить гамма-процентный ресурс машины с использованием вероятностных сеток.

Таблица 1 – Наработка до отказа, ч.

2. Построение графика статистической функции распределения

Удобным способом получить представление о распределении случайной величины Х является построение графика статистической функции распределения выборки (х₁, … х_n).

,                                                                                                                  (1)

где q – число опытов, в которых случайная величина Х принимала значения меньше х;

n – общее число произведенных опытов.

Статистический ряд перестраивается так, чтобы полученные числовые значения случайной величины располагались в возрастающем порядке (вариационный ряд).

Для возрастающего ряда эмпирическая функция распределения записывается:

 ,                                                   (2)

где i – порядковый номер опыта.

Таблица 2 – Вариационный ряд

3. Построение гистограммы наработок между отказами.

Другим видом обработки статистического ряда является поразрядная группировка полученных из опытов значений случайной величины. Из данных статистического ряда для непрерывных случайных величин строится гистограмма, изображающая статистическую плотность распределения. Она позволяет получить первое представление о виде распределения.

Высота разряда:

;                                                                                                                (3)

где P^*_i – частость;

l_i – длина разряда.

Частость:

,                                                                                                                (4)

где q_i – число наблюдений в разряде (частота).

Длина разряда:

,                                                                                                                (5)

где  - максимальное значение случайной величины в разряде i;

 -  минимальное значение случайной величины в разряде i.

Таблица 3.

4. Построение графика функции распределения на вероятностной сетке.

Экспоненциальное распределение

Функция распределения случайной величины:

.                                                         (6)

Логарифмируя это выражение, получим линейную зависимость:

.                                                   (7)

Переходим к построению вероятностных шкал F(t).

На горизонтальной оси наносится равномерная шкала t с масштабным фактором K_t :

.                                                               (8)

На вертикальной оси F(t) наносятся значения –ln [1-F(t)], а надписываются значения F(t). Шкала неравномерная. Наименьшее значение F(t) равно нулю. Наибольшее значение удобно принять равным 0,999. Тогда –ln (1-0,999)=6,908.

При длине вертикальной оси 150 мм,

.           (9)

Параметр распределения:

.                  (10)

Распределение Вейбула

Функция распределения случайной величины:

.                                                  (11)

Величина S_Fдля вертикальной шкалы длиной 150 мм определяется логарифмированием.

Примем F(t) крайние значения 0,001 и 0,999. Для этих значений размах составит 8,84  (от -6,91 до 1,93).

Числовые значения S_Fприведены в справочнике.

Горизонтальная шкала неравномерная и строится также, как для логарифмически нормального распределения.

Параметр b:

.                                                       (12)

Параметр а определяется точкой пересечения найденной прямой с осью t.

Графики на вероятностных сетках представлены на рисунке 3 и рисунке 4.

Полученная гистограмма наработок машины (см. рис. 2) позволяет сделать предположение о том, что наработки между случайными отказами подчиняются закону Вейбула, плотность распределения которого:

                                        (13)

где t – значение случайной величины T (наработка до отказа);

a – параметр масштаба (из графика а=990);

b – параметр формы ().

А закон распределения будет иметь вид:

                                                   (14)

Вероятность попадания случайной величины в i – й разряд:

              (15)

Подставляя в последнюю формулу а=990 ; b=1,8 и значения α_i и β_iдля каждого разряда из таблицы 3, получим теоретические значения p_i , приведенные в таблице 4.

Таблица 4.- Теоретические значения вероятностей.

Критерий Пирсона:

;                                         (16)

.

Распределение χ² зависит от параметра r, называемого числом степеней свободы.

r=k-(ψ+1)                                                              (17)

где k – число разрядов статистического ряда (k=5);

ψ – число параметров предполагаемого распределения (ψ=0).

r=5-0-1=4.

По справочнику для χ² =2,2 и r=4 находим вероятность p=0,7358 , что больше 0,05.

Следовательно предполагаемый закон Вейбула не противоречит опытным