Файл: Игнатов, В. А. Статистическая оптимизация качества функционирования электронных систем.pdf

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

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

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

Добавлен: 19.10.2024

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

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

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

При малых объемах выборок реальное распределение Кг имеет положительную асимметрию, а при увеличении числа экспериментов асимптотически стремится к нор­ мальному с параметрами (2-84). Поэтому приближенно целесообразно оценивать только нижнюю доверительную границу для Кг, так как ясно, что верхней доверительной границей служит 1.

Пример 2-8. Пусть т і= 1 0 _3 ч-1, аі=0,2-10”3 ч-1, г = 7, пі2 = 0,5 ч“1, а2=0,1 ч-1, у = 0,9, тогда ти=0,998003, о/г =1,78- ІО-3, К г м іш = т .к —l,65oft —0,995063.

Использование второго члена ряда Тейлора для

2

9

mh =

^ т у ' Дает незначительную поправку 3,97 X

Х Ю - 20/о.

Итак, в этом параграфе рассмотрены простейшие рандомизированные модели изменения надежности и по­ казано, что основной трудностью их использования в точ­ ных методах анализа является отыскание законов рас­ пределения нелинейных многомерных функций. Сущест­ венно упростить анализ и расчетные соотношения позволяет приближенный метод, использующий линеари­ зацию и асимптотическую нормальность функций от ве­ личин, имеющих конечные математические ожидания. Так как погрешности приближенного метода асимптоти­ чески убывают с ростом объема выборок и при числе экспериментов r> 10 не превышают 10—15%, этот метод можно успешно применять для исследования качества изделий с помощью рандомизированных моделей.

2-6. Инженерный метод анализа качества систем

Кроме моделей, рассмотренных в примерах предыду­ щих параграфов, в инженерном анализе качества функ­ ционирования обслуживаемых и необслуживаемых устройств может найти широкое применение еще одна простая рандомизированная модель. Эта модель имеет два работоспособных состояния устройства ( п = 2, каче­ ство функционирования оценивается по уровням прогно­ за и отказа), учитывает взаимосвязанное появление вне­ запных, постепенных и перемежающихся отказов и по­ строена на гипотезе, что АР и ПО состоят из двух (т = ='fe= 2) обобщенных операций с приведенными значе­ ниями интенсивностей ja и £. С помощью этой модели

62


получим аналитические соотношения для инженерной оценки характеристик качества устройств.

Дифференциальные уравнения, описывающие динами­ ку изменения качества, имеют вид:

(і)=

- Оі. +

А0)

