Прямой и непрерывный алгоритмы муравьиной колонии для идентификации параметров нечетких моделей

Изучение эффективности непрерывного и прямого алгоритма муравьиной колонии. Определение особенностей распределения ошибки в зависимости от размера архива, коэффициента испарения феромона. Решение задачи параметрической идентификации нечетких моделей.

Рубрика Математика
Вид статья
Язык русский
Дата добавления 19.01.2018
Размер файла 47,0 K

Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже

Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.

Размещено на http://www.allbest.ru/

ПРЯМОЙ И НЕПРЕРЫВНЫЙ АЛГОРИТМЫ МУРАВЬИНОЙ КОЛОНИИ ДЛЯ ИДЕНТИФИКАЦИИ ПАРАМЕТРОВ НЕЧЕТКИХ МОДЕЛЕЙ

Ходашинский И.А., д.т.н.,

Дудин П.А., профессор

Томский государственный университет систем управления и радиоэлектроники

e-mail: hodashn@rambler.ru

1. Постановка задачи

Рассматривается нечеткая модель типа синглтон, i-ое правило в которой имеет следующий вид:

IF x1=A1i AND x2=A2i AND … AND xn=Ani THEN y = ri,

где Aij - лингвистический терм, которым оценивается переменная xi;

ri - действительное число, которым оценивается выход y.

Нечеткая модель осуществляет отображение :

где x - входной вектор, R - число правил; n - количество входных переменных;- функция принадлежности (ФП), определяемая набором своих параметров, например, треугольная - тремя параметрами, трапециевидная - четырьмя, гауссова и параболическая - двумя.

Нечеткая модель может быть представлена как

y = f(x, и),

где и = ||и1,…, иN|| - вектор параметров, N = n (число параметров, описывающих одну функцию принадлежности) ·x - (число термов, описывающих одну входную переменную), y - скалярный выход модели.

Пусть дано множество обучающих данных {(xp; tp), p = 1,..., m}, тогда среднеквадратическая функция ошибки, которую необходимо минимизировать, будет иметь вид

Для решения проблемы минимизации предлагается использовать алгоритмы муравьиной колонии.

2. Непрерывный алгоритм муравьиной колонии

Классический алгоритм муравьиной колонии - процедура дискретной оптимизации [1], здесь при выборе очередной дуги, муравей руководствуется дискретным распределением вероятности.

В случае непрерывного алгоритма (НАМК) выбор, который делает муравей, не ограничен конечным множеством. Для этого дискретное распределение заменяется непрерывным, т.е. функцией плотности вероятности (ФПВ) [2].

В НАМК для описания используется ФПВ с гауссовым ядром. Под Гауссовым ядром Gi(x) понимается функция, основанная на взвешенной сумме нескольких одномерных гауссовых функций gli(x):

Каждому параметру ФП i соответствует свое гауссово ядро, i = 1, ..., N, N - число настраиваемых параметров нечеткой модели. Каждая функция Gi(x) описывается тремя векторами: - вектор весов, связанных с индивидуальными гауссовыми функциями, i - вектор математических ожиданий, и i - вектор среднеквадратичных отклонений. Количество элементов всех этих векторов равно числу функций Гаусса, составляющих гауссово ядро.

В НАМК вводится понятие архива решений. Архив решений представлен таблицей, в которой k строк. Каждая строка представляет собой найденное муравьем решение sl=||sl1, sl2,…, slN||, ошибку E и вес решения l. Решения упорядочены в архиве согласно их качеству.

Ниже приведен алгоритм работы НАМК для оптимизации параметров ФП:

Шаг 1. Задать начальные параметры.

Шаг 2. Сгенерировать популяцию муравьев в колониях.

Шаг 3. Сгенерировать несколько случайных решений.

Шаг 4. Найти значения вектора весов. Сделать текущим первого муравья первой колонии.

Шаг 5. Для текущего муравья текущей колонии вычислить номер l, используемой функции Гаусса. Определить параметры функций Гаусса. Сгенерировать вектор случайных величин ||и1,…, иN|| на основе полученных функций распределения.

Шаг 6. Вычислить ошибку E при найденных параметрах ||и1,…, иN||, если ошибка меньше текущей, то сохранить новые параметры.

