[论文解读] Simple pendulum dynamics: revisiting the Fourier-based approach to the solution

Riccardo Borghi|arXiv (Cornell University)|Feb 14, 2013

Control and Stability of Dynamical Systems被引用 6

一句话总结

本文通過引入一種迭代方法，重新探討基於傅里葉的簡諧擺非線性動力學解法，將擺的運動表示為由貝塞爾函數導出係數的傅里葉級數，提供了一種比橢圓函數更易理解的替代方案。該方法在振幅達 π/2 時仍能保持高精度的解析近似，最低階近似以 'besinc' 函數描述週期，誤差低於 0.2%。

ABSTRACT

The Fourier-based analysis customarily employed to analyze the dynamics of a simple pendulum is here revisited to propose an elementary iterative scheme aimed at generating a sequence of analytical approximants of the exact law of motion. Each approximant is expressed by a Fourier sum whose coefficients are given by suitable linear combinations of Bessel functions, which are expected to be more accessible, especially at an undergraduate level, with respect to Jacobian elliptic functions. The first three approximants are explicitely obtained and compared with the exact solution for typical initial angular positions of the pendulum. In particular, it is shown that, at the lowest approximation level, the law of motion of the pendulum turns out to be adequately described, up to oscillation amplitudes of $π/2$, by a sinusoidal temporal behaviour with a frequency proportional to the square root of the so-called "besinc" function, well known in physical optics.

研究动机与目标

• 開發一種更易理解的解析方法來處理簡諧擺的非線性動力學，避免依賴雅可比橢圓函數。
• 重新表述基於傅里葉的解法，使用更適合大學本科水平的貝塞爾函數。
• 提供一種迭代方案，生成越來越精確的擺的運動與週期的解析近似。
• 證明最低階近似能精確描述振幅達 π/2 時的擺的行為。
• 強調物理學中數學統一性，通過將擺動動力學與物理光學中熟知的 besinc 函數聯繫起來。

提出的方法

• 使用迭代傅里葉級數展開，將擺的角位置表示為基頻 ω 的奇次諧波之和。
• 通過涉及第一類貝塞爾函數的正交展開推導傅里葉係數，取代先前工作中使用的多項式近似。
• 在每一級迭代中，通過將當前近似代入非線性微分方程並利用傅里葉分析匹配係數來更新解。
• 透過強制滿足初始條件並使所得級數與運動方程匹配，自洽地確定頻率 ω。
• 使用貝塞爾函數的已知積分表達式（DLMF 公式 10.12.2 和 10.12.3），以貝塞爾函數表示 sin(ϑ(t)) 的傅里葉係數。
• 該過程迭代應用，以獲得逐級更高階的近似，並明確推導前三階近似，與精確解進行比較。

实验结果

研究问题

• RQ1 能否使用貝塞爾函數的傅里葉迭代方法，提供一種比橢圓函數更易理解的簡諧擺非線性方程求解替代方案？
• RQ2 最低階傅里葉近似在多大程度上能精確描述擺的運動與週期，特別是在大振幅情況下？
• RQ3 在振幅達 π/2 時，基於貝塞爾函數的係數與精確解在數值上相比如何？
• RQ4 能否利用擺動動力學與光學中 besinc 函數的聯繫，來增強教學理解？
• RQ5 基於貝塞爾函數的迭代近似在預測擺的週期與角位置方面有多準確？

主要发现

• 前三階近似明確使用貝塞爾函數推導，提供了擺運動越來越精確的表示。
• 在最低階近似下，擺的運動可良好地由與 besinc 函數平方根成比例的正弦函數描述。
• 對於振幅達 π/2 的情況，週期預測精度優於 0.2%，角位置預測誤差優於 0.1%。
• 該方法避開了雅可比橢圓函數，改用更易為大學本科生理解的貝塞爾函數。
• 與物理光學中熟知的 besinc 函數的聯繫，提供了跨不同物理現象的統一數學洞見。
• 該迭代方案系統且可擴展至任意高精度，每一步通過自洽的係數更新提升近似精度。