[论文解读] Pendulum Integration and Elliptic Functions

Pedro L. Garrido, Giovanni Gallavotti|arXiv (Cornell University)|Dec 12, 2008

Advanced Differential Equations and Dynamical Systems参考文献 11被引用 3

一句话总结

本文利用雅可比椭函数，严格推导了经典单摆在其不稳定平衡点附近的法向双曲坐标系，纠正了先前的错误，并建立了双曲固定点与椭圆固定点之间的直接对应关系。主要贡献在于以椭圆积分和模参数表示哈密顿量的闭式表达式，从而通过 $ g_0(x') $ 实现摆动运动的精确积分，同时给出了能量与作用量变量的显式幂级数和乘积展开式。

ABSTRACT

Revisiting canonical integration of the classical pendulum around its unstable equilibrium, normal hyperbolic canonical coordinates are constructed

研究动机与目标

• 使用椭函数理论重新推导摆在不稳定固定点附近的法向正则形式。
• 纠正早期处理中关于法向双曲坐标构造的不准确之处，特别是 [6, 附录 9] 中的问题。
• 建立不稳定平衡点附近双曲动力学与稳定平衡点附近椭圆动力学之间的精确对应关系。
• 提供以模参数表示的能量与作用量变量的显式幂级数和乘积展开式。

提出的方法

• 使用雅可比椭函数 $ \text{sn}, \text{cn}, \text{dn} $，以参数 $ u = gt/k $ 表示摆的时间演化。
• 引入模参数 $ x' = e^{-\frac{\tau}{\tau'}} $，其中 $ \tau = \frac{\text{K}(h')}{\text{K}(h)} $，以表达能量与作用量变量。
• 推导哈密顿量 $ G(pq) = \frac{d\tilde{\text{U}}}{d(pq)}(0) = g $，表明其依赖于乘积 $ pq $，从而确认法向双曲性。
• 通过包含 $ p', q' $ 的无穷级数构造生成函数 $ S_s' $ 和 $ R_s' $，其中 $ p' = \sqrt{x_s'} \cos(g_0^{(s)}t) $，$ q' = \sqrt{x_s'} \sin(g_0^{(s)}t) $。
• 应用变换 $ (p,q) = (a_s(x_s')p', a_s(x_s')q') $，其中 $ a_s^2(z) = -16I \frac{d}{dz}g_0^{(s)}(z) $，以保证正则结构。
• 利用乘积展开式 $ g_0(x') = g \prod_{n=1}^\infty \left( \frac{1 + x'^n}{1 - x'^n} \right)^2 $ 表示 $ g_0 $ 及其对数导数，与 theta 函数建立联系。

实验结果

研究问题

• RQ1 如何利用椭函数理论严格推导摆在不稳定平衡点附近的法向正则形式？
• RQ2 在 [6, 附录 9] 中所呈现的法向双曲坐标构造中，需要进行哪些修正？
• RQ3 能量 $ U $ 如何与模参数 $ x' $ 相关联？其幂级数或乘积展开的结构是什么？
• RQ4 通过解析延拓，双曲（不稳定）与椭圆（稳定）情形之间的精确变换关系是什么？
• RQ5 如何从椭函数参数化中显式构造作用量-角度变量？

主要发现

• 能量 $ U = 2g^2I / k^2 $ 以模参数 $ x' $ 表示，其中 $ g_0(x') = \frac{\pi}{2} \frac{g}{h' \text{K}(h)} $，并具有乘积展开式 $ g_0(x') = g \prod_{n=1}^\infty \left( \frac{1 + x'^n}{1 - x'^n} \right)^2 $。
• 证明能量 $ U(x') $ 满足 $ \frac{d}{dx'}U(x') = g_0(x') \frac{d}{dx'}(x' a(x')^2) $，其中 $ a(x')^2 = 8I \frac{d}{dx'}g_0(x') $，从而将能量变化与模导数联系起来。
• 对于稳定情形，能量为 $ U_s(x_s') = 32Ig_s^2 x_s' \prod_{n=1}^\infty \left( \frac{1 + x_s'^{2n}}{1 + x_s'^{2n-1}} \right)^8 $，其中 $ x_s' = -x' $，显示出与双曲情形的对偶结构。
• 稳定情形的法向哈密顿量为 $ \mathcal{U}_s(x) = 32Ig_s^2 W\left( \frac{x}{64Ig_s} \right) $，其中 $ W(z) = z(1 - 2z - 4z^2 - 20z^3 - \cdots) $，与已知的模形式一致。
• 从 $ (B,\beta) $ 到正则坐标 $ (p,q) $ 的变换通过 $ p' = \sqrt{x_s'} \cos(g_0^{(s)}t) $，$ q' = \sqrt{x_s'} \sin(g_0^{(s)}t) $ 构造，随后通过 $ a_s(x_s') $ 缩放以保证雅可比行列式为单位。
• 函数 $ g_0(x') $ 的对数导数为 $ 4 \sum_{n=1}^\infty \frac{n x'^{n-1}}{1 - x'^{2n}} $，其等于 $ \frac{1}{2} \frac{d^2}{dz^2} \log \theta_4(z, x') \big|_{z=0} $，从而与 theta 函数建立联系。