Тема: «Моделирование нечеткой системы инструментами нечеткой логики».

Цель работы: изучить метод построения нечеткой системы инструментами нечеткой логики.

В составе Matlab присутствуют пять основных средств графического интерфейса пользователя (ГИП), которые обеспечивают доступ к инструментам нечеткой логики: редакторы системы нечеткого вывода (СНВ), функции принадлежности, правил вывода, а также средства просмотра правил и поверхности вывода. Эти средства связаны между собой динамически и производимые изменения в одном из них влекут изменения в других.

Редактор СНВ предоставляет возможность формирования проектируемой системы на высоком уровне абстракции: количество входных и выходных переменных, наименование переменных.

Редактор функций принадлежности используется для определения формы функций принадлежности, ассоциированных с каждой переменной.

Редактор правил вывода применяется для редактирования списка правил, которые определяют поведение проектируемой системы.

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

Средство просмотра поверхности вывода используется для отображения зависимости выхода системы от одного или двух входов. Другими словами, оно генерирует и выводит карту поверхности вывода разработанной СНВ.

Редактор СНВ. Построение нечетких систем по Мамдани. Для построения создаваемой системы в командной строке основного окна Matlab необходимо набрать команду fuzzy. Окно редактора новой СНВ содержит входную, обозначенную input 1 и выходную — output 1 переменные. По умолчанию инструмент нечеткой логики предлагает создавать СНВ типа Мамдани.

Для того чтобы добавить новую переменную, необходимо выбрать в меню Edit соответствующий пункт (для входной переменной — Add input, для выходной — Add output). Изменение наименования переменной происходит по шагам.

Шаг 1. Отмечается переменная, которую необходимо переименовать.

Шаг 2. В поле редактирования изменяется наименование переменной по умолчанию на имя, предложенное пользователем.

Сохранение проектируемой системы в рабочее пространство среды MATLAB (в переменную) производится с помощью пункта меню File — Save to workspace as... В этом случае данные сохраняются до окончания сеанса работы с Matlab. Для сохранения данных на диске после окончания сеанса работы применяется соответствующий пункт того же меню — Save to disk as...

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

выбором в меню View пункта Edit Membership Functions...;

двойным щелчком мыши на изображении соответствующей переменной (входной или выходной);

набором в командной строке оператора mfedit.

С помощью редактора функций принадлежности можно отображать и редактировать любые функции принадлежности, ассоциированные (связанные) со всеми входными и выходными переменными разрабатываемой СНВ.

Связывание функции принадлежности с именем переменной происходит

следующим образом:

выбирается переменная по имени из набора графических объектов окна редактора функции принадлежности;

указывается диапазон изменения значений для базовой и видимый диапазон для текущей переменных;

в меню Edit выбирается пункт Add MFs... В появившемся окне выбирают вид функций принадлежности и их количество.

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

Таким образом, при построении СНВ необходимо с помощью редактора функций принадлежности определить соответствующие функции для каждой из входных и выходных переменных.

Редактор правил вывода. После того как указано количество входных и

выходных переменных, определены их наименования и построены соответствующие функции принадлежности, в СНВ необходимо включить правила вывода. Для этого в меню View выбирается пункт Edit Rules... или в командной строке Matlab набирается команда ruleedit.

Основываясь на описаниях входных и выходных переменных, определенных в редакторе функций принадлежности, редактор правил вывода формирует структуру правила автоматически. От пользователя требуется лишь связать значения входных и выходных переменных, выбирая из списка заданных ранее функций принадлежности и определить логические связки между ними. Также допускается использование логического отрицания (НЕ) и изменение весов правил в диапазоне от 0 до 1.

Правила вывода могут отображаться в окне в различных форматах, которые определяются путем выбора соответствующего пункта подменю Format меню Options. По умолчанию используется расширенный формат отображения правил вывода (verbose form):

If (input_lis [not]mf_lj1) <and, or>...

