感知与适应16 分钟

心率变异性来自哪里,能够反映哪些身体状态?

HRV描述正常心搏间隔怎样变化,但SDNN、RMSSD与频谱并不读取同一种变化。本文沿呼吸、自主神经、起搏工作点与电/光测量链解释这些指标,计算同分布不同顺序、增益与时标误差的最小例子,并用2025年原始补充数据复算指端和耳内设备的四组一致性结果。近期临床、运动员与训练研究表明,测量一致、身体状态推断与干预获益必须逐步建立联系;单个数值不能直接换算为压力、恢复或健康余额。

心跳每分钟相同,节拍仍可以不同

假设两段记录都包含300个正常心搏间隔,其中150个为980毫秒,另外150个为1020毫秒。第一段短、长交替;第二段先连续出现短间隔,再连续出现长间隔。两段的平均间隔相同,每种数值出现的次数也相同,但第一段每一搏都在改变,第二段只有中间转换一次。

通常所说的心率变异性,简称HRV,主要描述正常窦性心搏间隔的变化。它不是一个单独指标,更不只是“心跳不齐程度”。上述两段记录计算SDNN都是20.03毫秒,计算RMSSD却分别为40.00和2.31毫秒。

同样的间隔集合,顺序不同导致RMSSD不同
同样的间隔集合,顺序不同导致RMSSD不同

本库按定义计算的300搏示例,不是患者数据。SDNN的样本方差使用N−1作分母。

这个小例子已经说明:讨论“HRV升高”之前,需要知道究竟哪个指标升高、记录了多久、其中哪些心搏被纳入。少掉这些条件,同一个词可能在描述完全不同的现象。

不同指标分别保留什么

把相邻正常心搏之间的间隔记作$N_i$,平均值记作$\bar N$。常见的两个时间域指标为:

$$ \begin{aligned} \mathrm{SDNN}&=\sqrt{\frac{1}{M-1}\sum_{i=1}^{M}(N_i-\bar N)^2},\\ \mathrm{RMSSD}&=\sqrt{\frac{1}{M-1}\sum_{i=1}^{M-1}(N_{i+1}-N_i)^2}. \end{aligned} $$

SDNN看整段间隔围绕均值分散得多开,RMSSD看相邻间隔变了多少。前者不关心重排顺序,后者非常关心。把记录从五分钟延长到一天,又会纳入睡眠、活动与姿势等较慢变化,因此不同记录长度的SDNN不能直接当成同一把尺子。1

频率域指标换了一个问题:间隔变化中,多慢或多快的起伏各占多少功率?一个常见成人分析方案把0.04—0.15 Hz称为低频LF,把0.15—0.40 Hz称为高频HF。文献的下限和边界并不完全一致,例如某些生理实验采用更窄的频带。频带本身是分析约定,不能从名字直接推到交感或副交感神经。

还要分清“间隔中的高频”和“心电波形中的高频”。前者是在许多次心搏之间观察节拍调制,后者可能是在毫秒尺度观察波形细节;它们不是同一种频谱。

这里的“正常”同样有实质含义。漏检、重复计数、异位搏动与某些心律失常都能使原始RR间隔更分散。它们可以成为另一项心律分析的信号,但不能不作区分地并入自主神经HRV解释。插值和剔除能够处理部分问题,也会改变正在计算的变异。1

间隔为什么会随着呼吸和神经调节变化

心脏不是等待大脑逐次发出“现在跳一下”的指令。起搏组织自身持续活动,自主神经改变它的工作状态,呼吸、血压调节和其他过程又不断改变这些输入。最终的心搏间隔,是多条过程经过起搏系统后的输出。

其中一条比较容易理解的链是呼吸性窦性心律不齐,简称RSA。呼吸相关中枢活动、胸腔和血流变化、压力感受反射等环节,会共同改变心脏所受的调节。在常见静息条件下,心搏间隔随呼吸周期起伏,迷走神经介导的快速作用是其中一个关键环节。1

副交感和交感作用的时间响应也不同。前者能够较快改变起搏过程,交感作用包含较慢的受体和细胞内传递环节。因此,不同频率的输入到达心脏后,并不会以相同幅度保留。系统有自己的响应速度和增益。

