[论文解读] Mathematical pendulum and its variants

Chiş Oana, Dumitru Opriş|ArXiv.org|May 27, 2009

Quantum chaos and dynamical systems参考文献 9被引用 3

一句话总结

本文建立了一套数学框架，通过守恒律与等值面约束，将 R³ 中的欧拉陀螺系统与数学摆的运动联系起来。该框架进一步扩展至具有时滞、分数阶及随机扰动的两类系统，证明了在不变曲面上会涌现出摆动行为，并通过椭圆函数获得解析解，同时利用先进的积分格式进行数值验证。

ABSTRACT

In this paper we show that there are applications that transform the movement of a pendulum into movements in $\mathbb{R}^3$. This can be done using Euler top system of differential equations. On the constant level surfaces, Euler top system reduces to the equation of a pendulum. Those properties are also considered in the case of system of differential equations with delay argument and in the fractional case. Another aspect presented here is stochastic Euler top system of differential equations and stochastic pendulum.

研究动机与目标

• 通过哈密顿-泊松实现与守恒律，建立 R³ 中欧拉陀螺系统与数学摆之间的数学对应关系。
• 将此对应关系扩展至具有时间延迟、分数阶导数及随机扰动的系统。
• 推导在各种扩展情形下（包括椭圆与双曲函数）欧拉陀螺与摆的解析解。
• 开发并验证用于模拟分数阶与随机系统的数值方法，如 Adams-Moulton 法与 Milstein 法。
• 利用伊tô与斯特拉托诺维奇微积分分析欧拉陀螺与摆的随机动力学，包括概率密度演化对应的 Fokker-Planck 方程。

提出的方法

• 采用 R³ 中的欧拉陀螺常微分方程组：dx₁/dt = x₂x₃，dx₂/dt = -x₁x₃，dx₃/dt = x₁x₂，右侧为多项式表达式。
• 应用哈密顿-泊松实现，识别出三个守恒律（H₁, H₂, H₃）与三个卡西米尔函数（C₁, C₂, C₃），定义不变曲面。
• 将欧拉陀螺系统限制在这些守恒律的常数值曲面上，以恢复摆动动力学，包括标准、时滞、分数阶及随机变体。
• 引入沿 OZ 与 OX 轴的时滞欧拉陀螺系统，保持守恒律，并在曲面上导出带延迟的摆方程。
• 对欧拉陀螺系统应用 Caputo 分数阶导数，保持守恒律，并在曲面上约化为分数阶摆方程。
• 使用伊tô与斯特拉托诺维奇随机微分方程（SDE）建模随机欧拉陀螺与摆系统，噪声项为 √xᵢ dWᵢ，并推导相应的 Fokker-Planck 方程以描述概率密度演化。

实验结果

研究问题

• RQ1 欧拉陀螺系统在守恒律的常数值曲面上的限制，如何导出数学摆方程？
• RQ2 在能量常数相等（H = K）与不等（H ≠ K）的情形下，欧拉陀螺与摆系统的解析解是什么？
• RQ3 欧拉陀螺系统中的时滞如何影响摆动行为的出现及守恒律的保持？
• RQ4 分数阶欧拉陀螺系统是否能保持守恒律，并在不变曲面上约化为分数阶摆方程？
• RQ5 通过伊tô与斯特拉托诺维奇积分建模的随机扰动如何影响欧拉陀螺与摆的动力学？对应的 Fokker-Planck 方程是什么？

主要发现

• 当 H = K 时，系统存在精确的解析解，以双曲函数表示：x₁(t) = ±H√2 sech(H√2 t)，x₂(t) = ±H√2 tanh(H√2 t)，x₃(t) = ±H√2 sech(H√2 t)，代表异宿轨道。
• 当 H ≠ K 时，解以雅可比椭圆函数表示：x₂(t) = H√2 sn(H√2 t; √H/√K)，x₁(t) = H√2 cn(H√2 t; √H/√K)，x₃(t) = K√2 sn(H√2 t; √H/√K)，具有不同周期。
• 带有延迟参数的欧拉陀螺系统保持守恒律，其在曲面上的限制产生带延迟的数学摆，维持相同的动力学结构。
• 采用 Caputo 导数的分数阶欧拉陀螺系统保持守恒律，其在曲面上的限制产生分数阶摆方程。
• 随机欧拉陀螺与摆系统通过伊tô与斯特拉托诺维奇 SDE 建模，噪声项为 √xᵢ dWᵢ，其对应的 Fokker-Planck 方程已推导，用于描述概率密度演化。
• 采用 Milstein 格式（用于 SDE）与 Adams-Moulton 方法（用于分数阶系统）的数值模拟验证了理论动力学，所有情形下均观察到稳定收敛。