Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Конспект часть 1 (1)

.pdf
Скачиваний:
14
Добавлен:
07.03.2015
Размер:
397.42 Кб
Скачать

Согласно этому предположению значения случайной величины меньше х1 в

опытах не наблюдались, т.е. для x <

 

 

F * (x) = 0 . Далее при x =

 

происходит

x

x

1

1

скачкообразное увеличение значения функции распределения на величину Р*1,

которое сохраняется до x = x2 . При этом значении х функция распределения уве-

личивается до F * (x2 ) = P1* + P22 . Аналогичным образом график функции распре-

деления строится для остальных интервалов. Последний скачок функции на вели- чину Р*k происходит в точке x = xk . Здесь функция распределения достигает свое-

го максимального значения F(x) = 1. Это соответствует тому, что в опытах не наблюдались величины x > xk .

Статистическая плотность распределения представляется в виде гистограм-

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

ния прямоугольников нужно частоту каждого интервала разделить на его длину и полученное число взять в качестве высоты прямоугольника. Только в случае рав- ных по длине разрядов высоты прямоугольников пропорциональны соответству-

ющим частотам. Из способа построения гистограммы следует, что полная пло-

щадь ее равна единице. Пример построения гистограммы приведен на рис.4.2.

f(x)

0,06

0,05

0,04

0,03

0,02

0,01

0

 

 

 

 

x

 

 

 

 

 

 

 

 

 

0 10 20 30 40 50

Рис.4.2. Статистическая плотность распределения случайной

величины х

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

23

кривой, соответствующей графику плотности распределения случайной величины

Х.

В тех случаях, когда не требуется полной информации о законе распределе-

ния случайной величины, бывает достаточно ограничиться некоторыми числен- ными характеристиками. По аналогии с теорией вероятности, в математической

статистике для представления о поведении случайной величины наиболее часто

применяются характеристика положения (среднее арифметическое - аналог ма- тематическому ожиданию) и характеристики разброса значений случайной вели-

чины относительного среднего положения (выборочные дисперсия и среднеквад- ратичное отклонение - аналоги дисперсии и среднеквадратичного отклонения ге- неральной совокупности).

Среднее арифметическое наблюдаемой величины рассчитывается по форму-

ле

 

 

 

n

 

 

 

 

å xi

 

mx* =

 

=

i =1

.

(3.17)

x

 

 

 

 

n

 

Эту характеристику иногда называют статистическим средним случайной ве-

личины. При достаточно большом n величина x сходится по вероятности к мате-

матическому ожиданию. При ограниченном числе опытов статистическое среднее является случайной величиной, которая, тем не менее, связана с математическим ожиданием и может дать о нем известное представление.

Для того, чтобы различать дисперсию и среднее квадратичное генеральной

совокупности и подобные характеристики выборочной совокупности, обычно пер-

вые величины обозначают буквами D и σ, а выборочные характеристики - буквами

S2 и S.

Выборочная дисперсия рассчитывается по формуле

 

ån (xi

 

)2

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

S2 =

i =1

 

 

,

 

 

 

 

 

(3.18)

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

выборочное среднеквадратичное отклонение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ån

(xi

 

)2

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

S = S2 =

i =1

 

 

 

.

(3.19)

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

Теоретически показано, что средним значением S2 является не дисперсия ге-

неральной совокупности Dх, а несколько меньшая величина

n 1

Dx . Таким обра-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

24

зом, величина S2 будет сходиться по вероятности не к теоретическому значению

Dх, а к несколько меньшему значению nn-1Dx . При этом оценка величин D и σ по

значениям S2 и S дает некоторое смещение результатов, которое можно устра- нить, введя соответствующие поправки:

 

 

 

 

 

 

 

 

 

ån (xi -

 

)2

 

 

 

 

 

 

 

n

 

x

 

 

 

 

2 = S2

 

=

i =1

 

 

 

,

(3.20)

S

 

 

 

 

n -1

 

 

n -1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ån (xi -

 

)2

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

=

 

i =1

 

 

 

 

.

 

 

 

(3.21)

S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n -1

 

 

 

 

 

 

 

 

 

Величина S 2 называется несмещенной оценкой дисперсии или несмещенной выборочной дисперсией, величина S называется несмещенным средним квадра- тичным отклонением.

Расчет значения выборочной дисперсии удобно производить по формуле

 

 

 

 

n

 

1 æ

n

ö2

 

 

 

 

 

 

å xi2

-

 

çç

å xi ÷÷

 

 

 

 

 

 

n

 

 

 

 

2

=

i =1

 

èi =1

ø

,

(3.22)

S

 

 

 

n -1

 

 

 

 

 

 

 

 

 

