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

Методы оптимизации эксперимента в химической технологии

..pdf
Скачиваний:
39
Добавлен:
15.11.2022
Размер:
13.27 Mб
Скачать

В виде функции распределения можно задать распределение как

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

Как видно

из определения, F(x) есть неубывающая функция х; если

х ^ х 2, то

F(xy)^.F(x^ (рис. 2, а). Ордината этой кривой, соответствующая точке

хь

представляет собой вероятность того,

что

случайная

величина

X

при испытании

окажется

< х 1.

Разность

двух ординат,

соответ­

ствующая точкам

х} и х2,

дает

вероятность

того, что

значения

случайной величины будут лежать в интервале между х, и х2:

 

 

 

Р (Xj < X < х2) = F (х2) F (X!).

 

(1.12)

Значения функции распределения при предельных значениях аргумента соответственно равны 0 и 1:

F( — оо) =

О, F (-(- оо) = 1.

(1.13)

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

всегда

есть разрывная ступенчатая

функция, скачки которой происходят

в точках, соответствующих возможным значениям случайной величины,

и

равны вероятностям этих

значений (рис. 2, б). Сумма всех скач­

ков равна 1.

 

 

 

ся

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

производная функции

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

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

случайной величины X. Если F(x)

непрерывна и дифференцируема, то

 

 

/<*) =

а д

(1.14)

 

 

dx

Задание f(x) тоже полностью определяет случайную величину. Плотность распределения является неотрицательной функцией (рис. 3). Площадь, ограниченная осью х, прямыми х = х 1 и х=х2 и кривой плотности распределения, равна вероятности того, что случайная величина примет значения из интервала х, н-х2:

х*

 

 

Р (*! < X « х») = j f(x)Ax=F(xt) F(x1),

(1.15)

Х\

 

 

в частности

X

 

F (х) = Р (—оо < X < х) =

j f(x) dx.

(1.16)

— ОО

Отсюда же выводится еще одно важное свойство плотности рас­ пределения:

 

j f (x) d * = l ,

(1.17)

 

 

— ОО

 

 

гак как попадание случайной величины

 

в интервал -°°< Х< + 00 есть достовер­

 

ное событие.

 

 

2.

Числовые характеристики. Вместо,

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

Рис. 3. Плотность распределения не­

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

ны в виде законов распределения вероятностей в прикладных за­ дачах ее часто определяют при помощи числовых характеристик — чисел (вещественных), выражающих характерные особенности случай­ ной величины, называемых моментами случайной величины. Для дис­ кретной случайной величины начальный момент /t-го порядка опре­ деляется формулой

* = , . 2 .......

<1Л8)

1=1

 

для непрерывной случайной величины —формулой

 

оо

 

тц= J ** / (*) dJt.

(1.19)

— ОО

Начальный момент первого порядка (7t=l) называется математиче­ ским ожиданием (средним значением) случайной величины. Математи­ ческое ожидание принято обозначать различным образом:

М[Х], тх, т.

Для дискретных случайных величин

п

(1-20)

т1 — М (X) = 2 *lPi-

t=i

 

Для непрерывных случайных величин математическое ожидание выражается интегралом:

оо

 

тх = М[Х]= ^xf(x)dx.

(1.21)

— оо

 

Чаще, чем начальные моменты, применяются центральные момен­ ты. Центральный момент /t-го порядка для дискретной случайной величины определяется формулой

 

п

 

t*A=

mx)kPl,

(1.22)

 

1= 1

 

для непрерывной случайной величины —формулой

 

оо

 

НА = |

( x - m x)*f(x)dx.

(1.23)

— оо

Первый центральный момент всегда равен 0, pi=0. Второй централь­ ный момент называется дисперсией. Дисперсией случайной величины называется математическое ожидание квадрата отклонения случайной величины от ее математического ожидания, т. е.

D[X) = M[ ( X - m x)2].

(1.24)

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

 

