科普
编辑
分享
反馈
  • 动力系统在神经科学中的应用
  • 动手实践
  • 练习 1
  • 练习 2
  • 练习 3
  • 接下来学什么

动力系统在神经科学中的应用

SciencePedia玻尔百科
定义

动力系统在神经科学中的应用 是指利用数学模型将单个神经元及神经回路视为随时间演化的动力系统,并依据精确的数学定律研究其膜电位与状态的变化。该领域通过分析相空间几何结构和分岔现象,解释了动作电位、同步脑电波以及决策过程中的胜者全得逻辑。它结合了物理学与工程学中的最优反馈控制和共振等原理,为理解大脑的运动控制和感觉处理功能提供了强有力的理论框架。

关键要点
  • 动力系统提供了一个数学框架来模拟神经过程,从“漏电积分”神经元到复杂的网络振荡。
  • 动作电位、决策等神经现象可以被理解为系统相空间中的分岔和稳定状态等动力学事件。
  • 动力系统原理将神经科学与物理学、工程学和医学等领域深刻地联系起来,揭示了如共振、反馈控制和适应等共通概念。

引言

我们的大脑由数百亿个神经元构成,是宇宙中最复杂的系统之一。它如何产生思想、情感和行为?这一问题长期以来困扰着科学家。要理解这一庞大网络的集体行为,仅仅罗列其生物组成部分是远远不够的。我们缺少一种能够描述其动态演化和涌现特性的通用语言。本文旨在填补这一知识鸿沟,表明动力系统——一门研究事物如何随时间变化的数学分支——正是理解大脑运作逻辑的强大钥匙。

在接下来的篇章中,我们将踏上一段从理论到应用的探索之旅。首先,在“原理与机制”部分,我们将学习如何用简洁的数学方程来描述神经元、离子通道和突触等基本构件的动态行为。接着,在“应用与跨学科连接”部分,我们将看到这些模型如何解释学习、决策和运动控制等复杂功能,并揭示神经科学与工程学等领域的深刻联系。最后,通过一系列动手实践的练习,您将有机会亲自应用这些动力学概念来解决真实的计算神经科学问题。通过这次旅程,您将发现,看似神秘的大脑活动,实际上遵循着优美而深刻的动力学法则。

现在,让我们从最基本的单元开始,深入探索构成神经活动核心的原理与机制。

原理与机制

我们的大脑,这个由数百亿个神经元组成的复杂网络,是如何产生思想、情感和行动的呢?乍一看,这个问题似乎深不可测。但如果我们像物理学家一样,试图寻找其背后普适的原理,我们就会发现,大脑的许多核心功能都可以通过一个优美而强大的数学框架来理解——那就是动力系统​。

让我们抛开所有复杂的生物化学细节,把神经元想象成一个微小的计算设备。它的核心任务是什么?接收输入,进行处理,然后决定是否产生输出。这个过程,本质上就是一个随时间演化的动态过程。

神经元:一个会“漏水”的桶

想象一个最简单的神经元,它只有一个任务:整合来自其他神经元的输入信号。我们可以把它比作一个底部有小孔的桶。流进桶里的水是来自其他神经元的输入电流 Iin(t)I_{in}(t)Iin​(t),而桶里水位的高度就是神经元的膜电位 V(t)V(t)V(t)。因为桶会漏水,所以水位不会无限上涨;水漏出去的速度,正比于当前的水位 VVV。这个简单的比喻,恰好可以用一个简洁的微分方程来描述,即“漏电积分”模型:

τdVdt=−V+RmIin(t)\tau \frac{dV}{dt} = -V + R_m I_{in}(t)τdtdV​=−V+Rm​Iin​(t)

这个方程告诉了我们什么?左边的 dVdt\frac{dV}{dt}dtdV​ 是电压 VVV 随时间变化的速率,也就是我们的“水位”上升或下降的速度。右边的 −V-V−V 项代表着“漏电”——电压越高,泄露得越快,这使得神经元有遗忘过去的趋势。RmIin(t)R_m I_{in}(t)Rm​Iin​(t) 项则代表着输入电流带来的“充电”效应。

现在,让我们做一个小实验。假设这个神经元一开始处于静息状态 (V=0V=0V=0),然后我们突然给它一个恒定的输入电流 I0I_0I0​。会发生什么呢?就像往漏水的桶里持续注水一样,”水位“ V(t)V(t)V(t) 会开始上升,但由于泄漏,它上升的速度会越来越慢,最终会趋向于一个稳定的最大值 V∞=RmI0V_{\infty} = R_m I_0V∞​=Rm​I0​。有趣的是,这个过程不是瞬间完成的。它遵循一条优美的指数曲线。

