1. 课程导入与分岔现象初识

在工程与物理世界里，有一类现象反复出现：系统参数发生极其微小的连续变化，系统的行为却突然发生质的改变。一根细长压杆在轴向载荷缓慢增加时，起初只是轻微压缩，当载荷越过某个临界值，压杆会突然向一侧大幅弯曲；飞机机翼在气流速度逐渐提高时，原本平稳的升力会突然转化为越来越剧烈的自激振动；一个由运算放大器和电阻电容构成的双稳态电路，在输入电压缓慢扫过阈值时，输出电压会从一个稳定值瞬间跳变到另一个稳定值。这些现象的共同名字，就是分岔。

分岔（bifurcation）指的是：当系统参数连续变化并穿过某个临界值时，系统的定性行为——平衡点的数目、稳定性、周期解的存在性——发生突变。理解这句话需要抓住两个关键词。第一是“参数”，它是我们可以外部调节的量，例如载荷、流速、增益、阻尼系数；第二是“定性行为”，它不是数值大小的变化，而是结构的变化，比如原来只有一个稳定平衡点，越过临界值后变成两个稳定平衡点加一个不稳定平衡点。分岔分析要回答的核心问题正是：临界值在哪里？越过临界值后系统变成了什么结构？这种变化对工程安全意味着什么？

为了把这三个问题落到可操作的层面，本课程采用“现象—模型—判据—计算—验证”的主线。现象提供问题来源，模型把物理对象抽象为含参数的非线性微分方程或代数方程，判据给出识别分岔的数学条件，计算给出临界参数与分岔分支的数值方法，验证则通过数值延拓与相图、时程仿真确认结论。下面这张图概括了本课程的整体学习路径与各环节之间的依赖关系。

先看屈曲这个最经典的静力分岔。把一根两端铰支、长度为 L、抗弯刚度为 EI 的理想压杆作为对象，设轴向压力为 P。按小挠度理论，杆的横向挠度 w(x) 满足 EI·w'''' + P·w'' = 0。代入铰支边界条件 w(0)=w(L)=0、w''(0)=w''(L)=0，可解得非零解存在的条件为 P = n²π²EI/L²。取 n=1 得到临界载荷 P_cr = π²EI/L²。这个结果的物理含义非常明确：当 P < P_cr 时，只有 w≡0 的直线平衡态且稳定；当 P 越过 P_cr，直线平衡态失稳，同时出现两个对称的弯曲平衡分支。这正是叉式分岔的典型结构——一个稳定解分裂为两个稳定解加一个不稳定解。工程上，P_cr 就是设计许用载荷的依据，而分岔分析给出的不只是这个数值，还有“越过之后系统会走向哪一支”的结构信息。

再看颤振，它属于动态分岔。以典型二元机翼截面为例，取沉浮位移 h 与俯仰角 α 为广义坐标，气动升力和力矩随来流速度 U 变化，可写成二阶系统 M·q̈ + C·q̇ + K(U)·q = 0，其中 q = [h, α]ᵀ。把它化为一阶形式后，稳定性由系统矩阵的特征值决定。当 U 较小时，所有特征值实部为负，扰动衰减；随着 U 增大，某一对共轭复特征值的实部由负变正、在 U = U_f 处恰好穿过虚轴，系统从衰减振荡转为等幅振荡再转为发散振荡。这个穿越点就是Hopf 分岔点，U_f 即颤振临界速度。与屈曲相比，颤振的危险性在于分岔后出现的是极限环或发散振荡，属于动态失稳，必须通过特征值实部符号来判断，而不能只看平衡点个数。

双稳态电路则展示了另一种机制。以含正反馈的运算放大器电路为例，输出电压 u_o 与输入 u_i 之间满足一条 S 形的静态特性曲线，中间段斜率为负。当输入电压从低到高缓慢扫描时，工作点沿下支移动，到达上折点后无法继续，只能跳变到上支；反向扫描时则在上折点跳变回下支。两个折点对应鞍结分岔，而两条路径不重合形成的回线就是迟滞。这个例子说明：分岔不仅决定稳定性，还决定了系统对参数变化路径的依赖，即历史相关性。

三个例子放在一起，可以提炼出本课程要建立的三层能力。第一层是建模能力：把物理对象写成 ẋ = f(x, μ) 的形式，明确状态变量 x 与分岔参数 μ 各是什么。第二层是判据能力：会求平衡点、会算雅可比矩阵、会判断特征值是否穿越虚轴或零，从而区分鞍结、叉式、Hopf 等基本分岔类型。第三层是计算与验证能力：会用数值延拓追踪解分支，会画分岔图、相图和时程曲线，并能把数值结果与理论临界值对照，判断计算是否可信。