二维非线性系统
x˙=f(x,y),y˙=g(x,y)
很少能显式求解,但零增长线、局部线性化和相平面几何通常足以勾勒全局行为。
零增长线
- x-零增长线:f(x,y)=0,其上向量竖直;
- y-零增长线:g(x,y)=0,其上向量水平;
- 两者交点:不动点。
零增长线把平面分区。在每个区内判断 f,g 的符号,就知道向量指向右上、左上、左下还是右下。
局部线性化
设 x∗ 是不动点,令 η=x−x∗。Taylor 展开:
η˙=J(x∗)η+O(∥η∥2),
其中
J(x,y)=(fxgxfygy).
若雅可比的所有特征值实部都不为零,不动点称为双曲不动点。Hartman–Grobman 定理保证非线性系统在局部与线性化系统拓扑等价,所以可以沿用结点、焦点、鞍点分类。
一个完整例子
考虑
x˙=2x(1+x2−2y2),y˙=−y(1−4x2+3y2).
零增长线是
x=0或1+x2−2y2=0,
y=0或1−4x2+3y2=0.
交点为原点和 (±1,±1) 四点。雅可比矩阵
J=(2+6x2−4y28xy−8xy−1+4x2−9y2).
在原点为 diag(2,−1),所以是鞍点;在四个非零点处,迹负、行列式正、判别量负,所以都是稳定焦点。
这展示了考试标准流程:先解交点,再逐点代雅可比,不能用一个点的矩阵代表全部。
非双曲点为什么麻烦
如果存在零实部特征值,线性项不能控制全部方向。例如
x˙=−y−x3,y˙=x
在原点的线性化是中心,但
dtd(x2+y2)=2x(−y−x3)+2yx=−2x4≤0.
非线性阻尼使轨道缓慢向内,因此原点实际是稳定焦点式吸引子,而不是闭轨中心。
常见处理方法:
- 转极坐标看 r˙;
- 找守恒量或 Lyapunov 函数;
- 利用对称性、可逆性;
- 必要时使用中心流形或正规形理论。
极坐标公式
对 r2=x2+y2,
r˙=rxx˙+yy˙,θ˙=r2xy˙−yx˙.
极坐标能把“向内/向外”和“旋转”分开,尤其适合中心、焦点、极限环和 Hopf 分岔。
局部结论与全局相图
线性化只描述不动点附近。要补全全局相图,还应检查:
- 零增长线与各区域箭头;
- 不变直线、对称轴;
- 是否有轨道逃向无穷远;
- 是否存在闭轨或捕获域;
- 鞍点稳定/非稳定流形可能连接到哪里。
数值模拟的正确位置
数值相图适合检查草图、发现候选结构,但不能独立证明:
- 看似闭合的轨道可能只是缓慢螺旋;
- 有限时间图看不出极慢的不稳定;
- 步长过大会制造假极限环。
最稳妥的结构是“解析判据给结论,数值图作验证”。
考试要点
- 每个不动点都单独计算 J。
- 双曲点才可直接相信线性分类。
- 纯虚、零特征值、重根边界都要注明线性化失效。
- 零增长线上只有一个分量为零,不代表整条线都是不动点。