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

книги из ГПНТБ / Щукин, А. Н. Теория вероятностей и ее применение в инженерно-технических расчетах

.pdf
Скачиваний:
43
Добавлен:
19.10.2023
Размер:
4.5 Mб
Скачать

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

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

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

сложных

поверхностей,

обладающих большим числом

центров

отражения, при

движении этих

поверхностей,

т. е. амплитуды колебаний, отраженных

большинством

радиолокацнон н ых целей.

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

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

По определению (3.3) интеграл плотности вероятно­ сти р(х) в пределах ± оо равен единице, т. е. является конечной величиной. Кроме того, по своему физическому смыслу в большинстве случаев функция р(х) является непрерывной и гладкой. Функции, обладающие подоб­ ными свойствами, как известно, могут быть представ­

лены в

виде интеграла Фурье, т. е. выражены с

по­

мощью

функции ср(/) от новой переменной / таким

об­

разом,

что

 

<?(f) = +J p (x )e - i2*lxdx,

(4.13)

—00

 

+00

 

Р ( х )= J ? (Л е' 2"xfdf.

(4.14)

— 00

 

Если функция р(х) представляет собой сложное ко­ лебание (импульс), где x = t является временем,тофунк-

60

Ция cp(f) является частотным спектром функции р(х), выражающим зависимость плотности амплитуд Iф(/) | от частоты колебаний /. В области вероятностных явле­

ний функция ср(/)

называется

х а р а к т е р и с т и ч е ­

с к о й

для функции

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

плотности вероятно­

сти р(х). При этом

функции р{х)

и ср(/)

связаны друг

с другом зависимостями (4.13) и

(4.14).

Характеристи­

ческая

функция ср(/)

позволяет весьма

просто найти

среднее значение переменной х и среднее значение лю­ бой ее степени. Найдем производную

+00

/2т J хр (х) е~‘2к!х dx,

 

 

—00

 

и положим в ней / = 0. Тогда

 

 

 

+ оо

 

<р' (0) =

— /2т

j хр (х) d x = — /2тх,

ИЛИ

 

—00

 

^ — (//2т) <р' (0).

(4.15)

Точно так же

+00

 

 

 

 

 

= — 4т2

| х-р (х) e~'2r‘!xdx. ср” (0) =

— 4т2* 2,

ИЛИ

—00

 

 

 

 

 

 

£ 2= — ( 1 / 4 п 2) ф " ( 0 ) .

( 4 . 1 6 )

В общем случае

__

 

 

x k=[jk/(2л)'!]Фй(0).

(4.17)

При «свертывании» двух функций раопределения ха­ рактеристические функции, связанные с ними преобра­ зованием Фурье, перемножаются. Это очень важное положение может быть доказано следующим образом. Пусть имеются две непрерывные функции: pi{x) и рг{у). Их характеристическими функциями согласно выраже­ нию (4.13) являются

+ 00

 

 

+ 0 0

<Р, (!) =

f Р, {х) e~‘2r-,xdx,

 

% (/) =

J р2 (у) t~i2llUjdy.

— СО

 

 

— 0 0

Рассмотрим произведение

этих функций

 

+00 +00

 

 

 

< Р , if)

< Р *(f) = 5

J

Р Рг, М)

{Х+У) dxdy-

61

Обозначим

x-\-y z.

Тогда

y-— z x и предыдущая

формула приобретает вид

 

 

 

+ 00 + о о

 

% (f) <Ра (/) =

j

j Pi W

p2 (2 - x) e_/M* dxdz.

 

 

— 00 — 00

 

Перепишем это равенство как

 

+ 00

^

+ о о

?1 (/)

(/) =

е~'2ф dz

j р, (Л') р., (z х) dx (4.18)

 

— 00

 

— 00

и сравним его с формулой (4.13). Очевидно, что в фор-

+ 0 0

муле (4.18) интеграл J р,(х) р.-, (z—х) dx играет ту же

— 00

роль, какую в формуле (4.13) играет р(х), т. е. для это­ го интеграла, представляющего собой согласно (4.11) свертку функций pi(x) и Pz(z), характеристической функцией является произведение cpi(f) фа(/") - Пользуясь формулой (4.13) и преобразованиями, аналогичными вы­ воду формулы (4.18), можно показать, что

Pi М Pi (-*) = j e'2*x!df

[ ?, (ф) <р,(/ — f ) d<l>, (4.19)

— 00

— 00

т. е. характеристической функцией произведения двух функций pi(x) и рг(х) является свертка соответствую­ щих характеристических функций

£?i (Ф) *Ра (f — d<l>.

СО

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

т

(1-14)

1=1

62

Рассмотрим случай, когда средние значения Xi = 0, т. е. .гг2=Ог2, и сначала положим приближенно, что дис­

персии всех слагаемых

ху равны друг

другу:

аг=сг2/»г

Напишем выражение для характеристической

функции

одной из слагаемых случайных величин

(4.13)

 

 

СО

 

 

? i ( f ) =

ГРг(х)еЧ2к1хйх.

(4.13)

 

\)

 

 

СО

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

9г (f) = 9i

(0) + f?' (0) + -g- ?" (0) + -

Ввиду того, что при /= 0

еЧ2ф

+0О

= 1, а ^ р (х ) d x = 1,

— СО

пользуясь формулами (4.15) и (4.16), получаем, прини­

мая во внимание, что т = 0 и х ,-2= сГг2,

 

фг ( / ) = ! + 0 — 2 л Т о Ч т + . . .

(4 .2 0 )

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

р(х). Как было по­

казано раньше (4.11), плотность

вероятности

суммы

двух случайных величин является сверткой их плотно­ стей, а согласно (4.18) характеристической функцией свертки является произведение характеристических функций свертываемых функций.

'В данном случае мы отыскиваем плотность вероят­

ности не двух, а многих

 

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

и потому

характеристическая функция

ф(/)

плотности

вероятно­

сти р{х) 'будет равна

произведению т характеристиче­

ских функций фг(/),. т. е.

 

ф (/) ~ |фг(/) |,п. Если

т очень

велико, то согласно (4.20)

 

 

 

 

 

<? {f) « ( 1 -

2v-f2°2ltn)m=

(4.21)

- и потому на основании

формулы

(4.13)

 

Н-00

 

 

 

j Q- (2 ,^ 4 ^ n df:

р ( х )= ^ ( f ) e i2«x!d f =

—ОО

 

 

 

—00

 

 

- £ + с о ~ ( V 2\ c r - - g r - Y

 

20

Je

 

'

 

V2°>df.

 

 

] e

 

 

 

 

—OP

63

Введем новую переменную а — У 2 mf—jjc/аУ2 ; тогда df = da/m Y 2 ,

+ 00 -

р(л:) = Т^7 Г I е

 

 

 

 

—00

_

Значение

определенного

 

+ со

интеграла j" e~a*da. = У-к . По-

этому

 

 

 

— СО

 

 

 

 

 

 

 

р (х ) =

е- * /а7а 1/2^.

(4.22)

Формула

(4.22)

представляет

собой

н о р м а л ь н о е

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

в е р о я т н о с т е й

или закон рас­

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

(4.22)

мы считали среднее зна­

При выводе формулы

чение каждой слагаемой величины равным нулю (x; = 0). Это ограничение снимается, если при Хгф® ввести новые переменные yi = Xi—5;*, для которых справедливо условие

Уг=0. В таком случае

выражение для

плотности

нор­

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

вероятностей приобретает

вид

р (х - х ) = е~(х-~ху/2°УаУ§г,

(4.23)

Для того чтобы снять ограничение,

заключавшееся

в том, что дисперсии отдельных слагаемых были приня­ ты равными друг другу (ai2= a 2/m), рассмотрим, чему равна плотность распределения вероятности двух слу­ чайных величин, каждая из которых подчиняется нор­ мальному закону.

Характеристической функцией плотности нормально­ го распределения является, как мы видели, выражение (4.21). Поэтому для двух случайных величин с диспер­ сиями оч2 и сг22 характеристическими функциями будут

т,(/) = е - м - '' » ,, (f) =

Плотность вероятности суммы двух случайных вели­ чин равна их интегральной свертке (4.11). Характери­ стическая функция свертки равна произведению харак­ теристических функций свертываемых плотностей (4.18). Поэтому

? „ (/) = ?, (Г)«Р2Ш = е ^ (в,,+в1,,:

64

Характеристической функции cpi2(f) соответствует распределение плотности вероятности

-t»/2 (а.'Ч-оЩ (4.24)

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

Основываясь на этом положении и имея сумму весь­ ма 'большого числа случайных величин Х{, не равных друг другу, можно сгруппировать их таким образом, чтобы все они образовали несколько групп, внутри каж­ дой из которых слагаемые обладали бы приблизительно равными дисперсиями. Если число всех слагаемых Xi очень велико, то их число в каждой группе также будет достаточно большим. Поэтому плотность вероятности суммы случайных величин, входящих в каждую группу, будет подчиняться нормальному закону распределения вероятностей, а распределение плотности вероятности суммы этих групп также будет нормальным.

Последнее допущение, которое 'было сделано при вы­ воде выражения (4.22), заключалось в том, что, разла­ гая функцию q>(f) в ряд, мы ограничились только дву­ мя членами разложения: tp (f)~ (l—2nzfza2/m). Не при­ водя здесь доказательств справедливости такого допу­ щения, укажем, что оно тем более правильно, чем боль­ ше т.

Как показывает выражение (4.22), нормальное рас­ пределение плотности вероятностей представляет собой функцию, симметричную по отношению к прямой х—.г = 0 и простирающуюся от — оо до + оо.

Нормальное распределение имеет максимум при х—.т=0, и потому с его помощью могут быть приближен­ но выражены только те случаи распределения вероятно­ стей, которые имеют один максимум, совпадающий с ра­ венством нулю отклонения случайной величины от ее среднего значения. Свойством нормального распределе­ ния является также возможность неограниченно боль­ ших значений случайной величины х. Это свойство пока­ зывает, что в большинстве случаев, встречающихся на практике, функция нормального распределения лишь приближенно может характеризовать свойства реальных случайных явлений, так как в подавляющем большинст­ ве случаев физически возможные отклонения х ограни­

чены значительно

'бодее узкими пределами, чем ±'оо.

5-147

65

Вид кривых нормального распределения (4.22) для двух значений дисперсии .показан на рис. 4.1. Из выражения (4.22) и кривых на рис. 4.1 видно, что чем больше зна­ чение дисперсии, тем более пологой является кривая р(х).

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

ность нахождения х в пределах между х и x + dx была равна вероятности нахождения у в пределах между у и

y+dy т.

е. р (у)= ор (х). Принимая

это во внимание,

по­

лучаем

 

 

 

 

 

 

Р[У) = ’(1 /

е-г/:/2.

(4.25)

Вид

соответствующей

кривой

.показан

также

на

рис. 4.1, где значения абсцисс даны в долях ст. Интеграл функции (4.25) представляет вероятность того, что у ле­ жит в пределах от — оо до данного значения у

U _!f_

 

 

Р-ЛУ) = Р(у) = у = J e

dy.

(4.26)

—ОО

 

 

Кривая, изображающая p-i(y), показана

на рис. 4.2.

Интеграл

 

 

__ ОО

 

 

^ - ^ < r y'n dy,

 

(3.26)

У

66

которым мы уже пользовались в гл. 3 и значения кото­ рого приведены в табл. III. Приложения, представляет собой вероятность того, что х лежит за пределами ± у среднеквадратических отклонений; при этом очевидно,

что при у

 

 

 

 

 

Р(у) = 1—Ф (У)/2-

(4.27)

Значения р(у),

Р(у) и р'(у)

в зависимости от величины

у даны в табл. V приложения.

 

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

Как уже было сказано,

нормальное

характеризуется

величиной

дисперсии

о2 или средне­

квадратическим

отклонением а.

Иногда

применяется

также понятие

меры точности h,

связанной со средне-

квадратически'М отклонением

Л = 1 /о | / 2 ".

Кроме дисперсии и среднеквадратического или стан­ дартного отклонения, нормальное распределение можно также характеризовать с р е д н и м о т к л о н е н и е м хСр, которое равно

+“

---

Хср=77Ш

KlJ>dx— V

о

 

или хСр= 0,798а.

Наконец, часто употребляется термин вероятное или срединное отклонение, под которым подразумевается значение хр, при котором вероятность того, что х равно или меньше по абсолютной величине хр, равна 0,5, т. е.

+ х р

Р ( И < х р) = - ^ = j e ^ d x =

=

] / - ~ j

е dy — 0,5.

 

U

 

Из та'бл. III приложения

видно, что при -Р(|х|^Хр) =

= 1— Ф (у )~ 0,5

у = 0,6745,

т. е. хр=0,6745сг.

5*

67

Функции р(х)

(4.22)

и р(у)

(4.25)

неограниченно уво­

дят в обе стороны

от прямых х = 0 и г/= 0, Однако, как

можно видеть из табл. III, менее 5%

площади этих кри­

вых лежит за пределами

|х| > 2 о или

|г/| !>2, т. е. веро­

ятность

того, что

случайные

значения х или у будут

больше

по абсолютной

величине,

чем соответствено 2а

и 2, не превышает 0,05.

условно

считают максимально

На

практике часто

возможное отклонение равным утроенному среднек'вадратическому отклонению 3 }/гх2=За или учетверенному срединному отклонению 4хр.

Огибающие дискретных распределений: биномиаль­ ного (2.2) и Пуассона (4.'2) — приближаются к кривой нормального распределения плотности вероятности при больших значениях величин т и п в биномиальном рас­ пределении и величин -V п а: в распределении вероятно­ стей Пуассона. Для того чтобы сопоставить кривую нор­

мального распределения р{у) = ( 1/ 1^ 2 тг) е~и’!% с ди­ скретным распределением Пуассона Р(х) ='(хх/х\) е~х, последнюю следует провести к масштабу нормального распределения. Дисперсия распределения Пуассона мо­ жет быть получена из дисперсии биномиального распре­

деления

Об2= т р ( \ — р) при

тр = х\ т >-оо;

рИЗ, т. е.

х. ,

Поэтому изменению х на единицу в распределе­

нии Пуассона соответствует

изменение у на

1/з^ = 1/]/л:

в нормальном распределении.

Масштаб по оси ординат в распределении Пуассона

приводится

к

масштабу нормального

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

умножением на

V x так как ординате

р (0) =

1 /'|/’2тс в

нормальном

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

соответствует

ордината

Р (х— х) =

1/У'2'кх при

1 в распределении Пуассона.

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

величины х, ниже приводится

несколько примеров.

1.

ж=100; ж—*i = 10;

х х г = — 10; *i = 90;

*2= 110. Как у

было сказано выше, величина

(/= 1,0

в

нормальном

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

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

Пуассона

ап =У.т, т. е.

в рассматри­

ваемом случае оп =]/.г= 10. Поэтому

абсциссе х jci=10

соответст­

вует (/=1,0.

Принимая во внимание

коэффициент У*

при

переходе

от масштаба по оси ординат распределения Пуассона к нормально­ му, ордината, соответствующая *t= x — 10=90, может быть найдена

68

Мз Выражения

(4.3)

 

 

 

 

 

 

 

 

lg [Я (х,) K *j =

9 0 ( l g l0 0 - i g 90) —

 

 

 

0,434294-10 — 0,39908 —

(Ig 100 — lg 90) =

F,39925,

т.

е.

Я (х ,) К г =

0,25078.

 

 

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

соответствует

Р

Значению

(/=1,0

в нормальном

(у= 1,0) =0,24197.

Таким образом, ордината огибающей распреде­

ления

Пуассона

в

рассматриваемом

случае при xi=90 на +3,6%

больше соответствующем ординаты нормального распределения.

 

При .т=100 и х2=1Ю.

 

 

 

 

 

 

lg

 

(*,) Кxj =

110 (lg

ПО — lg 100) +

 

+ 0,434294,10 — 0,39908 — 0,5 (lg 110— lg 100)= 1,37023;

Я (х 2) Кж = 0,23456.

Сравнение с той же величиной Р ( у — 1) = 0,24197 показывает, что

Р (х2) ]Гх на —3,02% меньше ординаты нормального распределения. 2. .7=100; х—xi=20; 7—хг= —20; Xi=80; Х2= 420. В этом случае Я (х ,) Кг = 5,2035-ю - 2; Я ((/ = 2,0)=5,3990- 10_ 2,т. е. при х, =80

ординаты огибающей распределения Пуассона на —3,6% меньше

нормального распределения, а при х2=1'20,

Р(хг)1/г = 5,567 •10-2, т. е.

ордината

огибающей

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

Пуассона при Хг=120 на

+3,1% больше ординаты нормального

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

 

3. Расчеты, выполненные для случаев:

г =

10*;

х, = 9900; х2=

=

10 100, дают Р (х,) Кг =

0,24245

вместо 0,24197 по нормальному

распределению,

т.

е.

на

+ 0, 2%

 

больше,

а при х 2= 10 100,

Я (х2) Кг =

0,24125,

т. е. на —0,3%

меньше,

чем

по нормальному'

распределению.

 

 

 

 

 

 

 

 

 

5,39610~2, т. е.

на

4. Для 7 = 404; xi=9800; х2=Ю 200; Я(х,)У7 =

—0,06%

меньше,

чем

при

 

нормальном

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

а Р (х2) К.7=5,418-10-2, т. е. на +0,35%

больше,

чем по нормаль­

ному распределению.

увеличение

7

в

100

раз

(со

100 до 10 000)

 

Таким

образом,

уменьшает расхождение между ординатами огибающей распределе­ ния Пуассона и соответствующими ординатами нормального распре­ деления примерно на порядок, при величине расхождений на 7 = IО4 порядка нескольких десятых долей процента. При дальнейшем уве­ личении 7 расхождения становится еще меньше и огибающая рас­ пределения Пуассона практически совпадает с нормальным распре­ делением.

69

Соседние файлы в папке книги из ГПНТБ