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

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

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

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

Добавлен: 19.10.2024

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

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

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

16*

235

Рис. 7-1. Номограмма для определения доверительных границ

 

где i'ik — время ідо внезапного от­

 

каза

k-ro устройства из і-го состоя­

 

ния;

гі — число устройств, внезап­

 

но отказавших из і-го состояния,

 

 

г- 1

 

 

тх=

 

 

і—о

 

На рис. 7-2 показаны продол­

 

жительности ti и t'i, которые фик­

 

сируются для точечной и интерваль­

 

ной

оценки интенсивностей ухуд­

 

шения параметров и интенсивно­

что Г і + г ' і = Ышг п і . Кроме

стей внезапных отказов. Очевидно,

того, если

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

работоспособности устройств

после отказа, то в формулах (7-6) и

(7-7) t'i + г'і =ІѴИ.

 

 

Дисперсии и вариации оценок максимального правдоподобия для рассматриваемых интенсивностей по-прежнему находят по формулам

(7-2), (7-3) с учетом замены параметра г

на п или г ' і , интерваль­

ные оценки производят по формулам (7-4)

и (7-5) с учетом коэффи­

циентов

 

тгя - т ; и и я - щ

’ напр « еР'

22Г4

■ ’Пг

(7-8)

(Л^и ~772i) ?ih

(Nx —іПі) S tih

k-\

ft=i

Пример 7-1. Предположим, что на надежность испытывается пар­ тия приемно-усилительных ламп из Na= 20 шт. Выходным парамет­ ром служит крутизна X характеристики. Лампы работоспособны, если Х > х кѵ=0,5хт (хт — максимальное среднее значение крутизны характеристик для всей партии в начальный момент времени). До­ пустим, что учитывались две градации качества функционирования ламп: Х>0,7хт и 0,5xm<2f^;0,7xm. Имеются следующие результаты испытаний:

Го =14, г'0=6, т 0=0, Гі=9, г',=5, т і = 6,

и

 

6

 

 

Ц

іок = 11-Юз ч,

£

і'0к =

8,4.10*4,

ft=l

Ä=1

 

9

5

 

 

 

S

i 1K=8,4.10a ,,

S

^ , =

3,3.10«

k=l

k=i

 

Воспользуемся формулами (7-1)— (7-3) и получим оценки ма­ ксимального правдоподобия и их вариации

?)*„ =5г0,822-Ю-з

ч- \

V. = 0 ,2 8 9 ,

Х*0 0,178-ІО -3

ч - \

Ѵ\

=^0,5,

 

 

л0

 

236


vf, ^0,56 2 -1 0 -5 ч ~ \

V

Ö,3?8,

А.*, ^ 0 ,5 2 - 1 0 - 3 ч - \

^ = ^ 0 ,5 7 8 .

Выберем достоверность у=0,9 и доверительные вероятности уі = =0,95, у2=0,05. Используем номограмму рис. 7-1 для получения интервальных оценок для интенсивностей. Получим:

0,573- 10-3^ П о ^ 1 ,27Юг3;

0,928- 10-3s£Aos:0,375- 10“ 3;

0,351 • 10-3sgrii:=Sl,06- ІО-3;

0,247 -10-3^Х і^і1,17 - ІО“ 3.

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

Таблица 7-1

Состояние

Время изменения

Дата

Примечание

устройства

 

состояния

н

7

ч 05 мин

' 14/ѵі— 1969

х = 138

дб

Уг

20

ч 15 мин

18/ѴІ—1969

 

 

Вг

9

ч 30 мин

20/ѴІІ— 1969

 

 

НИ

14 ч 21 мин

5/ІХ— 1969

х=136

дб

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

В табл. 7-1 дан пример заполнения этой формы применительно к приемнику радиолокационной станции, для которого определяющим параметром служит чувствительность.

Для удобства записи основных состояний устройства применя­

ются

следующие сокращенные условные обозначения: Н — нормаль­

ная

работа; Уі — работа с t-м

ухудшенным значением

выходного

параметра (выходной параметр

находится в t-м интервале кванто­

вания); Ві — выполнение і-й операции

восстановления; ГК — выпол­

нение і'-й операции профилактического

обслуживания;

НИ — неис­

пользование. Сюда же в случае

необходимости можно

внести Т

транспортировка; X — хранение и др.

При заполнении формы необходимо использовать оптимальную периодичность контроля работоспособности устройства в зависимости

2 37


от дисперсии и интенсивности ухудшения выходных параметров (см. § 4-2). Так как записи в журнале учета статистики производят в ди­ скретные моменты времени контроля работоспособности или измене­ ния состояния, то исходные статистические данные получают при минимальных затратах времени обслуживающего персонала.

