第2章 · 2.5

2.5 稳定性与扰动

Stability and Disturbance
本节完成第 2 章由静态拓扑结构向动态行为评价的关键跃迁。首先从动力系统不变项公理出发,将平衡态形式化定义为代数变换的不动点;引入位移算子严格形式化外部扰动,建立渐近回返极限判据,区分渐近稳定、发散与中性平衡三种动力学形态;以相空间吸引域与分水岭流形刻画多稳态系统;将稳定性理论接回生理学,阐明坎农稳态的主动耗能本质与本质变量保护机制;证明部分平衡与整体平衡的充要等价性定理,提出子系统的“否决权原理”并指出局部稳定与全局稳定的本质分离;最后以稳态机为原型引出双层超稳定性架构,为后续适应性与调节理论提供数学准备。

2.5.1 不变项公理与代数平衡态

在 2.4 节的末尾,我们留下一道核心悬案:同样是在结构上构成闭环负反馈,瓦特的离心调速器能够平稳消除转速偏差,而高增益的伺服机构或存在神经病变的肢体运动回路却会爆发剧烈发散的振荡。要解开这一动力学分歧,科学分析必须首先建立关于稳定性(stability)的严格定义。正如著名应用数学家 Richard Bellman 曾针砭指出的那样,“稳定性”是现代科学中被滥用最为严重、歧义丛生的词汇之一(Ashby, 1956, S.5/1)。控制论的任务,正是拨开直觉与日常语汇的迷雾,将其完全还原为状态空间与变换算子的精确事实。

稳定性理论的原初公理基础是不变项(invariant):即在系统经历时间连续推移或相继离散变换演化的全过程中,某项关于系统状态的命题始终严格为真(Ashby, 1956, S.5/2)。设想一尊平稳放置在水平桌面上的实心立方体:若将其倾斜 \(5^\circ\) 后由静止释放,命题“其倾角在某一具体时刻恰好为 \(1^\circ\)”时真时假;然而命题“其瞬时倾角在任何时刻均不超过 \(6^\circ\)”则无论立方体如何左右摇摆均恒成立——这一界限构成了系统的一个不变项。与之相对,若将一只尖端朝下倒立的圆锥同样倾斜 \(5^\circ\) 释放,无论我们将角度边界放宽至何种数值,该界限均将在有限时间内被彻底击穿。“在时间演化中无法维持任何有界的空间不变项”,正是“不稳定”在物理上最原始、最直观的含义。

定义 2.7

平衡态(equilibrium state):设系统由封闭单值变换 \(T\) 支配。若系统处于某一特定状态 \(x^*\) 时,其在算子作用下的单步后继仍严格为其自身,即满足代数不动点方程:

\[ T(x^*) = x^* \]

则称状态 \(x^*\) 为系统在变换 \(T\) 作用下的平衡态(亦称不动点)。

对于由状态向量 \(\mathbf{x} = (x_1, \ldots, x_n)\) 刻画的复合系统,平衡态要求所有分量必须同时保持恒定;对于连续时间自主微分系统 \(\dot{\mathbf{x}} = \mathbf{f}(\mathbf{x})\),平衡态充要于所有状态变量的时间导数同时严格为零(\(\mathbf{f}(\mathbf{x}^*) = \mathbf{0}\))。在运动学图上,平衡态表现为自返有向环或无出边的陷阱节点(Ashby, 1956, S.5/3)。

在解析代数表象下,寻找决定性系统的平衡态等价于求解其不动点方程组。考察二维离散系统 \(U\):

\[ x' = 2x - y + 2, \qquad y' = x + y + 3 \]
(2.5-1) 二维向量系统 U 的单步离散规范表示方程;撇号表示经历一步离散演化后的取值。

在平衡态处,必须严格满足 \((x', y') = (x, y)\)。代入方程组 (2.5-1):

\[ \begin{cases} x = 2x - y + 2 \\ y = x + y + 3 \end{cases} \quad \Longrightarrow \quad \begin{cases} x - y = -2 \\ x = -3 \end{cases} \quad \Longrightarrow \quad (x^*, y^*) = (-3, -1) \]
(2.5-2) 第二方程直接解出 \(x^* = -3\);代入第一方程导出 \(y^* = -1\)。点 \((-3, -1)\) 是该系统全局唯一的代数平衡态。

代入核验验证:\(U(-3, -1) = \bigl( 2(-3) - (-1) + 2,\; -3 + (-1) + 3 \bigr) = (-3, -1)\),完全自洽。对于连续微分系统,求解过程完全同理,仅将单步后继替换为各分量的时间导数:

