CC BY
61
УДК 519.2
DOI: 10.17586/0021-3454-2016-59-3-197-201
ОЦЕНИВАНИЕ ХАРАКТЕРИСТИК НАДЕЖНОСТИ СЛОЖНОЙ СИСТЕМЫ ПРИ ОГРАНИЧЕННОМ ОБЪЕМЕ ЭКСПЕРИМЕНТАЛЬНЫХ ДАННЫХ
А. А. Ардашов, В. Н. Арсеньев, С. Б.Силантьев
Военно-космическая академия им. А. Ф. Можайского, 197198, Санкт-Петербург, Россия
E-mail: [email protected]
Рассматривается задача оценивания характеристик надежности сложной системы по ограниченному объему опытных данных. Для повышения точности оценок предлагается использовать метод приоритета опытной информации, позволяющий объединить все имеющиеся априорные и экспериментальные данные о характеристиках надежности.
Ключевые слова: сложная система, надежность, ограниченные испытания, метод приоритета опытной информации
Введение. Надежность является одним из основных свойств сложной системы, характеризующим ее способность выполнять целевые задачи в заданных режимах и условиях применения. Поэтому техническое задание на разработку системы всегда содержит требования к ее надежности. Для проверки соответствия сложной системы требованиям заказчика на всех этапах ее построения проводятся всесторонние исследования надежности. Получаемая при этом информация может быть условно разделена на априорную (накопленную до проведения испытаний опытных образцов системы) и опытную (полученную по результатам натурных испытаний). Число натурных испытаний, как правило, весьма ограничено, что не позволяет обеспечить высокую достоверность оценок характеристик надежности. Для повышения качества оценивания используется априорная информация.
В настоящее время круг методов, позволяющих производить совместную обработку априорных и опытных данных, достаточно ограничен [1—5]. Сложность применения таких методов связана главным образом либо с необходимостью выбора априорных распределений (методы, основанные на использовании формулы Байеса), либо с неопределенностью выбора весов априорной и опытной информации в апостериорных оценках характеристик надежности. Для решения данной задачи предлагается использовать метод приоритета опытной информации (ПОИ), который лишен этих ограничений [6, 7].
Постановка задачи. Полагается, что время безотказной работы сложной системы (слу-
л
чайное время наработки системы до отказа) 0 < X имеет гамма-распределение
Ф л (X; ц) = ф л (X; а, р) = 1 Xа exp {- X/р}, (1)
X X 1(а + 1) Р
T
где ц = [а Р] — вектор параметров распределения; а>-1; Р> 0; Г() — гамма-функция.
Параметры распределения (1) связаны с математическим ожиданием (средним временем
л
безотказной работы) M л и дисперсией D л величины X зависимостями
X X
Mл =(а + 1)Р, Dл =(а + 1)р2. (2)
X X
По результатам натурных испытаний Nc опытных образцов получены значения времени наработки системы до отказа: Xi, i = 1, Nc.
Имеется априорная информация, представленная оценками Mр и Dp математического ожидания MЛ и дисперсии DЛ , либо оценками а p , ß p параметров а, ß распределения (1).
Оценки ар, ßр и Mр, Dp связаны между собой соотношениями, аналогичными выраже-
Г ~Т
ниям (2), цр = [ар ßр .
Необходимо получить апостериорные оценки вероятностных характеристик времени
Л
безотказной работы X системы, учитывающие всю имеющуюся априорную и опытную информацию.
В соответствии с методом приоритета опытной информации на первым этапе решения задачи необходимо по опытным данным найти оценки максимального правдоподобия параметров распределения (1).
Определение опытных оценок. Функция правдоподобия, соответствующая распределению (1) и выборке X-, i = 1, Nc , имеет вид
N Nc 1 f Nc Ла f 1 Nc f
ПфЛ^;м) = ПФл (X;а,ß) = гN Nc(а„) ПX expf--Ix f. (3)
i=1 x i=1 x г c (а + 1)ß cV \ i=1 ) [ p i=i J
Оценки максимального правдоподобия а0, ßo параметров а, ß обеспечивают максимум функции (3), т. е.
Nc Nc
(а0,ß0) = arg max Пфл (Xi;а,ß) = arg max lnПфл (X;а,ß).
а>-1, ß>0 j=1 X а>-1, ß>0 X
Поскольку функция плотности распределения (1) является регулярной, то определение опытных оценок а0 и ß0 (или = [а0 ß0 ] ) осуществляется при необходимых условиях максимума функции (3) по параметрам а и ß :
Nc Nc
д ln Пфл (Х-; а, ß) д ln Пфл (X-; а, ß)
i=1 X
да
= 0;
i=1
X
а=а0 ,ß=ß0 dß
= 0:
а=а0 ,ß=ß0
что приводит к системе из двух нелинейных уравнении с двумя неизвестными:
1 "с 1 "с
( +1) = — 2X; У(а0 +1) + 1пР0 = — £ 1пX,, (4)
"о ,=1 "с ,=1
где ¥(•) — пси-функция.
Решение системы уравнений (4) может быть получено методом Ньютона. На основе опытных оценок а0 и (0 параметров а и ( в соответствии со вторым этапом метода ПОИ определяется вес априорной информации в апостериорных оценках и сами оценки.
Определение апостериорных оценок характеристик надежности. Функция правдоподобия (3) с учетом (4) может быть представлена в виде
Цц0;ц) =-1 АГ( ^ ехр|а"с Г^(а0 +1) + 1п(01-"с (ао +1)(о I.
Г "с (а +1)( "с (а+1) Р1 0 Г ^ ^ Ро 1 ( /
Поскольку близость результатов априорных оценок к опытным данным характеризуется отношением правдоподобия, для проверки гипотезы Н : а = ар, ( = (р :
V =■
Ц(Ро; р р)
о;Ро)
^ (ао + l)JЗ^c («о +1)
exp
N (ар - «о))Т («о +1) + 1п Ро ] + N (ао +1)
' Р ^
1 ^о
Р
р У
■ (а р + 1)Р N(а р +1)
в качестве величины, определяющей вес априорных данных в апостериорных оценках, ис-
*
пользуется значение V , которое может изменяться в диапазоне от 0 до 1.
Полагается, что опытные оценки ао и Ро параметров а и Р получены методом макси*
мального правдоподобия по некоторой гипотетической выборке объемом Nр =v из генеральной совокупности с распределением (1). Соответствующая функция правдоподобия имеет вид
Np (ар + 1)Рр
Ц (р; р )
1
ГNP (а + 1)Р
Np (а+1)
exp
а^
Т(а р + 1) + 1п Р;
Р
Тогда апостериорные оценки аа, Ра параметров а, Р определяются исходя из необхо-
димых условий максимума функции Ц р) = Ц(ро; р)р; р):
= 0, которые
Р=[аа Ра
приводят к системе уравнений
Т(а а + 1) + 1п Ра =
Т (ао +1) + 1п Ро + V* Гт (а р +1) + 1п Р;
(аа + 1)Ра =
1 + V
(«о + 1)Ро +v* (а р + 1)Р р
*
1 + V
(5)
Апостериорные оценки среднего времени безотказной работы системы и дисперсии, учитывающие результаты натурных испытаний и априорных исследований надежности, определяются по формулам, аналогичным (2):
Ма = (а а + 1)Ра =
Мо + V Мр
, = МаРа =
М2
(6)
1+ У "" " аа +1
Можно показать, что дисперсии оценок (6) меньше дисперсий соответствующих оценок ао и Ро, полученных без учета априорной информации. В частности, приближенный выигрыш в точности оценивания математического ожидания М Л , получаемый за счет учета ап-
х
риорной информации, определяется как
В[Ма] Оо (1 + ^)2
5 = -
(7)
О [Ма ] В0 +V*Dp '
откуда следует, что максимальный выигрыш 5 = 2 достигается при полном совпадении апри-
*
орных и опытных данных, поскольку в этом случае V = 1.
Пример. Рассмотрим сложную систему. Известно, что длительность интервала времени
Л
х между первым и последующими отказами системы имеет гамма-распределение.
Известны результаты испытаний 10 опытных образцов (Nc = 10): X = 1096 ч, Х2 = 981 ч, Х3 = 766 ч, Х4 = 1209 ч, Х5 = 993 ч, Х6 = 1240 ч, Х7 = 1166 ч, Х8 = 816 ч ;
2
Х9 = 718 ч; Хю = 924 ч, а также априорные оценки Мр = 1100 ч и Бр = 40000 ч математи-
л
ческого ожидания М л и дисперсии Б л величины X, на основе которых по формулам, ана-
X X
логичным (2), рассчитаны априорные оценки параметров распределения (1): а р = 29,250; (р = 36,364 ч .
Необходимо найти апостериорные оценки характеристик надежности системы. Оценки максимального правдоподобия параметров а, (, математического ожидания
л
М л и дисперсии Б л величины X, полученные по результатам испытаний опытных образ-
X X
цов, равны: а0 = 29,436; (0 = 32,554 ч, М0 =991 ч и Б0 =32256 ч2 .
Отношение правдоподобия V = 0,203 . Тогда апостериорные оценки характеристик надежности системы, полученные по формулам (5), (6), имеют следующие значения: аа = 28,037 ; (а = 34,756 ч ; Ма = 1009 ч; Ба = 35076 ч . Отсюда видно, что априорная оценка среднего времени безотказной работы системы является достаточно завышенной.
В соответствии с показателем (7) апостериорная оценка Ма точнее опытной оценки
М0 примерно в 1,16 раза.
Заключение. Применение метода приоритета опытной информации для оценивания характеристик надежности сложных систем по результатам испытаний ограниченного числа образцов позволяет улучшить качество оценивания этих характеристик и повысить достоверность принимаемых решений о соответствии системы требованиям технического задания.
СПИСОК ЛИТЕРАТУРЫ
1. Пугачев В. Н. Комбинированные методы определения вероятностных характеристик. М.: Сов. радио, 1973. 256 с.
2. Элементы теории испытаний и контроля технических систем / Под ред. Р. М. Юсупова. Л.: Энергия, 1978. 192 с.
3. Кринецкий Е. И., Александровская Л. И., Шаронов А. В., Голубков А.С. Летные испытания ракет и космических аппаратов. М.: Машиностроение, 1979. 464 с.
4. Щербаков П. С. Использование априорной информации для уточнения оценок параметров // Изв. АН СССР. Автоматика и телемеханика. 1988. № 5. С. 80—89.
5. Арсеньев В. Н. Новые методы принятия решений при ограниченных экспериментальных данных. СПб: ВКА им. А. Ф. Можайского, 1999. 90 с.
6. Арсеньев В. Н., Ардашов А. А., Зайцев М. А. Оценивание характеристик надежности сложных систем по результатам испытаний ограниченного числа опытных образцов // Сб. статей II Всерос. науч.-практ. конф. „Современные проблемы создания и эксплуатации вооружения, военной и специальной техники". СПб: ВКА им. А. Ф. Можайского, 2014. С. 190—191.
7. Арсеньев В. Н. Оценивание характеристик систем управления по ограниченному числу натурных испытаний. М.: Рестарт, 2013. 126 с.
Август Анатольевич Ардашов Владимир Николаевич Арсеньев
Сергей Борисович Силантьев
Сведения об авторах
— канд. техн. наук, доцент; ВКА им. А. Ф. Можайского
— д-р техн. наук, профессор; ВКА им. А. Ф. Можайского, кафедра бортовых информационных и измерительных комплексов;
E-mail: [email protected]
— канд. техн. наук, доцент; ВКА им. А. Ф. Можайского, кафедра автономных систем управления
Рекомендована кафедрой Поступила в редакцию
бортовых информационных 01.06.15 г.
и измерительных комплексов
Ссылка для цитирования: Ардашов А. А., Арсеньев В. Н., Силантьев С. Б. Оценивание характеристик надежности сложной системы при ограниченном объеме экспериментальных данных // Изв. вузов. Приборостроение. 2016. Т. 59, № 3. С. 197—201.
ASSESSMENT OF RELIABILITY CHARACTERISTICS OF A COMPLEX SYSTEM BY A LIMITED VOLUME OF EXPERIMENTAL DATA
A. A. Ardashov, V. N. Arseniev, S. B. Silantev
A. F. Mozhaysky Military Space Academy, 197198, St. Petersburg, Russia E-mail: [email protected]
The problem of evaluation of complex system reliability characteristics by a limited volume of experimental data is considered. To improve the estimates accuracy, it is proposed to apply the method of experimental data priority allowing for combining all available a priory and experimental data on the reliability characteristics.
Keywords: complex system, reliability, limited testing, method of experimental data priority
Data on authors
Avgust A. Ardashov — PhD, Associate Professor; A. F. Mozhaysky Military Space Academy Vladimir N. Arseniev — Dr. Sci., Professor; A. F. Mozhaysky Military Space Academy, Department
of Onboard Information and Measuring Complexes; E-mail: [email protected]
Sergey B. Silantev — PhD, Associate Professor; A. F. Mozhaysky Military Space Academy, Department of Independent Control Systems
For citation: Ardashov A. A, Arseniev V. N., Silantev S. B. Assessment of reliability characteristics of a complex system by a limited volume of experimental data // Izv. vuzov. Priborostroenie. 2016. Vol. 59, N 3. P. 197—201 (in Russian).
DOI: 10.17586/0021-3454-2016-59-3-197-201
