Экспериментальное и теоретическое исследование автогенераторных моделей нейронных систем тема автореферата и диссертации по физике, 01.04.03 ВАК РФ

Щапин, Дмитрий Сергеевич АВТОР
кандидата физико-математических наук УЧЕНАЯ СТЕПЕНЬ
Нижний Новгород МЕСТО ЗАЩИТЫ
2011 ГОД ЗАЩИТЫ
   
01.04.03 КОД ВАК РФ
Диссертация по физике на тему «Экспериментальное и теоретическое исследование автогенераторных моделей нейронных систем»
 
Автореферат диссертации на тему "Экспериментальное и теоретическое исследование автогенераторных моделей нейронных систем"

На правах рукописи

ЩАПИН Дмитрий Сергеевич

ЭКСПЕРИМЕНТАЛЬНОЕ И ТЕОРЕТИЧЕСКОЕ ИССЛЕДОВАНИЕ АВТОГЕНЕРАТОРНЫХ МОДЕЛЕЙ НЕЙРОННЫХ СИСТЕМ

01.04.03 - радиофизика

4852688

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата физико-математических наук

I 8 СЕН 2011

Нижний Новгород - 2011

4852688

Работа выполнена в Учреждении Российской академии наук Институт прикладной физики РАН (г. Нижний Новгород)

Научный руководитель: д. ф.-м. н., профессор

Некоркин Владимир Исаакович

Официальные оппоненты: д. ф.-м. н., профессор

Короновский Алексей Александрович

д. ф.-м. н., профессор Осипов Григорий Владимирович

Ведущая организация: Учреждение Российской академии наук

Институт радиотехники и электроники им. В.А. Котельникова РАН, г. Москва

Защита состоится « 28 » сентября 2011 г. в 15. часов на заседании диссертационного совета Д 212.166.07 при Нижегородском государственном университете имени Н И. Лобачевского (603950, Нижний Новгород, ГСП-20, проспект Гагарина, 23, корп. 1, ауд. 420).

С диссертацией можно ознакомиться в библиотеке Нижегородского государственного университета имени Н.И. Лобачевского.

-АЗ

Автореферат разослан «_!___» августа 2011 г.

Ученый секретарь диссертационного совета к.ф.-м.н., доцент ' / В.В. Черепенников

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

В последнее время наблюдается увеличивающийся интерес к изучению, методами нелинейной динамики, активных ансамблей, моделирующих поведение большого числа взаимодействующих нервных клеток — нейронов. С одной стороны, это объясняется современным прогрессом в развитии методов регистрации нейронной активности (мультиэлектродная регистрация, оптический нейроимиджинг и др.), позволяющих получить большое число новых экспериментальных данных о режимах активности достаточно больших нейронных сетей. Систематизация, обобщение и понимание закономерностей возникновения, распространения и исчезновения активности таких сетей стимулировали развитие так называемого нейроморфного моделирования. Модели этого направления учитывают морфологию индивидуальных нейронов и межнейронных связей, имеют архитектуру отражающую строение реальных нейронных систем. В нейроморфных моделях элементы-нейроны наделяются своей внутренней и, в частности, колебательной динамикой, которая в той или иной степени близка к динамике "живых" нейронов. С другой стороны, современные электронные технологии, в частности, так называемые программируемые кристаллы (программируемые логические интегральные схемы - ПЛИС, программируемые блоки аналоговых элементов - field-programmable analog array, и др.), допускают создание компактных электронных схем из большого числа активных единиц, которые, могут воспроизводить базовые свойства нейронов. В связи с этим, построение новых нетрадиционных высокоэффективных систем обработки и хранения информации, координации и управления движением, основанных на нейродинамических принципах, стало актуальной и практически реализуемой задачей.

В настоящее время нелинейно-динамический подход для исследования свойств нейронных сетей активно развивается как в России (Борисюк Г. Н.,

Борисюк Р. М., Безручко Б.П., Иваницкий Г.Р., Казанович Я. Б., Казанцев В.Б., Некоркин В. И., Осипов Г.В., Потапов А.Б., Смирнов Д. А., Яхно В. Г.), так и за рубежом (Abarbanel H.D.I., Afraimovich V. S.,'Atoara К., Arena P., Ayers J., Bazhenov M., Bilbaut J.M., Courbage M., Ennentrout G.B., Fortuna L., Ijspreert A. J., Izikevich E.M., Kurths J., Llinas R., Makarenko V. I., Rabinovich M. I., Rinzel J., Rulkov N. F., Sanjuan M. A . F., Tanaka G., Terman D. H., Varona P., Velarde M. G. и др.)

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

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

Научная новизна работы

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

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

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

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

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

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

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

Основные положения, выносимые на защиту

1. Разработанные автогенераторы обладают колебательными, возбудимыми и мультистабильными режимами.

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

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

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

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

Практическая и теоретическая значимость результатов

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

Апробация результатов работы

Основные результаты докладывались на следующих российских и международных конференциях: седьмой научной конференции по радиофизике (Нижний Новгород 2003), на конференции молодых ученых "Нелинейные волновые процессы" (Научная школа "Нелинейные волны" 2004, 2006, 2008, 2010 Нижний Новгород), на нижегородской сессии молодых ученых (естественнонаучные дисциплины 2004), на всероссийской конференции "Хаотические автоколебания и образование структур" (Саратов 2004, 2007), International Symposium "Topical Problems of Nonlinear Wave Physics" (Nizhny Novgorod, Russia 2003; St.-Petersburg - Nizhny Novgorod, Russia 2005; Nizhny Novgorod, Russia 2008), International Workshop on Nonlinear Dynamics of Electronic Systems, (Dijon, France 2006; Tokushima, Japan 2007), Chaos and dynamics in biological networks (Cargese, France 2010).

Личный вклад автора

Все эксперименты выполнены лично автором. В совместных работах, теоретический анализ и интерпретация полученных результатов были выполнены с научным руководителем. Результаты опубликованы в работах [1-22].

Структура и объем работы

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

КРАТКОЕ СОДЕРЖАНИЕ РАБОТЫ

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

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

Генератор ФитцХью-Нагумо. Схема генератора представлена на рис. 1а, где напряжение U играет роль мембранного потенциала, С - мембранной емкости, R=Rj+R2 - сопротивления утечки, Е - равновесного потенциала, контролирующего уровень деполяризации мембраны, F(U) - нелинейной проводимости мембраны, L=RtR2CL - индуктивности отвечающей за скорость изменения ионных токов, Евн - внешнего стимула. Нелинейная проводимость имеет форму кубической параболы (рис. 16), которая аппроксимируется функцией F(U)= aUfU-uo) (U+ и0), где коэффициент а имеет размерность (ОмВ2)"1. Динамика токов и напряжений данного генератора описывается системой

ах

(1)

Система (1) была исследована экспериментально и теоретически. При отсутствии внешнего стимула (Е„„(/)=0), построены теоретически и экспериментально бифуркационные диаграммы, разделяющие плоскость параметров (Я,Е) на области, соответствующие различным режимам генерации (рис. 2).

тС

Рис. 1. а - схема генератора ФитцХью-Нагумо; 6 - вольт-амперная характеристика нелинейного элемента, квадратами отмечены экспериментальные данные, линии соответствует кубическая функция Р(Ь')

«л

А ИВ)

\ 1 0,4 ■

с М-

—.н,

1 т 0-

