Навигация
Аналіз оптико-електронних ІІС для аналізу гемодинамічних показників
67501
знак
0
таблиц
36
изображений
1.2 Аналіз оптико-електронних ІІС для аналізу гемодинамічних показників
Реанімаційно-хірургічні монітори ЮМ-300 мають вбудовану систему автоматизованого кардіо- і реоаналізу. В основу аналізу покладений метод математичної обробки плетизмограми і кардіоінтервалів, зареєстрованих протягом визначеного часу.
У результаті аналізу ритму серця будуються такі графіки:
· РИТМОГРАМА - це послідовність вертикальних ліній, висота яких відповідає тривалості відповідного RR-інтервалу (в секундах). На осі абсцис відкладаються порядкові номери RR-інтервалів;
· ГІСТОГРАМА (варіаційна пульсометрія) - східчаста функція розподілу RR-інтервалів у досліджуваному ряді їхніх значень;
· СКАТТЕРГРАМА - послідовне нанесення на графік у прямокутній системі координат двох сусідніх RR-інтервалів. Скаттерграма особливо ефективна, при діагностиці аритмії.
· Крім того, проводиться спектральний аналіз ритму серця, будується відповідний графік (спектр ритму) і визначаються частотні складові спектра.
Система кардіоаналізу дозволяє визначати наступні параметри:
ДХ =RRmax - RRmin - варіаційний розмах; Мо - мода (значення RR-інтервалу, що найчастіше зустрічається); АМо - амплітуда моди (число реалізації (у відсотках) даної моди стосовно загального числа аналізованих кардіоінтервалів); RR - середня тривалість (мс) синусових RR-інтервалів за 5 хв: SDNN - середньоквадратичне відхилення від середньої тривалості всіх синусових RR-інтервалів (за 5 хв.); Cv - коефіцієнт варіації (%) - відношення SDNN до RR; SDANN - середньоквадратичне відхилення від середніх тривалостей синусових інтервалів, розрахованих на всіх 5-хвилинних інтервалах запису; PNN50 - частка сусідніх синусових інтервалів, що розрізняються більш ніж на 50 мс; RMSSD - середньоквадратичне відхилення між тривалістю сусідніх синусових інтервалів; ULF - потужність (енергія) складової спектра ультранизької частоти (0-0,04 Гц); LF - потужність низькочастотної складової спектра (0,04-0,15 Гц); HF - потужність високочастотної складової спектра (0,15-0,4 Гц); LF/HF - відношення потужностей низькочастотної і високочастотної складових спектра.
З аналізу плетизмографічних даних визначаються наступні параметри: t1 - тривалість RR-інтервалу; t2 - анакротична фаза; t3 - катакротична фаза; t4 - період швидкого кровонаповнення; t5 - період повільного кровонаповнення; t6 - період венозного відтоку; Vsf - сфігмографічна швидкість; 11 - дикротичний індекс; 12 - діастолічний індекс; 13 - індекс периферійного опору; 14 - інтегральний гідравлічний індекс; 15 - інтегральний артеріальний індекс; 16 - інтегральний венозний індекс;
Всі накопичені дані можна передати в комп'ютер і роздрукувати за допомогою спеціального програмного забезпечення [10].
2. ВИКОРИСТАННЯ ПЕРЕТВОРЕННЯ ФУР’Є ДЛЯ АНАЛІЗУ ПУЛЬСОВОЇ ХВИЛІ
Програма обчислення миттєвого значення частоти ударів пульсу входить у склад спеціалізованого програмного забезпечення (СПЗ) та проводить обробку аналогових сигналів, що надходять з датчиків пульсової хвилі різного типу (оптоелектронних, ємнісних, тензометричних і т.д.), з метою обчислення періоду пульсової хвилі та перерахунка миттєвого значення частоти ударів пульсу за хвилину. Задача ускладнюється тим, що пульсовій хвилі, як і іншим біомедичним сигналам, що повторюються, притаманний квазіперіодичний характер [5]. Це означає, що кожний наступний період сигналу лише приблизно відповідає попередньому, особливо за амплітудою відповідних ділянок (наприклад, екстремальних значень). Крім того, форма сигналу може різко змінюватись від періоду до періоду у зв'язку із загальним хвилюванням дослідженого хворого, що був поміщений у незвичну для нього обстановку. Тому неможливо використовувати прості методи обчислення періоду сигналу, що полягають у пошуку екстремальних точок з однаковою амплітудою. Хороші результати отримують при використанні перетворення Фур'є та аналізі періоду першої гармоніки розкладеної у спектр сигналу, однак дані методи потребують значних обчислювальних витрат за часом та об'ємом оперативної пам'яті.
В результаті моделювання запропонований достатньо простий швидкодіючий засіб обчислення миттєвого значення частоти пульсу, що використовує прості операції складання, віднімання та порівняння.Алгоритм обчислення миттєвого значення частоти ударів пульсу
Запропонований алгоритм обчислення миттєвого значення частоти ударів пульсу складається з таких етапів:
- виконується настройка таймера ТО ОМЕВМ на заданий дискрет часу t у нс роботи АЦП;
- виконується установка коефіцієнта перерахунка лічильника для коректування правильної роботи системних годинників;
- виконується установка максимального розміру місця у ОЗП, що відводиться під буфер відліків з АЦП N буфера;
- виконується установка програмного флага роботи АЦП, тобто дозволяється робота АЦП, що виконується програмою обробки переривання від таймера;
- відбувається циклічний аналіз стану програмного флага переповнення таймера ТО. Якщо підпрограма обробки переривання від таймера, то ще не заповнений весь буфер відліків у ОЗП, програмний флаг роботи зберігає одиничне значення. Таким чином, відбувається зупинка роботи основної програми до повного заповнення буферів відліку, коли підпрограма обробки переривання спрацьовує програмний флаг роботи та зупиняє запис відліків у буфер;
- відбувається перегляд всього буфера підрахунків та пошук максимального за величиною відліку Аmax, наповнення суми відліків та обчислення середнього за буфером відліку
A = sum [Aі, де i=0...Nбуф-1] / Nбуф
де і - номер відліку,
Аі - і-й відлік буфера,
Nбуф - розмір буфера відліків.
Якщо розмір буфера обліку Nбуф обраний кратним степеню двійки, тобто Nбуф = 2^j, де j=1,2,3..., то операція поділу накопичення суми замінюється на просте відкидання j молодших розрядів накопиченої суми відліків за буфером;
- відбувається обчислення величини порога Ф, який дозволяє виділити характерні фрагменти буфера відліку
Ф = A + (Amax - A)/2 = (Amax+A)/2 ;
- відбувається порогова обробка вихідного масиву відліку та занулення тих відліків, у яких величина менша порога, тобто порогових відліків Аі:
Аі, якщо Аі>P
Аі' =
0, якщо Аі<=P
Дана обробка дозволяє впустити лише ті фрагменти, які містять систолічний/діастолічний піки пульсової хвилі, а виключає невірні та випадкові викидає;
- відбувається пошук найближчого від початку буфера фрагмента A1 = {A1', де i=m...n} для якого А1'<>0 та у даному фрагменту визначається координата Хмах за величиною відліку А'(x) = max [A']. Якщо існує декілька розташованих поряд та рівних за величиною екстремальних відліків, тобто є у наявності горизонтальна "поличка" хвилі, відбувається обчислення ширини полички d та корекції координати Х на половину ширини полички Х1 = Х + d/2;
- координата Х1 записується як координата першої хвилі пульсу;
- відбувається обнуління знайденого фрагмента А1, тобто А1'=0 для всіх i=m...n;
- відбувається повторний пошук найближчого від початку буфера нульового фрагмента А11 = {Ai', i=k...1}, причому використовується одна й та ж підпрограма пошуку, що й у попередньому випадку;
-у фрагменту А11, як і раніше, визначається координата Х11максимального за величиною відліку з урахуванням можливої наявності горизонтальної полички. Координата Х11 запам'ятовується як координата наступної пульсової хвилі;
- вираховується період Т пульсової хвилі шляхом визначення різниці Х11-Х1 та множення його на тривалість одного такту роботи АЦП.
T = (X11 - X1) Xt;
- оскільки тривалість такту t вимірюється в одиницях мілісекунд, то миттєва частота ударів пульсу за хвилину F визначається як
F = 60 000 / T
Оскільки результат обчислень F представлений у загальному двійковому форматі, відбувається його перетворення у двійково-десяткову форму, зручну для людського сприйняття, та вивід результату обчислень на дисплей, тобто запис кодів Семи-сегментних індикаторів у буферний ЗП контролера клавіатури та дисплея.
Але під час виконання роботи був знайдений більш ефективний метод для аналізу пульсової хвилі – вейвлет-аналіз, якому і присвячений наступний розділ.
Похожие работы
... функції, що використовується. Ця ширина називається носієм функції. Якщо вікно досить вузьке, то говорять про компактний носій. Як побачимо надалі, ця термінологія особливо широко використовується в теорії вейвлет-перетворень. Часова інформація при ПФ відсутня. При ВПФ вікно має кінцеву довжину, накриває тільки частину сигналу, тому частотне розрізнювання погіршується. Отже, чим вужче вікно, тим ...
... залежить від віконної функції, що використовується при ВПФ або материнського вейвлета при вейвлет-перетворенні. 3. Апроксимуюча і деталізуюча компоненти вейвлет-аналізу Одна з основних ідей вейвлет-подання сигналу полягає в розбивці наближення до сигналу на дві складові: грубу (апроксимуючу) і витончену (деталізуючу), з подальшим уточненням ітераційним методом. Кожен крок такого уточнення ві ...
... регулирования движения судов: Отчет о НИР (промежуточный) // ХАИ. – 501-4/2002; – Харьков, 2002. – 30 с. АНОТАЦІЯ Жеребятьєв Д.П. Методи обробки динамічних сцен при впливі нестаціонарних завад у радіотехнічних системах супроводження надводних протяжних об'єктів. – Рукопис. Дисертація на здобуття наукового ступеня кандидата технічних наук за фахом 05.12.17 – радіотехнічні та телевізійні ...
... масштабу. (г) Многомасштабне градієнтное зображення Інші автори дотримуються підходу, при якому остаточна картина границь складається на основі аналізу градієнтних зображень від точних масштабів до не точних. При цьому, основними завданнями при такому підході є зменшення впливу шуму, до якого чутливі оператори градієнта малого розміру, і комбінування границь, отриманих на точних масштабах, із ...
0 комментариев