第 15 章：相平面分析、Routh–Hurwitz 条件、Jury 条件与三次方程根（Appendix Phase Plane Analysis, Routh–Hurwitz Conditions, Jury Conditions, Descartes' Rule of Signs and Exact Solutions of a Cubic）

附录 A 相平面分析（Phase Plane Analysis）

本附录以一段克制的方法论声明开篇：仅对如下形式的二阶自治常微分方程组做最简短的相平面讨论，并仅给出正文中所需要的基本结论；若需更深入的展开，读者可参考 Jordan and Smith (1999) 或 Guckenheimer and Holmes (1983) 的论著；若寻求一段简短而实用的常微分方程定性理论（含相平面方法）介绍，Odell (1980) 的章节是合适的选择。方程组形式为 dx/dt = f(x, y), dy/dt = g(x, y)。所谓相曲线或相轨线是 dx/dy = f(x, y)/g(x, y) 的解；除奇点 (x_s, y_s)（即 f(x_s, y_s) = g(x_s, y_s) = 0）之外，过任意点 (x_0, y_0) 恰有唯一相曲线。利用平移 x → x − x_s, y → y − y_s，可不失一般性地把奇点平移到原点：f(0, 0) = g(0, 0) = 0。

在 f 与 g 于原点解析的条件下，对 (A.2) 做 Taylor 展开并保留线性项，得到线性化形式 dx/dy = (ax + by)/(cx + dy)，对应的线性系统为 dx/dt = ax + by, dy/dt = cx + dy。矩阵 A 写作 |a b; c d| = |f_x f_y; g_x g_y| 在 (0, 0) 处的取值（偏导数均在原点处求值）；其特征值 λ_1, λ_2 由 det(A − λI) = 0 给出，即 λ = (1/2)(a + d ± [(a + d)^2 − 4 det A]^{1/2})。线性系统 (A.5) 的解可以写成 (x, y)^T = c_1 v_1 exp[λ_1 t] + c_2 v_2 exp[λ_2 t] 的形式，其中 c_1, c_2 为任意常数，v_1, v_2 是与 λ_1, λ_2 对应的特征向量（归一化形式为 v_i = (1 + p_i^2)^{−1/2} (1, p_i)^T，其中 p_i = (λ_i − a)/b, b ≠ 0, i = 1, 2）。式 (A.7) 的形式仅对相异特征值成立；若特征值相等（即 λ_1 = λ_2 = λ），则解形如 (c_1 + c_2 t) exp[λt]，即在 v 的线性组合基础上叠加上 t exp[λt] 的修正项。从 (A.7) 中消去 t 即得 (x, y) 平面上的相曲线（参数 t 仅作为辅助参数），其几何形态完全由 λ_1, λ_2 决定。

线性奇点按特征值情形分类（参见图 A.1 的典型相图与图 A.2 的 tr A − det A 判别图）。第一类为 λ_1, λ_2 实且相异：(a) 若 λ_1 与 λ_2 同号（典型情形如 λ_2 < λ_1 < 0），则 (A.7) 中 exp[λ_2 t] 相对于 exp[λ_1 t] 是高阶小量（o(exp[λ_1 t])），例如取 c_2 = 0, c_1 ≠ 0 时解形如 c_1 v_1 exp[λ_1 t]，在相平面上沿 v_1 方向趋于原点（c_1 > 0 时沿 PO, c_1 < 0 时沿 QO），并且 (A.7) 中所有解都最终沿 v_1 方向趋于原点——这是节点（node, Type I）奇点；当 λ_1 ≤ λ_2 < 0 时为稳定节点（所有轨线 t → ∞ 时趋于原点），而 λ_1 > λ_2 > 0 时为不稳定节点（(x, y) 在 t → −∞ 时趋于原点）。(b) 若 λ_1, λ_2 异号（典型情形如 λ_1 < 0 < λ_2），v_1 方向上 t → ∞ 时趋于原点，v_2 方向上 t → −∞ 时趋于原点，因此除严格沿 v_1 之外的小扰动都会指数增长——这是鞍点（saddle point），它始终是不稳定的。第二类为 λ_1, λ_2 = α ± iβ, β ≠ 0：解涉及 exp[αt] exp[±iβt]，对应振荡式趋近或离开原点。(a) α ≠ 0 时为螺旋（spiral）奇点，α < 0 时稳定（相轨迹从外向内盘旋趋于原点），α > 0 时不稳定。(b) α = 0 时相曲线为椭圆，称为中心（centre）；中心在通常意义下并不稳定——线性近似判定为中心时，必须考察比线性更高阶的项才能判定其真实稳定性，因为在线性近似下小扰动只是切到另一条相曲线而非回到原曲线。第三类为 λ_1 = λ_2 = λ：解中含有 t exp[λt] 修正项，只存在一个特征向量 v——称为节点（node, Type II）奇点（与 Type I 的区别在于收敛方向不再沿一个特征方向而更"对称"）；若解中不出现 t exp[λt] 项（特例情形），则得到星形（star）奇点，其稳定性取决于 λ 的符号。