\[ \frac{dx}{dt} = 2x - y^{2} = 0, \qquad \frac{dy}{dt} = xy - \tfrac{1}{2} = 0 \]
(2.5-3) 连续非线性动力系统平衡态的求解:联立方程解得平衡点 \((x^*, y^*) = (\frac{1}{2}, 1)\)。

平衡态在拓扑上的直接扩展是周期循环(cycle):代表点沿由有限个相异状态构成的闭合有向边链周而复始地循环往复(Ashby, 1956, S.5/4)。2.1 节已严格证明:在有限状态集合上,封闭单值动力系统的一切长程轨道归宿,必为且仅能为平衡态或周期循环两者之一。

2.5.2 位移算子与三种动力学结局

平置的立方体、桌面上的理想台球以及倒立于顶尖上的圆锥,在受力分析中均各自存在一个静态平衡态;然而,这三者在物理本质上的深刻鸿沟,只有在系统受到外部干扰后方能彻底暴露(Ashby, 1956, S.5/6)。

为了在纯代数框架下严格描述现实世界中无处不在的冲击与噪声,控制论将扰动(disturbance)形式化为一个瞬时作用于系统状态的位移算子(displacement operator) \(D: S \to S\)。扰动算子本身并不修改支配系统运转的固有变换法则 \(T\),它仅仅是在某一瞬时以“外力”将系统代表点从原平衡态 \(x^*\) 强行位移至新构型 \(D(x^*)\);此后,扰动撤除,系统重新完全遵循原变换 \(T\) 自主演化。

由此,关于系统对特定扰动的稳定性检验,被形式化为一个标准的极限判定问题:将系统置于扰动状态 \(D(x^*)\) 释放,考察其向前迭代轨道能否渐近回返至原平衡态:

\[ \lim_{k \to \infty} T^k\bigl(D(x^*)\bigr) = x^* \]
(2.5-4) 基于位移算子的渐近回返判据;表征系统在经历扰动 \(D\) 后自发消除偏差、收敛回到初始不动点的吸引性。

依据公式 (2.5-4),上述三个力学原型精确对应三种根本相异的动力学结局(Ashby, 1956, S.5/6):

  1. 稳定平衡(Stable Equilibrium):如立方体。位移算子将其推开 \(5^\circ\),算子 \(T\) 在后续演化中逐步耗散多余动能,使轨道单调或振荡地收敛回 \(0^\circ\),严格满足极限式 (2.5-4);
  2. 不稳定平衡(Unstable Equilibrium):如尖端倒立的圆锥。受扰位移后,算子 \(T\) 的相继迭代不仅无法将轨道带回原点,反而使代表点呈指数级迅速背离 \(0^\circ\),极限式 (2.5-4) 彻底失效;
  3. 中性平衡(Neutral Equilibrium):如光滑水平桌面上的理想台球。外力扰动将其推至新位置 \(D(x^*)\) 后,台球既不反向返回出发点,亦不自行进一步加速远离,而是永久停驻于扰动所达到的新坐标:
    \[ \lim_{k \to \infty} T^k\bigl(D(x^*)\bigr) = D(x^*) \neq x^* \]
    (2.5-5) 位移既不被衰减抵消,亦不被发散放大;系统在受扰后的位置处形成新的平衡态。
注记

控制论回返判据与现代李雅普诺夫(Lyapunov)稳定性的严格辨析:Ashby(1956)在此处建立的判据 (2.5-4) 偏重于动力系统的吸引性(attractivity)。在现代严格动力系统理论(Lyapunov, 1892;Tedrake, Underactuated Robotics)中,严格的稳定性分为两层:

  • 李雅普诺夫稳定(Lyapunov Stability):对任意给定的邻域有界误差 \(\varepsilon > 0\),均存在 \(\delta > 0\),使得只要初始扰动满足 \(\|x_0 - x^*\| < \delta\),对于所有未来时刻 \(t \ge 0\) 均有 \(\|x(t) - x^*\| < \varepsilon\)(轨道始终不离开 \(\varepsilon\) 管道);
  • 渐近稳定(Asymptotic Stability):系统同时满足李雅普诺夫稳定,且轨道具有局部吸引性(即满足极限式 (2.5-4))。中性平衡的台球在李雅普诺夫意义下是严格稳定的(取 \(\delta = \varepsilon\) 即可),但它不具备吸引性,因而不是渐近稳定的。控制论教材在使用“稳定”一词时,通常以渐近稳定为标准语义。

