约束与权衡16 分钟

相同的恢复速度,能告诉我们多少抗扰动能力?

恢复速度描述稳定态附近的长期衰减,有限扰动的去向还取决于耦合、方向和吸引域。本文构造局部特征值完全相同的双稳系统:耦合改变后,同样半径0.3的一次冲击可以从全部恢复变成约27.4%的方向转态。再用无耦合反例、原作者四稳态模型和纳米机械实验,区分短时放大、最终恢复、全过程约束与持续噪声下的逃逸。

恢复得一样快,为什么一次冲击后的结局不同

有两个系统,平时都停在一种稳定状态。轻轻推一下某个可见变量,偏差都会按同样的速度衰减。于是我们想知道:受到更大的冲击时,它们是否也一样容易恢复?

答案需要限定。轻推所测到的恢复率,确实说明某些局部扰动能够衰减;但它未必告诉我们另一个变量怎样把扰动传过来,更未必告诉我们多远之外存在另一种去向。本文把这些差别放进可以逐项核对的模型,再看真实实验提供了哪一部分证据。

这里的“抗扰动”先指固定参数下,受到一次指定冲击后,最终回到原来的状态。它与“整个恢复过程中都不越过允许范围”不同,也与“持续噪声下能待多久”不同。吸引域稳定性研究处理前一种去向,全过程存留研究则把沿途是否进入不允许区域单独定义。56

四个读数,分别在回答什么

在稳定态附近,把偏差记作向量$\mathbf u$,线性近似为$\dot{\mathbf u}=J\mathbf u$。若所有特征值的实部都小于零,足够小的偏差最终衰减。

读数 问题 还需要什么条件
局部谱与恢复率 很小的偏差长期怎样衰减? 被激发与观测的变量、线性近似范围
短时最大响应 最有利的扰动方向能先放大多少? 坐标尺度、范数、时间与输入方向
吸引域与返回比例 一次有限冲击后最终去哪里? 完整动力学、冲击范围及其概率分布
全过程存留 恢复途中是否越过允许范围? 允许区域、时间及损伤是否改变动力学

即使一个模型最终返回,途中也可能有很大的偏离;如果那种偏离会造成模型没有描述的损伤,单报“最终恢复”仍然不够。5

下面的自建模型均使用固定的无量纲$x,y$坐标,振幅按欧氏范数比较。半径相同表示这组坐标中的位移大小相同,不表示实物输入能量相同。不能把改变变量单位产生的数值变化当作机制改变。

最小过程:一个衰减的变量,还能把另一个变量推走

本库构造以下对照:

$$ \dot x=x-x^3+Ky,\qquad\dot y=-y,\qquad K\geq0. $$

$x$有左、右两种稳定状态;$y$是受到冲击后会自行衰减的附加变量。$K$决定$y$对$x$的作用。这是机制示例,不是纳米梁的标定方程,也不是新的普遍定律。

无论$K$取多少,稳定点都是$(-1,0)$与$(1,0)$,中间$(0,0)$是鞍点。在左态附近令$u=x+1$,有:

$$ \begin{pmatrix}\dot u\\\dot y\end{pmatrix} =J_K\begin{pmatrix}u\\y\end{pmatrix} +\begin{pmatrix}3u^2-u^3\\0\end{pmatrix}, \qquad J_K=\begin{pmatrix}-2&K\\0&-1\end{pmatrix}. $$

三组$K=0,2,8$的局部特征值完全相同,都是$-2,-1$。只在$x$方向施加小冲击时,$y$始终为零,三组的$x$恢复过程甚至完全相同。这种测量没有激发另一个输入通道。

换成$y$方向,虽然$y(t)=y_0e^{-t}$一直衰减,它仍会在消失之前持续推$x$。在线性范围内:

$$ e^{J_Kt}=\begin{pmatrix} e^{-2t}&K(e^{-t}-e^{-2t})\\0&e^{-t} \end{pmatrix}. $$