-—"г

4 V*

У и -ол-

я

И 5000 10000

15000 Л(Ом)

15000 Д(Ои)

Рис. 2. а - бифуркационная диаграмма системы (1); б - экспериментальная бифуркационная диаграмма

На рис. 2а кривые лил отвечают седло-узловым бифуркациям состояний равновесия и разделяют плоскость на две области, в которых система (1) име-

7

ет либо одно О, либо три состояния равновесия Оь 02 и 03, где 02 является седлом. Кривые Ах, А3 соответствуют бифуркациям Андронова-Хопфа состояний равновесия О, О/ и 0}, в результате которых происходит рождение устойчивых или неустойчивых предельных циклов. Кривые Яь Н2 отвечают нелокальным бифуркациям "малых" петель сепаратрис седла, а кривые Н21, Н\г (не указаны на рисунке) "больших" (охватывающих все три состояния равновесия) петель сепаратрис седла и соответствуют рождению неустойчивых предельных циклов при их разрушении. Кривая С отвечает бифуркации двукратных предельных циклов. Бифуркационные кривые разделяют плоскость (Я,Е) на 14 областей, отвечающих существованию колебательного, возбудимого, и различных мультистабильных режимов, переключение между которыми может осуществляться за счет выбора соответствующих параметров или за счет внешнего кратковременного воздействия £«„(?)#О (в случае мультистабильных режимов). Бифуркационная диаграмма полученная экспериментально, представлена на рисунке 26. Заметим, что она достаточно хорошо согласуется с теоретической. Однако, некоторые из бифуркационных кривых (например, бифуркация большой петли сепаратрис седла), находящихся близко друг к другу на плоскости (/?,£) на (рис. 2а) экспериментально различить не удалось.

Генератор, моделирующий динамику нейронов с подпороговой активностью. Генератор состоит из двух блоков, рис. За. Первый блок представляет генератор в колебательном режиме, а второй блок отвечает за генерацию импульсов возбуждения. Блоки соединены через повторитель, который обеспечивает однонаправленную связь между ними. Характерными режимами такого генератора являются генерация одиночных импульсов возбуждения (спай-ков) на пике подпороговых колебаний, а также генерация серий импульсов -беретов, на верхнем полупериоде подпороговых колебаний.

Г П 1 1—■■ 1

II И 1 1

1.4 1,6 С

а 6

Рис. 3. а - генератор, моделирующий нейрон с подпороговой активностью; б — режим перемежающейся синхронизации

Экспериментально была исследована динамика двух таких генераторов (далее обозначим их ¿V/, ЛУ, взаимодействующих через подавляющую обрат-

ную связь. Частоты блоков, отвечающих за генерацию подпороговых колебаний, различались на несколько процентов, а сами генераторы находились в режиме генерации одиночных спайков на фоне подпороговых колебаний. Генераторы связаны сопротивлением Я(ис), значение которого зависит от управляющего напряжения ис. Петля обратной связи изменяет сопротивление Щис) между генераторами следующим образом. Если N1 или М2 генерирует спайк, то сопротивление Я(ис) возрастает, т. е. связь между элементами становится малой на интервал времени г, где г является параметром, величину которого можно изменять. Показано, что в зависимости от значения параметра г в системе могут наблюдаться различные режимы синхронизации, начиная от полной до так называемой перемежающейся синхронизации. В режиме перемежающейся синхронизации существуют временные интервалы, в течение которых спайки синхронны, а вне этих интервалов синхронизации спайков нет. При этом интервалы синхронизации появляются во времени нерегулярно. Отметим, что такие режимы синхронизации наблюдаются и играют важную роль в формировании кластеров активности оливо-мозжечковой системы.

Ггнератор, моделирующий динамику нейронов с последеполяризацией. Блок-схема такого генератора приведена на рисунке 4а. Характерная особенность нейронов с последеполяризацией состоит в том, что после генерации потенциала действия возбудимость нейронной мембраны временно повышается (рис. 46), и для повторного возбуждения достаточно сравнительно небольшого стимула. Для реализации свойства последеполяризации был построен генератор с дополнительным блоком обратной связи. Если в генераторе возникает импульс возбуждения, то он активирует дополнительный блок генерации экспоненциального импульса, которой затем поступает на вход генератора, тем самым увеличивая его возбудимость. и

Рис. 4. а - генератор, моделирующий нейрон со свойством последеполяризации; б -последовательное воздействие трех внешних стимулов различной амплитуды

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

Во второй главе проведено экспериментальное исследование двух типов ансамблей взаимодействующих генераторов, представленных на рис. 5а,в. Изучена пространственно-временная динамика ансамбля из шестнадцати генераторов ФитцХью-Нагумо последовательно соединенных сопротивлениями связи Лс(рис. 5а) Показано, что в таком ансамбле, в зависимости от значений параметров, возможно распространение импульсов возбуждения и волновых фронтов переключения двух типов "кинков" и "антикинков". Экспериментально изучено взаимодействие волновых фронтов и импульсов возбуждения. Установлено, что волновые фронты переключения и импульсы возбуждения могут аннигилировать при взаимодействии друг с другом, либо проявлять частицеподобные свойства, т. е. отражаться при взаимодействии друг с другом и границами цепочки (рис. 56). Эффекты отражения волновых паттернов наблюдаются, когда параметры каждого генератора находятся вблизи бифуркационной границы С (см. рис. 2). Изучена пространственно-временная динамика импульсов возбуждения в ансамбле связанных автогенераторов,

Рис. 5. а - цепочка резистивно связанных генераторов ФитцХью-Нагумо; б - пространственно-временные диаграммы, в случае отражения фронтов переключения от границ и прохождение импульсов возбуждения друг через друга и отражение от границ ансамбля; в - цепочка генераторов с двумя нелинейными проводимостями; г - пространственно-временные структуры хаотической активности импульсов возбуждения

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

^0<"<>-

в(Ц) =

—и,и< о 3

_1_

(2)

и,и* о

У-1.2.....16, (/«,(0^ £/,(/), </„(0-£/„(/)

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

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

й = у,

■ у = су + Ф{и)-1+ (3)

а = -Е2-<%'(и)у,

Здесь Ф(м)=-Дм)+я(м), а функции .Дм) и ^и) соответствуют функциям Р(11) и С!(и), точкой обозначено дифференцирование по бегущей координате. Бифуркационный анализ системы (3) показал, что она имеет три состояния равновесия типа седло-фокус. Два из этих состояний равновесия и 03, имеют неустойчивое одномерное и устойчивое двумерное многообразия, а состояние равновесия 02, наоборот, имеет устойчивое одномерное и неустойчивое двумерное многообразия. Методом двумерных систем сравнения, локализующих в фазовом пространстве устойчивые и неустойчивые многообразия седловых состояний равновесия О/ и 03, было доказано существование гетероклинических траекторий (рис. 6). Этим траекториям соответствуют распространяющиеся с постоянной скоростью волновые фронты (кинки и антикинки) в цепочке генераторов. Показано, что кроме волновых фронтов в цепочке возможно распространение импульсов возбуждения, которым в системе (3) отвечают гомоклинические орбиты. На плоскости (1,с) установлены значения параметров, отвечающие бифуркации гомоклинических траекторий (множества Гь Г3), а также множества Н* и Я", отвечающие гетероклиническим орбитам, которые определяют зависимость скорости с, для волновых импульсов и фронтов переключения, от параметра деполяризации / (см. рис. 6).

