【python】空気抵抗を考慮した運動方程式の数値解析 - Qiita

Pythonを使って自由落下及びのシュミレーションをしてみようということでまとめました。 運動方程式の差分化 運動方程式の差分化 $$ F=ma $$ 高校物理でもおなじみの運動方程式である. 大学においては,運動方程式は微分を使った式で書かれることが多く, $$ F=m\frac{d^2x}{dx^2} $$ となる.上式より,以下の式が導き出される. $$ \begin{eqnarray} \frac{dv(t)}{dt} &=& a(t) \ \frac{dx(t)}{dt} &=& v(t) \end{eqnarray} $$ よって,この連立微分方程式をプログラムに落とし込めばシュミレーションが可能となる. ここで微分の定義を思い出してほしい $$ \frac{df(x)}{dx}=\frac{f(x+\Delta x)-f(x)}{\Delta x} $$ 但し,$\Delta