которая соответствует формуле (3.20).

25

Лекция 6. ЗАКОНЫ РАСПРЕДЕЛЕНИЯ НЕПРЕРЫВНЫХ СЛУЧАЙНЫХ ВЕЛИЧИН

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

функция распределения F(x). Для непрерывной случайной величины используют- ся интегральный закон, функция распределения F(x) и дифференциальная форма

плотность распределения f(x). На практике эти законы проявляются в виде реа- лизации определенных событий, причем данные о распределении случайной ве- личины существуют только в виде набора пар чисел: значение случайной величи-

ны частота появления этого значения.

Для того, чтобы проводить прогностические расчеты, необходимо иметь ана- литическую запись этих законов в виде определенных математических формул. Вид математического описания выбирается исходя из теоретических представле-

ний о распределении случайной величины или на основании экспериментальных данных. Соответствие аналитических зависимостей и экспериментальных данных

достигается путем подбора численных значений коэффициентов, входящих в ап-

проксимирующие формулы. Эти коэффициенты называются параметрами рас- пределения случайной величины. Для их подбора обычно применяется принцип максимального правдоподобия, частным случаем которого является метод наименьших квадратов.

Поскольку существует различие в описании распределения непрерывных и дискретных случайных величин, то и для их аппроксимации применяются различ-

ные формулы.

Для непрерывных случайных величин наиболее простым является равномер-

ное (прямоугольное) распределение с двумя параметрами {a, b}, соответствую-

щим границам интервала, внутри которого может находиться случайная величина. Это распределение встречается в основном в двух типовых ситуациях: во-первых,

когда в некотором интервале все значения случайной величины равновозможны и, во-вторых, при аппроксимации других непрерывных распределений в относи-

тельно малых интервалах.

Равномерное распределение на интервале [a, b] задается плотностью

1

 

f (x) = b a .

(4.1)

Функция распределения на этом интервале

26

x

1

x

x - a

 

F(x)= ò f (x)dx =

òdx =

.

 

 

a

b - a a

b - a

Математическое ожидание соответствует середине интервала:

mx = bò x f (x)dx =

b2 - a2

= a + b .

 

 

2(b - a)

 

 

a

 

 

2

 

 

 

 

 

 

 

 

Дисперсия распределения

 

 

 

 

 

 

Dx = bò(x - mx )2 f (x)dx =

1

bò(x - mx )2 d(x -mx )=

(b - a)2

.

 

12

a

b - a a

 

 

Графики функции и плотности распределения показаны на рис.4.1.

(4.2)

(4.3)

(4.4)

F(x)

f(x)

 

1

a)

x

б)

x

a

b

a

b

Рис.4.1. Графики функции а) и плотности б) распределения

случайной величины по прямоугольному закону

Некоторые случайные величины имеют экспоненциальное распределение с

одним параметром l0>0.

 

 

Функция распределения

 

 

F(x) = 1- e−λ0 x .

(4.5)

Плотность распределения

 

 

f (x)=

dF(x)

= λ e−λ0 x .

(4.6)

 

 

dx

0

 

 

 

 

Закон определен для случайной величины x ³ 0. Математическое ожидание

случайной величины

 

 

 

 

1

 

 

 

 

mx = ò x λ0e−λ0 xdx =

 

,

 

 

 

 

 

 

 

0

 

 

 

 

λ0

 

 

 

дисперсия

 

 

 

 

 

 

 

 

 

 

æ

1

ö2

λ0 e

−λ0 x

 

1

 

Dx = ò

ç

÷

 

 

ç x -

 

÷

 

 

 

dx =

 

.

λ0

 

 

 

λ02

0

è

ø

 

 

 

 

 

 

(4.7)

(4.8)

27

По этому закону, как правило, распределяется случайная величина продол-

жительность работы оборудования до отказа, если его отказы не связаны со ста- рением, т.е. износом оборудования, ухудшением определяющих его работоспо-

собность характеристик. Это распределение наблюдается, когда отказы вызыва- ются случайными внешними факторами, например, колебаниями тепловой, элек- трической или механической нагрузок, изменением качества топлива, особенно-

стями работы различных операторов или другими случайными воздействиями.

Графики функции и плотности экспоненциального распределения показаны на

рис.4.2.

1

F(x)

f(x)

 

 

 

a)

x

б)

x

 

Рис.5.2. Графики функции а) и плотности б) распределения случайной

величины по экспоненциальному закону

Наиболее общим является распределение Гаусса или нормальное распреде-

ление. Этот закон занимает особое положение. Доказано, что сумма достаточно большого числа независимых или слабо зависимых случайных величин, подчи-

ненных каким угодно законам распределения, приближенно подчиняется нор-