右上角表示耦合传递。单独激发$y$时,它对$x$的线性贡献在$t=\log2$达到$Ky_0/4$。这不是增长的本征模,而是两个衰减过程组合成的暂时响应。

按所有单位初始方向取最大振幅比:

$$ G_K=\max_{t\geq0}\|e^{J_Kt}\|_2. $$

本库算得$K=0,2$时$G_K=1$,$K=8$时$G_K=2.0801$,峰值时刻约$0.6451$。解析传播矩阵与独立矩阵指数计算的最大元素差约$5.0\times10^{-16}$。

耦合不对称使$J_K$在$K\ne0$时非正规,即$J_KJ_K^{\mathsf T}\ne J_K^{\mathsf T}J_K$。但非正规不保证振幅先增长。瞬时最大增长率由$(J_K+J_K^{\mathsf T})/2$的最大特征值决定;本例只有$K>2\sqrt2$时为正。$K=2$就是非正规、却没有欧氏范数放大的对照。

左:三组相同局部谱的最大线性振幅比。右:半径固定0.3、改变一次冲击方向后的最终去向。两图比较不同问题,左图增益不能直接换成右图概率。
左:三组相同局部谱的最大线性振幅比。右:半径固定0.3、改变一次冲击方向后的最终去向。两图比较不同问题,左图增益不能直接换成右图概率。

有限冲击要走完非线性过程

从左态出发,瞬间把$y$改为$+0.3$,随后停止输入。$K=2$最终回左态;$K=8$最终到右态。保持大小不变、改为$y=-0.3$,两组都回左态。只看冲击大小,遗漏了方向。

过程分成三步:附加变量保存短暂偏差;耦合把偏差传给$x$;当$x$的偏离进入非线性范围,原先的线性恢复近似不再负责决定去向,必须继续解完整方程。

为寻找门槛,令$z=Ky$。任意$K>0$的完整方程都变为:

$$ \dot x=x-x^3+z,\qquad\dot z=-z. $$

鞍点的稳定流形分隔两种去向。沿该边界有:

$$ \frac{\mathrm dx}{\mathrm dz}=-\frac{x-x^3+z}{z}, \qquad x(z)=-\frac z2-\frac{z^3}{32}+O(z^5). $$

本库从鞍点附近积分边界,找到它与$x=-1$相交于$z_*=1.60721816$,所以正$y$方向一次冲击的临界大小为:

$$ y_*(K)=\frac{1.60721816}{K}. $$

$K=2$的门槛约$0.80361$,$K=8$约$0.20090$,同样的$0.3$落在门槛两侧。另用完整方程从初态向前积分并搜索去向分界,24、40、56三个终止时刻的结果均与上述门槛相符,相对差小于$1.2\times10^{-11}$。这检验本例的数值一致性,不是一般证明。

再把冲击写成:

$$ (x_0,y_0)=(-1+0.3\cos\theta,\;0.3\sin\theta), \qquad\theta\sim\mathrm{Uniform}[0,2\pi). $$

720个等角方向的完整积分给出$K=0,2$全部恢复,$K=8$转态比例$27.50\%$。方向加密至1440个、时间步减半并延长观察,得到$27.4306\%$。独立求圆周与边界交点,转态角区间约$33.2283^\circ$至$131.9778^\circ$,连续角度下比例为$27.4304\%$。

这个比例属于固定半径、均匀方向、只冲击一次的任务,不能读成每单位时间的转态率。

把冲击后的状态分布记为$\mu$,原态吸引域记为$\mathcal B$,返回概率是:

$$ P_{\mathrm{return}}(\mu)=\int\mathbf1_{\mathcal B}(\mathbf x)\,\mathrm d\mu(\mathbf x). $$

分布也是任务的一部分。保持平均位移平方为$0.09$,把固定半径、均匀方向改成只冲击$y$、标准差为$0.3$的零均值高斯分布,本库按上面的门槛计算:$K=2$转态概率约$0.3695\%$,$K=8$约$25.1533\%$。这与圆周协议的$0\%$和$27.4304\%$不同。高斯协议有幅度尾部,也改变方向分布;相同均方大小没有匹配完整输入条件。