其中,参数 τ\tauτ 扮演了至关重要的角色,它被称为​膜时间常数。它衡量了神经元对输入的反应有多“迟钝”,或者说,它“记住”过去输入的时间有多长。一个非常有趣的问题是:电压达到其最终稳态值的一半需要多长时间?通过求解这个简单的微分方程,我们得到了一个令人惊讶的简洁答案:t1/2=τln⁡2t_{1/2} = \tau \ln 2t1/2​=τln2。这个时间只取决于 τ\tauτ 本身!这告诉我们,神经元整合信息的基本时间尺度,是由其自身的物理性质决定的。

感知之门:离子通道的随机舞蹈

神经元的“漏电”和“充电”从何而来?如果我们用更强大的显微镜去观察神经元膜,我们会发现成千上万个微小的蛋白质分子,它们像一扇扇门一样镶嵌在细胞膜上,控制着带电离子(如钠离子、钾离子)的进出。这些就是​离子通道。

这些门并不是简单地一直开着或关着。它们的状态在“关闭”(C)和“开放”(O)之间随机切换。我们可以用一个简单的模型来描述这个过程。从关闭到开放有一个速率 α\alphaα,从开放回到关闭则有另一个速率 β\betaβ。对于一大群这样的通道,我们可以问:在任意时刻 ttt,一个通道处于关闭状态的概率 PC(t)P_C(t)PC​(t) 是多少?

这个概率的演化同样遵循一个简单的动力学方程:

dPCdt=−(α+β)PC+β\frac{dP_C}{dt} = -(\alpha + \beta)P_C + \betadtdPC​​=−(α+β)PC​+β

这个方程的形式与我们之前看到的漏电积分模型惊人地相似!它同样描述了一个趋向于平衡的过程。无论初始状态如何,经过足够长的时间,系统都会达到一个稳态,此时关闭和开放的概率由速率 α\alphaα 和 β\betaβ 的比值决定。这揭示了一个深刻的原理:宏观上看似稳定的电学特性(如膜电导),实际上是大量微观离子通道随机“舞蹈”的统计平均结果。

静息态:一个动态的平衡点

到现在为止,我们已经接触了神经元的基本元件。现在,让我们把它们组装起来,构建一个更完整的神经元模型。真实的神经元比一个简单的RC电路要复杂得多,它的状态通常需要不止一个变量来描述。例如,除了膜电位 VVV 之外,可能还有一个或多个“恢复变量” uuu,它们代表了某些离子通道的缓慢激活或失活过程。一个著名的例子是Izhikevich模型,其简化版的动力学方程如下:

dvdt=0.04v2+5v+140−u+Idudt=a(bv−u)\begin{aligned} \frac{dv}{dt} &= 0.04v^2 + 5v + 140 - u + I \\ \frac{du}{dt} &= a(bv - u) \end{aligned}dtdv​dtdu​​=0.04v2+5v+140−u+I=a(bv−u)​

这是一个二维的动力系统。当神经元没有收到强烈刺激时,它会处于什么状态?它会处于“静息态”。在动力系统的语言里,静息态就是一个​不动点(Fixed Point)——在这个点上,所有的变化都停止了,即 dvdt=0\frac{dv}{dt}=0dtdv​=0 且 dudt=0\frac{du}{dt}=0dtdu​=0。

寻找这个不动点,就相当于解一个代数方程组。对于给定的参数,我们可能会找到一个唯一的解,例如 (v∗,u∗)=(−60,−12)(v_*, u_*) = (-60, -12)(v∗​,u∗​)=(−60,−12)。这个点 (v∗,u∗)(v_*, u_*)(v∗​,u∗​) 就是神经元在没有外界干扰时的“家”。所有微小的扰动最终都会衰减,系统会像滚回碗底的弹珠一样,自动返回到这个稳定的不动点。其他更简单的模型,比如一个由电压 vvv 和一个适应性变量 www 构成的线性系统,其不动点也是由系统内部的相互作用和外部输入 III 共同决定的。

神经元之舞:连接与同步

大脑的力量源于连接。单个神经元就像一个乐器,而整个大脑则是一支庞大的交响乐团。神经元之间是如何“交谈”和“协作”的呢?

直接的电学对话

最直接的交流方式是通过​电突触(或称间隙连接),它就像在两个相邻的细胞之间开了一扇窗,允许电流直接流过。想象两个相同的细胞,它们的电压分别是 v1v_1v1​ 和 v2v_2v2​。它们之间的电流正比于电压差 v2−v1v_2 - v_1v2​−v1​。这个系统的动力学可以这样描述:

