Структуры последовательной активности в нейронных сетях со случайными связями Текст научной статьи по специальности «Математика»

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ОПТИМАЛЬНОЕ УПРАВЛЕНИЕ

УДК 621.391.01

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

© 2010 г. ТА. Леванова, МА. Комаров, Е.Ю. Кадина, Г.В. Осипов

Нижегородский госуниверситет им. Н.И. Лобачевского levanova. tatiana@gmail .com

Посмупила в редасцию 15.10.2009

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

Ключевые слова: нейронная сеть, функциональная структура, модель Эрдёша - Реньи, модель Ват-тса - Строгаца, модель Барабаси - Альберта.

Введение

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

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

Активность одиночного нейрона

Для моделирования динамики нейрона была использована система Ходжкина - Хаксли [3]

Cm ~Г~ = gL (VL - Vi) + gNamlhi (VNa - Vi) + dt

+ gKnfV - Vi ) + Iext + lSyn, dm

' =am (1 - mi) -Pmmi , C1)

dt

~7~ = ah (1 - hi) - Phhi

dt

dni n \ n

— = an (1 - ni) -Pnni dt

где Vi - мембранный потенциал /-го нейрона; т Ь, п/ - воротные переменные /-го нейрона, определяющие долю открытых ионных каналов; ах , Рх - некоторые функции, зависящие от V/; gNa, gK, gL - постоянные, определяющие максимальную проводимость ионных каналов; VNa Гк, VI - равновесные потенциалы; 1х - внешний ток, подаваемый на клетки. Используемые параметры таковы, что при отсутствии связей каждый нейрон демонстрирует периодическую спайковую активность, например 1еха = -11. При наличии связей динамика каждого нейрона может быть очень разнообразной и варьируется от генерации периодических спайков до хаотической пачечной активности.

Конфигурация нейронной сети

Была рассмотрена нейронная сеть, состоящая из N элементов. Нейроны связаны ингибиторными синапсами согласно матрице связей G={gij}. Связи моделируются дополнительным внешним синаптическим током, добавляемым в первое уравнение системы (1), в виде

1?п =Еgjr] (0(Егеу - V-) . (2)

-

Здесь коэффициент связи gij обозначает максимальную проводимость канала в синапсе, соединяющем г-й нейрон с у-м нейроном; Егеу -потенциал реверсии. Доля открытых каналов /}(?) в синапсе описывается стандартным кинетическим уравнением dri

^ V- )(1 - Г )-РГ; ,

где

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

Был рассмотрим сначала мотив (малый ансамбль) из трех нейронов, связанных взаимными ингибиторными связями (см. рис. 1а).

Последовательные переключения активности (см. рис.1б) между нейронами в таком ан-

самбле возможны при наличии сильной асимметрии в коэффициентах синаптических связей, направленных по часовой стрелке и против часовой стрелки (£ и h соответственно). Последовательные переключения активности в мотиве из трех нейронов были подробно изучены в [4]. При рассмотрении двух типов связи (слабой

= 0) и сильной (£у =1)) можно утверждать, что последовательная активность в мотиве реализуется при сильной связи, направленной по часовой стрелке, и слабой связи, направленной против часовой стрелки (или наоборот).

Известно, что реальные биологические сети состоят из огромного числа нейронов, между которыми существует огромное число связей. Зададимся вопросом: какова должна быть организация связей между нейронами, чтобы в сети могли возникнуть последовательные переключения активности не только между отдельными нейронами, но также и между группами (кластерами) нейронов? Какое максимальное число мотивов можно получить одновременно в сети заданной размерности? Как зависит это число от архитектуры сети? Ответам на эти вопросы и посвящена данная статья.

Было показано [6, 7], что последовательные переключения пачечной активности в группе (ансамбле) из п нейронов w1,...,wn возникают в том случае, если между любой упорядоченной парой нейронов (кроме пар вида (мг, ^))

имеются симметричные ингибиторные связи, при этом пары вида (^ь ^) связаны

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

Рис. 1. Топология связей в мотиве из трех нейронов - а); последовательные переключения активности в мотиве из трех ингибиторно связанных нейронов (И > g) - б)

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

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

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