Шаг 7. Добавить в архив новое решение, ранжировать архив, удалить из архива худшее решение.

Шаг 8. Если в текущей колонии имеется следующий муравей, то сделать его текущим и перейти к шагу 5.

Шаг 9. Если имеется следующая колония, то сделать текущим первого муравья в этой колонии и перейти на шаг 5.

Шаг 10. Если условие окончание работы алгоритма выполнено, то закончить, иначе сделать текущим первого муравья первой колонии и перейти к шагу 5.

Условием окончания работы алгоритма является достижение заданного числа итераций, либо получение ошибки, меньше заданной.

3. Прямой алгоритм муравьиной колонии

В прямом алгоритме [3] муравей отвечает за вычисление значений закрепленного за ним параметра, поэтому муравьев в алгоритме столько, сколько параметров нечеткой модели. Каждый i-й муравей создает свое решение, генерируя нормально распределенное действительное число N(мi, уi). В алгоритме используются два вида ферамонов: первый связан с центрами нормальных распределений м = ||м1, …, мN||, второй с разбросом у = ||у1 …, уN||. Количество ферамона определяет значения параметров м и у. Для каждого параметра иj задан интервал изменения [aj, bj], где bj и aj - верхняя и нижняя граница параметра иj.

В качестве начальных значений для параметров м используются заданные случайным образом или другим способом значения параметров ФП. Начальные значения параметров у вычисляется по следующей формуле:

.

После того как муравьи нашли решения, определяется испарение ферамона. Для текущей t-й итерации испарение определяется следующим образом:

м(t) = (1- с) м(t-1)

у(t) = (1- с) у(t-1),

где с - эмпирический коэффициент испарения феромона, заданный на интервале [0, 1].

Далее происходит нанесение феромона:

м(t) = м(t) + си(t)

у(t) = у(t) + с|и(t) - м(t)|,

где и(t) - решение, найденное муравьиной колонией на текущей итерации, оно совпадает с глобальным лучшим решением.

Особенностью прямого алгоритма муравьиной колонии является включение в него простейшего локального поиска, состоящего из двух этапов: на первом значение параметра иj увеличивается с определенным шагом до значения иj + dj, на втором этапе значение параметра уменьшается с определенным шагом до значения иj - dj. Значение dj определяется по формуле:

dj = уj rand,

где rand - случайное равномерно распределенное число в интервале [0,1]. Шаг вычисляется по следующей формуле:

stj = dj / K,

где K - целое число, отвечающее за вычисление значения шага.

В результате локального поиска определяется новый вектор параметров и. Значения этих параметров передаются в нечеткую систему в качестве новых значений параметров ФП. Вычисляется ошибка и лучшее решение текущего шага. Глобальное лучшее решение запоминаются.

Для решения задачи параметрической идентификации необходима проверка изменения параметра иj на ограничения, накладываемые нечеткой системой. К таким ограничениям относится покрытие всей области определения термов, описывающих входную переменную, и упорядоченность пиков ФП.

Для преодоления локальных минимумов в алгоритме используется обновление параметров у. С этой целью введен параметр конвергенции, вычисляемый по следующей формуле:

Когда алгоритм приближается к локальному минимуму, коэффициент конвергенции cf приближается к 0. Как только коэффициент конвергенции становится меньше критического значения cfr , то вектор у возвращается в начальное состояние.

Ниже приведен пошаговый алгоритм:

Шаг 1.Сгенерировать N муравьев, инициализировать для них начальные параметры. Сделать текущим первого муравья.

Шаг 2. Для текущего муравья вычислить значение параметра dj. Вычислить шаг stj.

Шаг 3. Для текущего муравья, увеличить значение параметра иj на величину шага. Проверить полученное значение на ограничения, накладываемые нечеткой системой. Проверить возможное будущее значение параметра, полученное в результате применения формул испарения и нанесения феромона, на ограничения. Если предыдущие проверки допускают данное изменение, то вычислить ошибку.

Шаг 4. Если полученное значение не больше иj + dj, то перейти на шаг 3.

Шаг 5. Для текущего муравья, уменьшить значение параметра иj на величину шага. Проверить полученное значение на ограничения, накладываемые нечеткой системой. Проверить возможное будущее значение параметра, полученное в результате применения формул испарения и нанесения феромона, на ограничения. Если предыдущие проверки допускают данное изменение, то вычислить ошибку.