dv1dt=g(v2−v1)\frac{dv_1}{dt} = g(v_2 - v_1)dtdv1​​=g(v2​−v1​)
dv2dt=g(v1−v2)\frac{dv_2}{dt} = g(v_1 - v_2)dtdv2​​=g(v1​−v2​)

这里 ggg 是耦合强度。如果一开始它们的电压不同,比如 VAV_AVA​ 和 VBV_BVB​,最终会发生什么?通过巧妙的数学变换,我们可以发现,它们的电压差 v1−v2v_1 - v_2v1​−v2​ 会以指数形式衰减到零,而它们的电压之和 v1+v2v_1 + v_2v1​+v2​ 始终保持不变!这意味着,最终两个细胞的电压会达到一个共同的值,这个值恰好是它们初始电压的平均值:VA+VB2\frac{V_A + V_B}{2}2VA​+VB​​。这就像两个连通的装有不同水位的水池,最终水面会达到同一高度。耦合,天然地导致了​同步。

节律的共鸣

神经元更常见的交流方式是通过化学突触,这种相互作用更为复杂。但是,我们可以从一个更抽象的层面来理解同步。如果我们将每个神经元看作一个振荡器(比如一个钟摆),那么它的状态可以用一个相位 θ(t)\theta(t)θ(t) 来描述,表示它在其振荡周期中所处的位置。当两个这样的振荡器相互作用时,会发生什么?著名的Kuramoto模型给出了一个简洁的描述:

θ˙1=ω+Ksin⁡(θ2−θ1)θ˙2=ω+Ksin⁡(θ1−θ2)\begin{aligned} \dot{\theta}_1 &= \omega + K \sin(\theta_2 - \theta_1) \\ \dot{\theta}_2 &= \omega + K \sin(\theta_1 - \theta_2) \end{aligned}θ˙1​θ˙2​​=ω+Ksin(θ2​−θ1​)=ω+Ksin(θ1​−θ2​)​

这里 ω\omegaω 是每个神经元自身的振荡频率,KKK 是耦合强度。我们关心的是它们的相位差 ϕ=θ2−θ1\phi = \theta_2 - \theta_1ϕ=θ2​−θ1​ 如何演化。对 ϕ\phiϕ 求导,我们得到一个关于相位差的独立方程:ϕ˙=−2Ksin⁡(ϕ)\dot{\phi} = -2K \sin(\phi)ϕ˙​=−2Ksin(ϕ)。这个方程的不动点在哪里?在 ϕ=0\phi=0ϕ=0 和 ϕ=π\phi=\piϕ=π。通过稳定性分析,我们发现,当耦合强度 K>0K>0K>0 时,ϕ=0\phi=0ϕ=0 是一个稳定的不动点,而 ϕ=π\phi=\piϕ=π 是不稳定的。这意味着,两个神经元倾向于调整自己的节奏,直到它们的相位差为零,即达到“同相”振荡。这正是大脑中各种节律性活动(如脑电波)的数学基础。

可塑的连接:学习与记忆的本质

神经元之间的连接强度——即突触权重​——并不是固定不变的。它们会根据神经元的活动历史而改变。这被称为突触可塑性​,它被认为是学习和记忆的细胞基础。一个经典的法则是“赫布定律”:一起发放的神经元,连接会更强。

我们可以用一个简单的动力学方程来模拟一个突触权重 w(t)w(t)w(t) 的变化:

dwdt=−wτ+αC\frac{dw}{dt} = -\frac{w}{\tau} + \alpha Cdtdw​=−τw​+αC

这个方程描述了一场持续的“拔河比赛”。−wτ-\frac{w}{\tau}−τw​ 项代表着“遗忘”,即突触连接会随着时间自然衰减。αC\alpha CαC 项则代表着“学习”或“增强”,它源于突触前后神经元的同步活动。这场拔河比赛的最终结果是什么?系统会达到一个平衡点,此时遗忘和增强的速率恰好相等。这个平衡权重 weq=αCτw_{\text{eq}} = \alpha C \tauweq​=αCτ 就是一个被“储存”下来的记忆痕迹。它告诉我们,一个记忆的强度,取决于学习的速率、神经元活动的关联性以及遗忘的速度。

生命的火花:作为动力学事件的动作电位

我们终于来到了神经科学中最具标志性的事件——​动作电位​,即神经元发放的电脉冲。它不是一个简单的开关,而是一个壮观的、非线性的动力学过程。

穿越“相空间”的旅程

