沿圆周前进的两个人,即使始终在同一条圆周上,也可能相差半圈。周期运动的稳定性因此要先说明比较对象:是某一时刻的同相位点,还是整条闭合轨道?轨道稳定性采用后者,再用横截面检测偏离是否逐圈缩小。
形式陈述
设 f ∈ C 1 ( U , R n ) ,U ⊂ R n 开,ϕ 是它的局部流 理路 向量场的流 Flow of a vector field · Local flow 把光滑向量场积分为满足局部群律的一参数局部微分同胚族。 ,Γ = { ϕ ( s , p ) : 0 ≤ s < T } 是非平凡周期轨道。用欧氏距离 理路 度量空间 Metric space 用满足正定性、对称性与三角不等式的实值距离刻画点间远近的空间。 定义
d ( x , Γ ) = min z ∈ Γ ‖ x − z ‖ . 最小值存在,因为 Γ 紧。轨道称为稳定,若任意 ε > 0 都有 δ > 0 ,使 d ( x , Γ ) < δ 时解对所有 t ≥ 0 存在,且始终满足 d ( ϕ ( t , x ) , Γ ) < ε 。若还在某个邻域内有 d ( ϕ ( t , x ) , Γ ) → 0 ,称为渐近轨道稳定。
若存在 C ≥ 1 , α > 0 , δ > 0 ,使
(1) d ( ϕ ( t , x ) , Γ ) ≤ C e − α t d ( x , Γ ) , t ≥ 0 , d ( x , Γ ) < δ , 则称为局部指数轨道稳定。这里允许相位变化,右侧并非 ‖ ϕ ( t , x ) − ϕ ( t , p ) ‖ 。
横向收缩判据。 在 p 处取一个小横截面,记Poincaré回归映射 理路 Poincaré回归映射与返回时间导数 Poincaré return map · Poincare return map · 庞加莱回归映射 以横截条件构造周期轨道附近的局部返回时间,推导返回映射的投影导数,并证明横向乘子与单周期变分矩阵的关系。 为 P ,其导数在截面切空间上为 C 0 = D P ( p ) 。若 C 0 的所有特征值都严格位于单位圆内,则式(1)成立。等价地,单周期变分矩阵除去一个相位乘子1后,剩余乘子的模都小于1。
直觉
横截面像每圈经过一次的检测门。每次过门更接近周期轨道,还不足以直接写下连续时间结论:两次过门之间不能突然远离,过门的间隔也不能无限变长。证明分别控制这三件事。
找到真正收缩的范数
令 ρ ( C 0 ) < q < 1 。由Jordan形的矩阵幂估计 理路 Jordan 标准形 Jordan canonical form · Jordan normal form 在特征多项式分裂时,把有限维算子表示成 Jordan 块直和。 ,可先取 ρ ( C 0 ) < r < q ,得到 ‖ C 0 k ‖ ≤ K r k ;Jordan块产生的有限次幂因子被 r k 吸收。于是
‖ v ‖ ∗ := ∑ k = 0 ∞ q − k ‖ C 0 k v ‖ 收敛,且与原范数等价。移一项可得 ‖ C 0 v ‖ ∗ ≤ q ‖ v ‖ ∗ 。连续可微性使某个小凸坐标球内 ‖ D P ‖ ∗ ≤ κ < 1 。以不动点为中心,这个球被映回自身,并有
(2) ‖ P k ( y ) − p ‖ ∗ ≤ κ k ‖ y − p ‖ ∗ . 只看欧氏矩阵范数可能错失这个收缩;谱半径小于1不要求原坐标下的一步算子范数小于1。
从过门次数回到真实时间
缩小截面后,返回时间 τ ( y ) 连续且位于 [ T / 2 , 3 T / 2 ] 。轨道段 ϕ ( s , y ) 在 0 ≤ s ≤ 3 T / 2 上对 y 有共同Lipschitz界。因此从第 k 次过门到下一次过门,距 Γ 至多为式(2)右侧的一个固定倍数。
还须覆盖不是从门上出发的点。取截面局部坐标 v ↦ y ( v ) ,映射 ( s , v ) ↦ ϕ ( s , y ( v ) ) 在 ( 0 , 0 ) 的导数由流方向与截面方向组成,满秩。反函数定理 理路 逆函数定理 Inverse function theorem 导数可逆的光滑映射在该点邻域内存在光滑局部逆。 给局部流坐标。沿紧轨道取有限个这样的坐标盒,其映射和逆映射都有共同Lipschitz界;其中 v = 0 正好表示轨道。因此横向坐标大小与到轨道的距离可用统一常数互相控制,且到第一个门的时间有统一上界。
这些有限段运动连同式(2)始终留在定义域内的紧邻域,所以解能够不断延续。扣除第一次到门的有界时间后,经过 t 时间至少完成 2 t / ( 3 T ) 减去一个固定常数次返回。于是可取某个 C 和 α = − 2 log κ / ( 3 T ) > 0 得到式(1)。常数 C 吸收最初的时间段与坐标换算。
例子与边界
一个不需要隐藏常数的吸引圆
取
z ′ = ( 1 − | z | 2 ) z + 2 J z , J = ( 0 − 1 1 0 ) . 极坐标给 r ′ = r ( 1 − r 2 ) 、θ ′ = 2 。半径1的周期轨道周期为 π ,其横向乘子为 e − 2 π 。若 1 / 2 ≤ r ( 0 ) ≤ 3 / 2 ,此环带正向不变,并且
( r − 1 ) ′ = − r ( r + 1 ) ( r − 1 ) , r ( r + 1 ) ≥ 3 4 . 所以直接得到全时域证书
(3) d ( z ( t ) , Γ ) = | r ( t ) − 1 | ≤ e − 3 t / 4 | r ( 0 ) − 1 | . 这个环带上的统一界比局部线性衰减率2保守,但无需把线性近似当成精确解。
横向乘子1时,非线性项仍能决定成败
在 1 / 2 < r < 3 / 2 的环带考虑 θ ′ = 1 。对 r ′ = − ( r − 1 ) 3 ,记 u 0 = r ( 0 ) − 1 ,则
r ( t ) − 1 = u 0 1 + 2 u 0 2 t . 单位圆渐近稳定,收敛却只有代数速度,不满足式(1)。若改成 r ′ = + ( r − 1 ) 3 ,任何非零小径向扰动都会离开固定小环带,单位圆不稳定。两者沿圆的径向线性方程都为 u ′ = 0 ,横向乘子均为1。单位圆边界上的乘子不能替代非线性分析。
接近轨道不等于保持预先指定的相位
令 r ′ = 0 , θ ′ = r 。半径1附近的圆始终保持径向距离 | r − 1 | ,所以单位圆轨道稳定;但相位差以 ( r − 1 ) t 积累,同相位的两条轨道可逐渐走到圆的两侧。这里轨道不吸引,横向乘子仍为1。
推论与应用
平面中的一个标量检验
平面周期轨道只有一个横向乘子 ρ 。相位乘子为1,再用线性系统的Liouville行列式公式 理路 线性常微分方程组 Linear system of ordinary differential equations 形如 x′=A(t)x+b(t) 的向量值一阶线性方程组。 ,得到
(4) ρ = det M = exp ( ∫ 0 T div f ( ϕ ( t , p ) ) d t ) . 因此散度沿一圈积分为负时,轨道局部指数稳定。对上面的吸引圆,div f = 2 − 4 r 2 = − 2 ,式(4)给出 e − 2 π ,与精确回归映射一致。若积分为零,本页的严格收缩判据不作判断。
严格收缩还给相位误差的可求和性。由于 τ 为C1,| τ ( P j y ) − T | ≤ K κ j ‖ y − p ‖ ;逐圈的额外时间因此收敛到一个有限总偏移。结合两次返回之间的连续依赖,可为每个足够近的初值选一个极限相位,使它趋近对应的周期运动。这个结论使用了指数收缩;轨道稳定的定义本身并没有承诺相位锁定。
在偏置平面振子的闭轨证书 理路 Bendixson–Dulac判据与闭轨排除 Bendixson–Dulac criterion · Bendixson-Dulac theorem · Bendixson criterion · Dulac判据 用加权散度与周期边界零通量排除平面闭轨,明确单连通和严格符号条件,并在有孔环域上比较两条闭轨、证明扰动振子周期轨唯一。 中,可以先由紧环带证明周期存在,再以两条候选闭轨间的负散度排除多解。散度上界与角速度界还给 0 < ρ ≤ e − 15 π / 34 < 1 ,此时才使用本页判据取得局部指数轨道稳定。整个环带的omega集都趋于同一轨道,是另一份定性结论,并没有自动给出全环带共用的指数常数。
周期轨道的参数延续 理路 周期轨道的参数延续与周期敏感度 Continuation of periodic orbits · Period sensitivity · 周期轨道的相位约束 固定横截相位后,用非相位乘子不等于1证明局部参数延续;由带边框变分系统同时计算基点与周期导数,并解释退化和近退化边界。 只需横向乘子不等于1即可局部追踪轨道,允许某些乘子的模大于1。可延续性与吸引性是两项不同的判定,计算时应分别记录。
参考资料