这一解释有干预证据。Pomeranz等1985年研究了8名健康男性,其中6名完成第二天的逆序阻断比较。受试者按固定呼吸节律接受卧位、站立和自主神经阻断条件的比较。阻断副交感作用后,呼吸附近的高频功率大幅降低;低频的变化却依赖姿势:卧位与站立不能用同一种“低频等于交感”的解释概括。2

这项实验比单纯相关更能说明神经通路参与,但它仍是特定人群、姿势、呼吸频率和药理条件下的结果。不能把阻断后的百分比当成所有人的固定神经贡献比例,也不能将两个分支的作用简单相加。原研究分析的是瞬时心率功率谱,本文不把它的数值直接搬成毫秒间隔功率。

因此,较快的HRV变化在合适条件下能提供心脏副交感调节的信息;“当前迷走神经平均放电有多强”则是更进一步的问题。呼吸同步的调制幅度、平均神经活动与心脏受体的响应程度,并不天然一一对应。1

呼吸变慢,频带的解释也会改变

假如某种间隔起伏每4秒重复一次,频率是0.25 Hz,落入上述HF带;若改成每10秒重复一次,就变成0.10 Hz,落入LF带。仅这一步,已经可能改变LF/HF的比值,尚不需要假定交感和副交感的平均作用发生了相应的此消彼长。

Brown等1993年的9人实验也提供了具体依据:改变呼吸频率与潮气量会影响间隔功率,而平均间隔未出现相应的显著变化。本次仅取得该研究的原始摘要,因此只用它支持这一有限结论,不重建其完整实验与个体数据。3

尤其需要注意,RSA是按“与呼吸有关”定义的生理现象,HF是按频率边界定义的量。慢呼吸时,两者可以分开;说话改变呼吸方式,也可能同时改变HRV。如果一个压力实验把安静坐着与站着讲话比较,观察到的差异就混合了心理任务、姿势与呼吸变化。1

RMSSD也不是自动消除这些影响的办法。相邻差分对慢变化和快变化的响应不同;同样幅度的周期变化,换一个周期长度,就可能得到不同RMSSD。它通常比固定频带少一道“峰值跨出边界”的问题,但仍需要解释输入和记录条件。

同样的扰动,起搏系统可以放大出不同波动

另一层关系发生在心脏自身。Monfredi等2014年比较了人体、动物、离体心脏和单个窦房结细胞,讨论平均心率与间隔变异之间的生物物理关系。离体和单细胞记录中的变异,说明全部波动不能都归给完整身体内的自主神经输入;其大小也与起搏系统所处的工作状态有关。4

可以用一个很小的模型看清这一步。设起搏细胞从一个电位升到触发点需要补充固定电荷$Q$,平均向内电流的大小为$J$,每搏叠加扰动$\eta_i$;动作电位和其他固定阶段合并为$A$。则:

$$ N_i=A+\frac{Q}{J+\eta_i}. $$

这是假定爬升速度恒定、每搏状态简化到一个数的示意模型。扰动小时,间隔波动的尺度约为:

$$ \sigma_N\approx\frac{Q}{J^2}\sigma_\eta. $$

关键在$J^2$这个分母:相同电流扰动经过不同工作点,产生的时间变化可以不同。本库取自选的固定参数与同一组20000次随机扰动,将平均电流依次设为0.4、0.6、0.8 pA,间隔标准差得到50.27、22.21、12.47毫秒。由平均间隔换算的心率同时约为51.64、72.54、90.89次/分。

工作点改变波动增益,频率改变功率所属频带
工作点改变波动增益,频率改变功率所属频带

上图是恒定爬升模型及小扰动近似;下图是同幅度的均匀采样正弦间隔示意,仅用于展示频谱,不从这条重采样曲线计算RMSSD。两者均非人体拟合或原论文完整模型复现。

这个例子不说明人类HRV变化都能由平均心率解释。真实起搏过程包含多种离子电流、钙活动、非线性响应与相互作用,同一外部输入也未必对应模型中同一电流噪声。它说明的是:输出变异不仅依赖输入波动,还依赖把输入转换成心搏时间的系统。

因而,“先把心率影响全部回归掉,剩下的就是真实自主神经状态”也不稳妥。平均心率与HRV可能共同受神经调节影响,机械地去掉两者共享的变化,可能连想研究的部分一起去掉。2024年指导综述建议联合考察心搏间隔与RSA,并明确讨论校正假设,而非默认一种校正能还原神经真值。1

