1. Сравнительный анализ вероятностей

катастрофы летательного аппарата.

1.1Постановка задачи задания № 1

Летательный аппарат (ЛА) состоит из

-m двигателей с вероятностей отказа P₁, P₂,…P_m ;

-n дублирующих систем энергосбережения с вероятностей отказа

P_1Э , P_2Э ,... P_nЭ ;

-N c вероятностей отказа Р_с каждая.

Катастрофа наступает, если выходит из строя любая (r+1) и более двигателей, либо если все системы энергоснабжения, либо если хотя бы одна из N вспомогательных подсистем.

В случаи отказа любого r из m двигателей катастрофа наступает с вероятностью Р_D .

Определить вероятность катастрофы ЛА и сравнить ее с вероятностью катастрофы ЛА без дублирующих систем(один двигатель с вероятностью катастрофы P₁ , одна система энергосбережения с вероятностей отказа P_1Э и N вспомогательных подсистем с вероятностей отказа Р_с каждая ), предполагая, что все упомянутые выше системы и подсистемы ЛА функционируют независимо друг от друга.

В обоих случаях сравнить вероятности катастроф, связанных с отказом

- двигателей;

- систем энергосбережения;

- вспомогательных подсистем.

Дано

m = 5; Р₁ =6∙10^-4, Р₂ =5∙10^-4, Р₃=7∙10^-4, Р₄=2∙10 ^-4, Р₅=4∙10 ^-4

r=4 Р_D=0.1;

n=4 Р_1Э=3∙10^-4, Р_2Э=4∙10^-4, Р_3Э=10 ^-4, Р_4Э=6∙10 ^-4;

N=3∙10³ P_c=6∙10^-9.

Решение.

Математическая часть

Введем обозначение событий:

- D₁, D₂, D₃, D₄ - отказ 1-го, 2-го, 3-го и 4-го двигателей соответственно;

- В₁, В₂, В₃, - отказ 1-й, 2-й, и 3-й системы энергоснабжения соответственно;

• С_i - отказ i-ой вспомогательной подсистемы, i = 1,2,...,N;

• Е_к - катастрофа;

-E_kd,E_кэ, E_кc - катастрофы, связанные с отказом двигателей, систем энергоснабжения и вспомогательных подсистем соответственно.

А) Рассмотрим случай ЛА с дублирующими системами:

В этом случае:

Е_K=Е_KD+Е_KЭ+E_КС . (1.1)

Перейдем к противоположным событиям, будем иметь:

= (1.2)

Из равенства (1.2) в силу соотношения двойственности получим:

Е_K= ∙ ∙ ^(1.3)

Тогда вероятность катастрофы будет определяться по формуле:

P(E_K)=1- P( )=1-P( ∙ ∙ ) (1.4)

Из равенства (1.4) в силу независимости событий Е_KD , Е_KЭ, E_КС получим:

P(E_K)=1- P ∙ P( )∙ P(E_KC)=1-(1-P(E_KD))∙(1-P(E_KЭ))∙P(E_KC)). (1.5)

Рассмотрим структуру событий Е_KD, Е_KЭ, E_КС и найдем их вероятности, то есть вероятности катастроф, связанных с отказом

• двигателей Е_КD

• систем энергоснабжения Е_KЭ

• вспомогательных подсистем Е_KC

1) Рассмотрим структуру событий Е_KD и найдем P(E_KD)= P_KD

Так как событие Е_KD - это событие, состоящее в том, что катастрофа произошла из-за отказа двигателей, а по условию задачи катастрофа, связанная с отказом двигателей наступает, если выходят из строя любых (r+1) и более двигателей из m двигателей, а в случае отказа любого г из m двигателей катастрофа наступает с вероятностью Р_D . Значит:

Е_KD= Е_KDr+ Е_{KD≥ (r+1)}, где

Так как в нашем случае число двигателей m = 5, r = 4; то r + 1 = 4 + 1 = 5.

Значит:

Е_KD= Е_KD4+ Е_KD≥5 где:

Е_КD4 - событие, состоящее в том, что катастрофа произошла из-за отказа любого r =4 из m=5 двигателей;

Е_KD>5 - событие, состоящее в том, что катастрофа произошла из-за выходы из строя любых (r + 1) = 5 и более двигателей, а в нашем Е_KD>5= Е_KD5 - это событие, состоящее в том, что катастрофа произошла из-за отказа пяти двигателей. Из этого следует, что:

Е_KD≥5 = Е_KD5 = D₁∙ D₂∙ D₃∙ D₄∙ D₅ (1.6)

В свою очередь катастрофа, связанная с отказом ровно r = 4 двигателей (при работающих остальных), не обязательно влечет за собой катастрофу (ас вероятностью P_D ), значит

E_KD4=E_K∙ E_D4

Тогда:

E_KD= E_KD4+ E_KD≥5= E_K∙ E_D4+ E_KD≥5

Так как события E_KD4, и E_KD≥5 несовместны, то

P(E_KD)=P(E_KD4+ E_KD≥5)= P(E_KD4)+ P(E_KD≥5)=P(E_K∙ E_D4)+P(E_KD≥5)

