CC BY
45
пределенная система, охватывающая пункты управления (поисково-спасательные посты) различных уровней и аварийно-спасательные формирования в соответствии с зонами их ответственности, и обеспечивает информационно-техническое взаимодействие с поисково-спасательными службами других министерств и ведомств, в первую очередь ГМСКЦ Минтранса и МЧС России.
Литература
1. Организация аварийно-спасательного обеспечения на море // Морской флот. № 1. 2007. С. 4-12.
2. Положение о взаимодействии аварийно-спасательных служб министерств, ведомств и организаций на море и в водных бассейнах России // Российские вести. 1995. № 179. С. 5-6.
3. Справочник спасателя. Кн. 8: Надводные и подводные спасательные работы. М.: ФЦ ВНИИ ГОЧС, 2006. 204 с.
4. Абчук В.А., Суздаль В.Г. Поиск объектов. М.: Сов. радио, 1977. 336 с.
References
1. Organizatsiya avariyno-spasatelnogo obespecheniay na more [Organization emergency rescue at sea ensure], Morskoy flot [Marine fleet], 2007, no. 1, pp. 4-12.
2. Polozhenie o vzaimodeystvii avariyno-spasatelnykh sluzhb ministerstv, vedomstv i organizatsy na more i v vodnykh basseynakh Rossii [The position of the interaction of the emergency services of the ministries, departments and organizations at sea and in the river basins of Russia], Rossiyskie vesti [Russian Vesti], 1995, no. 179, pp. 5-6.
3. Spravochnik spasatelya. Kniga 8. Nadvodnye i podvodnye spasatelnye raboty [Rescuers Guide. Book 8. Abovewater and underwater rescue], Moscow, All-Russian Research Institute for Civil Defense and Emergencies, 2006, 204 p.
4. Abchuk V.A., Suzdal V.G., Poisk obyektov [Object search], Moscow, Sov. radio, 1977, 336 p.
УДК 519.718
МЕТОД ОБЕСПЕЧЕНИЯ ФУНКЦИОНАЛЬНОЙ НАДЕЖНОСТИ КОМПЛЕКСОВ СРЕДСТВ АВТОМАТИЗАЦИИ НА ОСНОВЕ ПРИМЕНЕНИЯ ИЗБЫТОЧНОСТИ РАЗЛИЧНОГО ВИДА
В.Л. Лясковский, д.т.н., профессор (Военная академия воздушно-космической обороны им. Маршала Советского Союза Г.К. Жукова, ул. Жигарева, 50, г. Тверь, 170022, Россия, l_vik_l@mail.ru); Р.В. Допира, д.т.н., зав. отделением (НИИ «Центрпрограммсистем»,
просп. 50лет Октября, 3а, г. Тверь, Россия, rvdopira@yandex.ru); А.Ю. Кабардинский, начальник кафедры; А.А. Догадов, инженер (Московский государственный технический университет им. Н.Э. Баумана, ул. 2-я Бауманская, 5, г. Москва, 105005, Россия, kabardinsky@pmstu.ru, ToxaDog@gmail.com);
А.А. Беглецов, ст. инженер (Военное представительство МО РФ, ул. 1-я Бухвостова, 12/11, г. Москва, 107076, Россия)
Предложен метод обеспечения требуемого уровня функциональной надежности комплексов средств автоматизации систем управления радиоэлектронными средствами на основе применения избыточности различного вида. В качестве показателя надежности рассмотрена вероятность достоверного решения задачи. Предложена формула для расчета вероятности достоверного решения задачи, учитывающая потоки случайных сбоев и отказов. Разработана целевая функция. Для применения методов бивалентного программирования целевая функция преобразована в форму, содержащую булевые переменные. Рассмотрены и реализованы в среде MATLAB точные и приближенные методы решения задачи выбора используемых средств и методов избыточности различного вида. Получены формулы для оценки среднего времени решения задачи методом полного перебора и эвристическими методами в зависимости от размерности задачи. Приведены оценки точности эвристических методов.
Ключевые слова: комплексы средств автоматизации, функциональная надежность, средства резервирования, методы временной избыточности, решение оптимизационных задач, эвристические методы.
A REDUNDANCY-BASED FUNCTIONAL RELIABILITY METHOD FOR AUTOMATION FACILITIES SETS
Lyaskovsky V.L., Ph.D., professor (Military Academy of the Aerospace Defence, Zhigareva St., 50, Tver, 170022, Russia, l_vik_l@mail.ru);
Dopira R. V., Ph.D., head of department (R&D Institute «Centrprogrammsystem», 50 let Oktyabrya Av., 3a, Tver, 170024, Russia, rvdopira@yandex.ru); Kabardinsky A. Yu., chief of department; Dogadov A.A., engineer (Bauman Moscow State Technical University, 2ndBaumanskaya St., 5, Moscow, 105005, Russia, kabardinsky@bmstu.ru, ToxaDog@gmail.com);
Begletsov A.A., senior engineer (Military Representation of the Ministry of Defence of Russia, 1st Bukhvostova St., 12/11, Moscow, 107076, Russia) Abstract. A redundancy-based method providing a required level of functional reliability of radio-electronic equipment control systems is given. The reliable task solution probability is chosen as the reliability index. The formula for calculating the probability of the reliable task solution is proposed. It considers soft errors and hard faults. The redundancy optimization
problem is formulated. An objective function for the optimization task is designed and converted to the form suitable for the bivalent programming. It includes boolean variables. Exact and approximate methods of optimization are realized in MATLAB. The formulas for calculating average time of the optimization using exhaustive and heuristic methods are given. Accuracy estimations of the heuristic methods are given.
Keywords: automation facilities set, functional reliability, backup, redundancy, transient error, breakdown, optimization, greedy algorithms, heuristic, brute-force search.
Одной из основных характеристик комплексов средств автоматизации (КСА) систем управления радиоэлектронными средствами является их функциональная надежность - способность выполнять требуемые функции в заданных режимах и условиях применения [1]. Актуальность решения задачи обеспечения надежности данной КСА связана с наличием большого количества задач, требующих высокой вероятности достоверного решения.
Цель работы - разработка метода обеспечения функциональной надежности КСА систем управления радиоэлектронными средствами на основе применения средств и методов избыточности различного вида.
Рассматривается задача обеспечения функциональной надежности КСА систем управления радиоэлектронными средствами на основе применения избыточности различного вида. Она заключается в выборе для каждого вычислительного модуля (ВМ) таких средств структурной избыточности (резервирования) и методов временной избыточности (программно-логических), с помощью которых для всех типов заявок, обслуживаемых КСА, вероятность достоверного решения была бы не ниже заданной, а суммарная стоимость применения средств и методов минимальной.
При постановке задачи зададим следующее: КСА состоит из множества ВМ, А={а1} (¡=1, ..., п) с производительностями 0={^Д, стоимостями С={с} и предназначен для обслуживания т типов заявок В={Ь¡}, (¡=1, ..., т), для которых задан вектор среднего объема вычислений 0={о,} и требуемая вероятность достоверного решения р*={р*}. Будем считать, что задачи, соответствующие заявкам, строго распределены по ВМ и существует матрица распределения задач по ВМ S={sy■}, где SiJ■=1, если задача по заявке /-го типа решается в ¡-м ВМ, иначе SiJ■=0.
Массогабаритные и энергетические характеристики КСА ограничены величинамиМ* и Я* соответственно. В зависимости от предъявляемых требований множество значений массогабаритных характеристик ВМ М={М¡} представляет собой вектор со значениями масс или объемов ВМ. Если на проектируемую систему наложены ограничения и по массе, и по объему, то вектор М={М¡} целесообразно разбить на вектор масс и вектор объемов ВМ. Множество значений энергетических характеристик Я={Я^ является вектором со значениями мощности, потребляемой от сети.
КСА характеризуется также потоками случайных отказов L0={X0j} и сбоев ЬС={ХС}} (рассматри-
ваются только критические сбои, которые не обнаруживаются системой контроля КСА и приводят к ошибкам при решении программ в КСА). Время восстановления работоспособного состояния ВМ задано вектором Т= {$.
При резервировании параллельно к основным ВМ а■ могут подключаться дополнительные ВМ в режиме нагруженного резервирования.
Все методы временной избыточности (программно-логические методы) повышения достоверности решения программ в КСА вызывают снижение интенсивности сбоев в Б={ё1} раз (1=1, ..., ф и увеличение времени выполнения задачи в Е={Е^} раз, а также характеризуются значениями стоимости реализации Срш = {с™ }.
Результатом решения поставленной задачи оптимизации являются вектор У={у} и матрица X={xft}, где у=т, если ¡-й ВМ резервируется т дополнительными, у=0 в противном случае; хц=1, если 1-й программно-логический метод применяется для комплекса задач по заявке /-го типа, хд=0 в противном случае.
Вектор У и матрица X ищутся такими, чтобы суммарная стоимость обеспечения функциональной надежности КСА была минимальной с учетом ограничений на массогабаритные и энергетические характеристики, а вероятность достоверного решения для каждой задачи не меньше заданной.
Формально для заданных исходных данных А, О, С, В, О, р*, £ М*, Я*, М, Я, Lo, Lc, Т, Б, Е, С™ требуется найти такие У', X', для которых
( п т д \
(7X') = argmin ^
ty
V j=1
j +ZZ<
=1 l=1
(1)
при ограничениях Р^У, X)>p*i, М(У)<М*, Я(У)<Я*, где Р,(У, X) - вероятность достоверного решения задачи по заявке /-го типа; при условии, что заявки /-го типа выполняются только в ¡-м ВМ, Р,(У, Х)=Ру(У, X); М(У) - суммарные массогабаритные характеристики КСА при применении резервирования; Я(У) - суммарные энергетические характеристики КСА при применении резервирования.
При условии, что заявки /-го типа выполняются только в ¡-м ВМ, вероятность достоверного решения задачи по заявке /-го типа [2] будем рассчитывать следующим образом:
( / / чЛЛ+Л
pj =
\ (
1 -I 1 -K. exp
ПЕ,))У
х exp
Е
(2)
tj П E
l у
где H - множество методов временной избыточ-
ности, применяемых для задачи по заявке /-го типа; Щ - коэффициент готовности ¡-го ВМ, 1/
^ = V
ti
среднее время выполнения
заявки /-го типа в ¡-м ВМ, t.. = —L.
В формуле (2) произведения по множеству Н могут быть вычислены с использованием компонент матрицы X:
П Е, = exp I Z хи ■ In (Е,)
П-т- = exp
Z х •ln
( T? \\
d,
V l //
(3)
(4)
Преобразуем исходную целевую функцию (1) в эквивалентную с булевыми переменными для применения методов решения, относящихся к задачам бивалентного программирования [2]. Так как у, может принимать значение 0 или целых положительных чисел, то у) можно свести к сумме двоичных переменных, разложив каждую из переменных у] по степеням 2:
у, = м,1 + 2п,2 + 22п,3 +... + 2к-1п^. (5)
Число 2к-1 является верхней границей для у, для ¡-го ВМ к находится из условия
2k -1 < min
f Л
M R
M R
V _ J _ _ j _ /
(6)
где [•] - ближайшее меньшее целое число.
Таким образом, целевую функцию можно преобразовать:
( п т д Л
(7X') = ш^т £с,у, с
¥-х I,=1 и
= argmin(c1n11 +... + 2к-1 с^ +...
и ,х
сРшх„
+cnunl +... + 2k-1 спил + сРГхп +... + cPLM
(7)
х )
qm qm '
при ограничениях Pi(U, X)>p*i, M(U)<M*, R(U)<R*, где U - вектор, полученный при представлении вектора Y={yj} в виде набора двоичных переменных.
Размерность приведенной задачи (7) определяется следующим образом:
w=mq+nk. (8)
Поставленная задача относится к задачам дискретной оптимизации. Методами ее точного решения являются, например, метод полного перебора и метод ветвей и границ [3]. Метод ветвей и границ, реализованный в пакете Optimization Toolbox среды MATLAB, позволяет решать задачи бивалентного программирования с одним линейным ограничением. Задача подобного типа рассмотрена в [4], где было установлено, что для компьюте-
ра AMD Athlon 64x2 Dual Core Processor 5000+ 2.61 ГГц, 2.00 Гб ОЗУ при размерности задачи w=200 не удается получить решение из-за переполнения памяти. Поэтому были исследованы два метода решения задачи (7): метод полного перебора и эвристические жадные методы (ЖМ). Поиск решений был проведен в среде MATLAB на компьютере AMD Athlon™ 64x2 Dual Core Processor 5000+ 2.61 ГГц, 2.00 Гб ОЗУ.
Решение задачи методом полного перебора предполагает последовательное нахождение целевой функции для 2w вариантов и выбор из всех полученных значений целевой функции экстремального.
В пакете Curve Fitting Toolbox среды MATLAB c помощью метода наименьших квадратов на интервале we [10; 20] с доверительной вероятностью 95 % было установлено, что экспериментально полученная зависимость среднего времени оптимизации методом полного перебора от размерности задачи (7) w описывается функцией
tex = a*2bw , (9)
где tex - среднее время оптимизации методом полного перебора (в сек.); ae[1,49-105; 2,60-10-5); be [0,96; 1,00].
Среднее время оптимизации для каждой размерности рассчитывалось по 20 измерениям.
Для решения поставленной задачи оптимизации можно предложить четыре ЖМ, основу каждого из которых составляют следующие эвристические правила выбора решения.
1. На каждом шаге выбор решения с наименьшим значением стоимости реализации.
2. На каждом шаге выбор решения, дающего наибольший прирост вероятности выполнения заявки.
3. На каждом шаге выбор решения с наибольшим отношением прироста вероятности выполнения заявки к стоимости реализации.
4. Комплексный ЖМ. Выбор лучшего решения из трех, полученных вышеперечисленными методами.
Для зависимости среднего времени решения задачи вторым ЖМ от размерности задачи w в пакете Curve Fitting Toolbox среды MATLAB c помощью метода наименьших квадратов на интервале we [10; 800] с доверительной вероятностью 95 % была найдена аппроксимирующая функция
tgr = a * w + b , (10)
где tgr - среднее время оптимизации ЖМ (в сек.);
ае[6,23-105; 6,66-105); Ье[-6,07-10-5; 1,71-Ю-4].
Среднее время оптимизации для каждой размерности рассчитывалось по 100 измерениям.
Экстраполяция полученных временных зависимостей на интервал ^е[20; 100] позволяет построить график зависимости времени выполнения задачи (7) / от размерности задачи w (см. рис.). (По оси / масштаб логарифмический.)
0
w
60
га
1E+25 1E+23 1E+21 1E+19 1E+17 о 1E+15 * 1E+13
(D
- 1E+11 1E+09 10000000 100000 1000 10 0,1 0,001 J
• Полный перебор----Второй "жадный" метод
В таблице приведено среднее значение точности ЖМ для решения задачи (7) с размерностью ^=20. Под точностью понимается отношение
-
W = *100%,
(11)
где - экстремальное значение целевой функции (7), полученное ЖМ; 2а - экстремальное значение целевой функции (7), полученное точным методом.
ЖМ Средний коэффициент точности W, % Количество задач, для которых не было найдено решение
Первый 74,6 7
Второй 48,9 0
Третий 90,0 9
Комплексный 96,7 0
Для каждого ЖМ проводилось 100 испытаний. Испытанием являлось решение задачи (7), составленной по первоначальной задаче (1), с п=1, т=1, Х°1=ХС1 =0,001 с-1, /ц =¿1 = 1 с. Векторы стоимостей, массогабаритных и энергетических характеристик задавались случайным образом.
На основании проведенных экспериментальных исследований установлено, что комплексный ЖМ обладает наибольшей точностью.
Таким образом, предложен метод обеспечения функциональной надежности КСА систем управления радиоэлектронными средствами на основе применения средств и методов избыточности различного вида. Для задач выбора средств избыточности, имеющих большую размерность, целесообразно применять, например, предложенный комплексный ЖМ или другие эвристические методы.
Литература
1. Малафеев С.И., Копейкин А.И. Надежность технических систем. Примеры и задачи. СПб: Лань, 2012. 320 с.
2. Коваленко К.А., Лясковский В.Л., Прохоров А.Г. К вопросу повышения надежности функционирования многомашинных вычислительных комплексов с использованием аппаратурных средств и программно-логических методов // Автоматика и телемеханика. 1997. № 3. С. 226-233.
3. Саати Т. Целочисленные методы оптимизации и связанные с ними экстремальные проблемы. М.: Мир, 1973. 304 с.
4. Лясковский В.Л., Догадов А.А. Об одном подходе к решению задачи обеспечения отказоустойчивости параллельных вычислительных систем // Технические и технологические системы: матер. IV междунар. науч. конф. Краснодар: Изд-во ВУНЦ ВВС ВВА, 2012. С. 243-247.
References
1. Malafeev S.I., Kopeykin A.I., Nadyozhnost tekhnicheskikh sistem. Primery i zadachi (Reliability of Technical Systems. Examples and Problems), St. Petersburg, Lan, 2012, 320 p.
2. Kovalenko K.A., Lyaskovsky V.L., Prokhorov A.G., Avtomatika i telemekhanika [Automatics and telemechanics], 1997, no. 3, pp. 226-233.
3. Saaty T., Optimization in Integers and Related Extremal Problems, McGraw-Hill, NY, 1970.
4. Lyaskovsky V.L., Dogadov A.A., Materialy 4 mezhdunar. nauch. konf. «Tekhnicheskie i tekhnologicheskie sistemy» (Proc. 4th Int. Sci. Conf. «Technic. and technolog. systems»), Krasnodar, Military and Training Research Center of the Air Force of the Military Air Acad., 2012, pp. 243-247.
20
40
80
100
120
УДК 004.652
ПРЕДСТАВЛЕНИЕ ДАННЫХ ПРИ РАЗРАБОТКЕ ПРОГРАММНОГО ОБЕСПЕЧЕНИЯ ПО ВЕДЕНИЮ ИНФОРМАЦИИ О СОСТОЯНИИ ПОИСКОВО-СПАСАТЕЛЬНЫХ ТЕХНИЧЕСКИХ СРЕДСТВ
Е.В. Реут, инженер-программист (НИИ «Центрпрограммсистем», просп. 50 лет Октября, 3а, г. Тверь, 1 70024, Россия,
теиХекаХеппа. cps.tver@gm.ail. сот)
В статье рассматривается форма представления в реляционной БД структурированных данных в виде древовидной зависимости с учетом специфики информации по состоянию поисково-спасательных технических средств. В данном случае древовидная зависимость имеет ограниченное количество уровней вложенности и различную сущность самих узлов дерева. Рассматриваются некоторые методы описания в реляционной БД структуры подобного типа с анализом каждого метода и рекомендациями по выбору формы представления такой структуры в реляцион-
