不变性, 吸引子/域和极坐标

dsdy Math 发布于 2025-12-27 67 次阅读


不变性(invariance set)

不变性分为不变集(invariant set)和对称不变性(symmetry)

不变域

如果(x(t0),y(t0))(x(t_0),y(t_0))在不变域SS中, 那么对于任意的 tt, 都有(x(t),y(t))S(x(t),y(t))\in S. 而不变域实际上是轨迹的集合. 对于特征值为纯虚数的平衡点centre, 每一个闭合的轨迹都是不变域.

一般来说, 不变域的描述是S={(x,y)|G(x,y)=0}S=\{(x,y)|G(x,y)=0\}. 在这个域上, 函数GG为一个常值, 即他对于tt的导数值为0. 就是 dGdt=Gxx˙+Gyy˙=0\cfrac{dG}{dt}=\cfrac{\partial G}{\partial x}\dot{x}+\cfrac{\partial G}{\partial y}\dot{y}=0.

对称不变性

对于方程x˙=f(x,y),y˙=g(x,y)\dot{x}=f(x,y),\dot{y}=g(x,y). 如果有映射x~=a(x,y,t),y~=b(x,y,t),t~=c(x,y,t)\tilde{x}=a(x,y,t),\tilde{y}=b(x,y,t),\tilde{t}=c(x,y,t)能使得x~˙=f(x~,y~),y~˙=g(x~,y~)\dot{\tilde{x}}=f(\tilde{x},\tilde{y}),\dot{\tilde{y}}=g(\tilde{x},\tilde{y}), 那么就称xa(x,y,t),yb(x,y,t),tc(x,y,t)x\to a(x,y,t),y\to b(x,y,t),t\to c(x,y,t)为系统的一个symmetry.

例如, 对于系统x˙=siny,y˙=sinx\dot{x}=\sin y,\dot{y}=\sin x, 不难发现yy,tty\to -y,t\to-txy,yxx\to y,y\to x是两对symmetry. 同时, 由于方程组的对称性, 还存在一个不变域E=cosycosxE=\cos y-\cos x, 其导数为y˙siny+x˙sinx=y˙x˙+x˙y˙=0-\dot{y}sin y+\dot{x}\sin x=-\dot{y}\dot{x}+\dot{x}\dot{y}=0.

同时, 一个具有对称不变性的系统可能会存在非线性中心(nonlinear centres)和异宿轨道(heteroclinic trajectories), 具体图像如下

吸引子(attractor) 和 吸引域(Basin of attraction)

吸引子是是动力系统中所有邻近轨道都趋向收敛的不变集合. 且这个不变域要足够小, 使得其中有且只能有一个吸引子.

每个吸引域对应一个吸引子. 任何以吸引域内的点作为初始值的(x,y)(x,y)都最终会在tt\to\infty时趋于该吸引子. 如图所示

极坐标

对于方程组x˙=f(x,y),y˙=g(x,y)\dot{x}=f(x,y),\dot{y}=g(x,y), 假设x(t)=r(t)cosθ(t),y(t)=r(t)sinθ(t)x(t)=r(t)\cos\theta(t),y(t)=r(t)\sin\theta(t). 就有r2(t)=x2(t)+y2(t)r^2(t)=x^2(t)+y^2(t). 于是可以得出r˙=xx˙+yy˙r\dot{r}=\cfrac{x\dot{x}+y\dot{y}}{r}. 对于x(t)=r(t)cosθ(t)x(t)=r(t)\cos\theta(t), 有x˙=r˙cosθrθ˙sinθ=xx˙+yy˙rxrθ˙y\dot{x}=\dot{r}\cos\theta-r\dot{\theta}\sin\theta=\cfrac{x\dot{x}+y\dot{y}}{r}\cdot\cfrac{x}{r}-\dot{\theta}y. 于是r2x˙=(xx˙+yy˙)xr2θ˙yr^2\dot{x}=(x\dot{x}+y\dot{y})x-r^2\dot{\theta}y. 由于r2(t)=x2(t)+y2(t)r^2(t)=x^2(t)+y^2(t), 所以(x2+y2)x˙=(xx˙+yy˙)xr2θ˙y(x^2+y^2)\dot{x}=(x\dot{x}+y\dot{y})x-r^2\dot{\theta}y. 仅需化简一下, 就有r2θ˙y=x2x˙+xyy˙x2x˙y2x˙=y(xy˙yx˙)r^2\dot{\theta}y=x^2\dot{x}+xy\dot{y}-x^2\dot{x}-y^2\dot{x}=y(x\dot{y}-y\dot{x}). 于是θ˙=xy˙yx˙r2\dot{\theta}=\cfrac{x\dot{y}-y\dot{x}}{r^2}.

仅需记住r˙=xx˙+yy˙r\dot{r}=\cfrac{x\dot{x}+y\dot{y}}{r}θ˙=xy˙yx˙r2\dot{\theta}=\cfrac{x\dot{y}-y\dot{x}}{r^2}即可. 如果r˙>0\dot{r}>0, 那么节点不稳定, 反之则稳定. 如果θ˙>0\dot{\theta}>0, 那么轨迹的方向就是逆时针(角度增大).

那么为什么要引入极坐标呢?

由于在非线性方程组的局部线性是由Hartman-Grobman理论给出. 其局限性在于实部为0时失效. 所以对于求出的平衡点为linear centre和non-isolated equilibrium point时, 我们需要借助极坐标来进一步判断.

注意此时我们仅对linear centre进行极坐标讨论. 其有可能变成spiral. 而对于non-isolated equilibrium point, 未看见方法.

此作者没有提供个人介绍。
最后更新于 2025-12-27