Спосіб оцінювання електричної альтернації серця

Реферат: Винахід належить до галузі медицини і стосується способу оцінювання електричної альтернації серця, заснований на аналізі порядку чергування елементів електрокардіограми з різними значеннями параметрів, що характеризують їх форму, зокрема, за амплітудою, тривалістю в часі та симетрією зубця Т. У способі оцінюють розподіл значень визначених параметрів на UA 101918 C2 (12) UA 101918 C2 вибраному інтервалі у вигляді гістограми з трьома градаціями значень, за формою якої приймають попереднє рішення про наявність електричної альтернації серця, а остаточне рішення про наявність електричної альтернації серця та оцінку рівня альтернації приймають на основі аналізу порядку чергування циклів двох типів, причому допускається, що в деякі моменти часу фаза чергування циклів змінюється. UA 101918 C2 5 10 15 20 25 30 35 40 45 50 55 60 Винахід належить до медицини, а більш точно до способів оцінювання електричної альтернації серця за сигналом, що несе інформацію про електричну активність серця (електрокардіограмою або магнітокардіограмою), та може бути застосована в медичних діагностичних системах для визначення ризику раптової серцевої смерті. За даними медичної статистики майже 50 % випадків смерті від серцево-судинних захворювань припадає на раптову серцеву смерть, коли хворий гине практично миттєво (від кількох секунд до години) після початку серцевого нападу [див. Zipes D. P., Wellens H. J. Sudden Cardiac death // Circulation. - 1998. - Vol. 98. - P. 2334-2351]. Тому визначення ризику раптової серцевої смерті є важливою характеристикою сучасних комп'ютерних систем медичної діагностики. Один з відомих предикторів раптової серцевої смерті, якій останнім часом став популярним в медичній практиці, є електрична альтернація серця, яка визначається за характерним чергуванням форми інформативних елементів електрокардіограми (ЕКГ), що відповідають періодам реполяризації шлуночків, зокрема чергування Т-зубців ЕКГ з високими та малими рівнями амплітуд [див. Liew R. Electrocardiogram-Based Predictors of Sudden Cardiac Death in Patients with Coronary Artery Disease // Clinical Cardiology.-2011. - Vol. 34. - Issue 8. - P. 466-473]. Відомо спосіб оцінювання електричної альтернації серця, який ґрунтується на аналізі морфології зубців Τ в частотній області з використанням алгоритму швидкого перетворення Фур'є [див. Rosenbaum D.S., Jackson L.E., Smith J.M. et al. Electrical alternans and vulnerability to ventricular arrhythmias // New England Journal of Medicine. - 1994. - No. 330. - P. 235-241]. Згідно з цим способом вважається, що електрична альтернація серця є на певному часовому інтервалі спостереження сигналу, якщо на усередненому спектрі, побудованому для послідовності серцевих циклів, що належать цьому інтервалу, з'являється характерний пік при частоті, яка дорівнює половині середньої частоти серцевих скорочень. Але цей спосіб передбачає, що під час спостереження у пацієнта за допомогою спеціальних засобів підтримується стабільна частота серцевих скорочень на високому рівні (не менш ніж 100 ударів в хвилину), що є небезпечним для певної групи пацієнтів. До того ж алгоритм швидкого перетворення Фур'є дає прийнятні результати лише, коли альтернація елементів ЕКГ змінюється за синусоїдальним законом. Відомо також інший спектральний спосіб оцінювання електричної альтернації серця [див. Патент Российской Федерации № 2371087 / Мельник О.В. Способ выявления альтернаций Тзубца электрокардиосигнала в режиме реального времени и устройство для его осуществления, RU 2008111542/14, заявлен 25.03.2008, опубликован 27.10.2009]. Спосіб передбачає розклад інтервалу реполяризації сигналу за базисними функціями Уолша з наступним аналізом значень коефіцієнтів такого розкладу. Але, як зазначає сам автор методу [див. Мельник О.В. Метод выявления альтернации Тзубца электрокардиосигнала // Биомедицинская радиоэлектроника, 2008. - № 7. - С. 56-59], цей метод дуже чутливий до точності визначення періоду реполяризації. Навіть незначна помилка у визначенні точок початку і закінчення зубця Τ може призвести до неприпустимої похибки в оцінці коефіцієнтів розкладу за функціями Уолша, що звужує область практичного використання методу в реальних умовах. Альтернативні способи аналізу електричної альтернації серця, на відміну від спектральних, засновані на аналізі альтернації Т-зубців ЕКГ безпосередньо в часовій області. Серед відомих способів аналізу електричної альтернації серця, заснованих на аналізі альтернації Т-зубців ЕКГ безпосередньо в часовій області, найбільш близьким до способу, що пропонується, є спосіб який описано [в патенті США № 6,169,919 / Nearing B.D., Verrier R.L. System and method quantifying alterbation in an electrocardiogram signal, US Appl. No. 09/335,680, File Date 18 June, 1999, Issue Date: Jan 2 Jan., 2001]. Цей спосіб передбачає, що вимірюють сигнал, який несе інформацію про електричну активність серця на вибраному інтервалі спостереження (ЕКГ), на послідовності серцевих циклів визначають елементи сигналу, що відповідають періодам реполяризації шлуночків (зубці Т), визначають амплітуди визначених зубців Τ окремо для парних та непарних серцевих циклів, визначають середнє значення амплітуд для парних та непарних циклів та оцінюють рівень електричної альтернації як різницю між визначеними середніми значеннями амплітуд окремо для парних та непарних циклів. Основним недоліком цього способу є те, що він розрахований лише на ідеалізований випадок, коли на парних циклах спостерігаються тільки зубці одного типу (умовно названим типом А), а на непарних - зубці другого типу (умовно названим типом В). Іншими словами, коли послідовність циклів може бути умовно позначена "словом" А В А В А В, в якому чітко чергуються "літери" А та В. 1 UA 101918 C2 5 10 15 20 25 30 35 40 45 50 55 60 У той же час клінічна практика показує, що при електричній альтернації фаза чергування зубців типу А та В у випадкові моменти часу може змінюватися, інакше кажучи, зубці одного типу можуть з'являтися як на парних, так і на непарних циклах [див. Nearing B.D., Verrier R.L. Modified moving average analysis of T-wave alternans to predict vetricular fibrilation with high accuracy // Journal of Applied Physiology. - 2002 - Vol. 92 - P. 541-549]. В цьому випадку морфологія ЕКГ описується "словом", в якому на сусідніх циклах з'являється одна і та ж "літера", наприклад, "словом" A B B А В А В …. В таких випадках результат, отриманий за відомим способом аналізу електричної альтернації серця, не є достовірним. Другим недоліком відомого способу є те, що він спрямований на аналіз альтернації зубців Τ лише за амплітудою, що не охоплює усіх практично можливих випадків альтернації форми цього елемента ЕКГ. В основу винаходу поставлена задача створення способу, який забезпечує підвищення достовірності та точності результату аналізу електричної альтернації серця в реальних умовах застосування. Спосіб, що пропонується, полягає в тому, що вимірюють сигнал, який несе інформацію про електричну активність серця на вибраному інтервалі спостереження, на послідовності серцевих циклів визначають елементи сигналу, що відповідають періодам реполяризації шлуночків (зубці Т), визначають значення параметрів, що характеризують форму визначених елементів сигналу на кожному з серцевих циклів, оцінюють розподіл значень визначених параметрів на вибраному інтервалі часу шляхом визначення відсотків попадання значень параметрів в кожен з трьох рівних інтервалів діапазону значень, приймають попереднє рішення про наявність електричної альтернації серця, якщо різниці між відсотками попадання значень параметрів в перший та третій інтервали відносно до попадання в другий інтервал перевищують заданий поріг, розподіляють усі цикли на два типи, до першого з яких відносять цикли, на яких визначені параметри перевищують медіану параметра, визначену на вибраному інтервалі часу, а до другого типу відносять цикли, на яких визначені параметри не перевищують визначену медіану, приймають остаточне рішення про наявність електричної альтернації на вибраному інтервалі, якщо цикли першого та другого типу чітко чергуються в часі або в деякі моменти часу фаза чергування циклів змінюється, але число таких моментів часу не перевищує задану величину. При цьому як параметри, що характеризують форму елементів сигналу, визначають амплітуди та (або) тривалості в часі та (або) параметри, що характеризують симетрію відповідних елементів на окремих циклах. Сукупність ознак, що відрізняють цей спосіб від відомих (а саме те, що додатково оцінюють розподіл значень визначених параметрів на вибраному інтервалі часу шляхом визначення відсотків попадання значень параметрів в кожен з трьох рівних інтервалів діапазону значень, приймають попереду рішення про наявність електричної альтернації серця, якщо різниці між відсотками попадання значень параметрів в перший та третій інтервали відносно до попадання в другий інтервал перевищують заданий поріг, розподіляють усі цикли на два типи, до першого з яких відносять цикли, на яких визначені параметри перевищують медіану параметра, визначену на вибраному інтервалі часу, а до другого типу відносять цикли, на яких визначені параметри не перевищують визначену медіану, та приймають остаточне рішення про наявність електричної альтернації на вибраному інтервалі, якщо цикли першого та другого типу чітко чергуються в часі або в деякі моменти часу фаза чергування циклів змінюється, але число таких моментів часу не перевищує задану величину, а як параметри, що характеризують форму елементів сигналу, що відповідають періодам реполяризації шлуночків, визначають амплітуди відповідних елементів та (або) тривалості в часі відповідних елементів та (або) параметри, що характеризують симетрію відповідних елементів на окремих циклах) разом з відомими ознаками дозволяє суттєво підвищити достовірність та точність аналізу електричної альтернації cepця в реальних умовах практичного застосування. На фіг. 1 наведено схему, яка ілюструє послідовність операцій, що забезпечують реалізацію способу. На фіг. 1 наведено такі позначення: 1 - блок реєстрації ЕКГ та аналого-цифрового перетворення; 2 - блок визначення параметрів форми Т-зубців ЕКГ; 3 - блок аналізу гістограми розподілу параметрів; 4 - блок поділу циклів на два типи; 5 - блок аналізу чергування циклів двох типів; 6 - блок визначення рівня альтернації; На фіг. 3 наведено фрагмент ЕКГ, яка несе інформацію про електричну активність серця. Показано три серцеві цикли, на кожному з яких позначено: 7 - зубець Р, який породжується електричним збудженням передсердя; 2 UA 101918 C2 5 10 15 20 25 30 35 40 45 50 55 8 - комплекс QRS, який відображає процес деполяризації шлуночків серця; 9 - зубець Т, форма якого відображує процес реполяризації шлуночків. На фіг. 3 наведено приклади позитивного (фіг. 3,а) і негативного (фіг. 3,б) зубців Т, на яких позначено: 10 - момент початку зубця Т; 11 - рівень ізоелектричної лінії відповідного циклу ЕКГ; 12 - максимального значення амплітуди зубця Τ при позитивному зубці; 13 - момент часу, коли досягається максимальне значення зубця Τ при позитивному зубці) або мінімальне значення зубця Τ (при негативному зубці); 14 - момент закінчення зубця Т. 15 - мінімальне значення амплітуди зубця Τ при негативному зубці. t1 - тривалість від початку 10 зубця Τ до моменту 13 часу, коли досягається максимальне значення зубця Τ (при позитивному зубці) або мінімальне значення зубця Τ (при негативному зубці); t 2 - тривалість від моменту 13 часу, коли досягається максимальне значення зубця Τ (при позитивному зубці) або мінімальне значення зубця Τ (при негативному зубці) до моменту 14 закінчення зубця Т. На фіг. 4 наведено приклади ЕКГ, зокрема - ЕКГ без альтернації зубця Τ (фіг. 4,а); - ЕКГ з чітким чергуванням амплітуд зубців Τ (фіг. 4,б); - ЕКГ зі зміною фази чергуванням амплітуд зубців Τ (фіг. 4,в); - ЕКГ зі зміною амплітуди зубців Τ в серединні інтервалу спостереження (фіг. 4,г). На фіг. 5 наведено приклади гістограм розподілу амплітуд Т-зубців ЕКГ при відсутності альтернації Т-зубців (фіг. 5,а) та при альтернації Т-зубців (фіг. 5,б). На фіг 6 наведені приклади ЕКГ та значення величин, що визначаються блоком 5 аналізу чергування циклів двох типів, а саме - ЕКГ з чітким чергуванням амплітуд зубців Τ (фіг 6,а); - ЕКГ зі зміною фази чергування на восьмому циклі (фіг 6,б); - ЕКГ, на якій двічі спостерігаються три послідовних цикли одного типу (фіг 6,в). Фіг. 7 ілюструє приклади ЕКГ, на яких спостерігається чергування позитивних та негативних зубців Τ (фіг. 7,а) та чергування зубців Τ з визначеною та "невизначеною" амплітудою (фіг 7,б). На фіг. 8 наведено приклади ЕКГ, на яких спостерігається чергування зубців Τ з різною тривалістю (фіг. 8,а) та ЕКГ, на яких відбувається чергування симетрії зубців Τ (фіг. 8,б). Спосіб оцінювання електричної альтернації серця, що пропонується, складається з таких дій. Блок 1 (фіг. 1) забезпечує реєстрацію сигналу, що несе інформацію про електричну активність серця, зокрема ЕКГ та перетворення в цифрову форму ЕКГ. Цифровий сигнал з виходу блока 1 надходить на вхід блока 2, який забезпечується автоматичне визначення параметрів, що характеризують форму зубців Τ на кожному серцевому циклі ЕКГ. Розглянемо спочатку деталі способу, що пропонуються, для випадку, коли як параметр, який характеризує форму зубця Т, використовується амплітуда цього зубця. За допомогою блока 2 на певному інтервалі спостереження визначається масив значень Y1, Y2,, YN амплітуд зубців Т кожного з N циклів ЕКГ, що належать цьому інтервалу. Амплітуда зубця Τ визначається як відхилення максимального значення 12 позитивного елемента (фіг. 3,а) відносно до рівня ізоелектричної лінії 11 відповідного циклу або відхилення мінімального значення 15 негативного елемента (фіг. 3, б) відносно до рівня ізоелектричної лінії 11 відповідного циклу. Слід зауважити, що в реальних клінічних умовах, коли на сигнал діють різного роду збурення, на окремих циклах спостерігаються певні зміни амплітуд зубців Τ, навіть при відсутності альтернації (фіг. 4,а). В цьому випадку амплітуди Y1, Y2,, YN можна розглядати як реалізації однієї випадкової величини, яка з обмеженою дисперсією варіює відносно деякого невідомого середнього значення. Коли ж спостерігається альтернація (фіг. 4,б), то амплітуди YA і YB зубців Τ на циклах різного типу будемо розглядати як реалізації двох випадкових величин, які з обмеженою дисперсією варіюють відносно "своїх" заздалегідь невідомих середніх значень, які нерівні один одному. В цьому випадку альтернацію зубців Τ буде характеризувати два непересічних інтервали амплітуд, коли 3 UA 101918 C2     min max min max YA  YA , YA , YB  YB , YB , , (1) та max min max min YA  YB або YB  YA . 5 10 15 (2) Будемо також розглядати загальний випадок, коли при альтернації порядок чергування зубців двох типів в випадкові моменти часу може бути порушений, але число таких моментів менше, ніж деяка задана величина. Приклад такої ЕКГ показано на фіг. 4, в, де на шостому циклі змінюється порядок чергування зубців: амплітуди зубців Τ на п'ятому та шостому циклах однакові. Враховуючи ці особливості, які спостерігаються на ЕКГ в реальних клінічних умовах, прийняття рішення про наявність або відсутність альтернації амплітуд зубців Т, згідно зі способом, що пропонується, приймається в два етапи за такою послідовністю операцій. Масив значень Y1, Y2,, YN надходить на вхід блока 3 (фіг. 1), за допомогою якого оцінюється розподіл цих параметрів у вигляді гістограми, що будується наступним чином. Спочатку за масивом Y1, Y2,, YN зубців Τ визначається впорядкований масив ~ ~ ~ ~ Y N  Y1, Y2,, YN , (4) наприклад, за зростанням амплітуд зубців Т. Далі визначається діапазон зміни амплітуд зубців Τ: 20 25 30 35 40 Q  max Y  min Y, i  1 ,N , i i , (5) ~ ~ де min Y  Y1 та max Y  YN . i i Наступним кроком є розбиття діапазону Q на M  3 рівних інтервали групування значень та ~ за впорядкованим масивом Y N визначається відсоток попадання елементів вихідного масиву Y1, Y2,, YN в кожен з цих трьох інтервалів. Числа q1, q2 , q3 визначають висоту відповідного стовпчика гістограми. При цьому, якщо будь який з елементів масиву Y1, Y2,, YN попадає на границю інтервалу групування, то при визначенні проценту попадання умовно вважається, що "половина" елемента належить лівому інтервалу, а друга "половина" - правому. Для ілюстрації на фіг. 5 показано приклади двох гістограм, побудованих згідно з наведеною вище послідовністю операцій. Перша гістограма (фіг 5, а) відповідає випадку, коли на інтервалі спостережень не відбувається альтернація зубців Τ (див. приклад на фіг. 4 а), а лише спостерігаються деякі зміни амплітуд відносно середнього рівня. Друга гістограма (фіг 5, б) відповідає випадку, коли на інтервалі спостережень відбувається альтернація зубців Τ (фіг 4, б), тобто зубці першого та другого типу варіюють відносно своїх середніх значень, які нерівні один одному. В цьому випадку на гістограмі відбувається "провал" в середньому інтервалі групування (фіг 5, б). Таким чином, про наявність альтернації зубців Τ можна судити на основі порівняння трьох рівнів гістограм q1, q2 , q3 і попереднє рішення про наявність альтернації приймається в тому випадку, коли виконуються співвідношення q1  q2  q0 та q3  q2  q0 , (6) де 45 q0 - параметр, що настроюється за експериментальними даними, наприклад встановлюється q0  20% . Слід зауважити, що умови (6) можуть виконуватись і в деяких екзотичних випадках, коли, наприклад, на перших N1  N / 2 циклах амплітуда зубців Τ варіює відносно деякого одного значення, а на інших циклах - відносно деякого другого значення (див. приклад на фіг. 5, г). Можна навести і деякі інші екзотичні приклади ЕКГ, які породжують гістограму, схожу на ту, що наведено на фіг 5, 6. 4 UA 101918 C2 5 Тому для прийняття надійного рішення про наявність або відсутність альтернації зубців Τ недостатньо лише перевірки співвідношень (6) і остаточне рішення приймається на другому етапі наступним чином. При виконанні умов (6) за допомогою блока 3 (фіг. 1) проводиться додатковий аналіз ~ порядку чергування циклів. Для цього за впорядкованим масивом Y N визначається порогове значення Y0 як медіана вхідного масиву Y1, Y2,, YN : ~N  Y[ ], якщо N  парне,  ,(7) Y0   2 ~ Y[(N div 2)  1], якщо N  непарне  10 де div - операція цілочисельного ділення. Порогове значення Y0 разом з масивом Y1, Y2,, YN надходять на вхід блока 4, який виконує операцію поділу циклів на два типи (тип А і В) та визначення масивів відповідних номерів циклів I A та IB . Операція поділу циклів на два типи здійснюється на основі співвідношень ТИП A, якщо Y  Y0 ; ТИП B, якщо Y  Y0 ;i  12,,N , i i , 15 20 25 (8) Отримані масиви I A та IB номерів циклів разом з масивом Y1, Y2,, YN надходять на вхід блока 5. Цей блок здійснює аналіз порядку чергувань циклів двох типів та прийняття остаточного рішення про наявність або відсутність альтернації зубців Τ на вибраному інтервалі спостереження. Альтернація підтверджується, коли на послідовності номерів i  12,,N циклів ЕКГ , спостерігається чітке чергування амплітуд зубців типів А та В або спостерігається не більш ніж задане число  0 змін фаз чергування зубців. В іншому випадку, навіть коли виконуються співвідношення (6), приймається остаточне рішення про відсутність альтернації. Для виконання такої перевірки за допомогою блока 5 визначається чи знайдуться на вибраному інтервалі спостереження три послідовних цикли ЕКГ одного типу (А або В). Така перевірка здійснюється на основі визначення двох величин , Z i  1  1 якщо i  I A ; , Z A    A i 0, якщо i  IB  (9) та 30 35 Z i  1  1, якщо i  IB ; , Z B    B i 0, якщо i  I A .  (10) Якщо задати початкові умови Z A 0  0 і ZB 0  0 , то при чіткому чергуванні циклів значення величин Z A  і ZB  не перевищують одиниці (фіг. 6,а). i i При зміні фази чергування, коли сусідні цикли одного типу (А або В), значення відповідної величини ( Z A  або ZB  ) буде дорівнювати двом (фіг. 6,б). i i Якщо ж на ЕКГ послідовно зустрічаються три або більше циклів одного типу, то значення однієї з величин Z A  або ZB  буде перевищувати два (фіг 6,в). i i Таким чином, виконання однієї з умов i Z A   2 або ZB   2 , i (11) 40 свідчить про відсутність альтернації зубців Τ, навіть при виконанні умови (6). Для остаточного прийняття рішення за допомогою блока 5 (фіг. 1) здійснюється також додаткова перевірка числа змін фаз чергування зубців. Для цього визначається величина   i  1  1 якщо Z A   2 або ZB   2, i  12,,N , i , i i , 5 (12) UA 101918 C2 при початковій умові 0  0 . Про відсутність альтернації свідчить виконання умови   0 , (13) 5 10 де  0 - задана величина, яку доцільно зв'язати з числом N циклів, що аналізуються, наприклад, визначив  0  0,1N . Таким чином, остаточне рішення про наявність альтернації зубців Τ на вибраному інтервалі спостереження приймається лише в тому випадку, коли виконуються умови (6), але не виконується жодна з умов (11) та (13). Якщо остаточне рішення про наявність альтернації прийнято, то масив Y1, Y2,, YN разом за масивами номерів циклів I A та IB надходять до блока 6 визначення рівня альтернації. Рівень альтернації визначається за співвідношенням  Y  YA  YB , (14) 15 де YA та YB - відповідно середні значення амплітуд зубців Τ циклів типу А і В на вибраному інтервалі спостереження, які визначаються за формулами YA  20 25 30 i  Y, iIA YB  40 (15) Співвідношення (14) дозволяє визначати рівень альтернації в загальному випадку, навіть коли при альтернації зубців Τ спостерігається чергування позитивних та негативних зубців Τ (фіг. 7,а) або чергування зубців Τ з визначеною та "невизначеною" амплітудою внаслідок низького рівня амплітуди (фіг. 7,б). Спосіб, який пропонується, без будь яких обмежень може бути застосований для аналізу альтернації зубців Τ не тільки за амплітудою, а й за іншими параметрами, що характеризують форму цих зубців. Наприклад, для визначення рівня альтернації зубців Τ з різною тривалістю в часі (фіг. 8,а) за допомогою блока 2 (фіг. 1) визначається масив 1, 2,, N тривалостей зубців Τ на кожному з N циклів інтервалу спостереження. Тривалість ,i  1 ,N зубця Τ на i-му циклі визначається i , за формулою   t E   t B  , i i i 35 i  Y , iIB (16) де t B , t E  - відповідно моменти 10 та 14 (фіг. 3) початку та закінчення зубця Τ на і-му циклі, i i які визначаються за моментами перетину сигналу рівня ізоелектричної лінії 11. Далі, згідно зі способом що пропонується, за допомогою блока 3 (фіг. 1) будується гістограма розподілу значень 1, 2,, N та перевіряється співвідношення (6) для прийняття попереднього рішення про наявність альтернації тривалостей зубців Τ на вибраному інтервалі спостереження. При виконанні умов (6) аналогічно тому, як це було описано вище, за допомогою блока 3 визначається порогове значення  0 як медіана масиву 1, 2,, N та за допомогою блока 4 здійснюється поділ циклів на два типи (тип А і В) та визначення масивів відповідних номерів циклів I A та IB за правилом, аналогічним правилу (8), а саме правилом: ТИП A, якщо   0 ; ТИП B, якщо   0 ;i  12,,N , i i , 45 (17) Масив 1, 2,, N разом з масивами номерів I A та IB , які визначені блоком 4, надходять до входу блока 5, який здійснює додаткову перевірку послідовностей чергування на основі співвідношень (9)-(13). 6 UA 101918 C2 Якщо за допомогою такої перевірки приймається остаточне рішення про наявність альтернації, то масив 1, 2,, N значень тривалостей зубців Τ разом з масивами номерів циклів I A та IB надходить до блока 6, за допомогою якою визначається рівень альтернації   таким чином: 5     A  B , (18) де  A та B - відповідно середні значення тривалостей зубців Τ циклів типу А і В на вибраному інтервалі спостереження, які визначаються за формулами A  i  , iIA B  i   , iIB (19) 10 15 Аналогічним чином здійснюється визначення рівня електричної альтернації серця за симетрією зубців Τ (фіг. 8,б). Для цього за допомогою блока 2 (фіг. 1) визначається масив 1,2,,N , елементи якого відповідають параметрам, що характеризують симетрії зубців Τ на кожному з N циклів інтервалу спостереження. Параметр ,i  1 ,N зубця Τ на і-му циклі згідно зі способом, що пропонується, i , визначається за формулою   i 20 25 30 t1 , i t 2  i (20) де t1 - тривалість від початку 10 зубця Τ (фіг. 3) до моменту 13 часу, коли досягається i максимальне значення 12 зубця Τ при позитивному зубці або мінімальне значення 15 зубця Τ при негативному зубці; t 2  - тривалість від моменту 13 часу, коли досягається максимальне значення 12 зубця Τ i при позитивному зубці або мінімальне значення 15 зубця Τ при негативному зубці до моменту 14 закінчення зубця Т. Далі, згідно зі способом, що пропонується, за допомогою блока 3 (фіг. 1) будується гістограма розподілу значень 1,2,,N та перевіряється співвідношення (6) для прийняття попереднього рішення про наявність альтернації симетрії зубців Τ на вибраному інтервалі спостереження. При виконанні умов (6) аналогічно тому, як це було описано вище, за допомогою блока 3 визначається порогове значення  0 як медіана масиву 1,2,,N , за допомогою блока 4 здійснюється поділ циклів на два типи (тип А і В) та визначення масивів відповідних номерів циклів I A та IB за правилом, аналогічним правилу (8), а саме правилом: ТИП A, якщо   0 ; ТИП B, якщо   0 ;i  12,,N , i i , 35 40 (21) Масив 1,2,,N разом з масивами номерів I A та IB , які визначені блоком 4, надходять до входу блока 5, який здійснює додаткову перевірку послідовностей чергування на основі співвідношень (9)-(13). Якщо за допомогою такої перевірки приймається остаточне рішення про наявність альтернації, то масив 1,2,,N значень параметрів симетрії зубців Τ разом з масивами номерів циклів I A та IB надходить до блока 6, за допомогою якою визначається рівень альтернації   таким чином:   A  B , 45 (22) де A та B - відповідно середні значення параметрів симетрії зубців Τ циклів типу А і В на вибраному інтервалі спостереження, які визначаються за формулами 7 UA 101918 C2 A  5 i  , iIA B  i   . iIB (23) Наведений вище детальний опис конкретної реалізація способу аналізу електричної альтернації серця дано лише з метою ілюстрації. Зрозуміло, що на практиці в реалізацію способу можуть бути внести деякі зміни і модифікації. Наприклад, як параметри, що характеризують форму зубців Τ можуть бути застосовані інші параметри, а замість зубця Τ застосовуватись інші інформативні елементи ЕКГ, наприклад S-Т-сегмент. Проте вважається, що такі або інші зміни і модифікації, які можуть бути виконані без суттєвих відхилень від даного способу, підпадають під дію даного винаходу. 10 ФОРМУЛА ВИНАХОДУ 15 20 25 30 35 40 45 50 1. Спосіб оцінювання електричної альтернації серця, у якому вимірюють сигнал, який несе інформацію про електричну активність серця на вибраному інтервалі спостереження, на послідовності серцевих циклів визначають елементи сигналу, що відповідають періодам реполяризації шлуночків, визначають значення параметрів, що характеризують форму визначених елементів сигналу на кожному з серцевих циклів, та оцінюють рівень електричної альтернації як різницю між середніми значеннями параметрів визначених елементів сигналу для двох типів циклів, який відрізняється тим, що оцінюють розподіл значень визначених параметрів на вибраному інтервалі часу шляхом визначення відсотків попадання значень параметрів в кожен з трьох рівних інтервалів діапазону значень, приймають попереднє рішення про наявність електричної альтернації серця, якщо різниці між відсотками попадання значень параметрів в перший та третій інтервали відносно до попадання в другий інтервал перевищують заданий поріг, розподіляють усі цикли на два типи, до першого з яких відносять цикли, на яких визначені параметри перевищують медіану параметра, визначену на вибраному інтервалі часу, а до другого типу відносять цикли, на яких визначені параметри не перевищують визначену медіану, та приймають остаточне рішення про наявність електричної альтернації на вибраному інтервалі, якщо цикли першого та другого типу чітко чергуються в часі або в деякі моменти часу фаза чергування циклів змінюється, але число таких моментів часу не перевищує задану величину. 2. Спосіб оцінювання електричної альтернації серця за п. 1, який відрізняється тим, що як параметри, що характеризують форму елементів сигналу, що відповідають періодам реполяризації шлуночків, визначають амплітуди відповідних елементів як відхилення максимального значення позитивного елемента відносно до рівня ізоелектричної лінії відповідного циклу або відхилення мінімального значення негативного елемента відносно до рівня ізоелектричної лінії відповідного циклу. 3. Спосіб оцінювання електричної альтернації серця за п. 1, який відрізняється тим, що як параметри, що характеризують форму елементів сигналу, що відповідають періодам реполяризації шлуночків, визначають тривалості в часі відповідних елементів на окремих циклах. 4. Спосіб оцінювання електричної альтернації серця за п. 1, який відрізняється тим, що як параметри, що характеризують форму елементів сигналу, що відповідають періодам реполяризації шлуночків, визначають параметри, що характеризують симетрію відповідних елементів на окремих циклах, які визначаються за формулою: t  1 , t 2 де t1 - тривалість від моменту часу початку визначеного елемента до моменту часу, коли досягається його максимальне значення при позитивному елементі або мінімальне значення при негативному елементі, а t 2 - тривалість від моменту часу, коли досягається максимальне значення визначеного елемента при позитивному елементі або мінімальне значення визначеного елемента при негативному елементі до моменту часу закінчення визначеного елемента. 8 UA 101918 C2 9 UA 101918 C2 10 UA 101918 C2 11 UA 101918 C2 Комп’ютерна верстка Л. Ціхановська Державна служба інтелектуальної власності України, вул. Урицького, 45, м. Київ, МСП, 03680, Україна ДП “Український інститут промислової власності”, вул. Глазунова, 1, м. Київ – 42, 01601 12