0.05 0.06

Рис. 6. Бифуркационные множества Н*, Н~ и Г\, Г3, отвечающие существованию соответственно гетероклинических и гомоклинических орбит в системе (3)

Показано, что бифуркационные множества Г\, Гъ, ИГ, Н, имеют несколько точек пресечения, которые отвечают существованию в фазовом пространстве системы (3) сложных гетероклинических контуров, образованных сепа-ратрисными многообразиями седло-фокусов О/ и 03 с положительной седло-вой величиной. В окрестности таких контуров, в частности, существует нетривиальное гиперболическое множество, содержащее бесконечное множество седловых периодических траекторий. Таким образом, для значений параметров, взятых в окрестности данных структур, система для бегущих волн демонстрирует чрезвычайно сложную динамику, что показывает сложное пространственно-временное поведение исходного ансамбля генераторов. В частности, это подтверждается экспериментально наблюдаемыми режимами хаотической импульсной активности (рис. 5г).

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

1,00,50,01

1,0

0,0 -04 -1,0

1

Рис. 7. а - эффект фазовой автопереустановки; б - шагающий робот

от начальной фазы), а определяется лишь интенсивностью внешнего импульса, рис. 7а. Изучена зависимость точности переустанавливаемой фазы от характеристик импульса. Показана возможность использования данного эффекта для импульсно управляемой синхронизации неидентичных автогенераторов. Реализован ансамбль с центральным генератором, относительно которого устанавливаются фазы всех остальных генераторов. В таком ансамбле возможно формирование фазовых паттернов, любой наперед заданной конфигурации. На основе такого ансамбля генераторов, обладающих эффектом фазовой автопереустановки, была построена система управления шагающим роботом, рис. 76. Экспериментально были установлены базовые шаблоны задающие движение робота (по плоской поверхности) "вперед", "влево" и "вправо". Система управления, используя базовые шаблоны движения и входной (сенсорный) сигнал с видеокамеры, формирует новые моторные шаблоны движения, позволяющие роботу следовать за целью.

В заключении приведены основные результаты и выводы диссертационной работы.

ОСНОВНЫЕ РЕЗУЛЬТАТЫ И ВЫВОДЫ

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

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

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

4. Реализован и исследован генератор, моделирующий динамику нейронов, обладающих свойством последеполяризации (понижение порога возбуждения после генерации спайка).

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

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

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

8. Экспериментально показано, что автогенератор ФитцХью-Нагумо обладает эффектом фазовой автопереустановки, который позволяет управлять фазой колебаний. При подаче на генератор прямоугольного импульса фаза колебаний переустанавливается к новому значению, которое не зависит от момента воздействия внешнего импульса (или от начальной фазы), а определяется лишь интенсивностью импульса, т. е. его амплитудой и длительностью.

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

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

СПИСОК РАБОТ ПО ТЕМЕ ДИССЕРТАЦИИ

1. Престунов В. А., Щапин Д. С., Казанцев В. Б., Некоркин В. И. Эффекты синхронизации и десинхронизации в системе двух динамических связанных автогенераторов. //Труды (седьмой) науч. конф. по радиофизике, посвященной 90-летию со дня рождения B.C. Троицкого / ред. А. В. Якимов. Н. Новгород: ТАЛАМ. 2003. С. 120.

2. Prestunov V. V., Shapin D. S., Kazantsev V. В., and Nekorkin V. I. The dynamics of two oscillators with pulse-controlled coupling. Abstr. Int Symp. "Topical Problems of Nonlinear Wave Physics". Nizhny Novgorod. Russia. 2003. P. 74-75.

3. Щапин Д. С., Престунов В. А., Казанцев В. Б., Некоркин В. И. Модель нейрона с подпороговыми колебаниями // Тезисы докл. Конф. молодых ученых "Нелинейные волновые процессы (XII науч. школа "Нелинейные волны-2004"). Н. Новгород, 2004. С. 92-93.

4. Престунов В. А., Щапин Д. С. Динамика модели нейрона с подпороговой активностью // Тезисы докладов IX Нижегородской сессии молодых ученых (естественнонаучные дисциплины), Н. Новгород, 2004. С. 120.

5. Дмитричев А. С., Щапин Д. С., Казанцев В. Б., Некоркин В. И. Сложная волновая динамика в ансамбле взаимосвязанных элементов ФитцХью-Нагумо и сепарат-рисные контура // Матер. VII междунар.школы "Хаотические автоколебания и образование структур - 2004". Саратов: Изд-во ГосУНЦ "Колледж", 2004. С. 119-120.

6. Престунов В. А., Щапин Д. С., Казанцев В. Б., Некоркин В.И. Управление фазой автоколебаний релаксационного генератора при помощи внешнего стимула // Матер. VII междунар. школы "Хаотические автоколебания и образование структур -2004". Саратов: Изд-во ГосУНЦ "Колледж", 2004. С.120-121.

7. Shapin D. S., Kazantsev V. В., and Nekorkin V. I. Synchronization and time binding in two neuronal oscillators with inhibitory feedback : abstr. // Int. Symp. "Topical Problems of Nonlinear Wave Physics", Nizhny Novgorod, 2005. P. 95-96.

8. Щапин Д. С. Фазовая автопереустановка и синхронизация релаксационных автоколебательных систем // Тезисы докл. X Нижегородской сессии молодых ученых (естественнонаучные дисциплины). Н. Новгород, 2005. С. 160-161.

9. Некоркин В. И., Дмитричев А. С., Щапин Д. С., Казанцев В. Б. Динамика модели нейрона со сложно-пороговым возбуждением // Математическое моделирование. 2005. Т. 17, №6. С. 75-91.

10. Щапин Д. С. Импульсное блокирование связи в сети нейроноподобных элементов в физическом эксперименте. Тезисы докладов конференции молодых ученых "Нелинейные волновые процессы". Н. Новгород. 2006. С. 169.

11. Kazantsev V. В., Dmitrichev A. S., Shapin D. S. and Nekorkin V. I. Multi-threshold excitability in-a nonlinear network of neuron-like units H The 14-th International Workshop on Nonlinear Danamics of Electronic Systems. Dijon. France. 2006. P. 69-72.

12. Kazantsev V. В., Shapin D. S. and Nekorkin V. I. Phase clusters in oscillatory neuronal network with inhibitory feedback // The 14-th International Workshop on Nonlinear Danamics of Electronic Systems. Dijon. France. 2006. P. 73-77.

13. Shapin D.. S. Phase control in nonlinear electronic circuit modeling oscillatory neurons // The 15th International Workshop on Nonlinear Dynamics of Electronic Systems. Tokushima. Japan. 2007. P. 194-196.

14. Щапин Д. С. Управление фазой колебаний в нелинейной электронной схеме модели нейрона // 8-я международная школа "ХАОС-2007": Тез. докл.- Саратов: Изд. СГУ, 2007. С. 100.

15. Некоркин В. И., Щапин Д. С., Дмитричев А. С. Сложная волновая динамика ансамбля нейроноподобных элементов со сложно-пороговым возбуждением // Изв. ВУЗов. Прикладная нелинейная динамика. 2007. Т. 15, №1. С. 3-22.