В этом параграфе рассмотрены особенности оценки параметров экспоненциальных законов распределений для марковских однород­ ных моделей. Особенностью оценки параметров неоднородных моде­ лей является необходимость решения не одного, а системы нелиней­ ных уравнений для определения нескольких параметров (см. § 7-3).

7-3. Оценка параметров математических моделей в условиях неполной определенности

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

Обозначим:

N — число испытываемых (эксплуатируемых) эле­

ментов (объем

партии); R(t) — число элементов, работоспособных

к моменту времени t\ d(t) — число элементов, отказавших к момен­

ту времени t.

Рассмотрим два основных вида неопределенности; неопределен­

ность

первого вида заключается в том, что N является случайной

величиной, а

R(t) и d(t) — случайными

процессами;

неопределен­

ность

второго

вида — в том, что известны

реализации

только неко­

торых из этих величин и процессов, например реализации d(t). По таким статистическим данным необходимо оценить интенсивность отказов Л(£), а также, если требуется, то и характеристики остав­ шихся неизвестных из числа /V, R (t) и d (t).

Ясно, что основную информацию о зависимости интенсивности отказов от времени несут процессы R(t) и d(t), поэтому целесооб­ разно рассматривать четыре основных типа задач. К первому типу относятся задачи, в которых известны реализации N и R(t), ко вто­

рому— задачи,

в которых известны реализации N и d(t),

к третье­

м у— задачи, в

которых известны только реализации R{t),

и, нако­

нец, к четвертому — задачи, в которых известны только реализации

d(t). В первых двух типах задач требуется определить характери­

стики

A(t) — моменты

и

интервальные оценки, в

следующих

двух — характеристики Л(£)

и N. На

практике

часто встречаются

задачи

. четвертого типа

,

решение

которых

является

наиболее

сложным.

 

 

 

 

 

Так как приближенное решение задач первых трех типов до­ вольно просто и требует привлечения всего лишь метода линеариза­ ции и нормальной аппроксимации одномерного распределения Л(£), мы остановимся на задаче четвертого типа.

Пусть реализации d(t) известны, требуется определить характе­ ристики A(t) и N. Так как реализации d(t) имеют разрывы первого рода, целесообразно использовать квантование времени и решение

искать на временной решетке £і = £А£, і=О, I, Л£=Т/І, где Т — про­ должительность интервала наблюдения. Для упрощения записи обо­ значим d(iAi)—di, A(iAt)=Ai. Величины di играют роль порядко­ вых статистик.

238


Р а с с м о т р и м д в а с о с е д н и х м о м е н т а в р е м е н и U и

д л я

них

 

Ш

I

— j A(jc) dx

 

dt = N \ \ - e

0

(7-9)

 

d j — d t

(7-10)

Л*=- {N-di)Lt '

Формула (7-10) для интенсивности отказов

является известной,

ее отличие в том, что di, dj и N являются случайными величинами.

Оодставляя в (7-10) значение di

из формулы '(7-9),

после пре­

образований получим:

Ш

 

 

 

 

 

 

г

у - И ) «

- I

_____

 

vt +

j Л (x)dx

е U

— vt = 0, i = l . l .

(7-11)

_

ibt

 

 

 

где vt = (di+1 — di) dt 1 .

Трансцендентное уравнение (7-11) является основным для опре­ деления характеристик интенсивности в том случае, когда известны только характеристики числа отказавших элементов. Так как полу­

чить

решение

выражения

(7-11) в аналитическом виде невозможно,

то приходится

применять

статистический подход — по

реализациям

d(t)

получать

реализации

A(f), а затем стандартными

методами их

обрабатывать. Отсюда же следует, что определение Л(/) для произ­ вольного закона распределения является серьезной и сложной за­ дачей.

Часто предполагают, что распределение име'ет удобный одноили двухпараметрический вид; например, принадлежит к семействам гамма-распределений, нормальных распределений, распределений

Вейбулла — Гнеденко и

др. Стандартные

методы точечной

оценки

параметров этих распределений — метод

максимального правдопо­

добия, метод моментов,

метод квантилей

и др. — в данном

случае

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

Иногда удобно отказаться от привычных распределений и при­ менять квазидетерминироваиные представления

п

п

л (0 =

11 A ttt, A (0 =

<?t_ °

.

(7-12)

 

;=о

 

 

 

где А і случайны. Такие

подходы развиваются в работах [Л.

64, 87].

Мы остановимся на

самых простых

задачах, в которых

с точ­

ностью-до параметров известен аналитический вид Л (аь аі, t), т. е. известен характер предполагаемого распределения. Тогда реализации

параметров а*, а і определяют из решения системы уравнений (7-11),

239