这些奇点类型由矩阵 A 的迹 tr A = a + d 与行列式 det A = f_x g_y − f_y g_x 共同决定（其中偏导数在奇点处求值），其判别汇总在图 A.2 中。Poincaré–Bendixson 定理是这一框架下的关键结论：若 (A.1) 存在一个约束集（即在某区域边界 ∂B 上向量 (dx/dt, dy/dt) 指向该区域内部）包围一个孤立的奇点，而该奇点是不稳定螺旋或不稳定节点，则任何相轨线既不能随 t → ∞ 收敛到奇点（因为奇点是不稳定的），也不能离开该约束集（因为约束集是封闭的）——因此 t → ∞ 时轨线必趋于某个极限环（limit cycle）解。这是该定理最简形式的应用：若唯一的奇点是鞍点，则极限环不可能存在（证明与更一般的应用可参见 Jordan and Smith 1999）。这一判据组合（约束集 + 不稳定奇点 → 极限环；唯一鞍点 → 极限环不存在）构成了第 7-8 章生物振荡器与 BZ 反应分析的几何基础。

附录 B.1 特征多项式、Routh–Hurwitz 条件与 Jury 条件（Characteristic Polynomials, Routh–Hurwitz Conditions and Jury Conditions）

本节从"线性稳定性的判定归结为特征多项式根的位置"这一基本事实出发。这种分析在种群相互作用模型与反应动力学系统（参见第 3 章与第 6 章）中反复出现，因此值得把所需的判据集中陈述。考虑向量形式线性系统 dx/dt = A x，其中 A 是非线性项在稳态处的 Jacobian 矩阵（生态学中称为群落矩阵，community matrix）。设 x = x_0 exp[λt]，代入后得到特征方程 det(A − λI) = 0；解 x = 0 的稳定性要求特征多项式所有根 λ 均落在左半复平面（即 Re λ < 0），此时 x 在 t → ∞ 时指数衰减到零。

对于 n 阶系统，特征多项式可取一般形式 P(λ) = λ^n + a_1 λ^{n−1} + ··· + a_n = 0，其中所有系数 a_i 为实数。若 a_n = 0 则 λ = 0 是平凡根，可等价地把多项式降到 n − 1 阶讨论。在 a_n ≠ 0 的假设下，Routh–Hurwitz 条件给出了一组关于 a_i 的充要判据，使得 P(λ) 的所有零点满足 Re λ < 0。这些条件与 a_n > 0 联合给出，构造方式如下：D_1 = a_1 > 0；D_2 = |a_1, 1; a_3, a_2| > 0；D_3 = |a_1, a_3, a_5; 1, a_2, a_4; 0, a_1, a_3| > 0；D_k 逐级递推，是奇下标系数 a_1, a_3, a_5, … 与偶下标系数 a_2, a_4, … 交叉排列所成行列式（式 (B.5) 中给出 D_k 的全 k × k 形式）。当多项式为三次 λ^3 + a_1 λ^2 + a_2 λ + a_3 = 0 时，Re λ < 0 的充要条件是 a_1 > 0, a_3 > 0, a_1 a_2 − a_3 > 0。作者坦言"老实说很难想象真的有人会把这些条件用到五次及以上的多项式上"——这暗示对高于三阶的特征多项式，更通常的做法是直接用根与系数之间的代数关系（即 Vieta 关系）：若 λ_1, …, λ_n 是 P(λ) 的非零相异根，则 Σ_{i=1}^n λ_i = −a_1, Σ_{i≠j} λ_i λ_j = a_2, …, λ_1 λ_2 … λ_n = (−1)^n a_n（式 (B.6)）。Lewis (1977) 给出了几种利用 Routh–Hurwitz 条件与网络概念直接对矩阵导出定性结论的方法。