16. Щапин Д. С. Экспериментальное исследование динамики модели нейрона с подпороговыми колебаниями // XIV научная школа "Нелинейные волны - 2008". Фундаментальные и прикладные задачи нелинейной физики. Конференция молодых ученых : тезисы докладов. Н. Новгород. 2008. С. 177.

17. Shapin D. S. The dynamics of two coupled neuron-like electronic elements with inhibitory feedback // Topical problems of nonlinear wave physics (NWP - 2008). Proc. of the 16th Int. Workshop on Nonlinear Dynamics of Electronic Systems (NDES - 2008): Niz-hny Novgorod, 2008. P. 87-88.

18. Nekorkin V. I., Shapin D. S., Dmitrichev A. S., Kazantsev V. В., Binczak S., Bil-bault J. M. Heteroclinic contours and self-replicated solitary waves in a reaction-diffusion lattice with complex threshold excitation // Physica D. 2008. V. 237. P.2463-2475.

19. Щапин Д. С. Динамика двух нейроноподобных элементов с подавляющей обратной связью // Радиотехника и электроника. 2009. Т. 54, №2. С. 185-195.

20. Щапин Д. С., Дмитричев А. С., Казанцев В. Б., Некоркин В. И. Генерация моторных паттернов и их перестройка по сенсорному сигналу в нейродинамической системе управления шагающего робота // XV Научная школа "Нелинейные волны -2010". Фундаментальные и прикладные задачи нелинейной физики. Конференция молодых ученых : тезисы докладов. Н. Новгород. 2010. С. 136.

21. Shchapin D.S., Dmitrichev A. S., Kazantsev V. В., Nekorkin V. I. Phase cluster based motor pattern generation in oscillatory network controller of hexapod robot walking gates // Chaos and dynamics in biological networks - 2010. Carges - France, 2010. P. 12.

22. Клиньшов В. В., Щапин Д. С., Некоркин В. И. Моделирование нейродинамической системы рабочей памяти // Радиотехника и электроника. 2010. Т. 55, №7. С. 812-817.

ЩАПИН Дмитрий Сергеевич

ЭКСПЕРИМЕНТАЛЬНОЕ И ТЕОРЕТИЧЕСКОЕ ИССЛЕДОВАНИЕ АВТОГЕНЕРАТОРНЫХ МОДЕЛЕЙ НЕЙРОННЫХ СИСТЕМ

Автореферат

Подписано к печати 6.07.2011 г. Формат 60x90 '/|б. Бумага офсетная. Усл. печ. л. 1,0. Тираж 100 экз. Заказ № 66(2011).

Отпечатано в типографии Института прикладной физики РАН, 603950 г. Нижний Новгород, ул. Ульянова, 46

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

Введение

1 Автогенераторы, моделирующие динамические свойства нервных клеток

1.1 Аналоговая модель автогенератора

ФитцХью-Нагумо.

1.2 Бифуркационный анализ модели.

1.2.1 Поглощающая область

1.2.2 Состояния равновесия системы (1.4).

1.2.3 Симметрия системы (1.4).

1.2.4 Функции Ляпунова. Существование в системе (1.4) гомоклинических траекторий.

1.2.5 Гомоклиническая орбита Г1.

1.2.6 Гомоклиническая орбита Г2.

1.2.7 "Большие" гомоклинические орбиты Гх^ и Г2д

1.2.8 Бифуркационные диаграммы и фазовые портреты

1.2.9 Экспериментальная бифуркационная диаграмма

1.3 Динамические режимы.

1.4 Генератор моделирующий динамику нейронов с подпороговой активностью.

1.5 Синхронизация двух автогенераторов с подпороговой активностью

1.6 Модель подавляющей обратной связи.

1.7 Генератор моделирующий динамику нейронов с последепо-ляризацией.

1.8 Выводы.

2 Ансамбли взаимодействующих генераторов

2.1 Распространение электрической активности в ансамбле резистивно связанных генераторов

ФитцХью-Нагумо.

2.2 Структуры хаотической спайковой активности в ансамбле генераторов с двумя нелинейными проводимостями.

2.3 Профили бегущих волн в ансамбле генераторов.

2.3.1 Система для бегущих волн.

2.3.2 Бифуркации состояний равновесия в системе для бегущих волн.

2.3.3 Нелокальное поведение многообразий.

2.3.4 Релаксационная динамика системы для бегущих волн

2.3.5 Гетероклинические траектории.

2.3.6 Построение систем сравнения.

2.3.7 Взаимное расположение многообразий.

2.3.8 Гомоклипические орбиты.

2.4 Выводы.

3 Фазово-управляемые колебания в ансамбле генераторов

3.1 Фазовая автопереустановка в автогенераторе ФитцХью-Нагумо

3.1.1 Описание процесса фазовой автопереустановки

3.2 Синхронизация неидентичных генераторов в ансамбле с центральным элементом.

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

3.4 Выводы.

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

Исследование систем, состоящих из большого числа взаимодействующих активных элементов, обладающих собственной колебательной динамикой, является актуальной проблемой современной радиофизики. Примерами таких систем в радиофизике являются массивы джозефсоновских контактов [1,2], сети связанных лазеров [3,4], ансамбли генераторов и сети фазовой синхронизации [5-7], фазированные антенные решетки [8] и т. д. В пространстве элементы могут быть расположены как произвольно, так и упорядоченно, например, находится в узлах пространственной решетки. Связь между элементами варьируется от простой линейной до сложной нелинейной, осуществляя как локальное, так и нелокальное взаимодействие между элементами. Фактически, ансамбли активных элементов представляют собой некоторые активные среды, изучение явлений и процессов в которых является фундаментальной задачей современной радиофизики. Для таких сред особенно интересны и важны режимы коллективной активности и распространения нелинейных волн, осуществляющие динамическую обработку информационных потоков.

В последнее время наблюдается увеличивающийся интерес к изучению, методами нелинейной динамики, активных ансамблей [9-15], моделирующих поведение большого числа взаимодействующих нервных клеток -нейронов. Это объясняется современным прогрессом в развитии методов регистрации нейронной активности (мультиэлектродная регистрация, оптический нейроимиджинг и др.), позволяющих получить большое число новых экспериментальных данных о режимах активности достаточно больших нейронных сетей. Исследования различных областей головного мозга показывают наличие колебательно-волной активности в различных нейронных системах в диапазоне от долей до сотен герц [16-25]. Например, для оливо-мозжечковой системы, отвечающей за контроль и координацию движений, значения регистрируемых частот составляет 812Гц [20,21,26-28], колебательные процессы в гиппокампе (с частотами 4-12Гц), связывают с процессами отвечающими за обработку сенсорной информации и за краткосрочную память [16,17,29-33]. Процессы в этих системах свидетельствуют о том, что нейронные системы способны формировать пространственно-временные кластеры синхронной колебательной активности за счет собственной динамики в ответ на поступающую сенсорную информацию. Таким образом, эксперименты подтверждают, что нейроны являются активными элементами,,которые генерируют сложные и даже хаотические колебания [34], а для нейронных сетей характерны такие явления как синхронизация и возбуждение в виде паттернов электрической активности [17,23,29], распространение уединенных волн - нервных импульсов [18,35-37], вынужденные колебания и синхронизация под внешним (сенсорным) воздействием [38] и др.

