Файл: Волков Е.Б. Основы теории надежности ракетных двигателей.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 11.04.2024

Просмотров: 206

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

между любыми ее компонентами

равен одному и тому же

числу Q, то

 

 

0ЛГ|1,2,...,ЛГ-1 — Q

N — 1

(1.94)

+ ( N — 2 )б

1

 

Существенным обстоятельством здесь является то, что в случае линейной регрессии (1.85) соотношения (1.90) — (1.92) строятся

с помощью р, и Hcffjll так, что знания вида закона распределения

не требуется.

 

 

 

 

важ­

Наряду с множественным коэффициентом корреляции

ное значение

имеет

частный

коэффициент корреляции

Qiv, jv—111,2...... n- 2

между

t.у и tN- 1

при фиксированных

t\ =

= .v'i;...; tN _ 2 — xN_2. Выражение для него имеет вид [1]

 

QjV,iV—1|1,2...,jV—2=

буу.лг—1|2,3,...,Л '-2 ~~ 6лГ|1,2.....IV—2 9/У -1|1,2,...,/У —2

(1 .9 5 )

1

(1

ejv|l,2....Дг—2) (1 — бдг—111,2,...,/V—2)

 

 

Соотношение (1.95)

является рекуррентным и позволяет, вычис­

лив частный коэффициент корреляции

 

 

Qn .n - 111'

Q/VJV-l

 

(1.96)

 

 

1 С1— 6Jv,l) 0 — 6w-l,l)

найти OjV, JV—I11,2, 6iV, JV—111, 2, з и т. д.

Интересно отметить, что если при рассмотрении двух случай­

ных величии t{ и tj

коэффициент корреляции qh^ [ —1,1], то в со­

вокупности tN= (t 1 ,

t2, ..., Cv) это

свойство может

не

выпол­

няться. Пусть, например, при N>2

коэффициенты

корреляции

между любыми двумя компонентами равны одному

и тому же

значению q, а матрица ||а,ц|| положительно определена.

Тогда

согласно работе [1]

р е [—1| (N—1),

1].

 

 

Множественный коэффициент корреляции может быть вычис­ лен с помощью частных коэффициентов корреляции по формуле

6l|2,3,...,/V =

* (1

6i,2)0 6?i3|2)-"(l

6i,/v|2... л'- i) -

(^-97)

Рассмотрим

теперь

некоторые типовые

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

слу­

чайных векторов.

1. 1.10. Непрерывные распределения

А. Многомерное нормальное распределение

Пусть непрерывный случайный вектор tN имеет вектор

сред­

них jlljv= (pi,..., щу) и ковариационную

матрицу ||a,-j||.

Тогда

распределение tN называется //-мерным

н о р м а л ь и ы м,

если

его функция распределения имеет вид

 

 

28


F { x n )= P (—оо< / ‘; < л',-, v i = \ , N ) =

---- СО

---- СО

П аУн

(1.98)

/-1

 

где

|crij’|

— определитель матрицы ||cr,’-i|l, обратной к ||а,ц||;

 

 

 

 

_

 

N

IV

 

 

 

 

 

 

 

 

Q{yiv) =

^ 2 я11(У1—ъ ) { У 1 —Ну);

 

 

 

 

 

 

 

 

i=l /=1

 

 

 

 

 

e'i — элементы

матрицы

||crij’||,

определяемые из

соотноше­

ния

(1.88).

 

 

 

 

 

 

 

 

 

 

 

Запишем выражение (1.98) в виде

 

 

 

 

 

F ( x JV)=

 

 

 

 

 

Xi- - *1.,

y i = T7n

) =

 

 

 

\

 

 

 

 

 

°i

 

'

 

 

 

= ^ { / ' n )= P ( — с о < 2 ; < Л ; , V i — l , N ) ,

 

 

где

Gi

полагается

большим нуля,

a

z*= (U\ц)!ай h{=

= (Xi—\ii)loi;

/ziY=

(Л ь Л2>. . . ,

Л лг).

 

 

 

 

 

Случайные

величины

 

имеют средние значения,

равные

нулю, и дисперсии, равные единице.

 

 

 

 

 

Найдем ковариацию Oij между г,- и 2 ,-. По определению

 

 

 

»/,= f

j

yiyjf{yi,yj)dyidyjt

 

 

где /(у,-,

ijj) — совместная плотность

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

tji

и у,.

 

Из соотношения (1.98)

следует, что

 

 

 

 

/ ( УпУ})-■

 

 

 

 

■ехр

 

й + у) ^аУлУ]

 

(1.99)

 

V

1

вЬ

 

2(1-0?,)

 

 

 

 

 

 

 

 

 

где Qij — коэффициент корреляции

между

t\ и tj,

и,

следова­

тельно,

 

 

 

 

 

 

 

 

 

 

 

 

<j

 

----о-

\

\

У1У]ЫР

у] + у) ~ ^QiiViUJ dy-Myj=

 

V 1-

Q?,

 

J

 

 

 

 

20 - е ? /)

 

 

 

 

 

J

— во — со

 

 

 

 

 

 

 

СО

 

 

 

2

 

 

 

 

 

 

 

 

= "^Г f

 

 

 

 

 

l-ef^+ Q /yiT /Jexp^—

 

 

 

111

\ У?ехР ( ----

e x p (-

dMyi+ 1 i-eJyX

 

— CO

29



х\ 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