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

книги из ГПНТБ / Волков Е.Б. Основы теории надежности ракетных двигателей

.pdf
Скачиваний:
55
Добавлен:
25.10.2023
Размер:
12.65 Mб
Скачать

х\ 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 ). В двумерном

случае (N1)

можно показать, что

справедливы соотношения [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 =

= 1tN. Тогда для случайного вектора

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

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