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

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

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

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

Добавлен: 19.10.2024

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

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

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

Ряд Грама—Шарлье для плотности вероятностей

Z(t) имеет вид [Л. 34]:

.о, [г(*)]=Ф<»>(г) - 4 ф(4) (г) +4ф<6)и-•••>

где А — коэффициент асимметрии; Е — коэффициент экс­ цесса распределения тоц функция Лапласа

2

2 du, ф т { г ) = * Ш .

О

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

Ü (О] =

Ф (z) - 4 Ф(8) (г) 4 " |-

Ф<4) (2) - -

Учитывая линейность преобразования (2-14), полу­

чим:

 

 

 

 

шИ*)1

1

 

Ф '( г ) - 4 Ф (4,(2) +

4 Ф (,' Р ) - . . . ] .

«1 (О

 

 

(2-17)

 

 

 

 

Из-за погрешностей квантования процесса Х(і) мате­ матическое ожидание процесса z(t) отличается от нуля, а его дисперсия не является единичной, поэтому с уче­ том квантования

А т3(0 . Е

m*(t) - з ,

(2-18)

°2 ( 0

4 ( о

 

где

(2-19)

(2- 20)

3 2 ( 0 =

Так как значения Ф ^г) табулированы [Л. 2І, 35], тб определение закона '<й [х(/)] сводится к определению ко­ эффициентов ряда (2-17) по формулам (2-18) — (2-21), известными величинами в которых являются оптималь­ ные квантованные значения х*гопт и вероятности Pi(t).

Рассмотрим марковскую модель перемежающихся отказов устройств. В этом случае процесс Х(1) является стационарным [Л. 22, 26] и нас могут интересовать сле­ дующие характеристики: вероятность нормального функ­ ционирования до первого отказа, плотность этой вероят­ ности, закон распределения определяющего параметра, коэффициент готовности устройства, закон распределе­ ния времени пребывания в неработоспособном состоя­ нии (длительность сбоя), математическое ожидание и корреляционная функция определяющего параметра.

Предположим, что интервал [х0, Хі] изменений опре­ деляющего параметра X известен. Выберем промежуточ­

ный уровень квантования х2< х і < х 0.

Обозначим интен­

сивности пересечения процессом Х(і)

уровней квантова­

ния Хо, Хі, х2 с положительными и отрицательными про­ изводными соответственно через Ѳо и г)о, Ѳі и т)і, Ѳг и r\z-

Вероятности событий x0<X(t), Xi<X(t)<Xo, x2< X (t)< < x i, X(>t)<x2 обозначим черезPo(t), Pi(;t), Pz(t) и Pz\t)-

Следовательно, вероятности Pi(t) и P2(t) характеризу­ ют вероятность пребывания устройства в момент време­

ни t в работоспособном состоянии, а вероятности

Po{t)

и Рз(і) — в неработоспособном.

уравнений, описываю­

Система дифференциальных

щих динамику перемежающихся

отказов, имеет вид:

 

P'At) = --n*P*{t) + K P M

 

р '1(0 =

Чо^о (0 — (Чі +

Ѳ0) Р1(0 +

Öl-Ра (t)’,

(2-22)

- Р \ (t) =

ъ р г (t) - (Ъ +

Ѳ.) Р 2 (t) +

02Р 3 (t);

 

p \{t) = %PAt) — ^ P 3{t)-

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

Используя изображения Р і ( s ) как производящие функции, можно простым способом найти среднюю дли­ тельность нормального функционирования

Г ъ ъ р г

Ъ (*І2 + Ql) — MoPj — ^lo hl + S0)P2 /р QO'

1 W A + Чо (Чі + Ѳо)(Ч2 + Ѳі) — '']оМД2+Ѳ1) - тЧоЧА

28


to среднюю продолжительность сбоя

