Это не официальный сайт wikipedia.org 01.01.2023

Нейронная сеть Кохонена — Википедия

Нейронная сеть Кохонена

(перенаправлено с «Векторное квантование»)

Нейронные сети Кохонена — класс нейронных сетей, основным элементом которых является слой Кохонена. Слой Кохонена состоит из адаптивных линейных сумматоров («линейных формальных нейронов»). Как правило, выходные сигналы слоя Кохонена обрабатываются по правилу «Победитель получает всё»: наибольший сигнал превращается в единичный, остальные обращаются в ноль.

По способам настройки входных весов сумматоров и по решаемым задачам различают много разновидностей сетей Кохонена[1]. Наиболее известные из них:

Слой КохоненаПравить

Базовая версияПравить

Слой Кохонена состоит из некоторого количества n   параллельно действующих линейных элементов. Все они имеют одинаковое число входов m   и получают на свои входы один и тот же вектор входных сигналов x = ( x 1 , . . . x m )  . На выходе j  го линейного элемента получаем сигнал

y j = w j 0 + i = 1 m w j i x i ,  

где:

  • w j i   — весовой коэффициент i  -го входа j  -го нейрона;
  • i   — номер входа;
  • j   — номер нейрона;
  • w j 0   — пороговый коэффициент.

После прохождения слоя линейных элементов сигналы посылаются на обработку по правилу «победитель забирает всё»: среди выходных сигналов выполняется поиск максимального y j  ; его номер j max = a r g max j { y j }  . Окончательно, на выходе сигнал с номером j max   равен единице, остальные — нулю. Если максимум одновременно достигается для нескольких j max  , то:

  • либо принимают все соответствующие сигналы равными единице;
  • либо равным единице принимают только первый сигнал в списке (по соглашению).

«Нейроны Кохонена можно воспринимать как набор электрических лампочек, так что для любого входного вектора загорается одна из них»[5].

Геометрическая интерпретацияПравить

 
Разбиение плоскости на многоугольники Вороного-Дирихле для случайно выбранных точек (каждая точка указана в своём многоугольнике).

Большое распространение получили слои Кохонена, построенные следующим образом: каждому ( j  -му) нейрону сопоставляется точка W j = ( w j 1 , . . . w j m )   в m  -мерном пространстве (пространстве сигналов). Для входного вектора x = ( x 1 , . . . x m )   вычисляются его евклидовы расстояния ρ j ( x )   до точек W j   и «ближайший получает всё» — тот нейрон, для которого это расстояние минимально, выдаёт единицу, остальные — нули. Следует заметить, что для сравнения расстояний достаточно вычислять линейную функцию сигнала:

ρ j ( x ) 2 = x W j 2 = W j 2 2 i = 1 m w j i x i + x 2  

(здесь y   — евклидова длина вектора: y 2 = i y i 2  ). Последнее слагаемое x 2   одинаково для всех нейронов, поэтому для нахождения ближайшей точки оно не нужно. Задача сводится к поиску номера наибольшего из значений линейных функций:

j max = a r g max j { i = 1 m w j i x i 1 2 W j 2 } .  

Таким образом, координаты точки W j = ( w j 1 , . . . w j m )   совпадают с весами линейного нейрона слоя Кохонена (при этом значение порогового коэффициента w j 0 = W j 2 / 2  ).

Если заданы точки W j = ( w j 1 , . . . w j m )  , то m  -мерное пространство разбивается на соответствующие многогранники Вороного-Дирихле V j  : многогранник V j   состоит из точек, которые ближе к W j  , чем к другим W k   ( k j  )[6].

Сети векторного квантованияПравить

Задача векторного квантования с k   кодовыми векторами W j   для заданной совокупности входных векторов S   ставится как задача минимизации искажения при кодировании, то есть при замещении каждого вектора из S   соответствующим кодовым вектором. В базовом варианте сетей Кохонена используется метод наименьших квадратов и искажение D   вычисляется по формуле

D = j = 1 k x K j x W j 2 ,  

где K j   состоит из тех точек x S  , которые ближе к W j  , чем к другим W l   ( l j  ). Другими словами, K j   состоит из тех точек x S  , которые кодируются кодовым вектором W j  .

Если совокупность S   задана и хранится в памяти, то стандартным выбором в обучении соответствующей сети Кохонена является метод K-средних. Это метод расщепления:

  • при данном выборе кодовых векторов (они же весовые векторы сети) W j   минимизацией D   находим множества K j   — они состоят из тех точек x S  , которые ближе к W j  , чем к другим W l  ;
  • при данном разбиении S   на множества K j   минимизацией D   находим оптимальные позиции кодовых векторов W j   — для оценки по методу наименьших квадратов это просто средние арифметические:
W j = 1 | K j | x K j x ,  

где | K j |   — число элементов в K j  .

Далее итерируем. Этот метод расщепления сходится за конечное число шагов и даёт локальный минимум искажения.

Если же, например, совокупность S   заранее не задана, или по каким-либо причинам не хранится в памяти, то широко используется онлайн метод. Векторы входных сигналов x   обрабатываются по одному, для каждого из них находится ближайший кодовый вектор («победитель», который «забирает всё») W j ( x )  . После этого данный кодовый вектор пересчитывается по формуле

W j ( x ) n e w = W j ( x ) o l d ( 1 θ ) + x θ ,  

где θ ( 0 , 1 )   — шаг обучения. Остальные кодовые векторы на этом шаге не изменяются.

Для обеспечения стабильности используется онлайн метод с затухающей скоростью обучения: если T   — количество шагов обучения, то полагают θ = θ ( T )  . Функцию θ ( T ) > 0   выбирают таким образом, чтобы θ ( T ) 0   монотонно при T   и чтобы ряд T = 1 θ ( T )   расходился, например, θ ( T ) = θ 0 / T  .

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

Самоорганизующиеся карты КохоненаПравить

Идея и алгоритм обученияПравить

Задача векторного квантования состоит, по своему существу, в наилучшей аппроксимации всей совокупности векторов данных k   кодовыми векторами W j  . Самоорганизующиеся карты Кохонена также аппроксимируют данные, однако при наличии дополнительной структуры в совокупности кодовых векторов (англ. codebook). Предполагается, что априори задана некоторая симметричная таблица «мер соседства» (или «мер близости») узлов: для каждой пары j , l   ( j , l = 1 , . . . k  ) определено число η j l   ( 0 η j l 1  ) при этом диагональные элементы таблицы близости равны единице ( η j j = 1  ).

Векторы входных сигналов x   обрабатываются по одному, для каждого из них находится ближайший кодовый вектор («победитель», который «забирает всё») W j ( x )  . После этого все кодовые векторы W l  , для которых η j ( x ) l 0  , пересчитываются по формуле

W l n e w = W l o l d ( 1 η j ( x ) l θ ) + x η j ( x ) l θ ,  

где θ ( 0 , 1 )   — шаг обучения. Соседи кодового вектора — победителя (по априорно заданной таблице близости) сдвигаются в ту же сторону, что и этот вектор, пропорционально мере близости.

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

Самоорганизующиеся карты Кохонена служат, в первую очередь, для визуализации и первоначального («разведывательного») анализа данных[7]. Каждая точка данных отображается соответствующим кодовым вектором из решётки. Так получают представление данных на плоскости («карту данных»). На этой карте возможно отображение многих слоёв: количество данных, попадающих в узлы (то есть «плотность данных»), различные функции данных и так далее. При отображении этих слоёв полезен аппарат географических информационных систем (ГИС). В ГИС подложкой для изображения информационных слоев служит географическая карта. Карта данных является подложкой для произвольного по своей природе набора данных. Карта данных служит заменой географической карте там, где географической карты просто не существует. Принципиальное отличие в следующем: на географической карте соседние объекты обладают близкими географическими координатами, на карте данных близкие объекты обладают близкими свойствами. С помощью карты данных можно визуализировать данные, одновременно нанося на подложку сопровождающую информацию (подписи, аннотации, атрибуты, информационные раскраски)[7]. Карта служит также информационной моделью данных. С её помощью можно заполнять пробелы в данных. Эта способность используется, например, для решения задач прогнозирования.

Самоорганизующиеся карты и главные многообразияПравить

Идея самоорганизующихся карт очень привлекательна и породила массу обобщений, однако, строго говоря, мы не знаем, что мы строим: карта — это результат работы алгоритма и не имеет отдельного («объектного») определения. Есть, однако, близкая теоретическая идея — главные многообразия (англ. principal manifolds)[8]. Эти многообразия обобщают линейные главные компоненты. Они были введены как линии или поверхности, проходящие через «середину» распределения данных, с помощью условия самосогласованности: каждая точка x   на главном многообразии M   является условным математическим ожиданием тех векторов z  , которые проектируются на x   (при условии x = P ( z )  , где P   — оператор проектирования окрестности M   на M  ),

x = E ( z | P ( z ) = x ) .  

Самоорганизующиеся карты могут рассматриваться как аппроксимации главных многообразий и популярны в этом качестве[9].