В настоящее время нелинейно-динамический подход для исследования свойств нейронных сетей активно развивается как в России (Безручко Б.П., Борисюк Г. Н., Борисюк Р. М., Иваницкий Г.Р., Казанович Я. Б., Казанцев В.Б., Некоркин В. И., Осипов Г.В., Потапов А.Б., Смирнов Д. А., Яхно В. Г. и др.), так и за рубежом (Abarbanel H.D.I., Afraimovich V. S., Aihara К., Arena P., Ayers J., Bazhenov M., Bilbaut J.M., Courbage M., Ermentrout G.B., Fortuna L., Ijspreert A. J., Izikevich E.M., Kurths J., Llinas R., Makarenko V. I., Rabinovich M. I., Rinzel J., Rulkov N. F., Sanjuan M. A. F., Tanaka G., Terman D. H., Varona P., Velarde M. G. и др.)

Существует достаточно большое число различных подходов построения моделей нервных клеток. Остановимся кратко лишь на двух из них, наиболее близких к теме диссертации. В основе первого подхода лежит детальное описание физико-химических процессов в клеточной мембране. Здесь одной из первых, и на сегодня общепризнанной, считается модель Ходжки-на и Хаксли [37], описывающая распространение электрической активности вдоль гигантского аксона кальмара. Модели такого типа, построенные на основе формализма Ходжкина-Хаксли (membrane conductance-based model - модель мембранной проводимости), базируются на описании динамики различных ионных трансмембранных токов [39] (токов протекающих через мембрану нервной клетки) являются высокоразмерными нелинейными системами дифференциальных уравнений и содержат большое число параметров, что делает довольно затруднительным их анализ. Существует и другой подход, который используется в данной работе, основанный на феноменологическом описании явлений. Здесь модели воспроизводят основные динамические свойства нервных клеток, не претендуя на количественное и детальное описание всех биофизических процессов. Например, большую популярность получила модель ФитцХью-Нагумо [40,41], в которой используется лишь две переменные - мембранный потенциал, и, так называемая, восстанавливающая переменная, характеризующая совокупное действие всех ионных токов. Широкое распространение таких моделей [42-50] связано с их относительно простой структурой (размерность систем не выше трех), и как правило, феноменологические модели обладают автоколебательными свойствами и могут быть адекватно описаны автогенераторными системами.

Сегодня аппарат нелинейной динамики широко применяется для исследования нейронных систем. Многие виды нейронной активности получили интерпретацию в терминах теории динамических систем. Такие понятия как регулярные и хаотические аттракторы, устойчивость, области притяжения, бифуркации [51-60] и др., прочно вошли в обиход целого направления, связанного с изучением различных аспектов нейронной активности методами нелинейной динамики. На данный момент, наиболее изучены такие явления как возникновение импульсов возбуждения и распространение их вдоль нервных волокон [61-72]. Изучение активности нейронных сетей, состоящих из большого числа нейронов, представляет собой сложную задачу . Такие сети содержат огромное количество нервных клеток (мозг человека содержит 1011 нейронов), взаимодействующих при помощи различных типов связей свойства которых, например, могут меняться во времени [74-77]. Кроме того нейронные сети, как правило обладают сложной пространственной организацией связей, обеспечивающей локальную и нелокальную связь одного нейрона со многими.

Интерес к исследованию нейронных систем обусловлен так же и прикладными аспектами, связанными с разработкой систем распознавания, передачи и обработки информации [39,78-81], построения систем управления и контроля, основанных на принципах нейродинамики [82-87]. Систематизация, обобщение и понимание закономерностей возникновения, распространения и исчезновения активности нейронных сетей являются важными задачами нелинейной динамики. Ключевую роль при решении этих задач играют, так называемые, нейроморфные модели. Модели учитывают морфологию индивидуальных нейронов и межнейронных связей, имеют архитектуру отражающую строение реальных нейронных систем. В нейроморфных моделях элементы-нейроны наделяются своей внутренней и, в частности, колебательной динамикой, которая в той или иной степени близка к динамике "живых" нейронов. Современные электронные технологии, в частности, так называемые программируемые кристаллы (программируемые логические интегральные схемы - ПЛИС, программируемые блоки аналоговых элементов - field-programmable analog array, большие интегральные схемы - БИС и др.), допускают создание компактных электронных схем из большого числа активных единиц, которые, могут воспроизводить базовые свойства нейронов. В связи с этим, построение новых нетрадиционных высокоэффективных систем обработки и хранения информации, координации и управления движением, основанных на нейродинамических принципах, стало актуальной и практически реализуемой задачей.

На сегодняшний день построено достаточно большое число электронных моделей автогенераторов [88-93], воспроизводящих динамику как отдельных нейронов, так и их ансамблей [94-98]. Существующие в настоящее время электронные нейроморфные модели, условно можно разделить натри основных типа: аналоговые системы [88-94,96-98] (оперирующие аналоговыми сигналами как и живые нейронные системы), аналого-цифровые [95] (модели нейронов являются аналоговыми, а информационные данные обрабатываются и передаются« цифровым образом) и цифровые (реализованные только на цифровых микросхемах). При реализации больших нейронных ансамблей используются различные полупроводниковые технологии построения интегральных схем [88,92,93,96,97], например, такие: КМОП (комплементарная логика на транзисторах металл-оксид-полупроводник, в англоязычной литературе CMOS - complementary metal-oxide-semiconductor) , МОП (металл - оксид - полупроводник, в англоязычной литературе MOSFET - metal-oxide-semiconductor field-effect transistor), BiCMOS (bipolar transistor and the CMOS) с использованием биполярных и КМОП-транзисторов на одном кристалле.

Отметим еще одно перспективное направление использования генераторных электронных моделей. Сравнительно недавно была показана возможность соединения реальных нейронов и "искусственных" [99-101]. Такие исследования могут помочь при диагностике и прогнозе динамических заболеваний (подавление эпилептических расстройств, связанных с нарушением колебательных ритмов в мозге и др.)

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

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

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

Научная новизна работы

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

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

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

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

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

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

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

Основные положения, выносимые на защиту

1. Разработанные автогенераторы обладают колебательными, возбудимыми и мультистабильными режимами.

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

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

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

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

Теоретическая и практическая значимость результатов

Реализованы аналоговые автогенераторы, демонстрирующие широкий спектр динамических режимов, характерных для колебательных нейронов. Реализованы ансамбли в виде цепочек резистивно связанных автогенераторов, в которых экспериментально получено распространение импульсов возбуждения, волновых фронтов переключения и самоподдерживающихся хаотических структур спайковой> активности. Построен ансамбль, на основе сети взаимодействующих автогенераторов со свойством фазовой автопереустаноки, позволяющий формировать фазовые паттерны синхронных колебаний, с заданным фазовым сдвигом между генераторами. На основе этого ансамбля разработана многопараметрическая система управления шагающим роботом. Разработанные автогенераторы могут быть использованы в качестве базовых единиц при построении информационно-вычислительных устройств нового поколения, способных осуществлять параллельное преобразование больших потоков информации. Результаты работы позволяют дать практические рекомендации по выбору параметров для существования требуемых динамических режимов. Результаты работы могут быть использованы в учебном процессе ВУЗов: ННГУ, СГУ, МГУ - при обучении студентов по специальностям радиофизического-профиля.

