Курсач Мой по Мобсу (1040401), страница 3
Текст из файла (страница 3)
Для исследования пространственных отношений электрической активности коры большого мозга был сконструирован прибор — электроэнцефалоскоп, позволяющий одновременно наблюдать электрическую активность 50 или 100 точек. С помощью электроэнцефалоскола был выявлен ряд закономерностей в пространственном протекании биоэлектрических процессов. Метод электроэнцефалоскопии был применен для анализа ЭЭГ здорового человека и больных с очаговыми поражениями мозга
Впервые вопрос о возможности применения корреляционного метода к ЭЭГ был поставлен Н. Винером в 30-х годах и разработан им в сотрудничестве с физиологами и физиками Массачусетского технологического института. Следует сказать, что независимо от американских исследователей, еще до их первых публикаций, корреляционный анализ ЭЭГ применили японские физиологи, проводившие весь корреляционный анализ с помощью очень трудоемкой обработки кривых вручную.
Частотный и автокорреляционный анализы в своем конечном выражении — спектре мощности — теоретически эквивалентны. Выбор типа анализа зависит от конкретных физиологических задач. Так, если ставится вопрос о временных соотношениях процессов, то целесообразно использовать корреляционный анализ. Если же ставится вопрос о наличии или отсутствии в процессе тех или иных частот, то можно пользоваться частотным анализом.
При корреляционном анализе производится исследование изменений процесса во времени [15].
Более того, корреляционный анализ ЭЭГ позволяет выявить и количественно оценить следующие стороны процесса: 1) наличие или отсутствие в ЭЭГ периодического (или квазипериодического) процесса; 2) период ритмических колебаний; 3) устойчивость периодики; 4) при кросскорреляции — определить временные отношения двух процессов по их фазовому сдвигу или сдвигу максимума кросскор-реляционной функции; 5) количественно оценить степень связи или сходства двух процессов ЭЭГ. Здесь необходимо отметить, что вычисление корреляционной функции ЭЭГ, в частности взаимной или кросскорреляции, а также, как это будет показано ниже, когерентности электрических колебаний потенциала не эквивалентно вычислению попарной корреляции амплитудных изменений в двух ЭЭГ, получаемой на основании данных энцефалоскопа.
При вычислении когерентности выявляются общие для двух ЭЭГ колебания, исходящие из единого источника, проводится оценка сходства по сравнению с остальными колебаниями даже совпадающей частоты, происходящими из другого источника. При корреляции амплитудных изменений двух ЭЭГ (по данным энцефалоскопа и других подобных приборов) сравниваются все колебания не только общей частоты, но и разного периода, возникшие на кривой ЭЭГ в данный момент измерения, т. е. может быть проведено сравнение увеличения амплитуды α-волны с падением δ-волны.
С применением корреляционного анализа ЭЭГ было показано, что α-ритм фоновой активности ЭЭГ и ритмика с частотой альфа, возникающая в форме разряда последействия в ответ на световое раздражение, развиваются независимо друг от друга, между ними нет фазовой связи.
Среди математических методов анализа ЭЭГ внимание исследователей привлекает метод вычисления спектров мощности, взаимных спектров и комплексной функции когерентности.
Спектр мощности отражает энергию каждой из частотных составляющих данной ЭЭГ. В каждой же ЭЭГ имеется, как правило, не одна, а несколько различно выраженных частотных составляющих. Вычисление спектра мощности позволяет выделить с большой точностью частотные составляющие каждого процесса ЭЭГ и показать веса разных частотных компонентов. Взаимные спектры дают количественные оценки связанности отдельных частотных составляющих ЭЭГ двух областей [13].
Частотный анализ.
Частотный анализ биопотенциалов мозга, в данном случае, производится при помощи электроэнцефалографического комплекса «Биофизприбор», который включал 8-канальный электроэнцефалограф, 8-канальный анализатор с широкополосовыми фильтрами, 8-канальный интегратор.
Анализатор имеет набор из 6 полосовых фильтров, в каждом 8 каналов с соответствующими значениями граничных частот: фильтр «δ» — 0,35—4,2 Гц; фильтр «θ» —4,6—7,4 Гц; фильтр «α» —7,7—13,4 Гц; два фильтра диапазона β-волн—«β1 » (низкий) — 13,7—21 Гц и «β2 » (высокий) — 21—30 Гц; фильтр «γ» — 31—70 Гц. Фильтры избирательно выделяют только те колебания ЭЭГ, частоты которых соответствуют полосе пропускания. Выходные сигналы фильтров характеризуют колебания в частотных полосах физиологических ритмов. Ширина полосы пропускания фильтром позволяет осуществлять анализ за короткие отрезки времени. Эпоха анализа сокращается до долей секунды (фильтр имеет небольшое время установления — постоянная времени составляет примерно 0,13 с). Фильтры обладают высокой селективностью.
Каждому каналу анализатора соответствует канал интегратора. Принцип его работы: измеряемые биопотенциалы подвергаются двухполупериодному выпрямлению, и выпрямленное напряжение поступает на интегратор. Последний измеряет суммарную биоэлектрическую активность пропорционально накопленному за эпоху анализа электрическому заряду. Интегратор с импульсным отсчетом отмечает определенные порции электрического заряда и соответственно этому посылает отметки на осциллограммы или на цифровой счетчик. Показания интегратора, отмечающиеся на цифровых электромеханических счетчиках, дают возможность определять средние значения амплитуд выделенных ритмов или ЭЭГ за любые более или менее длительные периоды времени, заданные заранее, и получать непосредственное цифровое выражение результатов анализа.
Показания интегратора приводятся в его относительных единицах, поскольку для оценки каждой ЭЭГ важно количественное соотношение выраженности ритмов разных диапазонов. При постоянном значении единицы интегратора и при одинаковом коэффициенте усиления имеется возможность сравнения средней активности ритмов в ЭЭГ разных отведений, а также в ЭЭГ разных людей. В ряде исследований вычислялись в процентах выраженности ритмов в каждой из исследованных областей, а также процентное отклонение в случае патологии того или иного ритма от среднего уровня в норме.
В одном из каналов анализатора предусмотрена возможность исследования частоты в любом из диапазонов ритмов с помощью периодметра. Действие периодметра основывается на измерении длительности отдельных периодов колебаний. На осциллограмме периодметр регистрирует соответственно каждой волне исследуемого диапазона импульс, высота которого пропорциональна длительности периода волны.
На рис. 10 приведены данные частотного анализа ЭЭГ здорового человека. Видно, что ЭЭГ содержит колебания всех диапазонов, выраженность которых различна; наибольшая выраженность — в диапазоне а-ритма.
При анализе ЭЭГ здоровых и больных людей эпоха анализа была выбрана длительностью в 5 с. Оценка средних значений фоновой ритмики производилась усреднением значений 5—10 замеров интегратора.
Для оценки изменений средней суммарной активности во время раздражений и в последействии вычислялось ее отношение к активности в фоне непосредственно перед раздражением, причем последняя принималась за 100%. Изменения показателей интегратора во время раздражения и в последействии оценивались в процентах [16].
Корреляционный анализ
Корреляционный анализ ЭЭГ производится с помощью автоматического коррелографа (конструктор Е. Н. Бурашников).
Принцип действия коррелографа состоит в предварительной записи напряжения биотоков на магнитную пленку и последующем автоматическом вычислении авто- или кросскорреляционной функции для выбранного отрезка записи.
Технические характеристики коррелографа: 1) число каналов записи и воспроизведения — два; 2) спектр частот анализируемых сигналов — от 0 до 1000 Гц; 3) максимальное время задержки τmax — 2,56 с; 4) среднеквадратическая ошибка измерений ±15% от значения корреляционной функции при τ=0, вычисление производится для дискретных значений аргумента в соответствия с формулой:
где Т – время интегрирования, ∆τ – шаг дискретизации по времени, κ·∆τ – аргумент корреляционной функции.
Результат интегрирования в виде постоянного напряжения, пропорционального значению корреляционной функции, поступает на электронный потенциометр и фиксируется в виде точки на диаграммной ленте. После этого величина аргумента автоматически изменяется на заранее заданную величину, и процесс вычисления повторяется.
Корреляционный анализ биотоков после предварительной фильтрации позволяет увеличить точность вычисления функций отдельных периодических и квазипериодических составляющих, имеющихся в ЭЭГ, и уменьшить время вычисления [20].
Шаг задержки, который применялся при обработке электрической активности мозга, был различным в зависимости от того, какой ритм исследовался. При анализе полной ЭЭГ применялся шаг задержки (∆τ), равный 10 мс. При обработке выделенных фильтрами диапазонов частот в зависимости от полосы ритмов шаг задержки был различным: при анализе α-ритма, так же как и при полной ЭЭГ, он составлял 10 мс; при анализе β-ритма — 5 мс; при обработке выделенного θ-диапазона — 20 мс; при обработке δ-диапазона — 40 мс. При определении шага задержки для обработки процесса исходили из того, чтобы на период колебания основной частоты приходилось по крайней мере 10 точек вычисленной корреляционной функции.
При корреляционном анализе весьма существенным является выбор длительности интервала наблюдения. Корреляционная функция может быть вычислена точно лишь в случае бесконечно большой величины интервала наблюдения. При исследовании ЭЭГ физиологи всегда стремятся сократить время наблюдения, так как особенности динамики нервных процессов, отраженные на ЭЭГ, имеют свои закономерности.
На основании проведенных исследований с различными интервалами наблюдения ЭЭГ человека, мы избрали интервал наблюдения, равный 10 с, который использовался как при обработке фоновой записи полной ЭЭГ и выделенных диапазонов ритмов, так и во время действия раздражений.
Для сравнения значения корреляционных функций разных по амплитуде ЭЭГ или электрограмм выделенных ритмов производилось нормирование: максимальное значение автокорреляционных функций всех анализируемых ЭЭГ и выделенных частот приводилось к одинаковой величине отклонения от нулевой линии на коррелограмме.
При корреляционном анализе можно вычислять автокорреляционную и кросскорреляционную (или взаимно-корреляционную) функции ЭЭГ. Использование корреляционного метода для анализа ЭЭГ обусловлено тем, что построение автокорреляционных функций является наиболее мощным средством выделения гармонических колебаний в случайном процессе. Вычисление кросскорреляционной функции двух ЭЭГ позволяет определить степень связи процессов в разных точках мозга и выделить периодические составляющие, общие для обеих ЭЭГ [12].
На рис. 11 в качестве примера разных форм автокоррелограммы (АКГ) приводятся: верхняя АКГ, вычисленная для ЭЭГ затылочной области обследуемого «без α-ритма», средняя — для ЭЭГ затылочной области с отчетливой перестройкой в ритме подаваемых мельканий 12,5 Гц, нижняя — для ЭЭГ затылочной области с хорошо выраженным α-ритмом.
При автокорреляционном анализе выделенных полосовыми фильтрами диапазонов физиологических ритмов форма АКГ изменяется соответственно соотношению в них случайных (непериодических) и квазипериодических колебаний.
Кросскорреляционный анализ ЭЭГ дает принципиально новые возможности оценки процессов двух точек мозга — позволяет количественно оценить степень сходства процессов или их связи, выявить общие компоненты и их соотношение, а также временные отношения разных ритмов. Вычисление кросскорреляционной функции двух ЭЭГ позволяет ближе подойти к раскрытию механизмов и путей функциональных связей между активностью разных отделов мозга.
Спектральный анализ