前面的$z=Ky$确实能把确定性方程中的$K$消去,却不能消去这个比较:原坐标中相同冲击,在新坐标中变成不同的$z$冲击。距离应同时变为$u^2+(z/K)^2$,噪声协方差与概率分布也必须一并变换。如果换变量后又统一用新的标准圆和标准噪声,实际上已经换了研究任务。

还有一个几何关键点:**固定自治系统的吸引域在自由演化下保持不变。**属于左态吸引域的初态,其后状态也属于左态吸引域。因此不能讲成“先在左态吸引域内,随后被自身瞬态放大甩到外面”。$K=8,y_0=0.3$时,冲击后的完整状态从一开始就属于右态吸引域;后来越过$x=0$只是投影读数改变。持续外部输入或噪声作用下,状态才可能穿过固定无输入系统的边界。

不同K的完整吸引域边界及冲击后的自由轨迹。实线对应正y冲击,点线对应负y冲击。右图只画x,不能据一次过零判断完整吸引域被穿越。
不同K的完整吸引域边界及冲击后的自由轨迹。实线对应正y冲击,点线对应负y冲击。右图只画x,不能据一次过零判断完整吸引域被穿越。

更强的替代解释:局部相同,全局门槛本来就不同

上述比较同时改变短时传递与完整吸引域几何,不能把全部差别只归给一个线性增益。完全没有附加变量,也能恢复率相同、转态门槛不同。考虑另一组自建对照:

$$ \dot x=-x\left(1-\frac xa\right)\left(1-\frac{x}{2a}\right),\qquad a>0. $$

稳定点为$0,2a$,分界点为$a$;两个稳定点的局部特征值均为$-1$,与$a$无关。对应势函数为:

$$ V_a(x)=\frac{x^2}{2}-\frac{x^3}{2a}+\frac{x^4}{8a^2}, \qquad V_a(a)-V_a(0)=\frac{a^2}{8}. $$

从$0$受到大小$0.5$的正$x$冲击,$a=0.25$时到另一稳定点,$a=1$时返回。分界距离差4倍、势垒差16倍,局部恢复率仍一样。这里改变了另一稳定点的位置和全局形状,只匹配局部恢复,不能宣称除了某个机制之外所有条件都相同。

在前一个二维模型中,$K=2$没有线性范数放大,却能被$y=0.9$的一次冲击转态。线性放大既不能单独给转态概率,也不是有限转态必须具备的条件。

原始研究在哪一步支持这条机制

2026年Troude与Sornette的研究尝试把保持局部谱的非正规耦合推进到随机逃逸理论。其出版方接受稿与补充材料已取得;补充说明明确把可分离二维模型视作便于隔离机制的最小情形,非可分离、高维推广列作后续工作。12

本库按固定版本作者代码中grad_U_main的漂移方程另写实现,检验其基础模型环节:

$$ \dot x=-\frac12x(x^2-1)+\frac{\kappa}{20}(1-y^2), \qquad\dot y=-\frac12y(y^2-1)+\frac{1}{20\kappa}(1-x^2). $$

作者用它研究带噪声的四稳态模型;本库关闭噪声,新增一次冲击任务。对100组状态、三个$\kappa$,另写漂移与实际执行的原始函数逐项一致。四个$(\pm1,\pm1)$确实都是稳定点:同号两态的谱为$-1\pm0.1$,异号两态为$-1\pm0.1i$,均不随$\kappa$变化。3

从$(1,1)$把$y$一次改为$1.3$,$\kappa=1,10$返回原态,$\kappa=50$到$(-1,1)$。两套积分容差给出相同去向。这把“谱相同,有限冲击去向不同”接回原始漂移模型;不是原论文随机逃逸率图的完整复现。

本篇没有采用“只测一个非正规性数值,就能普遍预测逃逸率”的更强结论。补充推导的有效噪声与势垒表达需要检查大耦合、弱噪声及快变量近似;控制参数、某一点的本征基条件数和完整路径作用量不能随意替换。2

