強制振動（ばね‐質点系） ： シミュレーション

振動を    初期位置に

外力： ${\displaystyle f(t)=B\sin \omega t}$       ${\displaystyle B}$〔m〕：振幅， ${\displaystyle \omega }$〔rad/s〕：角振動数

時刻 ${\displaystyle t}$〔s〕における物体の位置 ${\displaystyle x}$〔m〕，速度 ${\displaystyle v}$〔m/s〕

${\displaystyle x=A\cos ({\omega }_{0}t+\alpha )+\frac{B}{m({\omega }_0^2-{\omega }^2)}\sin \omega t}$

${\displaystyle v=\frac{dx}{dt}=-A{\omega }_0\sin ({\omega }_{0}t+\alpha )+\frac{B\omega }{m({\omega }_0^2-{\omega }^2)}\cos \omega t}$

${\displaystyle a=\frac{dv}{dt}=-A{\omega }_0^2\cos ({\omega }_{0}t+\alpha )+\frac{B{\omega }^2}{m({\omega }_0^2-{\omega }^2)}\sin \omega t}$

${\displaystyle {\omega }_0=\sqrt{\frac{k}{m}}}$      ${\displaystyle k}$〔N/m〕：ばね定数， ${\displaystyle m}$〔kg〕：物体の質量