МИНОБРНАУКИ РОССИИ
Санкт-Петербургский государственный
электротехнический университет
«ЛЭТИ» им. В.И. Ульянова (Ленина)
Кафедра БТС
отчет
по лабораторной работе №10
по дисциплине «ПСОБМД»
Тема: Методы ритмокардиографии
Вариант 4
Студентка гр. 7501 |
|
Фалевская А.А. |
Студент гр. 7501 |
|
Исаков А.О. |
Преподаватель |
|
Ровша М.В. |
Санкт-Петербург
2020
Цель работы: ознакомление с методами ритмокардиографии, построение графиков ритмограммы, гистограммы и скаттерограммы.
Основные теоретические положения:
Ритмокардиография – это метод оценки ритма сердца, основанный на графическом представлении последовательности значений длительности каридоцикла (RR-интервалов). Используются три вида графиков:
Ритмограмма – зависимость значений длительности RR-интервала от порядкового номера цикла измерения;
Гистограмма – относительное число RR-интервалов, относящихся к различным диапазонам значений их длительности;
Скаттерограмма – двумерное отображение ритма сердца, которое строится как совокупность точек, координаты каждой из которых на плоскости соответствуют величинам двух смежных RR-интервалов.
Для количественной оценки вариабельности сердечного ритма рассчитывается некоторые статистические показатели, в частности:
Среднее значение
Минимальное и максимальное значения
Вариационный размах
Мода
Амплитуда моды
Обработка результатов эксперимента
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
Рисунок 1 – Окно программы при выборе «нормального ритма»
Рисунок 2 – Окно программы при выборе «Желудочковая экстрасистолия»
Рисунок 3 – Окно программы при выборе «Мерцательная аритмия»
Выводы:
Исходные данные содержали три столбца значений RR-интервалов: в норме (первый столбец), при наличии желудочковых экстрасистол (второй столбец), и при мерцательной аритмии (третий столбец). Каждый из пунктов hMenu1 соответствует своему столбцу данных. Процедура построения графиков оформлена в виде специальной функции, размещенной в отдельном файле. Задача данной функции построение ритмограммы, гистограммы и скаттерограммы для выбранного при помощи меню записи сигнала. При выборе строки меню, соответствующей отсутствию сигнала, поля графиков очищаются. Для получения более наглядного графиков ритмограммы, гистограммы и скаттерограммы рассчитали верхний предел Ylim как 1,2 от максимального RR-интервала. В окне hAxes3 отображена нормированная гистограмма RR-интервалов, функция PH=H/SH*100; позволяет получить гистограмму в процентах. Рассчитали статистические параметры вариабельности сердечного ритма, описанные в разделе «Основные положения»: средний RR-интервал, минимальный RR-интервал, максимальный RR-интервал, вариационный размах, амплитуду моды (aMo) и индекс моды (iMo), моду (Mo).