心血管-自主神经（CVS-ANS）负反馈控制动力学的系统理论建模与稳态极限环分析：电路等效及状态空间方法

《Biomedical Signal Processing and Control》：Modeling and steady-state limit cycles of the cardiovascular-autonomic negative feedback control dynamics: A system theoretic approach

　　

心血管-自主神经（CVS-ANS）负反馈控制动力学涉及交感与副交感神经分支对心血管系统的快速调节，并在压力反射（baroreflex）与时滞非线性作用下产生有界振荡。本文提出一种面向系统理论的闭环模型：将心血管系统（CVS）分解为肺循环（PC）、左心室（LV）、体循环（SC）与右心室（RV）等模块，并采用电等效表示；心室收缩活动由非线性主动电压发生器及其收缩控制器（CVC）描述，心脏瓣膜由电压控制开关刻画阈值非线性。自主神经系统（ANS）通过压力反射的传入转导、静态非线性和具有传输延迟的传出通路调节心脏周期（T_heart）。研究人员构建了17^th阶状态空间模型，并在Simulink/MATLAB中仿真主要状态变量的时域演化，同时以相图进行数值稳定性分析。结果表明，在生理条件下系统轨迹收敛至稳定的极限环（limit cycle）；系统动脉压、心脏周期及自主传出变量呈现持续且有界的周期性振荡，邻近初始条件的轨迹收敛至同一闭合轨道。该研究提供关于极限环存在性、唯一性与稳定性的证据，并为CVS-ANS非线性动力学、稳态稳定性分析及病理情景预测提供了面向系统理论的in silico（计算机仿真）框架。

研究背景与问题。机体稳态调控由中枢神经系统、自主神经系统（ANS）和内分泌系统共同参与。ANS通过交感与副交感（迷走）分支调节非随意器官，并以交感-迷走平衡影响心率、心肌收缩性及血管阻力，使血压与心率在有界范围内稳定振荡。由于存在外源扰动和内源扰动，ANS仍能维持心血管系统（CVS）的稳态行为，因此定量评估交感-迷走作用对临床与生物医学工程均具重要意义。现有研究包括侵入性药理方法、以心率变异性（HRV）为间接指标的自主神经试验以及非侵入性数学模型；前两类方法常受解释复杂、个体差异和标准化不足限制。数学模型可研究调控过程及神经通路贡献，但既有CVS模型多采用液压表示并将CVS聚合成单一宏观模块，虽能再现总体血流动力学趋势，却掩盖分室动态，难以刻画心室强非线性；已有电等效模型也常依赖简化电路，难以准确表示主动、非线性心室。CVS-ANS的稳定生理振荡主要源于快速负反馈（以压力反射baroreflex为核心）和时滞（如ANS交感、副交感传递延迟），可用极限环（limit cycle）吸引子描述；这要求采用非线性动力学进行分析。论文因此旨在建立系统理论导向的模块化电等效模型，对CVS-ANS闭环调控进行状态空间建模与稳态稳定性分析。

研究人员开展的工作及结论意义。研究人员将CVS分解为PC、LV、SC、RV和CVC模块，采用定制电等效：心室收缩活动由非线性主动电压发生器表示，并由整合ANS信息的收缩控制器逐周期调控；心脏瓣膜由电压控制开关表示阈值驱动非线性；心房、体循环和肺循环按各自动态描述。进一步建立传入压力感受器动态、压力-放电频率转导、ANS的交感/副交感静态非线性及传出传输延迟，形成闭环17^th阶模型。仿真显示系统在稳态趋于稳定极限环，主要血流动力学变量在生理范围内周期波动，系统动脉压（x₆）约在80 mmHg舒张压至120 mmHg收缩压之间变化；T_heart总体振荡范围约46 ms，x₆范围约38 mmHg，得到范围型压力-周期增益约1.2 ms/mmHg。相空间结果表明，来自不同邻近初始条件的轨迹在瞬态后收敛至同一闭合轨道，提供稳定极限环局部吸引性的数值证据。其意义在于：该工作首次对完整CVS-ANS模型进行显式稳定性分析，说明闭环压力反射、时滞与非线性可产生自持周期节律，并为病理情景预测、参数探索及自主/循环功能障碍机制研究提供in silico平台。

