Свойства, обобщенные модели и задачи исследования тонкой структуры биоэлектрических сигналов

Размещено на http://www.allbest.ru/

электрокардиограмма помеха

Свойства, обобщенные модели и задачи исследования тонкой структуры биоэлектрических сигналов

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

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

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

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

Это обстоятельство ограничивает разрешающую способность любых методов измерения параметров.

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

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

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

Наибольшее влияние во всех без исключения исследованиях оказывают следующие виды помех:

- эффект поляризации электродов, приводящий к смещению нулевого уровня сигнала;

- квазигармонический процесс, представленный составляющими наводки напряжения промышленной частоты;

- артефакты смещения электродов, создающие выбросы случайной амплитуды и длительности;

- электрофизиологические помехи.

Реализация сигнала электрокардиограммы при наличии указанных видов помех приведена на рис. 1.

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

Степень компенсации помех определяется как характеристиками применяемых методов, так и соотношением амплитудных, спектральных и временных составляющих информационного и помехового процессов.

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

Рисунок 1 - Реализация кардиосигнала с помехами:

1 - шум ЭМГ; 2 -помеха промышленной частоты 50 Гц; 3 - смещение изоэлектрической линии вследствие эффекта поляризации электродов; 4 - артефакт смещения электрода

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

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

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

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

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

Для оценки соотношений между полезным сигналом и помехами в табл. 1 приведены значения амплитудно-временных параметров ЭКГ, соответствующих норме, а в табл. 2 - основные статистические параметры помех.

Соотношения частотных свойств полезного сигнала и помех представлены на рис. 2 в виде графиков СПМ различных составляющих.

Основная мощность QRS-комплекса сосредоточена в области частот 2-20 Гц с наличием максимума на частоте около 15 Гц.

Спектр ЭКГ-сигнала может изменяться в зависимости от морфологии сигнала. Спектр шумов от мышц является неоднородно распределенным и характеризуется значительной вариабельностью.

Таблица 1 Параметры нормальной ЭКГ

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

Таблица 2 Параметры аддитивных помех, сопровождающих регистрацию ЭКГ

Рисунок 2 - Характеристики относительной спектральной плотности мощности ЭКГ сигнала и шумов: 1 - ЭКГ-сигнал; QRS - комплекс; 3 - артефакты движения; Р-, Т-зубцы; 5 - напряжение поляризации; мышечные шумы (ЭМГ)

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

Актуальны также и задачи, которые на современном этапе целесообразно решать только с применением статистических методов обработки, в частности:

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

определение существования зубца Q, оценивание момента его начала;

оценивание момента окончания QRS-комплекса (точка J);

оценивание смещения сегмента ST относительно изоэлектрической линии.

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

Одним из наиболее значимых направлений исследований в данной области является изучение потенциалов замедленной деполяризации желудочков сердца. По предположению, подтверждаемому инвазивными методами исследований сердца (методы записи сигналов непосредственно с ткани миокарда), некоторые патологии сопровождаются возникновением слабых «всплесков» сигналов случайной формы в пределах зубца S и сегмента ST.

Амплитуда таких сигналов на поверхности тела весьма мала и обычно оказывается соизмеримой с уровнем шума ЭМГ даже при полном расслаблении мышц пациента. Задача исследования состоит в том, чтобы разрабатывать методы надежного обнаружения как можно более слабых компонент в ЭКГ при неинвазивном съеме.

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

С учетом изложенного в общем виде математическая модель ЭКС с помехами может быть представлена следующим образом:

, 1, 2, …, L,

где - высокоамплитудные составляющие сигнала ЭКГ (комплекс P-QRS-T), характеризующиеся вектором неизвестных параметров ; - низкоамплитудные составляющие сигнала ЭКГ, характеризующиеся вектором неизвестных параметров ; - реализация помехи на протяжении -го кардиоцикла.

Отметим, что представление ЭКС моделью общего вида (1), в принципе, неоднозначно при отсутствии дополнительных априорных сведений о свойствах , с учетом их малой амплитуды << . Распространено предположение, подтверждаемое инвазивными исследованиями, что низкоамплитудные составляющие сигнала ЭКГ являются сравнительно высокочастотными (ориентировочно диапазон частот 40-1000 Гц).

Для исследования поздних потенциалов желудочков сердца и любых низкоамплитудных сигналов вообще существует эффективный подход, основанный на свойстве повторяемости сигнала ЭКГ и реализующий принцип псевдосинхронного накопления.

В этом случае оцениваются все интервалы между кардиоциклами например, {, ? = 1, 2 …, L}, если за опорную точку считать максимум R-зубца QRS-комплекса, после чего вычисляется сигнал:

,

где - сигнал усредненной ЭКГ.

В предположении независимости отсчетов шума в усредняемых кардиоциклах (нормальное распределение) дисперсия шума в (2) уменьшается по сравнению с исходной в L раз, т. е. отношение сигнал/шум по напряжению улучшается в раз.

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

В процессе обработки приходится также решать задачи классификации форм комплексов с целью исключения атипичных комплексов и артефактов, которые приводят к значительным ошибкам при усреднении сигнала. Теоретически улучшать разрешающую способность по амплитуде можно неограниченно, однако время обследования пациента в расслабленном состоянии ограничивает количество доступных для анализа кардиоциклов величиной L ~ 500, помимо этого, суточные циклы изменения формы QRS-комплексов могут привести к неоднородности их формы в начале и конце записи, что исказит конечный результат.

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

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

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

Можно выделить два частных случая нестационарности проявления микропотенциалов (на основании изучения данных внутрисердечных исследований):

-микропотенциалы повторяются по форме от одного кардиоцикла к другому, но возникают в случайный момент времени, (т. е. «не синхронизированы» с высокоамплитудными составляющими ЭКГ-сигнала); эквивалентная модель кардиосигнала имеет вид

, 1, 2, …, L,

где - случайный момент появления низкоамплитудного сигнала ква-зидетерминированной формы в пределах кардиоцикла;

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

, 1, 2, …, L,

т.е. их параметры изменяются от одного кардиоцикла к другому (это отражает индекс ).

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

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

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

.

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

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

Детерминированные функции в ЭЭГ отсутствуют, фоновая активность (имеет характер случайного процесса) в этом случае рассматривается как помеха, а вызванные потенциалы могут иметь как стационарный, так и нестационарный характер проявления по времени появления или по форме.

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

Размещено на Allbest.ru