「心率变异性分析——ECG信号处理-第十二课」2025年4月1日

刘小小_算法工程师

已于 2025-04-01 17:17:14 修改

阅读量795

点赞数 14

文章标签：信号处理

于 2025-04-01 14:16:00 首次发布

本文链接：https://blog.csdn.net/qq_42936724/article/details/146904492

版权

心率变异性（Heart Rate Variability，HRV）是指逐次心跳周期差异的变化情况，它反映了心脏自主神经调节功能的动态变化。HRV分析在临床诊断、运动医学、心理生理学等领域具有重要的应用价值。本文将详细介绍心率变异性分析的流程，包括时间域分析、非线性分析、RR间期直方图分析以及频域分析，并结合正常参数范围和数据分析健康情况。

图1：时域与频域分析示意

一、HRV时域分析

时域分析是最直观的心率变异性分析方法之一，它直接对RR间期序列进行统计分析。RR间期是指连续两个R波之间的时间间隔，通常以毫秒（ms）为单位。时域分析的关键指标包括以下几个方面：

1.1 SDNN（标准差）

定义：SDNN表示RR间期的标准差，反映了心率变异性的总体水平。它衡量了RR间期的离散程度，值越大，说明心率变化越不稳定，自主神经调节功能可能越活跃。

计算公式：

其中，RRi表示第i个RR间期，RRˉ表示所有RR间期的平均值，N为RR间期的总数。

正常范围：50–150 ms（长期记录），20–50 ms（短期记录）。

健康评估：较高的SDNN值通常表明自主神经调节功能较好，个体处于较为放松的状态。较低的SDNN值可能表明自主神经调节功能受损，可能与心血管疾病、糖尿病等慢性疾病有关。

1.2 NN50（超过50ms的RR间期差值数量）

定义：NN50统计的是相邻RR间期差值大于50ms的次数。它反映了心率的快速变化情况，较高的NN50值表明心率在短时间内有较大的波动。

计算公式：

其中，RRi+1和RRi分别表示相邻的两个RR间期。

正常范围：短期记录约20–50次，长期记录可能达到数百次。

健康评估：较高的NN50值表明心率的快速变化能力较强，通常与较好的副交感神经功能相关。较低的NN50值可能表明副交感神经功能减弱，可能与应激状态、焦虑等有关。

1.3 pNN50（NN50占总RR间期数的百分比）

定义：pNN50是NN50占总RR间期数的百分比。它进一步量化了NN50的意义，能够更准确地反映心率的快速变化程度。

计算公式：

其中，N为RR间期的总数。

正常范围：短期记录约5%-20%。

健康评估：较高的pNN50值表明心率的快速变化能力较强，通常与较好的副交感神经功能相关。较低的pNN50值可能表明副交感神经功能减弱，可能与应激状态、焦虑等有关。

二、HRV频域分析

频域分析是通过将RR间期序列转换到频域，分析其功率谱密度（PSD），从而评估心率变异性的频率成分。频域分析的关键指标包括以下几个方面：

2.1 LF（低频功率）

定义：LF表示低频范围内的功率，通常频率范围为0.04-0.15 Hz。它主要反映了交感神经和副交感神经的共同作用，但以前者为主。

计算公式：

其中，P(f)表示功率谱密度。

正常范围：短期记录 200–800 ms²，长期记录 800–2500 ms²。

健康评估：较高的LF值可能表明交感神经活动较强，个体可能处于应激状态；较低的LF值可能表明交感神经活动较弱，通常与健康状态相关。

2.2 HF（高频功率）

定义：HF表示高频范围内的功率，通常频率范围为0.15-0.4 Hz。它主要反映了副交感神经的活动。

计算公式：

其中，P(f)表示功率谱密度。

正常范围：短期记录 150–1000 ms²，长期记录 500–2000 ms²。

健康评估：较高的HF值可能表明副交感神经活动较强，个体可能处于放松状态；较低的HF值可能表明副交感神经活动较弱，可能与应激状态、焦虑等有关。

2.3 TP（总功率）

定义：TP表示总功率，是LF和HF功率之和。它反映了心率变异性的总体水平。

计算公式：

正常范围：短期记录 1000–3000 ms²，长期记录 3000–10000 ms²。

健康评估：较高的TP值表明心率变异性较好，自主神经调节功能较强；较低的TP值可能表明自主神经调节功能受损，可能与心血管疾病、糖尿病等慢性疾病有关。

2.4 LF/HF（低频与高频的比值）

定义：LF/HF比值用于评估交感神经和副交感神经的平衡状态。较高的LF/HF比值可能表明交感神经占优势，而较低的比值则可能表明副交感神经占优势。

计算公式：

正常范围：健康成年人的LF/HF通常在0.5–2.0之间。

健康评估：较高的LF/HF比值可能表明交感神经占优势，可能与应激状态、焦虑等有关；较低的LF/HF比值可能表明副交感神经占优势，通常与健康状态相关。

三、非线性分析 Poincare图

Poincare图是一种直观展示RR间期动态变化的非线性分析方法。它通过绘制相邻RR间期的散点图，能够揭示心率变异性的复杂性。

图2：Poincare图示意

Poincare图的关键指标包括以下几个方面：

3.1 SD1（短轴标准差）

定义：SD1表示RR间期变化的短期波动性，反映了心率的快速变化情况。

计算公式：

其中，RRi+1和RRi分别表示相邻的两个RR间期。

正常范围：健康成年人的SD1通常在 20–70 ms之间。

