恢复得一样快,为什么一次冲击后的结局不同
有两个系统,平时都停在一种稳定状态。轻轻推一下某个可见变量,偏差都会按同样的速度衰减。于是我们想知道:受到更大的冲击时,它们是否也一样容易恢复?
答案需要限定。轻推所测到的恢复率,确实说明某些局部扰动能够衰减;但它未必告诉我们另一个变量怎样把扰动传过来,更未必告诉我们多远之外存在另一种去向。本文把这些差别放进可以逐项核对的模型,再看真实实验提供了哪一部分证据。
这里的“抗扰动”先指固定参数下,受到一次指定冲击后,最终回到原来的状态。它与“整个恢复过程中都不越过允许范围”不同,也与“持续噪声下能待多久”不同。吸引域稳定性研究处理前一种去向,全过程存留研究则把沿途是否进入不允许区域单独定义。56
四个读数,分别在回答什么
在稳定态附近,把偏差记作向量$\mathbf u$,线性近似为$\dot{\mathbf u}=J\mathbf u$。若所有特征值的实部都小于零,足够小的偏差最终衰减。
| 读数 | 问题 | 还需要什么条件 |
|---|---|---|
| 局部谱与恢复率 | 很小的偏差长期怎样衰减? | 被激发与观测的变量、线性近似范围 |
| 短时最大响应 | 最有利的扰动方向能先放大多少? | 坐标尺度、范数、时间与输入方向 |
| 吸引域与返回比例 | 一次有限冲击后最终去哪里? | 完整动力学、冲击范围及其概率分布 |
| 全过程存留 | 恢复途中是否越过允许范围? | 允许区域、时间及损伤是否改变动力学 |
即使一个模型最终返回,途中也可能有很大的偏离;如果那种偏离会造成模型没有描述的损伤,单报“最终恢复”仍然不够。5
下面的自建模型均使用固定的无量纲$x,y$坐标,振幅按欧氏范数比较。半径相同表示这组坐标中的位移大小相同,不表示实物输入能量相同。不能把改变变量单位产生的数值变化当作机制改变。
最小过程:一个衰减的变量,还能把另一个变量推走
本库构造以下对照:
$x$有左、右两种稳定状态;$y$是受到冲击后会自行衰减的附加变量。$K$决定$y$对$x$的作用。这是机制示例,不是纳米梁的标定方程,也不是新的普遍定律。
无论$K$取多少,稳定点都是$(-1,0)$与$(1,0)$,中间$(0,0)$是鞍点。在左态附近令$u=x+1$,有:
三组$K=0,2,8$的局部特征值完全相同,都是$-2,-1$。只在$x$方向施加小冲击时,$y$始终为零,三组的$x$恢复过程甚至完全相同。这种测量没有激发另一个输入通道。
换成$y$方向,虽然$y(t)=y_0e^{-t}$一直衰减,它仍会在消失之前持续推$x$。在线性范围内:
右上角表示耦合传递。单独激发$y$时,它对$x$的线性贡献在$t=\log2$达到$Ky_0/4$。这不是增长的本征模,而是两个衰减过程组合成的暂时响应。
按所有单位初始方向取最大振幅比:
本库算得$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$就是非正规、却没有欧氏范数放大的对照。

有限冲击要走完非线性过程
从左态出发,瞬间把$y$改为$+0.3$,随后停止输入。$K=2$最终回左态;$K=8$最终到右态。保持大小不变、改为$y=-0.3$,两组都回左态。只看冲击大小,遗漏了方向。
过程分成三步:附加变量保存短暂偏差;耦合把偏差传给$x$;当$x$的偏离进入非线性范围,原先的线性恢复近似不再负责决定去向,必须继续解完整方程。
为寻找门槛,令$z=Ky$。任意$K>0$的完整方程都变为:
鞍点的稳定流形分隔两种去向。沿该边界有:
本库从鞍点附近积分边界,找到它与$x=-1$相交于$z_*=1.60721816$,所以正$y$方向一次冲击的临界大小为:
$K=2$的门槛约$0.80361$,$K=8$约$0.20090$,同样的$0.3$落在门槛两侧。另用完整方程从初态向前积分并搜索去向分界,24、40、56三个终止时刻的结果均与上述门槛相符,相对差小于$1.2\times10^{-11}$。这检验本例的数值一致性,不是一般证明。
再把冲击写成:
720个等角方向的完整积分给出$K=0,2$全部恢复,$K=8$转态比例$27.50\%$。方向加密至1440个、时间步减半并延长观察,得到$27.4306\%$。独立求圆周与边界交点,转态角区间约$33.2283^\circ$至$131.9778^\circ$,连续角度下比例为$27.4304\%$。
这个比例属于固定半径、均匀方向、只冲击一次的任务,不能读成每单位时间的转态率。
把冲击后的状态分布记为$\mu$,原态吸引域记为$\mathcal B$,返回概率是:
分布也是任务的一部分。保持平均位移平方为$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$只是投影读数改变。持续外部输入或噪声作用下,状态才可能穿过固定无输入系统的边界。

