Лаба 10 / Лаба_10

Цель работы: ознакомление с методами ритмокардиографии, построение графиков ритмограммы, гистограммы и скаттерограммы.

Основные теоретические положения:

Ритмокардиография – это метод оценки ритма сердца, основанный на графическом представлении последовательности значений длительности каридоцикла (RR-интервалов). Используются три вида графиков:

•Ритмограмма – зависимость значений длительности RR-интервала от порядкового номера цикла измерения;

•Гистограмма – относительное число RR-интервалов, относящихся к различным диапазонам значений их длительности;

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

Для количественной оценки вариабельности сердечного ритма рассчитывается некоторые статистические показатели, в частности:

•Среднее значение

•Минимальное и максимальное значения

•Вариационный размах

•Мода

•Амплитуда моды

2

Обработка результатов эксперимента

close all clear all clc

global hMenu1 hA2 hA3 hA4 RR3

H = figure ('position', [250 50 1000 700]);

hMenu1 = uicontrol ('Style', 'popupmenu','position', [830 615 100 60]); hA2 = axes('Units', 'pixels','position', [50 420 900 220]);

hA3 = axes('Units', 'pixels','position', [50 60 420 320]); hA4 = axes('Units', 'pixels','position', [530 60 420 320]); set(hMenu1,'Callback','part2')

set(hMenu1, 'String', {'Выбрать сигнал', 'Нормальный ритм'...

'Желудочковая экстрасистолия', 'Мерцательная аритмия'}); RR3 = load('R10_04.txt');

function part2()

global hMenu1 hA2 hA3 hA4 RR3 RR

r = get(hMenu1,'Value'); if r == 1

axes(hA2) cla axes(hA3) cla axes(hA4) cla

else

RR = RR3(:, r-1); axes(hA2)

cla hold on

for i = 1:450 x(1) = i; x(2) = i; R(1) = 0; R(2) = RR(i); plot(x, R); hold on

end

maxRR = max(RR)*1.2; set(hA2, 'Ylim',[0 maxRR])

axes(hA3) dH = 0.05;

X = 0:dH:maxRR; H = histc(RR,X) SH = sum(H); PH= H/SH*100;

bar(X,PH, 'histc') maxh = max(PH)*1.2;

set(hA3, 'Ylim',[0 maxh])

axes(hA4)

NRR = length(RR);

plot( RR(1:NRR - 1), RR(2:NRR), 'o'); maxr = max(RR)*1.2;

set(hA4, 'Ylim',[0 maxRR],'Xlim',[0 maxRR])

RRsr = mean(RR); RRmin = min(RR); RRmax = max(RR);

dRR = RRmax - RRmin; [AMo, iMo] = max(PH); Mo = iMo*dH; axes(hA3)

xt = maxRR/20; yt = maxh;

dyt = maxh/15;

text (xt, yt - dyt, ['RRsr = ' num2str(RRsr) ' c']) text (xt, yt - 2*dyt, ['RRmin = ' num2str(RRmin) ' c']) text (xt, yt - 3*dyt, ['RRmax = ' num2str(RRmax) ' c']) text (xt, yt - 4*dyt, ['dRR = ' num2str(dRR) ' c']) text (xt, yt - 5*dyt, ['Mo = ' num2str(Mo) ' c'])

end return end

3

Рисунок 1 – Окно программы при выборе «нормального ритма»

Рисунок 2 – Окно программы при выборе «Желудочковая экстрасистолия»

4

Выводы:

Исходные данные содержали три столбца значений RR-интервалов: в

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

гистограммы и скаттерограммы рассчитали верхний предел Ylim как 1,2 от максимального RR-интервала. В окне hAxes3 отображена нормированная гистограмма RR-интервалов, функция PH=H/SH*100; позволяет получить гистограмму в процентах. Рассчитали статистические параметры вариабельности сердечного ритма, описанные в разделе «Основные положения»: средний RR-интервал, минимальный RR-интервал,

максимальный RR-интервал, вариационный размах, амплитуду моды (aMo) и

индекс моды (iMo), моду (Mo).

6