系统动力学状态 (Ecosystem State)
渐近稳定 (STABLE)
所有特征值均有负实部。微小扰动会被系统阻尼消化并恢复平衡。
最大特征值实部 $\max \text{Re}(\lambda)$
-0.500
s-1
理论近似主值 $\lambda_{\max}$
-0.485
s-1
物种生态位轴分布 (Resource Niche Axis $\xi$)
双击画布添加物种，选中后可拖拽圆点修改 $\mu_i$，右侧滑块修改其参数
生态群落矩阵 (Community Matrix $M$)
Symmetric
相互作用强度:
-1.0
0.0
特征值谱平面 (Complex Eigenvalue Spectrum)
10 Eigenvalues
稳定区 ($\text{Re}(\lambda) < 0$)
不稳定区 ($\text{Re}(\lambda) \ge 0$)
系统特征值
模型数学机理与学术背景 (Mathematical Context)
1. MacArthur-Levins 生态位重叠
物种 $i$ 和 $j$ 在一维资源生态位轴 $\xi$ 上的资源利用函数 $P_i(\xi), P_j(\xi)$ 通常建模为高斯分布：
$$P_i(\xi) = \exp \left[ -\frac{(\xi - \mu_i)^2}{2\sigma_i^2} \right]$$
两个物种由于共享资源而产生的竞争/重叠系数 $\alpha_{ij}$ 定义为其高斯资源分布卷积：
$$\alpha_{ij} = \sqrt{\frac{2\sigma_i \sigma_j}{\sigma_i^2 + \sigma_j^2}} \exp \left[ -\frac{(\mu_i - \mu_j)^2}{2(\sigma_i^2 + \sigma_j^2)} \right]$$
当所有物种生态位宽度恒定为 $\sigma_i = \sigma$ 且均匀相隔 $d$ 时，其重叠系数简化为：
$$\alpha_{ij} = \exp\left[ -\frac{d^2(i - j)^2}{4\sigma^2} \right]$$
这自然导出了一个对称的 Toeplitz（托普利茨）矩阵。
2. 群落动力学与稳定性判据
生态系统的 Lotka-Volterra 演化方程为：
$$\frac{d x_i}{d t} = x_i \left( r_i - s_i x_i + \sum_{j \neq i}^S A_{ij} x_j \right)$$
系统在平衡共存态附近的稳定性由雅可比矩阵（群落矩阵）$M$ 决定：
$$M_{ii} = -s_i, \quad M_{ij} = \pm \alpha_{\text{scale}} \alpha_{ij} \quad (i \neq j)$$
根据克莱普诺夫 (Lyapunov) 稳定性理论，系统渐近稳定的充分必要条件是 $M$ 的所有特征值的实部均小于零：
$$\max \text{Re}(\lambda(M)) < 0$$
若存在特征值实部 $\ge 0$，系统在平衡点将失去稳定性，扰动将被指数放大，导致系统发散或部分物种灭绝。
3. 论文核心结论与缩放规律
-
谱主特征值解析解：对于均匀间隔的 Toeplitz 群落矩阵，在连续极限下，其最大特征值可以解析逼近为：
$$\lambda_{\max} \approx -s - \text{type\_sign} \cdot \alpha_{\text{scale}} + \alpha_{\text{scale}} \frac{2\sigma\sqrt{\pi}}{d}$$这表明，生态系统稳定性主要由比值 $\sigma/d$ 调控，随着生态位重叠的增强（$\sigma$ 增大或距离 $d$ 缩短），系统显著向不稳定性转变。
- 异质分布的稳定效应：当物种分布中心 $\mu_i$ 在轴上呈指数密集堆叠（代表食物网低营养级物种密集）时，相比线性均匀分布，能够大幅降低大矩阵的谱半径，增加系统弹性。
- 人口统计丰度（Demographic Abundance）扰动：考虑共存平衡丰度异质性后，群落矩阵修正为 $M_X = X M$，其中对角矩阵 $X = \text{diag}(x_i^*)$ 引入的波动方差 $\sigma_X^2$ 会作为二阶微扰项，系统性地增大最大特征值实部，加速系统失稳。