. 6г Г7!g ("'ll —Ар) Ч~ воДг^Т — 9оѲіР2 — ѲоРі (Дг + 9і)1 ~Ь

9 г [ДоДг^о + До (Ді + Ѳо)(Дг + Ѳі) —

 

+ Дг [До9 0 Р 2 — Д о Д іЛ — 7]„ (уц + 6 0) Р г 1

(2-24)

— ДрѲр (Да + Ѳі) — Д о Д А І

 

где Рі и Р2— вероятности пребывания устройства в-со­ стояниях Si и S2 при t — 0. Такие характеристики особен­ но важны, например, для самолетных вычислительных устройств, используемых в навигационных системах.

Найдем характеристики качества устройства в со­ стоянии статистического равновесия: коэффициенты го­ товности Кѵ и простоя /Сп, закон распределения опреде­ ляющего параметра и корреляционную функцию процес­ са.

Решая систему (2-22) при условии Р'і (і)=0 и ис­ пользуя условие нормировки для предельных вероятно­ стей, получим:

Р .

__ Др^1^2

.ДоДЛ

. Д0Д1Д2

1

г

 

 

 

 

 

где

z = ѲоѲіѲ2+ тіоѲіѲг-ртіогрѲг+ Г|оГ]іТ]2.

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

Кг = Р. + Р,

Д о М г +

ДоД А

Kn= l —KT.

Ѳо0іѲ2 + 7]0ѲА

+

д„д А + Д„ДіД2

 

 

 

 

(2-25)

Для определения закона

распределения X(t) необхо­

димо найти оптимальные квантованные значения х*,-, і —

= 0,3. Выберем x*0=Xö+yÂXi, х*і = Хі-|-уАхі> х*2= х 2+

+уАх2, х*з— х2— уАх2, где Ахі— х0—хі, Ах2= хх—х2, и най­ дем оптимальное значение у из условия совпадения ма­ тематического ожидания тх и выборочного среднего т*х процесса (методом моментов), тогда

 

То

т*х — х0Я0хіР 1— х22+ Р 3)

(2-26)

 

&х> (Р0- Р , ) + Ах2

2— Р 3)

 

 

Из

выражения

(2-17), по значениям х*ІОлт и вероят­

ностям

Рі, і = 0,3

определяется закон распределения

м(х). Отличие заключается в том,

что его моменты не

29



зависят от времени, например

 

з

з

 

 

м-х=

3^=^] хіотРі т х .

(2-27)

 

1=0

і= о

 

Корреляционную функцию Кх (т) процесса X(t)

приб­

лиженно

найдем

следующим путем. Аппроксимируем

процесс

X(t) случайной последовательностью

прямо­

угольных импульсов с детерминированной амплитудой А = хМакс—хтш. Длительности импульсов выберем, рав­ ными длительностям положительных выбросов X(>t) за уровень тх, а интервалы между импульсами выберем равными длительностям отрицательных выбросов за уро­ вень тх. Функцию Kx(t) найдем с помощью обратного преобразования Фурье от спектральной плотности А(ш) случайной последовательности прямоугольных импульсов с постоянной амплитудой (Л. 35]

^

/уЛ__

2 (хмэкс Хмнн)2^

n c [ i - S i ( » ) ] [ i - . 6 i H ]

 

 

У ,

W2 (л + ^)

1 — И h И

где

(to)

и І2(со)— характеристические функции

дли­

тельностей положительного и отрицательного выбросов X(t) за уровень тх. Например,

СО

 

 

 

 

 

^ И = J eimt

ve-^dt =

 

 

 

 

о

 

 

 

 

 

где р — интенсивность

пересечения

процессом

уровня

тх с положительной производной.

 

 

 

 

После необходимых преобразований

 

 