мальному закону распределения. По этой причине нормальный закон распреде- ления наиболее часто встречается в природе.

Большинство встречающихся на практике случайных величин, таких, как, например, ошибки измерений, могут быть представлены как суммы большого чис- ла малых элементарных ошибок, каждая из которых вызвана действием отдель- ной причины, не зависящей от остальных. Каким бы законом распределения ни подчинялись отдельные элементарные ошибки, особенности этих распределений в сумме нивелируются и сумма оказывается подчиненной закону, близкому к нор-

мальному.

Нормальный закон распределения характеризуется плотностью вероятности:

f (x) =

1

 

 

(x m)2

 

 

 

e

2σ 2

.

(4.9)

 

 

σ

 

 

 

 

2π

 

 

28

 

Кривая плотности распределения по нор-

f(x)

мальному закону, приведенная на рис.4.3,

 

имеет симметричный холмообразный вид.

х

Максимальная ордината кривой, равная

1

 

m

σ 2π , соответствует точке x=m. По мере

Рис.4.3. График плотности рас-

удаления от точки x=m плотность распреде-

пределения случайной величи-

ления падает и при x→ ± ∞ кривая асимптоти-

ны для распределения Гаусса

чески прижимается к оси х.

 

Параметрами распределения являются коэффициенты m и σ. Величина m

здесь соответствует математическому ожиданию случайной величины; σ- среднее квадратичное отклонение величины х.

Если изменить величину m, то кривая распределения будет смещаться вдоль

оси абсцисс, не меняя своей формы. Центр рассеивания характеризует положе- ние распределения на оси абсцисс.

Параметр σ характеризует не положение, а саму форму кривой распреде-

ления. Это характеристика рассеяния. Наибольшая ордината кривой обратно про- порциональна величине σ, чем больше σ, тем ниже максимальное значение f(x). Влияние этих параметров показано на рис.4.4.

f(x)

f(x)

1

 

 

σ1

< σ2

< σ3

 

 

 

 

2

 

 

 

 

3

 

 

 

х

 

 

х

m1

m2 m3

 

 

 

 

 

Рис.4.4. Влияние параметров распределения на вид плотности распределения Гаусса

Во многих задачах, связанных с нормально распределенными случайными

величинами, необходимо определить вероятность попадания случайной величи- ны Х в интервал [α,β]:

P(α < X < β ) = F(β )F(α ) .

Здесь F(x) – функция распределения величины Х:

29

x

1

x

(x m)2

F(x) = òf (x)dx =

òe

 

2σ2 dx .

 

 

 

−∞

s 2p −∞

 

 

(4.10)

Вычисление значения функции распределения требует в каждом случае чис- ленного интегрирования. Чтобы избежать этого, используется нормированное нормальное распределение с m=0 и s=1. Для этого вводится новая переменная u:

u =

x - m

;

 

 

x = m +σu ; dx = σ du .

 

 

 

 

 

 

 

 

 

 

 

 

σ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Функция распределения принимает вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x m

 

 

u 2

 

 

 

 

 

 

u 2

 

F(x) =

 

 

 

1

 

 

 

 

σ

 

 

1

 

u

 

 

 

 

 

 

 

 

2 sdu =

2 du .

 

 

 

 

 

ò

e

 

 

òe

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s

 

2p

 

 

 

−∞

 

 

 

 

 

 

2p −∞

 

 

 

Обычно используется обозначение

 

 

 

 

 

 

 

Ф(u) =

1

 

 

 

u

t 2

 

 

 

 

 

 

 

 

 

 

 

 

 

2 dt .

 

 

 

 

(4.11)

 

 

 

òe

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2p −∞

 

 

 

 

 

 

 

 

 

 

 

 

 

t 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 dt не выражается

 

 

 

 

 

 

 

Интеграл òe

 

через элементарные функции. Это так

называемый интеграл вероятности. Для него составлены таблицы. Поскольку функция Ф(u) симметрична, то Ф(- u) = 1-Ф(u) .

Вероятность нахождения случайной величины X в диапазоне [a, b] определит-

ся выражением

 

 

 

 

 

 

 

 

 

 

 

 

æ

β - m ö

æ

α - m ö

 

 

 

 

 

 

P(α < X < β ) = Фç

 

÷

-Фç

 

÷ .

 

 

 

 

 

(4.12)

σ

σ

 

 

 

 

è

ø

è

ø

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

u

t 2

 

В ряде таблиц приводится значение интеграла

 

 

2 dt для u ³ 0 . При

 

 

òe

 

 

 

 

 

 

 

 

 

 

 

 

2p −∞

 

 

 

u=0 Ф(u)=0,5.