光学脉搏比心电晚到,不一定导致HRV不同

上一文讨论过,心电和光学脉搏经历不同传递过程。这里可以把逐搏时间关系写得更精确。

设第$n$次心电事件发生在$e_n$,对应光学标记发生在$p_n$,两者间的总延迟为$d_n$。这个延迟合并了电活动到收缩、血流和脉搏传播、测点波形及仪器处理等环节,不直接等同某一段血管的传导时间。对正确配对、共用时间尺度的事件,有:

$$ \begin{aligned} p_n&=e_n+d_n,\\ P_n&=R_n+(d_n-d_{n-1}). \end{aligned} $$

$P_n$和$R_n$分别是相邻光学标记与相邻心电事件的间隔。如果每一搏都晚200毫秒,固定延迟会相消,两组间隔完全相同。真正进入差异的是延迟怎样逐搏变化,以及检测和处理怎样改变标记。

变化的延迟既可能增加变异,也可能削弱变异。按方差关系:

$$ \operatorname{Var}(P)=\operatorname{Var}(R) +\operatorname{Var}(\Delta d) +2\operatorname{Cov}(R,\Delta d). $$

若最后的协方差为负,增加的传递过程并不必然增加最后的间隔方差。仅凭“PPG经过更多环节”不能推出固定的误差方向。

本库做了两个对照。第一例的真实间隔始终为1000毫秒,每个光学事件却带有标准差10毫秒的独立时标误差,计算出的RMSSD约为24.46毫秒。第二例特意构造延迟与原间隔变化的负向关系,得到的RMSSD恰好减半。

时标误差增加变异与特定延迟协变削弱变异的示意
时标误差增加变异与特定延迟协变削弱变异的示意

两个例子均为数学构造,10毫秒不是某款设备的实测误差,减半也不是某项实测低估的已证原因。第一例可解析得到RMSSD的总体尺度为$\sqrt{6}\times10$毫秒;数值结果与之相符。

这还说明,光学波形看起来更平滑,不足以证明逐搏间隔变异一定更小。波形滤波、峰位提取和间隔序列处理是相关但不同的步骤,要核对实际算法和数据。

用原始补充数据,看“相关很高”距离“数值相同”有多远

Parikh等2025年让30名参与者在尽量静止的条件下,分别完成自然呼吸和深慢呼吸记录,同时使用心电、指端与耳内设备。每种条件记录十分钟,并采用特定异常值处理流程。作者提供了处理后的逐人指标和排除标记。5

本库从这些数据重新计算RMSSD的一致性分析。按原排除标记,各比较保留27—30人;表中差值统一为心电RMSSD减去光学RMSSD。

设备与条件 纳入人数 平均差值(ms) 一致性界限(ms)
指端,自然呼吸 29 −0.47 −9.17至8.23
指端,深慢呼吸 30 0.36 −5.47至6.19
耳内,自然呼吸 28 −4.20 −19.42至11.02
耳内,深慢呼吸 27 −1.53 −14.46至11.40
按公开逐人处理结果重绘的四组一致性图
按公开逐人处理结果重绘的四组一致性图

每个点是一人在该条件下两种设备的十分钟指标比较。同一参与者可出现在多幅图,不把114个比较当成114名独立参与者。虚线为平均差值加减1.96倍差值样本标准差,是一致性界限,不是均值的置信区间。

四组结果与原表的舍入数值一致。这里可以直观看到:平均偏差接近零,并不意味着每个人误差都接近零。两种设备也可能把人与人的高低排序排得很像,仍留下不可忽略的个体差异。

按作者选用的相对一致性界限和30%判据,指端设备两种条件下的RMSSD达到其“可接受”标准,耳内未达到。这是对指定设备、十分钟记录和处理流程的有条件支持。该阈值不是一切医疗或训练用途的统一合格线;指端设备的这一结果也不直接证明手腕设备、短时读数或剧烈运动下同样可靠。5

本文取得的是已处理的逐人结果,重算验证了汇总与图表,没有从原始光学和心电波形重做峰值检测、插值与剔除。源头处理会不会抹掉有意义的变化,仍需要另一层验证。