真实装置的入口是2025年双模纳米机械实验:长50微米的氮化硅梁在持续驱动下,第一、第四振动模态各有高、低振幅状态;振动改变张力,张力又改变另一模态的频率。外加窄带噪声促成切换。作者用振幅、相位时间迹观察组合状态,以30秒记录中的驻留比例绘制概率图。4

它支持装置中存在耦合状态与噪声促成的跳态,却没有执行匹配完整局部谱的冲击对照。相近模态线宽也不能直接当作耦合、受驱动稳定态附近的全部恢复谱。作者留下一个未决:没观察到的中间组合,究竟真正不稳定,还是寿命短到当前时间分辨率看不见,需要降低噪声等进一步区分。这个实验提供物理过程与辨识边界,不能证明所有跳态都由非正规放大造成。4

噪声同样强,还要问它作用在哪里

对左态的局部线性近似,给$x$或$y$分别输入独立白噪声,协方差为$2D\operatorname{diag}(q_x,q_y)$。固定$q_x+q_y=1$,即固定这组坐标中的总强度。稳态协方差$C$满足:

$$ J_KC+CJ_K^{\mathsf T}+2D\operatorname{diag}(q_x,q_y)=0. $$

直接求解得:

$$ \operatorname{Var}(x)=\frac{Dq_x}{2}+\frac{K^2Dq_y}{6}, \quad\operatorname{Cov}(x,y)=\frac{KDq_y}{3}, \quad\operatorname{Var}(y)=Dq_y. $$

$D=0.01,K=8$时,同样总强度的噪声只输入$x$,$x$方差为$0.005$;只输入$y$则为$0.10667$,相差约21.33倍。解析表达与独立连续Lyapunov方程求解相符。只记录噪声大小还不够,输入必须投影到能被耦合传递的方向。

这仍是局部线性过程的方差,不能把它代入某个门槛就宣称得到非线性转态率。真实跨态还需噪声时间结构、完整路径、事件定义与观测窗口,分别留给后续课题。

左:局部恢复率相同的两个标量系统,全局势垒不同。右:同样总噪声强度输入不同方向,局部x方差不同。两图分别属于各自模型,不是实物能量与实测转态率。
左:局部恢复率相同的两个标量系统,全局势垒不同。右:同样总噪声强度输入不同方向,局部x方差不同。两图分别属于各自模型,不是实物能量与实测转态率。

怎样把候选机制分开

对照 区分什么 本篇证据的范围
固定全部局部特征值,改变耦合 局部谱是否足够,还是需要传递与全局几何 自建与原作者漂移模型均有不同去向
固定冲击半径,改变方向 扰动大小是否足够描述任务 强耦合时仅部分角区间转态
固定恢复率,改变标量全局形状 非正规放大是否必需 无耦合反例否定其必需性
固定总噪声强度,改变输入方向 噪声资源与投影是否混淆 局部方差改变;跨态率未算
同时记录最终状态与允许范围 最终恢复是否掩盖沿途损伤 两任务需分别定义;未给实物损伤门槛

只观测一个变量,可能没有激发隐藏模式;只记录占据比例,不能恢复全部转移速度;只看短时过零,不能确定最终到另一稳定态。两态系统的稳态占据比只约束两方向速率之比:把两个速率同时乘以同一倍数,驻留比例不变,等待时间却改变。

这篇回答到哪里

较稳的部分是:局部恢复谱、短时响应、有限冲击去向与全过程约束是不同读数;完整状态和输入方向不能省略。两个自建对照与原作者漂移模型给出可回查的具体因果链。

本文的整合判断是:**把恢复率用作有限抗扰动指标,需要额外核对扰动集合、输入通道和吸引域;关注沿途损伤时,还要增加允许区域。**这些信息充分时,局部恢复测量可以参与预测;仅有相同恢复率,一般不足以确定风险排序。

未解决的是:这些对照在具体纳米梁参数上对应多大冲击;部分观测下怎样重建边界;持续、带记忆噪声下的事件时间与路径怎样变化;高维系统中哪些低维分界还能保留。本篇不把数值比例搬成现实系统的通用风险。

