利用ECG信号进行HRV分析

HRV（HeartRateVariability,中文为心率变异度），在医学上有重要的意义。心脏除了本身的节律性放电引发的跳动之外，也受到自律神经系统所调控。过去二十年已有不少文献显示自律神经系统的调控与心血管疾病相关的死亡率有显著的关系，例如心因性垂死、高血压、出血性休克、败血性休克等。心率变异分析亦被发现可作为预测发生心肌梗塞后的死亡率的指标及预测末期肝癌病患的预后。

HRV的周期频率一般被分为三部分，通过时频域变换可以得到VLF、LF、HF分别对应超低频（<0.04Hz)、低频段（0.04-0.15Hz）、高频段（0.15-0.4Hz），每个频段的功率（或者说幅度）都代表了不同的生理信息，利用离散傅立叶变换将心跳间隔的时间序列转换为频域，以功率频谱密度（Power

spectral density）或是频谱分布(Spectral

distribution)的方式表现。一般心率变异讯号的频谱分析使用200至500连续心跳间期稳定记录表现，因此记录需要数分钟的时间。

一般的心跳间期频谱频率出现在1赫兹以下，在0到0.4赫兹的范围内可找到数个波峰。主要为高频区(0.15-0.40赫兹)及低频区(0.04-0.15赫兹)。高频区通常反映副交感神经的活性，低频区同时受到交感与副交感神经系统的调控。在这里摘录维基百科的一个表格，请看附录！

通常利用连续量测到的心电图波形，直接计算与分析其相连心跳间时间序列的关系，例如：

SDNN (Standard Deviation of Normal to

Normal)，全部正常心跳间距之标准差，单位为毫秒。

SDANN (Standard deviation of the averages of NN intervals in

all 5-minute segments of the entire

recording)，全程依五分钟分成连续的时段，先计算每五分钟心跳间期的标准差，再计算标准差的平均值，单位为毫秒。

NN50 count (Number of pairs of adjacent NN intervals differing

by more than 50 ms in the entire

recording)，心电图中所有每对相邻正常心跳时间间隔，差距超过50毫秒的数目。

pNN50 (NN50 count divided by the total number of all NN

intervals)，NN50数目除以量测之心电图中所有的正常心跳间隔总数。

HRV代表了心动周期的变化，所以我们首先要求出每个变化着的周期的具体值，一般方法是找到每个周期识别点，比如过零点、最大最小极值点以及一些能够容易检测到的点，在ecg信号中我们一般采用RR点之间的时间作为对应的周期，所以第一步是检测ecg信号的R点；第二步就是要计算出每个周期值，也就是RR点之间的时间值，这个在matlab中很好实现；因为HRV周期是随时间变化的函数，而且一般对应的周期点数有限，说以在fft之前进行插值计算是必要的，第三步，插值法；最后一步进行fft变换，得到HRV的频谱图，剩下的分析就交给医生了。

按照上面所说的四个步骤，下面的程序分为四段，有详细的说明，请读者自行阅读。

close all

clear

load ekg.mat%读入ecg信号

[map,r,delay]=pan_tompkin(ecg,fs,0)% 利用pan_tomkin算法找到R点

[a,l]=size(r)

for i=2:l

t(i-1)=r(i)-r(i-1)%求出R-R间的时间值，即使HRV

end

x=r(2:19)

y=interp1(x,t,r(2):1:r(19),'spline')%利用插值法求出以原ecg信号的采样率fs的拟合函数

plot(y)hold on,

scatter(r(2:19)-r(2),t(1:18))

N=length(y)

N1=20%确定频率轴的范围 每一单元为fs/N=0.06Hz

AF=fft(y)

AF=abs(AF)%求出傅里叶变换后的幅频特性

f=(0:N1-1)*fs/N

figure,plot(f,AF(1:N1))%这一步是对fft的频率轴进行一定的变化，但是不会影响到fft的具体信息，只是为了方

%便显示而已，具体的fft使用请参看上篇关于fft的介绍。

四、Matlab结果截图

通过pin_tomkin算法过程中产生的相应图片

通过'spline'逐段3次样条插值得到的函数曲线，圆心点为实际的周期点

最终的HRV的幅频图

本例的ecg信号来源于MIT的生理信号库，但是忘记了该信号属于哪一类，有什么对应生理特征，在这里作者知识想表达处理的过程。由fft图知：该ecg的HRV主要处于<0.1Hz的范围类，代表交感与副交感神经活性，应该是一个正常人的eck信号。

心律变异度频域分析测量指标、定义及临床意义

指标单位定义频谱范围临床意义

总功率

total power, TPms2全部正常心跳间期之变异数高频、低频、极低频的总和≤0.4Hz整体心律变异度评估

极低频范围功率

very low frequency power, VLFPms2极低频范围正常心跳间期之变异≤0.04Hz生理意义不明

低频范围功率

low frequency power, LFPms2低频范围正常心跳间期之变异数0.04-0.15Hz代表交感与副交感神经活性

高频范围功率

high frequency power, HFPms2高频范围正常心跳间期之变异数0.15-0.4Hz代表副交感神经活性

标准化低频功率

