第２８卷第５期 ２０１３年９月 数 据 采 集 与 处 理 Ｖｏ１．２８ Ｎｏ．５ Ｓｅｐ．２０１３ Ｊｏｕｒｎａｌ ｏｆ Ｄａｔａ Ａｃｑｕｉｓｉｔｉｏｎ ａｎｄ Ｐｒｏｃｅｓｓｉｎｇ 文章编号：１００４—９０３７（２０１３）０５—０６３８－０５ 运动前后心率变异性的高维时间不可逆性分析 侯凤贞 宁新宝 黄晓林。 （１．中国药科大学理学院，南京，２１０００９；２．南京大学电子科学与工程学院生物医学电子研究所，南京，２１００９３） 摘要：对心率变异性的分析已经广泛应用于心血管疾病临床实践和基础研究中，并逐渐开始在运动领域中流行 起来。本文利用时间序列的高维时间不可逆性（Ｈｉｇｈ－ｄｉｍｅｎｓｉｏｎａ１ ｔｉｍｅ ｉｒｒｅｖｅｒｓｉｂｉｌｉｔｙ，ＨＤＴＩ）分析方法，研究受 试者在跳绳运动前、中、后不同状态下短时心率变异信号的时间不可逆特性。结果表明，高维的时间不可逆性分 析特征参量的变化正好反映了运动前、中、后交感一迷走神经之间从平衡态到非平衡态、再从非平衡态逐渐恢复 到平衡态的过程；在分析快速变化的心电数据时，该特征参量与心率变化的趋势密切相关。研究结果显示ＨＤ— Ｔ１分析方法有望在运动领域中得到进一步的应用。 关键词：心率变异性；时间不可逆性；运动生理学 中图分类号：Ｒ３１８．Ｏ４ 文献标志码：Ａ Ｈｉ ｇｈ—Ｄｉｍｅｎｓｉｏｎａｌ Ｔｉｍｅ ＩｒｒｅＶｅｒｓｉｂｉｌｉｔｙ Ａｎａｌｙｓｉｓ ｏｆ Ｈｕｍａｎ Ｈｅａｒｔ Ｒａｔｅ Ｖａｒｉａｂｉｌｉｔｙ Ｂｅｆｏｒｅ ａｎｄ Ａｆｔｅｒ Ｅｘｅｒｃｉｓｅ Ｈｏｕ Ｆｅｎｇｚｈｅｎ ，Ｎｉｎｇ Ｘｉｎｂａｏ。，Ｈｕａｎｇ Ｘｉａｏｌｉｎ （１．Ｄｅｐａｒｔｍｅｎｔ ｏｆ Ｓｃｉｅｎｃｅ，Ｃｈｉｎａ Ｐｈａｒｍａｃｅｕｔｉｃａｌ Ｕｎｉｖｅｒｓｉｔｙ，Ｎａｎｊｉｎｇ，２１０００９，Ｃｈｉｎａ； ２．Ｉｎｓｔｉｔｕｔｅ ｏｆ Ｂｉｏｍｅｄｉｃａｌ Ｅｌｅｃｔｒｏｎｉｃ Ｅｎｇｉｎｅｅｒｉｎｇ，Ｓｃｈｏｏｌ ｏｆ Ｅｌｅｃｔｒｏｎｉｃ Ｓｃｉｅｎｃｅ ａｎｄ Ｅｎｇｉｎｅｅｒｉｎｇ，Ｎａｎｉｉｎｇ Ｕｎｉｖｅｒｓｉｔｙ，Ｎａ ｉｎｇ，２１００９３，Ｃｈｉｎａ） Ａｂｓｔｒａｃｔ：Ｈｅａｒｔ ｒａｔｅ ｖａｒｉａｂｉｌｉｔｙ（ＨＲＶ）ａｎａｌｙｓｉｓ ｈａｓ ｂｅｅｎ ｗｉｄｅｌｙ ｕｓｅｄ ｉｎ ｃｌｉｎｉｃａｌ ｐｒａｃｔｉｃｅ ａｎｄ ｂａｓｉｃ ｒｅｓｅａｒｃｈｅｓ ｏｆ ｃａｒｄｉｏｖａｓｃｕｌａｒ ｄｉｓｅａｓｅｓ．Ｍｏｒｅｏｖｅｒ，ｉｔ ｇｒａｄｕａｌｌｙ ｂｅｃｏｍｅｓ ｐｏｐｕｌａｒ ｉｎ ｔｈｅ ｆｉｅｌｄ ｏｆ ｓｐｏｒｔｓ．Ｔｈｉｓ ｐａｐｅｒ ｕｔｉｌｉｚｅｓ ｔｈｅ ｍｅｔｈｏｄ ｆｏｒ ｈｉｇｈ—ｄｉｍｅｎｓｉｏｎａｌ ｔｉｍｅ ｉｒｒｅｖｅｒｓｉｂｉｌｉｔｙ（ＨＤＴＩ）ａ— ｎａｌｙｓｉｓ ｔｏ ｓｔｕｄｙ ｓｈｏｒｔ—ｔｅｒｍ ＨＲＶ ｓｉｇｎａｌｓ ｄｅｒｉｖｅｄ ｆｒｏｍ １３ ｓｕｂｊｅｃｔｓ ｗｈｅｎ ｔｈｅｙ ａｒｅ ｂｅｆｏｒｅ，ｄｕｒｉｎｇ ａｎｄ ａｆｔｅｒ ｓｋｉｐｐｉｎｇ ｅｘｅｒｃｉｓｅ．Ｔｈｅ ｒｅｓｕｌｔｓ ｓｈｏｗ ｔｈａｔ ｔｈｅ ＨＤＴＩ ｍｅａｓｕｒｅｍｅｎｔ ｊ ｕｓｔ ｒｅｆｌｅｃｔｓ ｔｈｅ ｃｈａｎｇｅｓ ｏｆ ｃａｒｄｉａｃ ｓｙｍｐａｔｈｏ——ｖａｇａｌ ｍｏｄｕｌａｔｉｏｎｓ ｆｒｏｍ ｅｑｕｉｌｉｂｒｉｕｍ ｂｅｆｏｒｅ ｓｋｉｐｐｉｎｇ ｔｏ ｎｏｎ—．ｅｑｕｉ—— ｌｉｂｒｉｕｍ ｄｕｒｉｎｇ ｓｋｉｐｐｉｎｇ，ａｎｄ ｔｈｅｎ ｒｅｃｏｖｅｒｉｎｇ ｆｒｏｍ ｎｏｎ—ｅｑｕｉｌｉｂｒｉｕｍ ｔｏ ｅｑｕｉｌｉｂｒｉｕｍ ｇｒａｄｕａｌｌｙ ａｆｔｅｒ ｓｋｉｐｐｉｎｇ．Ｆｕｒｔｈｅｒｍｏｒｅ，ｉｔ ｉｓ ｆｏｕｎｄ ｔｈａｔ ｔｈｅ ＨＤＴＩ ｍｅａｓｕｒｅｍｅｎｔ ｉｓ ｌｉｎｅａｒｌｙ ｃｏｒｒｅｌａｔｅｄ ｗｉｔｈ ｔｈｅ ｔｒｅｎｄ ｏｆ ｈｅａｒｔ ｒａｔｅ ｗｈｅｎ ａｎａｌｙｚｉｎｇ ｔｈｅ ｄａｔａ ｗｈｉｃｈ ｃｈａｎｇｅ ｑｕｉｃｋｌｙ．Ｔｈｅ ｒｅｓｅａｒｃｈ ｓｕｇｇｅｓｔｓ ｔｈａｔ ｔｈｅ ＨＤＴＩ ｍｅｔｈｏｄ ｃａｎ ｂｅ ｆｕｒｔｈｅｒ ａｐｐｌｉｅｄ ｉｎ ｔｈｅ ｆｉｅｌｄ ｏｆ ｓｐｏｒｔｓ． Ｋｅｙ ｗｏｒｄｓ：ｈｅａｒｔ ｒａｔｅ ｖａｒｉａｂｉｌｉｔｙ；ｔｉｍｅ ｉｒｒｅｖｅｒｓｉｂｉｌｉｔｙ；ｓｐｏｒｔｓ ｐｈｙｓｉｏｌｏｇｙ 的研究也逐渐开始在运动领域流行起来，研究主要 引 目 集中在如下几个方面Ｌ６ ］：分析不同类型的运动对 人体的影响并应用ＨＲＶ对健康个体锻炼的效果 进行评价；利用ＨＲＶ去评定运动后副交感神经张 力的恢复能力；研究自主神经与运动表现的关系； 尝试应用ＨＲＶ来评定运动能力，这方面的一个典 心率变异性（Ｈｅａｒｔ ｒａｔｅ ｖａｒｉａｂｉｌｉｔｙ，ＨＲＶ）是 反映心脏自主神经功能状态的一种简便、可靠的无 创性技术，它已经广泛应用于心血管疾病临床实践 和基础研究中［ 。。。上世纪９０年代以来，对ＨＲＶ 型例子便是研究利用ＨＲＶ来测定青少年运动员 基金项目 ：国家自然科学基金（６１２７１０８２）资助项目；江苏省自然科学基金（ＢＫ２０１１５６５）资助项目；江苏省“青蓝工程” 资助项目 ０ 收稿日期 ：２０１３－０７—０６；修订日期：２０１３－０８—０９ 第５期 侯凤贞，等：运动前后心率变异性的高维时间不可逆性分析 ６３９ 无氧阈。１９８２年ＣＯＮＣＯＮＩ通过大量的实验证明 在功率自行车进行等级递增强度运动时，心率随功 率输出递增，某一点上心率的增加率下降，行成拐 点，这一心率拐点对应的速度即为无氧阚速度或无 氧阈心率。但ＣＯＮＣＯＮＩ法的一个不足是，心率 拐点并不是在所有人身上都能出现，即便是训练有 素的运动员，也不一定会出现明显的心率拐点 Ｊ。 因此，一些研究者开始尝试提出一些与心率拐点相 关的ＨＲＶ指标，如“ＨＲＶ第二峰值”［９ 等，作为心 率拐点的替代来评估无氧阈并取得了一定的成果。 对ＨＲＶ的分析既可以基于心跳间期（ＲＲ间 期）序列，也可以基于瞬时心率（ＲＲ间期的倒数） 序列。无论采用哪种序列，归根结底都是对时间序 列的分析研究。本文尝试利用时间序列的高维时 间不可逆性分析（Ｈｉｇｈ—ｄｉｍｅｎｓｉｏｎａｌ ｔｉｍｅ ｉｒｒｅｖｅｒｓ— ｉｂｉｌｉｔｙ ａｎａｌｙｓｉｓ，ＨＤＴＩ）方法ｌ＿１ ，研究受试者在不 同状态下（安静休息状态、５ ｍｉｎ的跳绳状态、跳绳 后的休息状态）短时的ＲＲ间期信号（约５００点）的 时间不可逆性。迄今为止，尚未见将ＨＤＴＩ方法 应用于运动前后ＨＲＶ分析的报道。这样做的目 的有三个：一是探索高维的时间不可逆指标与心电 信号何种生理学特性有关；二是研究在短时（１～ ５ ｍｉｎ）的ＨＲＶ分析中，作为一种统计学方法， ＨＤＴ１分析方法受哪些因素的影响，有没有其适用 的地方；三是看能否将心率变异信号的ＨＤＴ１分 析方法应用到运动领域中去。 １ 高维的时间不可逆性分析方法 考虑一个Ｎ点的时间序列 ｛ｚ ） １≤ｉ≤Ｎ 把该时间序列嵌入 维相空间［ｉｏ３，即对于每 一个ｚ ，取ｍ个点组成一 维矢量 ｘ （ｉ）＝［ ，ｚ斗Ｌ，…，Ｘ ＋（聃－１）Ｌ］ （１） 式中：ｍ为嵌人维数，Ｌ为延迟时间。取Ｌ＝１，那 么ｍ维矢量的个数就是Ｎ— ＋１。 然后分析各ｍ维矢量对应的坐标在 维相空 间的各个互不平行的侧面（ｚ ，３Ｃ ）上的投影，其 中，ｚ一１，２，…，ｍ一１。按式（２，３）定义二维侧面 （．ｚ ，ｚ ）上投影图的不对称性指标Ｐ （”）， Ｇ （ ）［ Ｐ ）一 。（２） Ｎ（ａｎｎＴ ∑ａＲＲ．＋（ ）ｚ Ｇ （ ）＝—— ———上 —————＝————一×１００ Ｎ（ａＲＲ：） Ｎ（△ＲＲ ） ∑△ＲＲ （ ） ＋∑ＡＲＲＴ（ ）。 （３） 式中：△ＲＲ 一ｚ ＋ 一３２ ，ＡＲＲ２表示△ＲＲ 小于０ 的点，△ＲＲ 表示△ＲＲ 大于０的点。 再按式（４，５）重定义指标Ｐ ，Ｇ 来表示嵌入 空间中ｍ维矢量的不对称性ｎ ０＿ Ｐ 一 ∑Ｔｎ －：：ｌ ５０一Ｐ％（ｎ）ｌ （４） Ｇ 一 ∑ｍ一 －：１ ５０一Ｇ％（ ）Ｉ（５） 最后，定义指标Ｄ 来综合Ｐ 和Ｇ 的值 Ｄ 一￣／Ｐ ２＋Ｇ （６） 由上述过程可见，Ｄ 可以反映 维矢量在ｍ 个投影面上关于主对角线的不对称程度的平均状 况，也反映了序列在逆转前后统计上的差异，也就 是序列的时间不可逆性。Ｄ 越大，代表着序列的 时间不可逆性越明显口。。。 ２ 实验验证 ２．１ 实验数据的采集 本文的实验对象是１３个健康年轻人（平均年 龄２５岁，年龄范围２Ｏ～３５岁，仅两人在３０岁以 上），其中８位女性，５位男性，均无心血管疾病史。 受试者从１～１３连续编号。利用增强型三通道 Ｈｏｌｔｅｒ记录器来采集体表心电图（Ｅｌｅｃｔｒｏ ｃａｒｄｉｏ— ｇｒａｍ，ＥｃＧ）信号，采样率为１２８ Ｈｚ，Ａ／Ｄ转换精 度为１２位。设备配有动态心电图分析软件，能自 动定位Ｒ波顶点并识别其是正常的窦性起搏点还 是异位起搏点。对于异位起搏点，以其前后两点时 刻的中间点来估计被抑制的窦性心率，并以此推算 其对应的ＲＲ间期值。 每位受试者连续采集了至少约１ ｈ的心电信 号。开始采集后先安静０．５ ｈ，但可进行一些日常 活动（后文将该状态表示为“跳绳前”），紧接着跳绳 ５ ｒａｉｎ（后文将该状态表示为“跳绳中”），跳绳结束 后每位受试者依个人情况安静散步几分钟至半小 时（后文将该状态表示为“跳绳后”）。记录每位受 试者开始跳绳和结束跳绳的时间（精确到秒）。 ２．２不同状态下ＨＲＶ的ＨＤＴ１分析结果 从跳绳前、跳绳中、跳绳后的ＲＲ间期序列中 分别截取了连续的５１２点来进行ＨＤＴ１分析。这 里需要说明的是，对于跳绳前，考虑到部分受试者 在跳绳开始前５ ｍｉｎ内有散步活动，因此，跳绳前 的数据从开始采样后１０ ｍｉｎ左右为首截取。而跳 绳中数据则从跳绳开始即刻为首截取，跳绳后的数 据是从跳绳停止即刻为首截取的。 ６４０ 数 据 采 集 与 处 理 第２８卷 本文考察的嵌入维数ｍ为２～１６。当嵌入维 为８和１６时，各受试者的Ｄ 参数值如图１（ａ，ｂ） 所示。从图１可见，几乎所有的受试者，其对应３ 种状态下的Ｄ 值大小排序都有：跳绳中最大，跳 绳后次之，跳绳前最小。这种现象在嵌入维ｍ值 大的图１（ｂ）中比嵌入维值稍小的图１（ａ）中表现得 更明显。成对样本Ｔ检验的结果表明，嵌入维为 １６时，跳绳前的Ｄ 参数值显著小于跳绳中（Ｐ一 １０ ）和跳绳后（ 一８×１０ ），另外，跳绳中的Ｄ 参数值也显著大于跳绳后（Ｐ一８×１０ ），对于其 他嵌入维ｍ，上述显著性差异依然存在，只是Ｐ值 稍大而已。 进一步计算了跳绳后的５个阶段的Ｄ 参数 值，结果如图２所示。第１个阶段对应跳绳后的第 １个５１２点，在图２中编号为状态３；第２个阶段对 应跳绳后的第２个５１２点，在图２中编号为状态 ４－．．…・；第５个阶段对应跳绳后的第５个５１２点， 在图２中编号为状态７。为了方便对照比较，将跳 绳前（编号为状态１）和跳绳中（编号为状态２）的结 果也列于图２中。 受试者编号 （ａ）ｍ＝８ 受试者编号 （ｂ）ｍ＝１６ 图１ 不同状态下（跳绳前、中、后）ＲＲ间期序列ＨＤＴ１ 分析结果 测 状态编号 图２不同状态下所有受试者ＨＤＴ１分析结果均值 从图２可以看出，在进行ＨＤＴ１分析时，嵌入维值 取５，８，１６时，Ｄ 参数值从静息时的状态１到跳绳 中的状态２显著上升；跳绳结束后，在状态２～３和 状态３～４内，Ｄ 参数值迅速下降，随后下降趋势 趋于缓慢，逐渐收敛于静息时值；状态４，５，６，７的 Ｄ 值之间并无显著性差异。运动中交感神经持 续兴奋而迷走神经受抑制，运动后迷走神经张力逐 渐恢复而趋向交感迷走平衡的现象在强度较大的 体育运动中普遍存在。就每个健康个体而言，如果 视安静清醒时的交感迷走活性状态为平衡态，那么 在运动中和运动后的短时内，其交感迷走活性之间 应该都是非平衡的，且运动后的非平衡性会随着恢 复时间的增加而逐渐减弱直至消失，重新达到平衡 态。这里Ｄ 值的变化正好反映了跳绳运动前、 中、后交感迷走神经之间的这种从平衡态到非平衡 态、再从非平衡态逐渐恢复到平衡态的过程。 仔细观察图２就可以发现另一个值得深究的 现象：在嵌入维 一５，８，１６时，在跳绳前静息态 下，三类Ｄ 参数值基本混叠在一起，成对Ｔ检验 结果表明两两之间并无显著性差异；在跳绳后交感 迷走接近静息平衡态，即状态４，５，６，７下，三 类Ｄ 参数值也基本混叠在一起，成对Ｔ检验结果 表明两两之间仅Ｄ 和Ｄ。之间在状态４下有显著 性差异（Ｐ一０．０３）；而在状态２下，Ｄ 显著大于Ｄｓ （单边Ｔ检验，Ｐ一２×１０ 。），Ｄ 显著大于Ｄ （ 一 ２×１０ ）；状态３下也是如此，Ｄ 显著大于Ｄ。（单 边Ｔ检验，Ｐ一７×１０ ），Ｄ８显著大于Ｄｓ（Ｐ一１× １０ ）。这种现象表明，在某些特定的状态下，随着 嵌入维数ｍ的增加，Ｄ 参数值也随之单调递增。 为了探究这一点的本质来源，在图３中绘出了受试 者１１在状态１、状态２、状态３和状态７下的ＲＲ 间期波形图。 第５期 侯凤贞，等：运动前后心率变异性的高维时间不可逆性分析 ６４１ 心跳编号 图３受试者１１在不同状态下的ＲＲ间期波形 从图３可以直观地发现：状态１下ＲＲ间期波形最 “粗糙”，展示了最大的ＨＲＶ，ＲＲ间期均值最大 （心率最／ｂ）也最平稳；在跳绳中，也就是状态２最 初的１００点内，ＲＲ间期呈线性趋势急速下降（心 率快速增加），而随后ＲＲ间期值下降速度（心率增 加速度）明显减慢，波形较光滑，ＨＲＶ较小，这与 现有报导一致；在跳绳结束后的第一个状态（状态 ３）初期，ＲＲ间期呈线性趋势增加（心率减慢），随 后ＲＲ间期增加的速度减慢，波形也逐渐变得“粗 糙”起来，ＨＲＶ逐渐增大，这一点也与现有报导一 致；在状态７下的波形并未显示出明显的上升或下 降趋势，具有一定的心率变异特性，但还未达到静 息时的水平。 总体来说，与其他状态不同的是，在状态２和 状态３的ＲＲ间期变化过程中，都存在着明显的快 速、线性变化趋势的部分，尤以状态２更为明显，而 这两种状态下的Ｄ 值也明显高于静息时的值。 因此，将数据按如下规则继续划分：跳绳前２．５ ｍｉｎ 一段，共考虑了７段；跳绳中１ ｍｉｎ一段，共考虑了 ４段；跳绳后１．５ ｍｉｎ一段，共考虑了１８段。然后 对上述２９段序列（序列点数不相同），分别进行 ＨＤＴ１分析。同时将每段ＲＲ间期序列转换成瞬 时心率序列，并对其做一阶线性拟合，以ＨＲｓｌｏｐｅ 表示线性拟合时的斜率 图４显示了２９段序列的 Ｄ 指标值和ＨＲｓｌｏｐｅ指标的绝对值（１３位受试者 的平均值），图中虚线代表跳绳前７段的指标均值。 线性相关性检验表明，１３位受试者的Ｄ 指标 均与ＨＲｓｌｏｐｅ指标的绝对值有着显著的正相关性 （户＜０．０５）。可以说，运动中和运动后恢复期中的 Ｄ 指标值更多地受心率变化趋势的影响而不是 心率变异性的影响。上述结果表明，在将ＨＤＴ１ 分析方法应用于极短的时间序列（约１５０￣２００点） 时，该方法极易受数据中单调变化的线性趋势的影 状态编号 （ａ）Ｉ－ＩＤＴ１分析，嵌入维ｍ＝１６ 鸯 状态编号 （ｂ）心率变化趋势分析 图４所有受试者在不同状态下２９段短时ＲＲ间期序 列的分析结果均值 响，此时Ｄ 指标值并不能准确地反映出ＨＲＶ的 大小；但是如果该趋势并非噪声干扰，而是由于内 在动力系统发生快速变化所致，而Ｄ 反映出来的 信息恰好可以很好地捕捉到这种变化，此时认为不 应该在进行ＨＤＴ１分析方法前去除该线性趋势。 ３ 结束语 通过对１３个健康年轻人在静息时、跳绳运动 早、中、晚期及恢复早、中、晚期进行ＲＲ间期序列 的ＨＤＴ１分析，本文认为高维的时间不可逆性分 析指标Ｄ 值的变化正好反映了跳绳运动前、中、 后交感一迷走神经之间从平衡态到非平衡态、再从 非平衡态逐渐恢复到平衡态的过程。而事实上，对 于不同的受试个体，Ｄ 值在从静息到运动中的上 升程度与从运动中到停止运动后短时内的下降程 度是很不相同的，甚至有些受试者如１２号受试者， 其Ｄ 值在停止运动后短时内不降反升。有趣的 是，１２号受试者主观反映跳绳后有胸闷、心悸不适 感，而其他受试者都没有反映有这种主观感觉。因 此，有如下设想，能否以从运动中到停止运动后短 时内的Ｄ 下降水平，与从静息到运动中的Ｄ 上 升水平的比值来衡量运动后迷走神经的恢复能力， 该比值越大代表恢复能力越强。但是这个设想是 否可行，需要大量的实验来验证。这将是未来的一 个研究方向。 本文的实验结果也表明，在分析运动中心率快 速变化的数据时，Ｄ 与心率变化的趋势，即变化 率密切相关，那么能否将它应用于无氧阈的评定问 ６４２ 数 据 采 集 与 处 理 第２８卷 题，成为心率拐点法的一个补充甚至替代，也有待 大量实验来验证。 最后，上述实验结果表明，在应用ＨＤＴＩ方法 来研究不同病理状态下的ＨＲＶ时，数据长度不宜 太短，因为ＲＲ间期是通过计算体表ＥＣＧ信号中 相邻ＲＲ波顶点之间的时间间隔得到的，对于较长 的ＲＲ间期序列（如５ ０００点，对应１ ｈ左右的 ＥＣＧ采集时间），一般不会有一个单调的递增或者 递减变化趋势的存在，这样计算出来的Ｄ 指标值 才能可靠地反映出ＨＲＶ的大小。 参考文献： ［１］ Ｍａｌｉｋ Ｍ，Ｂｉｇｇｅｒ Ｊ Ｔ，Ｃａｒｏｍ Ａ Ｊ，ｅｔ ａ１．Ｈｅａｒｔ ｒａｔｅ ｖａｒｉａｂｉｌｉｔｙ：Ｓｔａｎｄａｒｄ ｏｆ ｍｅａｓｕｒｅｍｅｎｔ，ｐｈｙｓｉｏｌｏｇｉｃａｌ ｉｎｔｅｒｐｒｅｔａｔｉｏｎ，ａｎｄ ｃｌｉｎｉｃａｌ ｕｓｅ［Ｊ］．Ｅｕｒｏｐ Ｈｅａｒｔ Ｊ， １９９６，１７（３）：３５４—３８１． ［２］ Ｋｌｅｉｇｅｒ Ｒ Ｅ，Ｓｔｅｉｎ Ｐ Ｋ，Ｂｉｇｇｅｒ Ｊ Ｔ．Ｈｅａｒｔ ｒａｔｅ ｖａｒｉ— ａｂｉｌｉｔｙ：Ｍｅａｓｕｒｅｍｅｎｔ ａｎｄ ｃｌｉｎｉｃａｌ ｕｔｉｌｉｔｙ［Ｊ］．Ａｎｎ Ｎｏｎｉｎｖａｓｉｖｅ Ｅｌｅｃｔｒｏｃａｒｄｉｏｌ，２００５，１Ｏ（１）：８８—１Ｏ１． ［３］ Ｔｈａｙｅｒ Ｊ Ｆ，Ａｈｓ Ｆ，Ｆｒｅｄｒｉｋｓｏｎ Ｍ，ｅｔ ａ１．Ａ ｍｅｔａ— ａｎａｌｙｓｉｓ ｏｆ ｈｅａｒｔ ｒａｔｅ ｖａｒｉａｂｉｌｉｔｙ ａｎｄ ｎｅｕｒｏｉｍａｇｉｎｇ ｓｔｕｄｉｅｓ：Ｉｍｐｌｉｃａｔｉｏｎｓ ｆｏｒ ｈｅａｒｔ ｒａｔｅ ｖａｒｉａｂｉｌｉｔｙ ａｓ ａ ｍａｒｋｅｒ ｏｆ ｓｔｒｅｓｓ ａｎｄ ｈｅａｌｔｈ［Ｊ］．Ｎｅｕｒｏｓｃｉ Ｂｉｏｂｅｈａｖ Ｒｅｖ，２０１２，３６：７４７－７５６． ［４］ Ｎｏｒｍａｎ Ｇ ｊ，Ｋａｒｅｌｉｎａ Ｋ，ｇｅｒｎｔｓｏｎ Ｇ Ｇ，ｅｔ ａ１． Ｈｅａｒｔ ｒａｔｅ ｖａｒｉａｂｉｌｉｔｙ ｐｒｅｄｉｃｔｓ ｃｅｌｌ ｄｅａｔｈ ａｎｄ ｉｎｆｌａｍ— ｍａｔｏｒｙ ｒｅｓｐｏｎｓｅｓ ｔｏ ｇｌｏｂａｌ ｃｅｒｅｂｒａｌ ｉｓｃｈｅｍｉａ［Ｊ］． Ｆｒｏｎｔ Ｐｈｙｓｉｏｌ，２０１２，３：１３１． Ｅ５］ Ｃｈａｎｄｒａ Ｐ，Ｓａｎｄｓ Ｒ Ｌ，Ｇｉｌｌｅｓｐｉｅ Ｂ Ｗ，ｅｔ ａ１．Ｐｒｅｄｉｃ— ｔｏｒｓ ｏｆ ｈｅａｒｔ ｒａｔｅ ｖａｒｉａｂｉｌｉｔｙ ａｎｄ ｉｔｓ ｐｒｏｇｎｏｓｔｉｃ ｓｉｇｎｉｆｉ— ｃａｎｃｅ ｉｎ ｃｈｒｏｎｉｃ ｋｉｄｎｅｙ ｄｉｓｅａｓｅ［Ｊ］．Ｎｅｐｈｒｏｌ Ｄｉａｌ Ｔｒａｎｓｐｌａｎｔ，２０１２，２７（２）：７００—７０９． ［６］ 宋淑华，刘坚，高春刚．ＨＲＶ指标在体育研究领域中 的应用（综述）［Ｊ］．山西师大体育学院学报，２０１０，２５ （５）：１２５—１２８． Ｓｏｎｇ Ｓｈｕｈｕａ，Ｌｉｕ Ｊｉａｎ，Ｇａｏ Ｃｈｕｎｇａｎｇ．Ａｐｐｌｉｃａｔｉｏｎ ｏｎ ｈｅａｒｔ ｒａｔｅ ｖａｒｉａｂｉｌｉｔｙ ｉｎ ｓｐｏｒｔｓ ｓｃｉｅｎｔｉｆｉｃ ｒｅｓｅａｒｃｈ （ｒｅｖｉｅｗ）口］．Ｊｏｕｒｎａｌ ｏｆ Ｐｈｙｓｉｃａｌ Ｅｄｕｃａｔｉｏｎ Ｉｎｓｔｉｔｕｔｅ ｏｆ Ｓｈａｎｘｉ Ｎｏｒｍａｌ Ｕｎｉｖｅｒｓｉｔｙ，２Ｏ１０，２５（５）：１２５— １２８． １－７］ 杜吟，李京诚．ＨＲＶ在运动领域应用研究的现状ＥＪ］。 首都体育学院学报，２０１１，２３（１）：８９—９６． Ｄｕ Ｙｉｎ，Ｌｉ Ｊｉｎｇｃｈｅｎｇ．Ｒｅｖｉｅｗ ｏｎ ｃｕｒｒｅｎｔ ＨＲＶ ａｐ— ｐｌｉｅｄ ｒｅｓｅａｒｃｈ ｉｎ ｓｐｏｒｔｓ ｆｉｅｌｄ［Ｊ］．Ｊｏｕｒｎａｌ ｏｆ Ｃａｐｉｔａｌ Ｉｎｓｔｉｔｕｔｅ ｏｆ Ｐｈｙｓｉｃａｌ Ｅｄｕｃａｔｉｏｎ，２０１ １，２３（１）：８９—９６． Ｅｓ］ 张立．Ｃ０ＮＣ０ＮＩ测试对足球运动员有氧能力评定的 可行性研究［Ｊ］．武汉体育学院学报，２００１，３５（２）：８３— ８７． Ｚｈａｎｇ Ｌｉ．Ｒｅｓｅａｒｃｈ ｏｎ ｆｅａｓｉｂｉｌｉｔｙ ｏｆ ｕｓｉｎｇ ＣＯＮＣＯＮＩ ｔｅｓｔ ｉｎ ａｓｓｅｓｓｉｎｇ ａｅｒｏｂｉｃ ｃａｐａｃｉｔｙ ｏｆ ｆｏｏｔｂａｌｌｅｒ［Ｊ］． Ｊｏｕｒｎａｌ ｏｆ Ｗｕｈａｎ Ｉｎｓｔｉｔｕｔｅ ｏｆ Ｐｈｙｓｉｃａｌ Ｅｄｕｃａｔｉｏｎ， ２００１，３５（２）：８３－８７． ［９］ Ｂｕｃｈｈｅｉｔ Ｍ，Ｓｏｌａｎｏ Ｒ，Ｍｉｌｌｅｔ Ｇ Ｐ．Ｈｅａｒｔ—ｒａｔｅ ｄｅ— ｆｌｅｃｔｉｏｎ ｐｏｉｎｔ ａｎｄ ｔｈｅ ｓｅｃｏｎｄ ｈｅａｒｔ－ｒａｔｅ ｖａｒｉａｂｉｌｉｔｙ ｔｈｒｅｓｈｏｌｄ ｄｕｒｉｎｇ ｒｕｎｎｉｎｇ ｅｘｅｒｃｉｓｅ ｉｎ ｔｒａｉｎｅｄ ｂｏｙｓＥＪ］． Ｐｅｄｉａｔｒ Ｅｘｅｒｃ Ｓｃｉ，２００７，１９（２）：１９２—２０４． ［１Ｏ］Ｈｏｕ Ｆｅｎｇｚｈｅｎ，Ｎｉｎｇ Ｘｉｎｂａｏ，Ｚｈｕａｎｇ Ｊｉａｎｊｕｎ，ｅｔ ａ１． Ｈｉｇｈ—．ｄｉｍｅｎｓｉｏｎａｌ ｔｉｍｅ ｉｒｒｅｖｅｒｓｉｂｉｌｉｔｙ ａｎａｌｙｓｉｓ ｏｆ ｈｕ— ｍａｎ ｉｎｔｅｒｈｅａｔ ｉｎｔｅｒｖａｌｓ ＥＪ］．Ｍｅｄ Ｅｎｇ Ｐｈｙｓ，２０１１， ３３：６３３—６３７． 作者简介：侯风贞（１９８０一），女，博士，讲师，研究方向：生物医 学信号处理；宁新宝（１９４１一），男，教授，研究方向：生物医学 电子学及生物医学信号处理；黄晓林（１９７７一），女，讲师，研究 方向：生物医学信号处理，Ｅ－ｍａｉｌ：ｘｌｈｕａｎｇ＠ｎｊｕ．ｅｄｕ．ｃｎ。