这一判据暗含两个具有深刻现实指导意义的方法论前提:

  • 微扰的物理持续性与自然选择:在真实的物理与工程世界中,微观扰动从未彻底停息(如电网热噪声、机械微震动、分子布朗运动)。因此,不稳定平衡态在自然界中根本无法持久存续,稍有扰动便瞬间瓦解;实践中能够被观察者长期稳定观测到的系统形态,必为稳定平衡态(Ashby, 1956, S.5/7);
  • 稳定性对扰动算子集合的相对性:脱离明确指定的扰动算子集合 \(D\),抽象谈论“系统是否稳定”在工程上是残缺的。一根平放于方形底座上的粗铅笔,对于倾斜 \(1^\circ\) 的微弱扰动是渐近稳定的,而对于倾斜 \(15^\circ\) 的大幅度冲击则是不稳定的;断言收音机放大电路“稳定”,默认的扰动集合是指微小的电网纹波,绝不包含雷击强电磁脉冲。在科学陈述中,**稳定性必须相对于特定的扰动集合与特定的相空间边界而确立**(Ashby, 1956, S.5/8)。

2.5.3 吸引流域与多稳态几何

在宏观复杂动力系统中,稳定性绝非仅仅围绕单一孤立平衡态展开。若一个特定的状态子集在变换作用下保持封闭,系统代表点一旦落入其中便永不溢出,该集合构成了相空间中的正向不变集(Ashby 称之为稳定区域stable region,S.5/5)。将所有能够最终收敛汇入某一特定极限吸引子(平衡态或周期极限环)的初始状态全体聚合在一起,便构成了该吸引子的吸引域 / 流域(basin of attraction)(Ashby, 1956, S.5/15)。

图 2.5-1 呈现了一个具有两个吸引流域的九状态离散动力系统运动学图:

流域 I 归宿:平衡态 b(单点自返环) a k n b 流域 II 归宿:周期循环 c → d → e → c c d e i g
图 2.5-1 一个九状态单值系统的离散相空间拓扑图。整个相空间被正交划分为两个互不相交的吸引流域:左侧流域 I 的全部轨道单向收敛至孤立平衡态 b;右侧流域 II 的全部轨道收敛至三态极限环 \(c \to d \to e \to c\)。受扰状态落入哪一个流域,直接决定其最终命运。

在连续时间多稳态系统中,不同吸引流域之间的相空间几何分界面,在数学上被称为分水岭(watershed)或分离流形 / 分界线(separatrix)(图 2.5-2):

状态分量 x 状态分量 y 分水岭(分离流形) 吸引子 A 吸引子 B P Q 扰动算子 D
图 2.5-2 连续二维动力系统的多稳态流线与分水岭流形。相空间包含两个稳定吸引子 A 与 B,二者的吸引盆地由高维鞍点稳定流形所构成的虚线“分水岭”完全分割。受扰后状态若落在分水岭左侧(点 P),轨道平滑收敛回到 A;若被推过分水岭落在右侧(点 Q),系统将不可逆地改投吸引子 B。

相空间几何揭示出控制论关于稳定性的另一项核心哲学反思:**稳定性本身是一项客观的中性动力学描述,绝不内禀道德或功能上的褒贬**。汽油蒸气一旦被电火花引燃,系统瞬间坍缩并永久锁定在“剧烈氧化燃烧”这一极其稳定的热力学平衡态,哪怕在外部风吹微扰下依然自持燃烧——这在消防工程师看来是灾难性的过度稳定(Ashby, 1956, S.5/11)。在临床医学中,Ragnar Granit 曾证明著名的灼性神经痛(causalgia):外周神经外伤导致传入感觉神经纤维与传出交感运动神经纤维发生异常的电接触耦联,在脊髓中枢形成一个自激性正反馈环路。该病理系统表现为典型的双稳态双阱系统:要么处于低激发的健康平稳态,要么落入高频放电的病理性剧痛吸引子,一旦强刺激将其推过分水岭,系统将被锁死在剧痛态中无法自拔。因此,控制论工程师与医学家关注的不仅是系统是否具备稳定性,更在于期望的目标状态究竟位于哪一个吸引流域之中。

习题 2.5-1

(相空间流域分解与平衡态上界定理证明题,源自 Ashby, 1956, 练习 5/3/7 与 S.5/15)

设 \(S\) 为有限状态集合,\(T: S \to S\) 为定义在 \(S\) 上的封闭单值变换。

(1) 严格证明系统全体吸引流域构成了状态空间 \(S\) 的一个数学划分(即互不相交且并集等于 \(S\));

(2) 证明系统吸引流域的总数严格等于相空间中极限循环的总数(将每个孤立平衡态视为长度为 1 的平凡循环);

(3) 系统的孤立代数平衡态的总数量,是否可能严格超过其吸引流域的总数量?请给出严密证明。

参考解答

