книги из ГПНТБ / Клемин А.И. Инженерные вероятностные расчеты при проектировании ядерных реакторов
.pdfД и ф ф е р е н ц и а л ь н ы м |
з а к о. и о м распределения |
не |
|||
прерывной случайной величины называется функция* |
|
||||
/ (л-) = |
dF (x)/dx, |
f |
(л-) > 0. |
(2.3) |
|
Очевидно, что |
|
|
|
|
|
j" |
/ (x)dx |
= |
F |
(x). |
(2.4) |
•—CO |
|
|
|
|
|
Функция f(x) численно равна вероятности попадания случайной
величины X в достаточно |
малый |
интервал единичной' длины, по |
||||
строенный около точки х(х |
< |
X < |
х + А, Д = 1). Поэтому диффе |
|||
ренциальный |
закон f(x) часто |
называют п л о т н о с т ь ю |
р а с |
|||
п р е д е л е н и я |
в е р о я т н о с т е й . |
|
||||
Чтобы найти дифференциальный закон f(x) опытным путем, на |
||||||
до произвести |
я |
наблюдений |
над |
случайной величиной X. |
Если |
|
результаты наблюдений нанести на числовую ось и пронумеровать по порядку xlt Хо, х„, а интервал [ху, х„] разбить на N= 5ч-12 не обязательно одинаковых подынтервалов шириной Ах, то по оп
ределению |
эмпирическая |
плотность вероятности |
|
|||||||
|
|
/Э |
(Л-) = |
Ал. (х)/пАх, |
|
|
(2.5) |
|||
где х — середина |
подынтервала; |
Ал |
(х) — число наблюдений, по- |
|||||||
павших в |
данный |
подынтервал |
|
X |
Дд; |
, X |
, |
ДдЛ |
. Обычно п и N |
|
|
2 |
+ |
2 |
|||||||
выбирают так, чтобы 5 ^ |
Ал(лг) < |
40. |
|
|
|
|||||
Из формулы (2.3) вытекает следующее свойство плотности рас |
||||||||||
пределения |
вероятностей: |
|
|
|
|
|
|
|
|
|
|
|
]f(x)dx |
= |
F(oo) = |
l. |
|
(2.6) |
|||
Для дискретной случайной величины X дифференциальный за кон распределения вырождается в ряд распределения: ряд чисел, равных вероятностям, что X примет конкретное значение xt:
Pi = Р(Х = ХІ] для всех і = 1,2, |
п, (2.7) |
где п — полное число возможных значений дискретной случайной величины (в случае счетного множества значений п = со ). Свой-
|
|
|
|
|
|
|
|
|
|
• п |
|
|
|
|
|
|
ство |
(2.6) в этом |
случае |
принимает |
вид |
^ Р г = |
1- |
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
і = 1 |
|
|
|
|
|
|
* З д е с ь |
п р е д п о л а г а е т с я , |
что |
ф у н к ц и я |
F |
(х) д и ф ф е р е н ц и р у е м а |
во |
всей |
|||||||||
о б л а с т и в о з м о ж н ы х |
з н а ч е н и й |
X; о б ы ч н о |
это |
в ы п о л н я е т с я |
в б о л ь ш и н с т в е |
|||||||||||
п р а к т и ч е с к и х |
з а д а ч . |
Е с л и F {х) не |
и м е е т |
п р о и з в о д н ы х в о т д е л ь н ы х |
т о ч к а х , то |
|||||||||||
э т о о г р а н и ч е н и е |
н е п р и н ц и п и а л ь н о , |
и б о д а ж е |
п р и |
нем |
п р а к т и ч е с к и |
в с е г д а |
||||||||||
у д а е т с я п о с т р о и т ь / |
(х), д и ф ф е р е н ц и р у я |
F |
(х) |
по |
о б л а с т я м , |
где |
с у щ е с т в у е т |
|||||||||
п р о и з в о д н а я , |
и |
и с п о л ь з у я н о р м и р о в к у |
(2.6). В |
п о с л е д н е м с л у ч а е |
/ (*) |
б у д е т |
||||||||||
и м е т ь |
р о в н о |
с т о л ь к о |
т о ч е к р а з р ы в а , |
в |
с к о л ь к и х т о ч к а х |
о т с у т с т в у е т |
|
dF/dx. |
||||||||
В е р о я т н о с т ь п о п а д а н и я с л у ч а й н о й |
в е л и |
||
ч и н ы в |
з а д а н н ы й |
и н т е р в а л а ^ Х < |
р находим |
по формуле |
|
|
|
Р |
{а < X < р} = |
F (Р) - F (а) = ) f (х) dx. |
(2.8) |
|
|
а |
|
Из нее вытекает интересное следствие для непрерывных случай ных величин. При а = р Р{Х = а) = 0, т. е. вероятность, что непрерывная случайная величина в результате опыта примет точ но заданное значение, равна нулю.
§ 2.2. Определение основных числовых характеристик случайной величины
Математическое ожидание. Среднее значение случайной вели чины X, или ее м а т е м а т и ч е с к о е о ж и д а н и е М(х), определяют следующим образом:
для дискретных величин
М (х) = 2 x l P u |
(2.9) |
дли непрерывных случайных величин
оо |
|
|
М (х) =— сJо xf |
(х) dx, |
(2.10) |
со |
|
|
если существует интеграл —j со |
dx. |
|
Из формул (2.9) и (2.10) вытекают правила для вычисления математических ожиданий случайных величин:
если с — константа, то
М(с)=с; |
(2.11) |
М (сХ) = сМ (X); |
(2.12) |
М (X + У) = М (X) + М (У); |
(2.13) |
М {X-Y) = М (X) М (У) + cov (ХУ), |
(2.14) |
где ковариация СОУ(ХУ) определяется ниже, см. выражение (2.25). Из последнего правила следует, что математическое ожидание про изведения независимых случайных величин равно произведению их математических ожиданий:
/ т \ |
т |
М[Пхк)= |
ПМ (Xh). |
(2.15) |
Медиана и мода. М е д и а и а |
численно равна |
такому значе |
||||
нию случайной величины Me(yY), для которого |
|
|
||||
|
Р {X > Me (X)} = Р {X < Me (X)} = |
1/2. |
|
|||
М о д а — наиболее вероятное |
значение случайной |
величины, |
||||
при |
котором |
функция дифференциального закона имеет макси |
||||
мум: /[Мо(Л/)1 = max, если X непрерывна, или Р{Х = |
Мо(Х)} = |
|||||
= max, если |
А' дискретна. |
|
|
|
||
Дисперсия |
и среднее |
квадратическое отклонение. |
Д и с п е р |
|||
с и я |
— мера |
разброса |
значений |
случайной величины около ее |
||
математического ожидания. Численно она равна математическому
ожиданию квадрата |
разности между самой |
случайной |
величиной |
и ее математическим |
ожиданием: |
|
|
D (X) = М [X — М (Х)Г- = М (X2) |
— 1М (Х)1\ |
(2.16) |
|
На практике для характеристики разброса значений случайной величины около ее математического ожидания используют также корень из дисперсии
|
|
|
a = VD, |
|
|
(2.17) |
||
называемый |
с р е д н и м |
к в а д р а т и ч е с к и м |
о т к л о |
|||||
н е н и е м. |
|
|
|
|
|
|
|
|
На основе |
формул (2.11) — (2.13) и |
(2.16) легко получить сле |
||||||
дующие правила для вычисления дисперсий: |
|
|
||||||
если с — константа, то |
|
|
|
|
|
|||
|
|
|
D (с) = |
0; |
|
|
(2.18) |
|
|
|
D (сХ) = c2D (X); |
|
(2.19) |
||||
D |
(X + |
Y) = D (X) |
+ D (Y) |
+ 2 cov (XY). |
|
(2.20) |
||
Из последнего правила следует, что для независимых |
случайных |
|||||||
величин |
|
( |
|
|
|
|
|
|
|
|
т |
\ |
т |
|
|
|
|
|
|
|
2 * J = |
2 |
|
|
(2.21) |
|
|
|
k=\ |
I |
k=\ |
|
|
|
|
Асимметрия и эксцесс. Эти две числовые характеристики |
зако |
|||||||
нов распределения |
случайных |
величин |
обычно используются |
для |
||||
оценки отличия кривой распределения /(х) от так называемого нормального закона (подробно о нем см. § 3.1). Как асимметрия, так и эксцесс являются безразмерными величинами, определяе
мыми соответственно |
по формулам: |
А = |
\М[Х-М{Х)\*^^; |
|
(2.22) |
Е = — М [Х—М (X)]"— 3 = -Ні- —3.
где ц т |
— центральный момент m-го порядка |
случайной величины |
X. Для |
всех симметричных распределений |
А = 0. Если А > 0 |
(положительная асимметрия), то длинная ветвь плотности рас
пределения |
лежит |
правее М(Х), при А < |
0 — левее. Эксцесс |
равен нулю |
только |
для нормального закона. |
Эксцесс Е > 0 по |
казывает, что кривая распределения имеет более острую и высо кую вершину, чем нормальный закон. Отрицательный эксцесс ука зывает на более плосковершинную кривую.
§ 2.3. Определение числовых характеристик корреляции и регрессии
Стохастическая связь. Между величинами Хн Y возможна связь двоякого рода: 1) функциональная, когда конкретному значению X = х соответствует одно единственное значение Y = у (много значные функции не будем рассматривать); 2) стохастическая или вероятностная, когда одному значению X = х соответствует мно жество возможных значений Y, распределенных по определенному
Р и с . 1. Ф у н к ц и о н а л ь |
н а я |
(1) |
и с т о х а с т и ч е с к а я (2) с в я з и д в у х |
|
в е л и ч и н X и |
Y |
(АВ |
—• л и н и я |
р е г р е с с и и ) . |
вероятностному закону f{ylx), |
зависящему |
от этого х (рис. 1). Эта |
||
связь проявляется в том, что при изменении х меняется закон рас
пределения |
f{y/x). |
|
|
|
Стохастическую связь |
характеризуют |
ф у н к ц и е й |
р е г |
|
р е с с и и |
Y по X: |
|
|
|
|
~~у{х) = |
М (Y/x) =]yf |
(у/х) dy. |
(2.23) |
Это неслучайная функция, в каждой точке х равная условному математическому ожиданию случайной величины У (на рис. 1 — кривая А В).
Большое значение для практики имеет случай (он встречается чаще других), когда функция регрессии линейна. Кстати, такой она будет, если случайные величины X и Y распределены по нор
мальному закону (закону |
Гаусса, |
подробнее о нем см. § 3.1). Ли |
||||||
нейная функция регрессии имеет следующий вид: |
|
|
||||||
|
~у (х) = М (Y) + |
г {у/х) їх — М (X)], |
|
(2.24) |
||||
где |
|
|
|
|
|
|
|
|
|
г |
Ш |
= р ^ |
^ |
Ш |
. ; |
(2.25) |
|
|
|
|
Ох |
|
°х |
|
и Y; |
|
°".v. |
— средние квадратические |
отклонения величин |
X |
|||||
соу(ХУ) = М{\Х — M(X))IY |
— М (К)]} — к о в а р и а ц и я,, |
раз |
||||||
мерная |
характеристика |
(или момент) |
связи случайных |
величин. |
||||
Ковариация тем больше, |
чем сильнее |
связь между X и Y. |
Если |
|||||
cov(XF) = 0, то случайные величины X и У некоррелированы.
Для независимых случайных величин всегда cov = 0, но обратное |
||||
не верно. |
|
|
|
|
К о э ф ф и ц и е н т , |
к о р р е л я ц и и |
(нормированная ко |
||
вариация) |
р = С 0 У |
^ ^ |
— безразмерная |
характеристика связи |
случайных |
величин |
X |
и У (—1 <; р ^ 1). Чем теснее и ближе к |
|
линейной |
эта связь, тем | р | ближе к единице. Для случайных ли |
||||||||
нейно связанных величин X и |
У = |
аХ + Ь, так как D(Y) = |
|||||||
= a-D{X), |
а СОУ(ХУ) = |
|
aD(X), |
|
|
|
|
|
|
|
" — |
aD |
(X) |
- |
а |
• = |
і . |
. . |
|
|
І a \ax a* |
|
(sign a) • 1, |
||||||
|
|
|
\a\ |
|
|
|
|||
т. e. p = 1 со знаком при a.
К о э ф ф и ц и е н т |
р е г р е с с и и |
У по |
X |
г(у/х) — |
это |
тангенс угла наклона |
к оси х прямой регрессии |
у(х) |
(2.24), т. |
е., |
|
грубо говоря, средней прямой, проходящей вдоль облака опытных точек* на плоскости ху через его центр тяжести. Для линейно-
зависимых случайных величин X и У = аХ + |
b г (у/х) — а. |
||||||
Зная r(ylx) |
и М(Х), М (У) — центр тяжести |
облака — можно |
|||||
прогнозировать |
(в среднем) величину У для заданного |
значения х. |
|||||
Зная р(х, у), |
можно |
судить, |
насколько связь между |
X |
и У близка |
||
к линейной |
(|p| - » - |
1) или |
же насколько эта связь |
вообще суще |
|||
ственна (р - > 0). По ковариации СОУ(ХУ) Ф 0 трудно судить о зна чимости связи, поскольку, изменив размерность, ковариацию мож но сделать сколь угодно большой.
* И м е ю т с я в в и д у т о ч к и , п о л у ч е н н ы е п у т е м н е п о с р е д с т в е н н о г о и з м е р в ' н и я н а о п ы т е в е л и ч и н ы Y п р и р а з л и ч н ы х з н а ч е н и я х X.
§ 2.4. Определение числовых характеристик и законов распределения функций случайных величин
Числовые характеристики функций случайных аргументов. Ес ли случайная величина X имеет дифференциальный закон рас пределения f(x), то для любой функции ф(Х) математическое ожи дание и дисперсия соответственно равны:
М [ф (X)] = Мщ = |
] Ф (*) / (*) dx; |
|
|
— о о |
(2.26) |
||
со |
оо |
||
|
|||
D [ Ф (X)] = J [<р (А-) — М ф 1 2 / (х) dx = |
J V (х) f (х) dx — M<р- |
|
|
Приближенное вычисление числовых характеристик методом линеаризации. Метод сводится к тому, что функцию Ф(Х) разла гают в ряд Тейлора в окрестности точки х° = М(Х) и сохраняют только члены первого порядка
Ф (X) ^ ср Ш (X)] + Ф' Ш (X)] -[X — М (X)].
К полученной линейной функции применяют известные правила (2.11), (2.12), (2.21) и получают:
М [ф (X)] =ё ф [М (X)]; |
|
|
| |
£ [ ф ( Х ) ] ^ { Ф ' Ш (X)]}2 D (X); а[ср (х)]^\ц>' Ш ( Х ) ] | а ( Х ) . | |
(2.27) |
||
Для функции многих случайных аргументов |
cp(Xlt |
Х 2 , |
Х„) |
имеем: |
|
|
|
М [ф ( Х ъ Х 2 , . . , Х„)] ~ ф [М (X,), М(Х2),М |
(Xn )]; |
j |
|
™ « 2 № ; ) D ™ + |
|
| |
( 2 , 8 , |
+Ї І І + 1 ( £ ) ( £ Н < * ' * А
где все производные йф°/дХй вычисляются в точке f*° - М ( Х ^ *° = М (Х2 ), .... = М (*„)},
на что указывает индекс 0; cov(Xf t X,) = pw 0(Xf e )a(X,) [см формулу (2.25)]. Если случайные величины Хх , Хп некоррелированы, т. е. р,а = 0, то дисперсия
d № ( х „ х „ x j ] = <TJ а 2 ( U ) ' C та =
(2.29)
k= I
15
Иногда в литературе эту формулу называют з а к о н о м |
н а к о п |
л е н и я п о г р е ш н о с т и . |
|
Закон распределения монотонной функции одного случайного |
|
аргумента. Пусть X — непрерывная случайная величина, |
распре |
деленная в интервале [а, Ь] по закону f(x). В случае монотонно воз
растающей или убывающей на [а, |
Ь] функции |
Y = ф(Х) диффе |
ренциальный закон для Y имеет вид |
|
|
dtp-1 |
(у) |
(2.30) |
|
I dyidx |
|
dy |
I |
|
где x — ф 1(y) — функция, обратная у = ф(х).
Закон распределения функции двух случайных аргументов.
Если функция Z = Ф(Х, Y), где ( X , Y) — система двух случайных величин, распределенная по закону f(x, у), то интегральный закон, согласно формуле (2.1),
G (г) = Р {Z < z) = Р {Ф (XY) < z) = Р {(XY) Є S).
Здесь область S представляет собой ту часть плоскости х, у, на которой Ф(ХУ) •< г. Отсюда
G (z) = J J / (х, г/) dwfy и g (г) = G' (г). |
(2.31) |
Величина г в выражение для G(z) входит неявно, через пределы ин
тегрирования. |
|
|
|
|
|
|
Композиция |
законов. |
Закон |
распределения суммы |
независи |
||
мых |
случайных |
величин |
называется к о м п о з и ц и е й |
з а к о |
||
н о в |
р а с п р е д е л е н и я |
слагаемых. Известно, что если |
X |
|||
и У независимы, то fix, |
у) — Д(л:) Д(г/) [1]. Подставляя |
это fix, |
у) |
|||
в формулу (2.31), получаем выражение для композиции двух зако нов, или плотность распределения вероятностей случайной вели
чины Z = |
X + Y: |
|
|
|
|
|
|
8 (г) = |
[ |
Д (*) Д (z — х) dx, |
|
|
|
|
—СО |
(2.32) |
|
|
|
|
|
||
|
|
|
|
оо |
|
|
или |
g{z) |
= |
l |
Д (z — у) Д (у) dy. |
Интегралы |
такого |
типа |
называются с в е р т к о й двух функций |
||
и часто обозначаются символом Д * Д.
Г л а в а 3.
О С Н О В Н Ы Е З А К О Н Ы Р А С П Р Е Д Е Л Е Н И Я С Л У Ч А Й Н Ы Х В Е Л И Ч И Н
§ 3.1. Законы распределения непрерывных случайных величин
Равновероятный закон. По этому закону распределяются по грешности измерений, связанные с округлением до ближайших целых значений (при считывании со шкал, при взвешивании в ус ловиях ограниченного количества мелких разновесок и т. п.). Он справедлив для отклонений параметров от номинальных значений в пределах полей допусков, связанных с погрешностями изготов ления, когда упомянутые отклонения вызываются следующими ус ловиями производства: а) наличием фактора, равномерно изменя ющегося во времени и доминирующего над совокупностью всех остальных (равномерная разлаДка оборудования, инструмента, равномерное ухудшение точности контрольных приборов и т. д.); б) использованием при изготовлении изделий большого числа раз личающихся по точности инструментов и оборудования или в слу чае смещения изделий многих партий, когда средние значения рас сматриваемого параметра для каждой партии равномерно распре деляются по заданному полю допуска [6].
Равновероятные дифференциальный и интегральный законы рас пределения имеют следующий вид:
О |
при х < а, |
|
/ (JC) = J—ї— |
при а ^ х ^ б , |
(3.1) |
Ь—а |
|
|
О |
при х > Ь; |
|
О |
при х ^ а, |
|
F(x) = \b-i |
при a^Zx^.b, |
(3.2) |
1 |
при х ^> Ъ. |
|
Математическое ожидание (2.10), дисперсия (2.16) и среднеквадра тичное отклонение (2.17) соответственно принимают вид
М (X) = (a+b)/2- D (X) = (6 - а)2 /12; о = (Ь — а)12 V3.
Экспоненциальный закон. По экспоненциальному закону рас пределяется неотрицательная случайная величина — время работы изделия до отказа при условии, что отказы не связаны со старе нием, износом (а вызываются, допустим, случайными внешними факторами, такими, как внезапные чрезмерные увеличения меха нической, электрической, тепловой и |Щгочих ^нагрузок). Сущест-
венным здесь является то, что вероятность возникновения отказов, за период ДЛдолжна определяться только длиной этого промежут ка и не зависеть от предшествующего периода времени (от преды стории). Это часто выполняется для сложных технических систем, содержащих большое количество элементов [7, 8]. Вообще экспо ненциальное распределение реализуется во всех случаях, когда рассматриваемая случайная величина допускает следующую ма тематическую интерпретацию. На положительную ось (или пло скость, объем) случайным образом бросаются точки так, что сред
Р и с . 2. Э к с п о н е н ц и а л ь н ы й з а к о н р а с п р е д е л е н и я , XY — Зк2.
нее число точек, попадающих на единицу длины, постоянно и рав но Я. Тогда случайная величина отрезка, на который не попало ни одной точки, будет распределена по экспоненциальному закону. Примером такой случайной величины может быть время между испусканием двух следующих друг за другом ос-частиц при радио активном распаде радия.
Дифференциальный закон для экспоненциально распределенной случайной величины X имеет вид (рис. 2)
(3.4).
(3.5)
Экспоненциальный закон обладает интересным свойством: для него математическое ожидание и среднее квадратическое откло нение совпадают:
М (X) = 0 = ПК; D (X) = 1/Я2. |
(З.бЬ |
Это обстоятельство часто помогает установить справедливость данного закона в конкретных практических случаях.
Нормальный закон (Гаусса). Большая распространенность нор мального закона на практике имеет под собой глубокую теорети ческую основу, которая выражена в центральной предельной тео реме А. М. Ляпунова. Сущность теоремы заключается в следующем. Если случайная величина Z может быть представлена в виде суммы достаточно большого числа п независимых (или слабо зависимых) случайных величин X,
Z= |
«=i |
(3.7) |
|
|
каждая из которых сравнительно мало влияет на сумму (т. е. сре ди слагаемых нет доминирующих), то для большого класса случай
ных величин Xt (практически для всех встречающихся |
в инженер |
|
ной практике) величина Z будет распределена по нормальному За |
||
гс |
п |
|
кону с параметрами M(Z) = 2 М(Х; ); |
D(Z) = ^D(Xi). |
Важным |
i= і |
£= і |
|
здесь является то, что результат оказывается не зависящим от за конов распределения величин Xt. Если же известно, что Xt рас пределены по нормальным законам с М(Хг) = УИг и D(Xt) — oi, то даже в случае зависимости Xt между собой закон распределения величины Z будет нормальным при любом п и любых Mi, ait при чем
|
M(Z)=2iMji;= і |
|
(3.8) |
|
|
|
|
|
|
D(Z)= |
І а ? + 2 " І ' |
І |
P t i a t a J t |
|
|
і = і |
/=і |
І=І+І |
|
где p;-j — коэффициент корреляции |
случайных величин X, и Xj |
|||
'[см. формулу (2.25)]. |
Из выражений |
(3.8) |
следует, что при компо |
|
зиции нормальных законов (т. е. при сложении независимых слу чайных величин, распределенных нормально) получится снова нор мальный закон. Законы, обладающие этим свойством; называются у с т о й ч и в ы м и .
Z, |
По нормальному |
закону, наряду с рассмотренными величинами |
в частности, распределяются погрешности экспериментов, вре |
||
мя |
восстановления |
некоторых ремонтируемых изделий, суммар |
ная наработка ряда восстанавливаемых изделий до капитального ремонта, погрешности изготовления (отклонения. параметров от
номинальных |
значений в пределах полей допусков [6, 9—14]). |
|||
Плотность |
вероятности, нормального распределения имеетвид |
|||
.(рис. 3): |
|
|
|
|
t, \ |
1 |
Г 1 fx— ЛП21 |
— о о < л : < о о , |
(3.9) |
|
|
|
||
