Файл: Игнатов, В. А. Статистическая оптимизация качества функционирования электронных систем.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— х2(Р2+ Р 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 г м состоя