CC BY
36
ИНФОРМАЦИОННЫЕ ТЕХНОЛОГИИ И СИСТЕМЫ
УДК 681.326
В. И. Сеньченков, И. Н. Некрасов
ОГРАНИЧЕНИЯ В ЗАДАЧАХ ПОСТРОЕНИЯ ОПТИМАЛЬНЫХ АЛГОРИТМОВ ОПРЕДЕЛЕНИЯ ТЕХНИЧЕСКОГО СОСТОЯНИЯ СИСТЕМЫ
Рассматриваются особенности построения алгоритмов определения технического состояния системы, оптимальных по заранее заданному критерию. Предложены формальные конструкции, описывающие как последовательность выполнения проверок контролируемых признаков, так и ограничения при оптимизации алгоритмов. Представлена оригинальная математическая постановка задачи синтеза алгоритма, оптимального по критерию максимума средней вероятности получения правильного решения о техническом состоянии системы, со всеми видами ограничений.
Ключевые слова: алгебра подмножеств, алгоритм, критерий, ограничения, оптимизация, показатель эффективности, проверка, техническое состояние.
Введение. Важнейшей составной частью диагностического обеспечения системы является алгоритм определения технического состояния. В зависимости от подходов к обработке выходных процессов системы при формировании диагностических моделей принципы построения алгоритмов могут существенно различаться. В работе [1] представлены решающие правила, являющиеся основой алгоритмов, в случае выполнения формализации выходных процессов на основе свойств пространства 1_2 (пространство измеримых функций, квадратично интегрируемых по Лебегу [2]). Возникает задача формирования всех возможных вариантов последовательности выполнения проверок 7 (; = 1, п) контролируемых признаков (т. е. всех возможных вариантов реализации алгоритма определения технического состояния) с целью нахождения оптимального по некоторому заранее заданному критерию.
Постановка задачи определения последовательности проверок. Задано:
— множество изображений
Е = {Ег | I = 0т} (1)
всех видов технического состояния системы, сформированных в результате обучения [3], где
Е| = (вп, в,2, ..., вг, ..., вт)т; (2)
— распределение вероятностей
р ={Р(Е,)| I=0т} (3)
на множестве видов технического состояния;
— множество проверок
П = {71; | 7 = М }, (4)
при этом изображения Ег заданных видов технического состояния попарно различимы между собой, т. е. выполняется условие
V Ег, Е е Е, I Ф/3 щ-: щ е П, в,- Ф в^. (5)
Условие (5) означает, что для любой пары рассматриваемых видов технического состояния системы на заданном множестве П имеется хотя бы одна проверка, которая дает разные исходы в этих состояниях.
Координатами в,- изображения (2) являются модельные (типовые) исходы выполняемых проверок щ-, соответствующие 1-му виду технического состояния. Значения индекса I используются в дальнейшем следующим образом. При контроле работоспособности и поиске отказов системы значение I = 0 указывает на работоспособное состояние, I = 1, т — на отказы функциональных элементов. При контроле правильности функционирования системы I = 0 обозначает состояние неправильного функционирования, а I = 1, т — отдельные режимы правильного функционирования.
На основе исходных данных (1)—(5) требуется найти для каждого /'-го вида технического состояния системы упорядоченное подмножество проверок П ^ П такое, что
П = {щ- | щ е П, V Ег- е Е, V Е/е Е\{Ег}, е1 - Ф вГ]}. (6)
Из соотношения (6) видно, что каждое подмножество содержит такие щ-, на которых все виды технического состояния наблюдаемы, т.е. изображения попарно различимы между собой.
Формальное описание последовательности выполнения проверок. Очевидно, что множество
т
п* = у п i
I=0
содержит необходимые проверки для распознавания всех т +1 видов технического состояния. Определение технического состояния производится следующим образом: выполняется
*
первая выбранная проверка щ- е П , анализируются ее результаты и на основании этого вы-
*
бирается вторая проверка щ- е П , т.е. процесс накопления информации о состоянии системы является рекуррентным.
Следовательно, диагностирование, связанное с выполнением ряда проверок, исходы которых заранее непредсказуемы, необходимо рассматривать как случайный эксперимент [4, 5]. Универсальной математической моделью случайного эксперимента является вероятностное пространство. Применительно к рассматриваемому эксперименту его модель задается множествами
М = (Е, Ае, Р, 0), (7)
где Е — множество изображений (исходов эксперимента); Ае = {К | К ^ Е} — алгебра подмножеств множества Е, в которой элементы К имеют смысл фазовых состояний моделируе-
мого процесса; Р =
Р(К)Р(Р.) = 2 Р(Ег), КеА
Е.-е К
вероятностная мера, заданная на ал-
гебре Ае; 0 = {0г | 0*: Е ^ Р^, I = 0, т } — множество операторов, описывающих переходы от начального фазового состояния К = Е к конечным К; = {Ег}, каждое из которых содержит единственный элемент Е;, соответствующий идентифицированному виду технического состояния системы.
Таким образом, К; может интерпретироваться как принятое решение о техническом состоянии, т. е. как элемент множества R возможных решений о наблюдаемом техническом состоянии системы. Естественно, что и всякое фазовое состояние К ^ Е, с одной стороны,
Алгеброй называется система множеств, замкнутая относительно конечного числа операций пересечения этих множеств [2].
означает подмножество предполагаемых видов технического состояния, в одном из которых находится система. С другой стороны, фазовые состояния могут рассматриваться и как подмножества множества Я решений о техническом состоянии (К ^ Я).
До начала диагностирования неизвестно, какой вид технического состояния системы соответствует реальному, а поэтому начальное фазовое состояние совпадает с множеством Е, т. е. К = Е. При выполнении первой и последующих проверок контролируемых признаков из множества Е исключаются те Ег- е Е, которые несовместимы с исходами выполняемых проверок. Таким образом, после каждой проверки размерность исходного фазового состояния К = Е сокращается. Процесс продолжается до получения конечного фазового состояния К-, содержащего единственный элемент Ег-, фазовые состояния, содержащие два и более элементов, в дальнейшем обозначаются через Къ
Отдельная проверка 7- при ее г-м исходе переводит процесс диагностирования из некоторого фазового состояния К е Е в состояние Яу ^ Къ т.е. реализует отображение
7--: Нк ^ Нгк;, (8)
где = {Ег- |Ег- е ^, 7- =7;, г = 1, ш- }; ю- — число исходов проверки 7-.
Для достижения конечного фазового состояния К- в общем случае требуется несколько раз выполнить отображение (8), используя разные проверки щ е П, образующие подмножество Пг- ^ П проверок, необходимых для идентификации наблюдаемого состояния системы с одним из заданных видов технического состояния. С учетом этого всякий оператор 9г- из множества 0, входящего в модель (7), формально может быть описан в виде композиции отображений (8), каждое из которых реализуется одной проверкой:
е- = П (7; : К к ^ Щ;). (9)
7 е П-
Так как первая проверка 7- е П применяется в начальном фазовом состоянии Ку = Е, а последняя приводит к получению конечного фазового состояния Ну- = К-, то правая часть выражения (9) может быть представлена в виде
П :Е ^К , I = 0т. (10)
7Iе Неизвестно [6], что отображение (10) может быть реализовано, если
^ &шкег7; = р| -1, (11)
где | — мощность индексного множества видов технического состояния системы (| р | = т + 1); ё1шкег 7 - — размерность ядра отображения (10), которое определяется выражением
кег7- = {Е- | Е- е Ку, 7- =7г-}.
Равенство (11) является условием достижимости конечного фазового состояния из начального К = Е. Оно задает правило остановки процесса диагностирования.
Целенаправленным выбором проверок с учетом их характеристик можно гибко изменять алгоритм диагностирования, придавая ему требуемые свойства. В связи с этим необходимо выбрать характеристики, от которых зависят доверие к результатам диагностирования, его оперативность и экономичность.
Характеристиками доверия являются вероятности ошибок 1-го (а-) и 2-го (|37) рода при выполнении проверок [5, 7]. Ошибка 1-го рода возникает, когда измеренное значение контролируемого признака не соответствует допустимому интервалу его изменения, а в результате
проверки принимается, что значение входит в допустимый интервал. Ошибка 2-го рода появляется, когда измеренное значение контролируемого признака попадает в допустимый интервал, а в результате проверки принимается, что не попадает. Указанные вероятности образуют множества
которые включаются в модель диагностирования в качестве ее элементов.
Для выполнения проверок, а следовательно и для определения технического состояния системы, необходимы ресурсы. На этапе эксплуатации системы определяющую роль играют динамические ресурсы, прежде всего временной и трудовой, от наличия которых непосредственно зависит получение информации о техническом состоянии. Поэтому данные ресурсы и учитываются в предлагаемой модели.
Выполнение произвольной проверки требует времени на подготовку и установку контрольно-измерительного средства, подачу тестового воздействия, измерение значения контролируемого признака, обработку и интерпретацию результатов измерения. Кроме этого, требуются и подготовительные работы на системе, от которых также зависит продолжительность проверки. Например, при виброакустическом диагностировании отдельных узлов с вращательным или возвратно-поступательным движением инерционных масс установка датчиков должна производиться в непосредственной близости от этих узлов с целью уменьшения влияния помех. Для этого на системе выполняются определенные демонтажные работы.
Выполнение проверок сопровождается также трудозатратами. Если величина трудозатрат указывается только для того, чтобы производить сравнение проверок по данному показателю, то могут использоваться условные единицы, которые отражают соотношения трудозатрат при выполнении различных проверок. Такие единицы и используются в дальнейшем.
Таким образом, в качестве элементов модели рассматриваются множества
позволяет придавать необходимые свойства процессу принятия решений о техническом состоянии системы.
Формирование критерия оптимальности алгоритма определения технического состояния системы. Задача оптимизации алгоритма должна быть инвариантной по отношению к конкретному виду показателя эффективности Р, который применяется для ее решения.
Идентификация I-го вида технического состояния предполагает выполнение в определенном порядке некоторого подмножества П проверок, обладающих заданными характеристиками. Целенаправленным выбором проверок контролируемых признаков для всех подмножеств
П (I = 0,т) в зависимости от их характеристик и сформулированных ограничений можно добиться экстремального значения показателя Р.
Очевидно, что оптимизационная задача в данном случае относится к классу комбинаторных. Каждое упорядоченное подмножество П можно рассматривать как элемент алгебры подмножеств А(П), заданной на множестве П. Далее в рассмотрение необходимо включить
множество т[А(П)] перестановок на элементах алгебры А(П). Это множество содержит все возможные упорядоченные сочетания проверок контролируемых признаков. Следовательно, экстремальное значение показателя Р необходимо искать именно на этом множестве альтернатив.
А = {а- | - = 1,п}, В = {р- | - = 1,п},
(12)
(13)
Оптимизационная задача формулируется следующим образом. Найти последовательность (П/)/ е р упорядоченных подмножеств проверок контролируемых признаков такую, что
Р(П/)/ер = ейг р{((П/ )/е р )р}. (14)
ре т[А(П)]
Ограничения на состав контролируемых признаков включают условие попарной различимости изображений (6) и условие достижимости конечного фазового состояния (11).
В каждом частном случае структура равенства (14) определяется тем, какое конкретное свойство принимается в качестве показателя эффективности процесса диагностирования. Пусть таким свойством является достоверность решения о техническом состоянии, т.е. его соответствие наблюдаемому техническому состоянию системы.
Тогда в качестве показателя эффективности может быть выбрана вероятность Б получения правильного решения о техническом состоянии. Оптимизация должна проводиться по критерию максимума данного показателя:
Б(П /)/ер = тах Б{((П /)/е р )р}, (15)
} } ре т[А(П)] ■> ■> Р
при этом должны выполняться ограничения (6) и (11) на состав контролируемых признаков, а также на временные и трудовые ресурсы.
Пусть система диагностирования зафиксировала 1-й вид технического состояния при условии, что наблюдаемое состояние системы также соответствует /-му виду состояния, т.е. одновременно произошли события Ег и К,. Следует заметить, что события К не могут считаться достоверным результатам из-за наличия погрешностей при выполнении проверок. Очевидно, что условная вероятность Р(Кг | Ег) представляет собой вероятность получения правильного решения для /-го вида технического состояния: = Р(Кг | Ег). Так как при идентификации /-го вида технического состояния реализуются проверки из упорядоченного подмножества Пг, можно записать, что = БПг.
В этом случае вероятность Б получения правильного решения в целом для алгоритма определяется как средняя, которая аппроксимирует математическое ожидание величины Бг:
т
Б = ЕБ(П /)/ =0-е =2 Р(Е/ )Р(К |Е/). (17)
/=о
Структура ограничений оптимизационной задачи на ресурсы формируется из аналогичных рассуждений. Если 2 ^(л) и 2 с(К]) — суммарные затраты временных и трудо-
К]еЦ- К ]еПг
вых ресурсов на выполнение проверок контролируемых признаков, необходимых для определения /-го вида технического состояния, то оценку математического ожидания затрат получим из выражений:
тт ЕТ = 2Р(Е/) ^ ЦК]); ЕС = 2Р(Е,) 2 с(к,). (18)
г=0 к ]еПг г=0 К ]еПг
Затраты (18) не должны превышать максимально допустимых величин М и Мс соответственно.
Для того чтобы оптимизировать алгоритм определения технического состояния системы, необходимо оперировать выражением (17) в явном виде. Следовательно, должны быть выведены соотношения для вычисления вероятностей Р(Кг | Е/) как правильных решений о техническом состоянии при г = / так и ошибочных (/ Ф/).
Для идентификации произвольного вида технического состояния необходимо выполнить ряд проверок контролируемых признаков. При этом вероятности ошибок 1-го и 2-го рода зависят только от метрологических характеристик контрольно-измерительных приборов и
допусков на контролируемые признаки. Поэтому вероятности возникновения этих ошибок для разных проверок статистически независимы, что позволяет использовать формулу умножения вероятностей при вычислении общей ошибки контроля и диагностирования на множестве Пг-, исходя из этого
Р(К|Еу) = П у-у(7;). (19)
7Iе Пг
Здесь сомножители у-у (7-) определяются на основе решающих правил [1] и вероятностей
ошибок 1-го и 2-го рода производимых проверок. При контроле функционирования системы сомножители для различных целевых задач вычисляются из следующих выражений:
1 -а-, если у--в-1 = ш!п{[Уу-ву|}, I = 1, т; У =1, т
1 -Р;, если [у— | = ш1^{[угу-ву|}, 1 = 0;
У=0, т
а;,если -ву-| = Р1!2-{| У; -в-|}, у = 0;
У=1, т
Р;, если угу - вГ} | = ш]П{угу - ву; |}, ■ = 0, У=1, т
где — текущее значение ;-го контролируемого признака, для которого принадлежность 1-му виду технического состояния установлена.
При контроле работоспособности и поиске отказов выражения для определения сомножителей у- у (7 - ) имеют вид:
1 -а-, если{у---ву| = ш^У/ -ву |}, / = 0;
У=0, т
У| У(7 -) =
(20)
у- У(7 -) =
(21)
1 -р-, если У--е-1 = щт.{|у---ву1}, / = 1,т;
У=1, т
а -, если |у-- - в-1 = ш]п{угу - в-1}, - = 0;
У=1, т
Р-, если ]уи -^^| = ш1П-{[yгy -в-| }, У = 0.
У=1, т
Таким образом, с учетом соотношений (19)—(21), а также ограничений (6), (11) и (18) задача оптимизации по критерию максимума средней вероятности получения правильного решения о техническом состоянии системы математически формулируется следующим образом. Найти упорядоченные подмножества П;- ^ П ( - = 0, т ) такие, что:
Г \
ЕБ = шах
ре ш[Л(П)]
Е Р(Ег)
-=0
П У-г(7-)
V71еП-
Пi = {7- | 7- е П, V Е- е Е, V Еу е Е \{Ег}, в- - * ву-}; Е кег7- = -1;
71еПг
(22)
Е р(е- ) Е '(7- ) * М; е Р(Е-) Е с(7- ) * М.
-=0
7 Iе П
-=0
7}еП
Заключение. Решение задачи (22) обеспечивает максимальную достоверность определения технического состояния системы. Наиболее рациональным для решения этой задачи является метод динамического программирования [8].
список литературы
1. Сеньченков В. И. Решающие правила в алгоритмах определения технического состояния системы // Изв. вузов. Приборостроение. 2013. Т. 56, № 3. С. 5—11.
2. Колмогоров А. Н., Фомин С. В. Элементы теории функций и функционального анализа. М.: Физматлит, 2009. 572 с.
3. Сеньченков В. И. Модели, методы и алгоритмы анализа технического состояния. Saarbrücken: LAP LAMBERT Academic Publishing, 2013. 377 с.
4. Дмитриев А. К., Юсупов P. M. Идентификация и техническая диагностика. М.: МО, 1987. 521 с.
5. Контроль и диагностика сложных технических систем / Под ред. В. Г. Беликова. М.: Энергоатомиздат, 1986. 104 с.
6. Ван дер Варден Б. Л. Алгебра. М.: Наука, 1979. 623 с.
7. Резчиков А. Ф., Кушников В. А., Твердохлебов В. А., Марков А. И. Информационно-измерительный комплекс для диагностирования дефектов геометрических параметров фюзеляжей вертолетов // Авиакосмическое приборостроение. 2012. № 4. С. 35—40.
8. Беллман Р., Калаба Р. Динамическое программирование и современная теория управления / Под ред. Б. С. Разумихина. М.: Наука, 1969. 118 с.
Сведения об авторах
Валентин Иванович Сеньченков — д-р техн. наук, профессор; Военно-космическая академия им. А. Ф. Можайского, кафедра специальных технических систем космических комплексов, Санкт-Петербург; E-mail: [email protected] — адъюнкт; Военно-космическая академия им. А. Ф. Можайского, кафедра специальных технических систем космических комплексов, Санкт-Петербург; E-mail: ponomarev igor [email protected]
Поступила в редакцию 29.03.14 г. космических комплексов
Игорь Николаевич Некрасов
Рекомендована кафедрой специальных технических систем
