Задача линейного разделения классов Персептрон

ЛАБОРАТОРНАЯ РАБОТА 1

МОДЕЛИ НЕЙРОНОВ И МЕТОДЫ ИХ ОБУЧЕНИЯ

Задача линейного разделения классов

Персептрон

Простой персептрон – это модель МакКаллока-Питса (рис.1). Весовые коэффициенты входов сумматора, на которые подаются входные сигналы , обозначаются , а пороговое значение - . Нелинейная функция активации персептрона является ступенчатой, вследствие чего выходной сигнал нейрона может принимать только два значения – 0 и 1 в соответствии с правилом

или

где обозначает выходной сигнал сумматора

.

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

Рис.1. Нейрон МакКаллока-Питса

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

Линейное разделение классов состоит в построении линейного решающего правила, т.е. такого вектора

- порог, что при

вектор

относится к первому классу, а при

- ко второму.

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

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

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

где - номер предыдущего цикла,  коэффициент обучения.

4. Если , а , то значения весов уточняются по формуле

.

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

1.2. Сигмоидальный нейрон

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

тогда как биполярная функция задается в виде

Параметр

влияет на крутизну графика функции

. При

сигмоидальная функция превращается в функцию ступенчатого типа, идентичную функции активации персептрона. На практике чаще всего используется значение

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

тогда как для биполярной функции

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

проводится в направлении отрицательного градиента целевой функции

, где

. Компонента градиента имеет вид

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

Значения весовых коэффициентов уточняются по формуле

где .

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

,

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

1.3. Модель нейрона Хебба
Структурная схема нейрона Хебба соответствует стандартной форме модели нейрона (рис.1). Д.Хебб предложил формальное правило, в котором вес

связи нейрона

с нейроном

изменяется пропорционально произведению их выходных сигналов

где - коэффициент обучения.

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

.

В каждом цикле обучения происходит суммирование текущего значения веса и его приращения :

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

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

2. Аппроксимация функций

2.1. Аппроксимация линейных функций. Нейрон типа «адалайн»

В нейроне типа «адалайн» (ADAptive LInear Neuron – адаптивный линейный нейрон) адаптивный подбор весовых коэффициентов осуществляется в процессе минимизации квадратичной ошибки, определяемой как

В связи с выполнением условия дифференцируемости целевой функции стало возможным применение алгоритма градиентного обучения. Значения весовых коэффициентов уточняются следующим способом

_.

2.2. Паде-нейрон
Паде-нейрон вычисляет произвольную дробно-линейную функцию вектора

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

Паде-нейрон может использоваться как обобщение нейрона типа «адалайн» в тех случаях, когда линейных функций становится недостаточно, в частности в задачах интерполяции эмпирических зависимостей.

В случае Паде-нейрона квадратичная ошибка определяется как

и значения весовых коэффициентов уточняются по следующим формулам

2.3. Нейрон с квадратичным сумматором

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

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

если , то объект принадлежит первому классу;

если , то объект принадлежит второму классу (при условии правильного выбора коэффициентов ).

Квадратичная ошибка здесь определяется как

.

Коэффициенты квадратичного сумматора уточняются по формулам

Недостаток такого классификатора – большое число настраиваемых параметров.

3. Нейроны типа WTA. Самообучение нейронов.
Нейроны типа WTA (Winner Takes All – Победитель получает все) имеют входной модуль в виде адаптивного сумматора. Выходной сигнал

-го сумматора определяется по формуле

По результатам сравнения сигналов отдельных нейронов победителем признается нейрон, у которого оказалось наибольшим. Нейрон-победитель вырабатывает на своем выходе состояние 1, а остальные (проигравшие) нейроны переходят в состояние 0.

Для обучения нейронов WTA учитель не требуется. На начальном этапе случайным образом выбираются весовые коэффициенты каждого нейрона, нормализуемые относительно 1 по формуле

.

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

определяется победитель этапа. Победитель переходит в состояние 1, что позволяет произвести уточнение весов его входных линий по правилу

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

Выходной сигнал -го нейрона может быть описан векторным отношением

.

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

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

Ниже представлена программа на языке системы математического моделирования MATLAB, описывающая решение простейшей задачи разбиения векторов на два класса методом центров масс. В качестве входных векторов нейрона рассматриваются двоичные векторы размерности 63, полученные из изображений букв А и В размером пикселей. Согласно методу центров масс, весовой вектор нейрона, различающего эти две буквы, вычисляется по формуле

где - вектор, задающий букву А, - вектор, задающий букву B. В результате, для буквы А нейрон дает выходной сигнал, равный 1, а для буквы B – выходной сигнал, равный –1.

function center

% СЧИТЫВАНИЕ ИЗОБРАЖЕНИЯ БУКВЫ А

A=imread('A.bmp');

% ВЫВОД ИЗОБРАЖЕНИЯ

subplot(1,2,1);

imshow(A)

% СЧИТЫВАНИЕ ИЗОБРАЖЕНИЯ БУКВЫ B

B=imread('B.bmp');

% ВЫВОД ИЗОБРАЖЕНИЯ

subplot(1,2,2);

imshow(B)

% РАСЧЕТ ВЕСОВЫХ КОЭФФИЦИЕНТОВ МЕТОДОМ ЦЕНТРОВ МАСС