近期的大样本研究增加了什么,又没有增加什么

2025年另一项临床研究在931名就诊者静坐记录中比较了胸部心电、上臂心电与上臂PPG。组平均心率接近,而上臂PPG得到的RMSSD与SDNN低于心电结果:例如胸部心电和上臂PPG的平均RMSSD分别为43.10与37.49毫秒,SDNN分别为77.23与64.10毫秒。6

它提供了具体配置下的重要负向验证,但不与前述指端结果构成逻辑矛盾:测点、设备、采样、算法、记录长度和人群都不同。组均值有差异,也没有单独给出任何用途都无法使用的结论。就诊者的疾病分组还可能重叠,药物等因素未充分控制,不能把组间差异直接归给单一疾病。

同一技术研究团队2026年发表的运动员研究,记录103名男性大学橄榄球运动员、三个赛季和12726组同步数据。其上臂PPG与心电的RMSSD相关系数达到约0.953,同时仍有明显的数值偏差。大量记录来自重复测量,不能作为12726名独立运动员理解。7

这篇还有一个值得细读的环节:作者比较两种读数识别“状态恶化”的速度,而所用参照事件,是心电HRV低于其个人滚动基线的某个标准化阈值。研究没有用独立测得的神经功能损害或运动损伤来定义这次恶化。因此它检验了不同读数追踪HRV阈值事件的差别,尚不能把结果直接换成“提前若干天发现真实自主神经损伤”。方法与图注对延迟截断天数的说明也有差异,本文不复述精确提前天数。

此外,两种序列各用自己的滚动均值和标准差时,固定线性偏差会被标准化消去。对正比例系数$a$,同一窗口有:

$$ z(aX+b)=\frac{aX+b-(a\mu+b)}{a\sigma}=z(X). $$

本库用滚动窗口实际核对了这个恒等关系。因此,若标准化阈值的跨越时间不同,仅有固定低估或固定偏移还解释不了;需要继续检查时变关系、噪声、清理与窗口处理。这个数学限定不否认原数据存在差异,而是限定差异的机制解释。

HRV到底能不能帮助判断压力、恢复与健康

可以提供信息,但要先说明判断的对象和证据层次。

判断一段生理调节是否发生变化,通常比判断“这个人有多少压力”更接近HRV直接支持的范围。若同一个人在相近时间、姿势、活动和测量条件下反复记录,设备与处理也稳定,持续偏离个人通常范围就值得解释。不过,它仍可能涉及呼吸、睡眠、训练负荷、疾病、药物或其他变化;HRV本身未把这些来源逐一分开。1

预测未来结局不要求指标只来自一种机制。Framingham的早期研究在736名合格老年参与者中,发现部分HRV指标与随后四年死亡风险存在调整后的统计关联。本次仅核到原始摘要,不评价完整模型的全部协变量。但这类研究说明HRV可以包含预后信息,同时完全不等于“主动把某个HRV数字提高,就能按同一比例降低风险”。8

把指标用于调整训练则需要比较实际策略效果。2022年一项低风险冠心病康复试验随机分配23人,21人完成六周方案。采用HRV指导与预定训练的两组,在周平均RMSSD变化上出现差异;其余所分析指标,包括最大摄氧量的变化,没有检测到组间显著差异。9

这个结果可以支持继续研究有条件的策略用途,却不能当成减少死亡的临床证明。它的日常指导指标来自手机PPG,实验室单次指标来自胸带,两种记录本就不同;样本小,而且完成者中的5名糖尿病患者全在预定训练组。没有显著差异不证明两策略等效,试验标题中的“死亡预测指标”也不是实际死亡结局。

所以,一种指标即使不完美,也可能有决策价值;但价值要由它在指定任务中增加了什么来证明。真正有区分力的后续研究应让状态参照尽量独立于HRV自身,在未参与拟合的人或时段检验,并比较加入HRV后是否优于已有的心率、活动、睡眠或自我感受信息。若声称干预获益,还要观察独立于该指标的功能或健康结局。

把读数放回产生它的过程

目前较稳的认识是:正常间隔变异包含自主神经与呼吸调节的信息,但也受起搏工作点、传递过程和数据处理影响;不同HRV指标和光学脉搏变异不能无条件互换。药理实验、设备验证与实际计算分别支撑这些环节。

