Стабилизация турбулентной спирально-волновой динамики возбудимых сред тема автореферата и диссертации по физике, 01.04.02 ВАК РФ
Высоцкий, Семен Андреевич
АВТОР
|
||||
кандидата физико-математических наук
УЧЕНАЯ СТЕПЕНЬ
|
||||
Москва
МЕСТО ЗАЩИТЫ
|
||||
2010
ГОД ЗАЩИТЫ
|
|
01.04.02
КОД ВАК РФ
|
||
|
На правах рукописи —-
004603988
Высоцкий Семен Андреевич
Стабилизация турбулентной спирально-волновой динамики возбудимых сред
Специальность 01.04.02 — теоретическая физика
АВТОРЕФЕРАТ
диссертации на соискание ученой степени кандидата физико-математических наук
1 7 «ЮН 2010
Москва - 2010 г.
004603988
Работа выполнена на физическом факультете
Московского государственного университета имени М.В. Ломоносова
Научный руководитель:
доктор физико-математических наук, профессор Лоскутов Александр Юрьевич
Официальные оппоненты:
доктор физико-математических наук, главный научный сотрудник Волков Евгений Израилевич доктор физико-математических наук, профессор Розов Николай Христович
Ведущая организация:
Институт космических исследований Российской академии наук
Защита состоится г^г^......... 2010 года в /.$.. ч.......м. на заседании
Диссертационного Совета Д 501.002.10 в Московском государственном университете по адресу: 119991, ГСП-2, Мос^в^^енинские горы, МГУ им. М.В. Ломоносова, физический факультет, ауд. ...
С диссертацией можно ознакомиться в библиотеке физического факультета МГУ им. М.В. Ломоносова.
Автореферат разослан «. ............2010 г.
Ученый секретарь Диссертационного совета Д 501.002.10
в Московском государственном университете __
доктор физико-математических наук, профессор Ю.В. Грац
Общая характеристика работы
Актуальность темы
Распределенные среды представляют собой достаточно широкий класс, включающий самые разные математические, физические, химические биологические и другие системы. Любой объект, если его компонента распределена в пространстве или времени по некоторому закону, является распределенным. К таким объектам относятся, например, жидкость или газ, популяция какого-либо вида, человеческий мозг или проводящая ткань сердца. Не удивительно, что уже очень давно данный класс систем привлекает исследователей из самых разных областей.
Наиболее сложная проблема теории распределенных сред — управление их динамикой и, в частности, стабилизация (подавление) сложных режимов поведения (квазипериодических или хаоса). Эта проблема возникла достаточно давно и связана она с тем, что такие сложные режимы, как правило, являются крайне нежелательными. В особенности это касается реальных физических, экологических, химических или биологических систем.
Под стабилизацией неустойчивого или хаотического поведения динамических систем обычно понимается искусственное создание и поддержание в этих системах устойчивых (как правило, периодических) колебаний посредством внешнего воздействия. Эта задача, не смотря на простоту формулировки, оказывается весьма сложной научной проблемой, особенно в приложении к распределенным средам.
Актуальность такой задачи вполне очевидна. Рассмотрим несколько реальных примеров. В приложении к сердечной ткани выведение системы на требуемый режим дает возможность влиять на сердечный ритм. Дело в том, что в настоящее время в теории возбудимых сред доминирует гипотеза, согласно которой возникновение фатальных сердечных аритмий — фибрилляций — есть следствие рождения в сердечной ткани большого количества автоволновых источников: спиральных волн или вихревых структур (т.е. пространственно-временного хаоса). Современные
методы стабилизации таких режимов с помощью одиночных электрических импульсов (в том числе от имплантируемых дефибрилляторов) являются весьма жесткими и далеко не всегда приводят к успеху. Однако исследования самого последнего времени открывают новые возможности. Оказывается, что турбулентный режим во многих возбудимых средах может быть стабилизирован достаточно слабым периодическим параметрическим, или силовым воздействием, приложенным к некоторой области среды. Для реакции Белоусова-Жаботинского такое воздействие позволяет создавать структуры нужного вида, получая таким образом системы, способные распознавать образы. Помимо перечисленных здесь идей, результаты настоящей работы можно использовать для кодирования информации, создания когерентных структур заданной геометрии, управления потоком частиц и других многочисленных приложений современной теоретической физики.
Цели работы.
— Построение математических моделей распределенных сред с различными граничными и начальными условиями, состоящих из возбудимых элементов.
— Исследование динамики системы в том числе со слабым почти точечным внешним воздействием.
— Анализ поведения системы в зависимости от параметров возбудимости среды, а также характеристик внешнего возбуждения.
— Стабилизация сложных режимов поведения, связанных со спирально-волновой турбулентностью.
— Изучение проблемы стабилизации турбулентной динамики силовым воздействием, применяемым ко всем точкам среды.
Научная новизна полученных в диссертации результатов
1. Исследована динамика математических моделей распределенной среды с внешним воздействием при различных граничных и начальных условиях и в широком диапазоне параметров как среды так и внешнего
воздействия.
2. На основе теории динамических систем и новых результатов теории распределенных сред показана возможность управления динамикой представленных моделей.
3. Показана принципиальная возможность выведения системы из состояния сложной (в том числе хаотической) динамики слабым почти точечным воздействием.
4. Исследована возможность одновременного воздействия на все элементы среды, указаны преимущества и недостатки такого подхода.
5. Предложен новый метод подавления спирально-волновой турбулентности с помощью движущихся ведущих центров, позволивший заметно повысить его эффективность. Получена зависимость эффективности подавления от количества ведущих центров и их характеристик.
Практическая ценность работы
1. Показана принципиально новая возможность подавления спирально-волновой турбулентности возбудимой среды внешним точечным воздействием малой амплитуды и выведения ее на периодический режим движения.
2. Разработаны практические методы, позволяющие вычислять наиболее предпочтительные для такой стабилизации частоты.
3. Найдено, что нестационарный пейсмекер, расположенный в среде, значительно повышает эффективность предложенного метода.
Структура и объем работы
Диссертация состоит из введения, общих сведений (обзора литературы), четырех глав, заключения и списка цитируемой литературы из 173 наименований. Диссертация содержит 126 страниц текста, включая рисунки, оглавление и список литературы.
Публикации
По результатам данной диссертационной работы опубликовано четыре статьи в реферируемых журналах (две — в международных и две —
в российских), а так же одна статья в сборнике, посвященном исследованию распределенных сред.
Апробация работы
Основные результаты исследования были доложены на на многочисленных международных конференциях. В их числе 17-й международный симпозиум в Швейцарии в июне 2009 года, международный симпозиум, посвященный сложным динамическим системам и приложениям в Индии в 2009 году. Дважды результаты докладывались на конференции РЬуэСоп в Санкт-Петербурге в 2005 году и в Германии в 2007 году.
Основное содержание диссертации
Во введении диссертации обоснована актуальность темы исследований и сформулированы цели работы.
Первая глава посвящена общим сведениям. Здесь содержится обзор литературы, формализуется проблема внешнего воздействия на динамическую систему и подавления нерегулярных режимов поведения при помощи такого воздействия. Далее приводятся некоторые существующие на сегодняшний день точные результаты из теории нелинейных сред с диффузией. Рассмотрены как дискретные, так и непрерывные модели распределенных сред. В частности приводятся простая модель ФитцХью-Нагумо, ее модификация Панфилова-Хогевега и более сложная трехкомпонеитная модель Фентона-Кармы.
Модель ФитцХыо-Нагумо обычно записывается в виде двухкомпо-нентной системы:
— = Аи - /(и) - у,
а (1)
= д{и,ь)(ки-у).
Здесь и(х, у, ¿) — концентрация активатора, у(х, у, £) — концентрация ингибитора. Легко видеть, способностью к диффузии обладает лишь активатор, он и служит прообразом потенциала действия, который используется в многочисленных приложениях.
Вид функций f(u) и д(и, и) обычно выбирают таким образом, чтобы
полученные при интегрировании профили и(х,у, ¿) максимально соответствовали профилям потенциала действия, полученным экспериментально. Например, в модификации Панфилова-Хогевега эти функции выбраны в кусочно-линейном виде:
С\и, и < щ, ГС?1> и < «2,
/(и) = < -С2и + а, и в [щ,и2), д(и,ь) = < С2, и>и2,
С3(и- 1), и>и2, [(?з, и<щ,и<и1.
(2)
В частности, значения параметров, приближающих систему к модели сердечной ткани и найденные эмпирически, таковы:
Сх = 20 «1 = 0.0026 а = 0.06
С2 = 3 и2 = 0.837 к = 3 (3)
С3 = 15 = 1.8 С2 = 1.0
Динамика активатора определяется функцией /, динамика ингибитора — функцией д. Параметры (?1 и бз являются управляющими и определяют относительный и абсолютный рефрактерный период соответственно.
Для количественного анализа возбудимых сред существуют более сложные модели, такие, например, как модель Фентона-Кармы. Она записывается в виде трехкомпонентной системы:
= - Лга(и, у) - Ми) -
Эр = в(ис - и)(1 - у)/т-(и) - в(и - ис)у/т+, (4)
= в(ис - и)(1 - ю)/т~ - в(и - ис)ю/т+.
Здесь и = , т~(и) = в (и - + 9(и„ - - пере-
менные состояния, 9(х) — стандартная функция Хевисайда. У^,,, и — токи натрия, калия и кальция, соответственно. Они выражаются следующим образом:
V
JNa{u,v) =--0(и - ис)(1 - и)(и - Ис),
^ 1 и 1
^к(и) = —0(ис - и) + —в(и - ис), (5)
•Ы«,«0 = -^(1 + Л[*(и - и?0)]).
2ГСа
Эта модель была предложена Фентоном (Fenton) и Кармой (Karma) в 1998 году. Ее цель не в том, чтобы точно копировать микроскопическую ионную сложность тех или иных процессов, а скорее в том, чтобы воспроизвести динамику потенциала действия на мезоскопическом уровне (на шкалах между внутриклеточным уровнем и целым органом, если речь идет о моделировании сердечной ткани). Модель построена так, чтобы ввести минимальный набор ионных мембранных токов, необходимых для воспроизведения экспериментально полученных кривых восстановления, и, таким образом, состоит из трех переменных.
Данные системы в диссертационной работе рассматриваются не только с точки зрения формального перечисления их характеристик и свойств, но также и с точки зрения возникающих в них решений. К таким решениям в частности относятся известные автоволновые структуры — спиральные волны.
Спиральная волна, как следует из названия, представляет собой самостоятельно вращающуюся в среде автоволновую структуру спиральной формы. Ядро спиральной волны (центральный сегмент квадрата на рис. 1) представляет собой разрыв фазы (или фазовую сингулярность). Фаза обычно определяется через значения функций и(х, у, t) и v(x, у, t) в точке и их среднего по всей области интегрирования в каждый момент времени и* nv":
Рассмотрим интеграл, через который определяется так называемый «топологический заряд»:
Если замкнутый контур С окружает точку фазовой сингулярности, то интеграл обратится в ±27г (знак определяет хиральность). В противном случае он будет равен 0 . Обозначим градиент фазы через вектор Уф = к. Понятно, что величина к обращается в 0 везде, где ф дифференцируема, и не равен 0 как раз в точках фазовых сингулярностей. Компоненты к
(б)
с
можно представить в виде конечных разностей на сетке:
кх[т,п] = ф[т+ 1 ,п] — ф[т, п], ку[т, п] = ф[т, п + 1] — ф[т, п].
Порядок суммирования градиента фазы по ячейкам показан на рис. 1:
Рис. 1 Порядок суммирования значений градиента фазы при поиске фазовых сингулярностей.
Интеграл (6) можно выразить через так называемый оператор конво-люции:
п = [к х V] ос Ух ® кх + ® ку
Здесь ® — оператор конволюции, а Ут и — ядра конволюции. Конкретный вид ядер зависит от того, каким методом будет аппроксимироваться интеграл и по какому количеству ячеек будет проводиться интегрирование. В диссертационной работе мы использовали ядра следующего вида:
" -1/2 0 1/2 V* = -10 1 _ -1/2 0 1/2
что соответствует суммированию по восьми соседним ячейкам.
Здесь необходимо отметить, что вовсе не обязательно любой фазовой сингулярности соответствует спиральная волна (обратное, однако, верно). Сингулярности могут появляться в среде также в виде небольших «обрывков» волнового фронта и исчезать.
Временная динамика фазовых сингулярностей впоследствии будет использована как критерий развития в среде пространственно-временного хаоса и подавления спирально-волновой турбулентности.
Главы со второй по пятую включительно содержат оригинальные результаты.
1/2 0 1/2 0 0 0 -1/2 0 -1/2
Во второй главе приводится аналитическое исследование системы (1). Рассматривается система с внешним воздействием. Обсуждаются оптимальные параметры такого воздействия и методики их подбора. В данной главе также показана принципиальная возможность выведения системы из состояния спирально-волновой турбулентности слабым почти точечным воздействием.
Рассмотрим вначале сосредоточенную систему, образованную из системы (1) без внешнего воздействия путем отбрасывания оператора Лапласа. Тогда получим следующие два простых уравнения:
й = -/(«) - V,
V = g(u,v)(3u — v).
Функции f(u) и д(и, и) имеют тот же вид, что и в соотношении (2), то есть состоят из трех линейных участков. Таким образом, фазовая плоскость системы (7) делится прямыми и = щ, и = щ и v = vi на три части, которые мы обозначим I, II и III. Целесообразно рассматривать исследуемую систему отдельно в каждой из указанных областей. Перепишем систему (7) в векторном виде:
(й, v) = F (и, v),
Рассчитаем значения дивергенции потока F в каждой из трех областей. Несложно показать, что
div(Fj) = -20 - G3 < 0, div(Fjj) = 3 - Gi > 0, div(F///) = -15 - G2 = -16 < 0.
Таким образом, в области III система консервативна, но, как видно из рис. 2, в целом она оказывается диссипативной, т.к. из области III фазовые траектории всегда уходят либо в область I либо в область II. Анализ на стационарные точки показывает, что в области I есть одна такая точка (0,0), в области II стационарных точек нет, и в области III — еще одна точка (5/6,5/2). Собственные значения матрицы линеаризации — действительные отрицательные числа: Ai,2|(o,o) = 1/2(—20 — G3±
\/(20 + бз)2 - 92С3) (где С3 > 0), и Л1>2|(5/б,5/2) » -14.8, -1.2 . То есть обе стационарные точки являются устойчивыми узлами.
В касательном пространстве системы имеется три независимых векторных поля скоростей (й, ь) и довольно сложно предсказать, как поведет себя фазовая траектория с начальными условиями, близкими к границам между областями 1,11 и III. Значения управляющих параметров были фиксированы: 61 = 0.01 и бз = 0.1. На представленном ниже рисунке показано векторное поле скоростей системы (7) (и, г>) в фазовом пространстве (и, у), и показана одна фазовая траектория с начальным
Рис. 2 Поле скоростей системы (7) в фазовом пространстве. Жирной, полужирной и тонкой линиями обозначены части поля I, II и III, соответственно. Узлы обозначены маленькими спиралями.
Добавим теперь к компоненте и системы (7) оператор Лапласа и перейдем, таким образом, от сосредоточенной к реакционно-диффузионной системе. Для того чтобы понять как именно появление лапласиана влияет на поведение среды, введем малый коэффициент а как:
^ = а£\и ~ f{u) - v,
^ = g(u,v)(ku - v).
Численное исследование этой системы при G1 = 0.01 и G3 = 0.1 показало, что при а ^ 0.01 распространение волн в среде невозможно -
волновые фронты затухают. При а = 0.2 наблюдается процесс дестабилизации и распада волнового фронта одиночной спиральной волны, т.е. волновой фронт претерпевает разрыв и постепенно распадается на небольшие куски. Они, в свою очередь, могут либо затухать, либо образовывать вторичные спирали, которые также постепенно разрушаются.
Вообще говоря, не вполне ясно, необходимо ли условие неустойчивости волнового фронта, чтобы можно было утверждать, что система находится в состоянии спирально-волнового хаоса. Простейшая неустойчивость появляется, когда спиральпые волны начинают блуждать при переходе (названный вторичной бифуркацией Хопфа) от вращающейся с одним периодом спиральной волны к двойной периодической системе, когда спираль вращается с одним периодом, а ее кончик (фазовая сингулярность) прецессирует с другим периодом. Следующая бифуркация имеет место при переходе от квазипериодического блуждания к хаотическому, при котором траектория фазовой сингулярности сильно нерегулярна. В этом случае появляется новый режим, когда непрерывно создаются и уничтожаются множественные вращающиеся волны.
Для того чтобы понять процесс возникновения пространственно-временного хаоса в исследуемой системе (1), необходимо сначала пояснить процесс возникновения спиральных волн. Существуют разные способы искусственного создания спиралей. Например, можно непосредственно задать ее в виде архимедовой спирали г = а/2р, где а — постоянный шаг спирали. В данной работе используется способ получения спирали из плоской полуволны. Принцип, лежащий в основе данной методики, очень простой. В свободную среду помещается прямой отрезок волнового фронта. Вследствие диффузии он начинает постепенно прорастать, закручиваться и принимать спиральную форму с обоих концов. Если теперь уменьшить параметр С?1, то волновые фронты постепенно будут терять устойчивость и со временем начинают разрушаться. Этот случай показам на рис. 3.
Что касается внешнего воздействия, то параметрическое воздействие в рассматриваемом классе систем не всегда осуществимо, т.к. оно пред-
полагает изменение управляющих параметров во всех точках среды одновременно. Поэтому мы изучали возможность точечного подавления путем подвода коротких импульсов к определенному участку среды, что — и это очень важно — легко реализовать на практике.
Рис. 3 Образование спирали из плоской полуволны с последующим ее разрушением и переходом к хаосу, О] = 0.01333, Сз = 1.
Нами было использовано несколько вариантов: синусоидальная форма, однофазные прямоугольные импульсы, бифазные прямоугольные импульсы, постоянное «напряжение», бифазный пилообразный импульс различных конфигураций и др. Как оказалось, однофазные импульсы не приводят к подавлению. Это связано с формой потенциала действия, а точнее с фазой реполяризации, которая частично находится в отрицательной области. По этой же причине не удалось достичь подавления с помощью постоянного сигнала, подаваемого в среду.
Отдельно мы исследовали бифазный пилообразный импульс различных конфигураций: импульс, обе фазы которого имеют одинаковую амплитуду А, а также импульс, амплитуда отрицательной фазы которого была равна А/2. Однако подавления во втором случае добиться не удалось при всех значениях т. Для импульса с одинаковыми амплитудами фаз подавления удалось достичь только при т = 1. Поэтому далее мы будем рассматривать только такой пилообразный импульс.
Для поиска частот эффективного подавления мы воспользовались следующей простой методикой. В сравнительно небольшом объеме с нулевыми начальными условиями и пейсмекером собственной частоты генерировались кольцевые волны. Поскольку среда нелинейна, то возникающие кольцевые волны будут иметь иную частоту. Обозначим ее шои(. В фиксированной точке, находящейся на значительном удалении от ведущего центра, снимались значения переменной и во времени. Далее по полученному графику и{Ь) выполнялось преобразование Фурье
и находилась частота основной гармоники wout- Далее строилась зависимость ojout(u>in). Очевидно, что частоты — кандидаты на наиболее эффективное подавление — существуют в окрестности глобального максимума этой зависимости. Здесь мы исходим из известного принципа, что в распределенных средах «выживает» возбуждение наибольшей возможной частоты.
Все проводившиеся численные исследования были выполнены с ограничением на максимальное время наблюдения в 10000 мс. Иными словами, мы считали подавление достаточно успешным если оно укладывалось в этот период.
Рассмотрим среду, для которой управляющие параметры суть G\ = 1/30, G3 = 1.0. Тогда оптимальные частота и амплитуда внешнего би-фазного прямоугольного импульса, соответствующие данным параметрам среды составят u>in — 0.48 и А = 6, т.к. глобальные максимумы амплитудно-частотных зависимостей отвечают именно этим значениям параметров импульса. Однако попытка подавить турбулентность при этих параметрах привела к неудаче. Спиральные волны убывают во времени очень медленно, их число остается практически постоянным в течении 7000 мс. (рис. 46). В то же время, при u)in = 1.2 и А = 6 результат оказался куда более важным. Ведущий центр полностью вытесняет спирали за весьма короткий интервал t « 1700 мс. (рис. 4а).
Рис. 4 Зависимость числа фазовых сингулярностей N в системе от времени при: = 1.2 (а), и иы = 0.48 (б) А = 6,т = 0.7,в1 = 1/30, = 1.0.
Многочисленные скачки на этих графиках вызваны особенностью алгоритма подсчета сингулярностей. Принципиальной роли они не играют,
N
так как нас интересует непосредственно сам момент времени, когда этот график касается оси абцисс. Тщательное исследование позволяет сделать вывод, что к этому моменту в среде действительно не остается ни одной спиральной волны, так что система переходит в состояние покоя.
При тех же параметрах среды для пилообразного импульса с амплитудой А = 7 зависимость w0Ut(u/¿n) получилась весьма любопытной. Она представляет собой совокупность линейных участков, отвечающих различным режимам синхронизации среды и пейсмейкера (1:1,1:2,1:3 и т.д.).
Рис. 5 Зависимость числа фазовых сингулярностей N в системе, от времени для пилообразного импульса с частотой и:!п = 1.82 и амплитудой А -= 7 при (?! = 1/30, Сз = 1.0.
Самый высокий максимум частотной характеристики приходится в этом случае на значение = 1.82. Пейсмекер с этой частотой вытеснил из среды все спиральные волны примерно за 2500 мс. (см. рис. 5) , что на 50% дольше, чем для прямоугольного импульса. Тем не менее, это весьма удачный результат.
Однако далеко не при всех параметрах среды результаты получились столь же оптимистичными. Это обстоятельство связано со значительной чувствительностью среды к начатыми условиям, и об этом еще будет сказано ниже.
Рассмотрим теперь известную систему Фентона-Кармы, которая была изначально предложена как упрощенная ионная модель распростра-
нения потенциала действия в сердечной ткани. В наших исследованиях внешнее воздействие представляло собой однофазные или бифазные прямоугольные импульсы с амплитудой А и частотой uiin и подавалось в область размером 2x2 узла. Точно так же как и для модели ФицХью-Нагумо, для определения отношения длительностей полуволн бифазного импульса здесь использовался параметр т.
Мы считали, что период внешнего импульса Т < 500 ms, это условие было необходимо, чтобы обеспечить достаточно высокую частоту импульсов стабилизации. Также необходимо было найти значения параметров А ит, обеспечивающих наиболее эффективное управление динамикой системы. Численно было найдено, что оптимальными значениями для параметра т являются следующие: г = 0.05 -г- 0.15 для монофазных импульсов и т — 0.25 -=- 0.3 для бифазных.
Вначале исследовалась возможность подавить сложную активность монофазным воздействием. Было обнаружено, что хотя подавление наблюдается для частот 3,13 Гц и 7 Гц, оно сильно зависит от начальных условий. Например, монофазное воздействие при ш¿п = 7 Гц приводит к стабилизации хаотической динамики, если начинается в момент времени 500 мс, но оно неудачно, если начало подавления приходится на 600 мс. И наоборот, стимуляция частоты 3,13 Гц, приложенная в момент 600 мс, привела к подавлению спирально-волновой турбулентности.
Затем исследовались бифазные импульсы внешнего воздействия. В отличие от монофазной стимуляции, бифазное воздействие приводит к стабилизации сложной динамики стимулами с частотой 3,13 Гц, прикладываемыми в момент времени 500 мс, и стимуляцией с частотой 7,25 Гц, начинающейся в момент 600 мс. Частотный интервал, соответствующий второму максимуму частотной зависимости довольно широкий, поэтому достаточно трудно выбрать корректное значение частоты стимуляции для эффективного подавления хаотической динамики.
Следует отметить также, что стимуляция с частотой = 7 Гц, начавшаяся при t = 600 мс (приводящая к восстановлению турбулентной активности при А = 10/¿Л/2), обеспечивает эффективное подавление
при удвоении амплитуды стимуляции (рис. 6).Однако утроение амплитуды не дает положительного результата. Следовательно, зависимость эффективности подавления спирально-волновой активности от амплитуды нелинейна.
Рис. 6 Число фазовых сингулярностей как функция времени во время бифазной стимуляции с и,п = 7 Гц и амплитудами А = 10дЛ/2 (жирная линия), А = 20цА/2 (тонкая линия) и Л = 30/хЛ/2 (пунктирная линия).
Таким образом из изложенного во второй главе материала можно сделать один крайне важный вывод. Развитый пространственно-временной хаос в среде, описываемой системой (1), можно подавить слабым почти точечным периодическим воздействием. Этот вывод имеет принципиальное значение для приложений, когда требуется стабилизировать турбулентную динамику возбудимой среды. Кроме того, наши исследования показали, что эффективность (или, иными словами, скорость) подавления хаоса в системе, можно существенно повысить.
В третьей главе изучается силовое воздействие на все точки среды одновременно, обсуждаются преимущества и недостатки такого выведения среды из хаотического состояния. Рассматриваются дополнительно некоторые любопытные аспекты поведения системы, которые были обнаружены в процессе нашего анализа.
В данной диссертационной работе была проведена целая серия исследований с воздействием на всю среда'. Первый вполне естественный вывод, к которому мы пришли, состоит в том, что характер граничных условий качественного влияния на эффективность подавления не оказы-
50
О
1000 2000 3000 4000 5000 6000
вает. Поэтому в большинстве случаев мы ограничились периодическими граничными условиями. Воздействие подавалось в среду с развитым пространственно-временным хаосом в виде однократного прямоугольного импульса длительностью 1 мс.
При первых же численных экспериментах нам сразу удалось достичь положительного результата для параметров среды Сг = 0.01, вз = 0.5 и амплитуды импульса А = 1.3.
Вполне очевидно, что эффективность описанного подхода при заданных параметрах среды должна зависеть от амплитуды подаваемого импульса. Действительно, при определенном пороговом значении амплитуды импульса А, турбулентность может полностью восстановиться. Например, для параметров среды С\ = 0.01 и бз = 0.5 пороговое значение амплитуды А = 1.29. То есть при использовании импульса с амплитудой А = 1.29 турбулентная спирально-волновая динамика полностью подавляется, а если воздействовать импульсом амплитуды А = 1.28, то небольшой участок среды остается в возбужденном состоянии и становится источником спиральных волн, которые постепенно заполняют всю область.
В дополнение к этим исследованиям мы также рассмотрели систему с одним ведущим центром при тех же условиях. В этом случае наблюдались некоторые новые интересные эффекты. В частности оказалось, что даже при больших амплитудах общего внешнего воздействия возрождение хаоса в такой системе вполне возможно. После подачи внешнего воздействия в каждый элемент среды фронты кольцевых волн от ведущего центра начинают распадаться на рефрактерных участках, и постепенно спирально-волновая турбулентность заполняет всю среду целиком. Происходит это при достаточно больших частотах волн, исходящих от ведущего центра.
Еще одно интересное наблюдение было сделано, когда мы еще более увеличили время рефрактерности. Оказалось, что в среде с длительным периодом рефрактерности турбулентная динамика исчезает сама собой еще до подачи общего воздействия. Это обстоятельство дает ключ к по-
ниманию миогих процессов, происходящих в реальных возбудимых средах.
Рис. 7 Зависимость амплитуды импульса от параметра Сех
Как уже упоминалось ранее, управляющие параметры С\ и (?з системы (1) - (3) обратно пропорциональны периодам абсолютной и относительной рефрактериости соответственно. Поэтому весьма любопытно изучить как порог возбудимости среды зависит от этих параметров. Для этих целей мы ввели новый совокупный параметр Сех = -¡¡г + ¿г-, и для каждого значения этого параметра было найдено значение пороговой амплитуды общего внешнего воздействия, при которой происходило полное подавление хаотической динамики в среде (рис. 7).
Хорошо видно, что максимальная амплитуда использовавшегося импульса равна 1.5, т.е. она по крайней мере в 4 раза меньше амплитуды точечного воздействия (см. Главу 2). Тем не менее, как уже упоминалось выше, применение этой методики в эксперименте может оказаться весьма затруднительным.
Таким образом, с точки зрения эффективности, метод общего одновременного силового воздействия на каждый элемент среды в целом можно считать предпочтительным с точки зрения эффективности. При правильно подобранной амплитуде всегда удается избавиться от турбулентности. Однако необходимо помнить, что если в системе еще дополнительно имеется точечный источник возбуждения, как, например, это имеет место для сердечной ткани, он может стать причиной достаточно
быстрого возрождения пространственно-временной хаотической динамики.
В четвертой главе рассматривается среда с несколькими неподвижными и движущимися ведущими центрами, обсуждаются преимущества и недостатки этого подхода, проводятся аналогии с результатами, полученными в предыдущих главах, производится оптимизация параметров внешнего воздействия, что позволяет обеспечить выведение из хаоса при любых начальных условиях.
Представим ситуацию, когда ведущий центр окружен несколькими спиралями таким образом, что пространство вокруг него занято волновыми фронтами. В этом случае пейсмекер будет полностью подавлен хаотичной средой, и динамика системы будет развиваться таким образом, как будто внешнее воздействие на нее практически не оказывается.
Для решения описанной проблемы, в диссертационной работе предложено два в принципе независимых направления дальнейших исследований:
• увеличение числа ведущих центров;
• переход от неподвижных ведущих центров к движущимся.
Самый простой способ — первый, т.е. увеличение числа областей, к
которым подводится внешнее возбуждение. Поэтому мы также изучали поведение системы с несколькими (от 2 до 8) внешними источниками возбуждения.
Из проведенных исследований мы можем сделать довольно любопытный вывод. В ряде случаев увеличение числа пейсмекеров может ускорить подавление турбулентной динамики. Однако так происходит далеко не всегда. Более того, иногда получалось, что увеличение числа ведущих центров давало обратный эффект, то есть время подавление увеличивалось. Связано это, по-видимому, с возникновением конкуренции между пейсмекерами.
Следующий этап работы был связан с анализом более общего случая, когда пейсмекеры не являются стационарными областями в среде, а движутся там с некоторой скоростью.
Система оказалась крайне чувствительной как к величине оу0(, так и числу ведущих центров и их расположению в области среды. В качестве закона движения для ведущих центров мы взяли проекцию синуса на горизонтальную либо вертикальную ось £ = зт(шго(£), то есть периодическое движение с амплитудой £о и частотой ал.о(. Но даже в этом случае мы столкнулись с очень сложным поведением. Система оказалась крайне чувствительной как к величине шго(, так и числу ведущих центров и их расположению в области среды. В целом наши исследования показали, что зависимость эффективности подавления хаоса в среде от расположения и скорости движения ведущих центров существенно нелинейна.
После всех проведенных расчетов было решено ограничиться двумя ведущими центрами, расположенными на вертикальной прямой, проходящей через центр области и медленно осциллирующими относительно нее по горизонтали. Частота и амплитуда колебаний ведущих центров была фиксированной: £0 = 10, и/Г0( = 2.5 • Ю-4. Иными словами, пейсме-кер оказывался локализованным в малой области среды. Тем не менее такой модификации пейсмекера было вполне достаточно, чтобы в системе произошли заметные изменения по сравнению со случаем неподвижных ведущих центров.
Глобальная задача, которая ставилась на данном этапе, состояла в следующем. Необходимо было найти такие параметры для пейсмекера, чтобы остаться в рамках концепции слабого почти точечного внешнего воздействия, но при этом полностью гарантировать подавление хаотической динамики в среде независимо от того, в каком состоянии находилась система в начальный момент времени (и в особенности область, окружающая ведущий центр). Основные характеристики воздействия уже фиксированы, однако, остался еще один параметр, который мы обошли стороной - отношение длительности положительной полуволны к отрицательной т. Было обнаружено, что с точки зрения эффективности целесообразно брать г = 0.18. В дальнейшем мы везде будем пользоваться именно этим значением т.
Рис. 8 Число фазовых сингулярностей в системе (1) с движущимися внешними пейсмекерами как функция времени. = 0.01, (?з = 0.5, = 0.65, г = 0.18.
При всех указанных выше параметрах внешнего импульса нам удалось решить задачу успешного подавления при всех использовавшихся начальных условиях. В качестве примера приведем здесь один из полученных графиков динамики фазовых сингулярностей (рис. 8).
В пятой главе рассматриваются и обсуждаются возможные механизмы подавления спиральных волн в исследуемой среде.
В целом можно сделать следующие выводы. Всего механизмов подавления спиралей в системе два: вытеснение на периферию (для нулевых граничных условий) и взаимное уничтожение фазовых сингулярностей с противоположными хиральностями. Соответственно, для случая периодических границ остается один механизм. При этом ключом к успешному подавлению спиральных волн является именно неустойчивость. Связано это, видимо, с тем, что неустойчивость создает благоприятные условия для образования пар сингулярностей с противоположными хиральностями. Как известно, ядра волн с противоположными хиральностями при наличии дрейфа могут сталкиваться и взаимно уничтожаться. Эти общие наблюдения также представляют исследовательский интерес, хотя, безусловно, нуждаются в математической формализации. Но этот аспект следует рассматривать как продолжение данного диссертационного исследования.
В заключении обсуждаются возможные приложения полученных в работе результатов, а также выносимые на защиту положения.
1. На базе модели ФицХыо-Нагумо показано, что стабилизация хаотической динамики внешним точечным воздействием вполне осуществима. Предложенный подход позволяет полностью вытеснить все спирали и восстановить регулярное поведение системы.
2. Разработанный метод стабилизации также опробован на примере модели Фентона-Кармы, продемонстрирована его эффективность.
3. Выявлены недостатки такой стабилизации и показано, что для получения успешного результата предпочтительно использование движущегося пейсмекера (пейсмекеров).
4. Продемонстрировано, что спирально-волновая турбулентность эффективно подавляется силовым воздействием, приложенным ко всей области среды, причем эффективность нелинейно связана с амплитудой этого воздействия.
5. Обоснована возможность использования предложенной методики подавления спирально-волновой динамики в различных областях физики сложных систем.
Список публикаций по теме работы
1. А.Ю.Лоскутов, Р.В.Черемин, С.А. Высоцкий. Стабилизация турбулентной динамики возбудимых сред внешним точечным воздействием.— ДАН, том 404, N 4, с. 1-4, (2005).
2. А.Ю.Лоскутов, С.А.Высоцкий. Новый подход к проблеме дефибрилляции: подавление спирально-волновой активности сердечной ткани,— Письма в ЖЭТФ, том 84, вып. 9, с. 616-621, (2006).
3. E.Zhuchkova, B.Radnayev, S.Vysotsky, A.Loskutov. Suppression of turbulent dynamics in models of cardiac tissue by weak local excitations.— In: Understanding Complex Systems, ed. S.K. Dana, P.K. Roy, J. Kurths, Springer Berlin/Heidelberg, pp. 89-105, (2009).
4. S.A.Vysotskiy, R.V.Cheremin, A.Loskutov. Suppression of spatio-temporal chaos in simple models of re-entrant fibrillations.— J. of Physics: Conference Series, v. 23, pp. 202-209, (2005).
5. S.A.Vysotsky, R.V.Cheremin and A.Loskutov. Suppression of spiral-wave turbulence by point weak excitations.— Proc. of 2005 Int. Conf. «Physics and Control», August 24-26, Saint Petersburg, Russia — IEEE, p.236-239, (2005).
6. S.A.Vysotskiy, R.V.Cheremin, A.Loskutov. Suppression of spiral-wave turbulence by point weak excitations.— IEEE XPlore, ISBN: 0-7803-9235-3, pp. 236-239, (2005).
7. A.Loskutov, S.Vysotskiy, S.Boccaletti. New Methods of Suppression of the Spiral Wave Activity in Cardiac Tissue.— Proceedings of 17th International Workshop on Nonlinear Dynamics of Electronic Systems, pp. 153-156, (2009).
8. S.A.Vysotsky, R.V.Cheremin and A.Loskutov. Suppression of spatio-temporal chaos in simple models of re-entrant fibrillations.— Proc. of Int. Conf. on Control and Synchronization of Dynamical Systems, October 4-7, 2005, Leon GTo., Mexico, p.20-21.
9. S.Vysotskyi, R.Cheremin and A.Loskutov. Stabilization of fatal cardiac rhythms.— Book of abstracts of XXV Dynamics Days Europe 2005, Berlin, Germany, July 25-28, 2005, p.229.
10. A.Loskutov and S.Vysotsky. A new strategy of the defibrillation: Suppression of the spiral wave turbulence by moving pacemaker(s).— Abstract Collection of the 3rd International IEEE Scientific Conference on Physics and Control (PhysCon 2007), p. 317, (2007).
11. A.Loskutov and S.Vysotsky. New Methods of Suppression of the Spiral Wave Activity in Cardiac Tissue.— Indo-Russian workshop «Complete networks and applications». Kolkata, India, December 1-2, 2009.
12. С.А.Высоцкий, Р.В.Черемин, А.Ю.Лоскутов. Исследование динамики пространственно-распределенных систем. Приложения к кардиологии. Международная конференция студентов, аспирантов и молодых ученых по фундаментальным наукам «Ломоносов-2005», Россия, Москва, (2005).
Подписано в печать: 12.05.10
Объем: 1,5 усл.печ.л. Тираж: 100 экз. Заказ № 256 Отпечатано в типографии «Реглет» 119526, г.Москва, пр-т Вернадского, 39 (495) 363-78-90; www.reglet.ru
I Общие свойства нелинейных систем
1.1 Подавление хаоса в динамических системах.
1.1.1 Системы с внешними возмущениями.
1.1.2 Мультипликативное и аддитивное возмущения.
1.1.3 Общие свойства периодически возмущаемых систем.
1.2 Нелинейные среды с диффузией.
1.3 Модели распределенных систем.
1.3.1 Дискретные модели.
1.3.2 Континуальная модель возбудимой среды. Реакционно-диффузионные системы.
1.4 Распространение возбуждения в активных средах.
1.5 Спиральные волны в возбудимых средах.
I.5.1 Фазовые сингулярности.
II Стабилизация турбулентной динамики точечным возбуждением
II. 1 Реакционно-диффузионные модели возбудимых сред.
И.2 Сосредоточенная система.
11.3 Формулировка проблемы стабилизации турбулентной динамики [122, 171-173].
11.4 Кинематический подход.
11.5 Рождение пространственно-временного хаоса [122,168,171-173].
11.6 Методика исследования [95,122,132].
11.7 Внешнее воздействие [95,122,132].
II.7.1 Форма импульса.
И.7.2 Оптимальные частоты
II.7.3 Амплитуда импульса.
11.8 Стабилизация пространственно-временного хаоса [95,132,133,172]
II.9 Модель Фентона-Кармы [171]
III Силовая стабилизация спирально-волновой турбулентности
III. 1 Преимущества и недостатки метода силовой стабилизации
III. 2 Результаты и обсуждение.
111.2.1 Возрождение хаоса
111.2.2 Система с ведущим центром.
IV Новые подходы к проблеме стабилизации турбулентной динамики возбудимых сред
IV. 1 Новый подход к проблеме [171].
1У.2 Нестационарные ведущие центры [171,173].
IV.2.1 Чувствительность среды.
IV.2.2 Оптимизация выбора параметров воздействия.
V Механизмы стабилизации спирально-волновой динамики
V.1 Резонансный дрейф спиральных волн.
V.2 Устойчивая среда.
V.3 Неустойчивая среда.
Распределенные системы представляют собой достаточно широкий класс, включающий самые разные математические, физические, химические и биологические системы. Любой объект, если его компонента распределена в пространстве или времени по некоторому закону, является распределенным. К таким объектам относятся, например, жидкость или газ, популяция какого-либо вида, человеческий мозг или проводящая ткань сердца. Не удивительно, что уже очень давно данный класс систем привлекает исследователей из самых разных областей.
Однако настоящий прорыв в этой области произошел в конце 80-х — начале 90-х годов двадцатого столетия, когда стала бурно развиваться компьютерная техника. Это связано в первую очередь с тем, что из-за чрезвычайной сложности реалистичные модели распределенных систем почти не поддаются теоретическому исследованию. Ситуацию усложняет еще и тот факт, что подавляющее большинство таких моделей существенно нелинейны и при определенных условиях могут демонстрировать весьма сложное (хаотическое и квазипериодическое) поведение. Развитый математический аппарат нелинейной динамики оказывается зачастую неприменимым к распределенным системам по ряду причин, и поэтому необходимо разрабатывать новые методы исследования.
В последние годы развитие теории динамических систем и компьютерных методов моделирования позволило сформировать комплексный подход к исследованиям сложных распределенных систем разнообразной природы. Один из таких методов — представить исследуемую среду как совокупность автоколебательных или возбудимых элементов, локально взаимодействующих друг с другом. Тогда использование комплексного подхода предоставит новые возможности для более глубокого понимания процессов и явлений, наблюдаемых в активных средах. Явления, к которым ранее был возможен только эмпирический подход, можно таким образом достаточно подробно исследовать теоретически.
Наиболее сложная проблема теории распределенных систем -управление их динамикой и, в частности, стабилизация (подавление) сложных режимов поведения (квазипериодических или хаоса). Эта проблема возникла достаточно давно, и связана она с тем, что такие сложные режимы, как правило, являются крайне нежелательными. В особенности это касается реальных физических, экологических, химических или биологических систем.
Под стабилизацией неустойчивого или хаотического поведения динамических систем обычно понимается искусственное создание и поддержание в этих системах устойчивых (как правило, периодических) колебаний посредством внешнего воздействия. Эта задача, не смотря на простоту формулировки, оказывается весьма сложной научной проблемой, особенно в приложении к распределенным средам.
Актуальность такой задачи вполне очевидна. Рассмотрим несколько реальных примеров. В применении к сердечной ткани выведение системы на требуемый режим дает возможность влиять на сердечный ритм. Дело в том, что в настоящее время в теории возбудимых систем доминирует гипотеза, согласно которой возникновение фатальных сердечных аритмий — фибрилляций — есть следствие рождения в сердечной ткани большого количества автоволновых источников: спиральных волн или вихревых структур (т.е. пространственно-временного хаоса). Современные методы стабилизации таких режимов с помощью одиночных электрических импульсов (в том числе от имплантируемых дефибрилляторов) являются весьма жесткими и далеко не всегда приводят к успеху. Однако исследования самого последнего времени открывают новые возможности. Оказывается, что турбулентный режим во многих возбудимых средах может быть стабилизирован достаточно слабым периодическим параметрическим, или силовым воздействием, приложенным к некоторой области среды.
Для реакции Белоусова-Жаботинского такое воздействие позволяет создавать структуры нужного вида, получив таким образом систему, способную распознавать образы.
Помимо перечисленных здесь идей, результаты данной работы можно применить для исследования самых разнообразных реальных систем, описывающих полимеризацию, структурообразование, гранулярные среды, фракталы и т.н.
Цели работы:
Построение математической модели распределенной среды с различными граничными условиями, состоящей из автоколебательных элементов.
Исследование динамики системы в том числе со слабым почти точечным внешним воздействием.
Анализ поведения системы в зависимости от параметров возбудимости среды, а так же характеристик внешнего возбуждения.
Стабилизация сложных режимов поведения, связанных со спирально-волновой турбулентностью.
Научная новизна полученных в диссертации результатов состоит в следующем.
1. Изучена динамика математической модели распределенной среды с внешним воздействием при различных граничных и начальных условиях и в широком диапазоне параметров как среды так и внешнего воздействия.
2. На основе теории динамических систем и новых результатов в области теории распределенных сред показана возможность управления динамикой представленной модели.
3. Показана принципиальная возможность выведения системы из состояния сложной (в том числе хаотической) динамики слабым почти точечным воздействием.
4. Исследована возможность одновременного воздействия иа все элементы среды, указаны преимущества и недостатки такого подхода.
5. Предложен новый метод подавления спирально-волновой турбулентности с помощью движущихся ведущих центров, позволивший заметно повысить эффективность стабилизации. Получена зависимость эффективности подавления от количества ведущих центров и их характеристик.
Практическая ценность работы вполне понятна и заключается в следующем:
1. Показана принципиально новая возможность подавления спирально-волновой турбулентности возбудимой среды внешним точечным воздействием малой амплитуды и выведения ее на периодический режим движения.
2. Разработаны практические методы, позволяющие вычислять наиболее предпочтительные для такой стабилизации частоты.
3. Найдено, что нестационарный пейсмекер, расположенный в среде, значительно повышает эффективность предложенного метода.
Работа состоит из введения, общих сведений (обзора литературы), четырех глав, посвященных результатам, заключения и списка литературы.
В первой главе, посвященной общим вопросам, содержится обзор литературы, вводятся необходимые математические понятия, и рассматриваются различные модели распределенных сред, как дискретные, так и континуальные, обсуждаются аспекты, касающиеся численного моделирования.
В следующей главе проводится предварительное теоретическое исследование моделируемой системы уравнений, приводятся результаты численного анализа системы с одним неподвижным внешним источником возбуждения и изучаются возможные режимы движения среды в зависимости от управляющих параметров. Подробно описан переход системы к режиму спирально-волновой турбулентности. Здесь же показана принципиальная возможность подавления сложной динамики слабым почти точечным воздействием, а также ставятся новые вопросы.
В третьей главе изучается силовое воздействие на все точки среды одновременно, обсуждаются преимущества и недостатки такого выведения среды из хаотического состояния, рассматриваются дополнительно некоторые любопытные аспекты поведения системы, которые были обнаружены в процессе исследования.
В четвертой главе рассматривается среда с несколькими неподвижными и движущимися ведущими центрами, обсуждаются преимущества и недостатки этого подхода, проводятся аналогии с результатами, полученными в предыдущих главах, производится оптимизация параметров внешнего воздействия, что позволяет обеспечить выведение из хаоса при любых начальных условиях.
В пятой главе рассматриваются и обсуждаются возможные механизмы подавления спиральных волн в исследуемой среде. Всего таких механизмов существует два: уничтожение спиральной волны при столкновении с границей области и взаимное уничтожение двух спиральных волн, обладающих противоположными хиральностями.
В заключении обсуждаются возможные приложения полученных в работе результатов, а также выносимые на защиту положения.
1. А.Ю.Лоскутов. Динамический хаос. Системы классической механики.— Успехи физ. наук, т. 177, N0 9, с. 989-1015, (2007).
2. Л.Больцман. Статьи и речи.— М., Наука, (1970).
3. Н.С.Крылов. Работы по обоснованию статистической физики.— М.-Л., Изд-во АН СССР, (1950).
4. Я.Г.Синай. К обоснованию эргодической гипотезы для одной динамической системы статистической механики.— Докл. АН СССР, т. 153, N0 6, с. 1261-1264, (1963).
5. В.И.Арнольд, А.Авец. Эргодические проблемы классической механики.— Редакция журн. «Регулярная и хаот. динамика», Ижевск, (1999).
6. Ю.Мозер. Лекции о гамилыпоновых системах — М., Мир, (1973).
7. А.Лихтенберг, М.Либерман. Регулярная и стохастическая динамика— М., Мир, (1984).
8. А.Н.Колмогоров. Об" энтропии на единицу времени как метрическом инварианте автоморфизмов.— ДАН СССР, т. 124, N0 4, с. 754-755, (1959).
9. Я.Г.Синай. О понятии энтропии динамической системы — ДАН СССР, т. 124, N0 4, с. 768-771, (1959).
10. Н.Мартин, Дж.Ингленд. Математическая теория энтропии.— М., Мир, (1988).
11. С.Смейл. Дифференцируемые динамические системы.— Успехи матем. наук, т. 25, вып.1, с. 113-185, (1970).
12. З.Нитецки. Введение в дифференциальную динамику — М., Мир, (1975).
13. Д.В.Аносов. Геодезические потоки на замкнутых римановых многообразиях отрицательной кривизны.— М., Наука, (1967).
14. A.Lasota, M.C.Mackey. Chaos, Fractals and Noise. Stochastic Aspects of Dynamics.— Springer, Berlin, (1994).
15. А.Б.Каток, Б.Хассельблат. Введение в современную теорию динамических систем — М., Факториал, (1999).
16. Динамические системы. Тома 1-9. Серия: Современные проблемы математики. Фундаментальные направления.— ВИНИТИ, (19851991).
17. Я.Г.Синай. Современные проблемы эргодической теории.— М., Наука, (1995).
18. Ргос. of the SPIE 1993 Annual Meeting «Chaos in Communications».— San Diego, California, 11-16 July, v. 2038, (1993).
19. A. Loskutov, V.M.Tereshko. Extraction of the prototypes encoded in a chaotic attractor. In: Artificial Neural Networks, eds. I.Alexander and J. Taylor.- Elsevier, North-Holland, pp. 449-452, (1992).
20. A.Loskutov, V.M.Tereshko. Processing information encoded in chaotic sets of dynamical systems.— SPIE, v. 2038, pp. 263-272, (1993).
21. S.Hayes, C.Grebogi, E.Ott. Communicating with chaos.— Phys. Rev. Lett., v. 70, No 20, pp. 3031-3034, (1993).
22. S.Hayes, C.Grebogi, E.Ott, A.Mark. Experimental control of chaos for communication.— Phys. Rev. Lett., v. 73, No 13, pp. 1781-1784, (1994).
23. H.D.I.Abarbanel, P.S.Linsay. Secure communications and unstable periodic orbits of strange attractors.— IEEE Trans. Circuits Systs., v. 40, No 10, pp. 643-645, (1993).
24. А.Ю.Лоскутов, С.Д.Рыбалко, А.А.Чураев. Система кодирования информации посредством, стабилизации циклов динамических систем,.— Письма в ЖТФ, т.ЗО, вып. 20, с. 1-7, (2004).
25. A.Loskutov, G.E.Thomas. On a possible mechanism of self-organization in a two-dimensional network of coupled quadratic maps — SPIE, v. 2037, pp. 238-249, (1993).
26. L.Bresler, G.Metcalfe, J.M.Ottino, T.Shinbrot. Isolated mixing regions: origin, robustness and control— Chem. Eng. Sci., v. 58, pp. 1671-1679, (1996).
27. T.Shinbrot, J.M.Ottino. A geometric method to create coherent structures in chaotic flows — Phys. Rev. Lett., v. 71, pp. 843-847, (1993).
28. А.Ю.Лоскутов, Г.Э.Томас. Хаос и дестохастизация в двумерной решетке сцепленных отображений— Вестн. Моск. ун-та, сер. Физ.-астр., т. 34, No 5, с. 3-11, (1993).
29. L.Glass. Cardiac arrhythmias and circle maps — A classical problem.— Chaos, v. 1, No 1, pp. 13-19, (1991).
30. A.Garfinkel, M.L.Spano, W.L.Ditto. Controlling cardiac chaos— Science, v. 257, pp. 1230-1235, (1992).
31. А.Ю.Лоскутов, С.Д.Рыбалко. О динамике отображения окружности при параметрическом воздействии.— Вестн. Моск. ун-та, сср. Физ.-астр., т. 34, No 4, с. 19-27, (1993).
32. А.Ю.Лоскутов. Нелинейная динамика и сердечная аритмия.— Прикладная нелинейная динамика, т. 2, No 3-4, с. 14-25, (1994).
33. A.Loskutov, V.M.Tereshko, K.A.Vasiliev. Predicted dynamics for cyclic cascades of chaotic deterministic automata.— Int. J. Neural Systems, v. 6, pp. 175-182, (1995).
34. R.V.Sole, L.Menendez de 1a. Prida. Controlling chaos in discrete neural networks.- Phys. Lett. A, v. 199, No 1-2, pp. 65-69, (1995).
35. S.J.Schiff, K.Jerger, D.H.Duong, T.Chang, M.L.Spano, W.L.Ditto. Controlling chaos in the brain.— Nature, v. 370, pp. 615-620, (1994).
36. G.Chen, X.Dong. From chaos to order — Perspectives and methodologies in controlling chaotic nonlinear dynamical systems— Int. J. Bifurcation and Chaos, v. 3, No 6, pp. 1363-1409, (1993).
37. T.Shinbrot. Progress in the control of chaos — Adv. Phys., v. 44, No 2, pp. 73-111, (1995).
38. R.Lima, M.Pettini. Suppression of chaos by resonant parametric perturbations.- Phys. Rev. A, v. 41, No 2, pp. 726-733, (1990).
39. J.Singer, Y-Z.Wang, H.H.Bau. Controlling a chaotic system— Phys. Rev. Lett,., v. 66, pp. 1123-1125, (1991).
40. L.Fronzoni, M.Geocondo, M.Pettini. Experimental evidence of suppression of chaos by resonant parametric perturbations — Phys. Rev. A, v. 43, pp. 6483-6487, (1991).
41. Y.Braiman, I.Goldhirsh. Taming chaotic dynamics with weak periodic perturbationsPhys. Rev. Lett., v. 66, pp. 2545-2548, (1991).
42. S.Rajasekar, M.Lakshmanan. Algorithms for controlling chaotic motion: application for the BVP oscillator.— Physica D, v. 67, No 1-3, pp. 282-300, (1993).
43. S.Bielawski, D.Derozier, P.Glorieux. Controlling unstable periodic orbits by a delayed continuous feedback.— Phys. Rev. E, v. 49, No 2, pp. 971-974, (1994).
44. Ph.V.Bayly, L.N.Virgin. Practical considerations in the control of chaos.— Phys. Rev. E, v. 50, Nol, pp. 604-607, (1994).
45. D.Vassiliadis. Parametric adaptive control and parameter identification of low-dimensional chaotic systems.— Physica D, v. 71, No 1-2, pp. 319-341, (1994).
46. B.Hiibinger, R.Doerner, W.Martienssen. Controlling chaos experimentally in systems exhibiting large effective Lyapunov exponents.— Phys.-Rev. E, v. 50, No 2, pp. 932-948, (1994).
47. R.Mettini, T.Kurz. Optimized periodic control of chaotic systems.— Phys. Lett. A, v. 206, No 5-6, pp. 331-339, (1995).
48. R.Chacon. Suppression of chaos by selective resonant parametric perturbations— Phys. Rev. E, v. 51, No 1, pp. 761-764, (1995).
49. T.Shinbrot. Chaos: Unpredictable Yet Controllable?— Nonlinear Sci. Today, v. 3, No 2, pp. 1-8, (1993).
50. T.Shinbrot, C.Grebogi, E.Ott, J.A.Yorke. Using small perturbations to control chaos — Nature, v. 363, pp. 411-417, (1993).
51. M.Pettini. Controlling chaos through parametric excitations. In: Dynamics and Stochastic Processes. Ed. R.Lima, L.Streit, R.Vilela Mendes — Springer, Berlin, pp. 242-250, (1990).
52. Yu.S.Kivshar, B.Rodelsperger, H.Benner. Suppression of chaos by non-resonant parametric perturbations— Phys. Rev. E, v. 49, pp. 319-324, (1994).
53. A.B.Corbet. Suppression of chaos in ID maps — Phys. Lett. A, v. 130, No 4-5, pp. 267-270, (1988).
54. E. A. Jackson, A. Hübler. Periodic entrainment of chaotic logistic map dynamics.— Physica D, v. 44, pp. 407-420, (1990).
55. E.A.Jackson. Control of dynamics flows with attractors — Phys. Rev. A, v. 44, pp. 4839-4853, (1991).
56. B.A.Huberman, E.Lumer. Dynamics of adaptive systems— IEEE Trans. Circ. Syst., v. 37, pp. 547-550, (1990).
57. K.Pyragas. Stabilization of unstable periodic and aperiodic orbits of chaotic systems by s elf-controlling feedback.-— Z. Naturforsch A, v. 48, pp. 629-632; (1993).
58. R.Chacon. Geometrical resonance as a chaos eliminating mechanism.— Phys. Rev. Lett., v. 77, pp. 482-485, (1995).
59. R.Mettin. Control of chaotic maps by optimized periodic inputs.— Int. J. Bifurcation and Chaos, v. 8, No 8, pp. 1707-1711, (1998).
60. A. Hübler, R. Georgii, M. Kuckler, W. Stelzl, E. Lüsher. Resonant stimulation of nonlinear damped oscillators by Poincaré maps.— Helv. Phys. Acta, v. 61, pp. 897-900, (1988).
61. J.F.Linder, W.L.Ditto. Removal, suppression, and control of chaos by nonlinear design — Appl. Mech. Rev., v. 48, No 12, pp. 795-807, (1995).
62. E.Ott, M.L.Spano. Controlling chaos — Physics Today, v. 48, No 5, pp. 34-40, (1995).
63. E.Ott, C.Grebogi, J.A.Yorke. Controlling chaos — Phys. Rev. Lett., v. 64, pp. 1196-1199, (1990).
64. F.J.Romeiras, E.Ott, C.Grebogi, W.P.Dayawansa. Controlling chaotic dynamical systems — Physica D, v. 58, pp. 165-192, (1992).
65. А.Ю.Лоскутов, А.И.Шишмарев. Об одном свойстве семейства квадратичных отображений при параметрическом воздействии — Успехи матем. наук, т. 48, вып. 1, с. 169-170, (1993).
66. A.Loskutov, A.I.Shishmarev. Control of dynamical systems behavior by parametric perturbations: an analytic approach.— Chaos, v. 4, No 2, pp. 351-355, (1994).
67. A.Loskutov. Non-feedback controlling complex behaviour: an analytic approach — In: Nonlinear Dynamics: New Theoretical and Applied Results. Ed. J.Awreicewicz.— Springer, Berlin, pp. 125-150, (1995).
68. A.N.Deryugin, A.Loskutov, V.M.Tereshko. Inducing stable periodic behaviour in a class of dynamical systems by parametric perturbations.— Chaos, Solitons & Fractals, v.7, NolO, pp. 1555-1567, (1996).
69. A.Loskutov, V.M.Tereshko, K.A.Vasiliev. Stabilization of chaotic dynamics of one-dimensional maps by a cyclic parametric transformation.— Int. J. Bif. and Chaos, v. 6, No 4, pp. 725-735, (1996).
70. А.Ю.Лоскутов. Хаос и управление динамическими системами.— В сб. Нелинейная динамика и управление. Вып. 1. Ред. С.В.Емельянов, С.К.Коровин. Москва, Физматлит, с. 163-216, (2001).
71. А.Н.Дерюгин, А.Ю.Лоскутов, В.М.Терешко. К вопросу о рождении устойчивого периодического поведения параметрически возбуждаемых динамических систем,— ТМФ, т. 104, No 3, с. 507512, (1995).
72. Z.Galias. New method for stabilization of unstable periodic orbits in chaotic systems.— Int. J. Bif. and Chaos, v. 5, No 1, pp. 281-295, (1995).
73. T.Shinbrot, E.Ott, C.Grebogi, J.A.Yorke. Using chaos to direct trajectories to targetsPhys. Rev. Lett., v. 65, pp. 3215-3218, (1990).
74. T.Shinbrot, C.Grebogi, E.Ott, J.A.Yorke. Using chaos to target stationary states of flows.— Phys. Lett. A, v. 169, pp. 349-354, (1992).
75. T.Shinbrot, E.Ott, C.Grebogi, J.A.Yorke. Using chaos to direct orbits to targets in systems describable by a one-dimensional map.— Phys. Rev. A, 1992, v. 45, No 6, pp. 4165-4168.
76. E.Kostelich, C.Grebogi, E.Ott, J.A.Yorke. Higher dimensional targetting- Phys. Rev. E, v. 47, pp. 305-310, (1993).
77. R.Meucci, W.Gadomski, M.Ciofini, F.T.Arecchi. Experimental control of chaos by weak parametric perturbations.— Phys. Rev. E, v. 49, No 4, pp. 2528-2531, (1994).
78. Ю.И.Неймарк. Динамические системы и управляемые процессы — М., Наука, (1978).
79. М.Розо. Нелинейные колебания и теория устойчивости.— М., Наука, (1971).
80. Е.Н.Дудник, Ю.И.Кузнецов, И.И.Минакова., Ю.М.Романовский. Синхронизация в системах со странным аттрактором.— Вестн. МГУ, сер. Физ.-астр., т. 38, No 4, с. 84-87, (1983).
81. Ю.И.Кузнецов, В.В.Милюлин, И.И.Минакова, Б.А.Сильнов. Синхронизация хаотических автоколебаний.— Докл. АН СССР, т.275, No 4-6, с. 1388-1391, (1984).
82. Ю.И.Кузнецов, П.С.Ланда, А.Ф.Ольховой, С.М.Перминов. Связь между амплитудным порогом синхронизации и энтропией в стохастических автоколебательных системах.— Докл. АН СССР, т. 281, No 2, с. 291-294, (1985).
83. Ю.И.Неймарк, П.С.Ланда. Стохастические и хаотические колебания.— М., Наука, (1987).
84. В.В.Алексеев, А.Ю.Лоскутов. Дестохастизация системы со странным аттрактором посредством, параметрического воздействия.— Вестник Моск. ун-та, сер. Физ.-астр., т. 26, No 3, с. 40-44, (1985).
85. В.В.Алексеев, А.Ю.Лоскутов. О возможности управления системой со странным аттрактором,— В сб. Проблемы экологического мониторинга и моделирования экосистем. Том VIII.— Ленинград, Гидрометеоизда.т, с. 175-189, (1985).
86. В.В.Алексеев, А.Ю.Лоскутов. Управление системой со странным аттрактором■ посредством периодического параметрического воздействия.—Л АН СССР, т. 293, вып. 6, с. 1346-1348, (1987).
87. P.So, E.Ott. Controlling chaos using time delay coordinates via stabilization of periodic orbits — Phys. Rev. E, 1995, v. 51, No -4, pp. 2955-2962.
88. G.A.Johnson, M.Locher, E.R.Hunt. Stabilized spatiotemporal waves in a convectively unstable open flow system: coupled diode resonators — Phys. Rev. E, v. 51, pp. 1625-1628, (1995).
89. D.Auerbach. Controlling extended systems of chaotic elements — Phys. Rev. Lett., v. 72, No 8, pp. 1184-1187, (1994).
90. M.Ding, E.Ott, C.Grebogi. Controlling chaos in a temporally irregular environment— Physica D, v.74, No 1-2, pp. 386-394, (1994).
91. J.E.S.Socolar, D.W.Sukow, D.J.Gauthier. Stabilizing unstable periodic orbits in fast dynamical systems — Phys. Rev E, v. 50, No 4, pp. 32453248, (1994).
92. A.Kittel, K.Pyragas, R.Richter. Prerecorded history of a system as an experimental tool to control chaos.— Phys. Rev. E, v. 50, Nol, pp. 262-268, (1994).
93. M.A.Matias, J.Guemez. Stabilization of chaos by proportional pulses in the system variables.— Phys. Rev. Lett., v. 72, Nol, pp. 1455-1458, (1994).
94. H.Zhang, Zh.Cao, N.-J.Wu, H.-P.Ying, G.Hu. Suppress Winfree Turbulence by Local Forcing Excitable Systems— Phys. Rev. Lett., v. 94, pp. 188301, (2005).
95. А.Ю.Лоскутов, Р.В.Черемин, С.А. Высоцкий. Стабилизация '. турбулентной динамики возбудимых сред внешним точечнымвоздействием — ДАН, том 404, N 4, с. 1-4, (2005).
96. G.Yuan, G.Wang , S.Chen. Control of spiral waves and spatiotemporal chaos by periodic perturbation near the boundary.— Europhys. Lett., v. 72, pp. 908-914, (2005).
97. G.Tang, M.Deng, B.Hu, G.Hu. Active and passive control of spiral turbulence in excitable media.— Phys. Rev., v. 77, pp. 046217, (2008).
98. Zh.Cao, P.Li, H.Zhang, F.Xie, G.Hu Turbulence control with local pacing and its implication in cardiac defibrillation.— CHAOS, v. 17, pp. 015107, (2007).
99. T.K.Shajahan, A.R.Nayak, R.Pandit Spiral- Wave Turbulence and Its Control in the Presence of Inhomogeneities in Four Mathematical Models of Cardiac Tissue.- PLoS ONE, v. 4, pp. e4738, (2009).
100. A.Loskutov. Dynamics control of chaotic systems by parametric destochastization — J. Phys. A, v. 26, Nol8, pp. 4581-4594, (1993).
101. A.Loskutov, S.D.Rybalko. Parametric perturbations and suppression of chaos in n-dimensional maps — Preprint ICTP IC/94/347, Trieste, Italy, November, (1994).
102. Е.Ф.Мищенко, В.А.Садовничий, А.Ю.Колееов, Н.Х.Розов. Автоволновые процессы в нелинейных средах с диффузией — Физматлит, Москва, 2005.
103. Г.Николис, И.Пригожин. Самоорганизация в неравновесных системах. От диссипативных структур к упорядоченности через флуктуации.— Мир, Москва, 1979.
104. N.Wiener, A.Rosenblueth. The mathematical formulation of the problem of conduction of impulses in a network of connected excitable elements, specifically in cardiac muscle.— Arch. Inst. Cardiologiade, Mexico, v. 16, pp. 205-265, (1946).
105. Т.С.Ахромеева, С.П.Курдюмов, Г.Г.Малинецкий, А.А.Самарский. Нестационарные структуры и диффузионный хаос.— М.: Наука, (1990).
106. Г.Г.Малинецкий, А.Б.Потапов. Современные проблемы нелинейной динамики.— (2002).
107. K.Nam, E.Ott, M.Gabbay, P.N.Guzdar. Spiral wave dynamics in the complex Ginzburg-Landau equation with broken chiral symmetry.— Physica D, v. 118, pp. 69-83, (1998).
108. R.FitzHugh. Impulses and Physiological States in Theoretical Models of Nerve Membrane.— Biophysical Journal, v. 1, (1961).
109. J.Nagumo, S.Arimoto, S.Yoshizawa. An active pulse transmission line simulating nerve axon— Proc IRE. 50: 2061-2070, (1964).
110. F.Fenton, A.Karma. Vortex dynamics in three-dimensional continuous myocardium with fiber rotation: Filament instability and fibrillation.— Chaos, v. 8, No. 1, pp. 20-47, (1997).
111. C.-H.Luo, Y.Rudy. A Dynamic Model of the Cardiac Ventricular Action Potential: I. Simulations of Ionic Currents and Concentration Changes.- Circulation Research 74, 1071-1096, (1994).
112. G.W.Beeler, H.Reuter. Reconstruction of the action potential of ventricular myocardial fibres.— J. Physiol., V.268. pp.177-210, (1977).
113. A. Huxley, A. Hodgkin. Movement of Radioactive Potassium and Membrane Current in a Giant Axon— Journal of Physiology (London), (1953).
114. H.Zhang, B.Hu, G.Hu, Q.Ouyang, J.Kurths. Turbulence control by developing a spiral wave with a periodic signal injection in the complex Ginzburg-Landau equation — Phys. Rev. E, v. 66, pp. 046303, (2002).
115. A.Panfilov, P.Hogeweg. Spiral breakup in a modified FitzHugh-Nagumo model— Phys. Lett. A, v. 176, pp. 295—299, (1993).
116. Р.В.Черемин. Стабилизация хаотических процессов в возбудимых средах.— Дипломная работа, Физический ф-т МГУ, (2004).
117. D.Allexandre, N.F.Otani, Preventing alternans-induced spiral wave breakup in cardiac tissue: An ion-channel-based approach — Phys. Rev. E, v. 70, pp. 061903 (2004).
118. N.H.Packard, S.Wolfram. Two-dimensional cellular automata— J. Stat. Phys., v. 38. pp. 901-946, (1985).
119. T.Hogg, B.A.Huberman. Artificial intelligence and large scale computation: a physics perspective — Phys. Rep., v. 156. pp. 227-310, (1987).
120. В.И.Кринский, К.И.Агладзе. Вихри с топологическим■ зарядом, 2, 3, 4 в активной химической среде.— ДАН СССР., т. 263. с. 335-337, (1982).
121. K.Agladze, Matthew W.Kay, V.Krinsky, N.Sarvazyan. Interaction between spiral and paced waves in cardiac tissue — Am. J. Physiol. Heart Circ. Physiol. 293: H503-H513, (2007).
122. С.А.Высоцкий. Анализ пространственно-временной динамики распределенных реакционно-диффузионных систем.— Дипломная работа, МГУ, (2005).
123. P.S.Hagan. Spiral waves in reaction-diffusion equations— SIAP v.41, No. 4, pp. 762-786, (1982).
124. S.Koga. Rotating spiral waves in reaction-diffusion systems — Progress of Theoretical Physics, v.67 No.l pp. 164-178, (1981).
125. В.Н.Ардашев, А.В.Ардашев, В.И. Стеклов. Лечение нарушений сердечного ритма — Медпрактика-М, Москва, (2005).
126. A.Pumir, V.I.Krinsky. How does an electric field defibrillate cardiac muscle?.— Physica D, v. 91, pp. 205-219, (1996).
127. J.P.Keener, A.V.Panfilov. A biophysical model for defibrillation of cardiac tissue.— Biophys. J., v. 71, pp. 1335-1345, (1996).
128. J.P.Keener, E.Cytrynbaum. The effect of spatial scale of resistive in-homogeneity on defibrillation of cardiac tissue — J. Theor Biol., v. 223, pp. 233-248, (2003).
129. M.G.Fishier. Syncytial heterogeneity as a mechanism underlying cardiac far-field stimulatioln during defibrillation-level shocks.— J. Cardiovascular Electrophysiology, v. 9, No., 4, pp. 384-394, (1998).
130. N.Trayanova, K.Skouibine, F.Aguel. The role of cardiac tissue structure in defibrillation— Chaos, v. 8, No.l, pp. 221-233, (1998).
131. A.Garfinkel, et. al. Quasiperiodicity and Chaos in Cardiac Fibrillation— Clin. Invest., v. 99, No. 2, pp. 305-314, (1997).
132. S.A.Vysotskiy, R.V.Cheremin, A.Loskutov. Suppression of spatiotemporal chaos in simple models of re-entrant fibrillations.— J. of Physics: Conference Series, v. 23, pp. 202-209, (2005).
133. S.A.Vysotskiy, R.V.Cheremin, A.Loskutov. Suppression of spiral-wave turbulence by point weak excitations.— IEEE XPlore, ISBN: 0-78039235-3, pp. 236- 239, (2005).
134. S.Takagi, A.Pumir, D.Pazo, I.Efimov, V.Nikolski, V.Krinsky, Unpinning and Removal of a Rotating Wave in Cardiac Muscle — Phys. Rev. Lett., v. 93, No. 5, pp. 058101, (2004).
135. M.-A.Bray, et. al. Experimental and Theoretical Analysis of Phase Singularity Dynamics in Cardiac Tissue.— Journal of Cardiovascular Elec-trophysilology, v. 12, No. 6, (2001).
136. M.-A.Bray, J.P.Wikswo. Use of Topological Charge to Determine Filament Location and Dynamics in a Numerical Model of Scroll Wave Activity — IEEE Transactions on Biomediaca.1 Engineering, v. 49, No. 10, (2002).
137. V.I.Krinsky, A.Pumir. Models of defibrillation of cardiac tissue.— Chaos, v. 8, No. 1, pp. 188-203, (1998).
138. R.Zou, J.Kneller, L.J.Leon, S.Nattel. Development of a computer algorithm for the detection of phase singularities and initial application to analyze simulations of atrial fibrillation.— Chaos, v. 12, No. 3, (2002).
139. В.А.Давыдов, В.С.Зыков, А.С.Михайлов, П.К.Бражник. Дрейф и резонанс спиральных волн в активных средах.— Изв. вузов сер. «Радиофизика», т. 31, с. 574-582, (1988).
140. К.И.Агладзе, В.А.Давыдов, А.С.Михайлов. Наблюдение резонанса спиральных волн в возбудимой среде.— Письма ЖЭТФ, т. 45, стр. 601-603, (1987).
141. S.Schmidt, P.J.Ortoleva. A new chemical wave equation for ionic systems— J. Chem. Phys., v 67, pp. 3771-3776, (1977).
142. H.Sevcikova, M.Marek. Chemical waves in electric field.— Physica D, v. 9, pp. 140-156, (1983).
143. Я.Б.Зельдович. Нелинейные волны. Распространение и взаимодействие — Москва, Наука, 1981.
144. В.С.Зыков, А.А.Петров. О роли неоднородности возбудимой среды в механизмах самоподдерживающейся активности.— Биофизика, т. 22, стр. 300-306, (1977).
145. В.С.Зыков, О.Л.Морозова. Скорость распространения возбуждения в двухмерной возбудимой среде— Биофизика, т. 24, стр. 717-722, (1979).
146. В.С.Зыков. Аналитическая оценка зависимости скорости волны возбуждения в двухм.ерной возбудимой среде от кривизны ее фронта.— Биофизика, т. 25, стр. 888-895, (1980).
147. P.Foerster, S.C. Viiller, В.Hess. Curvature and propagation velocity of chemical waves — Science, v. 241, pp. 685-687, (1988).
148. В.А.Давыдов, В.С.Зыков, А.С.Михайлов. Кинематика автоволновых структур в возбудимых средах— УФН, т. 161, стр. 45-85, (1991).
149. A.S.Mikhailov, V.A.Davydov, V.S.Zykov, Complex dynamics of spiral waves and motion of curves.— Physica D, v. 70, pp. 1-39, (1994).
150. M. Pascual Diffusion induced chaos in a spatial predator-prey system — Proc. R. Soc. Lond. B, 251, 1-7, (1993).
151. V. Rai, G. Jayaraman Is diffusion-induced chaos robust?.— Current Science, v. 84, N. 7, pp. 925-929, (2003).
152. И.М.Гельфанд, И.Л.Цетлин. О континуальных моделях управляющих систем— ДАН СССР, т. 131, стр. 1242-1245, (I960).
153. И.С.Балаховский. Некоторые режимы движения в идеальной возбудимой среде.— Биофизика, т. 10, стр. 1063-1067, (1965).
154. В.И.Кринский. Фибрилляции в возбудимых средах,— в сб. Проблемы кибернетики, N 20, Москва, Наука, стр. 59-80, (1968).
155. P.C.Fife. Propagator-conroller systems and chemical patterns— in Non-equilibrium dynamics in chemical systems, editor C.Vidal, A.Pacault, Springer-Verlag,Berlin, pp. 76-88, (1984).
156. П.К.Рашевский. Курс дифференциальной геометрии,— Москва, Физматгиз, (1956).
157. В.С.Зыков. Моделирование волновых процессов в возбудимых средах,— Москва, Наука, (1984).
158. В.С.Зыков. Исследование на ГВС-100 некоторых характеристик самоподдерживающейся активности в возбудимой среде,— В сб. Управление сложными системами, Москва, Наука, стр. 59-62, (1975).
159. В.А.Давыдов, А.С.Михайлов. Спиральные волны в распределенных активных средах.— в сб. Нелинейные волны. Структуры и бифуркации, Москва, Наука, стр. 261-279, (1987).
160. П.К.Бражник, В.А.Давыдов, В.С.Зыков, А.С.Михайлов. Вихревые кольца в распределенных возбудимых средах.— ЖЭТФ, т. 93, стр. 1725-1736, (1987).
161. D.P.Zipes, J.Jalife. Cardiac Electrophysiology from Cell to BedSide.— 4th ed. Philadelphia: W.B. Saunders, (2004).
162. Е.А.Жучкова. Исследование активных сред методами теории динамических систем,— Канд. дис., МГУ, (2006).
163. F.Fenton, A.Karma. Vortex dynamics in three-dimensional continuous myocardium with fiber rotation: Filament instability and fibrillation.— Chaos, v.8, pp. 20-47, (1997).
164. A.T.Stamp, G.V.Osipov, J. J.Collins. Suppressing arrhythmias in cardiac models using overdrive pacing and calcium channel blockers — Chaos, v. 12, pp. 931-940, (2002).
165. А.Ю.Лоскутов, А.С.Михайлов. Основы теории сложных систем.— М.-Ижевск: Институт компьютерных исследований,
166. П.К.Бражник, В.А.Давыдов, В.С.Зыков, А.С.Михайлов. Дрейф и резонанс спиральных волн в активных средах— Изв. вузов. Сер. Радиофизика, т. 31, № 5, с. 574-584, (1988).
167. К.И.Агладзе, В.А.Давыдов, А.С.Михайлов. Наблюдение резол,анса спиральных волн в возбудимой распределенной среде.— Письма ЖЭТФ, т. 45, с. 601-603, (1987).
168. A.S.Mikhailov, A.Loskutov. Foundation of Synergetics. Chaos and Noise.— Springer, Berlin, (1996).
169. F.Fenton, E.M.Cherry, H.M.Hastings, S.J.Evans. Multiple mechanisms of spiral wave breakup in a model of cardiac electrical activity. — Chaos, v. 12, pp. 852-892, (2002).
170. Z.Cao, P.Li, H.Zhang, F.Xie, G.Hu. Turbulence control with local pacing and its implication in cardiac defibrillation.— CHAOS, v. 17, pp. 015107, (2007).
171. А.Ю.Лоскутов, С.А.Высоцкий. Новый подход к проблеме дефибрилляции: подавление спирально-волновой активности сердечной ткани — Письма в ЖЭТФ, том 84, вып. 9, с. 616-621,
172. A.Loskutov, S.Vysotskiy, S.Boccaletti. New Methods of Suppression of the Spiral Wave Activity in Cardiac Tissue.— Proceedings of 17th International Workshop on Nonlinear Dynamics of Electronic Systems, pp. 153-156, (2009).2007).2006).