
книги из ГПНТБ / Волков Е.Б. Основы теории надежности ракетных двигателей
.pdfх\ j Ui ехр (—т~)Аехр (~т) dy‘dX•
—со — со
Здесь использованы замена Х= (г/,- — г/,-р,у)/ф 1— qj. и соотношеппе
У1 ~ 2QijUiyj + У)= ( 1 - 6;;) У? + (г/у - а,7г/,-)3.
Последний из двух двойных интегралов, входящих в выражение для CT,-j, равен нулю, так как
^[^,1 = 0 = - ^ ^ £/,ехр(—
— 00
Вследствие равенства
/ и И 1 = о 5 = - Д . | у ; е х р ( - ^ ) « < !, , = 1
получим Oij = Qij.
Таким образом, вместо выражения (1.98) может быть исполь зовано соотношение
|
- |
I |
ГТ77 |
hl |
!‘n |
|
|
N |
|
|
|
exp |
|
|
|
|
|||||
F u |
= F(hN)= - |
le j |
\ |
\ |
|
n |
^ - . d - 100) |
|||
|
|
(2л)Л/3 •> |
J |
|
|
|
|
|
||
где |
| q« | — определитель |
матрицы |
||д^||, |
обратной |
к |
корреля |
||||
ционной матрице ||д,7||; |
|
|
|
|
|
|
|
|||
|
,у |
N |
|
|
|
|
|
|
|
|
|
&(Улг) = V |
V егуУ/У/, |
//;= |
-Л;/~- — при i= T jV . |
||||||
|
jmJ |
JmU |
|
|
|
°i |
|
|
|
|
|
i=\ |
)~\ |
|
|
|
|
|
|
—► |
|
Вычисление вероятности попадания случайного |
|
|||||||||
вектора tN |
||||||||||
|
|
|
|
|
Л'г , |
xi\, |
т. е. определение величины |
|||
в выпуклый многоугольникП|х,-, |
||||||||||
|
|
|
|
|
/=1 |
|
|
|
|
|
|
|
|
Р {xi < t[ < х], v i = \ , N ) , |
|
|
|||||
осуществляется |
с |
помощью соотношений |
(1.98) |
или |
(1,100) |
|||||
и (1.57). |
интерес |
некоторые частные случаи |
функции |
|||||||
Представляют |
распределения F(hN). Пусть все значения /ii=(*i—pti)/сГг в вы ражении (1.100) равны между собой (hi = h), коэффициенты корреляции также одинаковы (Qa= Q)■ Тогда согласно работе [63] выражение (1.100) существенно упрощается:
30
F (л'лО — F ( h N ) = |
j" © (г /) F " |
У / |
q + |
h |
d у |
(I. ЮП |
||
|
|
|
|
l |
I — 0 |
|
|
|
или F (xN)= F (fiN)= |
1---- ^ |
[ F |
у / 2 ( 1 — Q)— h |
iJidy. |
(1.102) |
|||
|
|
|
|
|||||
|
1 я |
.) |
|
/ q |
|
|
|
|
Здесь |
|
—c© |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ll |
|
|
y ^ F 1* -1(]/ 2y) e - "a; |
a(z) = |
- ^ |
e--''72; |
F( / i ) = |
\ |
v\y)dz. (i. 103) |
||
|
|
у 2л |
|
|
|
J |
|
|
Функция плотности вероятности и функция одномерного нор мального распределения табулированы в работе [63]. Соотноше
ния (1.101) |
и (1.102), как и (1.100), определяют функцию рас |
||||
пределения |
наибольшей |
из |
компонентов вектора |
zx ={z\, |
|
z2, ..., zx ). В двумерном |
случае (N—1) |
можно показать, что |
|||
справедливы соотношения [46, 63]: |
|
|
|||
|
/Г (//1,Л2) = / Г (Л1)/7 (Л2)Ч -f 'b(h1,fi2,y)dy; |
(1.104) |
|||
|
|
|
о |
|
|
F (А,, /,,) = f <fl|> + f (|‘а) - |
T [7г., я».) - |
T(lh , а,„) - Ь; |
(1. 105) |
|
|
|
V , |
|
V! |
С;„ |
(1. 106) |
|
|
|
лАшА |
|
|
|
|
Здесь |
|
|
v= l |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
d F ( h l t ho) |
__ |
d^F ( h l t ho) |
|
|
У(/ >1, КУ) = |
|
д х {д х 2 |
|
||||
|
|
|
dQi о |
|
|
||
|
1 |
|
h \ + Ло — 2Л1/г2у |
(1. 107) |
|||
|
|
|
ехр |
2(1 -if-) |
|||
2л т^С— у2 |
|
||||||
Т |
1 |
Г |
е х р [— Л2 (1 + |
х 2)/2] dx — функция Оуэна |
|||
|
|
|
1 + а-2 |
|
|
|
|
(табулирована в работе [63]); |
|
|
|
|
|||
|
Но— |
/г iQ! о |
|
h\ — /ioQ12 |
|
||
|
|
|
ЛЛа |
h-2Т 1 |
QJ 2 |
|
|
|
hi \ |
|
1 — 0J2 |
|
|||
FW(li) — производная |
v-ro порядка |
от |
функции |
Лапласа |
|||
(1.103). В выражении (1.105) принимается |
Ь = 0, если |
h\ho>§, |
|||||
или если hlh2 = 0, |
но |
h\ + h2^ 0 . В |
противном случае |
прини |
|||
мается 6 = 0,5. |
|
|
|
|
|
|
|
31
Рассмотрим случай, |
когда каждая из |
величин U, |
входящих |
в выражение (1.100), |
представляет собой |
разность |
п, = Пг—Н* |
двух случайных величин П* и Нгсо средними П?-, Н,-, дисперсиями
ал |
>ан |
i |
11 коэффициентом корреляции Гл.п .. |
/ |
|
1 J |
Очевидно, что если Пг и Н, распределены нормально, то соот ношение (1. 100) для функции распределения случайного вектора
0v=(^i, U,---, |
tN) |
не |
изменится |
при |
подстановке |
значений |
|||||
|
1ц—--- |
X j- jU i- E j) |
|
= |
Х( _ „. |
j_[08) |
|||||
|
|
*1/ |
2 |
, |
2 |
' о |
|
|
|
|
> |
|
|
у |
сп. + |
|
°н.— ■‘гп.н/°п/ан/ |
|
|
|
|
||
где |
[i,- — ГГ,-— Н,-; |
ai ='j/rЗп. + |
Зн . ~ н ^ л ^ н . , |
|
|
||||||
и величин |
|
|
|
|
|
|
|
|
|
|
|
|
д |
= М [(и, - |
|
н-,){иj ~ ty)] = — |
гп и, + |
|
|||||
|
|
|
|
|
|
|
|
<w |
|
1 1 |
|
|
J Н ; |
Н ; |
Гн.Н- |
|
aiaj ' гП,Н; |
|
|
•гн.Пу) |
[. 109) |
||
|
|
|
|
ЧЧ |
|||||||
|
|
|
|
|
|||||||
где |
з .= |
л. з2н^ _ 2гпjHfiijHр |
|
|
|
|
|
||||
гл.,Пу, r H|.,Hyi гп.,н_,- |
и rH.uj — коэффициенты |
корреляции случай |
|||||||||
ных величин П; и П,-, Hj |
|
и Н,-, |
П,- |
и Hj, |
Ыг- |
и П3при |
/=1, N, |
||||
/ = TJV. |
|
|
|
|
|
|
|
|
|
|
|
Из соотношений (1.26) и (1.28) |
следует, |
что при всех q,; = 0, |
|||||||||
и всех оfj=1 соответственно |
|
|
|
|
|
|
|||||
|
|
|
* |
|
|
|
|
|
|
|
(1-110) |
|
F {Xn )\qu ~Q -= П F (III) И F (xN)\Qi r l= F (hm\ |
||||||||||
|
|
|
/-1 |
|
|
|
|
|
|
hm = m\Y[hi минимальное из Jti.
!</<//
Б. Многомерное экспоненциальное распределение
|
Пусть случайный вектор |
имеет компоненты tu каж |
|
дая из которых распределена экспоненциально: |
|
||
|
F (Xi) — |
= 1 —ехр(—%iXi), |
|
причем |
|
|
|
|
F(Xi) = Р (U>Xi) =ехр(—KiX{), |
(1.111) |
|
где |
— параметр экспоненциального распределения, |
равный |
|
— , а р, = М[/,]. |
|
|
|
У1 |
|
|
|
32
Тогда функция |
Fr (xN) определяется как [100] |
|
|||
F ( x n ) |
~ |
v <= l,;V)-~exp[ —2^max(x;, s,.)] , |
|||
|
|
|
|
ie i S| |
|
где {S} — множество векторов s= (sb s2, ..., s,\), |
а каждое s; = 0 |
||||
или 1, |
но (si, s2, ■■■, sN) Ф (0, 0 ,..., |
0). Для любого s g {5} вели |
|||
чина |
max (xit |
S;) |
— максимальное |
значение xit |
для которого |
Si=l-
Наиболее изучено к настоящему времени двумерное экспо ненциальное распределение. В работе [101] показано, что его функция плотности вероятности имеет вид
f { x ltx2)= |
Х|Хг |
ехр |
ГхjX[ |
Хц}\2 |
2 1-гб1а |
'‘•р'2 V^XV |
|
|
1— |
||||
(' — 012) |
1 — 01 2 |
|
||||
где qi г — коэффициент корреляции t\ и t2\ |
|
|||||
/„(■)— модифицированная |
функция |
Бесселя (табулирована |
||||
в работе [93]). |
|
|
|
|
||
В. |
Распределение порядковых статистик |
|
||||
Пусть из совокупности |
с функцией |
распределения F (х) = |
||||
= Р (О ^ Д ^х ) |
и |
функцией |
плотности |
вероятности f(x)=F'(x) |
извлечена выборка п наблюдений t\, t2, ..., tn случайной вели чины t. Расположим значения Г,- в порядке возрастания (в вариа ционный ряд):
|
*0) ■<tpi) |
/До- |
(1.112) |
Величина 7,-, где /=1, п, называется г-й п о р я д к о в о й |
с т а т и |
||
стикой. |
Величины Г, рассматриваются как взаимонезависимые |
||
случайные величины. Однако в ряду |
(1.112) вследствие |
упоря |
|
дочения |
величины ti оказываются |
зависимыми. Коэффициент |
|
корреляции Г(7„) и 1(h) при больших п равен [74] |
|
/,k 1 — m/n m 1 — k/n
Функция распределения in-го члена ряда согласно работе [74] определяется как
/=' W = P ( 0 < ^ ) < ^ ) = ( j J)ffl [ г/т_1(1 —y)ddy —
'А
= / f (x)(m,rf+l). |
(1.113) |
312 |
33 |
|
|
|
|
|
jY “ 1(1 — y)d dy |
|
|
|
|
||||||
Здесь |
J 0irn/Z-f 1)= |
—^--------------------- |
|
. |
(1.114) |
||||||||||
|
|
|
|
|
|
)' |
|
\ — y)i! <iy |
|
|
|
|
|||
|
|
|
|
|
|
b |
|
|
|
|
|
|
|
|
|
— неполная бета-функция; |
|
|
|
|
|
|
|
|
|
||||||
|
|
d = |
n — tn\ |
|
|
|
n\ |
|
|
|
|
||||
|
|
\ m |
J |
ini ( n — m)\ |
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|||||
Используем известные тождества: |
|
|
|
|
|
||||||||||
|
|
|
d |
|
|
|
|
|
|
а |
m — 1 -I- v |
|
|
||
J о (m, d + 1).- V |
l |
" ] fl«-»( l _ |
fl)v = |
fl« 2 |
Н-е)*. |
||||||||||
|
|
||||||||||||||
|
|
|
v=0 ' |
' |
|
|
|
|
v=0 |
|
|
(1.115) |
|||
первое из которых записывается также в виде |
|
|
|||||||||||||
|
|
|
|
||||||||||||
|
|
|
У о[т, |
d -}- 1 ) = |
B i (//, 0, d), |
|
|
|
|
||||||
где Bi (», 0, с/)= |
^ |
b (//, 0, v) |
и Ь\п, 0, v) — ^ п |
|
1 — 0)” — обо- |
||||||||||
значения. |
|
v= 0 |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Тогда |
|
|
/ r (A'm)='Bi[«, |
F(x), |
п — in] |
|
|
;i. ив) |
|||||||
и |
F ( x j = |
F [лT |
2 |
(m ~ l Г '’j 11 - F l*Г |
(1■1'17) |
||||||||||
|
|
|
|
|
|
v= 0 |
|
|
|
|
|
|
|
||
Соотношение (1.116) может быть использовано во многих |
|||||||||||||||
случаях. |
Пусть, |
например, |
случайный |
вектор |
tN<=R(x") |
имеет |
|||||||||
N компонентов, каждый из которых распределен равномерно на |
|||||||||||||||
отрезке (У, f], т. е. имеет функцию распределения |
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
х — Т |
|
|
|
(1.118) |
||
|
|
|
|
|
|
|
|
|
т— т |
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Без ограничения |
общности |
можно |
рассматривать |
отрезок |
|||||||||||
[ Л ? Ы 0,1]. |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Тогда F(x)=x. |
Пусть теперь t\, t2, . . . t N — независимые |
слу |
чайные величины, каждая из которых равномерно распределена на отрезке [0, 1]. Пусть далее выражение (1.112) представляет упорядоченные значения случайных величин U. Эти порядковые статистики разбивают отрезок [0, 1] на N+1 непересекающихся
интервалов, |
имеющих длины Ui = t{i)\ |
U2 = t(2)—i(ij;...; Un+i = |
= 1—tN. Тогда для случайного вектора |
Un+i на основе выраже |
|
ния (1. 116) |
можно получить [82] соотношение |
34
F u K „ ) = P ( U , > x r, V / = 1, W + 1)= ( 1 - - v ' * . Y \
/~1
В случае, когда рассматривается отрезок [О, Г],
|
Л Ч 1 |
, /V |
|
s |
* |
^(-vjVbl) = \ 1 |
/= 1 |
|
|
|
Основным в теории порядковых статистик является Д-мерное распределение Дирихле, подробно исследованное в работе [99].
1. 1. И. Дискретные распределения
А. Многомерное биномиальное распределение
Пусть рассматривается случайный |
вектор ТУеЛ?(ЛГ), каждая из |
||
компонент tj которого представляет |
возможное число отказов |
||
в «,• испытаниях |
(0^^-^/г,-), проводимых по схеме |
Бернулли |
|
[23]. Случайные величины t( независимы, причем |
|
||
Р {ti = |
v' )= ( * ) P"'_V/ (1- :р/ ^ ^ b I«/• Р/. v,\ |
(1.П9) |
а функция распределения /,• имеет вид
д';
*('?/, Р/, v,)=Bl(//,-, Р„л:,), (1. 120)
v=0
где Р,- —вероятность успешного исхода (события А{) .в одном (любом) из п, испытании, одинаковая в каждом испытании с номером 1, 2 , /г. При этом исходы испытаний считаются независимыми.
Функция распределения Д- согласно выражению (1.43) опре деляется как
Д(ГгЛ/)= Р(0 |
v i = l , W ) = V |
V ... V П A(«,.,P,.,v,.). |
|
|
.лай |
md |
=*0 /=>X |
|
=0 |
|
В некоторых случаях л-,- по условиям задачи оказываются связан ными между собой. Так, например, может потребоваться, чтобы
некоторая функция ф(л-,х) =0. Тогда
/V
(Т. 121)
2* |
35 |
где v;cfi; х'. — значение Х; удовлетворяющее условию cp(x,v) = = 0, a Q — множество значении \у, определяемое условием
ф(хл-) =0.
Наиболее естественно распределение (1.121) возникает при исследовании последовательных систем, элементы которых испы тываются отдельно от системы.
Б. Многомерное отрицательное биномиальное распределение
Пусть в условиях предыдущего раздела рассматривается слу чайная величина Б — число испытании до наступления d,- отка зов, имеющая отрицательное биномиальное распределение:
|
|
" 'о - р л |
|
d i + v . . . |
|
|
Тогда случайный вектор tN имеет функцию |
jV-мерного |
отрица |
||||
тельного биномиального распределения вида |
|
|
|
|||
|
.V |
N |
|
|
|
|
^ЧА'лг)— 2 |
V |
V/ — 1 |
PV |
-rf<( 1— P,-)rf', |
( . |
) |
|
п rf, — 1 |
1 122 |
||||
v, = rf, |
|
1=1 |
|
|
|
|
если /; независимы.
В. Многомерное гипергеометрическое распределение
d| |
Из совокупности N неразличимых предметов, в числе которых |
|||||
относится к первой категории, |
di — ко |
второй и т. |
д., dM |
|||
к М-й категории, а N— (d\ + d2 + ... +dN) — к |
нулевой |
катего |
||||
рии, извлекается выборка объема |
п. |
Вероятность P(/; = v;, i— |
||||
= |
1, М) того, что в выборке окажется |
г, предметов из |
первой |
|||
категории, мг — из второй и т. д., |
из M-i’i выражается как |
|||||
|
м |
|
( п - У а Л |
|
||
|
П |
(^1) 1 |
Л1 |
|
||
|
/=1 |
4 |
1=1 |
|
||
|
|
|
|
|||
|
|
|
N |
|
|
|
|
Функция Af-мерного гипергеометрического |
распределения |
||||
дается выражением |
|
|
|
|
|
\лг
36
Здесь if — случайная величина, равная числу предметов /-н кате гории в выборке п.
На основе метода производящих функций можно показать,
что справедливо тождество |
|
||
V .. .. |
м |
|
|
V м |
|
||
М |
|
jLd |
|
|
\)! i=1 |
j |
|
у |
T J |
|
|
|
|
|
Используя его, находим
где D = ' £ d i. i- 1
Г. Многомерное отрицательное гипергеометрическое распределение
Пусть теперь ti — число предметов, извлекаемых по одному из совокупности N до тех пор, пока не будет извлечено dj пред
метов /-й категории (/=1, М) |
(число предметов нулевой катего- |
м |
|
рии m0 — N — 2 ^ dj). Тогда |
распределение tM дается выраже- |
/=i |
|
нием [38] |
|
P(^/ = v/, |
v i = l , M ) = |
м
где V,- > d;■ k = V (vt. — di). i=i
Выражение дает возможность вычислить вероятность того, что для извлечения d{ предметов первой категории, d2— второй и т. д. потребуется лщ гг, • •. извлечений. Это эквивалентно собы
тию, состоящему в извлечении k предметов |
нулевой категории. |
Поэтому P(t=k) =P(^=V t,Y i'=l, М), где |
t — случайная вели |
чина, представляющая собой число извлечений до появления
k предметов нулевой категории. Функция распределения tM имеет вид
37
A-,., v / = l , T I ) = V P ( / = ^ ) -
A’ =0
•r/Vf
= 2 ••
i L=dL |
, .U“ d .U |
Л1 |
|
Xi-Cti). |
|
1.1.12. Полиномиальное распределение
Пусть проводятся п испытаний по схеме Бернулли, т. е. в каж дом из п независимых испытаний возможен один из двух исхо
дов: успех (событие Ао) |
пли отказ |
(событие |
Л0); вероятность |
||
Р = Р(Л0) = 1—q, где ^= Р(Ло) |
в каждом испытании одинаково. |
||||
|
|
_ |
ft |
т. е. отказ подраз- |
|
Но дополнительно к этому событие Л0=11Лг, |
|||||
деляется на k видов отказа (событие |
j - i |
___ |
причем |
||
Ai при |
г = 1, /г), |
||||
Л,- — несовместны. |
|
|
|
|
|
к |
|
|
|
|
|
Тогда Р -4-g= P-j- V |
q; = 1. |
Рассмотрим |
случайный |
вектор |
|
/~i |
|
|
|
|
|
th=(t), tk), где ti — возможное число отказов г-го вида
в п испытаниях. Компоненты вектора th линейно зависимы, так как t\-r ... +th + to = n, где to — возможное число успехов в п испытаниях. Тогда согласно работе [23]
____ |
к |
f ft |
' |
Р<Уо=''о'=Р(Л=Т;, ¥ / = 1 , / Н = « ! |
П <?/Pv |
П v/! |
|
|
/~1 |
0 |
, |
ft
где v0 -
/=1
Функция распределения tk имеет вид
F (**) = jU |
••• |
2 |
Я! П '7;,'PV° |
/ ( п V/!l > |
v=.\r0 |
VA=0 |
v^= 0 |
/ = 1 |
41=0 |
где x.'о =
i=1
Отсюда видно, что справедливо следующее интересное тож дество:
38
2 |
... 2 |
п 1?р-/ пV, |
|
|
VL к |
V. |
/ = 1 |
/ 1=0 |
|
2 |
v/ < D = « _ x 0 |
|
|
|
1= |
1 |
|
|
|
= |
2 f ” j Pv(l — P)"~v=Bi (/?., Р, /}). |
(1.123) |
||
|
v-Л’о |
|
|
|
При вычислении функции полиномиального распределения часто оказывается удобным известное соотношение
k |
k |
|
F (■**) < П р |
< х () = 11 |
Bi (л, Р„ jef). |
/=1 |
/=1 |
|
1. 1.13. Отрицательное полиномиальное распределение
Пусть теперь tQ— случайная величина, равная числу испыта ний до получения т0 успехов, a ti — число появляющихся при этом отказов — событий А {. Тогда отрицательным полиномиаль ным распределением называется выражение [38]
P(*/ = v„ v i= T 7 k )= P (t0 = m0)= Г(-л^ рт |
Y l |
я \ |
|||
|
|
^ ( т 0)П v;! |
/=1 |
|
|
ft |
|
ft |
/= 1 |
|
|
|
|
|
|
||
где п# = т0+ У ] V,-; |
Р = 1 — V q.\ Г(/г.^.)— гамма-функция. Это |
||||
i=i |
|
i=i |
того, |
что |
успешное |
выражение позволяет найти вероятность |
|||||
|
|
|
|
k |
\ |
испытание будет иметь место в /?г0-й раз при |
[ т0-j- ^ |
•ом |
|||
испытании. |
|
|
|
i=i |
|
|
|
|
|
|
|
Соответствующая функция распределения имеет вид |
|
||||
F (-*ft)=Р {t, < х„ |
v i = 1, k ) = Р (/„< х 0)= |
|
|||
-»Ч |
ЛА |
____ |
|
|
=2 p ^/ = v/. v i = l , k ) ,
»1”0
к
где х 0— т0-|- ^ х г •
1.1.14. Некоторые сведения из теории случайных функций
Пусть (Я, В, Р) — вероятностное пространство. Случайной функцией td ) называется конечная вещественная функция, кото-
39