Апробация результатов. Основные результаты докладывались на следующих российских и международных конференциях: седьмой научной конференции по радиофизике (Нижний Новгород 2003), на конференции молодых ученых "Нелинейные волновые процессы" (Научная школа "Нелинейные волны"2004, 2006, 2008, 2010 Нижний Новгород), на нижегородской сессии молодых ученых (естественнонаучные дисциплины 2004), на всероссийской конференции "Хаотические автоколебания и образование структур" (Саратов 2004, 2007), на международных симпозиумах International Symposium "Topical Problems of Nonlinear Wave Physics" (Nizhny Novgorod, Russia 2003; St.-Petersburg - Nizhny Novgorod, Russia 2005; Nizhny Novgorod, Russia 2008), International Workshop on Nonlinear Dynamics of Electronic Systems, (Dijon, France 2006; Tokushima, Japan 2007); Chaos and dynamics in biological networks (Cargese, France 2010).

Личный вклад автора. Все эксперименты выполнены лично автором. В совместных работах, теоретический анализ и интерпретация полученных результатов были выполнены с научным руководителем. Результаты опубликованы в работах [103-124].

Структура и объем работы. Диссертация состоит из введения, трех глав, заключения и списка цитируемой литературы. Общий объем составляет 153 страницы, включая 71 рисунок и список литературы.

 
Заключение диссертации по теме "Радиофизика"

3.4 Выводы

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

Заключение

Основными результатами диссертационной работы являются следующие:

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

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

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

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

5. Проведено экспериментальное и теоретическое исследование волно вых движений в двух типах ансамблей из шестнадцати резистивно связанных автогенераторов: в цепочке генераторов ФитцХыо-Нагумо и цепочке генераторов с двумя нелинейными проводимостями. В'таких ансамблях экспериментально установлено, распространение, импульсов возбуждения и волновых фронтов переключения. Показано, что в. зависимостиот значений параметров; волновые фронты и^ импульсы обладают различными: свойствами: либо импульсы возбуждения и волновые фронты переключения при столкновении друг с другом и границами ансамбля аннигилируют, либо демонстрируют частицеподобные свойства, т. е. отражаются при взаимодействии друг с другом, и от границ щепочки.

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

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

8. Экспериментально показано, что автогенератор ФитцХью-Нагумо обладает* эффектом фазовой автопереустановки, который позволяет управлять фазой колебаний. При подаче на генератор прямоугольного импульса фаза колебаний переустанавливается к новому значению, которое не зависит от момента воздействия внешнего импульса (или от начальной фазы), а определяется лишь интенсивностью импульса, т. е. его амплитудой и длительностью.

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

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

 
Список источников диссертации и автореферата по физике, кандидата физико-математических наук, Щапин, Дмитрий Сергеевич, Нижний Новгород

1. Strogatz S.H., Mirollo R.E. Splay states in globally coupled Josephson arrays: Analytical prediction on Floquet multipliers // Phys. Rev. E, 1993. V. 47. N.l. P. 220-227.

2. Лихарев К.К., Ульрих Б.Т. Системы с джозефсоновскими контактами. -М.: Изд-во МГУ. 1978. 446 с.

3. Otsuka К. Self-Induced Phase Turbulence and Chaotic Itenerancy in Coupled Laser Systems // Phys. Rev. Lett., 1990. V. 65. N. 3. P. 329332.

4. Winful H.G., Rahman L. Synchronized Chaos and Spatiotemporal Chaos in Arrays of Coupled Lasers // Phys. Rev. Lett., 1990. V. 65. N. 13. P. 1575-1578.

5. Капранов M.B. Взаимодействующие многосвязные СФС. Системы фазо-вой синхронизации / Под ред. Шахгильдяна В.В., Белюстиной JI.H. М.: Радио и связь. 1982. С. 55-73.

6. Afraimovich, V.S., Nekorkin, V.I., Osipov, G.V., Shalfeev, V.D. Stability, Structures and Chaos in Nonlinear Synchronization Networks. World Scientific. Singapore. 1995. 246 p.

7. Пиковский А., Розенблюм M., Курте Ю. Синхронизация. Фундаментальное нелинейное явление. М.: Техносфера. 2003. 496 с.

8. Самойленко В.И., Шишов Ю.А. Управление фазированными антенными решетками. М. Радио и связь. 1983. 408 с.

9. Llinas R., Yarom Y. Oscillatory properties of guinea-pig inferior olivary neurones and their pharmacological modulation: An in vitro study //J. Physiol. Lond., 1986. V. 376. P. 163-82.

10. Jensen O., Lisman J.E. An oscillatory short-term memory buifer model can account for data on the Sternberg task //J. Neurosci., 1998. V. 18. P! 10688.

11. Lisman J.E., Idiart M.A. Storage of 7 +/- 2 short-term memories in oscillatory subcycles // Science, 1995. V. 267. P. 1512.

12. Jensen O., Lisman J.E. Hippocampal sequence-encoding driven by a cortical multi-item working memory buffer // TRENDS in Neurosciences, 2005. V. 28. № 2. P. 67.

13. Haj-Dahmane S. Andrade R. Ionic Mechanism of the Slow Afterdepolarization Induced by Muscarinic Receptor Activation in Rat Prefrontal Cortex // J. Neurophysiol., 1998. V. 80 P. 1197-1210.

14. Raghavachari S., Lisman J. E., Tully M., Madsen J. R., Bromfield E. B. and Kahana M. J. Theta Oscillations in Human Cortex During a Working-Memory Task: Evidence for Local Generators // J. Physiol., 2006. V. 95(3). P. 1630-1638.

15. Lampl I., Yarom Y., Subthreshold oscillations and resonant behavior: two manifestations of the same mechanism // Neuroscience, 1997. V. 78. P. 325-341.

16. Murray J.D. Mathematical Biology. Springer-Verlag. Berlin. 1993. 767 p.

17. Winfree A.T. The geometry of Biological Time. Springer-Verlag. New-York. 1980.

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

19. Leznik E., Makarenko V.I., Llinas R. Electrotonically Mediated Oscillatory Patterns in Neuronal Ensembles: An In Vitro Voltage-Dependent Dye-Imaging Study in the Inferior Olive // J. Neurosci., 2002. V. 22. P. 28042815.

20. Абарбанель Г. Д. И., Рабинович М. И., Сильверстон А. , Баженов М. В., Хуэрта Р., Сущик М.М., Рубчинский JI.JI. Синхронизация в нейронных ансамблях // УФН, 1996. Т. 166, № 4. С. 363-390.

21. FitzHugh R. Impulses and physiological states in theoretical models of nerve membranes // Biophysical Journal, 1961. V. 1. P. 445-466.

22. Nagumo J., Arimoto S., Yoshizawa S. An active pulse transmission line simulating nerve axon // Proc. IRE, 1962. V.50 P. 2061-2070.

23. Morris C. and Lecar H. Voltage Oscillations in the barnacle giant muscle fiber // Biophysical Journal, 1981. V. 35(1). P. 193-213.

24. Rulkov N.F. Regularization of Synchronized Chaotic Bursts // Phys. Rev. Lett., 2001. V. 86. P. 183-186.

25. Izhikevich E.M. Neural excitability, spiking and bursting // Int. J. Bifurc. Chaos, 2000. V.10. N.6. P.1171-1266.

26. Hindmarsh J.L., Rose R.M. A model of neuronal bursting using three coupled first order differential equations // Proc. R. Soc. Lond. В Biol. Sci., 1984. V. 221. P. 87-102.

27. Клиныиов В.В., Некоркин В.И. Динамика системы с последеполяри-зацией и ингибиторной обратной связью // Изв. вузов. Радиофизика, 2005. Т. 48, № 3. С. 228-237.