来源与依据

  1. S1 原始理论与数值研究

    Phase transitions without instability from non-normal dynamics

    Virgile Troude; Didier Sornette · 2026

    出版方2026-07-25 Article in Press,未编辑接受稿;网页检索列表与发布日期有不齐,以该文章页为准 · 相关全文已核

    出版方接受稿:Box1;Model;Numerical simulations;Results中式10—25

    已取得13页接受稿,实际核对所列段落与数学表达,不声称逐页通读。用途为保持局部谱的耦合与随机逃逸框架;未复现随机速率图,也不采用通用相变/真实领域外推。Box1印刷符号与其Jacobian需分开核,本篇实际模型依据s3四稳态函数。

  2. S2 原始补充推导

    Supplementary Information: Phase transitions without instability from non-normal dynamics

    Virgile Troude; Didier Sornette · 2026

    MOESM2,随2026接受稿 · 相关全文已核

    Supplementary Note2,特别是S2.1可分离性、S2.3快恢复、S2.5双尺度近似;式S21—S50

    已取得23页,实际核相关段落;未将全部推导视作已逐行证明。大耦合、弱噪声与快恢复是不同条件;非可分离/高维推广在原文仍为进一步工作。本篇不依赖其有效温度公式给概率。

  3. S3 作者原始实现

    nonnormal-kramers: simulation and analysis code

    Virgile Troude · 2026

    固定提交3ab6bf7499393775625d010101dd1fbc79fc6fa0,MIT · 相关全文已核

    escape_rate.py的grad_U_main、grad_U_asym、compute_rate_for_k、simulation;measure_escape.py的estimate_escape_rates;README

    取得固定提交全部7个Python脚本与许可;实际阅读README、escape_rate.py相关实现、measure_escape.py全文及first_exp.py。只提取并实际执行grad_U_main定义,与独立重写在300组参数/状态上核对。使用该漂移新增无噪声一次冲击任务;未运行完整随机模拟程序。计数实现是x符号穿越,不自动等于完整吸引域转移;有限记录与速率处理留待T110。

  4. S4 原始纳米机械实验

    Multistability and noise-induced transitions in dispersively coupled nonlinear nanomechanical modes

    David Allemeier; I. I. Kaya; M. S. Hanay; K. L. Ekinci · 2025

    PRA112,023505;arXiv:2506.17026v2(2025-08-29,PDF内日期2025-09-01) · 相关全文已核

    作者版arXiv:2506.17026v2:II装置与TableI;IVB的Fig9/10及中间态未决;IVC概率的30秒分母;附随Supplemental中的SII、SIV噪声输入与读出带宽

    取得含附录和Supplemental的作者版;实际核上述段落、图注、条件和补充噪声校准。未取得原始时间迹,未复现完整实验图。驱动与模态线宽不等同于全部耦合稳定态恢复谱;中间态未见有真正不稳定/寿命过短两种解释。

  5. S5 原始动力系统理论与数值研究

    Survivability of Deterministic Dynamical Systems

    Frank Hellmann; Paul Schultz; Carsten Grabow; Jobst Heitzig; Jürgen Kurths · 2016

    Scientific Reports6,29654(2016-07-13) · 相关全文已核

    Introduction中吸引域稳定性与不良暂态的区别;survivability及初态概率测度μ定义;Methods定义入口

    全文HTML和PDF已取得,实际核定义及相关引论,不声称重算碳循环/神经/电网案例。用于区分最终返回与全过程允许区域、指定初态分布。

  6. S6 原始研究,本文仅摘要级使用

    How basin stability complements the linear-stability paradigm

    Peter J. Menck; Jobst Heitzig; Norbert Marwan; Jürgen Kurths · 2013

    Nature Physics9,89–92(2013) · 摘要已核

    出版方摘要与题录

    付费正文未取得。仅据原始摘要说明局部线性稳定与非局部吸引域量互补;不引用其图表、数值或具体推导。

最后更新 2026-10-08

来源