インテグレーターブロックを知ったので単振り子をモデル化してみる。θは鉛直方向下向きを0〔rad〕、左廻りをプラスとする。θが角度、θ'が角速度、θ''が角加速度。
θ'' = -(重力加速度 / 振り子の棒の長さ) * sin(θ)
https://github.com/ti-nspire/MATLAB/blob/main/pendulum_single.slx
↓ θの初期値を0.1〔rad〕(ほぼ真下)にした場合。理科の授業で習ったとおりほぼ単振動である。
↓ 棒の長さは変えずにθの初期値をpi * 0.99〔rad〕(ほぼ真上)にした場合。
