ПРОГРАММНЫЕ СИСТЕМЫ: ТЕОРИЯ И ПРИЛОЖЕНИЯ
ISSN 2079-3316 № 5(14), 2012, c. 71–79
УДК 519.8
Н. Ю. Жбанова
Моделирование процесса варки сахара разностной
нейронечёткой системой с переключениями
Аннотация. В статье рассматривается применение разностных нейронечётких переключаемых систем к моделированию процесса варки сахара.
Особое внимание уделяется учёту не только текущих, но и прошлых значений на входах модели. Предложен новый тип входных нечётких множеств, который позволяет снизить количество настраиваемых параметров нейронечёткой системы.
Ключевые слова и фразы: переключаемые нейронечёткие системы, разностные нечёткие модели, варка сахара.
Введение Важную роль в производстве сахара играет варка сахарного сиропа. Она осуществляется в вакуум-аппаратах – устройствах, поддерживающих условия, необходимые для образования и роста кристаллов сахара. Конструкция вакуум-аппарата позволяет воздействовать на три параметра варки посредством управляемых вентилей. Это (уровень сиропа), (вакуум) и (давление пара). Эти три параметра определяют четвёртый – (плотность сиропа). Именно оказывает наиболее существенное влияние на качество готового сахара.
Чтобы в итоге получился сахар нужного качества, параметр (а, значит, и параметры,, ) должен проследовать по некоторому графику. Графики параметров представляют собой ломаные линии, вид которых определяется несколькими контрольными точками.
Оператор выставляет эти точки, основываясь на опыте, а программа согласно им регулирует параметры варки.
Работа поддержана РФФИ, проект № 11-07-00580-а.
c Н. Ю. Жбанова, c Липецкий государственный технический университет, c Программные системы: теория и приложения, 72 Н. Ю. Жбанова Была поставлена задача разработки системы, которая в реальном времени рекомендовала бы оператору значения контрольных точек с учётом требуемого качества сахара, опираясь на снятые датчиками значения параметров,,.
Повторить рассуждения оператора на основе информации о значениях параметров варки можно с помощью нечёткой модели [1]. Модель с тремя входами, получая снятые датчиками значения,,, на выходе формировала бы контрольные точки.
Заметим, что в распоряжении автора есть данные по варке за сезон 2009 г., из которых для нечёткой модели можно сформировать обучающее множество. Это позволит избежать «ручной» настройки её параметров. Модели, по структуре аналогичные нечётким, и обладающие возможностью самонастройки, называются гибридными нейронечёткими системами [1].
Однако несколько входов, связанных на разных этапах варки разными зависимостями, сложно учесть в одной нейронечёткой модели.
Поэтому предлагается разбить варку на этапы, и моделировать каждый этап отдельной подмоделью с одинаковой структурой (три входа, один выход). Так от классической нейронечёткой модели мы перейдём к нейронечёткой модели с переключением.
1. Нечёткие и нейронечёткие модели с переключениями Нечёткие модели с переключениями – синтез переключаемых систем и нечётких моделей. Нечёткая переключаемая модель состоит из нескольких классических нечётких моделей (подмоделей) и переключающего сигнала [2, 3].
Переключаемую нечёткую Такаги –Сугено (TS) модель типа MISO (multiple input, single output) можно описать совокупностью правил вида:
: (1 = )... ( = ) = a u +.
(1) 1 Здесь = 1,..., – количество правил в каждой подмодели, – –– входные нечёткие множества, = 1,..., – количество входов, a =,...,, – параметры функций, содержащихся в – 1 заключениях правил, u = (1,..., ) – входные значения, = – (1, 2,..., ) –– переключающий сигнал. Каждая подмодель представляет собой классическую нечёткую TS-модель.
Выход TS-модели с переключением задается следующей формулой:
Моделирование варки сахара ) a u + / ( (2) =.
=1 = Уровень истинности предпосылки правила определяется по. Здесь – степень принадлежности входа формуле = = нечёткому множеству. Если используются нечёткие множества гауссовского типа, то степень принадлежности вычисляется так:
Здесь и — параметры (центр и ширина) гауссовского нечёткого множества. Нейронечёткие системы объединяют достоинства нейронных сетей и нечётких моделей. Существуют разные виды синтеза нечётких моделей и нейронных сетей, но самую широкую область применения нашли упомянутые выше гибридные нейронечёткие системы (ННС).
ННС –– это нечёткая модель, которая может обучаться, как нейронная сеть. Самый распространенный тип структуры ННС – ANFIS [1].
Функционально такая система эквивалентна нечёткой TS-модели.
По аналогии с переключаемыми нечёткими моделями, определим переключаемую нейронечёткую модель как систему, подсистемами которой являются собственно нейронечёткие модели.
В общем случае нейронечёткая переключаемая ANFIS-модель, как и нечёткая переключаемая модель Такаги–Сугено, имеет базу правил вида (1). Подмодели представляют собой ННС типа ANFIS с входами.
Реализуемая такой нейронечёткой переключаемой ANFIS-моделью функция имеет вид (2). Параметры заключений правил a, и нечетких множеств подбираются не вручную, а в процессе обучения.
В статьях нейронечёткие ANFIS-модели с переключениями не упоминаются (может быть потому, что они эквивалентны нечётким переключающимся TS-моделям, и отличаются от них лишь возможностью автоматической настройки параметров).
2. Разностная нейронечёткая модель с переключениями Большими возможностями по сравнению с моделями вида (1) обладают нечёткие модели Такаги–Сугено типа MISO, правила которых имеют следующий вид:
(4) Здесь a, b – векторы параметров заключений правил, – набор входных нечётких множеств.
Аргумент предпосылки правил такой базы X() задается в виде матрицы X() = [(), u(),..., ( + 1), u( + 1)] (+1), где u( + 1) 1 – вектор входных значений, ( + 1) 1 – выходное значение модели. Таким образом, в аргументе X() заложена информация не только о текущих значениях входов и выхода, но и о значениях в предшествующие моменты времени. В статьях такие модели называют разностными моделями, моделями с экзогенным входом, с глубиной памяти, или с запаздыванием [4].
Введя в разностную нечёткую TS-модель возможность настройки параметров заключений правил a, b, и параметров функций принадлежности, получим разностную ANFIS-систему.
Модели (4) учитывают большее количество информации о моделируемом объекте, и более точно отражают реальную ситуацию.
Далее будет рассматриваться частный случай модели (4) – модель – базой правил (5):
(5) –– векторы параметров заключений правил. Аргументом предпосылок правил является матрица U() = [u(),..., u( + 1)] размера, включающая в себя только входные значения (). Выход модели (5) определяется стандартно.
В развернутом виде база правил (5) имеет вид:
Рис. 1. Стандартная гауссовская функция принадлежности (а) и предлагаемая функция принадлежности (б) Степень принадлежности входа ( + 1) гауссовскому нечёткому множеству вычисляется по формуле (3). Множества гауссовского типа представлены на рис. 1 (а). Чтобы идентифицировать нечёткую модель (6), требуется определить параметры набора гауссовских нечётких множеств для каждого элемента входной матрицы ( + 1), = 1,...,, = 1,...,. Большое количество параметров, требующих идентификации, является недостатком моделей с запаздыванием.
Предлагается использовать функции принадлежности, которые фаззифицируют не одно входное значение ( + 1), а сразу вектор из элементов u () = [ (),..., ( + 1)]. При этом предполагается, что входной параметр изменяется во времени по линейному или близкому к линейному закону –– в этом случае вектор u () можно фаззифицировать предлагаемой трёхмерной функцией принадлежности, центр которой определяет прямая. Трёхмерная функция принадлежности определяет близость входного вектора к прямой так же, как гауссовская функция принадлежности определяет близость входного значения к своему центру.
На рис. 1 (б) показана предлагаемая трёхмерная функция принадлежности, оценивающая близость входного вектора u() к центрам нечётких множеств 1, 2. Степень принадлежности вектора u () нечёткому множеству вычисляется по формуле:
У предложенного нечёткого множества три параметра: 0 и 1 –– параметры прямой, определяющей центр множества, – пара- – метр, определяющий ширину нечёткого множества.
В случае использования множеств со степенью принадлежности, определяемой по формуле (7), упростится идентификация нечёткой модели –– для нее понадобится задать уже не ·, а наборов нечётких множеств. С учетом того, что у предложенных функций принадлежности на один параметр больше, чем у гауссовских, число настраиваемых параметров снизится в 2/3 раз.
Правила разностной нечёткой TS-модели в случае использования предложенных функций принадлежности имеют вид:
Здесь U() – матрица входных значений, () 1 – выход системы; a 1, b 1 – векторы коэффициентов системы.
Видно, что по сравнению с базой правил (6) функций принадлежности требуется меньше.
Предложенные функции принадлежности можно использовать и в сочетании с переключающимися нечёткими или нейронечеткими системами. База правил для переключающейся системы имеет вид:
(8) Ниже рассматривается применение модели (8) для описания процесса варки сахара.
3. Моделирование процесса варки разностной переключаемой ННС Заметим, что ход, результат процесса варки и контрольные точки зависят не только от конкретного значения параметров, и в каждый момент времени, но и от скорости их изменения. Это факт обосновывает применение для описания варки разностной нейронечеткой модели, учитывающей запаздывание входных параметров.
Как было сказано во введении, сложность процесса варки привела к необходимости разбить его на этапы и вместо классической нейронечёткой модели использовать модель с переключением. Всего было выделено 4 этапа варки, и каждому поставлена в соответствие отдельная подмодель.
Все подмодели представляют собой разностные ННС типа ANFIS.
Набор правил подмодели имеет вид:
Здесь переключающий сигнал = 1,..., 4 определяется значением уровня сиропа. Набор множеств для фаззификации каждого из трех входов состоял из двух нечётких множеств предложенного в п. 2 типа,. Количество правил в каждой подмодели = 1,..., 8.
Входная матрица U() 315 каждой подмодели составляется из значений трех параметров варки,, с одинаковой глубиной памяти = 15. Степень принадлежности входных векторов u1 () = [(),..., ((14)], u2 () = [(),..., ((14)], u3 () = [(),..., ((14)] нечётким множествам вычисляется по формуле (7). Угловые коэффициенты прямых, являющихся центрами нечётких множеств, вносят в подмодели информацию о скорости изменения параметров варки.
На выходе каждая подмодель формирует значение контрольной точки для своего этапа.
Обучающее множество представляет собой данные по варке сахара за сезон 2009 года, предоставленные Добринским сахарным заводом. Данные по варке – это значения 3 параметров варки, измеренные с интервалом 10 секунд, сведения о качестве полученного сахара, и значения контрольных точек. Для обучения модели были отобраны те циклы варки, в которых получился сахар хорошего качества.
Подмодели одинаковы по структуре и отличаются значениями параметров заключений правил a, b, и параметров функций принадлежности, которые подбираются в процессе обучения; при этом использование предложенных функций принадлежности (7) позволяет сократить число настраиваемых параметров ННС в 10 раз (по сравнению с гауссовскими функциями принадлежности).
4. Результаты моделирования Разностная переключаемая ННС была запрограммирована в MATLAB. В обучающее множество вошли 15 циклов варки, в тестовое – – 20 циклов. В таблице 1 приведены реальные и модельные значения контрольных точек для 4 циклов варки. Видно, что модель вычисляет значение контрольной точки, близкое к реальному. Это значит, что она правильно описывает варку. Отклонения модельных значений контрольных точек от реальных по всем циклам тестового множества не превышали 2 % (при абсолютных значениях контрольных точек примерно 80 brix).
Для подтверждения адекватности модели был вычислен коэффициент детерминации: 2 = 0, 917. Высокое значение говорит о достоверности построенной модели.
Список литературы [1] Пегат А. Нечёткое моделирование и управление. М. : Бином, 2009. – 800 c.
[2] Ojleska V., Stojanovski G. Switched Fuzzy Systems: Overview and Perspectives // 9th International PhD Workshop on Systems and Control: Young Generation Viewpoint. – Izola, 2008. [3] Котов К. Ю., Шпилевая О. Я. Переключаемые системы: устойчивость и проектирование (обзор) // Автомерия. Системы автоматизации в научных исследованиях и промышленности – Новосибирск, 2008. Т. 44, c. 71–87. [4] Кудинов Ю. И., Венков А. Г., Келина А. Ю. Моделирование технологических и экологических процессов: монография. Липецк : ЛЭГИ, 2001 ISBN 5– 90037–22–3. – 130 c. Рекомендовал к публикации Программный комитет Молодёжной школы-семинара Ассистент кафедры прикладной математики ЛГТУ.
Образец ссылки на эту публикацию:
Н. Ю. Жбанова. Моделирование процесса варки сахара разностной нейронечёткой системой с переключениями // Программные системы:
теория и приложения : электрон. научн. журн. 2012. T. 3, № 5(14), с. 71–79.
URL: http://psta.psiras.ru/read/psta2012_5_71-79.pdf N.Yu. Zhbanova. Switched Difference Neuro-Fuzzy Modelling of a Sugar-making Process.
Abstract. This article surveys switched difference neuro-fuzzy system application to sugar-making process modelling. Particular attention is given to system’s input vector, which consists of not only current but also past values. The new type of membership function, which reduces the amount of neuro-fuzzy system’s parameters, is suggested.
Key Words and Phrases: Switching Neuro-Fuzzy Systems, Difference Fuzzy Models, Sugar-making Process.