Routh–Hurwitz 条件不仅适用于连续时间系统，在离散时间系统（种群相互作用的离散模型，参见第 3 章）与单种群时滞模型（参见第 2 章）中，稳定性同样由特征多项式的根决定（参见第 3 章 3.10 节）。线性化后仍得到形如 (B.3) 的矩阵形式与形如 (B.4) 的多项式，但稳定性的判据不同：对于时滞差分方程（例如时滞 n 个时间步，u_{t+1} = b_1 u_t + ··· + b_n u_{t−n}，第 2 章式 (2.37) 即为一例），设 u_t ∝ λ^t 代入后仍得到 λ 的多项式；此时线性稳定性要求 |λ| < 1（与连续情形的 Re λ < 0 形成对比），这是因为 |λ| < 1 保证 u_t → 0 当 t → ∞。Jury 条件即为该情形下保证多项式所有根的模都小于 1 的充要判据。对于 P(λ) = λ^n + a_1 λ^{n−1} + ··· + a_n，Jury 条件的构造分三阶段递推：(i) 引入 b_n = 1 − a_2, b_{n−1} = a_1 − a_n a_{n−1}, …, b_1 = a_{n−1} − a_n a_1；(ii) 引入 c_n = b_n^2 − b_1^2, c_{n−1} = b_n b_{n−1} − b_1 b_2, …, c_2 = b_n b_2 − b_1 b_{n−1}；(iii) 类似地构造 d 序列（d_n = c_n^2 − c_2^2, …）、r 序列、s 序列，直到只剩三个 s 型元素 s_n = r_n^2 − r_{n−3}^2, s_{n−1} = r_n r_{n−1} − r_{n−3} r_{n−2}, s_{n−2} = r_n r_{n−2} − r_{n−3} r_{n−1}。Jury 条件（必要且充分）要求 |P(1)| > 0, (−1)^n P(−1) > 0, |a_n| > 1, |b_n| > |b_1|, |c_n| > |c_2|, |d_n| > |d_3|, …, |s_n| > |s_{n−2}|（式 (B.7)）。这些条件由 Lewis (1977) 详细描述与示例，他同时讨论了若干关于多项式根的大小与符号的解析与数值技巧。

附录 B.2 Descartes 符号法则（Descartes' Rule of Signs）

给定多项式 P(λ) = λ^n + a_1 λ^{n−1} + ··· + a_n，序列 {a_n, a_{n−1}, …, a_0} 中符号变化的次数（忽略零系数）记为 N。Descartes 法则断言：P(λ) 至多有 N 个正实根，且正实根的个数只能是 N, N−2, N−4, … 中的一个（即 N 与 N 同奇偶且 ≤ N）。该法则的实战价值在于：通过替换 ω = −λ 并再次应用，可以得到关于负实根个数的类似信息。两者合起来常常足以确定所有实根的符号（对稳定性问题而言这就是最关键的信息）。

以 λ^3 + a_2 λ^2 − a_1 λ + a_0 = 0（其中 a_i > 0, i = 0, 1, 2）为例：系数序列 {a_0, −a_1, a_2, 1} 中有两次符号变化（a_0 → −a_1 与 −a_1 → a_2），故正实根的个数为 0 或 2。再令 λ = −ω 代入得到 ω^3 − a_2 ω^2 − a_1 ω − a_0 = 0，序列中仅一次符号变化（−a_0 → ω^3），故 ω 的正实根至多一个。这意味着原方程 (B.8) 恰有一个负实根。

附录 B.3 一般三次多项式的根（Roots of a General Cubic Polynomial）

本节解释为何在数学文献中三次多项式是"我们真正想要写出显式根"的最大阶数：高阶多项式的根虽仍可求，但代价通常不抵其收益。Miura (1980) 与 Namias (1985) 给出了三次根的几种简洁推导。表 B.1 总结了一般三次 p(λ) = λ^3 + A λ^2 + B λ + C（A, B, C 实）的显式根。

引入简记：A ≡ 3a, B ≡ 3b（即将原系数除以 3），α ≡ a^2 − b, β ≡ 2a^3 − 3ab + C。根的形态分七种情形给出：

