Научная статья на тему 'Мезоморфная модель электрической активности кортекса'

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

CC BY
166
37
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
НЕЙРОН / СИСТЕМА ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ / NEURON / SYSTEM OF DIFFERENTIAL EQUATIONS

Аннотация научной статьи по математике, автор научной работы — Старченко И. Б., Борисова О. С.

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

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Mesomorphic Model of Cortex Electrical Activity

One of the mathematical models of mesoscopic cortex electrical activity was described. The system of 8 stochastic partial differential equations and 20 variables in dimemsionless form was used.

Текст научной работы на тему «Мезоморфная модель электрической активности кортекса»

воздействии.

Выводы

1. Показано, что появление звука при СКЭНАР-воздействии происходит по принципу работы электростатического громкоговорители. Роль гибкой полимерной пленки (мембраны) выполняет роговой (роговой+ блестящий) слой эпидермиса.

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

3. Выдвинуто предположение, что воздействие СКЭНАРа способствует изменению подкожной жидкостной структуры.

4. Рассмотрена гипотеза о дополнительной передаче информации по жидкостным каналам, включая информационные изменения, происходящие в жидкостных структурах при СКЭНАР-воздействии.

БИБЛИОГРАФИЧЕСКИЙ СПИСОК

1. Гринберг. Я.З. Об одном эффекте СКЭНАР-воздействия. Известия ТРТУ. Тематический выпуск. Материалы научно - технической конференции Медицинские информационные системы - МИС-2004». - Таганрог. - 2004, -№ 6(41), - С.100-105.

2. Гринберг Я.З. СКЭНАР-терапия и СКЭНАР-экспертиза. Некоторые аспекты. Рефлексология. - 2005. - №3(7). - С. 5-10.

3. Гринберг Я.З. СКЭНАР: построение, физические механизмы, основы эффективности. Нелекарственная медицина. - 2006. - №3(4). - С. 37-42.

4. Гринберг Я.З. Ещё раз об особенностях СКЭНАР-воздействия. Известия ТРТУ. Тематический выпуск. Материалы научно-технической конференции «Медицинские информационные системы - МИС -2006». - Таганрог. - 2006, - № 11(46), - С. 144-147.

5. БоголюбовВ.М., Пономаренко Г.Н. Общая физиотерапия. - М.: Медицина, 2003. - 432 с.

6. Финкельштейн А.В. Введение в физику белка: Курс лекций. 1999-2000 г. Цитируется по публикации в Интернете.

7. Зилов В.Г., Судаков К.В., Эпштейн О.И. Элементы информационной биологии и медицины: Монография. - М.: МГУЛ, 2000. - 248с.

8. Прохорцев И.В., Зимин Э.В. Смирнов А.П. Способ получения биоактивной воды. №2320545 от 03.27.2008.

9. Гринберг Я.З. Чрескожная электростимуляция: подход с позиции функционального континуума регуляторных пептидов //Рефлексотерапия. - 2002. - №1(1). - С. 29-32.

10. Зенин С.В. Исследование структуры воды методом протонного магнитного резонанса. // Докл. Акад. Наук, 332(3). 328-329 (1993).

11. Зенин С.В., Тяглов Б.В. Гидрофобная модель структуры ассоциатов молекул воды // Физ. Хим. - 1994. - № 68(4). - С. 636-641.

12. Розин И.Т. Исследования процессов энергоинформационного взаимодействия методами ИК спектроскопии // Сознание и физическая реальность. - 1999. - Т.4, - №2, - С. 60-64.

УДК 616.8-072.7(035)

И.Б. Старченко, О.С. Борисова

МЕЗОМОРФНАЯ МОДЕЛЬ ЭЛЕКТРИЧЕСКОЙ АКТИВНОСТИ

КОРТЕКСА

Для исследования электрической активности мозга человека были разработаны многочисленные модели [1, 2]. Существует математическая модель для описания поведения одиночного нейрона и сети нейронов. Чтобы смоделировать данные, полученные с одного