Была построена и изучена модель сети нейронов, представляющая собой ориентированный граф А(Ж, Е) (Ж - множество вершин, Е -множество дуг) с N вершинами, причем каждому нейрону из сети ставится в соответствие вершина графа, а каждому аксону /-го нейрона (где і изменяется от 1 до Щ, соединенного синаптическими связями с у-м нейроном, ставится в соответствие дуга (/, /). В графе А(Ж, Е) отсутствуют петли, т.к. нейрон не может быть синаптически связан сам с собой, и мультидуги, т.к. из одного нейрона в другой не может выходить более одного аксона.

Был рассмотрен также граф А'(№, Е'), являющийся подграфом графа А(Ж, Е), такой что:

1. Ж = и ... и ЖпЛ', Ж,' П Ж/ = 0,

/, / = 0...(п - 1);

2. |Ж'| = т, і = 0...(п - 1), т > 1;

3. V і = 0...(п- 1) V wi є Ж,

3 Wi+l є Ж+1 (W(i+1)mod q,Wi) = 0, (1)

4. V ij = 0.. .(п - 1), i ф j, j ф i-1, V w, є Wi,

З wj є Wj? (wi mod q, wjmod q) 1;

5. V w, w' є W', (w, w') = 0.

Таким образом, функциональным структурам в построенном графе будут соответствовать подграфы A', состоящие каждый из n кластеров с т не связанными друг с другом элементами в каждом кластере.

Для построения графа A(W, E) с указанными свойствами были использованы модели Эрдёша -Реньи, Ваттса - Строгаца и Барабаси - Альберта.

Алгоритм поиска функциональных структур

Для решения задачи поиска функциональных структур был выполнен переход к математической модели нейронной сети со случайными связями - графу A( W, E) с описанными выше свойствами.

На основе известного алгоритма поиска компонент связности и поиска мостов в графе [11] был разработан алгоритм, позволяющий в исходном графе, заданном матрицей смежности, отыскивать разбиение исходного графа на максимальное число неизолированных подграфов A' за время, существенно меньшее, чем отводимое на решение этой же задачи методом полного перебора. Известно, что полный пере-

г<пт г

бор занимает время , тогда как разрабо-

танный алгоритм требует время O(M+N) (где M - число ребер в графе, N - число вершин) для поиска компонент связности, O(M+N) для поиска и удаления мостов в каждой компоненте связности (оценка по худшему случаю, если компонента связности одна), O( CN / m) в худшем случае для формирования тхп-вершинных множеств и 0(1) для проверки их на соответствие искомому подграфу.

Рассмотрим алгоритм в общем виде:

1. На первом этапе с помощью алгоритма поиска в ширину [11] выделяются компоненты связности в графе. В рассмотрение не включаются компоненты связности, которые имеют меньше n вершин. Вершины, принадлежащие отброшенным компонентам связности, удаляются из графа вместе с инцидентными им ребрами.

2. Для каждой компоненты связности ищутся мосты с помощью известного алгоритма, основанного на поиске в глубину [11]. Из графа удаляются все найденные мосты и вершины, входящие в компоненты связности размерности меньше n, если таковые появились после удаления мостов.

3. В каждой из получившихся после предыдущих шагов компонент связности полным перебором ищутся неизолированные подграфы A'.

Исследования проводились для n = 3, 4 как для наиболее наглядных и простых случаев. Особенности топологии графа A' для случаев п = 3 и п = 4 учитывались при проверке условий из определения (1). Приведем пример адаптации определения графа A' для конкретных значений n = 3, m = 1:

A = A'(W, E), такой что

1. W=W0' и W1' и W2', W' n W'= 0, i, j = 0, 1, 2;

2. | Wi'i = m, i = 0, 1, 2;

3. (WW = 0, (W1,W2) = 0, (W2,W0) = 0, (W1,W)) = 1, (W2,W1) = 1, (W0,W2) = 1

или (W0,W1) = 1, (W1,W2) = 1,

(W2W0) = 1, (W1,W0) = 0, (W2,W1) = 0, (W0,W2) = 0.

Модель Эрдёша - Реньи (random graph)

В теории графов модель Эрдёша - Реньи [7] представляет собой модель для генерации случайного графа с постоянной вероятностью появления ребра между двумя вершинами незави-

симо от других ребер. Используется вариант модели Эрдёша - Реньи типа А^, р), согласно которой граф конструируется с помощью соединения пар вершин случайным образом. Каждое ребро включается в граф с вероятностью р независимо для каждой пары вершин. С ростом р от 0 до 1 граф становится все более насыщенным ребрами. При р = 0 имеем пустой граф (граф, не имеющий ребер), при р = 1 - полносвязный граф. При фиксированном количестве вершин в графе вероятность р - единственный входной параметр модели Эрдёша - Реньи.

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

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