更强的替代解释:局部相同,全局门槛本来就不同
上述比较同时改变短时传递与完整吸引域几何,不能把全部差别只归给一个线性增益。完全没有附加变量,也能恢复率相同、转态门槛不同。考虑另一组自建对照:
稳定点为$0,2a$,分界点为$a$;两个稳定点的局部特征值均为$-1$,与$a$无关。对应势函数为:
从$0$受到大小$0.5$的正$x$冲击,$a=0.25$时到另一稳定点,$a=1$时返回。分界距离差4倍、势垒差16倍,局部恢复率仍一样。这里改变了另一稳定点的位置和全局形状,只匹配局部恢复,不能宣称除了某个机制之外所有条件都相同。
在前一个二维模型中,$K=2$没有线性范数放大,却能被$y=0.9$的一次冲击转态。线性放大既不能单独给转态概率,也不是有限转态必须具备的条件。
原始研究在哪一步支持这条机制
2026年Troude与Sornette的研究尝试把保持局部谱的非正规耦合推进到随机逃逸理论。其出版方接受稿与补充材料已取得;补充说明明确把可分离二维模型视作便于隔离机制的最小情形,非可分离、高维推广列作后续工作。12
本库按固定版本作者代码中grad_U_main的漂移方程另写实现,检验其基础模型环节:
作者用它研究带噪声的四稳态模型;本库关闭噪声,新增一次冲击任务。对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$满足:
直接求解得:
$D=0.01,K=8$时,同样总强度的噪声只输入$x$,$x$方差为$0.005$;只输入$y$则为$0.10667$,相差约21.33倍。解析表达与独立连续Lyapunov方程求解相符。只记录噪声大小还不够,输入必须投影到能被耦合传递的方向。
这仍是局部线性过程的方差,不能把它代入某个门槛就宣称得到非线性转态率。真实跨态还需噪声时间结构、完整路径、事件定义与观测窗口,分别留给后续课题。

怎样把候选机制分开
| 对照 | 区分什么 | 本篇证据的范围 |
|---|---|---|
| 固定全部局部特征值,改变耦合 | 局部谱是否足够,还是需要传递与全局几何 | 自建与原作者漂移模型均有不同去向 |
| 固定冲击半径,改变方向 | 扰动大小是否足够描述任务 | 强耦合时仅部分角区间转态 |
| 固定恢复率,改变标量全局形状 | 非正规放大是否必需 | 无耦合反例否定其必需性 |
| 固定总噪声强度,改变输入方向 | 噪声资源与投影是否混淆 | 局部方差改变;跨态率未算 |
| 同时记录最终状态与允许范围 | 最终恢复是否掩盖沿途损伤 | 两任务需分别定义;未给实物损伤门槛 |
只观测一个变量,可能没有激发隐藏模式;只记录占据比例,不能恢复全部转移速度;只看短时过零,不能确定最终到另一稳定态。两态系统的稳态占据比只约束两方向速率之比:把两个速率同时乘以同一倍数,驻留比例不变,等待时间却改变。
这篇回答到哪里
较稳的部分是:局部恢复谱、短时响应、有限冲击去向与全过程约束是不同读数;完整状态和输入方向不能省略。两个自建对照与原作者漂移模型给出可回查的具体因果链。
本文的整合判断是:**把恢复率用作有限抗扰动指标,需要额外核对扰动集合、输入通道和吸引域;关注沿途损伤时,还要增加允许区域。**这些信息充分时,局部恢复测量可以参与预测;仅有相同恢复率,一般不足以确定风险排序。
未解决的是:这些对照在具体纳米梁参数上对应多大冲击;部分观测下怎样重建边界;持续、带记忆噪声下的事件时间与路径怎样变化;高维系统中哪些低维分界还能保留。本篇不把数值比例搬成现实系统的通用风险。