关键技术方法（不超过250字）。电等效建模（压力对应电压、流量对应电流，采用电阻-电感-电容（RLC）网络表示阻容感）；将CVS拆为PC、LV、SC、RV并以CVC闭合；以状态空间形式建立闭环17^th阶模型；以Hill、指数和sigmoid等静态非线性描述压力感受器转导及ANS分支；用Simulink Transport Delay实现传出时滞；在Simulink/MATLAB中以固定步长ODE3（Bogacki–Shampine）求解，采样时间为1×10^?3 s；以时域轨迹和相图进行极限环数值分析；参数取文献值。本文无样本队列来源。

研究结果。3.1 关键系统变量的时域行为：通过Simulink时域仿真，观察到T_heart与x₆的波动相互耦合，呈现反映压力反射闭环耦合的相干振荡节律；T_heart的总体振荡范围约46 ms，与文献中搏动间期变异性范围相符；以x₆总体范围约38 mmHg计算，得到约1.2 ms/mmHg的范围型压力-周期增益，但该模型估计不能直接等同于临床压力反射敏感指数。长期动力学进入状态空间中的不变集，表现为稳定周期轨道即极限环，提示存在由系统内部反馈结构支配的全局吸引振荡机制，并作为生理调节中的功能平衡。3.2 极限环的数值证据：通过有限维17维状态空间加传出通路传输延迟的数值仿真及相空间投影，观察到不同生理区域初始条件经短暂瞬态后向周期轨道收敛；在（x₆, x_eff1, x_eff2）及（x_eff1, x_eff2）平面的相图中，邻近初始条件的轨迹收敛至同一闭合轨道，符合不变性、周期性和局部收敛性。因此，结果与稳定极限环吸引子一致，提供极限环存在性和局部吸引性的数值证据；但作者明确指出，这属于数值证据，而非严格数学证明存在性、唯一性或全局吸引性。

讨论部分总结。讨论强调应采用homeodynamics（稳态动力学/动态稳态）视角：生理系统并非简单维持固定设定点，而是变量围绕特定轨迹稳定振荡。观察到的约1 s周期性主要反映心血管模型的搏动性质，与既有模型一致。该极限环代表压力反射闭环产生的自持周期节律，即稳态条件下连续心搏的稳定序列。模型优势在于模块化电等效、分室水平描述、显式心室非线性、瓣膜开关阈值非线性和ANS饱和/延迟机制；非线性延迟压力反射可促进振荡，而压力感受器和自主反应饱和（通过f_cs0、f_cs∞、f_p0、f_p∞、f_s0、f_s∞等限制）、瓣膜开闭及非线性心室收缩限制振幅，从而形成稳定闭合轨道。与液压集总模型及简化电模型相比，本模型以控制理论和状态空间形式统一描述CVS-ANS，并在模块层面显式纳入心室动力学、瓣膜非线性和自主反馈通路。讨论也指出局限：血管阻力、顺应性和传输延迟按文献取常数，以保持平稳状态空间表示；未显式建模窦房结电生理，ANS采用静态非线性；呼吸性窦性心律不齐等呼吸影响尚未纳入；当前工作主要为方法学贡献，尚非临床验证工具。

研究结论翻译。本研究提出了一种面向系统理论的途径，用于开发并实际实现CVS-ANS控制机制的闭环17^th阶模型。该模型能够模拟主要系统状态变量的时间轨迹，并突出显示这些轨迹在稳态下表现出与生理条件密切相似的振荡行为，其中变量并非保持恒定，而是在特定范围内自然且稳定地波动（homeodynamics）。所进行的状态空间分析首次为理解稳定极限环如何在CVS-ANS动力学控制系统的综合电等效模型中产生和维持提供了窗口。因此，本研究提供了一种新方法，通过利用信息工程领域方法（电等效和系统理论形式）的高维仿真来研究复杂自主-心血管控制系统中的非线性动力学。一方面，该途径为开发越来越详尽和全面的模型奠定基础，这些模型可通过控制理论原理充分研究；另一方面，它提供了一个实用的in silico平台，用于研究自主调节和心血管稳定性，并可能在情景预测、治疗计划和CVS-ANS病理生理学研究中应用。在这一视角下，未来扩展可纳入血管参数的主动自主调节，从而进一步丰富所提出框架的生理范围。