(1) 证明:依据 2.1 节已证明的抽屉原理,从任意状态 \(x \in S\) 出发的向前轨道,必在有限步内收敛进入且仅进入唯一的极限环或平衡态 \(C_x\)。定义二元等价关系 \(\sim\):若从状态 \(x\) 与 \(y\) 出发的轨道最终汇入同一个极限环(即 \(C_x = C_y\)),则记 \(x \sim y\)。

  • 自反性:\(C_x = C_x \implies x \sim x\);
  • 对称性:若 \(C_x = C_y\),显然 \(C_y = C_x \implies y \sim x\);
  • 传递性:若 \(x \sim y\) 且 \(y \sim z\),则 \(C_x = C_y = C_z \implies x \sim z\)。

因此,等价关系 \(\sim\) 将状态集合 \(S\) 严格划分为互不重叠的等价类集合。依据定义 2.2,每一个等价类恰好是对应极限环的吸引流域 \(B(C)\)。故流域构成 \(S\) 的完全划分:\(\bigcup_k B(C_k) = S\) 且 \(B(C_i) \cap B(C_j) = \varnothing\)(\(\forall i \neq j\))。

(2) 证明:每一个极限循环 \(C_k\) 本身包含的状态,在后续演化中永久停留在该循环内部,因而其自身必属于且仅属于其自身的吸引流域 \(B(C_k)\)。由于不同极限循环在集合上互不相交,且每个循环唯一指定了一个吸引流域,反之每个吸引流域依据定义唯一对应一个吸引核(极限环或平衡态),因此吸引流域的总数与系统中互不相同的极限循环(含单点平衡态)总数严格满足一一双射对应,二者在数量上必然严格相等。

(3) 证明:不可能。设系统包含 \(E\) 个孤立平衡态,包含 \(M\) 个长度 \(p \ge 2\) 的非平凡周期循环。依据不动点定义,每个平衡态在运动学图上为一个长度为 1 的自返自环,因而每一个平衡态均属于一个周期循环。因此极限循环的总数量为 \(N_C = E + M\)。依据第 (2) 问已证结论,吸引流域总数 \(N_B = N_C = E + M\)。由于周期循环数 \(M \ge 0\),显然有:

\[ E = N_B - M \le N_B \]

即代数平衡态的总数至多等于流域总数,绝对不可能超过吸引流域的总数量。以图 2.5-1 为例:系统仅有 1 个平衡态 \(b\)(\(E = 1\)),有 1 个三态循环(\(M = 1\)),流域总数恰为 \(N_B = 1 + 1 = 2\),显然 \(E = 1 \le 2\)。证毕。

2.5.4 稳态:主动的耗能调节与本质变量

在物理机械与生命有机体之间,稳定性呈现出质的差异。平放的立方体维持其平衡态是一种被动的、静态的几何持存:它不需要耗费一焦耳的内部能量即可永久停留在桌面上;然而在生命科学中,Cannon(1932)所阐明的生理稳态(homeostasis)则完全是另一幅图景:**稳态是一种主动的、动态的、持续做功消耗代谢自由能的高阶负反馈控制过程**。

以哺乳动物的核心体温调节为例(图 2.5-3):机体核心温度被死死锚定在约 37 ℃ 这一极度狭窄的生理调定点之上。这一平衡态的维持依赖于下丘脑调控中枢向下并行激发的两条方向完全相反的效应器支路:

  • 当机体产热过多导致体温偏高(误差 \(e = r - y < 0\))时:下丘脑抑制交感神经缩血管中枢,激活汗腺排汗蒸发散热,同时使外周皮肤血管大幅舒张,加速体核热量向外界对流辐射;
  • 当外界严寒导致体温受抑(误差 \(e = r - y > 0\))时:下丘脑迅速兴奋交感神经,使皮肤毛细血管剧烈收缩以减少散热,同时命令躯体运动神经元激发骨骼肌不自主的节律性高频抽搐(寒战产热),并刺激甲状腺素分泌提升基础代谢率。
生理调定点 r 核心温度 37 ℃ + Σ 误差 e 下丘脑中枢 (控制器) 体温偏高(e<0) 汗腺排汗 · 皮肤血管舒张 散热功率大幅增加 体温偏低(e>0) 骨骼肌寒战 · 血管收缩 代谢产热增加 · 限制散热 核心体温 y (核心受控变量) 外部扰动:气温剧变 · 剧烈运动 − 反馈测量回路:外周皮肤与中枢深部温度感受器实时回传 y
图 2.5-3 哺乳动物核心体温双向负反馈调节结构框图。当下丘脑比较器解算出偏差 \(e = r - y\) 后,分别驱动“散热支路”(排汗、扩血管)与“产热保温支路”(寒战、缩血管)双向拮抗做功,协同将体温拉回狭窄调定点。稳态的实质是依靠持续的负熵流与热力学做功对抗熵增。

