本节摘要:定义先行——多稳态是多个稳定平衡点在同一参数下并存的现象,最简形态为双稳;它带来三件线性世界没有的事:终态取决于初值(吸引域分界)、外部扰动可触发状态切换、加卸载路径不同则响应不同(迟滞)。本节以杜芬势阱与施密特触发器为例,把这三件事落到图与代码上。
本节是吸引子图鉴的第一站:不动点本尊。线性系统只有一个稳态,"初值决定命运"这句话在它那里没有着落;进入非线性世界,同一个参数下两个稳态并存,初值、扰动、历史全都开始参与决定结局。第二章的分界线概念在这里从画图技巧升级为核心机制。
把稳定平衡想成势能面上的坑最直观。杜芬振子 x'' + 0.4·x' − x + x³ = 0 的势能 V(x) = −x²/2 + x⁴/4 是一条著名的"双井曲线":x=±1 两个谷底(稳定平衡),x=0 一个山脊(鞍点)。停在任何一处的球都会滚进就近的井,但初始位置过了山脊,归属就翻盘——x=0 附近的邻域就是两个吸引域的分界,正是 2.3 节那条由鞍点稳定流形围出的分界线。
对照线性:线性弹簧 F = −kx 的势阱是一只独碗,碗底唯一,任何初值殊途同归。双井比独碗多出的不只是"第二个碗",而是三个新现象:

背景:电子工程里的施密特触发器是双稳最实用的化身——输入电压缓慢变化时,输出只有高、低两个稳态。用无量纲方程 dx/dt = x − x³/1.1 + u 写成"驱动 u 控制下的双井"(x 是输出状态,u 是输入偏置)。
操作:让 u 从 −0.4 缓慢升到 +0.4 再原路返回,每步积分到稳态,记录 x。
import numpy as np from scipy.integrate import solve_ivp def steady(u, x0): sol = solve_ivp(lambda t, s: [s[0] - s[0]**3 / 1.1 + u], [0, 50], [x0], rtol=1e-9) return sol.y[0][-1] up, down = [], [] for u in np.linspace(-0.4, 0.4, 41): # 加载路径:从负稳态出发 up.append(steady(u, -1.5)) for u in np.linspace(0.4, -0.4, 41): # 卸载路径:从正稳态出发 down.append(steady(u, 1.5)) print("u=+0.1 加载输出:", round(up[25], 3), "| 卸载输出:", round(down[15], 3)) print("u=-0.1 加载输出:", round(up[15], 3), "| 卸载输出:", round(down[25], 3))
结果:同一个输入 u=0.1,加载路径输出约 −0.3(仍在"负"支),卸载路径输出约 +1.1(已在"正"支)——输出不由输入唯一决定,而由路径决定。切换发生在 u≈±0.22 附近:那是井底被抬平、稳定解消失的鞍结分岔点(3.2 节的主角)。
解读:迟滞不是缺陷而是功能——两个切换阈值之间的"宽度"就是抗噪能力:噪声幅度小于回线宽度时永远无法误触发。数字电路的噪声容限、恒温器的温度死区、生物膜的开关响应,都是把这段宽度当成护栏来用。变式:把 1.1 调小(井更深)回线变宽、抗噪更强但切换更迟钝;调大则相反——单参数可调的"记忆与灵敏"权衡,正是多稳态器件的设计旋钮。
说"初值决定命运"容易,把命运地图画出来是另一门手艺。数值上的标准做法是网格初值法:在感兴趣的初值区域铺一张均匀网格,每个格点各积分一条轨道到足够长时间,按终点归入哪个吸引子给格点上色,色块的边界即吸引域分界。开销是网格数乘单轨时长,二维系统轻松,四维以上得靠延迟嵌入与抽样技巧。工程上还有一条省钱捷径:先用粗网格定位边界大致走向,再只在边界带内加密——梯度版地图与全网格版在工程精度下几乎无差。实验上的对应物是"扫初条件、等终态",施密特触发器的抗噪宽度正是这样量出来的。两个细节决定地图质量:积分时长必须超过最慢暂态(靠近分界线的点要绕很久才表态),边界附近的着色误差天然偏大——报告分界线位置时应附上这个不确定带,而不是给一条确定无疑的曲线。
迟滞是记忆现象最干净的动力学形态:同一输入对应两个输出,输出取哪个取决于历史路径——系统"记得"自己从哪支解来。但记忆不必依赖迟滞:时延系统(4.4 节)把记忆显式写进方程,慢变量系统用状态编码历史,都无需双稳。反过来,迟滞也不必依赖魔法:它只需要双稳加缓变参数,3.2 节的鞍结分岔负责换轨。分辨要点看机制——"输出与历史有关"是现象描述,"哪个变量、哪段时间、以什么方式参与演化"才是机制答案;画系统框图时把记忆载体标出来,是这类建模最值得养成的习惯。下一章 4.4 节的时延会与这里的迟滞再相遇:一个把记忆放进初值,一个把记忆放进解支选择。
稳态会并存,还会随参数生灭——下一节研究"结构改组"的分类学:分岔。