为了理解动作电位,我们需要在一个更抽象的空间——​相空间​(Phase Space)——中观察神经元的状态。对于一个由电压 VVV 和恢复变量 www 描述的二维系统,相空间就是一个以 VVV 和 www 为坐标轴的平面。神经元的每一个可能状态都对应这个平面上的一个点,而它的演化则是在这个平面上的一条轨迹。

通过分析相空间中的轨迹,我们可以解释一些非常奇特的现象。例如,“阳极中断兴奋”(anode break excitation):给神经元施加一个抑制性的(负)电流,当这个电流被突然撤销时,神经元反而发放了一个动作电位!这听起来有悖常理。但从相空间的视角看,一切都豁然开朗。在抑制性电流的作用下,系统被“推”到了相空间的一个特殊区域,在这里,电压很低,但恢复变量 www 也被拉到了一个很低的值。当电流被撤销时,系统从这个“蓄势待发”的位置开始演化,它的轨迹会划过一个巨大的环路,从而产生了一个动作电位,然后才返回静息态。这就像把一个弹弓向后拉得更远,松手后它会射得更远一样。

爆发的临界点:分岔

神经元是如何从静息状态突然转变为持续发放动作电位的活动状态的呢?这个转变通常不是平滑渐变的,而是“突变”的。在动力系统的语言中,这种由参数变化引起的系统质变被称为分岔​(Bifurcation)。

我们可以用一个描述振荡幅度 rrr 的方程来研究这个转变。当外部刺激(由参数 μ\muμ 控制)很弱时,唯一的稳定状态是静息态 (r=0r=0r=0)。随着刺激 μ\muμ 慢慢增强,到达一个临界点 μup\mu_{up}μup​ 时,静息态突然变得不稳定,系统会“跳跃”到一个大幅振荡的状态,代表神经元开始持续发放动作电位。有趣的是,如果我们再把刺激慢慢减弱,系统并不会在 μup\mu_{up}μup​ 就立刻停止发放。它会继续保持发放状态,直到刺激减弱到另一个更低的临界点 μdown\mu_{down}μdown​ 时,才会突然“跳回”到静息态。这种依赖于历史路径的现象被称为迟滞​(Hysteresis),它是由一种叫做​亚临界霍普夫分岔​的动力学机制造成的。这种机制使得神经元的反应具有“全或无”和“记忆性”的特征。

终极交响乐:霍奇金-赫胥黎模型

将所有这些概念集于一身的,是伟大的霍奇金-赫胥黎(Hodgkin-Huxley)模型。它是一个四维的动力系统,精确地描述了构成动作电位的各种离子电流的复杂相互作用。尽管方程复杂,但我们可以通过一个绝妙的近似来抓住其本质——快慢动力学分析​。

在这个模型中,膜电位 VVV 和钠离子通道的门控变量 (m,h)(m, h)(m,h) 是“快变量”,而钾离子通道的激活变量 nnn 是“慢变量”。我们可以把这个过程想象成:你在一个可以倾斜和移动的地板上行走。你(快变量)的动作非常迅速,总能很快在当前的地板上找到一个最稳的立足点。而地板本身(慢变量 nnn)则在缓慢地倾斜。

动作电位的产生过程就像一出戏剧:

  1. 静息​:你站在地板的“低处”稳定的一点上。
  2. 激发​:一个外部刺激(输入电流)像一只脚把你猛地踹到了地板的“高处”。你迅速在“高处”找到了一个新的立足点。
  3. 去极化与复极化:当你站在高处时,你脚下的地板(慢变量 nnn)开始因为你的存在(高电压)而慢慢向上倾斜。你为了保持稳定,不得不随着地板的倾斜而移动。
  4. 坠落​:你走到了高处平台的边缘——这个边缘在数学上是一个​鞍结分岔点​。在这里,你所站立的稳定立足点与一个不稳定的点碰撞并一起消失了!你脚下一空,瞬间“坠落”回地板的低处。
  5. 超极化​:在坠落后,你可能会落到比初始静息点更低的位置,然后慢慢滑回原来的静息点。

这趟从低到高,再从高处边缘坠落的旅程,其在电压 VVV 上的投影,就是我们观测到的那个标志性的、尖锐的动作电位波形。这个看似神秘的生物电信号,原来是一场在多维相空间中上演的、由分岔和快慢动力学主导的壮丽之舞。

从一个漏水的桶,到一场复杂的相空间之舞,动力系统为我们提供了一把钥匙,让我们能够打开神经科学的大门,窥见其背后普适、优美而深刻的数学原理。我们看到,看似复杂的生命现象,都可以被分解为不同时空尺度上相互作用的动态过程,它们共同谱写了思想与意识的交响乐。

应用与跨学科连接