健康评估：较高的SD1值表明心率的快速变化能力较强，通常与较好的副交感神经功能相关。较低的SD1值可能表明副交感神经功能减弱，可能与应激状态、焦虑等有关。

3.2 SD2（长轴标准差）

定义：SD2表示RR间期变化的长期波动性，反映了心率的长期变化趋势。

计算公式：

其中，std(RRi)表示所有RR间期的标准差。

正常范围：健康成年人的SD2通常在100-300 ms之间。

健康评估：较高的SD2值表明心率的长期变化能力较强，通常与较好的自主神经调节功能相关。较低的SD2值可能表明自主神经调节功能受损，可能与心血管疾病、糖尿病等慢性疾病有关。

3.3 SD1/SD2（SD1与SD2的比值）

定义：SD1/SD2用于评估心率变异性的复杂性。该比值越小，说明心率变异性越复杂，自主神经调节功能可能越活跃。

计算公式：

正常范围：健康成年人的SD1/SD2通常在 0.2-1.0 之间。

健康评估：较低的SD1/SD2值表明心率变异性较为复杂，自主神经调节功能较好。较高的SD1/SD2值可能表明心率变异性较为简单，可能与应激状态、焦虑等有关。

3.4 EA（椭圆面积）

定义：EA表示Poincare图中散点分布的椭圆面积，反映了心率变异性的总体水平。

计算公式：

正常范围：健康成年人的EA通常在3000–15000 ms²之间。

健康评估：较高的EA值表明心率变异性较好，自主神经调节功能较强。较低的EA值可能表明自主神经调节功能受损，可能与心血管疾病、糖尿病等慢性疾病有关。

3.5 VAI（角度变异指数）

定义：VAI表示相邻RR间期差值的角度变化情况，反映了心率变化的方向性。

计算公式：

正常范围：健康成年人的VAI通常在0.1-0.5 rad之间。

健康评估：较高的VAI值表明心率变化的方向性较强，通常与较好的自主神经调节功能相关。较低的VAI值可能表明自主神经调节功能受损，可能与心血管疾病、糖尿病等慢性疾病有关。

3.6 VLI（长度变异指数）

定义：VLI表示相邻RR间期差值的长度变化情况，反映了心率变化的幅度。

计算公式：

正常范围：健康成年人的VLI通常在 50–200 ms之间。

健康评估：较高的VLI值表明心率变化的幅度较大，通常与较好的自主神经调节功能相关。较低的VLI值可能表明自主神经调节功能受损，可能与心血管疾病、糖尿病等慢性疾病有关。

四、RR间期直方图

RR间期直方图用于展示RR间期的分布情况，它能够直观地反映心率的变化范围和频率。在直方图中，横轴表示RR间期的范围，纵轴表示每个RR间期范围内的概率。通过观察直方图，可以了解RR间期的集中趋势和离散程度。

图3：直方图示意

例如，如果直方图呈现出明显的单峰分布，说明心率变化相对稳定；如果直方图呈现出多峰分布，说明心率变化较为复杂。

单峰分布：表明心率变化较为稳定，通常与健康状态相关。

多峰分布：表明心率变化较为复杂，可能与应激状态、焦虑等有关。

五、算法实现与结果分析

5.1 算法实现过程

以下是基于 Pan-Tompkins 算法的心率变异性（HRV）分析的 Matlab 代码实现结果。实验使用了 MIT-BIH 心电数据库中的 221 号记录信号，采样频率为 360Hz。通过对该信号进行一系列处理和分析，我们得到了以下结果，并进行了相应的健康评估。

（1）从指定路径读取 ECG 信号数据，并进行预处理，包括去除基线漂移、肌电干扰和工频干扰。相关内容见我主页的如下链接：

「ECG信号处理-第十一课——心电信号中基线漂移、肌电干扰与工频干扰联合滤除」2025年3月31日-CSDN博客

预处理后的信号如图 4 和图 5 所示，可以看出，经过逐步处理，ECG 信号的噪声得到了有效去除。

图4：原始ECG预处理去噪结果

图5：去除的噪声（基线漂移、肌电干扰和工频噪声）

（2）使用 Pan-Tompkins 算法对预处理后的 ECG 信号进行 R 波检测，提取RR间期序列。相关内容见我主页的如下链接：

「ECG信号处理-第十课——心电信号特征提取——Pan-Tomkins算法（R波定位与检测）」2025年3月17日_心电r峰检测pantomking算法详解-CSDN博客

图6：Pan-tompkins算法R波定位结果

图7：提取的RR间期序列

（3）对检测到的 RR 间期序列进行时间域分析、频域分析、非线性分析（Poincaré 图）并绘制RR间期直方图。

图8：非线性分析 Poincaré 图（庞加莱图）

图9：RR间期直方图

图10：心率变异性量化指标

5.2 实验结果分析

(1) 患者记录的头文件信息

我们首先从222号记录的头文件中提取的信息提供了病人的基本临床背景：

病人基本信息：该患者为一位84岁的女性，正在接受**地高辛（Digoxin）和奎尼丁（Quinidine）**治疗。这两种药物常用于控制心律失常，尤其是房颤和房扑的治疗。
临床背景：根据记录，该患者经历了阵发性房颤/房扑（paroxysmal atrial flutter/fibrillation），这种心律失常通常伴随着窦房结逃逸（nodal escape beats）。房颤和房扑是由心房异常放电引起的，这导致了心率不规则并对自主神经系统造成较大的影响，进而影响心率变异性（HRV）。