多物種生態動力學：非線性互動、分岔與韌性

封存文章

模型參數及圖表主要為合成示範，不能直接推斷某個真實生態系統的承載量或保育門檻。

多物種系統的難處，在於競爭、捕食、互利和季節環境同時作用。一般形式可寫為

\dot N_i =N_i\left[ r_i(t) -\sum_j\alpha_{ij}N_j \right]+I_i(t)-H_i(t).

$\alpha_{ij}$ 的正負表示不同互動，但必須與實際機制一致；任意設定一個矩陣只會產生一個數學故事。

季節與平衡

若 $r_i(t)$ 週期變化，系統通常沒有固定平衡點。應尋找週期軌道或計算長時間週期平均：

\bar N_i =\frac1P\int_{t_0}^{t_0+P}N_i(t)\,dt.

在尚未消除暫態前取平均會產生假結果，因此要檢查多個初始條件及整合時間。

對自治模型的平衡 $N^*$ ，Jacobian

J_{ij}= \left.\frac{\partial f_i}{\partial N_j}\right|_{N^*}

的特徵值實部決定局部穩定性。沿環境或收穫參數延續平衡，可辨識鞍結、Hopf 或跨臨界分岔。只畫參數掃描散點不等於完成分岔分析；要追蹤解支及其穩定性。

擾動可改變物種數量、增長率或棲地容量。韌性至少有三個不同維度：

例如

\tau_\varepsilon =\inf\{t>t_p: \|N(t)-N_{\mathrm{ref}}(t)\|<\varepsilon\}.

若存在多穩態，單看最終總生物量會錯過群落組成已永久改變。

數值上要保持 $N_i\ge0$ 、做步長收斂及質量/能量來源檢查。科學上則需以獨立時間序列或操控實驗估計參數，檢查可辨識性，並對過度參數化採用正則化或階層模型。

模型較適合比較「哪一種機制會造成脆弱性」，而不是輸出單一精確族群值。保育方案應在參數後驗或多模型集合下測試，找出對假設最不敏感的措施。