28. Klinshov V.V., Nekorkin V.I. Activity clusters in dynamical model of the working memory // Network: Computation in Neural Systems, 2008. V.19, n. 2. P. 119-135.

29. Klinshov V.V., Nekorkin V.I. Working memory in the network of neuronlike units with noise // Int. J. Bifurcation Chaos, 2008. V. 18, № 1. P. 2743-2752.

30. Kazantsev V.B. , Nekorkin V.I., Makarenko V.I., Llinas R. Self-referential phase reset based on inferior olive oscillator dynamics // Proc. Natl. Acad. Sci. USA, 2004. V. 101(52). P. 18183-18188.

31. Velarde M.G., Nekorkin V.I., Kazantsev V.B., Makarenko V.A., Llinas R. Modeling Inferior Olive Dynamics // Neural Networks, 2002. V. 15(1). P. 5-10.

32. Баутин H. H., Леонтович E. А. Методы и приемы качественного исследования динамических систем на плоскости. М. Наука, 1976. 111 с.

33. Takens F. Singularities of vector fields // Publ. Math. IHES, 1974. V.43. P.47-100.

34. Богданов Р.И. Версальная деформация особой точки векторного поля на плоскости в случае двух нулевых собственных чисел // Труды семинара им. И.Г. Петровского, в.2. М.: Изд. МГУ, 1976. С.37-65.

35. Марсден Дж., Мак Кракен М. Бифуркация рождения цикла и её приложения. М. Мир. 1980.

36. Белых В.Н., Некоркин В.И. О качественном исследовании многомерной фазовой системы // Сибирский матем. журнал, 1977. Т. 18, №4. С. 723-735.

37. Belykh V.N. Homoclinic and heteroclinic linkages in concrete systems: nonlocal analysis and model maps // Amer. Math. Soc. Transl., 2000. V. 200. P. 51-62.

38. Митропольский Ю.А., Лыкова О.Б. Интегральные многообразия в нелинейной механике. М. Наука. 1973.

39. Рабинович М.И., Трубецков Д.И. Введение в теорию колебаний и волн. М. Наука. 1984.

40. Andronov А.А. Leontovich Е.А., Gordon I.I., Mayer A.G., Kuznetsov Y.A. Elements of Applied Bifurcation Theory. Springer-Verlag. Berlin. 1995.

41. Андронов A.A., Витт A.A., Хайкин С.Э. Теория колебаний. М. Физ-матгиз. 1959.

42. Васильев В.А., Романовский Ю.М., Яхно В.Г. Автоволновые процессы. М. Наука. 1987.

43. Морнев О.А., Асланиди О.В., Алиев P.P., Чайлахян JI. М. Солитон-ный режим в уравнении ФитцХыо-Нагумо: отражение сталкивающихся импульсов возбуждения. ДАН. 1996. Т. 347. С. 123-125.

44. Hayase Y. Collision and Self-Replication of Pulses in a Reaction Diffusion System // J. of the Physical Society of Japan, 1997. V. 66, N. 9. P. 25842587.

45. Hayase Y., Ohta T. Self replicating pulses and Sierpinski gaskets in exitable media // Phys. Rev. E, 2000. V. 62, N.5. P. 5998-6003.

46. Nekorkin V.I., Velarde M.G. Sinergetic phenomena in active lattices. Springer-Verlag. 2002.

47. Perez-Munuzuri V., Perez-Villar V., Chua L.O. Propagation failure in linear arrays of chua's circuits // Int. J. of Bifurcation and Chaos, 1992. V. 2(2). P. 403-406.

48. Keener J.P. Propagation and its failure in coupled systems of discrete excitable cells // SIAM Journal on Applied Mathematics, 1987. V. 47. P. 556-572.

49. Comte J.C., Morfu S., Marquie P. Propagation failure in discrete bistable reaction-diffusion systems: Theory and experiments // Phys. Rev. E, 2001. V. 64. P. 027102.

50. Binczak S., Bilbault J.M. Experimental propagation failure in a nonlinear electrical lattice // Int. J. Bifurcation and Chaos, 2004. V. 14(5). P. 18191830.

51. Nekorkin V.I., Chua L.O. Spatial disorder and wave fronts in a chain of coupled Chua's circuits // Int. J. Bifurcation and Chaos, 1993. V. 3. P. 1281-1297.

52. Максимов А.Г., Некоркин В.И. Гетероклинические траектории и фронты сложной формы модели ФитцХыо-Нагумо // Математическое моделирование, 1990. Т. 2, №2. С. 129-142.

53. Nekorkin V.I., Kazantsev V.B., Velarde M.G. Mutual synchronization of two lattices of bistable elements // Phys. Lett. A, 1997. V. 236. P. 505-512.

54. Некоркин В. И. Нелинейные колебания и волны в нейродинамике // УФН, 2008. Т. 178. № 3. С. 313-323.

55. Николлс Дж., Мартин Р., Валлас В., Фукс П. От нейрона к мозгу. Изд-во УРСС. 2003. 672 с.

56. Bennett M.V., Zukin R.S. Electrical coupling and neuronal synchronization in the mammalian brain // Neuron. 2004. V. 41(4). P. 495-511.

57. Kandel E.R., Schwartz J.H., Jessell T.M. (Eds.) Principles of Neural Science. Third Edition. Prentice-Hall Intern. Inc. 1991, 1135 p.

58. Korn H., Faure P. Is there chaos in the brain? II. Experimental evidence and related models. // C. R. Biologies, 2003. V. 326. P. 787-840.

59. Koch C. Biophysics of computation: information processing in single neurons. Oxford University Press. 1998.

60. Казанцев В.В., Некоркин В.И. Динамика колебательных нейронов. Информационные аспекты // Нелинейные волны 2002.- Нижний Новгород: ИПФ РАН. ред. А.В. Гапонов-Грехов, В.И. Некоркин. 2003.

61. Борисюк Г.Н., Борисюк P.M., Казанович Я.Б., Иваницкий Г.Р. Модели динамики нейронной активности при обработке информации мозгом -итоги "десятилетия"// УФН, 2002. Т. 172, №10. С. 1189.

62. Потапов А.Б., Али М.К. Нелинейная динамика обработки информации в нейронных сетях // Новое в синергетике: Взгляд в третье тысячелетие.-М.: Наука, ред. Г.Г. Малинецкий, С.П. Курдюмов. 2002.

63. Borisyuk R.M., Kazanovich Y.B. Oscillatory neural network model of attention focus formation and control // BioSystems, 2003. V. 71. P. 29-38.

64. Grillner S. Neural networks for vetebrate locomotion // Scientific American, 1996. V. 274 (1). P.64-69.

65. Ekeberg O. A., A combined neuronal and mechanical model of fish swiming // Biol. Cybern., 1993. V. 69. P. 363.

66. Казанцев В. В., Некоркин В. И., Принципы контроля и координации движений на основе динамики нейронов головного мозга // Известия ВУЗов Прикладная Нелинейная Динамика, 2001. Т. 9, № 1. С. 38-48.

67. Kazantsev V.B., Nekorkin V.I., Makarenko V.I., Llinas R. Olivo-cerebellar cluster-based universal control system // Procs. Natl. Acad. Sci. USA, 2003. V. 100(22). P. 13064-13068.