normalized LFP, nLFP标准化单位,n.u.LF/(TP-VLF)交感神经活性

定量指标

标准化高频功率

normalized HFP,nHFP标准化单位,n.u.HF/(TP－VLF)副交感神经活

性定量指标

低、高频功率的比值

LF/HF无单位低、高频功率的比值代表自律神经

活性平衡

心率变异性（HRV）是指逐次心动周期的时间变化及其规律，反映了心率连续的瞬时波动，这种波动受体内神经、体液的调控，为适应不同的生理状况或某些病理生理状态而做出的反应。（心率随机体状况和昼夜时间变化，这种心率的规则性变化可通过正常QRS波群周期性变化的变异数来反映）。HRV提供了一种新的非创伤性方法对自主神经系统的活动情况进行评价。交感神经活动使心率缓慢变化，迷走神经活动使心率快速变化，即心率加快或减慢在1至数个心搏便可发生。呼吸性窦性心律不齐受迷走神经活动调节，主要反应在高频带（HF）交感神经和迷走神经两者共同调节反应在低频带（LF），但以交感神经调节为主；超低频（VLF）也反映交感神经活动的调节。LF/HF的比值反映心脏交感神经和迷走神经活动的平衡向高交感神经活动的倾斜。简言之，HRV是指窦性心律不齐的程度，HRV大表示心律不齐的程度重、心率变化大或变化快；HRV小表示心律不齐轻、心率变化小或是变化慢。目前研究表明，某些心血管疾病，如心肌梗死、扩张型心肌病、糖尿病性心肌病、充血性心力衰竭等，其心率变异性降低。临床实践证明，HRV可作为心源性猝死危险性的一个独立预测指标，在评价心血管疾病的预后和预测急性心血管事件方面具有重要意义。

目录1 拼音2 英文参考3 概述4 心率变异性的别名5 心率变异性的适应证6 心率变异性的禁忌证7 准备8 *** 作方法和有关指标 8.1 1.时域分析指标8.2 2.频域分析指标8.3 3.图解法指标 9 注意事项 1 拼音

xīn lǜ biàn yì xìng

2 英文参考

heart rate variability

HRV

3 概述

心率变异性（HRV）即窦性心律不齐的程度，是判断自主神经功能的一组常用定量指标。

4 心率变异性的别名

HRV

5 适应证

心率变异性适用于：

1.急性心肌梗死  HRV在心肌梗死后呈动态演变，HRV恢复的快慢及心肌梗死后3周左右的HRV的大小有预后意义。

2.糖尿病  HRV是判断糖尿病患者是否伴有自主神经系统损害较为准确、敏感的指标，优于一些传统方法。

3.其他心血管疾病  对二尖瓣脱垂综合征、肥厚性心肌病、长QT间期综合征、扩张性心肌病、高血压病、心脏移植、心力衰竭等，HRV有一定的评估病情和预后意义。

4.非心脏疾病  胎儿宫内窒息、帕金森病、多发性硬化、吉兰巴雷综合征、血管迷走性晕厥及 *** 性低血压等疾病对HRV有影响。

6 禁忌证

除患者因精神病等原因不能配合、电极贴敷过敏、经费等原因外无禁忌证。

7 准备

向患者介绍本检查的目的、方法和注意事项，取得患者的合作。

8 *** 作方法和有关指标

多数的动态心电图记录仪中带有心率变异性的分析软件，心率变异性的 *** 作方法及程序详见动态心电图的 *** 作。心率变异性主要包括以下几项指标：

8.1 1.时域分析指标

（1）SD：正常RR间期标准差（ms）。

（2）ASD：连续5min正常RR间期标准差均值（ms）。

（3）SDA：连续5min正常RR间期均值的标准差（ms）。

（4）RMSSD：相邻正常RR间期差值均方根（ms）。

（5）PNN50：相邻正常RR间期超过50ms的百分比。

前4个指标分别被正常RR间期均数（秒）除，得如下各心率校正值：

（6）SDC：SD的心率校正值。

（7）ASDC：ASD的心率校正值。

（8）SDAC：SDA的心率校正值。

（9）RMSDC：RMSSD的心率校正值。

8.2 2.频域分析指标

（1）ULF：超低频成分，频谱范围0～0.0033Hz。

（2）VLF：极低频成分，频谱范围0.003 3～0.04Hz。

（3）LF：低频成分，频谱范围0.04～0.15Hz。

（4）HF：高频成分，频谱范围0.15～0.40Hz。

（5）LH：LF/HF低频与高频成分之比。

（6）TOT：总成分，频谱范围0～0.5 Hz。

8.3 3.图解法指标

（1）三角指数：正常窦性心搏（NN）间期的总个数除以NN间期直方图的高度。

（2）TINN：用最小方差法，求出全部NN间期的直方图近似三角形底边宽度（ms）。

9 注意事项

1.HRV时域分析以24h为宜，不宜取任意时间段分析。

2.各项时域指标不能相互取代，不能交叉比较。

3.HRV三角指数的计算结果与时间单位直接相关。

4.频域分析对短时程和长时程分析结果的意义有很大差别。

5.短时程分析采样过程中最好避免有早搏、漏搏等情况。

---