RU2326587C1 - Способ исследования вариабельности сердечного ритма человека - Google Patents
Способ исследования вариабельности сердечного ритма человека Download PDFInfo
- Publication number
- RU2326587C1 RU2326587C1 RU2007102732/14A RU2007102732A RU2326587C1 RU 2326587 C1 RU2326587 C1 RU 2326587C1 RU 2007102732/14 A RU2007102732/14 A RU 2007102732/14A RU 2007102732 A RU2007102732 A RU 2007102732A RU 2326587 C1 RU2326587 C1 RU 2326587C1
- Authority
- RU
- Russia
- Prior art keywords
- wavelet
- scale
- formula
- frequency
- power density
- Prior art date
Links
Images
Landscapes
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
Изобретение относится к медицине и предназначено для исследования вариабельности сердечного ритма человека. Регистрируют R-R интервалограмму. Проводят ее спектральный анализ методом непрерывного вейвлет-преобразования. Определяют мощность частоты кардиоинтервалограммы а в момент времени b по формуле
где W(a,b) - коэффициент вейвлетного преобразования; f(t) - анализируемая функция; ψ((t-b)/a) - анализирующий вейвлет. Осуществляют построение на основе вейвлетных коэффициентов скейлограмм на отрезке [bi,bj] по формуле
где V(al) - скейлограмма сигнала; N - количество коэффициентов; al - масштаб вейвлетного преобразования. Выделяют на скейлограммах физиологически значимые частотные диапазоны, исходя из расстояний между соседними локальными минимумами на кривой скейлограммы. Определяют значения вейвлетной плотности мощности U в каждом из частотных диапазонов и значения удельной вейвлетной плотности мощности. Предлагаемый способ позволяет повысить точность определения параметров вариабельности за счет учета индивидуальных реакций исследуемого пациента. 4 ил.
Description
Изобретение относится к медицине, а именно к кардиологии, нормальной физиологии, патологической физиологии.
Известен способ изучения вариабельности сердечного ритма (ВСР), который осуществляется с помощью методик, использующих режимы временного и частотного анализа (Л.И.Макаров. Холтеровское мониторирование. - 2000. - М: Медпрактика, - с.51-62). Временной анализ основан на вычислении ряда статистических параметров серии R-R-интервалов различной продолжительности (mean, SDNN, SDNN-i, SDANN-i, rMSSD, pNN50, SDSD, Counts) (Crowford MH, Bernstein SJ, Deedwania PC et al. AHA guidelines for ambulatory electrocardiography: a report of the American College of Cardiology / American Heart Association Task Force of Practice Guidelines (Commitee to Revise the Guidelines for Ambulatory Electrocardiography). J. Am. Coil Cardiol. 1999; 34:912-48).
Также известны так называемые геометрические методы анализа ВСР - построение интервальной гистограммы, дифференциальной гистограммы различий, корреляционной ритмограммы. Оценка результатов геометрических методов осуществляется с помощью измерения параметров построенных геометрических фигур, аппроксимации паттерна сердечного ритма через построение фигур и математическое преобразование, с последующей интерпретацией и непосредственное описание и интерпретация формы геометрических фигур сердечного ритма.
Общими недостатками вышеперечисленных методик является необходимость длительной записи сигнала, усреднение результатов анализа, что не позволяет использовать их для определения локализации и характеристик локальных всплесков и отсутствие определения частотных характеристик ВСР.
При частотном (спектральном) анализе с помощью авторегрессионного анализа или ряда модификаций быстрого преобразования Фурье производится разделение серии R-R-интервалов на частотные спектры различной плотности. Определяются показатели спектральной мощности низкочастотного (Lf) и сверхнизкочастотного (VLF) диапазонов, характеризующие влияния симпатического и высокочастотного (Hf) диапазона, отражающего парасимпатические влияния на регуляцию ВСР, соотношение Lf/Hf, а также общая мощность спектра колебаний R-R-интервалов (TF). Главным недостатком этого метода является невозможность его использования при анализе нестационарных сигналов, т.е. большинства паттернов ВСР. Метод Фурье в модификации Уэлча обеспечивает усреднение результатов анализа в пределах ширины окна, т.е. он не позволяет установить частотно-временную локализацию быстро затухающего процесса. Применяемое для оценки временной динамики оконное преобразование Фурье имеет одно и тоже разрешение по времени и частоте для всех точек плоскости преобразования (Н.М.Астафьева. Успехи физических наук, т.166, №11, 1996, с.1145-1170, с.1150), что делает этот метод математического анализа малоинформативным для изучения временной динамики ВСР при переходных процессах и, следовательно, и при тонкой диагностике скрытых изменений вегетативного гомеостаза.
Наиболее близким к предлагаемому способу является способ применения вейвлет-преобразования для анализа вариабельности сердечного ритма детей (RU №2241374, МПК А61В 5/0452, 5/02, 10.12.2004). Недостатком данного способа является его малая информативность и точность, невозможность выявить слабые по силе воздействия на сердце, так как выделение частот, в которых производится анализ, осуществляется по стандартным частотным диапазонам (Рекомендации «Рабочей группы Европейского Кардиологического Общества и Северо-Американского Общества Стимуляции и Электрофизиологии») без учета индивидуальных особенностей распределения частотных спектров.
Задачей предлагаемого изобретения является улучшение оценки состояния вегетативной регуляции вариабельности сердечного ритма человека.
Технический результат заключается в увеличении точности и информативности способа исследования вариабельности сердечного ритма человека.
Технический результат достигается тем, что способ исследования вариабельности сердечного ритма человека включает регистрацию R-R интервалограммы и дальнейший ее спектральный анализ методом непрерывного вейвлет-преобразования, в котором определяют мощность частоты кардиоинтервалограммы а в момент времени b по формуле
где W(a,b) - коэффициент вейвлетного преобразования;
f (t) - анализируемая функция;
ψ((t-b)/a) - анализирующий вейвлет;
построение на основе матрицы вейвлетных коэффициентов скейлограмм на отрезке [bi,bj] по формулам
где V(al) - скейлограмма сигнала;
N - количество коэффициентов;
al - масштаб вейвлетного преобразования;
выделение на скейлограммах физиологически значимых частотных диапазонов, согласно изобретению выделение на скейлограммах физиологических значимых частотных диапазонов осуществляется, исходя из расстояний между соседними локальными минимумами на кривой скейлограммы по формуле
Δа=аm-an,
где Δа - физиологически значимый диапазон,
am, an - соседние локальные минимумы на кривой скейлограммы;
определение значения вейвлетной плотности мощности (ВПМ) U в каждом из частотных диапазонов Δа=[аm, an] осуществляется по формуле
определение изменения вейвлетной плотности мощности во времени как U(t);
определение изменения частотных диапазонов во времени как Δа(t);
определение значения удельной вейвлетной плотности мощности U' во времени по формуле
U'=U(t)/Δa(t),
которая отражает динамику изменения тонуса симпатического и парасимпатического отделов ВНС на коротких промежутках времени.
Скейлограммы («энергетические» диаграммы) строятся на основе матрицы вейвлет-коэффициентов, заданные как среднее квадратов коэффициентов W(a,b) при фиксированном параметре а на отрезке [bi,bj]. Являясь функцией масштаба, скейлограмма отражает ту же информацию, что и спектральная плотность мощности Фурье, являющаяся функцией от частоты. Как известно, вейвлет-преобразование имеет преимущество прежде всего за счет свойства частотно-временной локализации вейвлетов. Вейвлет-преобразование, представляющее собой временную развертку спектра, позволяет получить и более локализованную во времени энергетическую информацию. Энергетические диаграммы (скейлограммы) строятся на кратковременных (порядка 2-3 секунд) отрезках, что позволяет отслеживать временную динамику процесса.
На скейлограммах выделяют локальные спектры и физиологически значимые частотные диапазоны Δа, которые рассчитывают, исходя из расстояний между локальными минимумами am, an, связанными с различными типами механизмов регуляции ВСР человека. При этом при выявлении трех наиболее значимых диапазонов определяются два наиболее выраженных минимума, при четырех - три и т.д.
Суммарное значение вейвлетной плотности мощности U отражает суммарную активность нервного центра и определяется в каждом из частотных диапазонов Δа=[аm,an].
Удельная вейвлетная плотность мощности U' характеризует удельную выраженность активного нервного центра и отражает процессы оптимизации частоты сердечных сокращений. Выделение физиологически значимых диапазонов между локальными минимумами на кривой скейлограммы, связанных с различными типами механизмов регуляции ВСР, и оценка данного параметра позволяют выявить даже слабые по силе воздействия вегетативной нервной системы на ВСР на различных этапах онтогенеза, в норме и патологии, как в покое, так и при переходных процессах, что качественным образом повышает информативность и точность способа оценки вариабельности сердечного ритма человека. Оценка динамики данного параметра во времени позволяет описать динамику изменения тонуса симпатического и парасимпатического отделов ВНС в покое на коротких промежутках времени.
На фиг.1 показаны скейлограммы (локальные вейвлетные спектры) двух разных участков кривой ВСР пациента К.
На фиг.2 показаны 2 различных увеличенных участка разделения скейлограмм по частотам в случае стандартного выделения физиологически значимых диапазонов (соответственно а и б).
На фиг.3 показаны соседние скейлограммы, построенные на основе кривой ВСР. Скейлограммы строились на трех последовательных пятисекундных отрезках (соответственно кривая 1, кривая 2, кривая 3).
На фиг.4 показана вейвлет-диаграмма кривой вариабельности сердечного ритма. По оси абсцисс - время, по оси ординат - масштаб (величина, обратная частоте). Светло-серые участки соответствуют максимумам, темно-серые - минимумам. Нижняя серая линия соответствует частотной границе между HF и LF диапазонами, средняя серая линия - границе между LF и VLF диапазонами.
Клинический пример.
Испытуемый К. 22 года. Рост 187 см, вес 78 кг. Регистрация вариабельности сердечного ритма осуществлялась в положении сидя, в покое, показатели спектрального анализа по Фурье ТР, мс^2=10377, VLF, мс^2=7695, LF, мс^2=1624, HF, мс^2=1057, LF norm, n.u.=60.6, HF norm, n.u.=39.4, LF/HF=1.54, %VLF=74.2, %LF=15.7, %HF=10.2.
Вейвлет-преобразование: непрерывное вейвлет-преобразование, вейвлет morlet, максимальный масштаб 100, временное усреднение при построении скейлограмм 10 кардиоциклов. Результаты эксперимента приведены на фиг.1-4. Значения ВПМ рассчитываются в частотных диапазонах между минимумами на каждой скейлограмме.
Данные, приведенные на фиг.1, получены за счет локализованного спектрального анализа, свидетельствуют об изменениях частотных составляющих ВСР и показывают вклад различных частот в общую картину ВСР.
На фиг.2а видно, что разделение частотных диапазонов симпатического (без выделения) и парасимпатического (серое выделение) отделов вегетативной нервной системы проведено успешно, т.к. граница проходит в области локального минимума. На фиг.2б разделение частотных диапазонов симпатического (без выделения) и парасимпатического (серое выделение) отделов вегетативной нервной системы проведено неудачно: граница прошла не по локальному минимуму, а вблизи локального максимума, тем самым искусственно разделив пик скейлограммы. Таким образом, из приведенных данных на фиг.2 видно, что разделение частотных диапазонов согласно стандартному выделению физиологически значимых диапазонов часто оказывается искусственным и не всегда совпадает с реальными границами генерации модулирующих влияний на ВСР симпатического и парасимпатического отделов вегетативной нервной системы.
Из данных, приведенных на фиг.3, видно, что частоты локальных максимумов и локальных минимумов постоянно изменяются и, следовательно, константная граница частотных диапазонов не отражает истинные характеристики процесса.
На фиг.4 видно, что при проведении разделения частотных диапазонов согласно предложенному изобретению диапазоны не постоянны, а изменяются в достаточно широких пределах.
Claims (1)
- Способ исследования вариабельности сердечного ритма человека, заключающийся в регистрации R-R интервалограммы и ее спектральном анализе методом непрерывного вейвлет-преобразования, включающим определение мощности частоты кардиоинтервалограммы а в момент времени b по формулегде W(a,b) - коэффициент вейвлетного преобразования; f(t) - анализируемая функция; ψ((t-b)/a) - анализирующий вейвлет;построение на основе вейвлетных коэффициентов скейлограмм на отрезке [bi,bj] по формулегде V(a1) - скейлограмма сигнала; N - количество коэффициентов; a1 - масштаб вейвлетного преобразования;выделение на скейлограммах физиологически значимых частотных диапазонов, отличающийся тем, что выделение на скейлограммах физиологических значимых частотных диапазонов осуществляется исходя из расстояний между соседними локальными минимумами на кривой скейлограммы по формулеΔа=аm-аn,где Δа - физиологически значимый диапазон, аm,an - соседние локальные минимумы на кривой скейлограммы,определение значения вейвлетной плотности мощности U в каждом из частотных диапазонов Δa=[am,an] осуществляется по формулеопределение изменения вейвлетной плотности мощности во времени как U(t);определение изменения частотных диапазонов во времени как Δa(t);определение значения удельной вейвлетной плотности мощности U' во времени по формулеU'=U(t)/Δa(t).
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU2007102732/14A RU2326587C1 (ru) | 2007-01-24 | 2007-01-24 | Способ исследования вариабельности сердечного ритма человека |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU2007102732/14A RU2326587C1 (ru) | 2007-01-24 | 2007-01-24 | Способ исследования вариабельности сердечного ритма человека |
Publications (1)
Publication Number | Publication Date |
---|---|
RU2326587C1 true RU2326587C1 (ru) | 2008-06-20 |
Family
ID=39637217
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
RU2007102732/14A RU2326587C1 (ru) | 2007-01-24 | 2007-01-24 | Способ исследования вариабельности сердечного ритма человека |
Country Status (1)
Country | Link |
---|---|
RU (1) | RU2326587C1 (ru) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2521729C1 (ru) * | 2012-12-07 | 2014-07-10 | Федеральное государственное бюджетное учреждение науки Институт проблем управления им. В.А. Трапезникова Российской академии наук | Бесконтактный радиоволновой способ измерения уровня жидкости в емкости |
RU2546103C2 (ru) * | 2013-03-06 | 2015-04-10 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Калининградский государственный технический университет" | Способ определения параметров вариабельности сердечного ритма |
IL289321B1 (en) * | 2021-12-23 | 2023-06-01 | Neteera Tech Ltd | Derivation of heartbeat interval hundreds of returns |
-
2007
- 2007-01-24 RU RU2007102732/14A patent/RU2326587C1/ru not_active IP Right Cessation
Non-Patent Citations (1)
Title |
---|
БОЙЦОВ С.А. и др., Анализ сигнал-усредненной ЭКГ (по данным вейвлет-преобразования) у здоровых и больных ИБС, Вестник Аритмологии. 2001, №23, с.32-35. АСТАФЬЕВА Н.М. Вейвлетный анализ: основы теории и примеры применения, Успехи физиол. наук, Т.166, 1996, №11, с.1145-1170. Heart Rate Variability Standarts of measurement, physiological interpretation and clinical use. Europ. Heart Journal. 1996, v.l7, p.354-381. * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2521729C1 (ru) * | 2012-12-07 | 2014-07-10 | Федеральное государственное бюджетное учреждение науки Институт проблем управления им. В.А. Трапезникова Российской академии наук | Бесконтактный радиоволновой способ измерения уровня жидкости в емкости |
RU2546103C2 (ru) * | 2013-03-06 | 2015-04-10 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Калининградский государственный технический университет" | Способ определения параметров вариабельности сердечного ритма |
RU2546103C9 (ru) * | 2013-03-06 | 2016-01-27 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Калининградский государственный технический университет" | Способ определения параметров вариабельности сердечного ритма |
IL289321B1 (en) * | 2021-12-23 | 2023-06-01 | Neteera Tech Ltd | Derivation of heartbeat interval hundreds of returns |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8233972B2 (en) | System for cardiac arrhythmia detection and characterization | |
EP3358485A1 (en) | General noninvasive blood glucose prediction method based on timing analysis | |
Zheng et al. | Computer-assisted diagnosis for chronic heart failure by the analysis of their cardiac reserve and heart sound characteristics | |
KR100493714B1 (ko) | 자율신경 검사장치 | |
US20150105666A1 (en) | Narrow band feature extraction from cardiac signals | |
US8280501B2 (en) | Systems and methods for quantitatively characterizing slow wave activities and states in sleep | |
CN110037668B (zh) | 脉搏信号时空域结合模型判断年龄、健康状态及恶性心律失常识别的系统 | |
US20220280096A1 (en) | Ssvep-based attention evaluation method, training method, and brain-computer interface | |
US8768442B2 (en) | Wakeful-state data generating apparatus and wakefulness degree determining apparatus | |
CN108992053B (zh) | 一种实时的无束缚检测心率和心跳间隔的方法 | |
EP1118945A1 (en) | Method and system for measuring heart rate variability | |
RU2326587C1 (ru) | Способ исследования вариабельности сердечного ритма человека | |
WO2018130897A1 (en) | Method and system for determining heart rate variability | |
US20220183631A1 (en) | Automatic test device and method for auditory brainstem response | |
Rohlén et al. | Optimization and comparison of two methods for spike train estimation in an unfused tetanic contraction of low threshold motor units | |
Hansson-Sandsten et al. | Multiple window correlation analysis of HRV power and respiratory frequency | |
CN108968946A (zh) | 一种基于hrv分析的女性内分泌管理系统 | |
Madhav et al. | Monitoring respiratory activity using PPG signals by order reduced-modified covariance AR technique | |
CN116616709A (zh) | 基于非静止状态下多通道视频的工作记忆负荷评估方法 | |
CN115969395A (zh) | 一种无创评估大脑皮层兴奋性的系统及方法 | |
CN114098645A (zh) | 一种睡眠分期方法及装置 | |
Chen et al. | A fast ECG diagnosis using frequency-based compressive neural network | |
RU2241374C2 (ru) | Способ исследования вариабельности сердечного ритма у детей | |
Fedotov et al. | Methods of mathematical analysis of heart rate variability | |
RU2328973C1 (ru) | Способ исследования регуляции микроциркуляторного русла человека |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
MM4A | The patent is invalid due to non-payment of fees |
Effective date: 20090125 |