[matplotlib animation] 14. 外サイクロイド(epicycloid)

はじめに
コード
解説
参考

はじめに

matplotlibのFuncAnimationを用いて、外サイクロイドのアニメーションを作成する。

コード

解説

モジュールのインポート

グリッドと軸ラベル、アスペクト比の設定

データの生成

外サイクロイドは極座標系の媒介変数でデータを生成する。
θを0〜4π,φを0〜2πとする。

外サイクロイド

中心円の半径をrc, 動円の半径をrm, 回転角を θとすると外サイクロイドは以下の式で表される。

$$x=(r_c + r_m)cos\theta -r_m cos( \frac{r_c + r_m}{r_m}\theta)\\ y=(r_c + r_m)sin\theta -r_m sin( \frac{r_c + r_m}{r_m}\theta)$$

外サイクロイドを表示するためのからのリスト、epicycloids =[]をつくる。
pointsは円周上の定点を表示するためのからのリスト。

動く円

回転する円をプロットするため、円を極座標表示のx,yで設定する。円をアニメーション表示するためにからのリスト、round_circles =[]をつくる。

円を移動させるための変数の設定。動円は半径=(rm+rc)の円運動をするので上記のようになる。

中心円

中心円はアニメーションで動かないのでここでplotしておく。

アニメーションの設定

アニメーションでは古いプロットを消す必要があるので、
グラフのリストにグラフが入っている場合は、pop().remove()で消す。

round_circle, = plt.plot(x_move[num]+xm, y_move[num]+ym, ‘b-‘, lw=2)
で、外側の円を円運動するアニメーションとして表示する。b-で青い線での表示となる。

epicycloid, = plt.plot(x[:num+1] ,y[:num+1], ‘g-‘, lw=2)
で、外サイクロイドの軌跡を表示する。g-で緑線での表示となる。

point, = plt.plot(x[num], y[num], ‘ko’, markersize=4)
でサイクロイドを描写する、外側の円上の定点を表示している。koで黒丸での表示となる。

それぞれ作成したプロットをからのリストに入れることでアニメーションを表示できる。

ax.set_title(theta_str + str(theta[num]/np.pi)[:4] + str(r’ $\pi$’))でタイトルにアニメーションで変化する回転角を表示する。πを表示するには、r’ $\pi$’とする。

アニメーションの表示

FuncAnimationでanimationの表示。updateを200step実行してアニメーションとする。intervalは100 ms なので20秒のアニメーションとなる。
HTML(ani.to_html5_video())でjupyter notebook上にアニメーションを表示できる。