п

 

D [X] = ц2= 2 (*1 - тх)г Pi,

(I • 25)

1= 1

для непрерывной

оо

(1.26)

D[X] = j (х —тх)3 f(x)Ax.

Другие обозначения для дисперсии Dx, сх2, а 2. Корень квадратный из второго центрального момента называется средним квадратичным отклонением (или стандартом):

о, = V D [X]= V 77 .

(1.27)

Третий центральный момент, разделенный на а*3, называется коэф­

фициентом асимметрии:

, о

(1.28)

 

Ti = М-з/< •

Через начальные моменты

Цз выражается следующим образом:

 

[Х3 = т 3 — 3/71^2 + 2т\ .

(1.29)

Четвертый центральный момент вычисляется по формуле

 

(х4 = т 4 — Агпггпъ+ 6т\ т2 — Ът\ .

( 1 .30)

Величина

Tf2 = (р4/а*) — 3

(1.31)

 

называется коэффициентом эксцесса.

На рис. 4 приведены примеры плотностей распределений с не­ нулевыми коэффициентами асимметрии и эксцесса. Для сравнения штриховой линией изображена кривая с тем же математическим ожиданием тх и дисперсией ст*, но с нулевыми значениями коэф­ фициентов эксцесса и асимметрии.

Моменты существуют, если соответствующие интегралы или ряды для дискретных величин сходятся. Для случайных величин, значе­ ния которых ограничены, моменты всегда существуют. Если у

случайной

величины

X

суще­

 

 

 

 

 

ствуют первый и второй момен­

 

 

 

 

 

ты, то можно построить норми­

fix)

 

1

 

 

рованную случайную eejimuny:

 

 

 

У

Т

____

р

\ * ' <0

 

 

 

 

 

Х0= ( Х - т х)/ах.

 

(1.32)

 

и

 

 

 

 

 

 

 

 

 

Для нормированной случайной

0

X

0

 

X

величины

 

 

 

 

 

 

 

 

М[Х0]= 0,

D [Х0] =

1.

(1.33)

fix)

 

 

 

 

 

Г С * * "

" "

 

 

 

 

 

 

 

 

 

 

Многие

таблицы

распреде­

•у

\ * г"

 

 

 

лений

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

 

 

 

 

 

 

 

 

нормированных случайных ве­

0

X

0

 

X

личин.

Существуют

следую­

 

 

 

 

 

щие соотношения между функ­

РИС; 4.

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

циями распределения, соответ-

выми коэффициентами асимметрии и эксцесса

ствующими нормированной величине Хо и ненормированной X:

 

 

(1.34)

fl (*о) = °xf (х) = ax f

+ ах *о),

(1.35)

F(x) = F 1(x0) = Fi y - -

х J,

(1.36)

Pi (*о) = F (х) = F (тх + ах х0).

(1.37)

Моменты являются общими (интегральными) характеристиками распре­ деления. Вторая группа параметров характеризует отдельные значения функции распределения. К ним относятся квантили. Квантилем хр распределения случайной величины X с функцией распределения F(x) называется решение уравнения

F(xp) = р,

(1.38)

т. е. хр есть такое значение случайной величины, что

 

Р ( Х с хр) = р.

(1.39)

Если известны два квантиля хр и xq, то

 

Р(хр < Х с х д) = q —p.

(1.40)

Наиболее важное значение имеет квантиль х% , называемый медианой распределения (рис. 5). Ордината медианы рассекает площадь между кривой плотности вероятности и осью абсцисс пополам. Если распре­ деление симметрично,

Xdt = mx

(1.41)

Квантили хр и хх_р называются симметричными. Для симметричного относительно нуля распределения всегда

хр = — *!-/>•

(1.42)

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

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

идисперсии случайных величин;

1.Математическое ожидание неслу­ чайной величины равно значению этой величины;

 

М[с] = с.

(1.43)

 