Упругие картыПравить

 
Визуализация набора данных по экспрессии генов в раке молочной железы с использованием упругих карт (b) и метода главных компонент (c). Классы точек показаны с использованием размера (ER — статус эстроген-рецептора), формы (GROUP — риск развития метастаз) и цвета (TYPE — молекулярный тип опухоли). На панели (a) показана конфигурация узлов двумерной упругой карты в проекции на первые три главные компоненты. Сравнивая (b) и (c), можно заметить, что базальный тип опухоли как кластер лучше отделён на нелинейной проекции (b).

Метод аппроксимации многомерных данных, основанный на минимизации «энергии упругой деформации» карты, погружённой в пространство данных, был предложен А. Н. Горбанём в 1996 году, и впоследствии развит им совместно с А. Ю. Зиновьевым, А. А. Россиевым и А. А. Питенко[7]. Метод основан на аналогии между главным многообразием и эластичной мембраной и упругой пластиной. В этом смысле он является развитием классической идеи сплайна (хотя упругие карты и не являются многомерными сплайнами).

Пусть задана совокупность входных векторов S  . Так же, как и сети векторного квантования и самоорганизующиеся карты, упругая карта представлена как совокупность кодовых векторов (узлов) W j   в пространстве сигналов. Множество данных S   разделено на классы K j  , состоящие из тех точек x S  , которые ближе к W j  , чем к другим W l   ( l j  ). Искажение кодирования D  

D = j = 1 k x K j x W j 2 ,  

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

На множестве узлов задана дополнительная структура: некоторые пары связаны «упругими связями», а некоторые тройки объединены в «рёбра жёсткости». Обозначим множество пар, связанных упругими связями, через E  , а множество троек, составляющих рёбра жёсткости, через G  . Например, в квадратной решётке ближайшие узлы (как по вертикали, так и погоризонтали) связываются упругими связями, а ребра жёсткости образуются вертикальными и горизонтальными тройками ближайших узлов. Энергия деформации карты состоит из двух слагаемых:

энергия растяжения U E = λ ( W i , W j ) E W i W j 2 ;  
энергия изгиба U G = μ ( W i , W j , W l ) G W i 2 W j + W l 2 ;  

где λ , μ   — соответствующие модули упругости.

Задача построения упругой карты состоит в минимизации функционала

U = D + U E + U G ;  

Если разбиение совокупности входных векторов S   на классы K j   фиксировано, то минимизация U   — линейная задача с разреженной матрицей коэффициентов. Поэтому, как и для сетей векторного квантования, применяется метод расщепления: фиксируем { W j }   — ищем { K j }   — для данных { K j }   ищем { W j }   — для данных { W j }   ищем { K j }   — … Алгоритм сходится к (локальному) минимуму U  .

Метод упругих карт позволяет решать все задачи, которые решают самоорганизующиеся карты Кохонена, однако обладает большей регулярностью и предсказуемостью. При увеличении модуля изгиба μ   упругие карты приближаются к линейным главным компонентам. При уменьшении обоих модулей упругости они превращаются в Кохоненовские сети векторного квантования. В настоящее время упругие карты интенсивно используются для анализа многомерных данных в биоинформатике.[10] Соответствующее программное обеспечение опубликовано и свободно доступно на сайте института Кюри (Париж)[11][12].

На рисунке представлены результаты визуализации данных по раку молочной железы. Эти данные содержат 286 примеров с указанием уровня экспрессии 17816 генов[13]. Они доступны онлайн как ставший классическим тестовый пример для визуализации и картографии данных[14].

Сети векторного квантования, обучаемые с учителемПравить

 
Пример возможного разделения классов, составленного с помощью разбиения Вороного-Дирихле.

Решается задача классификации. Число классов может быть любым. Изложим алгоритм для двух классов, A   и B  . Исходно для обучения системы поступают данные, класс которых известен. Задача: найти для класса A   некоторое количество k A   кодовых векторов W j A  , а для класса B   некоторое (возможно другое) количество k B   кодовых векторов W l B   таким образом, чтобы итоговая сеть Кохонена с k A + k B   кодовыми векторами W j A  , W l B   (объединяем оба семейства) осуществляла классификацию по следующему решающему правилу:

если для вектора входных сигналов x   ближайший кодовый вектор («победитель», который в слое Кохонена «забирает всё») принадлежит семейству { W j A }  , то x   принадлежит классу A  ; если же ближайший к x   кодовый вектор принадлежит семейству { W l B }  , то x   принадлежит классу B  .

С каждым кодовым вектором объединённого семейства { W j A } { W l B }   связан многогранник Вороного-Дирихле. Обозначим эти многогранники V j A  , V l B   соответственно. Класс A   в пространстве сигналов, согласно решающему правилу, соответствует объединению j V j A  , а класс B   соответствует объединению l V l B  . Геометрия таких объединений многогранников может быть весьма сложной (см. рисунок с примером возможного разбиения на классы).