Рис. 3. Математическое ожидание количества функциональных структур в зависимости от вероятности образования связи для случайного 100-вершинного графа по модели Эрдёша - Реньи, п = 3, т = 1 (красная линия) и п = 4, т = 1 (синяя линия)

_ s~im

Ц — CNCN - m

N!

N - (n-1)m

(m!)n(N - mn)!

/-tm где Cn —

N!

m!(N - m)!

. Согласно условию из

работы [5], для формирования функциональной структуры нужно п сильных и П2 — п = п2 — 2п слабых связей. Вероятность появления сильной связи р,, между двумя кластерами равна р_^ = = (1 - (1 — р)т)т, вероятность появления слабой

2

связи равна р№ = (1 — р)т . Вероятность того, что все нейроны внутри группы будут не связа-

а\т (т—1)

— р) 4 ,следова-

тельно, вероятность образования pf нужной топологии будет равна:

р/ = рП2—2пр: (I — р)тп( т—1) =

= ((1 - (1 - Р)m )т)”'-2п (1 - р)

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

< K >=НР/(п -1)! =

= N! ((1 - (1 - p)m )m )пП-2п (1 - p)n(2m2 -m)(n -1)!

= (m!)n(N - mn)! '

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

Результаты численных исследований, заключавшихся в проведении серии экспериментов и определении по полученной выборке числовых характеристик, где каждый эксперимент состоял в моделировании нейронной сети с помощью графа Эрдёша - Реньи и подсчете в графе числа функциональных структур, позволяют говорить о том, что наибольшее число неизолированных функциональных структур отыскивается в сетях, соответствующих графу модели Эрдёша - Реньи с вероятностью образования связи между двумя нейронами из интервала (0.25, 0.75)' Дисперсия числа структур в этом случае порядка 5-10 , что обусловливает достаточно большую выборку из 10000 экспериментов для каждого значения параметра.

Модель Ваттса - Строгаца (small-world graph)

)n(2m -m)

ции вершин (среднее значение по всем вершинам графа величины, характеризующей связность окрестности г-й вершины) и относительно небольшую среднюю длину пути (число вершин в простом полупути). Для генерации графа используются параметры N (число вершин), k (степень вершины) и в (специальный параметр: 0 < в < 1), причем 1 << 1^ << К << N .

