跳转至

第十一章：血液循环系统读书笔记

作者

本章作者为 James P. Keener（犹他大学数学系）与 James Sneyd（奥克兰大学数学系），合著 Mathematical Physiology, Vol. II: Systems Physiology（Springer, Interdisciplinary Applied Mathematics Vol. 8/II, 第二版, 2009）。本章是第二卷的开篇，承担着把第一卷建立的"细胞-分子-通道"建模工具推广到全身系统生理学的过渡任务：循环系统是所有器官系统（心、肺、肾、脑、内分泌）共用的"物流底座"，其压力-流量-容量关系的数学化是后续各章的前提。本章的写作风格与 Vol I 第一章一脉相承——从最简的 Poiseuille 流出发，逐层叠加顺应性（compliance）、微循环滤过、心泵功能、神经体液调节、胎儿分流及动脉脉搏传播，每一层都给出可解的方程组（含特征值、Hopf 分岔、激波形成的条件），并显式引用临床数据（动脉硬化、肺水肿、动脉导管未闭、动脉瓣闭锁不全等）以验证模型的预测能力。

内容概述

本章系统建立闭合血液循环系统的层次化数学模型：从单根血管的 Stokes 流（Poiseuille 定律）→ 顺应性血管的压力-面积关系 → 微循环滤过的 Starling 力学 → 心室的顺应性模型与 Starling 定律 → 多室（线性 Guyton）循环模型 → 局部自调节（autoregulation）与压力感受器反射（baroreflex）的反馈控制 → 胎儿循环的解剖学分流（ductus arteriosus、foramen ovale）→ 动脉脉搏传播（Windkessel 模型 + 线性波方程 + 非线性激波）。每一节都给出显式方程、典型参数值（如动脉顺应性 \(C_{sa}=0.0018\) L/mmHg、静脉顺应性 \(C_{sv}=0.22\) L/mmHg、心率 \(F=80\) bpm、心输出量 \(Q=5.6\) L/min），并以临床数据或动物实验数据（Guyton & Hall 1996、Sagawa et al. 1978、Berne & Levy 1993）作为模型校准基准。第 11.6.2 节的延迟微分方程 (11.108)–(11.115) 是本章最具数学深度的部分：它把心率-动脉压反馈环建模为 sigmoid 控制 + 反馈时滞 \(\tau\)，通过线性化分析得到特征方程并证明 Hopf 分岔产生约 5.8 秒周期的 Mayer 波振荡。

核心方程与概念

11.1 血流基础 — Poiseuille 定律

不可压缩 Navier-Stokes 方程 在忽略惯性项后化为 Stokes 方程 \(\mu\nabla^2 \mathbf{u} = \nabla P\)。对轴对称定常流，柱坐标下简化为 $\(\frac{1}{r}\frac{d}{dr}\left(r\frac{du}{dr}\right) = \frac{P_x}{\mu}\)$ 积分两次并施以无滑移边界条件 \(u(r_0)=0\)，得到抛物线速度剖面 $\(u(r) = -\frac{P_x}{4\mu}(r_0^2 - r^2) \qquad\qquad(11.6)\)$ 对截面积分得到总流量（Poiseuille 定律） $\(Q = -\frac{\pi P_x}{8\mu} r_0^4 = -\frac{P_x}{8\pi\mu} A_0^2 \qquad\qquad(11.7, 11.8)\)$ 关键结论：流量正比于半径的四次方，这是血管舒缩（vasodilation/vasoconstriction）成为循环系统首要调控手段的物理基础。Table 11.1 给出人体各段血管的直径、总截面积、平均压力、平均流速——动脉压差主要落在小动脉（arterioles）而非毛细血管，原因是小动脉半径虽小但数量相对少。

11.2 顺应性（Compliance）