электрода, необходимо воспроизвести поведение приблизительно 105 отдельных нейронов. К сожалению, используя физиологически точные математические модели, эту задачу нельзя выполнить. Поэтому была разработана математическая модель мезоскопической электрической активности кортекса. Мезоскопические модели в отличие от модели поведения единичного нейрона Ходжикина-Хаксли [1] аппроксимируют усредненную активность, или усредненное поле, совокупности клеток. Для развития этих моделей были определены выражения для усредненных свойств соседних нейронов. Полученные переменные описывают среднее значение соматического мембранного потенциала совокупности клеток или среднее значение подкорковой активности, полученные в объеме кортекса. В этих моделях наименьшая единица активности - совокупность клеток, а не один нейрон.

Мезоскопическая модель описывает электрическую активность столбца коры, а не отдельного нейрона. Имеется, по крайней мере, две причины для использования мезоскопической модели для обработки данных ЭКоГ человека. Во-первых, электрокортиограф записывает суммарную электрическую активность миллионов отдельных нейронов. Поэтому мезоскопическая модель и записи ЭКоГ воспроизводят моделируемые и реальные результаты, соответственно, в одинаковой пространственной области. Во-вторых, некоторые исследователи полагают, что в деятельности кортекса в значительной мере участвуют его слои, а не отдельные нейроны. Столбец (колонна) кортекса состоит изо всех тканей (нейроны, глия, аксоны) внутри цилиндрического объема площадью около 1 мм2 на поверхности коры, распространяющемся радиально вглубь через все слои коры.

В данной статье описана одна из математических моделей мезоскопической электрической активности нейронов коры.

Чтобы смоделировать электрическую активность мозга человека, использовалась система из 8-ми нелинейных стохастических частных дифференциальных уравнений и 20-ти переменных путем их преобразования в безразмерную форму [2]. Эта система представлена ниже:

% = 1 - И + г, (И0 - И )ІЄЄ +г (И0 - ие)!,, ді (1)

дИ = 1 - И. +Г (И0 - И. )І +Г (И0 - И. )І. 1 е ч е 1 ■> Є1 /V/ 1 ■> 11 ді (2)

1 д (—— + 1)2І = МРБ[И] + ф + Р +П> Ч™ — У ее е е\- е і те ее 1 Те ді (3)

1 д (—— +1)21 . = N [И ] + ф. + Р . + Г,> е. е е1- е J г. е. 2 Т е д (4)

1 д (- - +1)2 4 = N. [И. ] + Ріе +Г3, Ті ді (5)

1 д (- т; +1)21 и = N,в 8 [И, ] + Р. +Г4, Т. ді (6)

(-7- д+1)2 Фе =Л Ц-+(7- дт+[Ие ]. 7 ді 72 дх 7е ді (7)

(' д+ 1)2ф = а+ (' д+1)NЧ[И,]■ 7 ді ' 7 д х 7 ді . е е (8)

Восемь динамических переменных в этой системе (Ие, И, і,,, Іеі, і,, I й Фе, Ф) - функции безразмерного пространства х и времени і. Безразмерные переменные и параметры определены в табл. 1 и 2 соответственно. Безразмерные переменные (левая колонка) определены символами в средней колонке. Переменные описаны в правой колонке. Индексы е и і

соответствуют процессам возбуждения и торможения. Сделаны следующие условные

_Л _ ___ с тах_сша^сшах и гев^игевОиге^

упрощения: те-т,-т, ае -а, -а и д -д -п .

Таблица 1

Динамические переменные для безразмерной модели дифференциальных уравнений

Переменная Определение Описание

Ь-е.і Ье,і/Ьге5‘ Мембранный потенциал

Іее,іе Іее,іеїе/(0еехр(1)8шах) Совокупность е——е, 1——е клеток возбуждения

Іеі.іі Іеіііуі/(Оіехр(1)8тах) Совокупность е—1, 1—1 клеток торможения

Фе,і Феі/§™х Ряд е и 1

1 іУт Время

X х/(ху) Пространство

Таблица 2

Значения параметров для безразмерной макромодели

Переменная Определение Описание Типичные значения

е,і (Как индексы) совокупность клеток возбуждения или торможения

Ге,і ОеДехр(1)8тах/ТеД|Ье,Г-Ь.ге5‘| Влияние мощности на значение мембранного потенциала 1,42х10-3; 0,0774

Ье,° Ьеіге7Ь.ге!і‘ Обратный потенциал -0,643; 1,29

Те,і Тїе,і Константа 12,0; 2,6

^е,і ТУЛее,еі Безразмерная характеристика в инверсной шкале 11,2; 18,2

р . ее,іе р /стах г ее,іе'13 Подкорковый потенциал возбуждения 11,0; 16,0

р • •• еі,іі Реі,іі/стах Подкорковый потенциал торможения 16,0; 11,0

N а і>|е,і - Общее число синаптических окончаний от клеток возбуждения 4000,2000

N в і>|е,і - Общее число местных синаптических окончаний е и і 3034, 536

§е,і ёедЬге!!1 Сигмоидальный наклон в точке сгиба -19,6; -9,8

9е,і 9е>і/Ьге5‘ Точка изгиба сигмоидальной функции 0,857; 0,857

Уравнение (1) определяет эволюцию безразмерного напряжения возбудительных нейронов через само напряжение (кг), входной ток возбудительных нейронов (1ее) и входной ток тормозных нейронов (4). Уравнения (3) и (5) описывают динамическое поведение входных токов. Входной ток возбуждающих нейронов (1ее) устанавливается в соответствии с локальной мощностью (]^е £><&]), удаленной мощностью возбудительных нейронов коры (Фе), подкорковой мощностью (Рее) и стохастической мощностью (Г;). Член 8г[к] -безразмерная сигмоидальная передаточная функция