На начальном этапе генерации графа строится регулярная кольцевая цепочка с N вершинами, w0...wN-l, причем дуга ^г, Wj) будет добавлено тогда и только тогда, когда | г — j | < k'(mod К), | k'І є (1, У.2). Далее для каждой вершины wг от Wo до wN-l дуга ^г, Wj) (где і < j) с заданной вероятностью в заменяется дугой ^г, wl), где I выбирается из всех возможных значений I ^г, с дублированием дуг (дуга (^г, wl') не добавляется, если I' = I, т.е. в полученном графе будут отсутствовать мультидуги). Процедура замены дуг называется рандомизацией графа.

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

На рис. 4 отражены зависимости математического ожидания числа функциональных структур от параметров в и а при фиксированном k. С увеличением числа этапов рандомизации уже при небольших значениях в ускоряется рост числа функциональных структур, однако после достижения граничного значения (математическое ожидание не более 1500 структур для 100-вершинного графа) их число начинает медленно убывать. Это объясняется тем, что при рандомизации количество дуг в графе может уменьшиться за счет запрета на добавление дублирующих дуг, при увеличении в и/или а число выпавших из графа дуг медленно увеличивается. Соответственно, в графе с меньшим числом дуг в среднем будет находиться меньшее число функциональных структур.

Дисперсия числа функциональных структур для этого типа графов независимо от количества этапов рандомизации составляет порядка 3-10 , что существенно меньше, чем значение аналогичной характеристики для случайного графа по модели Эрдёша - Реньи.

Модель Барабаси - Альберта

Модель Ваттса - Строгаца [9] представляет Модель Барабаси - Альберта [10] представ-

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

О 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1

Рис. 4. Математическое ожидание количества функциональных структур в зависимости от вероятности изменения связи для сети, соответствующей случайному 100-вершинному графу по модели Ваттса - Строгаца, для п = 3, т = 1 и различного числа итераций рандомизации исходного графа: 1 - красная линия, 3 - синяя линия, 5 -зеленая линия, 7 - фиолетовая линия

0 0.05 0.1 0.15 0 2 0.25 0 3 0.35 0 4 0 45 0.5 0.55 0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1

Рис. 5. Математическое ожидание количества функциональных структур в зависимости от процентного соотношения входящих и исходящих дуг для случайного 100-вершинного графа по модели Барабаси - Альберта для п = 3, т = 1 и различных значений вероятности р образования дуги в начальном графе: р = 0.05 (бирюзовая линия), р = 0.1 (фиолетовая линия), р = 0.3 (красная линия), р = 0.5 (синяя линия), р = 1 (зеленая линия)

степеней вершин, т.е. доля вершин степени k в графе Р(к) ~№, где -3 < у < -2, k ^ ж. Степенной закон распределения степеней вершин обусловливает наличие в графе небольшого числа вершин с наибольшими степенями (так называемых хабов) и большого числа вершин с небольшими степенями.

Алгоритм построения графа по модели Ба-рабаси - Альберта состоит в следующем: к начальному графу с N вершинами, где N >2 и степень каждой вершины больше или равна единице, на каждом шаге добавляется одна новая вершина. Новая добавленная вершина соединяется с N существующими вершинами с вероятностью, пропорциональной числу уже существующих у вершины ребер. Таким образом, N увеличивается на 1 на каждом шаге алгоритма. Опыты показывают, что при таком способе построения граф почти всегда получается связным и среднее значение длины пути в таком графе не превышает 1.89.

В графах, построенных по модели Бараба-си - Альберта, среднее число функциональных структур много меньше, чем в других типах графов при самой большой дисперсии (порядка 106). Из-за этой особенности выборка была увеличена до 10000 опытов. Начальное задание дуг для графа-основы практически не влияет на количество функциональных структур в графе (за исключением граничных случаев: пустого и полносвязного начальных графов). В зависимости от процентного соотношения входящих и исходящих дуг в этом графе наблюдается различная динамика для разных значений параметра (см. рис. 5). Резкий, скачкообразный характер зависимости обусловлен большой дисперсией. Наибольшее число функциональных структур наблюдается при соотношении входящих и исходящих дуг 1:10.

Сравнение моделей. Выводы

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

Численные исследования, проведенные с помощью разработанного программного комплек-

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

Модель Эрдёша - Реньи при достаточно большом разбросе значений параметра (р принадлежит интервалу (0.25, 0.75)) позволяет конструировать искусственные сети с большим числом одновременно существующих в них функциональных структур.

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

Модель Ваттса - Строгаца, являясь по сути переходным звеном между регулярной кольцевой цепочкой и случайным графом модели Эрдёша - Реньи, позволяет конструировать сети с б.'.льшим числом функциональных структур, нежели в сетях Эрдёша - Реньи, однако при этом параметры сети Ваттса - Строгаца не должны выходить из довольно узкого диапазона значений.

Для подтверждения численных результатов был проведен следующий эксперимент: была сконструирована сеть по модели Эрдёша - Ре-ньи из 30 нейронов с параметром сети р = 0.28. Было выяснено, что в этой сети имеется 392 функциональные структуры, состоящие из п = = 3 кластеров по т = 3 нейрона в каждом. Далее рассматривались только 3 из всех обнаруженных функциональных структур. Варьировалось значение величины внешнего тока, подаваемого на некоторые нейроны. При наличии асимметрии связей между нейронами происходит переход от спайковой активности нейронов к пачечной активности, которая характеризуется временным интервалом, в течение которого наблюдается активность.

На рис. 6 представлены отклики сети со случайными связями в форме различных активных последовательностей кластеров на различные наборы сигналов. В промежуток времени (0, 750) были активированы нейроны из 3 кластеров Wl = (1, 3, 18), W2 = (20, 25, 27), Wз = (9, 21, 29); в промежуток времени (750, 1500) были активированы нейроны из 3 других кластеров Wl = (5, 12, 18), W2 = (6, 9, 14), Wз = = (7, 19, 21); в промежуток времени (1500, 2250) были активированы нейроны из кластеров W1 =

30'

I

25

20

с

3 15

С

і

10

і

5-

I I

""1 Ill'll 11 ■ , її , m

Illll Hill ІНШІ 1 mu mu

Illll Illll Illll IIIII in

' 1 i i mil і і і і

375 750 1125 1500

Time, ms

а)