Laplace 定律 描述血管壁张力 \(T\)、跨壁压 \(P\)、半径 \(r\) 与壁厚 \(M\) 的关系 $\(T = \frac{Pr}{M} \qquad\qquad(11.12)\)$ 线性顺应性模型 假设体积/面积与压力成正比 $\(V = V_0 + CP, \qquad A = A_0 + cP \qquad\qquad(11.13, 11.14)\)$ 静脉顺应性约为动脉的 24 倍（\(C_{sv}/C_{sa} \approx 122\)），使静脉成为血液储库。当跨壁压为负时血管会塌陷，全血（含红细胞）在动脉压低于约 20 mmHg 时停止流动，血浆在 5-10 mmHg 时停止。

将 (11.14) 代入 (11.8) 可得非线性流量-压差关系 $\(RQ = \frac{1}{3\gamma}\left[(1+\gamma P_0)^3 - (1+\gamma P_1)^3\right] = (P_0 - P_1)\left[1 + \gamma(P_0+P_1) + \frac{\gamma^2}{3}(P_0^2 + P_0P_1 + P_1^2)\right] \qquad\qquad(11.17, 11.18)\)$ 其中 \(R = 8\pi\mu L/A_0^2\)、\(\gamma = c/A_0\)。关键观察：高顺应性（\(\gamma\) 大）下驱动同一流量所需压差更小——这解释了为什么静脉压差远小于动脉。

11.3 微循环与滤过

Starling 滤过方程 描述毛细血管壁的液体通量 $\(\frac{dq}{dx} = K_f(-P_c + P_i + \pi_c - \pi_i) \qquad\qquad(11.24)\)$ 其中 \(P_c\) 为毛细血管内压、\(P_i \approx -3\) mmHg 为组织间液静水压、\(\pi_c \approx 28\) mmHg 为血浆胶体渗透压、\(\pi_i \approx 8\) mmHg 为组织间液渗透压、\(K_f\) 为滤过系数。

与连续性方程 \(dP_c/dx = -\rho q\)（\(\rho\) 为毛细血管阻力）联立，得到对称解 $\(q = Q_i \cdot \frac{\cosh\left(\frac{\beta}{2}(L - 2x)\right)}{\cosh(\beta L/2)}, \qquad \beta = \sqrt{K_f \rho} \qquad\qquad(11.27)\)$ 滤过率 \(Q_f = Q_i - q(L/2)\) 与 \(Q_i\) 的比值为 $\(\frac{Q_f}{Q_i} = 1 - \text{sech}\frac{\beta L}{2} \qquad\qquad(11.32)\)$ 关键结论：\(\beta L\) 是滤过的唯一无量纲参数；漏性大（\(K_f\) 大）、半径小（\(\rho \propto r^{-4}\) 大）的血管滤过更剧烈。

11.4 心输出量 — Starling 定律

心室被建模为时变顺应性腔 $\(V = V_0(t) + C(t)P \qquad\qquad(11.34)\)$ 收缩末期（end-systolic）\(V_{ES} = V_{\min} + C_s P_a\)（收缩顺应性 \(C_s\) 是心室的固有属性），舒张末期 \(V_{ED} = V_{\max} + C_d P_v\)。每搏输出量 $\(V_{\text{stroke}} = V_{\max} - V_{\min} + C_d P_v - C_s P_a \qquad\qquad(11.36)\)$ 心输出量 $\(Q = F V_{\text{stroke}} \qquad\qquad(11.37)\)$ 这就是 Starling 定律 的数学表达：心输出量随静脉回流压（\(P_v\)）线性增加，动脉负荷（\(P_a\)）影响较小。

11.5 循环系统整体模型

单室闭合循环（心 + 顺应性血管 + 单一阻力）的三个方程 $\(Q = F(V_h + C_d P_v - C_s P_a), \quad Q = \frac{P_0 - P_1}{R}\frac{1}{3\gamma}[(1+\gamma P_a)^3 - (1+\gamma P_v)^3], \quad V = V_0 + \frac{3V_0}{4} \cdot \frac{(1+\gamma P_a)^4 - (1+\gamma P_v)^4}{(1+\gamma P_a)^3 - (1+\gamma P_v)^3} \qquad(11.40, 11.41, 11.42)\)$ Guyton 线性两室模型 在小顺应性近似下得到闭式解 (11.55)–(11.57) $\(P_a = \frac{(1+FC_dR_s)V_t}{C_v + (1+FC_dR_s)C_a}, \quad P_v = \frac{V_t}{C_v + (1+FC_dR_s)C_a}, \quad Q = \frac{FC_dV_t}{C_v + (1+FC_dR_s)C_a}\)$ 关键定性结论：(1) 心率 ↑ → 动脉压 ↑ 至上限 \(V_t/C_a\)，静脉压 ↓；(2) 心衰时 \(P_a \downarrow\)、\(P_v \uparrow\)，最终 \(P_a = P_v\) 时无流动；(3) 系统阻力 \(R_s\) ↑ → \(Q \downarrow\)、\(P_a \uparrow\)、\(P_v \downarrow\)、血液从静脉系统向动脉系统转移。