Шаг 6. Если полученное значение не меньше иj - dj, то перейти на шаг 5.

Шаг 7. В качестве нового значения параметра иj выбирать значение, дающее наименьшую ошибку.

Шаг 8. Если есть следующий муравей, то сделать его текущим и перейти на шаг 2.

Шаг 9. Для всей колонии муравьев выполнить операции испарения и нанесения феромона.

Шаг 10. Если cf<cfr, вернуть вектор параметров у в начальное состояние.

Шаг 11. Если условие окончания работы алгоритма выполнено, то КОНЕЦ, иначе сделать текущим первого муравья и перейти на шаг 2.

Условием окончания работы алгоритма может быть выполнение заданного числа итераций Iter или достижение ошибки меньше заданной.

4. Эксперимент

Суть эксперимента заключалась в аппроксимации нечеткой системой следующей тестовой функции: f(x1, x2) = x1*sin(x2). Входные переменные описаны пятью нечеткими термами. База содержит 25 нечетких правил. Термы заданы треугольными ФП. Критерием аппроксимации является ошибка вывода E. Некоторые выводы по каждому из алгоритмов приведены ниже.

К сожалению, НАМК имеет свойство застревать в локальных минимумах, по этой причине необходимо говорить о лучших и худших решениях при одних и тех же параметрах алгоритма. Разные решения получаются в силу случайности, присутствующей в самом алгоритме, однако не всегда этот случайный выбор способен вывести алгоритм из локального минимума. Эксперимент показал, что худшее решение (большие ошибки) не зависит от числа итераций, лучшие решения

Рис. 1. Распределение ошибки в зависимости от размера архива

(малые ошибки) появляются уже при числе итераций, равном 75. Похожая зависимость наблюдается и при изменении числа муравьев. Время работы НАМК линейно зависит от числа итераций и от числа муравьев. На рисунке 1 приведено распределение ошибки при изменении размера архива решений, здесь каждая точка соответствует решению, найденному НАМК. Худшее решение не зависит от размера архива, наилучшие решения получены при размере архива 12-17 . Размер архива не оказывает существенного влияния на время вывода.

В прямом алгоритме муравьиной колонии исследовались три параметра алгоритма: коэффициент испарения феромона с, параметр K и число итераций алгоритма Iter. Средняя оценка ошибки определялась по результатам 25 прогонов при заданном сочетании параметров алгоритма.

Параметр K тестировался в диапазоне от 4 до 20. При изменении параметра от 4 до 15 ошибка уменьшается от 4,36892E-4 до 2,5423E-7, дальнейшее увеличение K не приводит к существенному уменьшению ошибки.

Параметр Iter тестировался в диапазоне от 3 до 15. Изменение параметра от 3 до 12 позволяет уменьшить ошибку с 2,42398E-4 до 1,58104E-7, дальнейшее увеличение числа итераций не приводит к уменьшению ошибки.

Параметр с протестирован в диапазоне от 0,1 до 0,9999 на трех различных базах правил. Первая база содержала 25 нечетких правил, обе входные переменные описаны пятью нечеткими термами. Вторая база содержала 36 нечетких правил, входные переменные описаны шестью нечеткими термами. Третья база содержала 14 правил, входная переменная x1 задана двумя нечеткими термами, переменная x2 описана семью нечеткими термами. На рис.2 видно, что средняя оценка ошибки уменьшается, когда с увеличивается от 0,1 до 0.96. При дальнейшем увеличении с ошибка резко увеличивается.

Рис. 2. Распределение ошибки в зависимости от коэффициента испарения феромона

Заключение

В работе рассмотрены два алгоритма муравьиной колонии для решения задачи параметрической идентификации нечетких моделей. Результаты эксперимента показали, что эффективность непрерывного алгоритма несколько выше по сравнению с прямым алгоритмом. Эффективность и скорость сходимости непрерывного и прямого алгоритмов много выше классического алгоритма муравьиной колонии, в котором переход от дискретного типа оптимизации к непрерывному выполнен путем построения полного ориентированного графа, количество вершин в котором определяется точностью нахождения значений параметров нечеткой модели [4].