2. Неслучайную величину можно вы­

 

носить за знак математического ожидания:

Рис. 5. Медиана распределения

М[сХ] = сМ[Х].

(1.44)

14

3. Математическое ожидание суммы случайных величин равно сумме математических ожиданий этих случайных величин:

М [ Х г + Х 2 + . . . + Х Я] = М [ Х 1] + М 1 Х %) + . . . + М [ Х Я].

(1.45)

4. Математическое ожидание произведения независимых случайных величин равно произведению математических.ожиданий сомножителей:

М [ Х г- Х% . . . Х п] = М [ Х 1] М [ Х 2'\ . . . М [ Х п].

(1.46)

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

Свойства дисперсии:

1. Дисперсия неслучайной величины равна нулю:

D [с] = 0.

(1.47)

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

D [сХ] = c2D [ X ] .

(1.48)

3. Дисперсия случайной величины равна математическому ожида­ нию квадрата случайной величины минус квадрат ее математического ожидания:

D [Х\ = М [ Х * \ - т 2х .

(1.49)

4. Дисперсия суммы независимых случайных величин равна сумме дисперсий этих величин:

0 [ Х г + Х 2 + . .. + X n] = D [ X 1] + D [ X 2] + ... + D [ X n).

(1.50)

Используя свойства математического ожидания и дисперсии, покажем, что для нормированной случайной величины справедливо утверждение (L33), т. е. если Х0=(Х — тх)/ ах, то

 

М[Х0) = 0,

D [Х0] =

1,

М 1Хв1 = М | Х а

т* j =

(X — т,) =

|л1 (X) — тх 1 = 0,

( ^

) ^

0 ( Х _

^

к н = 0 и и , .

В | Х Л - 0 ' ^ = х

 

 

 

 

Пример 2. В результате испытаний двух расходомеров установлена вероятность наблюдения помех, оцениваемых по двухбалльной системе:

Уровень помех,

Вероятность

наблюдения

помех данного уровня

балл

 

 

 

 

расходомер

1

расходомер 2

1

0,20

 

0,03

2

0,065

 

0,15

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

Р е ш е н и е . Обозначим через X случайный уровень помех расходомера. Определим средний уровень помех для каждого расходомера по формуле (1.20):

Ml [XJ = 0,20.1 +0,065-2 = 0,33,

М2[Х[ = 0,03-1 +0,15-2 = 0,33.

Таким образом, средний уровень помех у обоих расходомеров одинаков и по этому показателю нельзя выбрать лучший прибор. Определим устойчивость показаний, для этого по формуле (1.25) посчитаем дисперсии уровня помех для каждого расходомера:

Dx [X] = (1 -0,33)2 0,2 + (2— 0,33)*.0,065 = 0,11;

D2[X] = ( I— 0,33)2-0,03 + (2 — 0,33)2-0,15 = 0,43.

Следовательно, лучшим является первый расходомер.

Рис. 6. Плотность вероятности

Рис. 7. График функции F(x) рав-

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

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

4. Равномерное распределение. Определим основные числовые характеристики одного из простейших непрерывных распределений — равномерного распределения. Равномерным распределением называется распределение, для которого плотность вероятности постоянна в определенных пределах и равна нулю вне этих пределов (рис. 6). Плотность f(x) постоянна и равна с на отрезке [а, Ь\ вне этого отрезка она равна нулю:

с при а < х<Ь,

0 при х < а или х>Ь.

Так как площадь, ограниченная кривой распределения, равна едини­ це: c(b - а) = 1, то с = \/(Ь - а), и плотность распределения f(x) имеет вид

I

/ (*) = 7----- при а < х < Ь,

оа

(1.51)

f (х) = 0 при х с а или х > Ь.

Функция распределения выражается площадью кривой распределе­ ния, лежащей левее точки х. Следовательно,

 

0 при х < а ,

 

 

F(x) =

х — а

х < Ь,

(1.52)

------ при а <

 

Ьа

 

 

I при х > Ь.

График функции F(x) приведен на рис. 7. ‘Математическое ожидание случайной величины X, имеющей равномерное распределение на отрезке [а, b], равно

С х f

а+Ь

 

,1' -

— -

<, ю )

В силу симметричности равномерного распределения медиана вели­ чины X также равна + Ь)/2.

По

формуле

(126)

определим

дисперсию случайной величины X:

 

 

 

 

 

ь

 

Ф — д)*

 

 

 

 

 

 

 

 

(1.54)

 

 

 

 

 

 

 

12

 

 

 

 

 

а

 

 

 

 

 

 

 

 

 

 

Коэффициент

асимметрии Yie H3/ax

равен нулю (распределение

симметрично).

 

 

 

 

 

 

Для определения коэффициента эксцесса найдем четвертый

центральный момент:

 

 

 

 

 

 

 

 

 

 

 

— а)4

 

 

 

 

 

 

 

 

80

 

отсюда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

- 1, 2 -

 

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

величины

X

на

отрезок [ a , р ],

представляющий

собой часть от­

резка

[а,

Ь]

(рис. 8),

определяется

отношением

длины отрезка

[ a, Р ] к длине всего отрезка [а, Ь]:

 

 

 

 

 

 

 

 

Р( а<Х<$ ) =

 

(1.55)

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

—т )*

f ( x )

где тх и —математическое ожи­ дание и дисперсия случайной вели­ чины X.

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

 

 

 

( х - т ) *

1

с

е

х dx (1.57)

f w = • ■-__—

\

У2" «*

-t

 

 

оо<дг<оо), (1.56)

Рис. 8. Определение вероятности попадания равномерно распреде­ ленной случайной величины на заданный участок

Нормальное распределение наиболее часто встречается на практике и теоретически наиболее полно разработано. Нормальный закон при некоторых условиях является предельным законом для суммы боль­ шого числа п независимых случайных величин, каждая из которых подчинена какому угодно закону распределения (теорема Ляпунова). Основное ограничение состоит в том, чтобы все слагаемые играли в общей сумме относительно малую роль. Множество событий происходит случайно вследствие воздействия на них большого числа независимых (или слабо зависимых) возмущений. У таких явлений закон распределения близок к нормальному. Нормальный закон распределения широко используется при обработке наблюдений. Пользуясь методами теории информации, можно показать, что нор­ мальное распределение содержит минимум информации о случайной величине по сравнению с любыми распределениями с той же дисперсией. Следовательно, замена некоторого распределения экви­ валентным нормальным не может привести к переоценке точности наблюдений. График плотности нормального распределения называет­ ся нормальной кривой или кривой Гаусса (рис. 9).

Рис. 10. График функции. Fo(x) стандартного нормального рас­ пределения

Нормальное распределение нормированной случайной величины называется стандартным. Его функция распределения имеет вид

X

(I.5S)

График этой функции представлен на рис. 1Q Для такой величины

Р (*! < Х0<

х2) = F0 (х2) — F0(*i).

(1.59)

Функция

 

 

Ф (X) = Г, (*) -

(1.60)

называется функцией Лапласа:

X

 

 

 

Ф (х) = Г0 (*) -

Г0 (0) = — != f е~х'п Ах.

(1.61)

 

У2тГ Jj

 

Значения этой функции приводятся в табл. 1 приложения. Функция Лапласа —нечетная функция, т. е.

Ф(— *) = — Ф(*),

(1.62)

поэтому таблицы значений Ф(х) составлены лишь для х>0.

Для нормированной случайной величины, учитывая (1.59) и (1.60), имеем

Р (*i < Х0< х2) = Ф (х2) + V* — Ф (хг) — V2 = Ф (х2)— Ф (*,). (1.63) В общем случае

Р(х1< Х < х 2) = Р (— ^ тх < Х В< — ~ т* } =

Во многих практических задачах х, и х2

симметричны относитель­

но математического ожидания, в частности

в задаче об абсолютном

отклонении. Абсолютным отклонением называется величина

 

АХ = |Х — тх\.

(1.65)

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

Р (АХ < е) = Р (тх —с < х < т л. + е).

(1.66)

В частности, для нормированной случайной величины

 

Р (ДХ0 < е) = Я ( - Е< Х0 < + «) = Ф(0 - Ф ( - «) = 2Ф (г).

(1.67)

Для случайной величины, имеющей нормальное распределение с параметрами тх И GX ,

Р (ДХ < 0 = Р^ДХ0 < J - J = 2Ф

.

(1.68)

Обозначив е/вх = К получим из (1.68) следующее соотношение:

 

Р (АХ < k?x) = 2Ф (k),

 

(1.69)

отсюда

 

 

Р (АХ < ах) = 2Ф (1) = 0,6826,

 

 

Р (АХ < 2ах) = 2Ф (2) = 0,9544,

 

(1.70)

Р (АХ < Зох) = 2Ф (3) = 0,9973.

 

 

Таким образом, отклонения больше чем утроенный стандарт (средне­ квадратическое отклонение) практически невозможны. Нормальное распределение обладает свойством линейности: если независимые случайные величины Х} и Х2 имеют нормальные распределения, то для произвольных чисел а и Р величина

К = аХ1 + рХ2

(1.71)

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

ту =

amXl + $mXt

 

(1.72)

и

 

 

 

 

Яу — | /

а2

+ Р*°х,

'

0*73)

Пример 3. Толщину керамической глазурованной плитки И можно считать нормально

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

со

стандартом

а = 0,3

мм. Какова вероятность

брака, если бракуются плитки, толщина которых отклоняется от номинала (математи­ ческого ожидания) более чем на 0,5 мм.

Р е ш е н и е . Определение вероятности брака сводится к решению задачи об абсо­ лютном отклонении для случайной величины И—толщины плитки. Необходимо опре­ делить Р(АИ** 0,5). Найдем вероятность противоположного события по формуле (1.68):

Р h < 0,5) = 2Ф (0,5/0,3) = 2Ф (1,67) = 2 0,4525 = 0,905,

откуда

Р(Д/t^0,5) = 1-0,905 = 0,095.

6.Системы случайных величин. В практических приложениях теории вероятностей и математической статистики очень часто при­ ходится иметь дело с задачами, в которых результат эксперимента описывается не одной случайной величиной, а двумя или более случайными величинами, образующими систему.

Свойства системы нескольких случайных величин не исчерпывают­ ся свойствами отдельных величин ее составляющих: помимо этого, они определяются также взаимными связями (зависимостями) между случайными величинами. Информация о случайных величинах содер­ жится в законах распределения. Рассмотрим систему из двух случай­ ных величин (Xf У).

Функцией распределения системы двух случайных величин (X, У)

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

F(x, у) = Р ( Х с х у Ус у ) .

(1.74)

В § 1 приведены основные свойства функции распределения

F(x)

для одной случайной величины. Сформулируем аналогичные свойства

для функции распределения системы (X, У).

неубывающая

функция

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

F(x, у) есть

обоих аргументов, т. е.

 

 

 

при *2 > хг

F (х2, у) > F (х1у

у);

(1.75)

при у2 >Ух

F (ху y2) > F ( x y ух).

 

2. Всюду на —00 функция распределения равна нулю:

 

F (х, оо) = F (—ос, у) = F (— оо, —оо) = 0.

(1.76)

3. При одном из аргументов, равном +°°, функция распределения системы превращается в функцию распределения случайной величины, соответствующей другому аргументу:

F(x, + oo) = F1(x),

(1.77)

F(+oot y) = F2(y)t

 

где Fy(x), F2(y)~ соответственно функции распределения случайных величин X и У.