到现在为止,我们已经领略了描述事物变化的数学语言——动力系统的基本原理和机制。你可能会觉得这很有趣,但这些抽象的方程和相图与我们头颅中那个重达三磅、构成我们思想、情感和行动源泉的复杂器官——大脑——又有什么关系呢?

答案是:关系重大。事实上,动力系统正是我们理解大脑内在逻辑和统一之美的最强大的透镜之一。它让我们能够从纷繁复杂的生物细节中抽丝剥茧,抓住神经活动的本质。现在,让我们走出理论的殿堂,踏上一段激动人心的旅程,去看看这些思想如何在神经科学的广阔天地中开花结果,并与其他学科交相辉映。这不仅仅是理论的应用,更是一次发现之旅,我们将看到,从单个神经元的电脉冲到我们做出决策和控制身体运动,背后都贯穿着同样的动力学原理。

单个神经元的动态生命

让我们从最基本的构成单位——单个神经元——开始。一个神经元并非一个简单的开关,它是一个充满活力的、复杂的动力学实体。

积分器、谐振器与记忆的萌芽

神经元最基本的功能之一是整合它从成千上万个其他神经元那里接收到的输入信号。我们可以用一个非常优雅的“漏桶”模型来描绘这个过程。想象一下,输入信号是流入桶里的水流,而神经元膜的自然漏电特性就像是桶壁上的一个小洞,水会从中慢慢渗出。神经元的膜电位(它的“激活水平”)就是桶里的水位。这个过程可以用一个简单的一阶线性微分方程来描述,即“泄漏积分器”(leaky integrator)模型。

当一个短暂的输入脉冲(比如一个想法的火花)传来时,桶里的水位会迅速上升,然后,即使输入停止,水位也会缓缓下降,而不会立即消失。这短暂的活性“余晖”便构成了短期记忆的最基本形式。这个简单的模型告诉我们,记忆并非存储在某个静态的“文件”里,而是动态地“活”在神经元的持续活动中。

更有趣的是,我们不仅可以用这个模型理解大脑的自然活动,还能用它来指导我们与大脑互动。在一种名为“光遗传学”的前沿技术中,科学家可以利用光来精确控制神经元的活动。如果我们用一束节律性闪烁的光来照射神经元,会发生什么呢?我们的泄漏积分器模型可以精确地预测神经元膜电位的反应,它会随着光的节奏起舞,但其振幅和相位会依赖于光的频率和神经元自身的时间常数。这不仅是一项强大的实验工具,也再次印证了神经元对动态输入的处理能力。

然而,有些神经元远比简单的“漏桶”更为精妙。它们是天生的“谐振器”。就像我们在恰当的时刻推一把秋千,能让它越荡越高一样,这些神经元对特定频率的输入信号反应最为强烈。我们可以用一个受驱阻尼振子的二阶微分方程来描述这种行为。这种谐振特性是产生大脑中各种节律(如theta波和gamma波)的关键,这些节律被认为在注意力、感知和记忆中扮演着至关重要的角色。

适应:习惯的力量与内在的恒温器

大脑不仅仅被动地响应世界,它还在不断地适应世界。你是否注意到,刚走进咖啡店时浓郁的香气,几分钟后就几乎闻不到了?这就是感觉适应。我们可以用一个简单的双变量动力系统来捕捉这个现象的精髓。在这个模型中,神经元对刺激的响应(RRR)会受到一个内部“适应”变量(AAA)的负反馈抑制。当刺激持续存在时,适应变量会缓慢积累,从而逐渐减弱神经元的响应。这个简单的反馈回路,根据参数的不同,既可以平滑地衰减至适应状态,也可以在适应过程中产生优美的阻尼振荡。

除了适应外部世界,神经元还必须维持自身的内部稳定,这被称为“稳态可塑性”。一个神经元不能过于“兴奋”,也不能过于“沉默”,它需要将自己的平均放电率维持在一个健康且信息丰富的“设定点”(r0r_0r0​)附近。这就像一个内在的恒温器。一个简洁的动力学模型可以描述神经元如何通过动态调整其放电阈值(θ\thetaθ)来实现这一目标。如果神经元放电太频繁,它就提高门槛;如果太安静,就降低门槛。这是一个美妙的负反馈机制,确保了神经系统在面对持续变化时仍能保持稳定和功能。

噪声的奇迹:随机共振

我们通常认为噪声是信号的敌人,是需要被滤除的杂质。但在神经系统的非线性世界里,噪声有时却能扮演一个出人意料的建设性角色。这种现象被称为“随机共振”。