2 (Хмакс ' ‘ Хмин)а Д*

 

eia>t d(o=

 

+

р)[со2+ (7 +

(J.)2]

 

 

(Хмаке

Хмиң)2 рХ

 

 

 

О + р )2

 

 

 

 

Уточним Хмакс и xMHH с помощью

числовых характе­

ристик (2-27).

 

 

 

 

 

Обозначим Хмако==tnx-\-kax, хмш= тхkax,

восполь­

зуемся условием Д х ( 0 ) = <г2х, тогда

 

 

 

(2-28)

k=m?=l’

 

 

 

 

Следовательно, для приближенной оценки качества при перемежающихся отказах устройств достаточно по

30


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

Итак, в данном параграфе рассмотрены некоторые основные типы марковских моделей появления отказов элементов и необслуживаемых устройств, построенные с помощью квантования по уровню выходных случайных функций. Достоинство этих моделей в том, что они по­ зволяют связать изменение качества устройств с измене­ нием их параметров, учесть взаимосвязанное появление внезапных и постепенных отказов, оценить различные характеристики качества устройств при перемежающихся отказах. При построении моделей на вид случайных функций, характеризующих изменение параметров эле­ ментов и необслуживаемых устройств, не накладывается никаких ограничений. Они могут быть линейными, вы­ пуклыми вверх или вниз, иметь точки перегиба и т. п. С учетом точности исходных данных сложность структу­ ры моделей можно оптимизировать по требуемой точно­ сти расчетов характеристик качества.

Отличительная особенность рассмотренной марков­ ской аппроксимации заключается в том, что вероятност­ ные характеристики качества элементов и необслуживае­ мых устройств оценивают по минимуму статистики о времени пересечения случайными функциями фиксиро­ ванных уровней квантования. Вопросы получения и об­

работки статистических данных будут рассмотрены в

гл.

7.

2-3.

Модели технического обслуживания систем

Техническое обслуживание обычно изучают метода­ ми теории массового обслуживания и теории восстанов­ ления [Л. 3, 6]. Эти методы позволяют при некоторых упрощающих анализ предположениях определять веро­ ятностные характеристики длительности обслуживания: вероятность завершения технического обслуживания за фиксированное время, плотность этой вероятности, ма­ тематическое ожидание, дисперсию времени обслужива­ ния и др.

Наряду с этими характеристиками значительный ин­ терес представляют и те, которые получают с учетом

31

влияния обслуживания на характеристики качества об­ служиваемого устройства. Очевидно, что в таких иссле­ дованиях необходимо применять вероятностные модели, в которых фиксируют значения определяющих парамет­ ров устройств до и после обслуживания.

Кроме того, большинство реальных видов техничес­ кого обслуживания, как правило, включает не детерми­ нированное, а случайное число операций, которое зави­ сит от технического состояния устройства, поступающе­ го на обслуживание, вида обслуживания, характера отказов, наличия запасных частей и необходимого инст­ румента, квалификации обслуживающего персонала и т. п. Поэтому необходимо привлекать и математические модели со случайным числом операций. Рассматривая это число как параметр рандомизации [Л. 27, 28], в дальнейшем такие модели будем называть рандомизи­ рованными моделями технического обслуживания, а са­ ми процессы — рандомизированными.

В этом параграфе мы рассмотрим простую вероятно­ стную модель технического обслуживания, которая, учи­ тывает изменение определяющего параметра устройства в процессе обслуживания [Л. 15], и рандомизированную модель [Л. 27, 28]. Предполагая длительности отдельных операций экспоненциально распределенными величина­ ми, с помощью решения дифференциальных уравнений А. Н. Колмогорова, описывающих динамику обслужива­ ния, найдем аналитические соотношения для вероятно­ стных характеристик процессов обслуживания.

Рассмотрим восстановление работоспособности уст­ ройства после отказа, включающее т последовательных операций: отыскание неисправного узла или элемента, приготовление инструмента и необходимых деталей для ремонта, разборка неисправного узла, ремонт, сборка отремонтированного узла, настройка и т. п. Предполагая интенсивности проведения отдельных операций извест­ ными, найдем вероятность К(т) восстановления устрой­ ства за фиксированное время т, плотность н(т) этой ве­

роятности, математическое ожидание Тѵ и дисперсию времени восстановления.

Обозначим через S* такое состояние отображающей марковской системы, когда выполнено і операций вос­

становления і = 0,т. Дифференциальные уравнения для вероятностей Р%{х) пребывания системы в 5 г м состоя­