
книги из ГПНТБ / Щукин, А. Н. Теория вероятностей и ее применение в инженерно-технических расчетах
.pdfсама суммой Весьма большого числа случайных вели чин, каждая из которых вносит малую долю в оконча тельный результат.
Нормальный закон распределения дает возможность совершенно строго выразить такие величины, как мгно венные значения электрического тока, обусловленные его дискретной структурой из отдельных электронов, и ши роко используется при анализе шумов радиоэлектронной аппаратуры.
Нормальному закону подчиняются мгновенные значе ния величины электрических колебаний, отраженных от
сложных |
поверхностей, |
обладающих большим числом |
|
центров |
отражения, при |
движении этих |
поверхностей, |
т. е. амплитуды колебаний, отраженных |
большинством |
радиолокацнон н ых целей.
Нормальным законом распределения часто пользуют ся для определения случайных отклонений снарядов от центра рассеяния, а также отклонений многих других величин от их средних значений. Следует, однако, отме тить, что здесь нормальный закон является только при ближением, более или менее справедливым лишь в тех случаях, когда среди множества причин, вызывающих отклонения, нет доминирующих.
Прежде чем переходить к формулированию закона нормального распределения, остановимся предваритель но на некоторых общих свойствах функций распределе ния плотности вероятности.
По определению (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