непрерывный муравьиный колония параметрический

Литература

1. Dorigo M., Maniezzo V., Colorni A. Ant System: Optimization by Colony of Cooperating Agents // IEEE Transaction Systems, Man and Cybernetics. Part B. - 1996. - Vol. 26. - P. 29-41.

2. Socha K., Dorigo M. Ant Colony Optimization. Faculte Libre de Bruxelles. -- Bruxelles, 2005. - 34p. or Continuous Domains/ Technical Report TR/IRIDIA/2005-037, Universit

3. Kong M., Tian P. Application of ACO in Continuous Domain / Ed. By L. Jiao et al.. ICNC-2006, LNCS 4222, Part II. - Berlin, Springer-Verlag, 2006. - P. 126-135.

4. Ходашинский И.А., Дудин П.А. Параметрическая идентификация нечетких моделей на основе гибридного алгоритма муравьиной колонии // Автометрия. -2008. -Том 44, № 5. - С. 24-35.

Размещено на Allbest.ru


Подобные документы

  • Механизмы реализации эвристических алгоритмов муравьиной колонии. Основная идея - использование механизма положительной обратной связи, помогающего найти наилучшее приближенное решение в сложных задачах оптимизации. Области применения алгоритма муравья.

    реферат [361,6 K], добавлен 07.05.2009

  • Математическая теория нечетких множеств, история развития. Функции принадлежности нечетких бинарных отношений. Формирование и оценка перспективного роста предприятия оптовой торговли. Порог разделения ассортимента, главные особенности его определения.

    контрольная работа [22,3 K], добавлен 08.11.2011

  • Функция принадлежности в форме трапеции, ее представление. Составление проекта бюджета. Сумма и разность нечетких переменных. Операция нечеткого выбора. Порядок вычисления бюджета. Решение задачи с использованием трапециевидной функции принадлежности.

    презентация [32,5 K], добавлен 15.10.2013

  • Математическая теория нечетких множеств и нечеткая логика как обобщения классической теории множеств и классической формальной логики. Сферы и особенности применения нечетких экспертных систем. Анализ математического аппарата, способы задания функций.

    презентация [1,0 M], добавлен 17.04.2013

  • Оптимальная настройка параметров "алгоритма отжига" при решении задачи коммивояжера. Влияние начальной температуры, числа поворотов при одной температуре и коэффициента N на результат. Сравнение и определение лучшей функции для расчётов задачи.

    контрольная работа [329,9 K], добавлен 20.11.2011

  • Значения коэффициента регрессии (b) и сводного члена уравнения регрессии (а). Определение стандартной ошибки предсказания являющейся мерой качества зависимости величин Y и х с помощью уравнения линейной регрессии. Значимость коэффициента регрессии.

    задача [133,0 K], добавлен 21.12.2008

  • Нахождение полинома Жегалкина методом неопределенных коэффициентов. Практическое применение жадного алгоритма. Венгерский метод решения задачи коммивояжера. Применение теории нечетких множеств для решения экономических задач в условиях неопределённости.

    курсовая работа [644,4 K], добавлен 16.05.2010

  • Описание метода потенциалов Математическая постановка задачи об оптимальных перевозках. Метод решения задачи об оптимальных перевозках средствами Ms Excel. Постановка параметрической транспортной задачи, ее математическое и компьютерное моделирование.

    курсовая работа [802,5 K], добавлен 21.10.2014

  • Построение модели множественной регрессии теоретических значений динамики ВВП, определение средней ошибки аппроксимации. Выбор фактора, оказывающего большее влияние. Построение парных моделей регрессии. Определение лучшей модели. Проверка предпосылок МНК.

    курсовая работа [352,9 K], добавлен 26.01.2010

  • Числовые характеристики непрерывных величин. Точечные оценки параметров распределения. Статистическая проверка гипотез. Сравнение средних известной и неизвестной точности измерений. Критерий Хи-квадрат для проверки гипотезы о виде распределения.

    курсовая работа [79,0 K], добавлен 23.01.2012

Работы в архивах красиво оформлены согласно требованиям ВУЗов и содержат рисунки, диаграммы, формулы и т.д.
PPT, PPTX и PDF-файлы представлены только в архивах.
Рекомендуем скачать работу.