{"id":7157,"date":"2023-12-20T16:54:41","date_gmt":"2023-12-20T07:54:41","guid":{"rendered":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/?page_id=7157"},"modified":"2023-12-20T16:56:35","modified_gmt":"2023-12-20T07:56:35","slug":"python-%e3%81%a7%e6%a5%95%e5%86%86%e8%bb%8c%e9%81%93%e4%b8%8a%e3%81%ae%e6%99%82%e5%88%bb%e3%81%94%e3%81%a8%e3%81%ae%e4%bd%8d%e7%bd%ae%e3%82%92%e6%95%b0%e5%80%a4%e7%9a%84%e3%81%ab%e6%b1%82%e3%82%81","status":"publish","type":"page","link":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/%e3%82%b3%e3%83%b3%e3%83%94%e3%83%a5%e3%83%bc%e3%82%bf%e6%bc%94%e7%bf%92\/%e6%a5%95%e5%86%86%e8%bb%8c%e9%81%93%e4%b8%8a%e3%81%ae%e6%99%82%e5%88%bb%e3%81%94%e3%81%a8%e3%81%ae%e4%bd%8d%e7%bd%ae%e3%82%92%e6%95%b0%e5%80%a4%e7%9a%84%e3%81%ab%e6%b1%82%e3%82%81%e3%82%8b%e6%ba%96\/python-%e3%81%a7%e6%a5%95%e5%86%86%e8%bb%8c%e9%81%93%e4%b8%8a%e3%81%ae%e6%99%82%e5%88%bb%e3%81%94%e3%81%a8%e3%81%ae%e4%bd%8d%e7%bd%ae%e3%82%92%e6%95%b0%e5%80%a4%e7%9a%84%e3%81%ab%e6%b1%82%e3%82%81\/","title":{"rendered":"Python \u3067\u6955\u5186\u8ecc\u9053\u4e0a\u306e\u6642\u523b\u3054\u3068\u306e\u4f4d\u7f6e\u3092\u6570\u5024\u7684\u306b\u6c42\u3081\u3066\u30a2\u30cb\u30e1\u3092\u3064\u304f\u308b"},"content":{"rendered":"
\n
<\/div>\n
\n
\n

\u300c\u6955\u5186\u8ecc\u9053\u4e0a\u306e\u6642\u523b\u3054\u3068\u306e\u4f4d\u7f6e\u3092\u6570\u5024\u7684\u306b\u6c42\u3081\u308b\u6e96\u5099<\/a>\u300d\u3067\u4e0b\u3054\u3057\u3089\u3048\u3057\u305f\u5f0f\u3092 Python \u3067\u6570\u5024\u7684\u306b\u89e3\u304f\u3002<\/p>\n<\/div>\n<\/div>\n<\/div>\n

<\/p>\n

\n
<\/div>\n
\n
\n

\u5fc5\u8981\u306a\u30e2\u30b8\u30e5\u30fc\u30eb\u306e import<\/h3>\n<\/div>\n<\/div>\n<\/div>\n
\n
\n
In\u00a0[1]:<\/div>\n
\n
\n
\n
import<\/span> numpy<\/span> as<\/span> np<\/span>\r\nimport<\/span> matplotlib.pyplot<\/span> as<\/span> plt<\/span>\r\nfrom<\/span> matplotlib.animation<\/span> import<\/span> FuncAnimation<\/span>\r\nfrom<\/span> scipy.integrate<\/span> import<\/span> solve_ivp<\/span>\r\n\r\n# \u4ee5\u4e0b\u306f\u30b0\u30e9\u30d5\u3092 SVG \u3067 Notebook \u306b\u30a4\u30f3\u30e9\u30a4\u30f3\u8868\u793a\u3055\u305b\u308b\u8a2d\u5b9a<\/span>\r\n%<\/span>config<\/span> InlineBackend.figure_formats = ['svg']\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n
\n
<\/div>\n
\n
\n
\u7121\u6b21\u5143\u5316\u3055\u308c\u305f\u5f0f<\/h3>\n
\u300c\u6955\u5186\u8ecc\u9053\u4e0a\u306e\u6642\u523b\u3054\u3068\u306e\u4f4d\u7f6e\u3092\u6570\u5024\u7684\u306b\u6c42\u3081\u308b\u6e96\u5099<\/a>\u300d\u3067\u4e0b\u3054\u3057\u3089\u3048\u3057\u305f\u5f0f\uff08\u7121\u6b21\u5143\u5316\u6e08\u307f\uff09\u3002<\/p>\n
scipy.integrate.solve_ivp()<\/code> \u306f $\\displaystyle \\frac{dy}{dt} = f(t, y)$ \u306e\u5f62\u306e\u5f0f\u3092\u89e3\u304f\u306e\u3067\uff0c\u5909\u6570\u3092\u63c3\u3048\u308b\u3002<\/p>\n
$ \\phi \\rightarrow y, \\ T \\rightarrow t$<\/p>\n
\\begin{eqnarray}
\n\\frac{dy}{dt} &=& f(t, y) = 2 \\pi \\frac{(1+e \\cos y)^2}{(1-e^2)^{3\/2}} \\tag{1}\\\\
\nR(y) &\\equiv& \\frac{r}{a} = \\frac{1-e^2}{1+e\\cos y} \\tag{2}\\\\
\nX(y) &\\equiv& \\frac{x}{a} = R(y) \\cos y \\tag{3}\\\\
\nY(y) &\\equiv& \\frac{y}{a} = R(y) \\sin y \\tag{4}
\n\\end{eqnarray}<\/p>\n<\/div>\n<\/div>\n<\/div>\n
\n
\n
In\u00a0[2]:<\/div>\n
\n
\n
\n
# \u5fae\u5206\u65b9\u7a0b\u5f0f\u306e\u53f3\u8fba \u72ec\u7acb\u5909\u6570\u306f t, y \u6c7a\u3081\u6253\u3061<\/span>\r\ndef<\/span> f<\/span>(<\/span>t<\/span>,<\/span> y<\/span>):<\/span>\r\n    global<\/span> e<\/span>\r\n    return<\/span> 2<\/span>*<\/span>np<\/span>.<\/span>pi<\/span>*<\/span>(<\/span>1<\/span>+<\/span>e<\/span>*<\/span>np<\/span>.<\/span>cos<\/span>(<\/span>y<\/span>))<\/span>**<\/span>2<\/span>\/<\/span>((<\/span>1<\/span>-<\/span>e<\/span>**<\/span>2<\/span>)<\/span>**<\/span>(<\/span>3<\/span>\/<\/span>2<\/span>))<\/span>\r\n\r\ndef<\/span> R<\/span>(<\/span>y<\/span>):<\/span>\r\n    global<\/span> e<\/span>\r\n    return<\/span> (<\/span>1<\/span>-<\/span>e<\/span>**<\/span>2<\/span>)<\/span>\/<\/span>(<\/span>1<\/span>+<\/span>e<\/span>*<\/span>np<\/span>.<\/span>cos<\/span>(<\/span>y<\/span>))<\/span>\r\n\r\ndef<\/span> Xdaen<\/span>(<\/span>y<\/span>):<\/span>\r\n    global<\/span> e<\/span>\r\n    return<\/span> R<\/span>(<\/span>y<\/span>)<\/span>*<\/span>np<\/span>.<\/span>cos<\/span>(<\/span>y<\/span>)<\/span>\r\n\r\ndef<\/span> Ydaen<\/span>(<\/span>y<\/span>):<\/span>\r\n    global<\/span> e<\/span>\r\n    return<\/span> R<\/span>(<\/span>y<\/span>)<\/span>*<\/span>np<\/span>.<\/span>sin<\/span>(<\/span>y<\/span>)<\/span>\r\n\r\n# \u96e2\u5fc3\u7387\u3002\u304b\u306a\u308a\u5927\u304d\u3081<\/span>\r\ne<\/span> =<\/span> 0.9<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n
\n
<\/div>\n
\n
\n
solve_ivp()<\/code> \u3067\u5e38\u5fae\u5206\u65b9\u7a0b\u5f0f\u3092\u6570\u5024\u7684\u306b\u89e3\u304f<\/h3>\n<\/div>\n<\/div>\n<\/div>\n
\n
\n
In\u00a0[3]:<\/div>\n
\n
\n
\n
# T \u306e\u521d\u671f\u5024 <\/span>\r\nt0<\/span> =<\/span> 0<\/span>\r\n# T \u306e\u7d42\u4e86\u5024<\/span>\r\nt1<\/span> =<\/span> 1<\/span>\r\nt_span<\/span> =<\/span> [<\/span>t0<\/span>,<\/span> t1<\/span>]<\/span>\r\n# phi \u306e\u521d\u671f\u5024<\/span>\r\ny_ini<\/span> =<\/span> [<\/span>0<\/span>]<\/span>\r\n\r\n# \u5206\u5272\u6570<\/span>\r\nNdiv<\/span> =<\/span> 36<\/span>\r\n\r\nt_list<\/span> =<\/span> np<\/span>.<\/span>linspace<\/span>(<\/span>t0<\/span>,<\/span> t1<\/span>,<\/span> Ndiv<\/span>+<\/span>1<\/span>)<\/span>\r\nanswer<\/span> =<\/span> solve_ivp<\/span>(<\/span>f<\/span>,<\/span> t_span<\/span>,<\/span> y_ini<\/span>,<\/span> \r\n                   t_eval<\/span> =<\/span> t_list<\/span>,<\/span> \r\n                   # \u8a08\u7b97\u7cbe\u5ea6\u306e\u8a2d\u5b9a<\/span>\r\n                   rtol<\/span> =<\/span> 1.e-13<\/span>,<\/span> atol<\/span> =<\/span> 1.e-15<\/span>)<\/span>\r\n\r\n# \u8a08\u7b97\u7d50\u679c\u306f .y \u3068\u3057\u3066\u53d6\u308a\u51fa\u3057\u307e\u3059\u3002<\/span>\r\nphii<\/span> =<\/span> answer<\/span>.<\/span>y<\/span>[<\/span>0<\/span>]<\/span>\r\n\r\nX<\/span> =<\/span> Xdaen<\/span>(<\/span>phii<\/span>)<\/span>\r\nY<\/span> =<\/span> Ydaen<\/span>(<\/span>phii<\/span>)<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n
\n
<\/div>\n
\n
\n
\u3061\u3087\u3046\u30691\u5468\u671f $T = T_{\\rm end}$ \u305f\u3066\u3070 $\\phi = 2 \\pi$ \u306b\u306a\u3063\u3066\u3044\u308b\u306f\u305a\u3060\u304b\u3089\uff0c\u6570\u5024\u8aa4\u5dee\u3092\u78ba\u8a8d\u3059\u308b\u3002<\/p>\n
$$\\phi(T_{\\rm end}) – 2 \\pi$$<\/p>\n<\/div>\n<\/div>\n<\/div>\n
\n
\n
In\u00a0[4]:<\/div>\n
\n
\n
\n
# [-1] \u306f\u6700\u5f8c\u306e\u9805<\/span>\r\n\r\nphii<\/span>[<\/span>-<\/span>1<\/span>]<\/span> -<\/span> 2<\/span>*<\/span>np<\/span>.<\/span>pi<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n
\n
\n
\n
Out[4]:<\/div>\n
\n
-7.275335889289636e-11<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n
\n
<\/div>\n
\n
\n
\u4e00\u5b9a\u6642\u9593\u9593\u9694\u3054\u3068\u306e\u4f4d\u7f6e\u3092\u30b0\u30e9\u30d5\u306b<\/h3>\n
\u4e00\u5b9a\u6642\u9593\u9593\u9694 $\\displaystyle \\Delta T = \\frac{1}{N_{\\rm div}}$ \u3054\u3068\u306e\u4f4d\u7f6e\uff1a<\/p>\n<\/div>\n<\/div>\n<\/div>\n
\n
\n
In\u00a0[5]:<\/div>\n
\n
\n
\n
fig<\/span> =<\/span> plt<\/span>.<\/span>figure<\/span>(<\/span>figsize<\/span>=<\/span>[<\/span>6<\/span>,<\/span>6<\/span>])<\/span>\r\n# \u7e26\u6a2a\u306e\u30a2\u30b9\u30da\u30af\u30c8\u6bd4\u3092\u7b49\u3057\u304f<\/span>\r\nplt<\/span>.<\/span>axes<\/span>()<\/span>.<\/span>set_aspect<\/span>(<\/span>'equal'<\/span>)<\/span>\r\n# \u8ef8\u306e\u8a2d\u5b9a<\/span>\r\nplt<\/span>.<\/span>axis<\/span>([<\/span>-<\/span>2<\/span>,<\/span> 1<\/span>,<\/span> -<\/span>1<\/span>,<\/span> 1<\/span>])<\/span>\r\n\r\nplt<\/span>.<\/span>scatter<\/span>(<\/span>X<\/span>,<\/span> Y<\/span>,<\/span> c<\/span>=<\/span>\"blue\"<\/span>,<\/span> s<\/span>=<\/span>20<\/span>);<\/span>\r\n\r\n# x\u8ef8 y\u8ef8<\/span>\r\nplt<\/span>.<\/span>axhline<\/span>(<\/span>0<\/span>,<\/span> color<\/span>=<\/span>'black'<\/span>,<\/span> dashes<\/span>=<\/span>(<\/span>3<\/span>,<\/span> 3<\/span>),<\/span> linewidth<\/span>=<\/span>0.6<\/span>)<\/span>\r\nplt<\/span>.<\/span>axvline<\/span>(<\/span>0<\/span>,<\/span> color<\/span>=<\/span>'black'<\/span>,<\/span> dashes<\/span>=<\/span>(<\/span>3<\/span>,<\/span> 3<\/span>),<\/span> linewidth<\/span>=<\/span>0.6<\/span>);<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n
\n
\n
\n
<\/div>\n
\n
<\/p>\n
\"\"<\/p>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n
\n
<\/div>\n
\n
\n
\u554f\u984c<\/h3>\n
    \n
  1. \u96e2\u5fc3\u7387 $e$ \u3092\u5909\u3048\u3066\u6570\u5024\u8a08\u7b97\u3092\u884c\u3046\u3002<\/li>\n
  2. \u4e00\u5b9a\u6642\u9593\u9593\u9694\u3054\u3068\u306e\u4f4d\u7f6e\u3092\u30a2\u30cb\u30e1\u306b\u3059\u308b\u3002\n
      \n
    1. \u30b7\u30fc\u30f31\uff1a\u8ecc\u9053\u304c\u6955\u5186\u306b\u306a\u308b\u69d8\u5b50\u306e\u30a2\u30cb\u30e1\u3002<\/li>\n
    2. \u30b7\u30fc\u30f32\uff1a\u6955\u5186\u4e0a\u306b\u4e00\u5b9a\u6642\u9593\u9593\u9694\u3054\u3068\u3092\u4f4d\u7f6e\u3092\u793a\u3059\u30a2\u30cb\u30e1\u3002<\/li>\n
    3. \u30b7\u30fc\u30f33\uff1a\u4e00\u5b9a\u6642\u9593\u9593\u9694\u3054\u3068\u306b\u6383\u304f\u6247\u5f62\u3092\u5857\u308a\u3064\u3076\u3059\u30a2\u30cb\u30e1\u3002<\/li>\n<\/ol>\n<\/li>\n<\/ol>\n<\/div>\n<\/div>\n<\/div>\n
      \n
      <\/div>\n
      \n
      \n
      \u30a2\u30cb\u30e1\u30fc\u30b7\u30e7\u30f3\u4f5c\u6210<\/h3>\n
      \u6ed1\u3089\u304b\u306b\u3059\u308b\u305f\u3081 Ndiv \u306f\u5927\u304d\u3081\u306b\u3068\u3063\u3066\u6570\u5024\u7a4d\u5206<\/h4>\n<\/div>\n<\/div>\n<\/div>\n
      \n
      \n
      In\u00a0[6]:<\/div>\n
      \n
      \n
      \n
      # Ndiv = 720 \u3068\u3057\u3066\u6570\u5024\u8a08\u7b97<\/span>\r\n\r\n# T \u306e\u521d\u671f\u5024 <\/span>\r\nt0<\/span> =<\/span> 0<\/span>\r\n# T \u306e\u7d42\u4e86\u5024<\/span>\r\nt1<\/span> =<\/span> 1<\/span>\r\nt_span<\/span> =<\/span> [<\/span>t0<\/span>,<\/span> t1<\/span>]<\/span>\r\n# phi \u306e\u521d\u671f\u5024<\/span>\r\ny_ini<\/span> =<\/span> [<\/span>0<\/span>]<\/span>\r\n\r\n# \u5206\u5272\u6570<\/span>\r\nNdiv<\/span> =<\/span> 720<\/span>\r\n\r\nt_list<\/span> =<\/span> np<\/span>.<\/span>linspace<\/span>(<\/span>t0<\/span>,<\/span> t1<\/span>,<\/span> Ndiv<\/span>+<\/span>1<\/span>)<\/span>\r\nanswer<\/span> =<\/span> solve_ivp<\/span>(<\/span>f<\/span>,<\/span> t_span<\/span>,<\/span> y_ini<\/span>,<\/span> \r\n                   t_eval<\/span> =<\/span> t_list<\/span>,<\/span> \r\n                   # \u8a08\u7b97\u7cbe\u5ea6\u306e\u8a2d\u5b9a<\/span>\r\n                   rtol<\/span> =<\/span> 1.e-13<\/span>,<\/span> atol<\/span> =<\/span> 1.e-15<\/span>)<\/span>\r\n\r\n# \u8a08\u7b97\u7d50\u679c\u306f .y \u3068\u3057\u3066\u53d6\u308a\u51fa\u3057\u307e\u3059\u3002<\/span>\r\nphii<\/span> =<\/span> answer<\/span>.<\/span>y<\/span>[<\/span>0<\/span>]<\/span>\r\n\r\n# \u5f8c\u3005\u306e\u3053\u3068\u3092\u8003\u3048\u3066\u30ea\u30b9\u30c8\u306b\u5909\u63db\u3057\u3066\u304a\u304f<\/span>\r\nX<\/span> =<\/span> Xdaen<\/span>(<\/span>phii<\/span>)<\/span>.<\/span>tolist<\/span>()<\/span>\r\nY<\/span> =<\/span> Ydaen<\/span>(<\/span>phii<\/span>)<\/span>.<\/span>tolist<\/span>()<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n
      \n
      <\/div>\n
      \n
      \n
      \u307e\u305a\u306f\u9014\u4e2d\u307e\u3067\u306e\u8ecc\u9053\u3092\u63cf\u304f<\/h4>\n
      \u6570\u5024\u8a08\u7b97\u306e\u7d50\u679c\u306f X<\/code> \u3068 Y<\/code> \u306b\u683c\u7d0d\u3055\u308c\u305f\u306e\u3067\uff0c\u9069\u5b9c\uff08\u9593\u5f15\u304d\u3057\u305f\u308a\u3057\u3066\uff09\u30b0\u30e9\u30d5\u306b\u3059\u308b\u3002<\/p>\n<\/div>\n<\/div>\n<\/div>\n
      \n
      \n
      In\u00a0[7]:<\/div>\n
      \n
      \n
      \n
      # \u30b0\u30e9\u30d5\u306e\u30b5\u30a4\u30ba\u306e\u8a2d\u5b9a<\/span>\r\nfig<\/span> =<\/span> plt<\/span>.<\/span>figure<\/span>(<\/span>figsize<\/span>=<\/span>[<\/span>6<\/span>,<\/span>6<\/span>])<\/span>\r\nax<\/span> =<\/span> fig<\/span>.<\/span>add_subplot<\/span>()<\/span>\r\n\r\n# \u5916\u67a0\u3068\u76ee\u76db\u3092\u975e\u8868\u793a\u306b<\/span>\r\nax<\/span>.<\/span>axis<\/span>(<\/span>'off'<\/span>)<\/span>\r\n\r\n# \u30b0\u30e9\u30d5\u306e\u7e26\u6a2a\u306e\u30a2\u30b9\u30da\u30af\u30c8\u6bd4\u3092 equal \u306b\u3002<\/span>\r\nax<\/span>.<\/span>set_aspect<\/span>(<\/span>'equal'<\/span>)<\/span>\r\n\r\n# \u8ef8\u306e\u8a2d\u5b9a<\/span>\r\nplt<\/span>.<\/span>axis<\/span>([<\/span>-<\/span>2<\/span>,<\/span> 1<\/span>,<\/span> -<\/span>1<\/span>,<\/span> 1<\/span>])<\/span>\r\n\r\ni<\/span> =<\/span> 1<\/span>\r\n# i*20 \u756a\u76ee\u307e\u3067\u306e\u8ecc\u9053<\/span>\r\nxi<\/span> =<\/span> X<\/span>[:<\/span>i<\/span>*<\/span>20<\/span>]<\/span>\r\nyi<\/span> =<\/span> Y<\/span>[:<\/span>i<\/span>*<\/span>20<\/span>]<\/span>\r\nplt<\/span>.<\/span>plot<\/span>(<\/span>xi<\/span>,<\/span> yi<\/span>)<\/span>\r\n\r\n# i*20 \u756a\u76ee\u306e\u4f4d\u7f6e\u306b\u8d64\u4e38<\/span>\r\nplt<\/span>.<\/span>plot<\/span>(<\/span>X<\/span>[<\/span>i<\/span>*<\/span>20<\/span>],<\/span> Y<\/span>[<\/span>i<\/span>*<\/span>20<\/span>],<\/span> \"or\"<\/span>,<\/span> markersize<\/span> =<\/span> 3<\/span>)<\/span>\r\n\r\n# x\u8ef8 y\u8ef8<\/span>\r\nplt<\/span>.<\/span>axhline<\/span>(<\/span>0<\/span>,<\/span> color<\/span>=<\/span>'black'<\/span>,<\/span> dashes<\/span>=<\/span>(<\/span>3<\/span>,<\/span> 3<\/span>),<\/span> linewidth<\/span>=<\/span>0.6<\/span>)<\/span>\r\nplt<\/span>.<\/span>axvline<\/span>(<\/span>0<\/span>,<\/span> color<\/span>=<\/span>'black'<\/span>,<\/span> dashes<\/span>=<\/span>(<\/span>3<\/span>,<\/span> 3<\/span>),<\/span> linewidth<\/span>=<\/span>0.6<\/span>);<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n
      \n
      \n
      \n
      <\/div>\n
      \n
      <\/p>\n
      \"\"<\/p>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n
      \n
      <\/div>\n
      \n
      \n
      \u30b7\u30fc\u30f31\uff1aFuncAnimation()<\/code> \u3067\u30a2\u30cb\u30e1\u30fc\u30b7\u30e7\u30f3\u4f5c\u6210<\/h5>\n
      \u7279\u5b9a\u306e i<\/code> \u306e\u5024\u306e\u3068\u3053\u308d\u306e\u30b0\u30e9\u30d5\u304c\u3067\u304d\u305f\u3089\uff0c\u5f53\u8a72\u30b3\u30fc\u30c9\u3092 def func(i):<\/code> \u306e\u4e2d\u306b\u30b3\u30d4\u307a\u3057\u3066 FuncAnimation()<\/code> \u3092\u547c\u3079\u3070\uff0c\u30d1\u30e9\u30d1\u30e9\u30a2\u30cb\u30e1\u3092\u3064\u304f\u308b\u3053\u3068\u304c\u3067\u304d\u308b\u3093\u3060\u3063\u305f\u3088\u306d\u3002<\/p>\n<\/div>\n<\/div>\n<\/div>\n
      \n
      \n
      In\u00a0[8]:<\/div>\n
      \n
      \n
      \n
      # \u30b0\u30e9\u30d5\u306e\u30b5\u30a4\u30ba\u306e\u8a2d\u5b9a<\/span>\r\nfig<\/span> =<\/span> plt<\/span>.<\/span>figure<\/span>(<\/span>figsize<\/span>=<\/span>[<\/span>6<\/span>,<\/span>6<\/span>])<\/span>\r\nax<\/span> =<\/span> fig<\/span>.<\/span>add_subplot<\/span>()<\/span>\r\n\r\n# \u30b0\u30e9\u30d5\u306e\u7e26\u6a2a\u306e\u30a2\u30b9\u30da\u30af\u30c8\u6bd4\u3092 equal \u306b\u3002<\/span>\r\nax<\/span>.<\/span>set_aspect<\/span>(<\/span>'equal'<\/span>)<\/span>\r\n\r\ndef<\/span> func<\/span>(<\/span>i<\/span>):<\/span>\r\n    # \u524d\u306e frame \u3092\u6d88\u3059<\/span>\r\n    plt<\/span>.<\/span>cla<\/span>()<\/span>\r\n\r\n    # \u5916\u67a0\u3068\u76ee\u76db\u3092\u975e\u8868\u793a\u306b<\/span>\r\n    ax<\/span>.<\/span>axis<\/span>(<\/span>'off'<\/span>)<\/span>\r\n\r\n    # \u8ef8\u306e\u8a2d\u5b9a<\/span>\r\n    plt<\/span>.<\/span>axis<\/span>([<\/span>-<\/span>2<\/span>,<\/span> 1<\/span>,<\/span> -<\/span>1<\/span>,<\/span> 1<\/span>])<\/span>\r\n    \r\n    # i*20 \u756a\u76ee\u307e\u3067\u306e\u8ecc\u9053<\/span>\r\n    xi<\/span> =<\/span> X<\/span>[:<\/span>i<\/span>*<\/span>20<\/span>]<\/span>\r\n    yi<\/span> =<\/span> Y<\/span>[:<\/span>i<\/span>*<\/span>20<\/span>]<\/span>\r\n    plt<\/span>.<\/span>plot<\/span>(<\/span>xi<\/span>,<\/span> yi<\/span>,<\/span> linewidth<\/span> =<\/span> 0.5<\/span>)<\/span>\r\n\r\n    # i*20 \u756a\u76ee\u306e\u4f4d\u7f6e\u306b\u8d64\u4e38<\/span>\r\n    plt<\/span>.<\/span>plot<\/span>(<\/span>X<\/span>[<\/span>i<\/span>*<\/span>20<\/span>],<\/span> Y<\/span>[<\/span>i<\/span>*<\/span>20<\/span>],<\/span> \"or\"<\/span>,<\/span> markersize<\/span> =<\/span> 3<\/span>)<\/span>\r\n\r\n    # x\u8ef8 y\u8ef8<\/span>\r\n    plt<\/span>.<\/span>axhline<\/span>(<\/span>0<\/span>,<\/span> color<\/span>=<\/span>'black'<\/span>,<\/span> dashes<\/span>=<\/span>(<\/span>3<\/span>,<\/span> 3<\/span>),<\/span> linewidth<\/span>=<\/span>0.6<\/span>)<\/span>\r\n    plt<\/span>.<\/span>axvline<\/span>(<\/span>0<\/span>,<\/span> color<\/span>=<\/span>'black'<\/span>,<\/span> dashes<\/span>=<\/span>(<\/span>3<\/span>,<\/span> 3<\/span>),<\/span> linewidth<\/span>=<\/span>0.6<\/span>);<\/span>\r\n\r\n    plt<\/span>.<\/span>title<\/span>(<\/span>\"\u30b1\u30d7\u30e9\u30fc\u306e\u7b2c1\u6cd5\u5247\uff1a\u6955\u5186\u8ecc\u9053\"<\/span>)<\/span>\r\n\r\nframes<\/span> =<\/span> 36<\/span>\r\nani<\/span> =<\/span> FuncAnimation<\/span>(<\/span>fig<\/span>,<\/span> func<\/span>,<\/span>  \r\n        interval<\/span> =<\/span> 100<\/span>,<\/span> \r\n        frames<\/span> =<\/span> range<\/span>(<\/span>frames<\/span>))<\/span>\r\n\r\nani<\/span>.<\/span>save<\/span>(<\/span>\"Pd2anim01.mp4\"<\/span>,<\/span> dpi<\/span>=<\/span>300<\/span>)<\/span>\r\nani<\/span>.<\/span>save<\/span>(<\/span>\"Pd2anim01.gif\"<\/span>,<\/span> dpi<\/span>=<\/span>150<\/span>)<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n
      \n
      \n
      \n
      <\/div>\n
      \n
      <\/p>\n
      \n