想象一个粒子在对称的双势阱中运动,两个阱底代表神经元的静息态和预备放电态,而“放电”则对应于粒子从一个阱跳到另一个阱。如果一个非常微弱的、本身不足以让粒子越过势垒的周期性信号作用于系统,粒子将无动于衷。但现在,如果我们加入一些随机噪声,就像轻轻地晃动整个系统,情况就大不相同了。恰到好处的噪声能够时不时地“助推”粒子一把,使其更容易在周期性信号的“鼓励”下越过势垒。结果是,原本看不见的微弱信号,在噪声的帮助下,竟能被系统清晰地“探测”出来。这个惊人的原理表明,大脑中无处不在的“噪声”可能并非缺陷,而是一种巧妙的设计,用以增强其对微弱感觉信号的探测能力。

神经网络的交响乐

单个神经元固然神奇,但大脑真正的力量源于它们组成的庞大网络。动力系统为我们提供了一把钥匙,去开启理解这些神经元如何“合唱”的奥秘。

动态突触:会思考的连接

神经元之间的连接点——突触——并非一成不变的电线。它们是动态的、可塑的,其连接强度会随着近期的活动历史而改变。这种“短期可塑性”主要有两种形式:

  • 突触抑制(Synaptic Depression):当一个突触被高频率使用时,它可能会“疲劳”,其传递信号的效率会暂时下降。这可以用一个简单的动力学模型来描述,其中可供释放的神经递质资源(xxx)会因高放电率(RRR)而被消耗,其消耗速度超过了补充速度。这就像一支笔的墨水会越写越少。

  • 突触易化(Synaptic Facilitation):与之相反,有些突触在被反复激活后会“热身”,其效率会暂时增强。这通常与突触前末梢的“残余钙离子”有关。

这两个看似简单的规则,为神经计算增添了极其丰富的动态。它们使得神经元不仅能响应“现在”的输入,还能响应输入的“变化”和“历史”。更令人兴奋的是,这些数学模型并非空中楼阁。例如,“残余钙离子假说”可以通过实验来验证。当科学家向突触中注入一种能快速吸收钙离子的化学物质BAPTA时,他们观察到突触易化现象被大大削弱甚至消失了。这完美地展示了理论模型与实验验证之间的强大协同作用。

电路基序:思想的构件

在复杂的神经网络中,科学家们发现了一些反复出现的、小型的连接模式,被称为“电路基序”(circuit motifs)。它们就像是思想的“构件”或音乐中的“乐句”。

  • 竞争与决策​:我们如何在两个选项之间做出选择?大脑中一个普遍存在的机制是“侧向抑制”或“赢者通吃”网络。我们可以用一个简单的模型来描述两个相互抑制的神经元群体。当外部输入同时驱动这两个群体时,它们之间的相互抑制会展开一场“竞争”。最终,活动稍微占优的那个群体会进一步抑制另一个,最终“赢得”竞争,代表着决策的做出。这个简单的动力学原理是构成更复杂决策过程的核心。

  • 时机与精度:前馈抑制​:在神经系统中,时机就是一切。一个被称为“前馈抑制”(Feedforward Inhibition, FFI)的电路基序在精确控制时间上起着关键作用。在这个结构中,一个兴奋性信号不仅会激活下游目标,还会通过一个中间的抑制性神经元,向同一目标发送一个略微延迟的抑制信号。这个兴奋信号和紧随其后的抑制“阴影”,共同创造了一个非常短暂的“机会窗口”,只有在这个窗口内到达的其他信号才能被有效处理。这就像一位敏锐的编辑,为神经信号的传递设定了严格的截止时间,从而大大提高了信息处理的时间精度。

跨越边界:连接行为、疾病与工程

动力系统的视角不仅限于解释微观的神经元和电路,它还能连接到宏观的行为、疾病的根源,甚至与工程学等领域产生深刻的共鸣。

脑环路与节律:从运动控制到疾病模型

大脑中的神经元组成了宏大的反馈环路,例如连接大脑皮层和基底节的环路,后者在运动控制和行为选择中至关重要。一个简化的动力学模型可以揭示,这样的环路结构天然地倾向于产生振荡。这些振荡在正常情况下是协调运动所必需的,但当它们变得异常强大或无序时,就可能导致病理状态,例如帕金森病患者的震颤。有趣的是,模型中的一个常数项(DDD),可以被看作是多巴胺等神经调质的抽象代表。改变这个参数就能极大地影响整个环路的动态,这为我们理解神经调质如何调节大脑状态以及相关药物为何有效提供了深刻的洞见。

神经药理学:设计更智能的药物

