6.2 生物与生命科学：种群、神经元与心律

6.2 生物与生命科学：种群、神经元与心律

本节摘要：定义先行——生命系统里非线性不是"扰动项"而是主干：种群依靠密度依赖的增益与饱和维持存续，神经元依靠可激发性与阈值完成信号编码，心脏依靠起搏细胞同步与螺旋波动力学协调泵血。本节沿"个体数量—单个细胞—器官节律"的尺度阶梯，各解剖一个模型，并回答同一句话：线性版本的这个生命系统，为什么根本活不下去。

本节上承 2.2 的逻辑斯蒂对照与 3.3 的极限环：种群模型是逻辑斯蒂的正宫用法，神经元是可激发介质的原型，心律是同步与螺旋波的主战场。对照轴：线性模型里的生命是常数解，非线性模型里的生命才是过程。

种群：同一个 r，两种政策

2.2 节已证明连续与离散逻辑斯蒂模型"同名不同命"，这里补上管理含义。连续版 dx/dt = r·x·(1 − x/K)：r 是内禀增长率，K 是承载量，任何 r>0 都平滑收敛到 K——最大可持续产量在 K/2 处（该点增长速率最大），这是渔业与林业配额理论的数学底座。离散版 x(n+1) = r·x(n)(1 − x(n/K))：世代不重叠的物种（许多昆虫、一年生植物）r 越过 3 后进入周期与混沌——r≈3.7 的种群年际波动看似随机，其实是确定性混沌。

管理含义的对照清单：连续世代物种，配额定在 K/2 附近且波动可平滑预测；离散世代物种，r 一旦越过临界值，"平均而言稳住"是幻觉——单年崩塌可以在没有任何外因的情况下发生（3.5 节瀑布的生态现场）。同一个 r 字母，两套模型，两种政策：先判世代结构，再谈配额数字。

神经元：阈值与可激发，编码的物理基础

神经元的放电不是对输入的线性放大，而是阈值触发：静息时无论多小的刺激都衰减回静息（稳定平衡点管着），超过阈值则爆发一个固定形状的动作电位，然后回位。这就是"可激发"——3.1 节双稳与 3.3 节极限环的组合拳。经典的 FitzHugh-Nagumo 模型两行抓本质：

• v' = v − v³/3 − w + I（快变量：膜电位）
• w' = (v + a − b·w)/τ（慢变量：恢复通道）

恒定输入 I 在零与阈值之间时，系统静默；I 越过阈值，轨迹冲出去绕一大圈（发放一个尖峰）再回来。线性化的判稳（2.4 节）在静息点有效，但"发放"这个事件本身是全局轨道行为——线性框架里根本不存在"越线即全有"的语言。神经编码的关键事实——刺激强度编码在发放频率里而非幅度里（全或无）——正是非线性动力学的直接推论。

心律：同步救了命，同步也能要命

心脏的节律来自窦房结起搏细胞群的同步放电（4.3 节同步的生理版），电信号作为行波扫过心房心室，驱动协调收缩。病理发生在波的组织方式崩溃时：心房颤动对应多条杂乱的小波；心室颤动则与螺旋波有关——行波在某处断裂后卷成自持旋转的螺旋，绕着不可激发的核心打转，泵血节律彻底瓦解。螺旋波动力学是 3.3 节旋转波在可激发介质里的化身：它的钉扎、漂移与失稳，决定了一块心肌是"抖一下就恢复"还是"颤到底"。

除颤器的工作原理在这个框架里一目了然：不是"电击 restart"这么含糊，而是用一次强场把全部可激发态同时重置、抹平断裂的波前，让行波重新组织——用一次受控的强扰动，消灭失控的自持波。对照线性处置（比如"加点阻尼"式的弱干预）在这里完全无效，因为问题不是振荡太强，而是波的拓扑组织错了。

生命系统三个尺度对照表

展开案例：给 FitzHugh-Nagumo 数一遍发放频率

背景：验证"刺激强度编码在频率里"：恒流输入 I 从阈值下方扫到上方，统计发放频率。

操作：取 a=0.7、b=0.8、τ=12.5，I 取几档，各积分两百个时间单位，数 v 的尖峰数：

import numpy as np
from scipy.integrate import solve_ivp

a, b, tau = 0.7, 0.8, 12.5
def fhn(I):
    def rhs(_t, s):
        v, w = s
        return [v - v**3 / 3 - w + I, (v + a - b * w) / tau]
    sol = solve_ivp(rhs, [0, 200], [0.0, 0.0], rtol=1e-9, max_step=0.05,
                    dense_output=True)
    vs = sol.sol(np.linspace(50, 200, 3000))[0]      # 丢弃暂态
    spikes = np.sum((vs[1:] > 1.0) & (vs[:-1] <= 1.0))
    return spikes / 150                              # 每单位时间发放数

for I in (0.2, 0.34, 0.5, 1.0, 2.0):
    print(f"I={I}: 发放频率 ≈ {fhn(I):.3f}")
​

结果（典型输出）：I=0.2 时频率为零（静默区）；I=0.34 起跳到约 0.5；I=0.5、1.0、2.0 依次约 0.7、0.9、1.2——幅度不涨、频率渐涨，全或无加频率编码同时现形。

解读：对照线性放大器的输出（幅度随输入线性涨、无阈值、无饱和），可激发系统的输入—输出曲线是一张"频率阶梯"。变式：把 I 改成周期输入并扫频率，会看到 3.3 节预告的锁频（1:1、2:1 拖带）——深脑刺激与心脏起搏拖带用的是同一条数学通道。

⚠️ 常见坑：把种群混沌、神经发放、心律失常都归因于"环境扰动"。三者分别由密度依赖、可激发性、波组织失效内生驱动——归因错位会导致管理动作全部打空（比如往混沌种群年度波动里找气候原因，而 r 早已越线）。

本节要点回顾

• 世代结构决定政策：连续收敛对离散分岔，配额逻辑完全不同；先判世代，再谈数字。
• 可激发 = 阈值 + 全或无 + 频率编码：线性放大器给不了这三种性质，而它们是神经编码的地基。
• 心律失常是波的拓扑病：螺旋波失稳是组织问题不是强度问题，除颤的原理是重置波前。
• 同一种非线性，三种生命形态：密度依赖、阈值可激发、同步行波——工具箱一以贯之。

生命之后是社会与市场：那里做不了实验，只剩观测数据——下一节看确定性阶梯的下一级。