(input_i is [not ]mf_ij...<and, or>

(input_n is [not]mf_njn) then

(output_l is [not]mf_n + 1 j n + 1<and,or>...

(output_k is [not]mf_k + njk + n) <and,or>...

(output_m is[not]mf_m + njm + n) (w),

где i — номер входной переменной;

ji — номер функции принадлежности i-й переменной;

k— номер выходной переменной;

n — количество входных переменных;

т— количество выходных переменных;

w — вес правила.

Круглые скобки заключают в себе обязательные параметры, квадратные — необязательные, а угловые — альтернативные параметры (один на выбор).

Кроме формата по умолчанию, существуют еще два вида форматов отображения правил: символьный (symbolic form) и индексный (indexed form). Символьный формат имеет следующий вид:

(input_l<~=,==>mf_lj1) <&,| >...

(input_i<~=,==>mf_ij1)...<&, |>

(input_n<~=,==>mf_nj п) =>

(output_l<~=,==>mf__n + ljn +1)...<&, | >

(output_k<~=,==>mf__k + njk + n)<&,|>...

(output__m<~=,==>mf_m + njm n) (w)

Отличие символьного формата от расширенного состоит в том, что вместо словесной интерпретации связок используется символьная (символы «&» и «|» — соответственно определяют логическое И и логическое ИЛИ, символ «~» — логическое отрицание, а символ «=>» является разделителем условной и заключительной частей правила (антецендента и консеквента).Общее описание правила вывода в индексном формате может быть представлено в следующем виде: [-]1j1…[-]iji…[-]njn+1jn+1..[-]k+njk+1…[-]m+njm+n(w):<1,2>.

Здесь порядок следования чисел соответствует очередности вводимых переменных, причем символ «,» разделяет правило на условную и заключительную части. До двоеточия записывается порядковый номер соответствующей функции принадлежности, после двоеточия — вид логической связки («1» — логическое И, «2» — логическое ИЛИ). Логическое отрицание задается символом «—».

После определения правил вывода в одноименном редакторе можно утверждать, что СНВ полностью создана.

Пример 10. Создание СНВ. Рассмотрим следующую ситуацию. Необходимо оценить степень инвестиционной привлекательности конкретного бизнес-проекта на основании данных о ставке дисконтирования и периоде окупаемости.

Шаг 1. Вызываем редактор для создания СНВ, набирая в командной строке fuzzy. Добавляем входную переменную посредством выбора в меню Edit пункта Add input. В результате получаем следующую структуру СНВ: два входа, механизм нечеткого вывода по Мамдани, один выход. Объявляем первую переменную как discont, а вторую — period, которые соответственно будут представлять ставку дисконтирования и период окупаемости бизнес-проекта. Наименование выходной переменной, на основании которой принимается решение о степени инвестиционной привлекательности бизнес-проекта, задается как rate. Сохраним создаваемую модель под именем Invest.

Шаг 2. Каждой входной и выходной переменной поставим в соответствие набор функций принадлежности. Данная процедура реализуется в редакторе функций принадлежности. Для discont определяем диапазон базовой переменной от 5 до 50 (единица измерения — проценты). Такой же диапазон выбираем для ее отображения. Добавим три функции принадлежности, тип которых — trimf Последовательно выделяя мышью отдельные функции принадлежности, присвоим наименования — small, middle, big соответственно небольшой, средней и большой ставке дисконтирования. Переменной period диапазон базовой переменной определен равным [3, 36] (единица измерения — месяцы), поставлены в соответствие три функции принадлежности типа gaussmf с наименованиями — short, normal, long. Таким образом, переменная срока окупаемости бизнес-проекта будет принимать следующие значения: короткий, обычный и длительный срок окупаемости. Наконец, для переменной rate определяем: базовая переменная изменяет значение в диапазоне [0, 1], семантика описывается тремя функций принадлежности типа trimf с наименованиями: bad, normal, good.

Шаг 3. Заключительным этапом построения СНВ является определение набора правил, которые задают связь входных переменных с выходными. Для этого в редакторе правил вывода определим:

ЕСЛИ discont = small И period = short ТО rate = good ЕСЛИ discont — HE small И period = long TO rate = bad ЕСЛИ discont = middle И period = normal TO rate = normal ЕСЛИ discont — big И period = short TO rate = normal

Рисунок 10 – Окно редактора системы нечеткого вывода

Рисунок 11 – Окно редактора функций принадлежности

В расширенном формате отображения указанные правила вывода представляются следующим образом:

Рисунок 12 – Окно редактора правил вывода

if(discont is small) and (period is short) then (rate is good) (1)

if (discont is not small) and (period is long) then (rate is bad) (1)

if(discont is middle) and (period is normal) then (rate is normal) (1)

if(discont is big) and (period is short) then (rate is normal) (1)

При изменении формата на символьный правила вывода будут иметь вид:

(discont == small) & {period == short) => (rate == good) (1)

(discont ~= small) & {period —= long) => (rate == bad) (1)

{discont == middle) & (period == normal) => (rate == normal) (1) (discont ==

big) & (period == short) => (rate == normal) (1)

Наконец, то же самое, но в индексном формате:

 11,3(1): 1 -13,1(1): 1 22,2(1): 1 31,2(1): 1.

Средство просмотра правил вывода. Данное средство просмотра правил вывода позволяет отобразить процесс нечеткого вывода и получить результат. Главное окно средства просмотра состоит из нескольких графических окон, располагаемых по строкам и столбцам. Количество строк соответствует числу правил нечеткого вывода, а количество столбцов — числу входных и выходных переменных, заданных в разрабатываемой СНВ. Дополнительное графическое окно служит для отображения результата нечеткого вывода и операции де- фазификации. В каждом окне отображается соответствующая функция принадлежности, уровень ее среза (для входных переменных) и вклад отдельной функции принадлежности в общий результат (для выходных переменных).

В нижней части главного окна можно отобразить номера правил вывода в различных форматах вывода, отмечая их мышью. Для изменения формата в меню Options выбирается пункт Rule display format.

Изменение значений входных переменных допустимо двумя способами:

путем ввода в поле Input записи входного вектора, размерность которого равна количеству входных переменных;

щелчком мыши в любом графическом окне, которое относится к входной переменной.

Входной вектор в каждом из этих вариантов определения исходных данных будет задавать набор красных вертикальных прямых.

Для СНВ, рассмотренной в примере П10, при входном векторе [15 10] (ставка дисконтирования — 15 %, период окупаемости бизнес-проекта — 10 месяцев) результат (степень инвестиционной привлекательности) будет составлять 0,585.

Средство просмотра поверхности вывода. Средство просмотра поверхности вывода позволяет строить трехмерную поверхность как зависимость одной из выходных переменных от двух входных.

Рисунок 13 – Окно средства просмотра правил вывода

Выбор входных и выходных переменных осуществляется посредством ниспадающих меню главного окна рассматриваемого программного средства. Количество выводимых линий по осям Хи Yопределяется в полях ввода Xgrids, Ygrids. Поверхность вывода, соответствующая правилам вывода примера П10.

Построение нечетких систем типа Суджено. Рассмотрим построение СНВ двумя редакторами — СНВ и функций принадлежности. Для построения СНВ типа Суджено необходимо в меню File выбрать пункт New Sugeno FIS. Количество входных и выходных переменных определяется так же, как и при построении СНВ типа Мамдани.

Редактор функций принадлежности. Для СНВ типа Суджено изменения касаются только схемы определения функций принадлежности для выходных переменных. Инструмент нечеткой логики в среде Matlab позволяет разрабатывать два вида нечетких моделей. Первая модель — это нечеткая модель Суджено нулевого порядка. Нечеткое правило вывода имеет следующий вид:

if х is A and у is В then z — k,

где А и В — нечеткие множества антецедента;

к — четко заданная константа консеквента.

Для построения такой модели при добавлении функции принадлежности необходимо выбрать тип — константа (constant) и задать в качестве параметра функции принадлежности численное значение соответствующей константы. Вторая модель — нечеткая модель Суджено первого порядка.

Рисунок 13 ­– Окно просмотра поверхности решений

Рисунок 14 – Окно просмотра правил вывода (вывод по Суджено)

Для нее нечеткое правило вывода записывается следующим образом:

if х as A and у is В then z = р • х + q • у + r, где р, q и r — константы.

В данном случае тип функции принадлежности — линейная зависимость (linear). Для определения параметров функции принадлежности необходимо ввести вектор, элементы которого соответствуют численным значениям констант консеквента.

Работа с редактором правил вывода, а также со средствами просмотра правил и поверхности вывода выполняется аналогично случаю построения СНВ по Мамдани.

Пример нечеткого вывода по Суджено с использованием нечеткой модели нулевого порядка и правил вывода, определенных выше (выходная переменная имеет три значения: bad, normal, good, которые задаются соответственно тремя константами — 0, 0,5, 1).

Вывод: в данной лабораторной работе был изучен метод построения нечеткой системы инструментами нечеткой логики в MATLAB.

Практическая работа №3