Правила обучения сети онлайн строится на основе базового правила обучения сети векторного квантования. Пусть на вход системы подаётся вектор сигналов x  , класс которого известен. Если он классифицируется системой правильно, то соответствующий x   кодовый вектор W  слегка сдвигается в сторону вектора сигнала («поощрение»)

W n e w = W o l d ( 1 θ ) + x θ ,  

Если же x   классифицируется неправильно, то соответствующий x   кодовый вектор W  слегка сдвигается в противоположную сторону от сигнала («наказание»)

W n e w = W o l d ( 1 + θ ) x θ ,  

где θ ( 0 , 1 )   — шаг обучения. Для обеспечения стабильности используется онлайн метод с затухающей скоростью обучения. Возможно также использование разных шагов для «поощрения» правильного решения и для «наказания» неправильного.

Это — простейшая (базовая) версия метода[15]. Существует множество других модификаций.

ПримечанияПравить

  1. How many kinds of Kohonen networks exist? Internet FAQ Archives. Online Education  (неопр.). Дата обращения: 31 августа 2008. Архивировано 11 мая 2008 года.
  2. Hecht-Nielsen, R. (1990), Neurocomputing, Reading, MA: Addison-Wesley, ISBN 0-201-09355-3.
  3. Kohonen, T. (1989/1997/2001), Self-Organizing Maps, Berlin — New York: Springer-Verlag. First edition 1989, second edition 1997, third extended edition 2001, ISBN 0-387-51387-6, ISBN 3-540-67921-9
  4. Kohonen, T. (1988), Learning Vector Quantization, Neural Networks, 1 (suppl 1), 303.
  5. Уоссермен, Ф. Нейрокомпьютерная техника: Теория и практика = Neural Computing. Theory and Practice. — М.: Мир, 1992. — 240 с. — ISBN 5-03-002115-9. Архивированная копия  (неопр.). Дата обращения: 1 сентября 2008. Архивировано из оригинала 30 июня 2009 года.
  6. Real time interactive Voronoi and Delaunay diagrams with source code  (неопр.). Дата обращения: 1 сентября 2008. Архивировано 1 сентября 2008 года.
  7. 1 2 3 Зиновьев А. Ю. Визуализация многомерных данных. — Красноярск: Изд. Красноярского государственного технического университета, 2000. — 180 с.
  8. Диссертация T. Хасти: Hastie T., Principal curves and surfaces Архивная копия от 21 февраля 2017 на Wayback Machine, Ph.D dissertation, Stanford Linear accelerator center, Stanford University, Stanford, California, US, November 1984. А также на сайте PCA Архивная копия от 7 ноября 2018 на Wayback Machine. С этой работы началось изучение главных многообразий.
  9. Yin H. Learning nonlinear principal manifolds by self-organising maps Архивная копия от 6 марта 2019 на Wayback Machine, In: Gorban A. N. et al (Eds.), LNCSE 58, Springer, 2007. ISBN 978-3-540-73749-0
  10. Gorban A. N., Kegl B., Wunsch D., Zinovyev A. Y. (Eds.), Principal Manifolds for Data Visualisation and Dimension Reduction, Series: Lecture Notes in Computational Science and Engineering 58, Springer, Berlin — Heidelberg — New York, 2007, XXIV, 340 p. 82 illus. ISBN 978-3-540-73749-0 (а также онлайн Архивная копия от 16 марта 2019 на Wayback Machine).
  11. VIMIDA: a Java-applet for VIsualization of MIcroarray DAta  (неопр.). Дата обращения: 6 сентября 2008. Архивировано 9 октября 2008 года.
  12. ViDaExpert: a software for multidimensional vectorial data visualization  (неопр.). Дата обращения: 6 сентября 2008. Архивировано 26 апреля 2012 года.
  13. Wang Y., Klijn J.G., Zhang Y., Sieuwerts A.M., Look M.P., Yang F., Talantov D., Timmermans M., Meijer-van Gelder M.E., Yu J. et al. Gene-expression profiles to predict distant metastasis of lymph-node-negative primary breast cancer. Lancet 365 (2005), 671—679.
  14. Principal manifolds for data cartography and dimension reduction, Leicester, UK, August 2006. A web-page with test microarrays datasets provided for participants of the workshop Архивная копия от 24 сентября 2008 на Wayback Machine.
  15. DLVQ Fundamentals  (неопр.). Дата обращения: 7 ноября 2018. Архивировано 19 декабря 2018 года.

См. такжеПравить