(1) α > 0 且 β = 0：λ_1 = −a, λ_2 = (3α)^{1/2} − a, λ_3 = −(3α)^{1/2} − a（这是判别式为零的退化情形，其中两复根退化为实根）。

(2) α > 0 且 |β| ≤ 2α^{3/2}：三个实根由正弦参数化给出，λ_1 = 2α^{1/2} sin φ − a, λ_2 = −2α^{1/2} sin(π/3 + φ) − a, λ_3 = 2α^{1/2} sin(π/3 − φ) − a，其中 φ = (1/3) sin^{−1}[β/(2α^{3/2})]，φ 落在 [−π/6, π/6] 区间内——这是三实根的"标准三角度"型表达。

(3) α > 0 且 β > 2α^{3/2}：一个实根与一对共轭复根，λ_1 = −2α^{1/2} cosh ψ − a，ψ = (1/3) cosh^{−1}[|β|/(2α^{3/2})]，λ_2 = −α^{1/2} cosh ψ − a + i(3α)^{1/2} sinh ψ, λ_3 = −α^{1/2} cosh ψ − a − i(3α)^{1/2} sinh ψ。

(4) α > 0 且 β < −2α^{3/2}：对称情形，与 (3) 类似但 ψ 由 cosh^{−1} 给出，结构相同。

(5) α = 0，β 任意：λ_1 = −β^{1/3} − a, λ_2 = β^{1/3}/2 − a + (3i/4) β^{2/3}, λ_3 = β^{1/3}/2 − a − (3i/4) β^{2/3}。

(6) α < 0：引入 (−α)^{1/2}，由双曲正弦参数化给出，θ = (1/3) sinh^{−1}[β/(2(−α)^{3/2})]，λ_1 = −2(−α)^{1/2} sinh θ − a, λ_2 = (−α)^{1/2} sinh θ − a + i(−3α)^{1/2} cosh θ, λ_3 = (−α)^{1/2} sinh θ − a − i(−3α)^{1/2} cosh θ。

这套显式根（连同 Routh–Hurwitz 充要条件一起）使得三次特征方程情形在生物数学应用中具有特殊地位——可读、可算、可解释，而高于三次时一般退而求其次，仅用根的代数关系做定性判断。

本章个人批注

