主稳定函数与可同步性

对于相同的振荡器（如相同的神经元模型），完全同步态总是存在的，问题是：它稳定吗？Pecora 与 Carroll（1998，PRL）提出主稳定函数（master stability function, MSF）方法。设

$$\dot{\mathbf x}_i=\mathbf F(\mathbf x_i)-\sigma\sum_j L_{ij}\,\mathbf H(\mathbf x_j)$$

$\mathbf x_i\in\mathbb R^n$，$\mathbf F$ 是单个节点的动力学，$\mathbf H$ 是耦合函数（决定“通过哪些变量耦合”），$\sigma$ 是全局耦合强度，$L$ 为拉普拉斯矩阵，特征值 $0=\mu_1<\mu_2\le\cdots\le\mu_N$（连通无向网络）。因为 $L$ 行和为零，同步流形 $\mathbf x_1=\cdots=\mathbf x_N=\mathbf s(t)$（$\dot{\mathbf s}=\mathbf F(\mathbf s)$）是不变集。

核心思想

在同步态附近线性化，扰动可以沿 $L$ 的特征向量分解成 $N-1$ 个互不耦合的横向模式，第 $k$ 个模式满足

$$\dot{\boldsymbol\eta}_k=\big[D\mathbf F(\mathbf s)-\sigma\mu_k\,D\mathbf H(\mathbf s)\big]\boldsymbol\eta_k .$$

所有模式都是同一个方程，只是参数 $\alpha=\sigma\mu_k$ 不同。于是定义主稳定函数 $\Lambda(\alpha)$：主方程 $\dot{\boldsymbol\zeta}=[D\mathbf F(\mathbf s)-\alpha D\mathbf H(\mathbf s)]\boldsymbol\zeta$ 的最大李雅普诺夫指数。

MSF 的三种类型

按 $\Lambda(\alpha)<0$ 的区域形状（Boccaletti 等 2006 的归纳）：

1. 无稳定区：$\Lambda(\alpha)>0$ 处处成立，任何网络、任何耦合强度都不能使同步稳定（耦合方式选错了变量）；
2. 无界稳定区 $(\alpha_1,\infty)$：只需 $\sigma\mu_2>\alpha_1$，$\mu_2$ 越大越好，所需最小耦合 $\sigma_{\min}=\alpha_1/\mu_2$；
3. 有界稳定区 $(\alpha_1,\alpha_2)$：需要 $\sigma\mu_2>\alpha_1$ 且 $\sigma\mu_N<\alpha_2$。存在这样的 $\sigma$ 当且仅当 $$\frac{\mu_N}{\mu_2}<\frac{\alpha_2}{\alpha_1}.$$ 耦合太弱，慢模式不稳定；耦合太强，快模式（短波长模式）反而失稳。

例子

例 1（可手算，第 2 类）：标量节点 $\dot x=ax$（$a>0$，单独时不稳定），耦合 $H(x)=x$。主方程 $\dot\zeta=(a-\alpha)\zeta$，$\Lambda(\alpha)=a-\alpha$，稳定区为 $(a,\infty)$。取 $a=1$：3 节点路径 $\mu_2=1$，需要 $\sigma>1$；三角形 $\mu_2=3$，只需 $\sigma>1/3$；$N=5$ 的路径 $\mu_2\approx0.382$，需要 $\sigma>2.6$。

例 2（数值，第 3 类）：Rössler 振子（$a=b=0.2,\ c=7$，混沌），只通过 $x$ 分量耦合（$D\mathbf H=\mathrm{diag}(1,0,0)$）。本站代码对每个 $\alpha$ 积分变分方程并求最大李雅普诺夫指数，得到 $\Lambda(0)\approx+0.07$（即孤立 Rössler 振子的最大李雅普诺夫指数，混沌），稳定区约为 $(0.14,\,4.7)$，$\alpha_2/\alpha_1\approx35$（有限积分时长，数值仅作示意；参数不同时区间会变化）。于是 $N=100$、平均度 6 的规则环（$p=0$ 的 WS 网络）$\mu_N/\mu_2\approx156$，无论 $\sigma$ 取多少都无法同步；重连概率 $p=0.05$ 时降到约 34，处于边缘；$p=0.2$ 时约 19，可以同步。

由此得到的结论

• 小世界结构利于同步：在环形规则网络中加入少量捷径，$\mu_N/\mu_2$ 急剧下降（Barahona 与 Pecora 2002，PRL）。长程连接让远处的节点能直接“对表”。
• 异质性悖论：Nishikawa 等（2003，PRL）发现，度分布越异质的网络，尽管平均路径更短，按特征值比衡量却更难同步——枢纽抬高 $\mu_N$（$\mu_N\ge k_{\max}+1$），而大量低度节点压低 $\mu_2$。对耦合按度加权（如每个节点的输入除以其度）可以缓解这一问题（Arenas 等 2008 第 4 章）。
• 注意区分：度异质性使 Kuramoto 模型更早开始部分同步（$K_c$ 低），却可能使相同振荡器的完全同步更难稳定。两个结论回答的是不同问题。

完整推导见数学推导与算法 C 节，数值算法见算法 2。