На практике часто встречается задача вычисления вероятности попадания

нормально распределенной случайной величины на участок, симметричный отно-

сительно центра рассеяния m.

Вероятность попадания величины на участок x = m ± σ :

u

β

= m + σ - m

= 1 ;

u = m -σ - m

= -1 ;

 

σ

 

α

σ

 

 

 

 

 

 

30

 

Ф(1)=0,8413;

 

 

Ф(-1)=1-Ф(1)=1-0,8413=0,1587;

 

 

 

P±σ

= P(−σ + m < X < σ + m) = 0,8413 0,1587 = 0,6826 .

 

 

Вероятность того, что случайная величина не попадет на этот интервал

 

 

P±σ

= 1P±σ

= 10,6826 = 0,3174 .

 

 

 

 

 

 

 

Вероятность попадания на интервал m ± 2σ :

 

 

 

 

 

 

P(2) = 0,9772 ,

 

P(2) = 10,9772 = 0,0228 ;

 

 

 

 

 

P±2σ

= 0,9772 0,0228 = 0,9544 ; P±2σ = 10,9544 = 0,0456 .

 

 

Вероятность попадания на интервал m ± 3σ :

 

 

 

 

 

 

P(3) = 0,9986 ,

 

P(3) = 10,9986 = 0,0014 ;

 

 

 

 

 

P±3σ

= 0,9986 0,0014 = 0,9972 ; P±3σ = 10,9972 = 0,0028 .

 

 

 

Таким образом, вероятность того, что случайная величина попадает на интер-

вал m ± 3σ , равна 99,72%. Вероятность того, что она не попадет на этот интер-

вал, равна 0,3%. Отсюда следует так называемое «правило трех сигм», соглас-

но которому с практической достоверностью можно считать, что диапазон изме-

нения случайной величины ограничен пределами m ± 3σ .

 

 

 

 

 

Логарифмически нормальное распределение имеет место, когда нормально

распределена не сама случайная величина, а ее логарифм, т.е. в том случае, ко-

гда случайная величина есть произведение большого числа независимых случай-

ных величин, среди которых нет преобладающих по влиянию на произведение.

f(x)

 

 

 

 

 

Плотность

распределения

 

 

 

 

 

определяется выражением

 

0,04

 

 

 

 

 

 

0,03

 

 

 

 

 

f (x )=

1

 

(ln xm )2

 

0,02

 

 

 

 

 

 

e

2σ

2

 

 

 

 

 

 

2

π

 

 

0,01

 

 

 

 

 

xσ

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

параметр

нормированного

рас-

0

 

 

 

 

 

 

 

 

 

100 x

пределения

 

 

 

 

0

20

40

60

80

 

 

 

 

 

u = ln x m

.

Рис.5. Кривая плотности логарифмически

σ

нормального распределения

 

 

Закон определен для случайной величины x > 0. График плотности распреде-

ления показан на рис. 5. Это распределение позволяет описать распределение случайной величины, имеющую положительную асимметрию.

31

Распределение Вейбулла было предложено в 1939 г. профессором Королев-

ского технологического института Швеции В. Вейбуллом, изучавшим проблему усталости материалов. Функция распределения имеет вид

æ x-m ön

 

÷

 

ç

x0

÷

F(x )=1- e è

ø

Плотность распределения

 

dF (x )

 

 

 

n-1

 

æ x-m ön

 

æ x - m ö

n

 

 

÷

 

x

 

f (x )=

 

ç

 

÷

 

= ç

 

÷

 

 

e è

 

0

ø

dx

x0

 

x0

 

 

è

ø

 

 

 

 

.

Закон определен при x>m, x0>0. Распределение имеет три параметра: m, x0 и n. Параметр m соответствует характеристике положения случайной величины, па-

раметр x0 нормирует значение случайной величины и определяет степень рассе- яния. Параметр n определяет вид кривой плотности распределения и также ока- зывает влияние на рассеяние случайной величины. При 0<n£ 1 график плотности

представляет собой L- образную кривую, при n>0 кривая принимает колоколооб-

разный вид. При m=0 и n=1 распределение Вейбулла переходит в экспоненциаль- ное распределение. Характерный вид кривых показан на рис. 6.

В отличие от нормального распределение Вейбулла не имеет математического обоснования и является только достаточно хорошей аппроксимацией опытных

данных. Вейбулл отмечал, что единственным достоинством этого распределения

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

ше описывает некоторые наблюдения, чем другие известные функции.

Для представления непрерывных случайных величин используются также другие виды распределений, в частности c2 Пирсона, t-распределение Стьюдента, F-распределение Фишера и т.д. Значения функций распределения приводятся в табличном виде в справочной литературе.

32

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]