Ashby 在此完成了生物学生存概念向控制论状态集合代数的严格映射:在系统状态空间的所有可能取值中,观察者可以划定一个狭窄的子集,称为系统保持结构完整与生存所必须维持的物理边界;在此边界内部的状态变量被称为本质变量(essential variables)(Ashby, 1956, S.10/4)。生命的存续,在数学形式上充要于:在外界连续涌入的扰动序列冲击下,机体的负反馈机制能够始终将本质变量的相空间轨迹约束在生存耐受集合内部(图 2.5-3)。

进一步,控制论必须彻底澄清关于反馈极性与稳定性关系的常见误区:**“负反馈必然稳定,正反馈必然发散”这一经验直觉在数学上是不严密的**。系统的最终稳定性完全取决于转移算子幂迭代的全局特征,而非回路单步局部的正负号速记。习题 2.5-2 给出了一个极其优雅的反例:系统内部两变量完全呈现正向相互促进(严格正反馈),然而由于开环增益衰减,整体系统依然渐近收敛至唯一平衡点。

习题 2.5-2

(正反馈系统的渐近稳定性反例,源自 Ashby, 1956, S.5/10)

考察由两个相互作用的连续变量构成的离散线性动力系统:

\[ x' = \tfrac{1}{2}y, \qquad y' = \tfrac{1}{2}x \]

(1) 从初始状态 \((10, 10)\) 与 \((10, 12)\) 出发各推演一步,严格验证两变量之间的交互影响具有正的极性(正反馈);

(2) 运用代数递推与矩阵特征值,证明对相空间中的任意初始状态 \((x_0, y_0) \in \mathbb{R}^2\),当 \(k \to \infty\) 时均有 \(T^k(x_0, y_0) \to (0, 0)\),且原点为全局唯一渐近稳定的代数平衡态;

(3) 结合本题计算结果,深入评述“正反馈系统必然发散失稳”这一流行工程直觉的逻辑漏洞。

参考解答

(1) 单步计算:

  • 从基准点 \((10, 10)\) 出发:\(T(10, 10) = (\frac{1}{2}\times 10, \frac{1}{2}\times 10) = (5, 5)\);
  • 若将初始 \(y\) 从 10 增加至 12(取点 \((10, 12)\)):下一步状态为 \(T(10, 12) = (\frac{1}{2}\times 12, \frac{1}{2}\times 10) = (6, 5)\)。可见,\(y\) 的正向增加导致下一步 \(x'\) 由 5 提升至 6(偏导数 \(\frac{\partial x'}{\partial y} = \frac{1}{2} > 0\));
  • 同理,保持 \(y=10\) 将 \(x\) 从 10 增至 12,下一步 \(y'\) 亦从 5 升至 6(\(\frac{\partial y'}{\partial x} = \frac{1}{2} > 0\))。

两变量之间均互相施加同向正效应,回路增益为正(正反馈)。

(2) 系统的状态转移矩阵为 \(A = \begin{pmatrix} 0 & 1/2 \\ 1/2 & 0 \end{pmatrix}\)。两步演化有:

\[ T^2(x, y) = A^2 \begin{pmatrix} x \\ y \end{pmatrix} = \begin{pmatrix} 1/4 & 0 \\ 0 & 1/4 \end{pmatrix} \begin{pmatrix} x \\ y \end{pmatrix} = \left( \frac{x}{4},\ \frac{y}{4} \right) \]

数学归纳法可知:对于偶数步 \(2k\),\(T^{2k}(x, y) = \frac{1}{4^k}(x, y)\);对于奇数步 \(2k+1\),\(T^{2k+1}(x, y) = \frac{1}{2 \cdot 4^k}(y, x)\)。由于分母中 \(4^k \to \infty\)(当 \(k \to \infty\)),对相空间中任意初态 \((x_0, y_0)\),其向前轨道均严格满足:

\[ \lim_{k \to \infty} T^k(x_0, y_0) = (0, 0) \]

令 \(T(x^*, y^*) = (x^*, y^*) \implies \frac{1}{2}y^* = x^*\) 且 \(\frac{1}{2}x^* = y^* \implies x^* = \frac{1}{4}x^* \implies x^* = 0, y^* = 0\)。因此原点 \((0, 0)\) 是唯一的全局渐近稳定平衡态。

(3) 该例无可辩驳地推翻了“正反馈必发散”的简单直觉:虽然回路的符号极性是正向相互激发的,但由于闭环回路的单循环总传递增益(loop gain)为 \(g = \frac{1}{2} \times \frac{1}{2} = \frac{1}{4} < 1\),能量沿环路循环传播一圈后几何收缩至原先的四分之一。动力系统的最终收敛性严格由系统转移算子谱半径(特征值模长均严格小于 1)决定,局部的“正/负反馈”标签仅代表导数符号,绝不能替代严格的算子谱半径稳定性判据。