S [h ] =--------1---------, (9)

л eJ 1 + exp[-ge(he -e,)] которая необходима для преобразования локального напряжения возбудительных нейронов в среднюю скорость работы возбудительной популяции (непрерывно-импульсное преобразование). Мощность нейронов удаленного кортекса, которая всегда возбудительная, устанавливается в соответствии с уравнением распространения (7) диффузного типа. Для вычисления неизвестных подкорковых мощностей добавляется последний член Г1 в (3), который определяется как безразмерная стохастическая мощность

Г =а^Лу[ x, t ], (10)

где - источник Гауссова белого шума с нулевым значением. Для вычислений аппроксимируем как

йм=RM-. (ii)

V Ax At

где x=mAx и t=nAt, (m, n - целые числа) определяют пространственные и временные координаты в двухмерном пространстве Ax и At соответственно. Похожие уравнения определены для динамики нейронов торможения (2), местных токов для нейронов торможения (4) и (6) и отдаленных возбудительных токов для нейронов торможения (8).

Самая важная переменная в модели - переменная he - усредненный по поверхности соматический мембранный потенциал возбудительных клеток кортекса. В уравнениях (1)-(8) безразмерный параметр he = h У(-70 мВ). Для того чтобы вызвать приступоподобное состояние, в модели необходимо поменять два параметра, связанные с процессом возбуждения в модели: увеличить параметр Рее - силу тока возбуждения в подкорковом веществе - до 2000 % и уменьшить параметр Ге - влияние возбудительных токов на значение мембранного потенциала сомы возбудительных и тормозных клеток - до 30 %. Затем вычисляются значения уравнений (1)-(8) для одномерного пространства и определяется активность бегущей волны у he.

Данная модель позволяет проводить динамический анализ эволюции комплексной электрической активности.

БИБЛИОГРАФИЧЕСКИЙ СПИСОК

1. Жадин М.Н. Биофизические механизмы формирования электроэнцефалограммы. - М.: Наука, 1984. - 197 с.

2. Mark A. Kramer, Andrew J. Szeri, James W. Sleigh. Mechanisms of seizure propagation in a cortical model. (2006).

УДК 616.12-073.97

С.В. Хашева, А.Н. Калиниченко МЕТОДЫ ФИЛЬТРАЦИИ ЭКГ МАТЕРИ И ПЛОДА

В ходе наблюдения за протеканием беременности, особенно беременности с осложнениями, возникает необходимость анализа и мониторинга сердечной деятельности плода. Возможны два варианта для съема ЭКГ плода: инвазивный и неинвазивный.

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

i Надоели баннеры? Вы всегда можете отключить рекламу.