【要約】単振り子の運動方程式を無限級数で解く ―― Taylor解とSymPy数値解の比較 [Qiita_Trend] | Summary by TechDistill
> Source: Qiita_Trend
Execute Primary Source
// Problem
物理学者が、振幅の大きい単振り子の挙動を正確に計算しようとする際、非線形性の問題に直面する。
- ・$\sin x \approx x$ の近似では、大きな振幅において誤差が無視できない。
- ・$\sin x$ を直接展開しても、非線形方程式が残り、解の構成が困難である。
- ・厳密解であるJacobi楕円関数は、計算が複雑になる。
// Approach
著者は、非線形項を近似せずに解を構成するため、補助変数を用いた多項式系への変換手法を提案した。
- ・$v=\dot x, s=\sin x, c=\cos x$ を導入し、方程式を1階の多項式系へ変換する。
- ・係数間の関係を、行列形式の畳み込み漸化式として定式化する。
- ・$Y_{n+1}=\frac{1}{n+1}\sum_{k=0}^{n}A_k\,Y_{n-k}$ により、係数を逐次計算する。
// Result
数値実験を通じて、提案手法が収束半径内において極めて高い精度を持つことを示した。
- ・$N=40$ のとき、$t \le \pi/2$ の範囲で誤差 $10^{-11}$ 以下を達成した。
- ・打ち切り次数 $N$ を増やすほど、収束半径内での精度が向上することを確認した。
- ・収束半径を超えると、次数を増やしても誤差が急増することを実証した。
Senior Engineer Insight
> 数値計算の精度と計算コストのトレードオフを理解する上で、非常に示唆に富む内容である。多項式系への変換による漸化式の構築は、計算の安定性と実装の簡潔さを両立している。ただし、収束半径の制約は実用上の課題となる。実システムへの投入には、解析接続や周期性の利用といった、時間軸を跨ぐための戦略的な実装が不可欠である。