CC BY
109
ким кругом пользователей, не специализирующихся в области вычислительной техники и программирования.
Изменяя значения весовых коэффициентов тех или иных показателей, можно создавать требуемые оценочные предпочтения, а также учитывать логические условия получения результата.
Определение состава показателей и их весовых коэффициентов осуществляется в основном субъективными экспертными методами. С учётом этого штабная коэффициентная методика с достаточным количеством промежуточных показателей, вводимых в целях их использования при анализе получаемых результатов, представляет собой простейшую экспертную систему, так как содержит присущие любой экспертной системе основные элементы - подсистему вывода знаний, реализованную на основе суммирования значений показателей с умножением их на экспертно определяемые значения весовых коэффициентов, и подсистему объяснения в виде совокупности значений промежуточных показателей, получаемых при конкретных значениях исходных показателей.
Простота построения коэффициентной методики, отсутствие сложного математического аппарата вычислений результатов оценки, а также возможность применения программной реализации на персональных электронновычислительных машинах (ПЭВМ) для решения задач оценки (в том числе рейтинговой) объектов, относящихся к самым разнообразным предметным областям, обеспечивает возможность её использования значительным кругом пользователей.
УДК 681.53
А. О. Беляев
МИКРОКОНТРОЛЛЕРНЫЙ АЛГОРИТМ ОБРАБОТКИ ЭКС
Электрокардиограмма (ЭКГ) - графическое представление разности потенциалов, возникающей во время работы сердца на поверхности тела, регистрируемой аппаратом под названием электрокардиограф в процессе электрокардиографии. ЭКГ является одним из основных методов диагностики сердечно-сосудистых заболеваний. По ней можно оценить источник (так называемый водитель) ритма, регулярность сердечных сокращений, их частоту. Все это имеет большое значение для диагностики различных аритмий. По продолжительности различных интервалов и зубцов ЭКГ можно судить об изменениях сердечной проводимости.
Электрокардиограмма представляет собой периодическую последовательность кардиоциклов. В типичном кардиоцикле (рис. 1) выделяют несколько элементов: Р-волна, ря^-комплекс и Т-волна.
РГ?-интервал |<-------Н
—Н |<—(ЗРЭ-интервал
ОТ-интервал ---------------------►)
Рис. 1. Общий вид сегмента кривой ЭКГ
Отправной точкой ряда современных методик компьютерной электрокардиографии является выделение положения QRS-комплекса, которое определяется позицией своего максимума - R-зубца. Полученная последовательность позиций R-пиков используется для сегментации кардиоцикла, измерения амплитуд и длительностей отдельных его элементов, а также в методиках анализа ВСР - вариабельности сердечного ритма.
Системы реального времени, к которым относятся и программные ЭКГ-системы, осуществляющие съем, обработку и анализ ЭКГ, предъявляют повышенные требования к производительности алгоритмов. Поэтому необходим разумный компромисс между точностью метода и его производительностью. Современные аппаратные и программные средства позволяют достичь высокой чувствительно -сти без заметной потери производительности.
Описание метода. Данная методика выделения R-зубца основана на анализе фильтрованного кардиосигнала, прошедшего предварительную аналоговую и цифровую обработку. Поскольку в общем случае на снимаемый с пациента кардиосигнал могут влиять различные электрические факторы (сетевая помеха 50 Гц, импульсные помехи, постоянные наводки на кабелях отведений), то существует необходимость по возможности снизить их влияние на сигнал, используемый для анализа алгоритмом выделения R-зубца. Для удобства рассмотрения представим алгоритм в виде последовательности функционально законченных модулей. В таком случае последовательность обработки кардиосигнала будет выглядеть, как показано на рис. 2.
Рис. 2. Функциональная последовательность обработки сигнала
На вход функционального модуля поступает сигнал фц от АЦП. Поскольку
кардиосигнал расположен и в положительной, и в отрицательной полуплоскостях, то для удобства АЦ-преобразования, при помощи аналоговой цепи, его необходимо сдвинуть в положительную полуплоскость, при этом уровень изолинии поднимается обычно до середины диапазона от 0 Вт до Vref . Имея ввиду, что возможен дрейф изолинии, а также смещение уровня аналоговой земли (рис. 3), после АЦ-преобразования необходимо убрать низкочастотную (постоянную) составляющую с частотой обычно менее 1 Гц.
Как видно из рисунка изолиния отклоняется от своего нормального положения (П-^о) до минимального значения П-5от^п и максимального - П-&'отах . Размах колебаний изолинии в различных случаях может составлять от 5 до 20% от размаха полезного сигнала, возможны и большие значения, но они обычно характеризуются плохим контактом электродов с кожей пациента или неисправностями во входных аналоговых цепях.
Блок ФВЧ представляет собой стандартный фильтр верхних частот Баттер-ворта 1-го порядка с частотой среза 1 Гц.
На выходе ФВЧ имеем сигнал фЦ, форма которого представлена на рис. 4.
Блок ФНЧ выполнен в виде фильтра нижних частот Баттерворта 2-го порядка с частотой среза 20 Гц и передаточной функцией.
На выходе ФНЧ-сигнал (фЦ) изменяется относительно нулевого значения и
не содержит в себе сетевых и импульсных помех. В таком виде он пригоден для выделения QRS комплекса и R-зубца в частности.
Функциональные узлы ФВЧ и ФНЧ могут быть программно реализованы различными методами в зависимости от разрядности и вычислительных возможностей используемого процессора. Единственное условие - время выполнения функций-фильтров должно удовлетворять требованиям режима реального времени.
Функциональный модуль КФ представляет собой сложную функцию, отслеживающую крутизну фронтов (скорость изменения) входного сигнала, математически описываемую выражениями:
( -X--2 )2 (-(т -Х{_2т-2 )2
Рі = Рі-1 -
К
и = и і-і +
К
(Рі - и-і)
м
edgei = Р- - П-.
где X - входные отсчеты кардиосигнала, К и М - масштабные коэффициенты, индекс V является константой, указывающей на смещение отсчета кардиосигнала относительно 0-го элемента, нулевым считается центральный отсчет во временном окне из последних (2^+3) отсчетов, данный коэффициент зависит от частоты дискретизации сигнала, например, для частоты съема ЭКГ 200 Гц коэффициент V равен 10.
Данным функциональным узлом отслеживаются фронты интервалов QR и RS, а собственно R-зубец характеризуется резким отрицательным фронтом выходного сигнала модуля КФ ( ), как показано на рис. 5.
Отслеживая величину сигнала edgei, возможно регистрировать R-зубец на кривой кардиосигнала.
Общий алгоритм выделения R-зубца представлен на рис. 6, в него не входят цифровые фильтры Ф1 и Ф2, как было сказано выше, их реализация зависит от разрядности и вычислительной мощности используемой аппаратной платформы, данный алгоритм описывает работу только узла Ф3.
Вычисление переменной Р, производится по следующей формуле:
(тВцТ/егс-т-2 -тВи//егс-т)2 (тВцТ/егс+т-2 -тВи//егс+а)2
Pi = Pi-1 -
2
2а
где inBuffer - массив отсчетов кардиосигнала, индексы с, w являются константами, с - центральный элемент массива, w - коэффициент, зависящий от частоты дискретизации сигнала, для частоты съема ЭКГ 200 Гц, коэффициент w равен 10.
Блок фиксирования R-зубца содержит пользовательскую процедуру, в которой может выполняться измерение RR-интервалов, расчет ЧСС (частота сердечных сокращений) и статистических параметров ЭКГ для анализа ВСР (вариабельности сердечного ритма).
Подстройка граничных параметров заключается в динамической корректировке порогового значения threshold, отражающего порог фиксации R-зубца. Такая подстройка необходима для адаптации алгоритма к различным амплитудам кардиосигнала, таккак сложно и невозможно установить среднее значение размаха сигнала ЭКГ для различных пациентов, также амплитуда сигнала может изменяться и у одного того же пациента с течением времени.
Адаптивный алгоритм подстройки показан на рис. 7.
Суть алгоритма заключается в определении величины фронта входного сигнала и установки граничного значения для фиксации R-зубца в процентах от этого значения, в данном случае порог установлен на 20% от максимальной величины фронта. С течением времени максимальная зафиксированная величина фронта декрементируется до тех пор, пока не будет зафиксировано новое максимальное значение. Значение декремента в данном случае выбрано в 1% от максимальной зафиксированной величины фронта.
На рис. 8 представлен результат моделирования работы алгоритма в пакете Simulink системы Matlab.
Параметры модели:
• длительность временного окна кардиосигнала - 20 с (на графике для удобства выведен интервал всего 4 с);
• частота съема (дискретизации) кардиосигнала - 200 Герц;
• разрядность исходного кардиосигнала - 32 бита.
С
НАЧАЛО
)
1 Г
Опрос АЦП
1 г
ФВЧ
1 г
ФНЧ
1 Г
Запись текущего отсчета ЭКГ в переменную: іпБи//егі
1 г
Вычисление переменной Рі
1 г
Вычисление переменной и
1 Г
edgei = Pi - Ut
Подстройка граничных параметров
Адаптация граничных (пороговых) значений
зависимости от изменения уровня сигнала
Рис. 6. Общий алгоритм выделения R-зубца
1
Рис. 8. Кардиосигнал и сигнал на выходе блока выделения Я-зубца
Из-за небольшого временного интервала моделирования в настоящей модели не представлена работа адаптивного алгоритма подстройки.
На основании полученного на выходе алгоритма выделения R-зубца сигнала можно проводить дальнейший анализ электрокардиограммы. Сам процесс фиксации R-зубца представляет собой отслеживание фронта сигнала на выходе алгоритма. Данная методика позволяет с достаточной точностью измерять временные интервалы между R-зубцами (RR-интервалы), строить КИГ (кардиоинтервало-граммы) проводить статистический анализ ВСР (вариабельности сердечного ритма).
Конечная реализация представленного алгоритма зависит от производительности используемой вычислительной платформы и её разрядности, требований к энергопотреблению. Как уже было сказано выше, существует возможность исключить из алгоритма блоки цифровой фильтрации, перенеся их функции на аналоговые фильтры.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Рябыкина Г.В., Соболев А.В. Вариабельность ритма сердца/ Монография. -М.: Изд-во «Оверлей», 2001.
2. Хэмптон Дж. Р. Основы ЭКГ. -М.: Изд-во: «Медицинская литература», 2006.
3. Дэбни Дж. Б., Харман Т. Л. Simulink 4. Секреты мастерства/ пер. с англ. Симонова М. Л. - М.: БИНОМ. Лаборатория знаний.
4. Pan J. A real time QRS detection algorithm”, IEEE Trans. Biomed. Eng., vol. BME-32, pp. 230-236, 1985.
5. Wan H., Cammarota J.P, Akin A., and others Comparison of QRS peak detection algorithms in extracting HRV signal. Proceedings - 19th International Conference - IEEE / EMBS Oct. 30 -Nov. 2, 1997 Chicago, IL, USA.
6. Chouhan V.S., Mehta S.S. Detection of QRS Complexes in 12-lead ECG using Adaptive Quantized Threshold”/ IJCSNS International Journal of Computer Science and Network Security, VOL.8 No.1, January 2008.
7. Sahambi J. S., Tandon S. N., Bhatt R. K. P. A New Approach for On-Line ECG Characterization.
8. Moody G.B., MarkR.G. Development and evaluation of a 2-lead ECG analysis program. Massachusetts Institute of Technology, Cambridge, Mass., and Beth Israel Hospital, Boston, Mass. USA.
УДК 621.396.218
Е.Н. Остроух, М.Ю.Чегодарь
ОБУЧЕНИЕ НЕЙРОННОЙ СЕТИ С ИСПОЛЬЗОВАНИЕМ ГЕНЕТИЧЕСКИХ АЛГОРИТМОВ
Большой ряд задач, в том числе социально-экономических, задач проектирования технологических процессов, моделируется нейронными сетями, решение которых существенно упрощается при использовании эволюционных (генетических) алгоритмов. При этом получается большая точность при уменьшении времени, затраченного на обучение нейронной сети. Выбор структуры НС осуществляется в соответствии с особенностями и сложностью задачи. Наиболее популярные и изученные архитектуры - это многослойный персептрон, нейронная сеть с общей регрессией, нейронные сети Кохонена и другие. Если же задача не может быть сведена ни к одному из известных типов, разработчику приходится решать сложную проблему синтеза новой конфигурации [1].
На следующем этапе следует "обучить" выбранную нейронную сеть, т. е. подобрать такие значения ее весов, чтобы она работала нужным образом. Обучить