2.5.5 部分平衡、整体平衡与否决权原理

当多个子系统相互耦合构成高维复杂网络时,稳定性在“部分”与“整体”之间呈现出极其优美的代数结构。设整机由 \(n\) 个部件组成,全局状态为向量 \(\mathbf{x} = (x_1, \ldots, x_n)\)。整体处于平衡态,意味着整个向量在时间步进下恒定不动,即每一个分量均保持恒定。就此基本事实,Ashby 推导出了一条极具哲学深度的基本定理(Ashby, 1956, S.5/12):

定理 2.3

复合系统平衡态充要等价定理:一个由若干子系统耦合而成的复合系统处于全局平衡态,**当且仅当每一个组成部件在其相邻其余部件所提供的当前输入条件下,均各自处于局部的平衡态**。

定理 2.3 的逻辑证明极为直接:

  • 充分性(局部平衡 \(\implies\) 全局平衡):若每个部件在其接收到的当前输入参数下均已达到平衡态,则在下一时间步,没有任何部件会主动发生状态改变。由于所有部件状态未变,提供给相互之间的输入参数亦未发生任何变动,于是系统在下一步仍停留在原构型,整体处于平衡;
  • 必要性(全局平衡 \(\implies\) 局部平衡):若整机处于全局平衡态,则所有部件的状态向量恒定。这意味着每个部件所接收到的由其他部件决定的外部参数全部保持恒定;在该恒定参数下,该部件的状态从当前时刻到下一时刻保持未变,这恰好满足了带输入的机器在固定参数下的局部平衡态定义。

这一定理引出了控制论中著名的否决权原理(power of veto)(Ashby, 1956, S.5/13):设想系统正沿相空间轨道演化,若部件 \(A\) 在某一瞬时恰好抵达了它在当前外部条件下的平衡点(\(A\) 倾向于停在此处),然而与其耦联的伙伴部件 \(B\) 此刻尚未达到平衡(\(B\) 必须在下一步发生位移);在下一步中,\(B\) 的移动瞬间改变了它向 \(A\) 输送的输入参数,这导致 \(A\) 的原有平衡态位置发生漂移,\(A\) 不得不在新的参数强迫下被迫打破平衡、重新起步。在控制论意义上,**由于部件 \(B\) 未处于平衡态,它对 \(A\) 的停滞状态行使了单方“否决权”**。任何全局平衡态的达成,必须获得网络中所有组成子系统的全票一致赞成;任何一个处于非平衡动态的子构件,都拥有破坏整体平衡的绝对否决权。

更为关键的推论在于稳定性的动力学分离:各子系统在孤立状态下各自具有渐近稳定性,绝不保证耦合之后的整体系统必然稳定!这一事实直接宣告了简单层级还原论的破产(参见习题 2.5-3)。

习题 2.5-3

(子系统局部稳定与整机失稳的分离反例,改编自 Ashby, 1956, S.5/13)

考察两个具有外控输入参数的换能器部件:

\[ P:\ x' = \tfrac{1}{2}x + \beta,\qquad Q:\ y' = \tfrac{1}{2}y + \alpha \]

其中 \(\beta\) 为部件 \(P\) 的输入参数,\(\alpha\) 为部件 \(Q\) 的输入参数。现通过全双工交叉反馈接口进行耦合:设定 \(\alpha = x\),\(\beta = y\)。

(1) 在输入参数冻结恒定的前提下,分别求解部件 \(P\) 与 \(Q\) 的局部平衡态,并分析它们在受到扰动后的孤立稳定性;

(2) 写出耦合后闭环整机的规范方程组,并从初始状态 \((10, 10)\) 与 \((10, -10)\) 出发计算整机前 3 步的数值轨迹;

(3) 求解整机在原点的平衡态,并判定其稳定性;该结果是否与定理 2.3 产生矛盾?结合本例阐明系统工程中“逐层独立验收合格无法确保整机稳定”的深刻机理。

参考解答