68. Allen I. Selverston, Attila Szucs, Ramon Huerta, Reynaldo Pinto and Marcelo Reyes Neural mechanisms underlying the generation of the lobster gastric mill motor pattern // Frontiers in Neural Circuits, 2009. V. 3. Article 12.

69. Kohno T., Aihara K. A MOSFET-based model of a class 2 nerve membrane // IEEE Trans. Neural Networks, 2005. V. 16. P. 754-773.

70. Binczak S., Kazantsev V.B., Nekorkin V.I., Bilbault J.M. Experimental study of bifurcation in modified Fitzhugh-Nagumo cell // Electron. Lett., 2003. V. 39. P. 961-962.

71. Wagemakers A., Sanjuan M.A.F., Casado J. M., Aihara K. Building electronic bursters with the Morris-Lecar neuron model // Int. J. Bifurcation and Chaos, 2006. V. 16, n. 12. P. 3617-3630.

72. Linares-Barranco B., Sanchez-Sinencio E., Rodriguez-Vazquez A., Huertas J. L. A CMOS Implementation of FitzHugh-Nagumo Neuron Model // IEEE Journal of solid-state circuits, 1991. V. 26(7). P. 956.

73. Jacobo D. Sitt and J. Aliaga. Versatile biologically inspired electronic neuron // Phys. Rev. E, 2007. V.76. P. 051919.

74. Simoni M.F., Cymbalyuk G.S., Sorensen M.E., Calabrese R.L., DeWeerth S.P. A Multi-Conductance Silicon Neuron with Biologically Matched Dynamics // Biomedical Engineering, IEEE Transactions, 2004. P. 342354.

75. Van Schaik A. Building blocks for electronic spiking neural networks // Neural Networks, 2001. V. 14. P. 617-628.

76. Douence V., Renaud-Le Masson S., Saighi S. and Le Masson G. A Field-Programmable Conductance Array IC for Biological Neurons Modeling // Bio-Inspired Applications of Connectionism 6th International Work

77. Conference on Artificial and Natural Neural Networks. IWANN 2001 Granada. Spain. Proceedings. Part II. J. Mira and A. Prieto (Eds.) Lecture Notes in Computer Science. Springer-Verlag Berlin. 2001. V. 2085/2001. P. 31-38.

78. Le Masson S., Laflaquiere A., Bal T., Le Masson G. Analog circuits for modeling biological neural networks: Design and applications // IEEE Trans. Biomed. Eng., 1998. V. 46. P. 638-645.

79. Arthur J.V., Boahen K. Silicon Neurons that Inhibit to Synchronize // Circuits and Systems; 2007. ISCAS 2007. IEEE International Symposium. P. 1186.

80. Marquie P., Comte J.C., Morfu S. Analog simulation of neural information propagation using an electrical FitzHugh-Nagumo lattice // Chaos, Solitons and Fractals, 2004. V. 19. P. 27-30.

81. Szucs A., Varona P., Volkovskii A. R., Abarbanel H. D. I., Rabinovich M. I., Selverston A. I. Interacting biological and electronic neurons generate realistic oscillatory rhythms // Neurorep., 2000. V. 11. P. 563-569.

82. Zeck G., Fromherz P. Noninvasive neuroelectronic interfacing with synaptically connected snail neurons immobilized on a semiconductor chip // Proc. Natl. Acad. Sci. USA, 2001. V. 98. P. 10457-10462.

83. Sorensen M., DeWeerth S., Cymbalyuk G., and Calabrese R. L. Using a Hybrid Neural System to Reveal Regulation of Neuronal Network Activity by an Intrinsic Current // The Journal of Neuroscience, 2004. V. 24(23). P. 5427-5438.

84. Cole K. S. and Baker R.F. Longitudinal impedance of the squid giant axon //J. Gen. Physiol., 1941. V. 24. P. 771.

85. Prestunov V. V., Shapin D. S., Kazantsev V. В., "and Nekorkin V. I. The dynamics of two oscillators with pulse-controlled coupling. Abstr. Int. Symp. "Topical Problems of Nonlinear Wave Physics". Nizhny Novgorod. Russia. 2003. P. 74-75.

86. Щапин Д. С., Престунов В. А., Казанцев В. В., Некоркин В. И. Модель нейрона с подпороговыми колебаниями // Тезисы докл. Конф. молодых ученых "Нелинейные волновые процессы (XII науч. школа "Нелинейные волны-2004"). Н. Новгород, 2004. С. 92-93.

87. Престунов В. А., Щапин Д. С. Динамика модели нейрона с подпоро-говой ак-тивностью // Тезисы докладов IX Нижегородской сессии молодых ученых (естествен-нонаучные дисциплины), Н. Новгород, 2004. С.120.

88. Shapin D. S., Kazantsev V. В., and Nekorkin V. I. Synchronization and time bind-ing in two neuronal oscillators with inhibitory feedback // abstr. Int. Symp. "Topical Prob-lems of Nonlinear Wave Physics", Nizhny Novgorod, 2005. P. 95-96.

89. Щапин Д. С. Фазовая автопереустановка и синхронизация релаксационных автоколебательных систем // Тезисы докл. X Нижегородской сессии молодых ученых (естественнонаучные дисциплины). Н. Новгород, 2005. С. 160-161.

90. Некоркии В. И., Дмитричев А. С., Щапин Д. С., Казанцев В. Б. Динамика модели нейрона со сложно-пороговым возбуждением // Математическое моделирование, 2005. Т. 17, №6. С. 75-91.

91. Щапин Д. С. Импульсное блокирование связи в сети нейроноподоб-ных элементов в физическом эксперименте. Тезисы докладов конференции молодых ученых "Нелинейные волновые процессы". Н. Новгород. 2006. С. 169.

92. Kazantsev V. В., Dmitrichev A. S., Shapin D. S. and Nekorkin V. I. Multi-threshold excitability in a nonlinear network of neuron-like units // The 14-th International Workshop on Nonlinear Danamics of Electronic Systems. Dijon. Prance. 2006. P. 69-72.

93. Kazantsev V. В., Shapin D. S. and Nekorkin V. I. Phase clusters in oscillatory neuronal network with inhibitory feedback // The 14-th International Workshop on Nonlinear Da-namics of Electronic Systems. Dijon. France. 2006. P. 73-77.

94. Shapin D. S. Phase control in nonlinear electronic circuit modeling oscillatory neurons // The 15th International Workshop on Nonlinear Dynamics of Electronic Systems. Tokushima. Japan. 2007. P. 194-196.

95. Щапин Д. С. Управление фазой колебаний в нелинейной электронной схеме модели нейрона // 8-я международная школа "ХАОС-2007": Тез. докл.- Саратов: Изд. СГУ, 2007. С. 100.

96. Некоркин В. И., Щапин Д. С., Дмитричев А. С. Сложная волновая динамика ансамбля нейроноподобных элементов со сложно-пороговым возбуждением // Изв. ВУЗов. Прикладная нелинейная динамика, 2007. Т. 15, т. С. 3-22.

97. Щапин Д. С. Динамика двух нейроноподобных элементов с подавляющей обратной связью // Радиотехника и электроника, 2009. Т. 54, №2. С. 185-195.

98. Клинынов В. В., Щапин Д. С., Некоркин В. И. Моделирование нейродинами-ческой системы рабочей памяти // Радиотехника и электроника, 2010. Т. 55, т. С. 812-817.