三室（系统动脉 + 系统静脉 + 肺循环）线性化模型 给出的饱和心输出量 (11.86) $\(Q_\infty = \frac{2V_e}{R_{sv}(2C_{sa}+C_{sv}) + R_{sa}C_{sa} + C_p R_p}\)$ 灵敏度分析 (Table 11.2) 显示循环系统对 \(C_{sv}\)（静脉顺应性）和 \(R_{sv}\)（静脉阻力）高度敏感（\(\sigma \approx -85\%\)、\(-76\%\)），对 \(C_{sa}\) 和 \(R_{sa}\) 不敏感——这解释了为何体液调节主要通过静脉系统而非动脉系统进行。

11.6 心血管调节

11.6.1 自调节（Autoregulation）

组织自调节 假设小动脉阻力与静脉血氧浓度线性相关 $\(R = R_0(1 + A[O_2]_v) \qquad\qquad(11.91)\)$ 结合 \(Q = (P_a - P_v)/R\) 与代谢平衡 \(([O_2]_a - [O_2]_v)Q = M\)，解得 $\(Q = \frac{1}{1 + A[O_2]_a}\left(MA + \frac{P_a}{R_0}\right) \qquad\qquad(11.92)\)$ 当 \(A = 0\) 时退化为无调节的欧姆定律。典型参数：\(A \approx Q^*/(3M^* - Q^*[O_2]_a^*)\)，\(R_0 = P^*(3M^* - Q^*[O_2]_a^*)/(2Q^*M^*)\)，与 Fig. 11.10/11.11 的实验数据吻合。

11.6.2 压力感受器反射（Baroreflex）— 延迟微分方程

