3. 一维离散动力系统的平衡点与稳定性

连续系统用微分方程描述状态随时间的连续变化，而离散动力系统用映射描述状态在离散时刻之间的跳变。一维离散系统的一般形式为

x_n+1 = f(x_n)，其中 f 为实数区间上的连续可微函数。

本节要解决三个问题：平衡点如何定义、如何用蛛网图和线性化判断稳定性、以及这些工具在人口模型与混沌映射中的实际表现。

平衡点的定义与求解

若某点 x* 满足 f(x*) = x*，则称 x* 为映射 f 的平衡点（不动点）。一旦系统到达 x*，后续所有时刻都停留在该点。求解平衡点就是解方程 f(x) − x = 0，通常得到若干候选点，再结合实际问题背景（如人口非负）筛选出有意义的平衡点。

需要注意，离散系统的平衡点与连续系统不同：连续系统平衡点要求 f(x*)=0，而离散系统要求映射值等于自身。这一差别决定了后续稳定性判据的形式。

线性化判据：导数决定局部稳定性

设 x* 为平衡点，在 x* 附近作一阶展开：

x_n+1 − x* ≈ f′(x*)(x_n − x*)

记偏差 e_n = x_n − x*，则 e_n+1 ≈ λ e_n，其中 λ = f′(x*) 称为乘子。由此得到判据：

• |λ| < 1：平衡点局部渐近稳定，偏差按几何级数衰减；
• |λ| > 1：平衡点不稳定，偏差被放大；
• |λ| = 1：线性化失效，需借助高阶项或图形方法判断，可能出现临界情形。

λ 的符号还决定趋近方式：0 < λ < 1 时偏差单调衰减，系统从同侧逼近；−1 < λ < 0 时偏差正负交替，系统在平衡点两侧振荡收敛。λ 越接近 0，收敛越快。

蛛网图：把迭代过程画出来

蛛网图（cobweb plot）是分析一维映射最直观的工具。在同一坐标系中画出曲线 y = f(x) 与对角线 y = x，从初值 x₀ 出发：竖直走到曲线上得到 f(x₀)，水平走到对角线上把该值转为新的横坐标，再竖直走到曲线……如此反复，形成阶梯状折线。折线与对角线的交点即平衡点；折线是螺旋收敛、阶梯收敛还是向外发散，直接反映稳定性。

蛛网图的价值在于它同时展示了平衡点位置、稳定性类型和收敛速度，尤其适合 |λ| 接近 1 或线性化失效的情形。

案例一：人口模型（Logistic 映射）

经典离散人口模型为 Logistic 映射：

x_n+1 = r x_n(1 − x_n)，x_n ∈ [0,1]，r > 0 为增长率参数。

令 f(x) = rx(1−x)，解 f(x) = x 得两个平衡点：x*₁ = 0（灭绝态）与 x*₂ = 1 − 1/r（正平衡态，仅当 r > 1 时落在 (0,1) 内）。

求导得 f′(x) = r(1 − 2x)。在 x*₁ = 0 处 λ = r，故 r < 1 时灭绝态稳定，r > 1 时不稳定。在 x*₂ = 1 − 1/r 处 λ = r(1 − 2(1 − 1/r)) = 2 − r，因此：

• 1 < r < 3 时 |2 − r| < 1，正平衡点稳定；其中 1 < r < 2 时 λ > 0 单调收敛，2 < r < 3 时 λ < 0 振荡收敛；
• r > 3 时 |2 − r| > 1，正平衡点失稳，系统进入周期振荡乃至混沌。

以 r = 2.8 为例，λ = 2 − 2.8 = −0.8，属于振荡收敛。取初值 x₀ = 0.2，手算前几步：x₁ = 2.8×0.2×0.8 = 0.448，x₂ = 2.8×0.448×0.552 ≈ 0.692，x₃ ≈ 0.597，x₄ ≈ 0.674，x₅ ≈ 0.615……数值在 x* = 1 − 1/2.8 ≈ 0.643 两侧来回摆动并逐步收拢，与 λ < 0 的预测一致。若把 r 调到 3.2，λ = −1.2，摆动幅度反而增大，蛛网图呈向外螺旋，说明平衡点已失稳。

案例二：混沌映射与稳定性判据的边界

Logistic 映射在 r = 4 时化为混沌映射 x_n+1 = 4x_n(1 − x_n)。此时正平衡点 x* = 3/4，λ = 2 − 4 = −2，|λ| > 1，平衡点不稳定。蛛网图上迭代轨迹不再收敛到任何单点，而是在 [0,1] 内密集游走，对初值极其敏感——两个相差 10⁻⁶ 的初值经过数十步后轨迹完全分离。这正是“局部线性化判据只能判断平衡点邻域行为”的典型体现：判据告诉我们 x* 不稳定，但系统最终去向（周期、混沌还是发散）需要借助分岔图、Lyapunov 指数等工具进一步分析。

实操要点与常见误区

实际分析时，建议把“解析求平衡点—计算乘子—画蛛网图验证”作为固定动作。蛛网图不仅能验证结论，还能暴露初值选取不当导致的误判：例如初值落在不稳定平衡点附近时，前几步看似收敛，随后才发散，只看数值序列容易误判，而蛛网图能一眼看出轨迹走向。

常见误区有三：其一，把连续系统的 f(x*)=0 误用于离散系统；其二，忽略 λ 的符号，只判断收敛与否而漏掉振荡特征；其三，在 |λ|=1 时仍强行用线性化下结论。对 Logistic 映射而言，r=3 正是 λ=−1 的临界点，此时线性化失效，系统由稳定平衡点转为周期 2 振荡，必须借助分岔分析才能说清。

掌握平衡点求解、乘子判据与蛛网图这三件工具，就具备了对一维离散系统做局部稳定性分析的基本能力，也为后续高维离散系统与分岔理论打下基础。