(0 +

вд (t) + т (0 + *Рь (<);

 

Р'і (t) =

loPo (0 -

(А, + 1], +

Ѳ+

V) Д (0 + Р-Рз (f);

 

P\{t) =

X0p0t + (А, +

7і,) р, (/) -

(f);

(2-85)

 

 

 

 

 

 

 

P'» (0 =

(0 -

2р.рз (0 ; Д 4(0 =

ѵд (0 6р4(0 ;

 

Для определения характеристик отдельных процессов (ухудшения качества, АР, ПО) необходимо рассматри­ вать переходные режимы. Начнем с определения харак­ теристик ухудшения качества. В этом случае в системе (2-85) необходимо положить ѵ, |, р, Ѳ= 0, тогда в ней останутся первые три уравнения, решение которых при самых общих начальных условиях позволяет найти ха­ рактеристики качества необслуживаемого устройства.

Вероятность безотказной работы

 

 

 

,

(2-86)

где <Хі=Аг+

 

 

 

 

плотность этой вероятности

 

 

 

/ (0 — *«Pe

е

+ а. (Pi + ЪРо

е ~ ^ \

(2-87)

интенсивность

отказов

 

 

 

 

Л (t) = f(t) [p{t)]~l;

 

(2- 88)

среднее время безотказной

работы

 

 

 

 

 

 

(2-89)

дисперсия времени

 

 

 

2

 

 

 

 

t

 

 

 

 

вариация времени

 

 

(2-90)

 

 

( 2- 91)

 

Vt = otT0- 1.

 

63


Чтобы найти закон распределения и числовые харак­ теристики определяющего параметра, методом наимень­ ших квадратов вычислим оптимальные квантовые значе­ ния х*і из условия

д

 

— л:в

2+\m x(t =

T0) - x l]^ =

дх*і

 

где

 

= 0,

t = ÖT2,

 

 

(2-92)

 

х*0ра(t) +

х*,р, (t) +

х*грг(f);

 

тх {t) =

 

Po = P o ^ ; Рг (О = ^

ао

+

 

 

 

 

 

а1

 

 

+ ( Р, +

 

е"“*'.

Р-2 ( 0 = 1 -

Po (0 -

Л (О-

Математическое ожидание и дисперсия выходного

параметра

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

тх (0 = Ц

0П1 р,-(0,

з* (0 =

 

 

 

і—0

 

 

 

 

 

'

= s

(-^*1 out)2 Рі (0

- л £ ( 0 -

(2-93)

 

i=0

 

 

 

 

 

Вычислим коэффициенты асимметрии и эксцесса вспо­

могательного

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

 

 

 

 

 

л = =>п,(0 ; ң __^

з

 

(2-94)

 

 

"г (О

 

°г4 (0

 

 

 

где

 

 

 

 

 

 

 

л ,( о = 5 ] [ —

-тх

 

/МО;

 

 

 

 

 

 

 

і=0

 

 

 

 

 

,W = ] і]

 

 

 

 

(2-95)

 

 

1=0

 

 

 

 

 

 

 

2

тх (0

 

 

 

 

 

 

 

 

 

 

 

(=0

]*/М 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

(2'96)

 

S [

 

 

 

 

і=а

 

 

 

'

 

64


и найдем плотность распределения выходного параметра

ю W1 =

[ ф' (z) “ Ж фі & + #

Ф5 (Z) - ..." •

 

 

(2-97)

Таким образом, по известному диапазону допустимо- )

го изменения

определяющего параметра

и статистиче-1

ским путем найденным интенсивностям пересечения слу- | чайной функцией выбранных уровней квантования с по- 1 мощью соотношений (2-86) — (2-97) можно определить » все основные характеристики качества функционирова­ ния элементов и необслуживаемых устройств.

Перейдем к определению характеристик АР и ПО. Положим в системе (2-85) ло, ль ta, Хь ѵ, I, 0 равными нулю, тогда, учитывая, что вероятность восстановления работоспособности устройства

V(t) = l —p2(t)—p3(t),

 

из оставшихся уравнений получим:

 

V (^) =

1 — (1 -j-

(2-98)

плотность вероятности

восстановления

 

 

V (t) =

(2-99)

интенсивность проведения АР

 

jx±(f )

= 1 X ^ ( 1 — p i ) - 1.

( 2 - 1 0 0 )

Аналогично, полагая в системе (2-85) т]о, Ль ^о» ^ь Ць

V, Ѳ равными нулю,

для вероятности завершения ПО по­

лучим:

 

 

 

Я(*) = І - - ( 1 + 8 ) < Г 6*,

(2-101)

плотность этой вероятности

 

интенсивность ПО

v{t) = \ 4 e ^ ,

(2-102)

 

 

 

 

 

+

(2-103)

Как следовало ожидать, время АР и время ПО име­

ют гамма-распределения с параметрами

 

Tpf= 2 p ._1,

ар — 2р-~2,

=

=

 

Ѵр = V. = 0,5 У ~ 2 .

(2-104)

5—385

65


Чтобы получить характеристики качества функциони­ рования обслуживаемого устройства в установившемся режиме эксплуатации и в дежурном режиме, рассмотрим состояние статистического равновесия. Полагая в систе­ ме (2-85) p 'i(t)— 0 и учитывая условие нормировки, най­ дем:

р

___ ____________ __________ і ц .

( 8 - f

v +

0 , 5 c t i ) ________________________

 

 

I [I* ( Ѳ + ѵ

+ а 0 +

0 . 5 * !

-

0 , 5 Х 0)

+

1 . 5 А 0 +

ѵ ) + ~ " *

 

 

 

+jl .5*0*,] +

2v(fc ( а о- 0

, 5 А о) ;

 

(2-105)

n

 

6 ІР +

 

( а о

 

0 . 5 Х 0)

 

 

 

 

 

1 “

V + а0 +

0,5а1 -

 

0 ,5 \0)]+ 1 ,5 [Х0 (8+ ѵ )+ ^

 

 

 

+ “ оа і] + 2ѵ (х ( а 0 — 0 , 5 Л 0)

!

 

( 2 - 1 0 6 )

р

I

 

Х0| К “і +

Хо (ѵ +

8)1

 

 

 

2 ~

[ * о « 1 + Х . ( Ѵ + Ѳ ) ] ( 1 , 5 \ 0 + | І . ) + | * (Ѵ)0+ 0 , 5 Х 0)

[ 6

( Х 0-

* , ) + 2 ѵ Х 0] '

Коэффициент готовности устройства

 

 

(2-107)

 

 

 

 

__ I* { h o +

0.5Хо) (Х„ — *,) +

К « 1

+

Х„ (у +

9)]} Р

Л г~

X, [*,*, + X, (ѵ +

0)]

 

 

(2-108)

коэффициент простоя на АР

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/Ср=1,5 Р 2,

 

 

 

 

(2-109)

коэффициент простоя на ПО

 

 

 

 

 

 

 

 

 

 

 

2ѵ[Д (-По+ 0 . 5 Х о)

 

 

 

(2- 110)

 

 

 

і;[«о“і+ х0'.(ѵ +

8)1

 

 

 

 

 

 

 

 

Характеристики тх, ах, (о(х) определяют по форму­ лам (2-93) — (2-97). Особенностью является то, что они не зависят от времени. Коэффициент надежности обслу­ живаемого устройства, используемого в дежурном режи­ ме, рассчитывается, как обычно:

А(0 = КХ7 7 ‘ JP W dx,

(2-111)

t

 

где характеристики Р(т), Т0 и Кг определяют соответст­ венно по формулам (2-86), (2-89), (2-108).

Полученные в этом параграфе соотношения позволя­ ют по относительно простой исходной статистике оценить наиболее часто применяемые в инженерном анализе ха-

66