W=(double(A)-double(B))/2;

% ТЕСТИРОВАНИЕ НЕЙРОНА

S=sum(sum(W.*double(A)));

yA=sgn(S)

S=sum(sum(W.*double(B)));

yB=sgn(S)

%-----------------------------------------------------

% ФУНКЦИЯ АКТИВАЦИИ НЕЙРОНА

function y=sgn(x)

if (x>=0)

y=1;

else

y=-1;

end
Пример 2. Обучение персептрона распознаванию букв (класс1) и цифр (класс 2)

function perceptron

% СЧИТЫВАНИЕ ИЗОБРАЖЕНИЙ

L(:,:,1)=imread('A.bmp');

L(:,:,2)=imread('B.bmp');

L(:,:,3)=imread('C.bmp');

L(:,:,4)=imread('D.bmp');

D(:,:,1)=imread('0.bmp');

D(:,:,2)=imread('1.bmp');

D(:,:,3)=imread('2.bmp');

D(:,:,4)=imread('3.bmp');

%ВЫВОД ИЗОБРАЖЕНИЙ

subplot(2,4,1); imshow(L(:,:,1));

subplot(2,4,2); imshow(L(:,:,2));

subplot(2,4,3); imshow(L(:,:,3));

subplot(2,4,4); imshow(L(:,:,4));

subplot(2,4,5); imshow(D(:,:,1));

subplot(2,4,6); imshow(D(:,:,2));

subplot(2,4,7); imshow(D(:,:,3));

subplot(2,4,8); imshow(D(:,:,4));

%МЕТОД ЦЕНТРОВ МАСС

% расчет вектора весовых коэффициентов

W=0;

for i=1:4

W=W+double(L(:,:,i))-double(D(:,:,i));

end

W=W/8

% тестирование нейрона на изображениях букв

for i=1:4

S=sum(sum(W.*double(L(:,:,i))));

yL(i)=sgn(S);

end

yL

% тестирование нейрона на изображениях цифр

for i=1:4

S=sum(sum(W.*double(D(:,:,i))));

yD(i)=sgn(S);

end

%АЛГОРИТМ ОБУЧЕНИЕ ПЕРСЕПТРОНА

training=1;

epoch=0; % число итераций

while (training > 0)

% этап обучения

training=0;

epoch=epoch+1;

alpha=1/epoch; % коэффициент обучения

for i=1:4

delta=(1-yL(i));

if (delta~=0)

training=1;

W=W+alpha*delta*double(L(:,:,i));

end

delta=(-1-yD(i));

if (delta~=0)

training=1;

W=W+alpha*delta*double(D(:,:,i));

end

W % вывод весовых коэффициентов

%этап тестирования

% тестирование на буквах

for i=1:4

S=sum(sum(W.*double(L(:,:,i))));

yL(i)=sgn(S);

end

% тестирование на цифрах

for i=1:4

S=sum(sum(W.*double(D(:,:,i))));

yD(i)=sgn(S);

end

end % of while

epoch % число итераций

%-------------------------------------------

% ФУНКЦИЯ АКТИВАЦИИ НЕЙРОНА

function y=sgn(x)

if (x>0)

y=1;

else

y=-1;

end

ВАРИАНТЫ ЗАДАНИЙ

В системе моделирования MATLAB реализовать следующие модели нейронов:

Персептрон. Решить задачу линейного разделения входных данных на два класса (класс букв и класс цифр) обучением по всему задачнику при значениях коэффициента обучения , где  номер итерации.
Сигмоидальный нейрон с униполярной функцией активации. Методом наискорейшего спуска решить задачу линейного разделения входных данных на два класса: класс букв и класс цифр. Обучать по отдельным примерам.
Сигмоидальный нейрон с униполярной функцией активации. Методом наискорейшего спуска решить задачу линейного разделения входных данных на два класса: класс букв и класс цифр. Обучать по всему задачнику.
Сигмоидальный нейрон с биполярной функцией активации. Методом обучения с моментом решить задачу линейного разделения классов. Обучать по отдельным примерам.
Сигмоидальный нейрон с биполярной функцией активации. Методом обучения с моментом решить задачу линейного разделения классов. Обучать по всему задачнику.
Нейрон Хебба. Решить задачу линейного разделения классов. Обучать по отдельным примерам.
Нейрон Хебба. Решить задачу линейного разделения классов. Обучать по всему задачнику.
Пусть для некоторого набора действительных чисел определены значения функции . Требуется построить наилучшую линейную аппроксимацию функции , т.е. найти такой вектор , чтобы приближение давало наименьшую ошибку. Решить задачу аппроксимации линейной функции нейроном типа «адалайн».
Паде-нейрон. Решить задачу аппроксимации функции .
Паде-нейрон. Решить задачу аппроксимации функции .
Нейрон с квадратичным сумматором. Решить задачу аппроксимации квадратичной функции .
На плоскости задать два класса объектов с нормальным распределением. Используя нейрон с квадратичным сумматором построить на плоскости оптимальную кривую, разделяющие классы.
Самообучение распознаванию объектов на изображении. Задать наборы изображений объектов «прямоугольник», «треугольник», «крест» различных размеров. Обучить слой нейронов типа WTA распознаванию этих объектов.