本文的整合判断是,HRV更适合作为在已说明条件下观察身体调节的一个窗口。要把它解释成压力、疲劳或恢复,需要相应状态的独立依据;要将它用于选择行动,还要看实际策略是否改善了目标结局。没有一个通用的“越高越健康”刻度能够代替这条证据链。

仍未解决的是,哪些长期个体模式能稳定区分不同状态,何种设备和处理能保留这些模式,以及这种额外信息究竟能带来多少实际帮助。上述研究没有提供跨设备、跨人群、跨任务的统一答案。

本篇的复算脚本保存在仓库中,结果记录保存模型参数、数值输出与原数据重算口径。生理过程解释与自选模型明确分开,本文不构成个人诊断或训练处方。

来源与依据

  1. S1 领域指导综述;用于机制与测量概念,不是本库投稿验收规范

    Publication guidelines for human heart rate and heart rate variability studies in psychophysiology

    Karen S Quigley; Peter J Gianaros; Greg J Norman; J Richard Jennings; Gary G Berntson; Eco J C de Geus · 2024

    · 相关全文已核

    §3相关平均量定义,§4心脏神经/呼吸/姿势/工作点与§5.2、Table1指标及处理段

    芝加哥机构库公开原PDF;2024发表,2026机构记录更新不当新论文。实际读相关段,不称63页全部精读。HF与RSA、相位调制与平均迷走活动分开;不继承固定LF/HF神经平衡或通用心率残差校正。

  2. S2 原始人体阻断与姿势实验

    Assessment of autonomic function in humans by heart rate spectral analysis

    Bruce Pomeranz et al. · 1985

    · 相关全文已核

    H151–H153全文,研究方法、Table1及讨论

    原PDF公开副本来自psicolibra.it并核原发表题录。8名22–36岁健康男性,第二天6名完成;固定呼吸,瞬时心率谱0.04–0.12与0.224–0.28Hz。不能直接混同NN功率或普遍神经贡献百分比;未复现生理实验。

  3. S3 原始呼吸条件实验的摘要

    Important influence of respiration on human R-R interval power spectra is largely ignored

    T E Brown; L A Beightol; J Koh; D L Eckberg · 1993

    · 摘要已核

    Europe PMC原始摘要与PMID8307890元数据

    9名健康参与者,呼吸频率/潮气量影响谱,平均RR未显著变化。APS PDF403、直接PubMed验证页并非摘要,已改名failed;不声称读过正文或重算数据。

  4. S4 原始跨准备实验与生物物理模型

    BIOPHYSICAL CHARACTERISATION OF THE UNDER-APPRECIATED AND IMPORTANT RELATIONSHIP BETWEEN HEART RATE VARIABILITY AND HEART RATE

    Oliver Monfredi; Alexey E Lyashkov; Anne-Berit Johnsen; Shin Inada; Heiko Schneider; Ruoxi Wang; Mahesh Nirmalan; Ulrik Wisloff; Victor A Maltsev; Edward G Lakatta; Henggui Zhang; Mark R Boyett · 2014

    · 相关全文已核

    Methods;Results人体/离体/细胞与简化模型;Discussion;PMC4326239

    人体11、活体大鼠11、兔离体心脏58、大鼠离体8、兔SAN细胞67等不同层级,不能当独立同质人群汇总。主张工作点与输入共同决定输出;不继承普遍指数规律或全部HRV都由HR解释。本库恒定爬升例子受其思路启发,非重现全文ODE或拟合参数;HTML扁平公式不作符号指控。

  5. S5 原始同步设备验证与公开处理后逐人数据

    Evaluation of In-Ear and Fingertip-Based Photoplethysmography Sensors for Measuring Cardiac Vagal Tone Relevant Heart Rate Variability Parameters

    Ankit Parikh; Gwyn Lewis; Hamid GholamHosseini; Usman Rashid; David Rice; Faisal Almesfer · 2025

    · 相关全文已核

    Methods处理/排除/统计,Table1–3,Discussion/Limitations;DataS1.xlsx、DataS2.xlsx、DataS4.txt

    30人,每条件10min。S1的180行含3设备×2条件,按discard配对保留指端29/30、耳内28/27,共114对,不当114人。本库从S1重算四组RMSSD均差、样本SD一致界限,并对照S2差值与S4排除;与原表至0.01ms一致。LOA%按最大绝对界限/全部配对均值重现23.27/14.20/47.22/33.49;不是逐人的相对误差或统一合格线。未重跑原波形清洗/检测或ICC,非神经调节/刺激疗效验证。主文一句“no data excluded”与前后计数冲突,以S1/S4和表中实际人数。其余附件取得未逐项核;补充ZIP527717字节、SHA记录。

  6. S6 原始临床ECG与上臂PPG同步比较

    Pulse rate variability is not the same as heart rate variability: findings from a large, diverse clinical population study

    Kantrowitz, Allen B. ; Ben-David, Kfir ; Morris, Michael ; Wittels, Harrison L. ; Wishon, Michael J. ; McDonald, Samantha M. ; Renaghan, Eric J. ; Feigenbaum, Luis A. ; Wittels, S. Howard · 2025

    · 相关全文已核

    Methods与Table1/各病组及Discussion/Limitations

    2024-05至09单中心931就诊者,静坐5–7min,上臂PPG与上臂/胸部ECG。主文Table1平均值用于描述;不自行重算配对CI,不继承所有PPG无用、光学低通必削弱间隔变异或错误神经递质描述。pNN50表头ms不合定义、年龄上下限/成人措辞有差异,本篇均不复述。疾病分组重叠,药物等缺口保留。原件完整,首次article节点仅导航的抽取已改为全文,不能用初次2段抽取冒充全文。

  7. S7 原始运动员重复同步记录与阈值比较

    Pulse Rate Variability Is Not the Same as Heart Rate Variability: Implications for Sports Performance and Injury Prevention

    Kfir Ben-David; Allen B Kantrowitz; Michael Morris; Eric J Renaghan; Luis A Feigenbaum; Kyle Bellamy; Joe Girardi; Harrison L Wittels; Michael J Wishon; Samantha M McDonald; S Howard Wittels · 2026

    · 相关全文已核

    Methods采集/处理/21日阈值,Table3–4,Fig2–4图注,Discussion/Limitations

    2026-06-12发表,年度合集Dec不是将来发表。103名男性单队,2022–24三个27周赛季、398session、12726记录。红光RMSSD r=.953与bias−23.004ms不当一致。ECG与PRV各用自己滚动基线,ECG阈值定义参照事件,不是独立神经损害/运动伤。Methods延迟封顶10天、Fig3写7天;检测比例的正文/摘要也不一致,不重复精确天数或泛化伤害预防。重复记录CI未确认处理个体聚类,不作新推断;race表重叠计数不归因。与s6同技术团队非独立设备重复,研究人群不同。

  8. S8 原始前瞻性队列摘要

    Reduced heart rate variability and mortality risk in an elderly cohort. The Framingham Heart Study

    H Tsuji et al. · 1994

    · 摘要已核

    Europe PMC原始摘要与PMID8044959元数据

    736名合格者均龄72,前2h动态ECG,4y74死亡;调整关联仅据原摘要,不列未经核对全部协变量或断言独立于HR。全文未得,不能由观察关联推出提高HRV降低死亡。

  9. S9 原始小样本随机康复策略比较

    Heart Rate Variability-Guided Training for Improving Mortality Predictors in Patients with Coronary Artery Disease

    Agustín Manresa-Rocamora; José Manuel Sarabia; Silvia Guillen-Garcia; Patricio Pérez-Berbel; Beatriz Miralles-Vicedo; Enrique Roche; Néstor Vicente-Salar; Manuel Moya-Ramón · 2022

    · 相关全文已核

    Methods日常PPG/实验室胸带;Results流程、Table1–2;Discussion

    2018–19实施,2022发表。25招募、23随机、21完成10指导11预定,不继承摘要21随机。周平均RMSSD来自手机PPG;实验室单次来自PolarH7,未混为ECG金标准。同一日常指标参与指导与部分结局,周均变化p=.034、其他组间未显著,不当等效或死亡结局。5糖尿病均在预定组,作者以部分糖尿病者改善排除混杂的解释未接受。未取SI或重分析原个体临床数据。

最后更新 2026-09-11

来源