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