非线性相平面与局部线性化

Views: --

二维非线性系统

x˙=f(x,y),y˙=g(x,y)\dot x=f(x,y),\qquad \dot y=g(x,y)

很少能显式求解,但零增长线、局部线性化和相平面几何通常足以勾勒全局行为。

零增长线

  • xx-零增长线:f(x,y)=0f(x,y)=0,其上向量竖直;
  • yy-零增长线:g(x,y)=0g(x,y)=0,其上向量水平;
  • 两者交点:不动点。

零增长线把平面分区。在每个区内判断 f,gf,g 的符号,就知道向量指向右上、左上、左下还是右下。

局部线性化

x\boldsymbol x^* 是不动点,令 η=xx\boldsymbol\eta=\boldsymbol x-\boldsymbol x^*。Taylor 展开:

η˙=J(x)η+O(η2),\dot{\boldsymbol\eta} =J(\boldsymbol x^*)\boldsymbol\eta +O(\|\boldsymbol\eta\|^2),

其中

J(x,y)=(fxfygxgy).J(x,y)= \begin{pmatrix} f_x&f_y\\ g_x&g_y \end{pmatrix}.

若雅可比的所有特征值实部都不为零,不动点称为双曲不动点。Hartman–Grobman 定理保证非线性系统在局部与线性化系统拓扑等价,所以可以沿用结点、焦点、鞍点分类。

一个完整例子

考虑

x˙=2x(1+x22y2),y˙=y(14x2+3y2).\dot x=2x(1+x^2-2y^2),\qquad \dot y=-y(1-4x^2+3y^2).

零增长线是

x=01+x22y2=0,x=0\quad\text{或}\quad1+x^2-2y^2=0, y=014x2+3y2=0.y=0\quad\text{或}\quad1-4x^2+3y^2=0.

交点为原点和 (±1,±1)(\pm1,\pm1) 四点。雅可比矩阵

J=(2+6x24y28xy8xy1+4x29y2).J= \begin{pmatrix} 2+6x^2-4y^2&-8xy\\ 8xy&-1+4x^2-9y^2 \end{pmatrix}.

在原点为 diag(2,1)\operatorname{diag}(2,-1),所以是鞍点;在四个非零点处,迹负、行列式正、判别量负,所以都是稳定焦点。

这展示了考试标准流程:先解交点,再逐点代雅可比,不能用一个点的矩阵代表全部。

非双曲点为什么麻烦

如果存在零实部特征值,线性项不能控制全部方向。例如

x˙=yx3,y˙=x\dot x=-y-x^3,\qquad \dot y=x

在原点的线性化是中心,但

ddt(x2+y2)=2x(yx3)+2yx=2x40.\frac{d}{dt}(x^2+y^2)=2x(-y-x^3)+2yx=-2x^4\le0.

非线性阻尼使轨道缓慢向内,因此原点实际是稳定焦点式吸引子,而不是闭轨中心。

常见处理方法:

  • 转极坐标看 r˙\dot r
  • 找守恒量或 Lyapunov 函数;
  • 利用对称性、可逆性;
  • 必要时使用中心流形或正规形理论。

极坐标公式

r2=x2+y2r^2=x^2+y^2

r˙=xx˙+yy˙r,θ˙=xy˙yx˙r2.\dot r=\frac{x\dot x+y\dot y}{r}, \qquad \dot\theta=\frac{x\dot y-y\dot x}{r^2}.

极坐标能把“向内/向外”和“旋转”分开,尤其适合中心、焦点、极限环和 Hopf 分岔。

局部结论与全局相图

线性化只描述不动点附近。要补全全局相图,还应检查:

  • 零增长线与各区域箭头;
  • 不变直线、对称轴;
  • 是否有轨道逃向无穷远;
  • 是否存在闭轨或捕获域;
  • 鞍点稳定/非稳定流形可能连接到哪里。

数值模拟的正确位置

数值相图适合检查草图、发现候选结构,但不能独立证明:

  • 看似闭合的轨道可能只是缓慢螺旋;
  • 有限时间图看不出极慢的不稳定;
  • 步长过大会制造假极限环。

最稳妥的结构是“解析判据给结论,数值图作验证”。

考试要点

  • 每个不动点都单独计算 JJ
  • 双曲点才可直接相信线性分类。
  • 纯虚、零特征值、重根边界都要注明线性化失效。
  • 零增长线上只有一个分量为零,不代表整条线都是不动点。

评论