CC BY
15В. Р. Цибульский, О. И. Сергейчик, М. В. Абрамов
МЕТОДЫ ПОЛУЧЕНИЯ ПЕРЕДАТОЧНОЙ ФУНКЦИИ НА ОСНОВЕ ЧАСТОТНЫХ ХАРАКТЕРИСТИК ЭКГ
Предложены варианты математической модели функционирования сердца и сердечно-сосудистой системы человека и варианты получения передаточной функции с последующей оценкой устойчивости моделей.
Одной из главных причин высокой смертности в развитых странах являются заболевания сердечно-сосудистой системы. Рост количества людей в среднем и пожилом возрасте с заболеваниями сердца — важная проблема в современном здравоохранении. Выявление таких заболеваний на ранних стадиях затрудняется тем, что основным его методом по сей день остается визуальный анализ электрокардиограммы врачом-специалистом.
Применение компьютерных вычислительных пакетов для оценки оцифрованного сигнала ЭКГ по какому-либо параметру [11] позволяет уменьшить время и повысить точность анализа результатов обследования. Затруднение при этом заключается в отсутствии математической модели сердечной деятельности, которая охватывает определяющие процессы в сердце.
Основной моделью, описывающей работу сердца с учетом ряда возможных патологий, можно считать четырехкамерную динамическую модель сердца, связывающую более двадцати параметров, основными из которых являются давление и объем в камерах сердца, кровоток между отдельными камерами [12]. При этом учитываются закон сохранения массы крови, перемещаемой между желудочками и предсердиями, с использовыванием принятых в физиологии формул Пуазейля для определения кровотока; состояния клапанов сердца, и, таким образом, может быть построена общая динамическая модель кровеносной системы [12].
Представляет интерес и модель, описывающая напряженно-деформированное состояние стенок сердца [1]. Она позволяет оценить нагрузку на ту или иную часть сердца, определить вероятность превышения нагрузки и возможной дисфункции при тех или иных состояниях человеческого организма.
Более сложная модель, описывающая взаимодействия всех элементов сердца [6], обусловлена системой закономерностей, фундаментальных в физиологии. Эта модель учитывает ряд патологий и включает в себя уравнения основных процессов сердечной деятельности: механики сердечной мышцы, сокращения и возбуждения отделов сердца, гомеометрической саморегуляции, геометрии сердца, работы клапанов сердца. Причем при таком большом объеме закономерностей дифференциальные уравнения модели имеют довольно простой вид. Автором дано описание программной реализации модели.
Охарактеризованные модели не позволяют достаточно просто оценить, например, устойчивость функционирования сердца. В данной статье предлагается получить модель на основе передаточной функции [2, 3]. Ее можно представить, если рассматривать сердце как объект управления, находящийся в режиме самовозбуждения [7]. Передаточная функция системы (сердца) может быть получена через частотную характеристику.
В общем виде передаточная функция Ж(р) определяется как отношение изображений выходной характеристики У(р) к входной (управляющей) и(р):
Нормированный вид передаточной функции:
л 2 n ,
1 + a • p + a2 • p +... + an • p
(2)
где к — коэффициент усиления; аи Ьг- — коэффициенты передаточной функции; р — оператор дифференцирования.
Важным свойством передаточной функции является взаимосвязь между ней и временной (или импульсной переходной) характеристикой объекта, выраженная следующей формулой:
где h(t) — импульсная переходная характеристика.
Вместе с тем существует ряд объектов управления, определить передаточную функцию которых затруднительно. Примером такого объекта является сердце человека. Идентификация передаточной функции сердца классическим путем, т. е. со входом и выходом, невозможна, так как в данном случае сложно определить, что является входом (U(p)).
Предлагается несколько вариантов определения передаточной функции работы сердца:
1. Предположеним, что при остановке сердца у больного дается короткий и мощный импульс, после чего работа сердца возобновляется и наблюдается ЭКГ. Используя импульсную переходную характеристику и частотную характеристику ЭКГ пациента, перенесшего остановку и запуск сердца, получаем передаточную функцию на основании формулы (3).
2. Переход от модели сердца к моделям кардиологических процессов, которые вызывают отклонения или, проще говоря, заболевания. В этом случае за вход U(p) можно принять ЭКГ [4] и ее спектр у пациента в норме, а за выход Y(p) — ЭКГ и ее спектр у больного того же возраста [9]. Частотная переходная функция будет представлена отношением частотных характеристик выхода и входа. При использовании нормальной ЭКГ получаем модель болезни и оценка устойчивости улучшения состояния будет другая, может быть некая среднеквадратическая характеристика, стремящаяся к нулю.
3. Использование методики определения передаточной функции через модель генератора, например RC [8]. В качестве исследуемой берется одна из известных RCL-моделей. Отличие будет лишь в том, что в нормальном состоянии амплитуда, форма ЭКГ и частота стабильны, значит, для модели больного, возможно, следует ввести нелинейную обратную связь [8].
Рассмотрим первый метод подробнее.
Суть метода заключается в получении на выходе объекта временных (или импульсных переходных) характеристик при подаче ступенчатого или импульсного сигнала на входе. Неизвестные параметры передаточной функции определяются либо исходя из связи передаточной функции и импульсной переходной характеристики (3), определяющей передаточную функцию как изображение по Лапласу импульсной переходной характеристики, либо с помощью методов идентификации параметров [10].
W(p) = L[h(t)],
(3)
Помимо данных методов, для выходных сигналов, имеющих периодический вид, возможен метод получения передаточной функции на основе частотных характеристик ЭКГ, полученных методом Фурье [2, 5, 11]. Основываясь на последнем, можно сформулировать следующий алгоритм:
1. Получение спектра сигнала при помощи алгоритма быстрого преобразования Фурье (периодограмма).
2. Получение аналитической записи сигнала — ряд Фурье.
3. Получение передаточной функции посредством преобразования Лапласа, примененного к данному ряду Фурье.
Аналитическая запись сигналов в виде ряда Фурье будет иметь вид:
а ^
х(*) = “2 + 2 (ап С0<Пп ■*) + ЬП • *)) . (4)
2 п=1
Для простоты получения первого результата выберем линейную модель, для которой справедлив принцип суперпозиции и любые возмущения можно свести к одному входу.
Преобразование Лапласа для формулы (4) с учетом удаленного среднего будет представлять собой сумму дробей вида
£[*)] = 2 ап- Р + Ьп 2“п . (5)
Я Р2 +пп
Основываясь на формуле (3), при условии что х(*) — импульсная переходная характеристика, можно говорить о том, что выражение (5) является передаточной функцией Ж (р).
Поскольку количество частот может быть бесконечным, о чем говорит верхний предел суммирования, то главным вопросом при использовании данного метода является определение набора частот, которые адекватно определяют систему.
Пусть верхний предел равен М, тогда выражение (5) можно представить в следующем виде:
А0 • р2 ■ (М-1) + А1 • р2 ■ (М-2) +... + А2 (М 2) • р2 +1 Ж(р) = К• -----------—-----------2-&-2) г------ (6)
ггуу) 1V • N 2• (М-1) , , г> 2 . л ■
р + р + ... + Вм-1• р + 1
Выражение (6) имеет тот же вид, что и выражение (2). Степенью данной передаточной функции является удвоенное значение количества частот, выбранных для описания системы моделью (5).
Получив передаточную функцию системы, можно использовать различные методы для анализа ее устойчивости. К ним можно отнести критерий Ляпунова, Гурвица, Найквиста и др. [3].
В качестве иллюстрации применения метода предлагается пример получения передаточной функции реальной ЭКГ больного.
На рис. 1 приведен вид электрокардиограммы больного, перенесшего остановку и запуск сердца (данные из Тюменского кардиологического центра). Этот сигнал можно считать выходной характеристикой сердца.
Рис. 1. Электрокардиограмма больного
Для него получен спектр по алгоритму БПФ, для 500 частот. Амплитудночастотная характеристика (для нижних частот) представлена на рис. 2.
Рис. 2. Амплитудно-частотная характеристика электрокардиограммы
Для упрощения определения передаточной функции предложенным методом выберем всего пять частот, ограничив спектр на уровне 15 % от максимальной амплитуды. На рис. 2 эти частоты отмечены более жирными линиями, значения их амплитуд приведены в табл. 1.
Параметры передаточной функции на основе первой выборки частот
і Юі а Ьі
1 0,012 0,016529072 -0,013810079
2 0,015 -0,01749109 -0,004379352
3 0,018 0,002004919 0,018348789
4 0,019 0,009302966 0,015070977
5 0,029 0,001980545 -0,016824831
Передаточная функция сердца при данной выборке определяется выражением
ч 7,7 • 1016 • р9 - 6,42 • 1014 • р8 +1,2 • 1014 • р7 + 2,75 • 1011 • р6 + 6,47 • 1010 • р5 +
Ж (р) = —---------------------------------------------------------------—
6,25 • 1018 • р10 +1,18 • 1016 • р8 + 8,05 • 1012 • р6 + 2,52 • 109 • р4 + (7)
+ 2,5 • 108 • р4 +1,53 • 107 • р3 +1,51 • 104 • р2 +1,42 • 103 • р - 4,19.
+ 3,67 105 • р2 +19,9
Основываясь на утверждении о том, что реальные части корней характеристического уравнения устойчивой системы лежат в отрицательной полуплоскости [3], найдем решения данного уравнения для полученной модели.
Полагая систему замкнутой, запишем характеристическое уравнение [3]:
1 + Ж (р) = 0 (8)
Решения уравнения (8) для данной выборки приведены в табл. 2.
Таблица 2
Корни характеристического уравнения первой выборки
№ Р
1 -0,01±0,023і
2 -2,995 10"3±7,845-10"3і
3 1,205 10-4±0,018і
4 1,076 10-3±0,015і
5 5,66110-3±0,026і
Поскольку есть реальные части корней характеристического уравнения, которые лежат в положительной полуплоскости, мы не можем говорить о том, что работа сердца устойчива.
Вместе с тем при выборе других частот можно получить совсем иной результат. Выберем частоты, относительно равномерно распределенные по всему диапазону (табл. 3).
Таблица 3
Параметры передаточной функции на основе второй выборки частот
і Юі а Ьі
1 0,007 0,006378 -0,0006
2 0,022 0,015886 -0,00201
3 0,031 0,012162 -0,002
4 0,043 0,006931 0,000316
5 0,048 0,008406 -0,00127
Передаточная функция сердца при данной выборке определяется выражением
, 1,94-1016 • р9 - 6,161013 • р8 + 8,95 • 1013 • р1 - 2,7110й- р6 + 1,3210й- р5 -
Ж (р) =--------------------------------------------------------------------
3,9 • 1017 • р10 + 2,2 1015- р8 + 4,3 • 1012 • р6 + 3,37-109 • р4 + (9)
- 3,79-108 • р4 + 6,61107 • р3 -1,66 -105 • р2 + 6,94-103 • р - 9,89 + 9,29 105 • р2 + 37,9
Решения характеристического уравнения (8) для новой выборки приведены в табл. 4.
Таблица 4
Корни характеристического уравнения второй выборки
№ Р
1 -0,028
2 -6,851 -10-3
3 -4,13510-3±0,011 /'
4 -1,81410-3±0,041 /'
5 -1,082 10"3±0,028/
6 -4,11-10"4±0,046г
Поскольку реальные части корней характеристического уравнения отрицательны, мы можем говорить о том, что после запуска сердце работает устойчиво.
На основании приведенных примеров можно говорить о том, что применение данного метода возможно, но основной трудностью является правильный выбор частот, использование которых позволяет описать адекватную модель. Для более уверенной оценки необходимо исследование всех методов получения передаточной функции, предложенных выше.
ЛИТЕРАТУРА
1. Бегун П. И., Афонин П. Н. Моделирование в биомеханике: Учеб. пособие. — М.: Высш. шк., 2004. — 390 с.
2. Бендат Дж., Пирсол А. Прикладной анализ случайных данных / Пер. с англ. д-ра физ.-мат. наук В. Е. Привальского; Под ред. И. Н. Коваленко — М.: Мир, 1989. — 540 с.
3. Бесекерский В. А., Попов Е. П. Теория систем автоматического управления. — Изд. 4-е, перераб. и доп. — СПб.: Профессия, 2003. — 752 с.
4. Вартак Ж. М. Интерпретация электрокардиограммы. — М.: Медицина, 1978. — 152 с.
5. Марпл К. (мл.). Современные методы спектрального анализа // ТИИЭР. — 1981, ноябрь. — Т. 69, № 11.
6. Мосткова Е. В. Математическая модель сердца // Клиническая физиология кровообращения. — 2006. — № 4. — С. 26.
7. Парашин В. Б., Иткин Г. П. Биомеханика кровообращения: Учеб. пособие. —
М.: Изд-во МГТУ им. Н. Э. Баумана, 2005. — 224 с.
8. Ройтман М. С. Принципы построения, вопросы теории и практики прецизион-
ных источников переменных напряжений: Автореф. ... дис. д-ра техн. наук. — Томск, 1971.
9. Рычкова А. Ю., Цибульский В. Р., Сергейчик О. И., Копылова Л. Н. Спектральный анализ нормальной ЭКГ // Вестн. аритмологии. — 2004. — № 35. — С. 52.
10. Сарычев А., Искандеров Г. Методы идентификации объектов управления на
основании результатов активного эксперимента // Материалы XXXI науч. конф. по результатам работы ППС / Северо-Кавказский гос. техн. ун-т. — Режим доступа: http://www/ncstu.ru/Science/conf/past/2001/xxxi/14/19.pdf, свободный.
11. Сергейчик О. И. Модели и алгоритмы спектрального анализа обработки кардиологических временных рядов: Автореф. ... дис. канд. техн. наук. — Тюмень, 2007.
12. Холодов А. С. некоторые динамические модели внешнего дыхания и кровообращения с учетом их связности и переноса веществ // Компьютерные модели и прогресс медицины / Под ред. О. М. Белоцерковского. — М.: Наука, 2001. — 300 с.
V. R. Tsibulsky, O. I. Sergeychik, M. V. Abramov
METHODS OF OBTAINING TRANSFER FUNCTION BASING ON FREQUENCY CHARACTERISTICS OF ELECTROCARDIOGRAM
The article suggests alternatives of a mathematical model with respect to man’s heart and cardiovascular system together with variants of obtaining transfer function with a subsequent evaluation of stability of these models.
