1. 课程导入与动力系统建模基础

在工程、物理、生物与经济领域，我们经常遇到一类问题：系统状态会随时间演化，而我们关心的往往不是它在某一瞬间的具体数值，而是它长期会走向哪里、会不会停在某个位置、受到扰动后能否回到原处。研究这类问题的数学工具就是动力系统。本课程的核心对象只有一个：平衡点及其稳定性。理解了这个对象，就掌握了分析大量实际系统长期行为的主线。

一、什么是动力系统

动力系统可以理解为“状态 + 演化规则”的组合。设系统的状态用一个向量 x 描述，演化规则写成微分方程：

dx/dt = f(x),   x ∈ R^n

其中 t 是时间，f 是给定的向量场。给定初始状态 x(0)，方程的解 x(t) 就描述了系统状态随时间的轨迹。这里需要先明确一个约定：本课程统一采用时间正向演化，即 t 从 0 增大，所有“稳定”“收敛”均指 t → +∞ 时的行为；涉及一维位置量时，坐标轴正方向按图中标注为准，全文与图示保持一致。

动力系统分为连续时间（微分方程）与离散时间（差分方程 x_{k+1} = f(x_k)）两类。本课程以连续时间为主，离散情形在方法上完全平行。

二、平衡点的定义与物理意义

平衡点（也称不动点、奇点）是使系统“停住”的状态：若 x* 满足

f(x*) = 0

则从 x* 出发的系统永远停留在 x*，即 x(t) ≡ x*。物理上，平衡点对应系统的静态工作点：单摆的竖直下垂位置、种群数量的环境容纳量、电路的工作偏置点、控制系统的设定值。但“能停住”不等于“停得稳”。真实系统总有扰动，因此真正重要的是：平衡点附近的轨迹是被吸回 x*，还是被推离 x*。这就是稳定性问题，也是本课程后续所有章节要回答的问题。

三、从实际对象建立模型：两个典型案例

案例一：单摆。设摆长 l、质量 m、重力加速度 g，摆角 θ 从竖直向下位置量起，忽略阻尼。由力矩平衡得

m l² d²θ/dt² = -m g l sinθ
即  d²θ/dt² + (g/l) sinθ = 0

令 x₁ = θ、x₂ = dθ/dt，化为一阶系统：

dx₁/dt = x₂
dx₂/dt = -(g/l) sin x₁

令右端为零，得 x₂ = 0 且 sin x₁ = 0，即平衡点为 θ = kπ。其中 θ = 0（下垂）与 θ = π（倒立）在物理直觉上截然不同：前者受扰后会摆动回来，后者一碰就倒。这个直觉正是稳定性概念的来源。取 g = 9.8 m/s²、l = 1 m（示例参数），则 g/l ≈ 9.8，在 θ = 0 附近用 sinθ ≈ θ 线性化得 d²θ/dt² + 9.8θ = 0，特征根为纯虚数 ±j·3.13，对应无阻尼等幅振荡——这是“临界情形”，加入阻尼后才会真正收敛。

案例二：种群增长。设种群数量 N，内禀增长率 r，环境容纳量 K，采用逻辑斯蒂模型：

dN/dt = r N (1 - N/K)

令右端为零得 N* = 0 与 N* = K。取示例参数 r = 0.8 /单位时间、K = 1000：N* = 0 表示灭绝（不稳定，少量个体即可增长），N* = K 表示饱和（稳定，数量偏离后会回归）。在 N = K 处线性化，令 N = K + u，得 du/dt ≈ -r u，特征值为 -0.8 < 0，故稳定；在 N = 0 处线性化得 dN/dt ≈ rN，特征值 +0.8 > 0，故不稳定。

四、本节知识结构总览

下图把“实际对象 → 微分方程模型 → 平衡点 → 稳定性判据”的建模链条与两个案例的具体结果放在一起，便于对照理解。

五、建模时容易踩的坑

• 平衡点求全。令 f(x)=0 时容易漏解，尤其含三角函数、乘积项时。单摆的 θ = kπ 有无穷多个，只取 θ = 0 会漏掉倒立点。
• 先化一阶系统。二阶及以上方程必须先引入状态变量降阶，否则无法统一用向量场语言讨论平衡点。
• 线性化只在平衡点附近有效。在 θ = 0 用 sinθ ≈ θ 合理，在 θ = π 附近应令 θ = π + u，此时 sinθ ≈ -u，线性化方程变为 u'' - 9.8u = 0，特征值 ±3.13，一正一负，正对应“一碰就倒”。
• 参数要带单位与来源。示例中的 g=9.8、l=1、r=0.8、K=1000 均为教学假设值，实际建模须用实测或标定参数。