Ottesen 1997 模型把循环系统化简为 lumped 两室 + 压力感受器 sigmoid 控制 + 副交感时滞 \(\tau\) $\(c_a \frac{dP_a}{dt} = -\frac{P_a}{R} + V_s F, \qquad \frac{dF}{dt} = \alpha g(P_a(t-\tau)) - \beta(1 - g(P_a)) \qquad\qquad(11.108, 11.109)\)$ 其中 sigmoid 函数 \(g(P_a) = \mu^7/(\mu^7 + P_a^7)\)（\(g'(P) < 0\)）。

稳态：\(\bar{P} = \mu(\alpha/\beta)^{1/7}\)，\(\bar{F} = \bar{P}/(V_s R)\)。线性化后特征方程 $\(c_a \lambda^2 + \frac{1}{R}\lambda - V_s g'(\bar{P})(\alpha e^{-\lambda\tau} + \beta) = 0 \qquad\qquad(11.115)\)$ 设 \(\lambda = iz\) 并分离实虚部得 $\(c_a z^2 - b^2 = a^2 \cos z\tau, \quad \frac{z}{R} = a^2 \sin z\tau, \qquad a^2 = -\alpha V_s g'(\bar{P}),\ b^2 = -\beta V_s g'(\bar{P})\)$ 两式平方相加消去 \(\tau\)（11.118），用 Ottesen 1997 参数（\(c_a=1.55\) ml/mmHg, \(R=1.05\), \(V_s=67.9\) ml, \(\mu=93\), \(\alpha=0.84\), \(\beta=1.17\)）得 \(z=1.09\)、\(\tau=1.15\) s，对应振荡周期 \(\approx 5.8\) s，振幅约 20 mmHg——与 Mayer 波（约 0.1 Hz，10 s 周期）的实验观测定性一致。

11.7 胎儿循环

胎儿肺未通气，肺阻力 \(R_p\) 极高，动脉导管（ductus arteriosus）和卵圆孔（foramen ovale）形成右-左分流。关键方程组 (11.119)–(11.130) 给出 6 条流路的解。出生后第一口呼吸使 \(R_p\) 骤降，\(P_{pv} > P_{sv}\)，卵圆孔关闭；脐带结扎使 \(R_{sa}+R_{sv}\) 骤升，动脉导管内流反向并逐渐闭合。

病理情况： - 动脉导管未闭 (PDA)：约 1/5500 新生儿，主动脉-肺动脉持续分流； - 房间隔缺损 (ASD)：左→右分流，\(Q_p/Q_s = C_{rd}/C_{ld}\)，肺血流量大于体循环； - 肺动脉环缩术（banding）可使 ASD + PDA 患儿 \(Q_p \approx Q_s\)，但仅在 \(R_d \neq \infty\) 时有效。

11.8 动脉脉搏

11.8.1-11.8.2 Windkessel 模型

质量守恒 \(A_t + (Au)_x = 0\)（11.152），动量守恒 \(\rho(u_t + u u_x) = -P_x\)（11.155）。令 \(\rho \to 0\) 则 \(P_x = 0\)，即 Windkessel 模型：压力只随时间变化。广义顺应性 \(\theta(P) = \int_0^L A_0 c(x)/(A_0 + c(x)P) dx\) 给出 ODE $\(Q(t) = \theta(P) P_t + \frac{P}{R} \qquad\qquad(11.161)\)$ 心脏射血时 \(Q \uparrow \Rightarrow P \uparrow\)（血管膨胀）；心脏停射时 \(Q=0 \Rightarrow P_t = -P/(\theta(P)R)\)（血管回弹驱动外周流动）——动脉系统本质上是"风箱"。

11.8.3 线性波方程

忽略非线性项后 $\(\rho u_t + P_x = 0, \quad c P_t + A_0 u_x = 0 \qquad\qquad(11.162, 11.163)\)$ 交叉微分得经典波方程 \(P_{tt} = (A_0/(c\rho)) P_{xx}\)，波速 \(s = \sqrt{A_0/(c\rho)} \approx 4\) m/s（与颈动脉-踝动脉脉搏波到达时差一致）。D'Alembert 解 \(P(x,t) = f(t - x/s) + g(t + x/s)\)（11.165），前向波 + 反向波。

11.8.4 激波形成

特征线 \(C_\pm\) 满足 \(dx/dt = u \pm K(P)\)，\(K(P) = \sqrt{A(P)/(\rho c)}\)。黎曼不变量 \(u \pm 2K(P)\) 分别沿 \(C_\pm\) 守恒。\(u(0,t) = 2K(P) - 2K(P_0)\) 全场成立，因此 \(C_+\) 斜率 \(s(t) = (3/2)u(0,t) + K(P_0)\) 随 \(u(0,t)\) 单调变化。当 \(u(0,t)\) 先升后降时，\(s(t)\) 也先升后降，特征线收敛形成包络（envelope）——即激波。

激波首次形成位置 (11.183)-(11.184) $\(t_i = \frac{s(t)}{s'(t) + t}, \quad x_i = \frac{s^2(t)}{s'(t)}\)$ \(s'(t) > 0\) 时 \(x_i\) 有限，对应激波形成于距心脏约 40 cm 处（Anliker et al. 1971 数值模拟）。临床意义：主动脉瓣闭锁不全时 \(u(0,t)\) 振幅大、上升陡，导致 pistol-shot 现象（股动脉听诊可闻的"枪击声"）。

关键结论

  1. Poiseuille 定律 \(Q \propto r^4\) 决定了血管舒缩成为循环系统首要调控手段，血管半径改变 16% 即可使流量翻倍。
  2. Starling 定律（\(Q = F(C_d P_v - C_s P_a)\)）是心-血管耦合的线性近似，在生理范围内成立；偏离线性（高 \(P_v\) 或高 \(F\)）时出现饱和。
  3. Guyton 线性两室模型的闭式解显示：(a) 心率 ↑ → \(P_a \uparrow\)、\(P_v \downarrow\)；(b) \(R_s \uparrow\) → \(Q \downarrow\)、\(P_a \uparrow\)、\(P_v \downarrow\)，血液从静脉向动脉转移；(c) 心衰时 \(P_a \downarrow\)、\(P_v \uparrow\)，最终无流动。
  4. 三室模型的灵敏度分析（Table 11.2）显示系统对 \(C_{sv}\)（\(\sigma = -85\%\)）和 \(R_{sv}\)（\(\sigma = -76\%\)）高度敏感，对 \(C_{sa}\)（\(\sigma = -2.8\%\)）和 \(R_{sa}\)（\(\sigma = -1.5\%\)）不敏感——这定量解释了为何体液调节（血容量、血管紧张素）主要作用于静脉系统。
  5. 微循环滤过仅由无量纲参数 \(\beta L = L\sqrt{K_f \rho}\) 决定；漏性高、半径小的血管滤过更剧烈。
  6. 动脉导管未闭约 1/5500，主动脉-肺动脉持续分流导致肺血流量过大，肺阻力上升时可反转为右-左分流。
  7. 压力感受器反射的延迟微分方程模型（Ottesen 1997）在 \(\tau = 1.15\) s 处发生 Hopf 分岔，振荡周期 5.8 s、振幅 20 mmHg，定性重现 Mayer 波（10 s 周期、0.1 Hz 的动脉压自发振荡）。
  8. 动脉脉搏的 Windkessel 模型揭示动脉系统作为"风箱"的功能本质——射血期膨胀储血、舒张期回弹驱动外周流动。
  9. 激波形成的特征线分析给出首次激波位置 \(x_i = s^2/s'\)——陡峭脉搏（\(s'\) 大）在距心脏 40 cm 内即可形成激波，这定量解释了主动脉瓣闭锁不全患者的 pistol-shot 现象。
  10. 胎儿循环的 6 流路模型显示出生时 \(R_p\) 骤降导致 \(P_{pv} > P_{sv}\)，卵圆孔功能性闭合；脐带结扎使 \(R_{sa}+R_{sv}\) 骤升，动脉导管反向并最终关闭。

挑战和开放性问题

  1. Mayer 波与 Traube-Hering 波的统一机制 — 11.6 节末尾明确指出"there's a confusing plethora of oscillations"且"no unified model"。不同作者对 Mayer 波（≈ 0.1 Hz, 10 s）和 Traube-Hering 波（呼吸频率）的分类不一致；延迟微分方程只能解释约 5.8 s 周期但实测约 10 s 周期。Julien (2006) 和 Cohen & Taylor (2005) 提供了多机制并存的统计模型，但反馈延迟是否为核心机制仍有争议。
  2. 微循环的连续介质假设失效 — 11.1 节明确承认 Stokes 方程在毛细血管不适用（红细胞直径接近毛细血管直径），但本章未给出非连续介质修正。Pries & Secomb (2000, 2005) 给出了 Fåhræus-Lindqvist 效应的经验修正，但缺乏从第一性原理推导的微循环流变学模型。
  3. 心室顺应性的真实非线性 — 11.4 节的线性 \(V = V_0 + CP\) 是简化假设，真实心肌的应力-应变关系高度非线性（被动成分的指数硬化 + 主动收缩的钙依赖性激活）。将心脏力学耦合到本章的集总参数模型仍是开放问题——Guyton 1963 之后 60 年仍未有公认的标准模型。
  4. 脉搏波在动脉网络中的反射与分支效应 — 11.8 节只处理了单根均匀直管，忽略了动脉分叉、锥度（tapering）和弯曲。临床观察的"第二波"（常用于评估动脉硬化）需要完整的网络模型，但目前没有解析框架能处理整个体循环的 10^9 根血管。
  5. 胎儿-新生儿过渡期的血流动力学 — 11.7.1 节列出了 PDA、ASD、closed foramen ovale 等病理，但出生第一口呼吸与脐带结扎的时间差（通常 30-60 s）如何影响脑血氧供应仍不清楚。新生儿缺氧缺血性脑病的数学预测模型目前依赖经验阈值。
  6. 血液流变学的非牛顿行为 — 11.1 节明确指出含红细胞血浆的 Stokes 方程描述"is suspect"，但本章未引入 shear-thinning、屈服应力、Fåhræus 效应等非牛顿修正。剪切率依赖的表观黏度在大血管（高剪切率）和小血管（低剪切率）下差异可达 10 倍，对微循环和静脉回流有显著影响。
  7. 动脉导管关闭的分子机制 — 11.7 节提到"For reasons that are not completely understood (probably because of the increased concentration of oxygen)"，氧分压如何触发导管平滑肌收缩的信号通路在数学模型中仍完全缺失——这是胎儿-新生儿循环的"黑盒"。
  8. 自调节的非线性饱和 — 11.6.1 节的线性关系 \(R = R_0(1 + A[O_2]_v)\) 在极低氧时失真，缺氧驱动的最大血管舒张有明确上限但模型未描述。Exercise 9 建议研究 \(R = R_0(1 + a[O_2]_v)/(1 + b[O_2]_v)\) 形式，但未给出临床验证。

个人反思与批判性分析

本章的最大优点在于层次的递进性和方程组的可解性：从 Poiseuille 出发，每引入一个新物理量（顺应性、滤过、心泵、调节、分流、脉搏）都给出一个新方程，并展示如何从解析解或特征值分析中提取临床相关结论。11.5.3 节的三室模型 + 灵敏度表 (Table 11.2) 尤为精彩——它把"循环系统主要通过静脉调控"这一临床共识用 12 个偏导数定量证伪。延迟微分方程 (11.108)–(11.118) 则是生理学 ODE 建模的范例：1 个 sigmoid + 1 个时滞就足以产生 Hopf 分岔和生理相关周期的自发振荡。

然而本章的主要缺陷在于过于依赖线性化：

  • 11.5 节的单室闭合循环明确承认"these are not uniquely determined by this equation, so that this model does not have a unique solution"，转而依赖经验假设"\(\rho \to 0\)"或"\(C_s \to 0\)"补足方程。这种为可解性而牺牲物理性的做法在系统生理建模中需要明示其代价：所有闭式解都是有限参数范围内的近似，超出范围（如极高血压或极低血容量）时数值发散。
  • 11.6.1 节的自调节模型 \(R = R_0(1 + A[O_2]_v)\) 在 \([O_2]_v \to 0\) 时 \(R \to R_0\)（无调节），但实际上缺氧应触发最大血管舒张（\(R\) 应远小于 \(R_0\)）——线性假设的代价是无法描述低氧饱和度的临床真实行为。
  • 11.8.4 节的激波分析假设 \(u\) 足够小使 \(C_-\) 斜率为负，但主动脉瓣闭锁不全时 \(u\) 振幅可达正常 3 倍，此时 \(C_-\) 可能变正，特征线交点可能不存在——模型在"pistol-shot"的关键临床场景下外推失效。

对我自己研究血管生物力学（vascular biomechanics）的启发：(1) 集总参数模型（lumped parameter）的真正价值不在于预测精确数值，而在于揭示参数之间的灵敏度排序——Table 11.2 的 \(-85\%\) vs \(-2.8\%\) 比任何实验都更直观地告诉我们该调控哪个参数；(2) 临床上"动脉硬化"症状常被归因于 \(C_{sa}\) 变化，但本章显示 \(C_{sa}\) 的灵敏度仅 \(-2.8\%\)，真正的元凶是 \(R_{sa}\)（+68\%）和 \(R_{sv}\)（\(-48\%\)）——动脉硬化实际上是外周阻力增加而非顺应性降低，这一区分对药物靶点选择有直接意义；(3) 11.7 节展示的"脐带结扎时序"问题，与我们组研究的"血管手术后血流再分布"在数学上同构——一个外科操作的时间窗决定了组织灌注的瞬态响应。

值得指出，本章的 Starling 滤过模型 (11.24) 是所有微循环教材的标准起点，但它忽略了 glycocalyx（内皮表面层）的关键作用——2000 年后的大量实验显示，有效滤过压主要由 glycocalyx 决定而非传统 Starling 力的代数和。这一修正需要把 \(\pi_c\) 从血浆整体渗透压替换为sub-glycocalyx 渗透压，在数学上引入跨 glycocalyx 的渗透压阶跃（类似半导体 PN 结的能带阶跃），目前文献中尚未见到完整的 ODE/PDE 形式。

重要参考文献

[X1] Guyton, A. C. (1963). Circulatory Physiology: Cardiac Output and Its Regulation. W. B. Saunders. — 本章循环系统整体模型和临床参数的主要参考，11.5.2 的两室线性模型直接基于此书。 [X2] Guyton, A. C., & Hall, J. E. (1996). Textbook of Medical Physiology (9th ed.). W. B. Saunders. — 11.6.1 节自调节的实验数据 (Fig. 11.10, 11.11) 出处；Table 11.1 的血压、流速、截面积数据来源。 [X3] Pedley, T. J. (1980). The Fluid Mechanics of Large Blood Vessels. Cambridge University Press. — 11.1 节血流基础和 11.8 节脉搏分析的参考书；Poiseuille 流的实验验证和动脉脉搏的非线性现象。 [X4] Lighthill, J. (1975). Mathematical Biofluiddynamics. SIAM. — 11.8.3 节波方程解和反射波分析的理论基础。 [X5] Sagawa, K., Maughan, W. L., Suga, H., & Sunagawa, K. (1978). Cardiac Contraction and the Pressure-Volume Relationship. Oxford University Press. — 11.4 节压力-容积环 (Fig. 11.5) 的实验数据出处；收缩末期斜率 \(C_s\) 的实验定义。 [X6] Ottesen, J. T. (1997). Modelling of the baroreflex-feedback mechanism with time-delay. Journal of Mathematical Biology, 36(1), 41-63. — 11.6.2 节压力感受器反射的延迟微分方程模型原型；Table 11.3 的参数值来源；Mayer 波振荡的 Hopf 分岔分析。 [X7] Berne, R. M., & Levy, M. N. (1993). Physiology (3rd ed.). Mosby. — Fig. 11.4 微循环示意、Fig. 11.16 动脉脉搏狗实验数据出处。 [X8] Hoppensteadt, F. C., & Peskin, C. S. (2001). Modeling and Simulation in Medicine and the Life Sciences (2nd ed.). Springer. — 11.5 节多室循环模型（Fig. 11.6-11.9）和 11.7 节胎儿循环模型 (Fig. 11.15) 的原型；Fig. 11.5C 压力-容积环示意来源。 [X9] Pries, A. R., & Secomb, T. W. (2000). Modeling the microcirculation. In Modeling and Simulation in Medicine and the Life Sciences (2nd ed., Springer). — 微循环滤过和血流-阻力关系的现代综述；Fåhræus-Lindqvist 效应的处理方法。 [X10] Anliker, M., Histand, M. B., & Ogden, E. (1971). Dispersion and attenuation of small artificial pressure waves in the canine aorta. Circulation Research, 28(5), 572-591. — 11.8.4 节激波形成的数值模拟实验数据；主动脉瓣闭锁不全时激波在 40 cm 内形成的依据。 [X11] Cohen, M. A., & Taylor, J. A. (2002). Short-term cardiovascular oscillations in man: Measuring and modelling the physiologies. The Journal of Physiology, 542(3), 669-683. — Mayer 波与 Traube-Hering 波的区分与综述。 [X12] Danielsen, M., & Ottesen, J. T. (2004). Exact linearization and boundary control of a model of the baroreceptor-reflex. Journal of Mathematical Biology, 48(2), 175-195. — Ottesen 1997 模型的精确线性化与边界控制；sigmoid 函数 \(g(P_a) = \mu^7/(\mu^7 + P_a^7)\) 的实验依据。