Экспоненциальное (показательное) распределение
Рассмотрим семейство распределений, широко используемое при принятии управленческих решений и других прикладных исследованиях - семейство экспоненциальных распределений. Проанализируем вероятностную!! модель, приводящую к таким распределениям. Для этого рассмотрим «поток событий», т.е. последовательность событий, происходящих одно за другим в какие-то моменты времени. Примерами могут служить: время безотказной работы компьютерной системы, интервал между последовательными поступлениями автомобилей к стон-линии перекрестка, поток обращений клиентов в отделение банка; поток покупателей, обращающихся за товарами и услугами; поток вызовов на телефонной станции; поток отказов оборудования в технологической цепочке и т.д.
В теории потоков событий справедлива теорема суммировании потоков событий. Суммарный поток состоит из большого количества независимых частных потоков, ни один из которых не оказывает преобладающего влияния на суммарный поток. Так, поток вызовов, поступающих на телефонную станцию, состоит из большого числа независимых потоков вызовов, исходящих от отдельных абонентов. В случае, когда характеристики потоков не зависят от времени, суммарный поток полностью описывается одним числом X - интенсивностью потока. Для суммарного потока функция распределения случайной величины X - длины промежутка времени между последовательными событиями имеет следующий вид:
Это распределение называется экспоненциальным (показательным) распределением. В данную функцию иногда вводят параметр сдвига с.
Экспоненциальное распределение имеет только один параметр, который и определяет его характеристики. Плотность распределения имеет следующий вид:
где X - постоянная положительная величина.
График функции /(х) представлен на рис. 9.12.
Рис. 9.12.
На рис. 9.13 представлен график плотности экспоненциального распределения при разных параметрах X.
Экспоненциальное распределение характеризует распределение времени между независимыми событиям, появляющимися с постоянной интенсивностью. Экспоненциальный закон характерен для распределения случайных величин, изменение которых обусловлено влиянием какого-то доминирующего фактора. В теории надежности это распределение описывает распределение внезапных отказов, так как последние являются редкими событиями. Экспоненциальное распределение служит также для описания
Рис. 9.13. Плотность экспоненциального распределения при разных параметрах X
наработки сложных систем, прошедших период приработки, и для описания времени безотказной работы системы с большим числом последовательно соединенных элементов, каждый из которых не оказывает большого влияния на отказ системы.
Теоретические частоты для экспоненциального закона распределения определяют по формуле
где N - объем совокупности; 1г к - длина интервала; е - основание натурального логарифма; X - условные отклонения середин классов:
Рассмотрим выравнивание эмпирического распределения (табл. 9.4) по экспоненциальному закону.
Таблица 9.4
Эмпирические частоты для выравнивания распределения по экспоненциальному закону
Имеем N = 160; Ь к = 41; х = 54,59. Расчет величин условных отклонений середин классов, вспомогательных величин е _1 и теоретических частот произведен в табл. 9.5.
Таблица 95
Выравнивание эмпирических частот по экспоненциальному закону
Эмпирические данные, х |
Эмпирическая частота, т |
Теоретические частоты |
|||
Эмпирические и теоретические частоты экспоненциального распределения изобразим графически на рис. 9.14.
Показательное распределение представляет собой частный случай распределения Вейбулла - Гнеденко (соответствующий значению параметра формы b = 1).
Экспоненциальный закон распределения называемый также основным законом надежности, часто используют для прогнозирования надежности в период нормальной эксплуатации изделий, когда постепенные отказы еще не проявились и надежность характеризуется внезапными отказами. Эти отказы вызываются неблагоприятным стечением многих обстоятельств и поэтому имеют постоянную интенсивность. Экспоненциальное распределение находит довольно широкое применение в теории массового обслуживания, описывает распределение наработки на отказ сложных изделий, время безотказной работы элементов радиоэлектронной аппаратуры.
Приведем примеры неблагоприятного сочетания условий работы деталей машин, вызывающих их внезапный отказ. Для зубчатой передачи это может быть действием максимальной нагрузки на наиболее слабый зуб при его зацеплении; для элементов радиоэлектронной аппаратуры - превышение допустимого тока или температурного режима.
Плотность распределения экспоненциального закона (рис. 1) описывается соотношением
f (x ) = λe −λ x ; (3)
функция распределения этого закона - соотношением
F (x ) = 1− e −λ x ; (4)
функция надежности
P (x ) = 1− F (x ) = e −λ x ; (5)
математическое ожидание случайной величины Х
дисперсия случайной величины Х
(7)
Экспоненциальный закон в теории надежности нашел широкое применение, так как он прост для практического использования. Почти все задачи, решаемые в теории надежности, при использовании экспоненциального закона оказываются намного проще, чем при использовании других законов распределения. Основная причина такого упрощения состоит в том, что при экспоненциальном законе вероятность безотказной работы зависит только от длительности интервала и не зависит от времени предшествующей работы.
Риc. 1. График плотности экспоненциального распределения
Пример 2. По данным эксплуатации генератора установлено, что наработка на отказ подчиняется экспоненциальному закону с параметром λ=2*10 -5 ч -1 . Найти вероятность безотказной работы за время t =100 ч. Определить математическое ожидание наработки на отказ.
Р е ш е н и е. Для определения вероятности безотказной работы воспользуемся формулой (5), в соответствии с которой
Математическое ожидание наработки на отказ равно
Как было сказано ранее, примерами распределений вероятностей непрерывной случайной величины Х являются:
- равномерное распределение вероятностей непрерывной случайной величины;
- показательное распределение вероятностей непрерывной случайной величины;
- нормальное распределение вероятностей непрерывной случайной величины.
Дадим понятие равномерного и показательного законов распределения, формулы вероятности и числовые характеристики рассматриваемых функций.
Показатель | Раномерный закон распределения | Показательный закон распределения |
---|---|---|
Определение | Равномерным называется распределение вероятностей непрерывной случайной величины X, плотность которого сохраняет постоянное значение на отрезке и имеет вид | Показательным (экспоненциальным) называется распределение вероятностей непрерывной случайной величины X, которое описывается плотностью, имеющей вид |
где λ – постоянная положительная величина |
||
Функция распределения | ||
Вероятность попадания в интервал | ||
Математическое ожидание | ||
Дисперсия | ||
Среднее квадратическое отклонение |
Примеры решения задач по теме «Равномерный и показательный законы распределения»
Задача 1.
Автобусы идут строго по расписанию. Интервал движения 7 мин. Найти: а) вероятность того, что пассажир, подошедший к остановке, будет ожидать очередной автобус менее двух минут; б) вероятность того, что пассажир, подошедший к остановке, будет ожидать очередной автобус не менее трех минут; в) математическое ожидание и среднее квадратическое отклонение случайной величины X – времени ожидания пассажира.
Решение. 1. По условию задачи непрерывная случайная величина X={время ожидания пассажира} равномерно распределена между приходами двух автобусов. Длина интервала распределения случайной величины Х равна b-a=7, где a=0, b=7.
2. Время ожидания будет менее двух минут, если случайная величина X попадает в интервал (5;7).
Вероятность попадания в заданный интервал найдем по формуле:
Р(х 1 <Х<х 2)=(х 2 -х 1)/(b-a)
.
Р(5 < Х < 7) = (7-5)/(7-0) = 2/7 ≈ 0,286.
3. Время ожидания будет не менее трех минут (т.е. от трех до семи мин.), если случайная величина Х попадает в
интервал (0;4). Вероятность попадания в заданный интервал найдем по формуле:
Р(х 1 <Х<х 2)=(х 2 -х 1)/(b-a)
.
Р(0 < Х < 4) = (4-0)/(7-0) = 4/7 ≈ 0,571.
4. Математическое ожидание непрерывной, равномерно распределенной случайной величины X – времени ожидания пассажира, найдем по формуле: М(Х)=(a+b)/2 . М(Х) = (0+7)/2 = 7/2 = 3,5.
5. Среднее квадратическое отклонение непрерывной, равномерно распределенной случайной величины X – времени ожидания пассажира, найдем по формуле: σ(X)=√D=(b-a)/2√3 . σ(X)=(7-0)/2√3=7/2√3≈2,02.
Задача 2.
Показательное распределение задано при x ≥ 0 плотностью f(x) = 5e – 5x. Требуется: а) записать выражение для функции распределения; б) найти вероятность того, что в результате испытания X попадает в интервал (1;4); в) найти вероятность того, что в результате испытания X ≥ 2 ; г) вычислить M(X), D(X), σ(X).
Решение. 1. Поскольку по условию задано показательное распределение , то из формулы плотности распределения вероятностей случайной величины X получаем λ = 5. Тогда функция распределения будет иметь вид:
2. Вероятность того, что в результате испытания X попадает в интервал (1;4) будем находить
по формуле:
P(a < X < b) = e −λa − e −λb
.
P(1 < X < 4) = e −5*1 − e −5*4 = e −5 − e −20 .
3. Вероятность того, что в результате испытания X ≥ 2 будем находить по формуле:
P(a < X < b) = e −λa − e −λb при a=2, b=∞.
Р(Х≥2) = P(1< X < 4) = e −λ*2 − e −λ*∞ = e −2λ − e −∞ =
e −2λ - 0 = e −10 (т.к. предел e −х при х стремящемся к ∞ равен нулю).
4. Находим для показательного распределения:
- математическое ожидание по формуле M(X) =1/λ = 1/5 = 0,2;
- дисперсию по формуле D(X) = 1/ λ 2 = 1/25 = 0,04;
- среднее квадратическое отклонение по формуле σ(Х) = 1/λ = 1/5 = 1,2.
где λ – постоянная положительная величина.
Из выражения (3.1), следует, чтопоказательное распределение определяется одним параметром λ.
Эта особенность показательного распределения указывает на его преимущество по сравнению с распределениями , зависящими от большего числа параметров. Обычно параметры неизвестны и приходится находить их оценки (приближенные значения) разумеется, проще оценить один параметр, чем два или три и т. д . Примером непрерывной случайной величины, распределенной по показательному закону , может служить время между появлениями двух последовательных событий простейшего потока.
Найдем функцию распределения показательного закона .
Итак
Графики плотности и функции распределения показательного закона изображены на рис. 3.1.
Учитывая, что получим:
Значения функции можно находить по таблице.
Числовые характеристики показательного распределения
Пусть непрерывная случайная величина Χ распределена по показательному закону
Найдем математическое ожидание , используя формулу её вычисления для непрерывной случайной величины:
Следовательно:
Найдем среднее квадратическое отклонение , для чего извлечем квадратный корень из дисперсии:
Сравнивая (3.4), (3.5) и (3.6), видно, что
т. е. математическое ожидание и среднее квадратическое отклонение показательного распределения равны между собой.
Показательное распределение широко применяетсяв различных приложениях финансовых и технических задач, например, в теории надежности.
4. Распределение «хи-квадрат» и распределение Стьюдента.
4.1 Распределение «хи-квадрат» (- распределение)
Пусть Χ i (ί = 1, 2, ..., n)-нормальные независимые случайные величины , причем математическое ожиданиекаждой из нихравно нулю , а среднее квадратическое отклонение - единице .
Тогдасумма квадратов этих величин
распределена по закону с степенями свободы , если же эти величины связаны одним линейным соотношением, например , то число степеней свободы
Распределение хи-квадрат нашло широкое применение в математической статистике.
Плотность этого распределения
где - гамма-функция, в частности .
Отсюда видно, чтораспределение хи-квадрат определяется одним параметром - числом степеней свободы k.
С увеличением числа степеней свободыраспределение хи-квадрат медленно приближается к нормальному.
Хи-квадрат распределение получается, если в законе распределения Эрланга принять λ = ½ и k = n /2 – 1.
Математическое ожидание и дисперсия случайной величины, имеющей хи-квадрат распределение, определяются простыми формулами, которые приведем без вывода:
Из формулы следует, что при хи-квадрат распределение совпадает с экспоненциальным распределением при λ = ½ .
Интегральная функция распределения при хи-квадрат распределенииопределяетсячерез специальные неполные табулированные гамма-функции
На рис.4.1. приведены графики плотности вероятности и функции распределения случайной величины, имеющей хи-квадрат распределениепри n = 4, 6, 10.Рис.4.1. а )Графики плотности вероятности при хи-квадрат распределении
Рис.4.1. б)Графики функции распределения при хи-квадрат распределении
4.2 Распределение Стьюдента
Пусть Z – нормальная случайная величина, причём
а V – независимая от Z величина, которая распределена по закону хи-квадрат с k степенями свободы.Тогда величина:
имеет распределение, которое называют t -распределением или распределением Стьюдента (псевдоним английского статистика В. Госсета),
с k = n - 1 степенями свободы (n - объём статистической выборки при решении задач статистки).
Итак , отношение нормированной нормальной величинык квадратному корню из независимой случайной величины, распределенной по закону «хи квадрат» с k степенями свободы , деленной на k, распределено по закону Стьюдента с k степенями свободы.
Плотность распределения Стьюдента:
Определение. Показательным (экспоненциальным) называется распределение вероятностей непрерывной случайной величины Х, которое описывается плотностью
где l - положительное число.
Найдем закон распределения.
Графики функции распределения и плотности распределения:
f(x) F(x)
Найдем математическое ожидание случайной величины, подчиненной показательному распределению.
Результат получен с использованием того факта, что
Для нахождения дисперсии найдем величину М(Х 2).
Дважды интегрируя по частям, аналогично рассмотренному случаю, получим:
Тогда
Итого: Видно, что в случае показательного распределения математическое ожидание и среднее квадратическое отклонение равны.
Также легко определить и вероятность попадания случайной величины, подчиненной показательному закону распределения, в заданный интервал.
Показательное распределение широко используется в теории надежности .
Допустим , некоторое устройство начинает работать в момент времени t 0 =0 , а через какое - то время t происходит отказ устройства.
Обозначим Т непрерывную случайную величину - длительность безотказной работы устройства.
Таким
образом
, функция распределения F(t) = P(T
Вероятность противоположного события (безотказная работа в течение времени t ) равна R(t) = P(T>t) = 1 - F(t).
Определение. Функцией надежности R(t) называют функцию, определяющую вероятность безотказной работы устройства в течение времени t .
Часто на практике длительность безотказной работы подчиняется показательному закону распределению.
Вообще говоря , если рассматривать новое устройство, то вероятность отказа в начале его функционирования будет больше, затем количество отказов снизится и будет некоторое время иметь практически одно и то же значение. Затем (когда устройство выработает свой ресурс) количество отказов будет возрастать.
Другими словами , можно сказать, что функционирование устройства на протяжении всего существования (в смысле количества отказов) можно описать комбинацией двух показательных законов (в начале и конце функционирования) и равномерного закона распределения.
Функция надежности для какого- либо устройства при показательном законе распределения равна:
Данное соотношение называют показательным законом надежности .
Важным свойством , позволяющим значительно упростить решение задач теории надежности, является то, что вероятность безотказной работы устройства на интервале времени t не зависит от времени предшествующей работы до начала рассматриваемого интервала, а зависит только от длительности времени t .
Таким образом , безотказная работа устройства зависит только от интенсивности отказов l и не зависит от безотказной работы устройства в прошлом.
Так как подобным свойством обладает только показательный закон распределения, то этот факт позволяет определить, является ли закон распределения случайной величины показательным или нет.
2.8 Распределение «Хи-квадрат»
Пусть X i (i=1,2,…,n) - нормальные независимые случайные величины, причем математическое ожидание каждой из них равно нулю, а среднее квадратическое отклонение - единице. Тогда сумма квадратов этих величин
распределена по закону («Хи-квадрат») с k=n степенями свободы; если же эти величины связаны одним линейным соотношением, например , то число степеней свободы k=n-1.
Плотность этого распределения
где -Гамма-функция; в частности,
Отсюда видно , что распределение «Хи-квадрат» определяется одним параметром - числом степеней свободы k. С увеличением числа степеней свободы распределение медленно приближается к нормальному.
2.9 Распределение Стьюдента
Пусть Z -нормальная случайная величина, причем M(Z)=0, s(Z)=1, а V- независимая от Z величина, которая распределена по закону с k степенями свободы. Тогда величина
имеет распределение, которое называют t- распределением или распределением Стьюдента, k степенями свободы. Итак отношение нормированной нормальной величины к квадратному корню из независимой случайной величины, распределенной по закону
«Хи-квадрат» с k степенями свободы , деленной на k, деленной на k распределено по закону Стьюдента с k степенями свободы. . С увеличением числа степеней свободы распределение медленно приближается к нормальному.
2.9 Нормальный закон распределения
Определение. Нормальным называется распределение вероятностей непрерывной случайной величины, которое описывается плотностью вероятности
Нормальный закон распределения также называется законом Гаусса .
Нормальный закон распределения занимает центральное место в теории вероятностей. Это обусловлено тем, что этот закон проявляется во всех случаях, когда случайная величина является результатом действия большого числа различных факторов. К нормальному закону приближаются все остальные законы распределения.
Можно легко показать , что параметры и , входящие в плотность распределения являются соответственно математическим ожиданием и средним квадратическим отклонением случайной величины Х.
Найдем функцию распределения F(x) .
График плотности нормального распределения называется нормальной кривой или кривой Гаусса .
Нормальная кривая обладает следующими свойствами:
1 ) Функция определена на всей числовой оси.
2 ) При всех х функция распределения принимает только положительные значения.
3 ) Ось ОХ является горизонтальной асимптотой графика плотности вероятности, т.к. при неограниченном возрастании по абсолютной величине аргумента х , значение функции стремится к нулю.
4 ) Найдем экстремум функции.
Т.к. при y’ > 0 при x < m и y’ < 0 при x > m , то в точке х = т функция имеет максимум, равный .
5 ) Функция является симметричной относительно прямой х = а , т.к. разность
(х - а ) входит в функцию плотности распределения в квадрате.
6 ) Для нахождения точек перегиба графика найдем вторую производную функции плотности.
При x = m + s и x = m - s вторая производная равна нулю, а при переходе через эти точки меняет знак, т.е. в этих точках функция имеет перегиб.