振り子のシミュレーション

質問

剛体振り子のシミュレーションをしたのですが、減衰係数は無し、初期角度は10°、力は重力のみです。結果、周期は一定なのですが振幅角は減少していきます。つまり、エネルギー保存で、同じ高さまであがるはずなのにそうなりません。これはなぜなのでしょうか？よろしくお願いします。

回答No.1

ojisan7 · ベストアンサー

方程式は正しいですか？また、数値計算では、誤差の評価もしなければなりませんが、どうでしょうか。

プログラムのミスでなければ、数値計算で反復計算をするなかで「桁落ち」の累積が原因であることも考えられます。「桁落ち」によって、期待通りに動作しない（誤動作）ことはよくあることです。プログラムはそんなに長くないと思いますので、もう一度見直しをしてください。

お礼

速い解答、大変ありがとうございました。
「なるほど」と自分でつぶやいてしまいましたｗ。
桁落ちが原因とは思いもよらずお恥ずかしい限りです。
ありがとうございました。