а для нашего случая и учитывая (1.6), получим:

P(E_KD)=P(E_KD4+ E_KD≥5)= P(E_KD4)+ P(E_KD≥5)=P(E_K∙ E_D4)+P(E_KD≥5)= P(E_K∙ E_D4)+P(E_KD5) = P(E_K∙ E_D4)+ P(D₁∙ D₂∙ D₃∙ D₄∙ D₅)

С другой стороны, катастрофа, связанная с отказом ровно r=4 двигателей при работающих остальных из пяти имеющихся у ЛА по условию задачи, есть следующее событие:

E_D4 = D₁∙ D₂ ∙ D₃∙D₄ + D₁∙ D₂ ∙ D₃∙ ∙ D₅ + D₁∙ D₂ ∙ ∙D₄∙D₅ +

+ D₁∙ ∙ D₃∙D₄∙D₅ + ∙ D₂ ∙ D₃∙D₄ ∙ D₅ (1.8)

то есть не работают 5-й, 4-й, 3-й, 2-й, 1-й двигатели из пяти, имеющихся у ЛА.

Замечание.

Тот факт, что события E_KD4 и E_KD≥5 несовместны, можно доказать следующим образом:

E_KD4∙ E_KD≥5 =< согласно (1.7) >= E_K∙ E_D4∙ E_KD≥5=< согласно (1.6) >= E_K∙ E_D4∙ Е_KD5 = =< согласно (1.6) и (1.8) = E_K(D₁∙ D₂ ∙ D₃∙D₄ + D₁∙ D₂ ∙ D₃∙ ∙ D₅ + D₁∙ D₂ ∙ ∙D₄∙D₅ + D₁∙ ∙ D₃∙D₄∙D₅ + ∙ D₂ ∙ D₃∙D₄ ∙ D₅) ∙ D₁∙ D₂∙ D₃∙ D₄∙ D₅ = E_K ((D₁∙ D₂∙ D₃∙ ∙D₅ D₁∙ D₂∙ D₃∙ D₄∙D₅)+( D₁∙ D₂∙ ∙ D₄∙D₅ ∙ D₁∙ D₂∙ D₃∙ D₄∙D₅)+( D₁∙ ∙ D₃∙ D₄∙D₅ ∙ D₁∙ D₂∙ D₃∙ D₄∙D₅)+( ∙ D₂∙ D₃∙ D₄∙D₅ ∙ D₁∙ D₂∙ D₃∙ D₄∙D₅)+( D₁∙ D₂∙ D₃∙ D₄∙ ∙D₁∙ D₂∙ D₃∙ D₄∙ )=

= E_K((D₁∙ D₁)∙( D₂ ∙D₂)∙( D₃∙ D₃)∙( D₄∙ D₄) ∙( D₅ ∙ ) + (D₁∙ D₁)∙( D₂∙ D₂)∙( D₃∙ D₃)∙( D₄∙ )∙( D₅ ∙ D₅)+(D₁∙ D₁)∙( D₂∙ D₂)∙( D₃∙ )∙( D₄∙ D₄) ∙( D₅ ∙ D₅) +(D₁∙ D₁)∙( D₂∙ )∙( D₃∙ D₃)∙( D₄∙ D₄) ∙( D₅ ∙ D₅)+( D₁∙ )∙( D₂∙ D₂)∙( D₃∙ D₃)∙( D₄∙ D₄) ∙( D₅ ∙ D₅)

Используя тот факт, что A∙A = A и A∙ =Ø, получим

E_KD4∙ E_KD≥5 =E_K((D₁ ∙D₂ ∙ D₃ ∙D₄∙ Ø) + (D₁ ∙D₂ ∙D₃∙ Ø∙ D₅ )+ (D₁ ∙D₂ ∙ Ø ∙ D₄∙D₅) + (D₁ ∙ Ø ∙ D₃ ∙D₄∙D₅) + (Ø ∙D₂ ∙D₃ ∙D₄∙ D₅)) = Ø

А как известно, что, если произведение двух событий равно невозможному событию (пустому множеству), то такие события являются несовместными.

По определению условной вероятности имеем:

P(E_KD)=P(E_K / E_D4)∙P(E_D4)+P( )

а в силу независимости событий D_i , i= , далее имеем:

P(E_K / E_D4) ∙ P(E_D4)+ P( )

Используя (1.7) и несовместимость его (E_D4) слагаемых

P(E_K / E_D4)∙(P(D₁∙ D₂ ∙ D₃∙D₄ ) + P(D₁∙ D₂ ∙ D₃∙ ∙ D₅ ) + P(D₁∙ D₂ ∙ ∙D₄∙D₅)+ P(D₁∙ ∙ D₃∙D₄∙D₅) + P( ∙ D₂ ∙ D₃∙D₄ ∙ D₅))+ )

В силу всех независимых событий D_i , i= и потому, что

P( )=1-P(D_i), получим далее:

P(E_K / E_D4)∙ [(P(D₁)∙P(D₂)∙P(D₃)∙(P(D₄) ∙(1-P(D₅ ))+ (P(D₁)∙P(D₂)∙P(D₃)∙(1-P(D₄)) ∙P(D₅ )+P(D₁)∙P(D₂)∙(1-P(D₃) )∙P(D₄) ∙P(D₅) +P(D₁)∙(1-P(D₂) )∙P(D₃)∙P(D₄) ∙P(D₅) +(1-P(D₁)∙ P(D₂)∙ P(D₃)∙P(D₄) ∙P(D₅)]+ )

Так как P(D_i)=P_i , i= и P(E_K / E_D4)=P_D , имеем

P(E_KD)=P_D∙[ P₁∙ P₂∙ P₃ ∙P₄∙(1-P₅)+P₁∙ P₂∙ P₃ ∙(1- P₄)∙P₅ + P₁∙ P₂∙(1- P₃)∙ P₄∙P₅ + P₁∙(1- P₂)∙ P₃∙ P₄∙P₅ +(1- P₁)∙ P₂∙ P₃∙ P₄∙P₅]+ P₁∙ P₂∙ P₃ ∙ P₄ ∙P₅=P_D∙[ P₁∙ P₂∙ P₃∙ P₄+ P₁∙ P₂∙ P₃∙P₅+ P₁∙ P₂∙ P₄∙P₅+ P₁∙ P₃∙ P₄∙P₅+ P₂∙ P₃∙ P₄∙P₅]∙(1-5P_D)∙ P₁∙ P₂∙ P₃ ∙ P₄∙P₅≡P_KD;

Если выполняется условие P«P_D для всех i= (1.9)

и учитывая, то что значение вероятности случайного события есть величина, меньшая единицы, то

P₁∙ P₂∙ P₃ ∙ P₄∙ P₅→0

А значит тоже

(1-5P_D)∙ P₁∙ P₂∙ P₃ ∙ P₄∙ P₅→0

И тогда имеем

P(E_KD)≡P_KD≈P_D∙( P₁∙ P₂∙ P₃∙ P₄+ P₁∙ P₂∙ P₃∙P₅+ P₁∙ P₂∙ P₄∙P₅+ P₁∙ P₃∙ P₄∙P₅+ P₂∙ P₃∙ P₄∙P₅) (1.10)

Подставив значения, данные из условия задания, получим

P(E_KD)≡P_KD≈P_D∙( P₁∙ P₂∙ P₃∙ P₄+ P₁∙ P₂∙ P₃∙P₅+ P₁∙ P₂∙ P₄∙P₅+ P₁∙ P₃∙ P₄∙P₅+ P₂∙ P₃∙ P₄∙P₅)=

=0.1∙(6∙10^-4∙5∙10^-4∙7∙10^-4∙2∙10^-4+6∙10^-4∙5∙10^-4∙7∙10^-4∙4∙10^-4+6∙10^-4∙5∙10^-4∙2∙10^-4∙4∙10^-4+6∙10^-4∙7∙10^-4∙2∙10^-4∙4∙10^-4+5∙10^-4∙7∙10^-4∙2∙10^-4∙4∙10^-4)=

=0.1∙10^-16∙(420+840+240+336+280)=21.16∙10^-16 (1.10)

2) Рассмотрим структуру событий Е_кэ и найдем P(E_КЭ)=P_КЭ

E_КЭ≡ B₁∙ B₂∙ B₃∙ B₄ - катастрофа, связанная с отказом всех трех систем энергоснабжения (п= 4 по условию задачи).

В силу независимости всех событий B_i , i= имеем

P(E_КЭ) ≡P(B₁∙B₂∙B₃∙ B₄)=P(B₁) ∙P(B₂) ∙P(B₃) ∙P(B₄) =P_1э∙P_2э∙P_3э∙P_4э (1.12)

Подставив значения, данные из условия задания, получим

P(E_КЭ)≡P(B₁∙B₂∙B₃∙ B₄)=P(B₁) ∙P(B₂) ∙P(B₃) ∙P(B₄)=P_1э∙P_2э∙P_3э∙P_4э =

=3∙10^-4∙4∙10^-4∙10^-4∙6∙10^-4=120∙10^-16 (1.13)

3) Рассмотрим структуру события е_кс и найдем P(е_кс) = P_кс.

Событие Е_кс наступает, если отказывает хотя бы одна из вспомогательной подсистемы, значит

е_кс≡C₁+C₂+…+C_N=

В силу закона двойственности

е_кс≡ = ∙ ∙…∙ =

в силу независимости событий , i= получим

P ( ) ≡P( =P( ) ∙ P( )∙…∙ P( )= = 1-P(C_i))

Так как P(C_i)=P_{c ,} i= получим

P ( )= = 1-P_с)=(1-P_c)^N

тогда

P(е_кс)=(1- P ( )=1-(1-P_c)^N≡P_KC

Если выполняется NP_C<<1=>

P ( )=( 1-P_c)^N=1-NP_C+ P_C²-…(-1)^N P_c^N ≈ 1-NP_C (1.14)

Подставив значения, данные из условия задания, получим

P(е_кс) 1-1+NP_C=NP_C=3∙10³∙6∙10^-9=18∙10^-6 (1.15)