本章与第 14 章风格上截然不同：第 14 章是"分形理论的生物学应用与滥用"——一篇有立场、有方法论警惕、有大量生物学实例（视网膜细胞、亚细胞膜表面、肺血流、空间填充曲线）的论述型章节；本章则是两个附录的汇编——附录 A 是"相平面分析"的方法学速成，附录 B 是"特征多项式稳定性判据"的工具箱。作者在两节附录中的共同意图是：为读者在阅读正文章节时遇到的数学工具提供一份可即时查阅的参考，而不必为每个工具翻找外部教材。这一意图在 A 节开篇就明确宣示："仅对……做最简短的相平面讨论，并仅给出正文中所需要的基本结论"；在 B 节中体现为"老实说很难想象真的有人会把这些条件用到五次及以上的多项式上"这一类坦率评论。我把本章定位为"工具箱附录"而非"叙事章节"——它的内容与前 14 章的具体应用相连接，但自身并不推进任何新的生物学议题。技术细节方面，我特别关注以下几点：(1) 附录 A 的核心是 2D 自治系统奇点的分类：节点（I 型与 II 型）、鞍点、螺旋、中心、星形——所有这些都可由矩阵 A 的迹与行列式唯一决定，其汇总在图 A.2 中。(2) 节点（I 型）vs 节点（II 型）的区别在于：I 型对应两个不同特征值 λ_1 ≠ λ_2 同号，解沿较慢方向趋于原点；II 型对应 λ_1 = λ_2 出现 t exp[λt] 修正项，结构上有本质不同。(3) 中心奇点（α = 0, λ = ±iβ）在线性近似下稳定但实际未必——这是 Poincaré–Lindstedt 思想的雏形：必须考察高阶项才能判定真实稳定性。(4) Poincaré–Bendixson 定理给出极限环存在的最简充分条件：约束集内仅有一个不稳定螺旋或不稳定节点奇点；这一条件在 BZ 反应、心脏电生理、捕食者-猎物等若干实例中得到应用（参见第 7-8 章）。我注意到附录 A 末尾明确说明"若唯一的奇点是鞍点则极限环不可能存在"——这一对偶判据在排除某些"假想极限环"时非常实用。(5) 附录 B.1 的核心是 Routh–Hurwitz 充要条件族（D_k 行列式 > 0, k = 1, …, n）与 Jury 条件族（|a_n| > 1, |b_n| > |b_1| 等递推不等式），分别针对连续时间（Re λ < 0）与离散时间（|λ| < 1）系统。(6) Vieta 关系（Σ λ_i = −a_1 等）作为高阶多项式定性判断的辅助工具：作者坦诚"五次及以上用 Routh–Hurwitz 实在麻烦"，此时 Vieta + Descartes 法则往往已经足够。(7) Descartes 法则（N = 符号变化次数，正实根数 ∈ {N, N−2, …}）是一种"宽松但快速"的实根数上界估计；与 Vieta 联用可确定实根符号（这是稳定性判别的全部所需信息）。(8) 三次方程显式根（表 B.1）的七种情形分类依据两个组合不变量 α ≡ a^2 − b 与 β ≡ 2a^3 − 3ab + C：α > 0 时对应三实根或一实两复；α = 0 时为临界情形；α < 0 时对应一实两复。这一分类本质上是"判别式 Δ = 18abc − 4b^3 + b^2 c^2 − 4c^3 − 27a^2"在归一化变量下的不同符号情形。我个人对本章作为"工具箱"的印象比较深：它不是新内容的载体，而是把分散在正文章节中用到的数学工具集中起来——相平面法是第 2 章性别决定、第 3 章相互作用种群、第 6 章反应动力学、第 7 章生物振荡器所共用的几何语言；Routh–Hurwitz 是第 3 章 Lotka–Volterra 与第 6 章反应系统稳定性分析的基础；Jury 条件是第 2 章离散时滞模型与第 3 章差分方程稳定性分析的工具；Descartes 与三次显式根则是当 Routh–Hurwitz 在三次情形下需要数值化、或者只需要"几个实根"时的高效替代。这种"工具分散使用 + 附录集中参考"的结构，使得附录本身短而精——每个工具的呈现都聚焦于"判定什么 + 怎么算"，不展开证明与应用案例（那些都在正文章节中）。

与上下章的衔接（一段话）

第 14 章是第三部分"波与空间中的图式"中关于分形几何的方法论收束——它处理的是"在越来越小的尺度上重复自身的几何对象如何与生物学测量相关"；本章则跳出具体生物学议题，转向"工具箱"性质的数学附录——附录 A 提供相平面方法（奇点分类、轨迹几何、Poincaré–Bendixson 极限环存在性），附录 B.1 提供连续/离散系统的多项式稳定性判据（Routh–Hurwitz 与 Jury 条件），附录 B.2 提供快速判断实根个数的 Descartes 法则，附录 B.3 给出三次特征方程的显式根。从结构上看，本章是 Murray 第一卷正文章节（1-14 章）与全书参考文献之间的桥梁——它不复述任何具体模型，而是把正文中反复使用的几种数学工具以紧凑形式集中呈现，供读者在阅读第 2 章性别决定（奇点稳定性）、第 3 章相互作用种群（特征方程）、第 6 章反应动力学（极限环与稳定性）、第 7 章生物振荡器（中心与螺旋奇点）、第 8 章 BZ 反应（极限环实例）等章节时即时查阅。位置上看，本章与第 14 章并列位于全书末尾的"附加材料"部分，区别于第 14 章对单一议题（分形）的方法论阐述，本章是跨议题的工具汇总。从上下衔接看：本章向上承接第 14 章——后者是 Murray 第一卷正文章节中最后一个有完整叙事结构的章节（其个人批注中我曾特别关注 Murray 对"机制 vs 现象"的方法论警示）；本章向下开启参考文献（Bibliography），从 [1] AIDS 到 [432] 及之后共约 432 条文献，覆盖了从经典（Euler 1760、Malthus 1798、Bernoulli 1760）到现代（Nelson–Perelson–Murray 2000 的 HIV 延时模型）的所有主要参考文献。Murray 在本章末尾用"This page intentionally left blank"分隔附录 B 与参考文献，是 Springer 风格出版物的标准排版——这一分隔也提示我们附录 B 是正文的"附属技术材料"，参考文献则是独立的引文索引。