1875 2250

Рис. 6. Отклики сети на различные наборы сигналов. (а) Пространственно-временная диаграмма: последовательные переключения активности между различными группами как отклики на различные сигналы. Фиолетовые вертикальные линии соответствуют моментам изменения сигналов, подаваемых на нейроны. (б) Среднее значение мембранного потенциала V группы нейронов

= (2, 16, 26), W2 = (9, 15, 21), W3 = (14, 17, 29). В экспериментах внешние стимулы в различных случаях незначительно отличались. Тем не менее существование синхронной активности нейронов внутри одной группы хорошо видно из рис. 6б, на котором представлена зависимость среднего значения мембранного потен-

1 ы

циала нейронов =-,—в каждой группе

| уу1\ г=1

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

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

Работа выполнена при финансовой поддержке РФФИ („ранты 08-02-92004, 08-02-97049, 10-0200940) и ФЦП «Научные и научно-педагогические кадры инновационной России» на 2009—2013 гг. (ГК № П2308, № 02.740.11.5138, № П2018, № П15, № П942, № 02.740.11.5188).

Список литературы

1. Hahnloser R.H.R., Kozhevnikov A.A., Fee M.S. An ultra-sparse code underlies the generation of neural sequences in a songbird // Nature. 2002. V. 419. P. 65-70.

2. Jones L.M., Fontanini A., Sadacca B.F. et al. Natural stimuli envoke dynamic sequences of states in sensory cortical ensembles // PNAS. 2007. V. 104 (47). P. 18772-18777.

3. Hodgkin A.L., Huxley A.F. A quantitative description of membrane current and its application to conduction andexcitation in nerve // J. Physiol. 1952. V. 117. P. 500-544.

4. Nowotny T., Rabinovich M.I. Dynamical origin of independent spiking and bursting activity in neural microcircuits // Phys. Rev. Lett. 2007. V. 98(12): 128106.

5. Rabinovich M.I., Varona P., Selverston A.I., Ab-arbanel H.D. I. Dynamical principles in neuroscience // Reviews of Modern Physics. 2006. V. 78(4): 1213.

6. Komarov M.A., Osipov G.V., Suykens J.A.K., Rabinovich M.I. Numerical studies of slow rhythms

emergence in neural microcircuits: bifurcations and stability // Chaos. 2009. V.19(1):015107.

7. Komarov M.A., Osipov G.V., Suykens J.A.K. Sequentially activated groups in neural networks // EPL. 2009. V. 86. P. 6006.

8. Chung F., Lu L. Complex graphs & networks // CBMS Regional Conference Series in Mathematics. American Mathematical Society. 2006. V. 107.

9. Watts D.J., Strogatz S.H. Collective dynamics of ‘small-world’ networks // Nature. 1998. V. 393 (6684). P. 440-442.

10. Barabasi A-L., Albert R. Emergence of scaling in random networks // Science.1999. V. 286. P. 509.

11. Новиков Ф.А. Дискретная математика для программистов. СПб.: Питер, 2002.

SEQUENTIAL ACTIVITY STRUCTURES IN NEURAL NETWORKS WITH RANDOM CONNECTIONS

T.A. Levanova, M.A. Komarov, E. Yu. Kadina, G. V. Osipov

The effect of sequential bursting activity in neural networks with different connection topology is studied. We consider the networks with random connectivity between neurons, small-world networks and scale-free networks with power law degree distribution. The dependence is established of the number of neural structures (in which the sequential bursting activity is possible in response to external stimuli) on the parameters of each type of network.

Keywords: neural network, functional structure, Erdôs - Rényi model, Watts - Strogatz model, Barabasi - Albert model.