Skip to content

定理Theorem

周期轨道的轨道稳定性

Orbital stability of periodic orbits · Transverse stability of a limit cycle · 周期轨道横向稳定性

用到整条轨道的距离定义稳定性,将横截面严格收缩转成连续时间指数界,并以中性径向例子划清乘子判据的边界。

沿圆周前进的两个人,即使始终在同一条圆周上,也可能相差半圈。周期运动的稳定性因此要先说明比较对象:是某一时刻的同相位点,还是整条闭合轨道?轨道稳定性采用后者,再用横截面检测偏离是否逐圈缩小。

形式陈述 ​

设 f∈C1(U,Rn),U⊂Rn 开,ϕ 是它的局部流,Γ={ϕ(s,p):0≤s<T} 是非平凡周期轨道。用欧氏距离定义

d(x,Γ)=minz∈Γ‖x−z‖.

最小值存在,因为 Γ 紧。轨道称为稳定,若任意 ε>0 都有 δ>0,使 d(x,Γ)<δ 时解对所有 t≥0 存在,且始终满足 d(ϕ(t,x),Γ)<ε。若还在某个邻域内有 d(ϕ(t,x),Γ)→0,称为渐近轨道稳定。

若存在 C≥1,α>0,δ>0,使

(1)d(ϕ(t,x),Γ)≤Ce−αtd(x,Γ),t≥0,d(x,Γ)<δ,

则称为局部指数轨道稳定。这里允许相位变化,右侧并非 ‖ϕ(t,x)−ϕ(t,p)‖。

横向收缩判据。 在 p 处取一个小横截面,记Poincaré回归映射为 P,其导数在截面切空间上为 C0=DP(p)。若 C0 的所有特征值都严格位于单位圆内,则式(1)成立。等价地,单周期变分矩阵除去一个相位乘子1后,剩余乘子的模都小于1。

直觉

横截面像每圈经过一次的检测门。每次过门更接近周期轨道,还不足以直接写下连续时间结论:两次过门之间不能突然远离,过门的间隔也不能无限变长。证明分别控制这三件事。

找到真正收缩的范数 ​

令 ρ(C0)<q<1。由Jordan形的矩阵幂估计,可先取 ρ(C0)<r<q,得到 ‖C0k‖≤Krk;Jordan块产生的有限次幂因子被 rk 吸收。于是

‖v‖∗:=∑k=0∞q−k‖C0kv‖

收敛,且与原范数等价。移一项可得 ‖C0v‖∗≤q‖v‖∗。连续可微性使某个小凸坐标球内 ‖DP‖∗≤κ<1。以不动点为中心,这个球被映回自身,并有

(2)‖Pk(y)−p‖∗≤κk‖y−p‖∗.

只看欧氏矩阵范数可能错失这个收缩;谱半径小于1不要求原坐标下的一步算子范数小于1。

从过门次数回到真实时间 ​

缩小截面后,返回时间 τ(y) 连续且位于 [T/2,3T/2]。轨道段 ϕ(s,y) 在 0≤s≤3T/2 上对 y 有共同Lipschitz界。因此从第 k 次过门到下一次过门,距 Γ 至多为式(2)右侧的一个固定倍数。

还须覆盖不是从门上出发的点。取截面局部坐标 v↦y(v),映射 (s,v)↦ϕ(s,y(v)) 在 (0,0) 的导数由流方向与截面方向组成,满秩。反函数定理给局部流坐标。沿紧轨道取有限个这样的坐标盒,其映射和逆映射都有共同Lipschitz界;其中 v=0 正好表示轨道。因此横向坐标大小与到轨道的距离可用统一常数互相控制,且到第一个门的时间有统一上界。

这些有限段运动连同式(2)始终留在定义域内的紧邻域,所以解能够不断延续。扣除第一次到门的有界时间后,经过 t 时间至少完成 2t/(3T) 减去一个固定常数次返回。于是可取某个 C 和 α=−2log⁡κ/(3T)>0 得到式(1)。常数 C 吸收最初的时间段与坐标换算。

例子与边界

一个不需要隐藏常数的吸引圆 ​

取

z′=(1−|z|2)z+2Jz,J=(0−110).

极坐标给 r′=r(1−r2)、θ′=2。半径1的周期轨道周期为 π,其横向乘子为 e−2π。若 1/2≤r(0)≤3/2,此环带正向不变,并且

(r−1)′=−r(r+1)(r−1),r(r+1)≥34.

所以直接得到全时域证书

(3)d(z(t),Γ)=|r(t)−1|≤e−3t/4|r(0)−1|.

这个环带上的统一界比局部线性衰减率2保守,但无需把线性近似当成精确解。

横向乘子1时,非线性项仍能决定成败 ​

在 1/2<r<3/2 的环带考虑 θ′=1。对 r′=−(r−1)3,记 u0=r(0)−1,则

r(t)−1=u01+2u02t.

单位圆渐近稳定,收敛却只有代数速度,不满足式(1)。若改成 r′=+(r−1)3,任何非零小径向扰动都会离开固定小环带,单位圆不稳定。两者沿圆的径向线性方程都为 u′=0,横向乘子均为1。单位圆边界上的乘子不能替代非线性分析。

接近轨道不等于保持预先指定的相位 ​

令 r′=0,θ′=r。半径1附近的圆始终保持径向距离 |r−1|,所以单位圆轨道稳定;但相位差以 (r−1)t 积累,同相位的两条轨道可逐渐走到圆的两侧。这里轨道不吸引,横向乘子仍为1。

推论与应用

平面中的一个标量检验 ​

平面周期轨道只有一个横向乘子 ρ。相位乘子为1,再用线性系统的Liouville行列式公式,得到

(4)ρ=det⁡M=exp(∫0Tdivf(ϕ(t,p))dt).

因此散度沿一圈积分为负时,轨道局部指数稳定。对上面的吸引圆,divf=2−4r2=−2,式(4)给出 e−2π,与精确回归映射一致。若积分为零,本页的严格收缩判据不作判断。

严格收缩还给相位误差的可求和性。由于 τ 为C1,|τ(Pjy)−T|≤Kκj‖y−p‖;逐圈的额外时间因此收敛到一个有限总偏移。结合两次返回之间的连续依赖,可为每个足够近的初值选一个极限相位,使它趋近对应的周期运动。这个结论使用了指数收缩;轨道稳定的定义本身并没有承诺相位锁定。

在偏置平面振子的闭轨证书中,可以先由紧环带证明周期存在,再以两条候选闭轨间的负散度排除多解。散度上界与角速度界还给 0<ρ≤e−15π/34<1,此时才使用本页判据取得局部指数轨道稳定。整个环带的omega集都趋于同一轨道,是另一份定性结论,并没有自动给出全环带共用的指数常数。

周期轨道的参数延续只需横向乘子不等于1即可局部追踪轨道,允许某些乘子的模大于1。可延续性与吸引性是两项不同的判定,计算时应分别记录。

参考资料
关系图谱16 个相邻概念 · 2 类关系

拖动节点调整位置。

显示关系

显示:依赖

  1. 前置三跳
  2. 前置二跳
  3. 前置一跳
  4. 当前条目
  5. 后续一跳
  6. 后续二跳
  7. 后续三跳
文字版关系按与当前条目的最短距离分组
类型化关系