药物如何影响大脑?动力系统可以帮助我们定量地理解和预测。以一类被称为“使用依赖性”的通道阻断剂为例,这类药物的特殊之处在于,它们更倾向于阻断那些正在被频繁使用的神经通道。这意味着它们能更精准地靶向那些过度活跃的神经元(例如在癫痫发作中),而对正常神经元的影响较小。通过建立一个描述药物结合与解离过程的动力学模型,药理学家可以预测在不同放电频率下药物的效果,从而指导设计出更有效、副作用更小的“智能”药物。

大脑,一位卓越的工程师:最优反馈控制

我们旅程的最后一站,将神经科学与一个看似遥远的领域——工程控制论——联系起来。当你伸手去拿一个杯子时,你的手臂是如何运动得如此平滑、精准和高效的?一个激动人心的理论——“最优反馈控制”(Optimal Feedback Control)——认为,大脑的运动系统在某种意义上是在解决一个工程优化问题。

这个理论假设,大脑在选择运动指令(utu_tut​)时,会寻求一个最佳平衡:一方面要最小化运动误差(即状态xtx_txt​与目标的偏差),另一方面也要最小化能量消耗或“努力程度”(即控制指令utu_tut​的成本)。这是一个经典的控制论问题,而其解的形式却惊人地简单:最优的控制策略是一个线性的反馈法则,ut=−Kxtu_t = -K x_tut​=−Kxt​。这意味着,运动指令的大小与当前状态的偏差成正比。更妙的是,这样一个看似复杂的计算,在生物学上可能仅仅通过一个具有特定突触权重(KKK)的神经连接就能实现。这个理论不仅为我们理解运动的优雅与高效提供了一个全新的数学框架,也预示着大脑可能在用一套深刻的优化算法来与世界互动,而这正是现代机器人学和人工智能所追求的目标。

结语

回顾我们的旅程,我们从单个神经元的电火花,走到了由亿万神经元组成的、能够决策、行动和适应的宏伟网络,最终甚至触及了大脑作为一名“最优工程师”的深刻思想。贯穿始终的,是动力系统这门描述变化的语言。它如同一条金线,将物理学、数学、生物学、医学和工程学等不同学科的智慧编织在一起。

动力系统不仅为我们提供了分析大脑的工具,更重要的是,它揭示了一种看待大脑的新视角:不再是一个固定的、由无数零件组成的复杂机器,而是一个流动的、充满节律和模式的动态整体。在这背后,隐藏着一种深刻的数学之美与内在逻辑。我们的探索才刚刚开始,而这门美丽的语言,无疑将继续引领我们去发现大脑那“附魔的织机”上更多未知的奇迹。

动手实践

练习 1

大脑的连接并非一成不变,而是会随着经验而改变,这一过程被称为突触可塑性。本练习将探讨一种最基本的可塑性形式——突触抑制,即连接强度因持续的突触前活动而减弱。通过求解一个简单的一阶微分方程,您将亲身体验如何为基本的学习规则建模,并理解在突触层面记忆痕迹是如何随时间衰退的。

问题​: 在计算神经科学中,一个简单的突触衰减模型描述了两个神经元之间的连接强度如何因持续的突触前活动而随时间减弱。考虑一个单一突触,其强度由一个无量纲权重 w(t)w(t)w(t) 表示。该权重的变化率由以下微分方程决定:

dwdt=−ηrprew\frac{dw}{dt} = -\eta r_{pre} wdtdw​=−ηrpre​w

在此,rprer_{pre}rpre​ 是突触前神经元的发放速率,单位为赫兹 (Hz),而 η\etaη 是一个无量纲的学习率常数,决定了突触减弱的速度。

假设一个感觉神经元受到一个恒定不变的背景刺激,使其以 rpre=5.0r_{pre} = 5.0rpre​=5.0 Hz 的恒定速率发放。连接该神经元与其下游神经元的突触初始权重为 w(0)=w0=0.9w(0) = w_0 = 0.9w(0)=w0​=0.9。该突触的学习率为 η=0.08\eta = 0.08η=0.08。

计算在 t=3.0t = 3.0t=3.0 秒后的突触权重 w(t)w(t)w(t)。将最终答案四舍五入至三位有效数字。

显示求解过程
练习 2

在理解了突触之后,我们来看看神经元本身是如何整合数千个输入信号的。本练习将使用计算神经科学的基石模型——“漏电积分-发放” (LIF) 模型来探索突触整合过程。该模型描述了神经元的膜电位 V(t)V(t)V(t) 如何在累加兴奋性和抑制性输入的同时不断“泄漏”电荷,通过对模拟为狄拉克-德尔塔函数 δ(t)\delta(t)δ(t) 的离散突触事件的响应进行分析,您将深入了解单个神经元执行的基本计算。