(1) 孤立稳定性分析:

  • 冻结输入 \(\beta = \beta_0\) 时,令 \(x^* = \frac{1}{2}x^* + \beta_0 \implies x^* = 2\beta_0\)。受扰 \(x = 2\beta_0 + \delta\) 后,一步演化为 \(x' = \frac{1}{2}(2\beta_0 + \delta) + \beta_0 = 2\beta_0 + \frac{1}{2}\delta\),偏差每步衰减一半(特征值 \(\lambda_P = \frac{1}{2} < 1\))。因此 \(P\) 在任意固定输入下均全局渐近稳定;
  • 同理,冻结输入 \(\alpha = \alpha_0\) 时,部件 \(Q\) 的平衡态为 \(y^* = 2\alpha_0\),偏差每步衰减一半(\(\lambda_Q = \frac{1}{2} < 1\)),\(Q\) 在固定输入下亦全局渐近稳定。

(2) 耦合整机方程与轨迹计算:代入 \(\alpha = x, \beta = y\),导出闭环整机方程:

\[ \begin{cases} x' = \tfrac{1}{2}x + y \\ y' = x + \tfrac{1}{2}y \end{cases} \]

数值轨迹演化:

  • 从初态 \((10, 10)\) 出发:
    • 第 1 步:\(x_1 = 5 + 10 = 15,\ y_1 = 10 + 5 = 15\);
    • 第 2 步:\(x_2 = 7.5 + 15 = 22.5,\ y_2 = 15 + 7.5 = 22.5\);
    • 第 3 步:\(x_3 = 11.25 + 22.5 = 33.75,\ y_3 = 22.5 + 11.25 = 33.75\)。
    状态沿正对角线以每步 \(1.5\) 倍的公比呈指数爆炸级发散失稳;
  • 从初态 \((10, -10)\) 出发:
    • 第 1 步:\(x_1 = 5 - 10 = -5,\ y_1 = 10 - 5 = 5\);
    • 第 2 步:\(x_2 = -2.5 + 5 = 2.5,\ y_2 = -5 + 2.5 = -2.5\);
    • 第 3 步:\(x_3 = 1.25 - 2.5 = -1.25,\ y_3 = 2.5 - 1.25 = 1.25\)。
    状态沿反对角线以每步 \(-0.5\) 的比率迅速收敛至原点。

(3) 整机稳定性与原理解读:

  • 联立方程 \(x = \frac{1}{2}x + y\) 与 \(y = x + \frac{1}{2}y\),唯一解为 \((x^*, y^*) = (0, 0)\),原点为唯一平衡态;
  • 整机状态转移矩阵为 \(A = \begin{pmatrix} 1/2 & 1 \\ 1 & 1/2 \end{pmatrix}\)。其特征多项式为 \(\det(\lambda I - A) = (\lambda - 1/2)^2 - 1 = 0 \implies \lambda - 1/2 = \pm 1\),解得特征值: \[ \lambda_1 = 1.5,\qquad \lambda_2 = -0.5 \] 由于存在模长严格大于 1 的不稳定正实特征值 \(\lambda_1 = 1.5 > 1\),原点在相空间中呈现典型的鞍点(saddle point)失稳拓扑,按李雅普诺夫判据是严格不稳定的;
  • 与定理 2.3 完全自洽:定理 2.3 仅判定平衡态的存在性,而不判定稳定性。在原点处,当 \(y=0\) 时 \(P\) 确平衡于 0,当 \(x=0\) 时 \(Q\) 确平衡于 0,双方在局部均满足静态平衡条件,故原点确实是整机的平衡态;然而,各个部件在孤立条件下的局部渐近稳定性,完全建立在“外部输入被人工冻结为常数”的前提之下。一旦闭环耦合,输入随状态动态位移,原先成立的稳定性前提荡然无存。这揭示出复杂系统工程的至高纪律:**局部子系统的孤立验收合格绝不能替代闭环整机的全局动态校核,系统失稳的灾难往往滋生于部件交汇的接口反馈拓扑之中**。

2.5.6 超稳定性:当一级反馈不足以自持

普通的固定参数负反馈回路(如体温自主调节)能够有效吸收和抵御日常的微弱扰动;然而,一旦外界施加的环境扰动强度超过了分水岭临界点,系统代表点将被强行击出原有的吸引流域(图 2.5-2),固有的负反馈回路将彻底丧失纠偏能力。面对这一生存危机,生命系统与高级自动化机器开辟了高阶演化路径:**当既定回路内部的偏差持续无法消除时,系统必须启动二级调节机制,主动修改主导回路本身的参数结构,重置动力学变换,以搜寻新的稳定模式**。

Ashby(1952, 1956)将这种设计思想实体化,发明了控制论历史上著名的机械心智原型——稳态机(Homeostat):

  • 子系统 A(一级内部动力学单元):由四组相互电磁耦合的双向指针表构成,代表机体的基础内部状态向量。每一指针的平衡位置取决于其余指针传递的电流及其自身的内部电位器参数。当环境发生外来强迫变动时,指针通常偏离中心零位;
  • 子系统 B(二级自适应换能器):由四个多掷步进开关(uniselectors)与继电器构成,共有 \(25^4 = 390,625\) 种离散参数组合(若包含继电器状态共 \(781,250\) 种状态)。步进开关串联接入系统 A 的电磁线圈之中,其旋转物理更迭了 A 的耦合增益参数;
  • 二级闭环连接准则:当且仅当子系统 A 的所有指针维持在中心安全阈值(本质变量安全区)时,流经继电器的电流被切断,子系统 B 处于静止平衡;一旦任意指针触碰危险边缘,继电器通电闭合,驱动步进开关在极高频率下随机跳步旋转。

在稳态机的运行中,否决权原理展现出震撼的力量:只要子系统 A 的指针未曾回归中心安全区,子系统 B 便永远处于非平衡动态之中,持续行使其“否决权”,驱动参数开关持续跳步;在盲目随机的参数搜寻中,一旦命中某个能够使子系统 A 建立全局渐近稳定的参数构型,四根指针自发收敛至中心,B 随即断电静止。整机由此依靠纯机械机制,自主实现了对未知剧烈扰动的自愈式自适应重构。

这种**受本质变量偏差驱动、动态更迭系统底层动力学参数的两级嵌套反馈组织**,被 Ashby 命名为超稳定性(ultrastability)(Ashby, 1952, 1956, 序言;本教程将在 4.3 节与第 5 章系统推导其数学模型)。

至此,我们完成了机制理论的全部正向构造:状态、变换、决定性机器、耦合、反馈、稳定性与吸引域。然而,科学研究在绝大多数情形下必须面对逆向工程的严峻挑战:面对一个完全封盖、内部物理结构不可直接透视的复杂客体,研究者仅能通过有限的外部输入激励与输出读数进行观测——我们究竟如何从离散的输入输出时序数据中,逆向重构出系统的规范表示与状态变换?这正是 2.6 节“黑箱方法”所要确立的认识论大纲。

关键术语

不变项 (invariant)
在系统动力学演化中保持严格为真的性质、代数关系或空间边界;稳定性的公理化核心基石。
平衡态 / 不动点 (equilibrium state / fixed point)
在代数变换算子作用下映射为其自身的相空间状态,\(T(x^*) = x^*\);连续系统导数全为零的拓扑驻留点。
位移算子 (displacement operator)
表示外部瞬时扰动的代数算子 \(D\),其将系统瞬间置换至偏离平衡态的新构型,而不改变系统自身的内在规律 \(T\)。
渐近稳定 (asymptotic stability)
系统在受到有限位移扰动后,其向前时间轨道始终有界并最终收敛回原平衡态的动力学属性。
中性平衡 (neutral equilibrium)
系统在受扰后既不返回原不动点亦不进一步发散,而是以位移后的终点作为新平衡态持存的物理相态。
吸引流域 (basin of attraction)
相空间中所有正向轨道最终汇入同一吸引子(平衡态或极限环)的初始状态全体构成的开集区域;流域构成相空间的完全划分。
分水岭 / 分离流形 (watershed / separatrix)
在多稳态系统中分割不同吸引流域的拓扑临界超曲面;跨越分水岭的微扰将导致系统最终命运的灾变性倒转。
本质变量 (essential variables)
在系统调节与生存模型中被定义表征生命存续或系统完整性的核心状态变量;稳态的本质在于守住本质变量的有界区域。
否决权原理 (power of veto)
复合系统的全局平衡态充要于每个子部件在相邻条件下同时平衡;任一未处于平衡态的子构件对整体平衡具有单方否决破坏力。
超稳定性 (ultrastability)
当主导负反馈回路无法消除本质变量偏差时,由二级机制自动改变系统底层参数与动态结构的自适应双层反馈机制。

参考文献与延伸阅读

  1. Ashby WR. 1956. An Introduction to Cybernetics. London: Chapman & Hall / New York: John Wiley & Sons. [尤其是第5章:Stability, The Invariant, Basin, Power of Veto, The Homeostat].
  2. Ashby WR. 1952. Design for a Brain. London: Chapman & Hall. [超稳定性的奠基性专著].
  3. Cannon WB. 1932. The Wisdom of the Body. New York: W. W. Norton & Company.
  4. Maxwell JC. 1868. On Governors. Proceedings of the Royal Society of London 16:270–283.
  5. Lyapunov AM. 1892. The General Problem of the Stability of Motion. 博士学位论文, 卡尔可夫大学 (英译本: Taylor & Francis, 1992).
  6. Granit R. 1947. Sensory Mechanisms of the Retina. Oxford University Press. [灼性神经痛病理反馈实验].
  7. Tedrake R. 2023. Lyapunov Analysis. In: Underactuated Robotics: Algorithms for Walking, Running, Swimming, Flight, and Manipulation. MIT Course Notes.