盯住一条周期轨道附近的每一瞬间,容易把相位差与真正的远离混在一起。选一扇横切轨道的小门,只记录每次穿门的位置,就把连续运动变成了少一维的离散映射。但门必须横着放,返回时间也不是固定常数。
形式陈述
设 f ∈ C 1 ( U , R n ) ,U ⊂ R n 开、n ≥ 2 ,ϕ ( t , x ) 是其局部流 理路 向量场的流 Flow of a vector field · Local flow 把光滑向量场积分为满足局部群律的一参数局部微分同胚族。 。设 p 在非平凡周期轨道上,最小正周期为 T ,故 ϕ ( T , p ) = p 且 f ( p ) ≠ 0 。
取 C 1 函数 h ,在 p 附近定义小截面
(1) Σ = { x : h ( x ) = 0 } , h ( p ) = 0 , ℓ = ∇ h ( p ) , ℓ T f ( p ) ≠ 0. 最后一个条件称横截。缩小截面后,存在唯一接近 T 的 C 1 返回时间 τ ( x ) ,使
τ ( p ) = T , ϕ ( τ ( x ) , x ) ∈ Σ . 定义局部Poincaré回归映射
靠 近 (2) P ( x ) = ϕ ( τ ( x ) , x ) , x ∈ Σ 靠近 p . 它在截面内为局部微分同胚,P ( p ) = p 。它的其他不动点对应附近周期轨道,而非原流的平衡点。
记 M = D x ϕ ( T , p ) 。对截面切向量 v ∈ T p Σ = ker ℓ T ,有
(3) D τ ( p ) v = − ℓ T M v ℓ T f ( p ) , (4) D P ( p ) v = ( I − f ( p ) ℓ T ℓ T f ( p ) ) M v . 括号是沿流方向投到截面切空间的投影,不一定是正交投影。
直觉
固定时间 T 的流先把扰动向量送到 M v 。终点一般偏离小门,必须再沿流走一小段时间,把终点挪回截面;式(4)减掉的正是这条流向分量。若省略返回时间导数,就会把相位偏差算进横向变化。
横截为何给出光滑返回时间
令 H ( t , x ) = h ( ϕ ( t , x ) ) 。由解对初值的可微依赖 理路 ODE变分方程与参数敏感度 Variational equation for ODEs · Smooth dependence of ODE solutions · 常微分方程的初值导数 在共同存在窗口中证明非线性ODE对初值和参数可微,构造齐次与非齐次变分方程,并给出一阶近似的余项控制与失效边界。 ,H 为C1,且
H ( T , p ) = 0 , ∂ t H ( T , p ) = ℓ T f ( p ) ≠ 0. 隐函数定理 理路 隐函数定理 Implicit function theorem 当相关偏导块可逆时,方程组局部可把部分变量表示为其余变量的函数。 因此能把 H = 0 局部唯一解成 t = τ ( x ) 。时间邻域选在 T > 0 附近,与起点的零时间相交分开。
这个论证首先保证“接近T的那次返回”。若还要称为第一正返回,可以将截面取成足够小的局部片:横截条件排除起点附近的第二次穿越;对远离0、T的一段紧时间,最小周期性使参考轨道与p保持正距离,再靠连续依赖排除邻近轨道提前碰到小片。于是没有更早的正返回。不能把任意大的截面或多次绕回的时间也无条件叫第一返回。
求导时必须把时间变化一起带入
对恒等式 h ( ϕ ( τ ( x ) , x ) ) = 0 使用链式法则 理路 链式法则 Chain rule 复合映射的导数等于各层导数按计算顺序组成的线性映射复合。 ,在p处得到
0 = ℓ T { M v + f ( p ) D τ ( p ) v } . 解出时间导数就是式(3),再对式(2)求导即得式(4)。横截分母越小,返回时间对扰动可能越敏感;分母为零时这套计算没有保证。
相位乘子怎样从矩阵中分离
自治流恒等式给 M f ( p ) = f ( p ) 。用截面切空间加流方向分解
R n = T p Σ ⊕ span { f ( p ) } , 则 M 的矩阵呈块形
(5) ( D P ( p ) 0 ∗ 1 ) . 所以 M 的谱按代数重数恰为 D P ( p ) 的谱加一个1。若还有其他1,它们仍留在横向谱中,不能全都当作相位删除。
D P ( p ) 可逆也直接可证:若式(4)为零,则 M v = c f ( p ) 。由于 M 可逆且 M f ( p ) = f ( p ) ,有 v = c f ( p ) ;再用 v ∈ T p Σ 和横截条件得 c = 0 。有限维同维单射可逆,隐函数定理的局部逆版本便给截面内局部微分同胚。
例子与边界
固定时间流有剪切,返回映射却是恒等
在环带 1 / 2 < r < 3 / 2 上用极坐标写
r ′ = 0 , θ ′ = r . 每个半径各自绕圈,周期为 2 π / r 。截面取正横轴附近的一段,径向坐标就是r。于是
τ ( r ) = 2 π r , P ( r ) = r , τ ′ ( 1 ) = − 2 π , P ′ ( 1 ) = 1. 在基点p=(1,0)、固定时间T=2π处,以径向、切向为基,
M = ( 1 0 2 π 1 ) . 固定时间流让径向扰动额外积累切向位移;式(4)投影后才恢复截面的导数1。这个例子展示返回时间修正的真实作用,不是一个可以省略的记号。
一圈后的径向收缩
对 r ′ = r ( 1 − r 2 ) 、θ ′ = 2 ,令 q = r − 2 ,则 q ′ = − 2 q + 2 。因此在正横轴截面附近,返回时间为π,且
P ( r ) = [ 1 + ( r − 2 − 1 ) e − 2 π ] − 1 / 2 , P ′ ( 1 ) = e − 2 π < 1. 半径1是不动点。周期轨道稳定性 理路 周期轨道的轨道稳定性 Orbital stability of periodic orbits · Transverse stability of a limit cycle · 周期轨道横向稳定性 用到整条轨道的距离定义稳定性,将横截面严格收缩转成连续时间指数界,并以中性径向例子划清乘子判据的边界。 将把这份截面收缩变成整个连续轨道的稳定结论,还需控制两次穿门之间的运动。
相切时,返回可能消失或不光滑
取水平流 x ′ = 1 , y ′ = 0 ,曲线 h ( x , y ) = y − x 2 = 0 在原点与流相切。邻近水平轨道的交点时间含 ± y :y > 0 有两个交点,y < 0 没有交点,靠近零时导数无界。局部唯一性仍成立,失败的是横截条件,而不是ODE解不存在。
推论与应用
改变横截面时,可沿同一流把一个小门上的点送到另一个小门,得到局部可逆映射 H 。两套回归满足 P ~ ∘ H = H ∘ P ;在固定点求导,两个横向矩阵相似。因此局部坐标或门的位置改变,不应改变横向乘子的稳定结论。
时间周期的非自治系统也可用固定时间T的映射,但其理由不同:把时间模T加入状态空间后,截面是时间坐标为零的小门,返回时间恒为T。仍必须保证初值在整个T窗口内存在,不能把固定时间映射无条件定义在全空间。
本页的局部返回接口从一条已存在周期轨出发。若起初只有一个紧正向不变平面区域,Poincaré–Bendixson定理 理路 Poincaré–Bendixson定理与平面闭轨存在 Poincaré–Bendixson theorem · Poincare-Bendixson theorem · 庞加莱–本迪克松定理 证明平面紧正向轨道的无平衡点omega集必为一条周期轨,用流盒、横截穿越次序和Jordan分离闭合证明,并构造非径向振子的环带存在证书。 可以在区域无平衡点时先认证闭轨存在;其流盒只构造局部穿门投影,不需要预先知道一圈的返回时间。周期轨的唯一性仍要另外证明,不能从本页在一个门附近的隐函数结论扩成全区域唯一。
参数延续 理路 周期轨道的参数延续与周期敏感度 Continuation of periodic orbits · Period sensitivity · 周期轨道的相位约束 固定横截相位后,用非相位乘子不等于1证明局部参数延续;由带边框变分系统同时计算基点与周期导数,并解释退化和近退化边界。 把 P λ ( x ) − x = 0 作为有限维隐式方程。是否可解取决于 D P − I ,而不是错误地对必含相位乘子1的 M − I 直接求逆。
参考资料