ИНФОРМАТИКА
УДК 0.48+628.16.087+631.171:636.5
ОПТИМИЗАЦИЯ ФУНКЦИОНИРОВАНИЯ НЕЧЕТКИХ КОГНИТИВНЫХ КАРТ С ИСПОЛЬЗОВАНИЕМ НЕЙРОННЫХ СЕТЕЙ (НА ПРИМЕРЕ УПРАВЛЕНИЯ ПРОЦЕССАМИ ВОДООЧИСТКИ)
В. Н. ШТЕПА
Учреждение образования «Полесский государственный университет», г. Пинск, Республика Беларусь
Введение
Анализ современных комбинированных систем водоочистки (КСВ), включающих несколько способов воздействия на водные растворы, показал сложность построения эффективных систем управления, поскольку такие объекты характеризуются размытостью и неопределенностью входной информации, что может иметь место, например, при залповых сбросах загрязнителей и/или действии чрезвычайных ситуаций техногенного и природного характера [4], [5].
Существующие системы автоматизации технологических процессов водоочистки не обеспечивают быстрого комплексного реагирования на оперативные изменения ситуационного поведения объектов управления, которые зависят от многих факторов технологического и организационного характера. Улучшить ситуацию можно за счет использования сценариев управления на основе когнитивно-сценарных моделей, а также алгоритмов с применением интеллектуальных механизмов [1], [2]. Поэтому разработка систем многоцелевого управления технологическими процессами водоочистки на основе сценарного подхода будет способствовать повышению производительности, уменьшению удельных потерь и расходов ресурсов и сырья, улучшению экологической безопасности.
При этом нейросетевой подход к оптимизации функционирования когнитивной карты позволяет использовать типовые нейросетевые алгоритмы и концепции обучения, что повышает эффективность получения соответствующих моделей и создает предпосылки для адаптации параметров в режиме реального времени.
Последовательность задач для реализации предлагаемого подхода может быть определена следующим образом [3]: выявление базовых целей, характеризующих исследуемый процесс или систему; выявление ключевых элементов процесса или системы; выявление взаимовлияния факторов; построение когнитивной модели; построение возможных сценариев развития; определение критериев оценки сценариев; оценка сценариев и выявления лучшего с точки зрения выбранных критериев.
Задача оптимального обобщения экспертных мнений с целью построения адекватных нечетких когнитивных карт (НКК) является актуальной, поскольку большинство фактической информации по объекту управления, в нашем случае КСВ, во многом носит субъективный характер (экспертные оценки), причем мнения экспертов по одному и тому же вопросу могут существенно, иногда принципиально, отличаться.
Целью исследования является обоснование и создание нейросетевого блока оптимизации функционирования нечетких когнитивных карт управления комбинированными системами водоочистки с реализацией интерактивного обобщения экспертных
мнений при формировании матрицы значений концептов с целью обеспечения эффективной работы технологического оборудования в режиме реального времени, включая условия действия нештатных ситуаций.
Основная часть
Для достижения цели сформулированы следующие задачи исследований:
- создание структуры НКК для сценарного управления КСВ;
- разработка методики оптимизации работы НКК КСВ с использованием нейронных сетей;
- синтез нейронной сети Кохонена поддержки работы НКК КСВ;
- создание методики использования вероятностной нейронной сети для выбора в режиме реального времени значений коэффициентов НКК КСВ.
Структура нечеткой когнитивной карты комбинированных систем водоочистки. При разработке структуры НКК системы управления комбинированными установками водоочистки исходим из того, что она представляет собой сложную организационно-техническую систему, которая состоит из кортежа:
< D(t), S(t), Y(t), E(t), t >, (1)
где D - действия системы водоочистки; S - влияние факторов внешней среды; Y -выходные показатели работы установки; E - множество концептов и взаимосвязей между ними; t - время функционирования.
В нашем случае для системы управления КСВ задача НКК - определение сценариев работы оборудования.
Указанная система характеризуется большой неопределенностью элементов, которые входят в ее состав (человеческий, экономический и другие факторы), и для моделирования развития такой системы невозможно получить ее точную математическую модель. Поэтому модель комбинированной установки представим в виде обобщенной НКК.
Соответствующая структура нечеткой когнитивной карты разрабатывалась исходя из экспериментальных исследований [4] и объектно-ориентированного анализа предприятий [5]. Исходя из технологических особенностей комбинированных систем водоочистки [2], [4], сценарное планирование выполняется раз в 5 мин.
Экспертно выбраны следующие элементы матрицы взаимовлияния НКК:
1. Промежуточные концепты:
- Е1 - техническое и технологическое оснащение установок водоочистки;
- Е2 - степень использования оборудования;
- Е3 - управление производством.
2. Входные действия:
- Х1 - цены на энергоносители;
- Х2 - расход воды;
- Х3 - рН входной воды;
- Х4 - концентрация взвесей (входная вода);
- Х5 - концентрация нитратов (входная вода);
- Х6 - концентрация фосфатов (входная вода);
- Х7 - биологическая потребность кислорода (БПК) входной воды;
- Х8 - концентрация поверхностно-активных веществ (ПАВ) (входная вода);
- Х9 - концентрация взвесей (входная вода);
- Х1о - температура входной воды;
- Х11 - рН выходной воды;
- Х12 - концентрация взвесей (выходная вода);
- Х13 - концентрация нитратов (выходная вода);
- Х14 - концентрация фосфатов (выходная вода);
- Х15 - БПК выходной воды;
- Х16 - концентрация ПАВ (выходная вода). 3. Исходные действия:
- У1 - финансовые расходы;
- У2 - энергоэффективность.
Модель процессов в комбинированной установке водоочистки представлена в виде соответствующего орграфа (нечеткой когнитивной карты), которая приведена на рис. 1 и иллюстрирует множественные связи и характер взаимодействия факторов.
Рис. 1. Модель комбинированной установки водоочистки в виде нечеткой когнитивной карты
Концепт энергоэффективности (У2) является интегральным показателем, который объединяет расход электроэнергии и качество водоочистки.
Методика формирования, оптимизации и адаптивного использования блока расчета значений весовых коэффициентов НКК КСВ, полученных экспертной оценкой. Формирование значений весовых коэффициентов на основе экспертных оценок должно решить проблему невозможности оперативного опроса экспертов при изменении параметров процессов водоочистки.
Поэтому эксперты (принято количество экспертов равное трем) дают оценку не только типовым, полученным в результате экомониторинга, диапазонам параметров качества воды (рис. 2), но и их потенциально возможным значениям в случае действия нештатных ситуаций (залповые сбросы загрязнителей, порывы и другие непредвиденные ситуации).
При предложенной последовательности (рис. 2) экспертные оценки сводятся в одну таблицу с экспериментальными данными. С целью увеличения гибкости оценки каждый эксперт дает не одно возможное значение, а три наиболее вероятные, по его мнению, значения весовых коэффициентов.
Потом данные группируются для определения единого значения весовых коэффициентов на основе мнения экспертов - используется нейронная сеть.
Рис. 2. Блок-схема алгоритма формирования и использования модуля определения значений весовых коэффициентов концептов НКК (при использовании экспертной оценки)
На следующем этапе формируется адекватная НС, задача которой - расчет значений параметров НКК КСВ исходя из информации, полученной от датчиков ИИК (информационно-измерительного комплекса).
Синтезируемое программное обеспечение включает реализацию работы обоих модулей, что поддерживает условие адаптивности: при выходе системы за пределы установленной эффективности возможно переобучение системы.
Периодическое (по времени) изменение параметров модулей группировки и выбора значений коэффициентов выполняться не будет, поскольку без использования экспертов такие действия могут со временем (при высокой частоте коррекции) внести искажения в работу системы.
Среди всего перечня концептов экспертно определяются: Е\ - техническое и технологическое оснащение установок водоочистки; Е2 - степень использования оборудования; Е3 - управление производством.
Для группировки экспертных оценок (кластеризации) и определения единых значений используем самоорганизационные карты Кохонена.
Модель Кохонена относится к классу алгоритмов векторного кодирования. Она обеспечивает топологическое отображение, оптимально размещает фиксированное число векторов во входном пространстве более высокой размерности, обеспечивая, таким образом, сжатие данных [3].
Алгоритм самоорганизации состоит из следующих этапов [2]:
1) инициализация синаптических весов в сети (с использованием датчика случайных чисел);
2) конкуренция (competition): для любого входного образа и для всех нейронов сети вычисляется значение дискриминантной функции, это является основой конкуренции; нейрон с максимальным значением дискриминантной функции становится победителем;
3) кооперация: нейрон-победитель определяет пространственное расположение соседних возбужденных нейронов;
4) настройки весовых коэффициентов (адаптация): значение дискриминантной функции возбужденных нейронов увеличивается для данного образа путем настройки весовых коэффициентов. При адаптации отзыв нейрона-победителя на близкий входной образ увеличивается.
Размерность входного пространства (данных): x = [x1, ..., xm]T - входной образ. Вектор синаптических весовых коэффициентов для j-го нейрона:
ш j = [ш И m ]T,
(2)
где l - число нейронов сети. Найдем нейрон-победитель:
i(х) = arg max ш^;
i,i
max wTx о min х - w. j j " '
i(x) = argmin||x -ш j|, j = 1, l. (3)
Выражение (3) описывает процесс конкуренции. Непрерывное пространство входных образов отражается на дискретное пространство нейронов в процессе конкуренции между нейронами сети.
Нейрон-победитель определяет центр группы (окрестность) нейронов, участвующих в кооперации. Нейрон-победитель сильнее влияет на топологически близких соседей, чем на более отдаленные нейроны.
Пусть da - латеральное расстояние между данным нейроном и возбужденным
нейроном. Топологическая окрестность - это унимодальная функция латерального расстояния, которая удовлетворяет двум требованиям:
- hj - симметричность относительно центра;
- hj - амплитуда монотонно уменьшается с увеличением латерального расстояния и при dj ^ 0.
Этим условиям удовлетворяет функция Гаусса:
h (х)=exp
' d2 ^
2g2
d2 =1Г - j2, (4)
где г- позиция возбужденного нейрона; у, т{ - дискретная позиция нейрона-победителя.
Обе координаты определяются в дискретном выходном пространстве. Ширина топологической окрестности со временем уменьшается. Если п - дискретное время, то экспоненциальное падение обеспечивает зависимость:
ст(п) = о0е"т1, п = 0,1,2,..., (5)
где т1 - временная константа, а0 - начальное значение. Функция окрестности:
. (-)(п) = ехР
г Я ^
2С2(П) ;
п = 0,1,2,.... (6)
Впоследствии окрестность сужается, так как его ширина уменьшается. Нейроны окрестности будут участвовать в адаптации весовых коэффициентов. Цель кооперации -скоррелировать изменение весовых коэффициентов нейронов одной окрестности.
Адаптация заключается в изменении весового коэффициента wj в зависимости
от входного вектора х. Она базируется на постулате обучения Хебба: правильные связи усиливаются, а ложные ослабевают.
Однако в случае самоорганизации это правило неприменимо, поскольку не известен целевой выход. Если связи будут модифицироваться только в сторону усиления, то вскоре все они достигнут насыщения.
Модификация весовых коэффициентов вычисляется по формуле
=^}х - %(у. , (7)
где "л - коэффициент скорости обучения.
Выберем линейную функцию:
У.) = Л У..
Тогда (7) будет иметь вид:
= ЛУ.(х - ) = л. х)(х - ). (8)
Тогда при переходе от момента времени п к п +1 получим:
о.(п +1) = ш.(п) + л(п)к] ях)(п)(х - о.(п)), . =1,1. (9)
Таким образом, модифицируются весовые коэффициенты всех нейронов в окрестности нейрона-победителя I. Значение весового вектора wi нейрона-победителя I приближается к х. Векторы синаптических весовых коэффициентов отслеживают распределение входных векторов в соответствии с выбором окрестности, обеспечивая тем самым топологическое составление карты признаков во входном пространстве.
Синтез нейросетевого блока расчета значений весовых коэффициентов НКК КСВ, полученных экспертной оценкой. Используя пакет прикладных математических программ «^айБйса» на основании экспериментальных исследований [4] синтезируем нейронную сеть Кохонена (рис. 3). Для каждого набора экспериментальных данных создается своя НС.
Эксперт 1 Эксперт 2 Эксперт 3
Рис. 3. Архитектура самоорганизационной нейронной сети Кохонена определения числовых значений параметров НКК (экспертные оценки)
Адекватность НС подтверждена оценкой относительной среднеквадратической ошибки (рис. 4): для учебной выборки - 0,83 %, для контрольной выборки - 0,03 %.
Относительная ошибка,
2,6 ,—
2,4 ■
2,2 ■
2,0 ■
1,8 ■
1,6 ■
1,4 ■
1,2 -
1,0 -
0,8 -
0,6 ■
0,4 ■
0,2 ■
0,0 ■
-0,2 ■
-0,4 --200
Рис. 4. Проверка адекватности самоорганизационной НС Кохонена определения числовых значений параметров НКК (экспертные оценки)
С использованием синтезированной топологической карты НС Кохонена (рис. 5) и предложенной методики рассчитываются единичные значения для всех весовых коэффициентов концептов НКК КСВ:
- значения коэффициентов влияния концепта Е\ наУ —0,77 и У2 - 0,79;
- значения коэффициентов влияния концепта Е2 на У1 —0,78 и У2 - 0,80;
- значения коэффициентов влияния концепта Е3 на У1 —0,89 и У2 - 0,90.
После расчета с помощью самоорганизационной НС Кохонена единичных экс-пертно определенных значений весовых коэффициентов НКК КСВ при разных комбинациях входных параметров необходимо решить задачу их выбора в режиме реального времени в зависимости от значений с датчиков информационно-измерительного комплекса (ИКК). Для этого предлагается использовать вероятностные нейронные сети (РКЫ), которые являются частным случаем Байесовских сетей. Такая сеть не требует обучения в том смысле, которое необходимо для сетей типа персептрона, радиаль-но-базисной функции и других архитектур, так как все параметры РКЫ-сети (число элементов и значения весов) определяются учебными данными [2].
Рис. 5. Плато прогноза самоорганизационной НС Кохонена определения числовых значений параметров НКК (экспертные оценки)
В основу классификации в сети РМК положена идея, что для каждого образца можно принять решение на основе выбора наиболее вероятного класса из тех, которым мог бы принадлежать этот образец. Такое решение требует оценки функции плотности вероятности для каждого класса. Эта оценка устанавливается в результате рассмотрения учебных данных. Правилом является то, что класс с плотным распределением в области неизвестного образца будет более многочисленным, по сравнению с другими классами. Точно так же будет иметь преимущество и класс с высокой априорной вероятностью или высокой ценой ошибки классификации.
Для двух классов А и В, согласно этому правилу, выбирается класс А, если
ИлсЛ!Л( х) = ИБсБ/Б( х) (10)
где И - априорная вероятность; с - цена ошибки классификации; Дх) - функция плотности вероятностей.
В контексте нашей задачи интересует не столько дискретная классификация, сколько значение выхода слоя добавления, который высчитывает плотности распределения вероятности для совокупности экземпляров соответствующего образа. То есть на выходе этого слоя сможем отслеживать динамику изменения (потенциальной смены) коэффициентов.
Для этого традиционно используется форма функции Гаусса, куда включен квадрат евклидового расстояния от неизвестного образа к элементу слоя образцовых образов:
( II || 2 Л
" 1-1 |х - X 8(х) = ^ ехр
1=1
а2
(11)
где а - параметр, который задает ширину функций Гаусса; х - неизвестный входной образец.
Выход нейронной сети - номер класса (значения коэффициентов), к которому относится полученный набор входных величин. При эксплуатации предложенного блока оптимизации работы НКК КСВ информация подается на отдельные РМЫ-сети, а они реализуют выбор экспертных значений (рис. 6) в зависимости от данных с ИКК.
Рис. 6. РМЫ-сеть определения значений весовых коэффициентов НКК КСВ (структура использования)
Заключение
Использование нейронных сетей Кохонена (оптимизация экспертных заключений) и вероятностных нейронных сетей (выбор значений коэффициентов) позволяет повысить функциональную эффективность НКК КСВ путем обеспечения возможности адаптации в режиме реального времени ее параметров в зависимости от влияния техногенных и природных факторов: качества воды, ее расхода, состояния технологического оборудования.
Литература
1. Лисенко, В. П. Спещальш роздши вищо'1 математики (Неч^ю множини) / В. П. Лисенко, Б. В. Кузьменко. - К. : НАУ, 2004. - 83 с.
2. Люггер, Дж. Ф. Искусственный интеллект: стратегии и методы решения сложных проблем : пер. с англ. / Дж. Ф. Люггер. - М. : Вильямс, 2005. - 864 с.
3. Рассел, С. Искусственный интеллект: современный подход : пер. с англ. / С. Рассел, П. Норвиг. - 2-е изд. - М. : Вильямс, 2006. - 1408 с.
4. Штепа, В. М. Обгрунтування алгоритму експериментально-анал^ичних досль джень режимiв електротехшчно! очистки спчних вод агропромислових об'ект1в з метою побудови енергоефективних систем управлшня / В. М. Штепа // Енергети-ка i автоматика. - 2012. - 01 (11). - Режим доступа: http://nbuv.gov.ua/ j-pdf7eia_2014_2_10.pdf. - Дата доступа: 10.09.2015.
5. Штепа, В. М. Обгрунтування арх^ектури системи управлшня комплексними методами очистки спчних вод промислових об'екпв / В. М. Штепа // Вюник Харюв-ського нацюнального техшчного ушверситету сшьського господарства iменi Петра Василенка. - Харюв : ХНУТСГ, 2014. - Вип. 154. - С. 48-50.
Получено 02.12.2015 г.