问题​: 一位计算神经科学家正在研究一种名为“漏电积分-发放”(leaky integrate-and-fire)的神经元简化模型。在此模型中,神经元的膜电位 V(t)V(t)V(t) 随时间的变化遵循以下微分方程:

dVdt=−Vτ+Isyn(t)\frac{dV}{dt} = -\frac{V}{\tau} + I_{syn}(t)dtdV​=−τV​+Isyn​(t)

其中,VVV 的单位是毫伏 (mV),时间 ttt 的单位是毫秒 (ms)。项 −Vτ-\frac{V}{\tau}−τV​ 代表“漏电”,即电位以特征膜时间常数 τ\tauτ 自然衰减回 0 mV 的静息态。项 Isyn(t)I_{syn}(t)Isyn​(t) 代表来自输入突触信号的电流。

假设神经元从其静息电位开始,因此 V(0)=0V(0) = 0V(0)=0 mV。然后,它在不同时间接收到两个瞬时突触输入。第一个是时间 t1t_1t1​ 时的兴奋性输入,第二个是时间 t2t_2t2​ 时的抑制性输入。这些输入可以使用狄拉克 (Dirac) delta 函数 δ(t−t0)\delta(t-t_0)δ(t−t0​) 来建模,该函数表示在时间 t0t_0t0​ 的一个无限尖锐的脉冲。总突触输入由下式给出:

Isyn(t)=w1δ(t−t1)+w2δ(t−t2)I_{syn}(t) = w_1 \delta(t-t_1) + w_2 \delta(t-t_2)Isyn​(t)=w1​δ(t−t1​)+w2​δ(t−t2​)

该特定神经元的参数如下:

  • 膜时间常数:τ=10.0\tau = 10.0τ=10.0 ms
  • 第一次输入的时间:t1=5.0t_1 = 5.0t1​=5.0 ms
  • 第一次输入的权重:w1=15.0w_1 = 15.0w1​=15.0 mV
  • 第二次输入的时间:t2=12.0t_2 = 12.0t2​=12.0 ms
  • 第二次输入的权重:w2=−10.0w_2 = -10.0w2​=−10.0 mV

计算在时间 t=20.0t = 20.0t=20.0 ms 时神经元的膜电位 VVV。最终答案以 mV 为单位表示,并四舍五入到三位有效数字。

显示求解过程
练习 3

单个神经元和突触是基本构建模块,而像记忆这样的复杂认知功能则源于大量神经元群体的集体活动。本练习将研究一个循环连接的神经群体速率模型,这是工作记忆理论的核心。通过分析该系统的非线性动力学,您将学习如何运用分岔分析来确定系统产生多种稳定状态的临界条件,从而将动力系统的抽象概念与大脑的记忆功能直接联系起来。

问题​: 一个描述神经元群体平均活动动力学的简化模型,对工作记忆理论至关重要,由一个一维速率方程给出。该群体的平均活动水平,由无量纲变量 A(t)A(t)A(t) 表示,随时间演化,遵循以下微分方程: dAdt=−A+tanh⁡(wA+I)\frac{dA}{dt} = -A + \tanh(wA + I)dtdA​=−A+tanh(wA+I) 在此模型中,AAA 代表活动,其取值范围为-1到1。参数 www 是一个正常数,表示群体内循环连接的强度;III 是一个恒定的外部输入。对于特定的参数值,该系统可以表现出多个稳定稳态,这一特性被认为是记忆保留的基础。稳态是使 dAdt=0\frac{dA}{dt} = 0dtdA​=0 的 AAA 值。

给定循环连接强度为 w=3w=3w=3,确定系统首次支持超过一个可能的稳态活动水平时的最小正外部输入值 IcI_cIc​。

从以下选项中选择 IcI_cIc​ 的正确表达式。

A. 6−ln⁡(3+2)\sqrt{6} - \ln(\sqrt{3}+\sqrt{2})6​−ln(3​+2​)

B. ln⁡(3+2)−6\ln(\sqrt{3}+\sqrt{2}) - \sqrt{6}ln(3​+2​)−6​

C. 3−ln⁡(3+12)\sqrt{3} - \ln\left(\frac{\sqrt{3}+1}{\sqrt{2}}\right)3​−ln(2​3​+1​)

D. 6−ln⁡(3)\sqrt{6}-\ln(\sqrt{3})6​−ln(3​)

E. 3−ln⁡(2)\sqrt{3}-\ln(\sqrt{2})3​−ln(2​)

显示求解过程
接下来学什么
动力系统
尚未开始,立即阅读
在电路中的应用
模型神经元的放电率