CC BY
64
ВЕСТНИК ТГГПУ. 2008. №4(15)
УДК 6128:6815
МОДЕЛЬ РОСТА НЕЙРОННОЙ СЕТИ
© Ф.А.Галимянов, Ф.М.Гафаров, Н.Р.Хуснутдинов
В статье представлена модель, описывающая динамическое образование нейронных сетей. Целью работы является моделирование процесса возникновения нейронной сети. Образование нейронной сети в нашей модели осуществляется путем соединения кончика аксона с телом нейрона. Все нейроны по своим формам, свойствам, и поведению абсолютно одинаковы. Положения нейронов фиксировано, и рост аксонов единственный способ образование сети. Предложенная модель применяется как для трехмерного, так и двухмерного случаев.
Ключевые слова: математическая модель, нейронная сеть, аксон, численный анализ
Математическая модель
Координаты центра нейрона располагаются в вершинах прямоугольника или куба (в зависимости размерности). Положение кончика аксона каждого нейрона описывается радиус-вектором Г (і), и в начальный момент равно координате
центра нейрона. Нейроны обладают свойством называемым активностью, которая подчиняется дифференциальному уравнению:
&З, (і) ( * ^
&
1=1
Здесь т, (і) коэффициент затухания ак-
тивности и сама активность і -го нейрона в момент времени і; /(х~) = хв(х) , где в(х) - ступенчатая функция, м ,і - веса, которые могут
принимать три значения "1", "0" или "-1". Веса бывают положительными, отрицательными или же "нулевыми" и определяют тип связи между нейронами. Типы связей (веса) зависят от активности і -го нейрона, с которым соединился аксон 1 -го нейрона, если активность і -го нейрона больше чем пороговая Зпор (величина безразмерная), то связь отрицательная и наоборот [1], [2]:
Г-1, при З{ (і)> 0,51 м. . = ^
і,1 | 1, при З. (і)< 0,51
Вначале все связи (веса) считаются "нулевыми", то есть связи считаются не появившимися и остаются нулевыми до образования связи. Связь появляется, если конец аксона 1 -го нейрона приблизится к центру і -го нейрона на расстояние, которое меньшее чем половина радиуса аксона. ЗГ1 (і) - сигнал входящий в систему извне.
В зависимости от активности З. (і) нейрон
испускает в межклеточное пространство вещество, которое влияет на аксоны нейронов находя-
щихся в системе, а также на собственный аксон. Аксон не растет, в случае если активность его собственного нейрона больше порогового значения Зпор . Связь не может появится у нейрона с
его же аксоном то есть если і = 1, то всегда
^1 = °.
Воспользуемся стандартным уравнением теплопроводности для описания распространения вещества в среде. Решая уравнение:
- В2 Ас. - ке. = 3 (г,і),
ді ' ' Л1’ г
можно найти концентрацию вещества выделившегося из і -го нейрона в момент времени і .
Здесь к - коэффициенты деградации. Решение этого уравнения описывает концентрацию вещества в точке Г в момент времени і, и оно
хорошо известно:
е. (г. - г.,і) = Г О, (г - г,і) е. (г )&?.+
1 \ 1 і / ^ & \ J 1 / 1 V і / і
і
+| &ік І О& ( - гі,і - ік) (^л Н,
0
где (г,і) - функция Грина:
О& (г,і):
1
2
ехр
Г 2 Л
- кі --
4іБ
(піБ2 )
Здесь Б2 - коэффициент теплопроводности, а плотность тока 3І (г,і) ,
3 (г,і) = а5{() ( - г{) (і),
выражается через активность З. (і) .
Решая эти уравнения совместно, получаем уравнение для полной концентрации:
і
ег ( - Г,,і) = аІ&ікОа ( - гі,і - ік) (ік) .
0
Уравнение движения аксона 1 -го нейрона, описываемого радиус - вектором §і , записывается в следующем виде:
^ = АР(X,)£ Vс. (|, - г„.,) .
Ш .=1
Здесь А - коэффициент, описывающий чувствительность и подвижность аксона,
Р (Я) = 0( Бпор - Я) , где Я пор - это пороговое значение активности.
Численный анализ
Далее мы численно решаем уравнения движения аксонов и вычисляем активности. Для решения дифференциального уравнения используем метод Эйлера, программной средой является язык С++, погрешность вычисления допускается при шаге итераций А1 = 20 сек. и менее.
Параметры, используемые при решении уравнений, имеют следующие значения: ё = 3*10~3 см - диаметр сомы нейрона; г = 0,5 мм - расстояния между нейронами;
Япор = 0,51 - пороговое значение активности;
к = 10-3 1/сек - коэффициенты деградации;
О2 = 6*10-7 см2 сек - коэффициент теплопроводности [3];
6 2 I
А=4*10- см /сек*нМ - коэффициент, описывающий чувствительность и подвижность аксона.
Далее задаем последовательность и время действия внешних сигналов. Девятый нейрон получает внешний сигнал с нулевой секунды и до 800 секунды; первый нейрон с 800 до 1200; восьмой с 1200 до 3500; на первый нейрон вновь начинает поступать внешний сигнал с 3500 секунды и до конца эксперемента, т.е. до 40000 секунды. См. Рис.1.
40000 0
Рис.1. Рост аксонов на плоскости. Кружками обозначены месторасположения нейронов. В нижнем правом углу показано время в секундах. Рядом с нейронами указаны их активности в данный момент времени
В трехмерного случае расчет производим аналогичным образом, меняется лишь число нейронов и время подачи внешнего сигнала на них. На первый нейрон сигнал подается с началом отсчета и действует до 2000 секунды; на восьмой с 2000 по 4000; на третий с 4000 по 6000; на семнадцатый с 6000 по 9000; на одиннадцатый с 9000 до конца эксперимента, т.е. 40000 секунды. См. Рис.2.
о п ■» lb / 0.0
r U і] 0.0
/ 1\ l'w 0.0
----------------- • /л / ^ \ V
Ф Хо ф у о.о
і 0.0 ф 0.0 ф 0.0
40000.0
Рис.2. Рост аксонов в пространстве. Сферами обозначены месторасположения нейронов. В нижнем правом углу показано время в секундах. Рядом с нейронами указаны их активности в данный момент времени
Суммируем результаты работы. В статье предложена реальная модель роста нейронной сети, основанная на уравнении диффузии вещества, управляющего ростом аксона. Используя экспериментальные значения параметров нейронной сети, показано, что картина роста сети соответствует наблюдаемым данным. Наблюдаемая топология растущей нейронной сети управляется активностями нейронов.
1. Dickson B.J. Molecular Mechanisms of Axon Guidance // Science. - 2002. - 298. - P.1959-1964.
2. Gafarov F. Self-wiring in neural nets of point-like cortical neurons fails to reproduce cytoarchitectural differences // J. Int. Neurosc. 2006. - 5. - P.159-169.
3. Gomez T.M. and Zheng J.Q. The molecular basis for calcium-dependentaxon pathfinding // Nature Neurosci. - 2006. - 7. - P. 115-125.
ФИЗИКА
THE MODEL OF DYNAMIC GROWTH OF NEURON NETWORK
F.A.Galimyanov, F.M.Gafarov, N.R.Khusnutdinov
Here we present the model of dynamic growth of the neuron network. The network is created by connection of growth cone with neuron's soma. All neurons are equivalent to each other. The neurons have fixed positions, and the growth of axons is the only way of creation of the network.
We use the model for two and three dimensional cases with neuron's soma as a sphere.
Key words: mathematical model, neuron network, axon, numerical analysis
Галимянов Фанис Анисович - аспирант кафедры информационных технологий в образовании
Гафаров Фаиль Мубаракович - кандидат физико-математических наук, доцент кафедры теоретической физики Татарского государственного гуманитарно-педагогического университета
Хуснутдинов Наиль Рустамович - доктор физико-математических наук, заведующий кафедрой информационных технологий в образовании Татарского государственного гуманитарно-педагогического университета
E-mail: ff@tggpu.ru
