{"id":3922,"date":"2022-09-15T17:02:55","date_gmt":"2022-09-15T08:02:55","guid":{"rendered":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/?p=3922"},"modified":"2024-02-22T11:44:07","modified_gmt":"2024-02-22T02:44:07","slug":"gnuplot-%e3%81%a7%e3%82%bd%e3%83%ac%e3%83%8e%e3%82%a4%e3%83%89%e3%82%92%e8%9e%ba%e6%97%8b%e7%8a%b6%e3%81%ab%e3%81%9d%e3%82%8c%e3%82%89%e3%81%97%e3%81%8f%e6%8f%8f%e3%81%8f","status":"publish","type":"post","link":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/3922\/","title":{"rendered":"gnuplot \u3067\u30bd\u30ec\u30ce\u30a4\u30c9\u3092\u87ba\u65cb\u72b6\u306b\u305d\u308c\u3089\u3057\u304f\u63cf\u304f"},"content":{"rendered":"

\u300c\u5186\u96fb\u6d41\u56de\u8def\u306e\u91cd\u306d\u5408\u308f\u305b\u3067\u30bd\u30ec\u30ce\u30a4\u30c9\u3092\u63cf\u304f<\/a>\u300d\u3067\u306f\uff0c\u5186\u96fb\u6d41\u56de\u8def\u306e\u91cd\u306d\u5408\u308f\u305b\u3068\u3057\u3066\uff08\u624b\u3092\u629c\u3044\u3066\uff09\u63cf\u3044\u3066\u307f\u305f\u304c\uff0c\u305d\u3046\u3067\u306f\u306a\u304f\uff0c1\u672c\u306e\u96fb\u7dda\u3092\u87ba\u65cb\u72b6\u306b\u5bc6\u306b\u5dfb\u3044\u305f\u3082\u306e\u3068\u3057\u3066\u30bd\u30ec\u30ce\u30a4\u30c9\u3092\u305d\u308c\u3089\u3057\u304f\u63cf\u304f\u3002\u30bd\u30ec\u30ce\u30a4\u30c9\u306e\u5185\u90e8\u78c1\u5834\u3082\u63cf\u304f\u3002\u5b9f\u969b\u306e\u78c1\u5834\u306e\u8a08\u7b97\u306e\u969b\u306f\uff0c\u30bd\u30ec\u30ce\u30a4\u30c9\u3092\u5186\u96fb\u6d41\u56de\u8def\u306e\u91cd\u306d\u5408\u308f\u305b\u3068\u3057\u3066\u8a08\u7b97\u3059\u308b\u306e\u3067\u3059\u3088\u3002<\/p>\n

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

\u30bd\u30ec\u30ce\u30a4\u30c9\u3092\u87ba\u65cb\u72b6\u306b\u305d\u308c\u3089\u3057\u304f\u63cf\u304f<\/h3>\n

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

\n
\n
In\u00a0[1]:<\/div>\n
\n
\n
\n
%gnuplot<\/span> inline<\/span> svg<\/span> size<\/span> 640<\/span>,<\/span>640<\/span> fixed<\/span> enhanced<\/span> font<\/span> 'Noto Sans CJK JP,14'<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n
\n
\n
In\u00a0[2]:<\/div>\n
\n
\n
\n
# \u30bd\u30ec\u30ce\u30a4\u30c9<\/span>\r\nreset<\/span>\r\nunset<\/span> xtics<\/span>\r\nunset<\/span> ytics<\/span>\r\nunset<\/span> ztics<\/span>\r\nunset<\/span> border<\/span>\r\nset<\/span> xyplane<\/span> 0<\/span>\r\nset<\/span> zeroaxis<\/span>\r\nset<\/span> xrange<\/span> [<\/span>-2<\/span>:<\/span>2<\/span>]<\/span>\r\nset<\/span> yrange<\/span> [<\/span>-2<\/span>:<\/span>2<\/span>]<\/span>\r\nset<\/span> samples<\/span> 1000<\/span>\r\nset<\/span> view<\/span> 60<\/span>,<\/span>20<\/span>,,<\/span>\r\n\r\n# reds.pal \u6539<\/span>\r\nset<\/span> palette<\/span> defined<\/span> (<\/span> \\\r\n  0<\/span> '#FCBBA1'<\/span>,<\/span>\\\r\n  1<\/span> '#FCBBA1'<\/span>,<\/span>\\\r\n  2<\/span> '#FC9272'<\/span>,<\/span>\\\r\n  3<\/span> '#FB6A4A'<\/span>,<\/span>\\\r\n  4<\/span> '#EF3B2C'<\/span>,<\/span>\\\r\n  5<\/span> '#CB181D'<\/span>,<\/span>\\\r\n  6<\/span> '#99000D'<\/span>,<\/span>\\\r\n  7<\/span> '#600000'<\/span> \\\r\n)<\/span> \r\nset<\/span> key<\/span> inside<\/span> sample<\/span> 2<\/span>\r\nset<\/span> title<\/span> \"\u30bd\u30ec\u30ce\u30a4\u30c9\"<\/span>\r\nunset<\/span> colorbox<\/span>\r\n\r\na<\/span> =<\/span> 1<\/span>\r\npitch<\/span>=<\/span>0.01<\/span>\r\nN<\/span> =<\/span> 20<\/span>\r\nuini<\/span> =<\/span> -<\/span>pi<\/span>\/<\/span>2<\/span>\r\nuend<\/span> =<\/span> 2<\/span>*<\/span>pi<\/span>*<\/span>N<\/span>+<\/span>pi<\/span>\/<\/span>2<\/span>\r\n\r\n# y > 0 \u5074\u306e\u30c7\u30fc\u30bf<\/span>\r\nset<\/span> print<\/span> \"solep.txt\"<\/span>\r\n  dx<\/span> =<\/span> 0.1<\/span>\r\n  do<\/span> for<\/span> [<\/span>i<\/span>=<\/span>0<\/span>:<\/span>30<\/span>]{<\/span>\r\n    print<\/span> -3<\/span>+<\/span>dx<\/span>*<\/span>i<\/span>,<\/span> -<\/span>a<\/span>,<\/span> pitch<\/span>*<\/span>uini<\/span>\r\n  }<\/span>\r\n  Ndiv<\/span> =<\/span> 200<\/span>\r\n  Nwin<\/span> =<\/span> 20<\/span>\r\n  du<\/span> =<\/span> (<\/span>uend<\/span>-<\/span>uini<\/span>)<\/span>\/<\/span>(<\/span>Ndiv<\/span>*<\/span>Nwin<\/span>)<\/span>\r\n  do<\/span> for<\/span> [<\/span>i<\/span>=<\/span>0<\/span>:<\/span>Ndiv<\/span>\/<\/span>2<\/span>]{<\/span>\r\n    u<\/span> =<\/span> uini<\/span> +<\/span> du<\/span> *<\/span> i<\/span>\r\n    print<\/span> a<\/span>*<\/span>cos<\/span>(<\/span>u<\/span>)<\/span>,<\/span> a<\/span>*<\/span>sin<\/span>(<\/span>u<\/span>)<\/span>,<\/span> pitch<\/span>*<\/span>u<\/span>\r\n  }<\/span>\r\n  do<\/span> for<\/span> [<\/span>i<\/span>=<\/span>0<\/span>:<\/span>Ndiv<\/span>*<\/span>Nwin<\/span>]{<\/span>\r\n    u<\/span> =<\/span> uini<\/span> +<\/span> du<\/span> *<\/span> i<\/span>\r\n    if<\/span> (<\/span>sin<\/span>(<\/span>u<\/span>)<\/span> >=<\/span> -0.02)<\/span> {<\/span>print<\/span> a<\/span>*<\/span>cos<\/span>(<\/span>u<\/span>)<\/span>,<\/span> a<\/span>*<\/span>sin<\/span>(<\/span>u<\/span>)<\/span>,<\/span> pitch<\/span>*<\/span>u<\/span>}<\/span>\r\n      else<\/span> {<\/span>print<\/span> \"      \"<\/span>;<\/span>}<\/span>\r\n  }<\/span>\r\n  do<\/span> for<\/span> [<\/span>i<\/span>=<\/span>0<\/span>:<\/span>30<\/span>]{<\/span>\r\n    print<\/span> -3<\/span>+<\/span>dx<\/span>*<\/span>i<\/span>,<\/span> a<\/span>,<\/span> pitch<\/span>*<\/span>uend<\/span>\r\n  }<\/span>\r\n  print<\/span> \"          \"<\/span>\r\n  print<\/span> \"          \"<\/span>\r\n  print<\/span> 0<\/span>,<\/span> -<\/span>a<\/span>,<\/span> 0<\/span>\r\nset<\/span> print<\/span> \r\n\r\n# y < 0 \u5074\u306e\u30c7\u30fc\u30bf<\/span>\r\nset<\/span> print<\/span> \"solem.txt\"<\/span>\r\n  dx<\/span> =<\/span> 0.1<\/span>\r\n  do<\/span> for<\/span> [<\/span>i<\/span>=<\/span>0<\/span>:<\/span>30<\/span>]{<\/span>\r\n    print<\/span> -3<\/span>+<\/span>dx<\/span>*<\/span>i<\/span>,<\/span> -<\/span>a<\/span>,<\/span> pitch<\/span>*<\/span>uini<\/span>\r\n  }<\/span>\r\n  Ndiv<\/span> =<\/span> 500<\/span>\r\n  Nwin<\/span> =<\/span> 20<\/span>\r\n  du<\/span> =<\/span> (<\/span>uend<\/span>-<\/span>uini<\/span>)<\/span>\/<\/span>(<\/span>Ndiv<\/span>*<\/span>Nwin<\/span>)<\/span>\r\n  do<\/span> for<\/span> [<\/span>i<\/span>=<\/span>0<\/span>:<\/span>Ndiv<\/span>*<\/span>Nwin<\/span>]{<\/span>\r\n    u<\/span> =<\/span> uini<\/span> +<\/span> du<\/span> *<\/span> i<\/span>\r\n    if<\/span> (<\/span>sin<\/span>(<\/span>u<\/span>)<\/span> <=<\/span> 0.02)<\/span> {<\/span>print<\/span> a<\/span>*<\/span>cos<\/span>(<\/span>u<\/span>)<\/span>,<\/span> a<\/span>*<\/span>sin<\/span>(<\/span>u<\/span>)<\/span>,<\/span> pitch<\/span>*<\/span>u<\/span>}<\/span>\r\n      else<\/span> {<\/span>print<\/span> \"      \"<\/span>;<\/span>}<\/span>\r\n  }<\/span>\r\n  do<\/span> for<\/span> [<\/span>i<\/span>=<\/span>0<\/span>:<\/span>30<\/span>]{<\/span>\r\n#    print -3+dx*i, a, pitch*uend<\/span>\r\n  }<\/span>\r\n  print<\/span> \"          \"<\/span>\r\n  print<\/span> \"          \"<\/span>\r\n  print<\/span> 0<\/span>,<\/span> a<\/span>,<\/span> 0<\/span>  \r\nset<\/span> print<\/span> \r\n\r\nset<\/span> print<\/span> \"sole-vec.txt\"<\/span>\r\n  print<\/span> -1<\/span>.5<\/span>,<\/span> -1<\/span>,<\/span> pitch<\/span>*<\/span>uini<\/span>,<\/span> 1<\/span>,<\/span> 0<\/span>,<\/span> 0<\/span>\r\n  print<\/span> -0<\/span>.5<\/span>,<\/span> 1<\/span>,<\/span> pitch<\/span>*<\/span>uend<\/span>,<\/span> -1<\/span>,<\/span> 0<\/span>,<\/span> 0<\/span>\r\nset<\/span> print<\/span>\r\n\r\nsplot<\/span> \\\r\n  \"solep.txt\"<\/span> u<\/span> 1<\/span>:<\/span>2<\/span>:<\/span>3<\/span>:<\/span>(<\/span>-<\/span>$<\/span>2<\/span>)<\/span> w<\/span> l<\/span> lw<\/span> 6<\/span> lc<\/span> palette<\/span> notitle<\/span>,<\/span> \\\r\n  \"solem.txt\"<\/span> u<\/span> 1<\/span>:<\/span>2<\/span>:<\/span>3<\/span>:<\/span>(<\/span>-<\/span>$<\/span>2<\/span>)<\/span> w<\/span> l<\/span> lw<\/span> 6<\/span> lc<\/span> palette<\/span> notitle<\/span>,<\/span> \\\r\n  \"sole-vec.txt\"<\/span> using<\/span> 1<\/span>:<\/span>2<\/span>:<\/span>3<\/span>:<\/span>4<\/span>:<\/span>5<\/span>:<\/span>6<\/span>:<\/span>(<\/span>-<\/span>$<\/span>2<\/span>)<\/span> \\\r\n    w<\/span> vec<\/span> filled<\/span> head<\/span> lw<\/span> 6<\/span> lc<\/span> palette<\/span> title<\/span> \"\u96fb\u6d41\"<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n
\n
\n
\n
<\/div>\n
\n
\"\"<\/p>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n
\n
\n
In\u00a0[3]:<\/div>\n
\n
\n
\n
set<\/span> output<\/span> \".\/sole-rasen-kairo.svg\"<\/span>\r\nreplot<\/span>\r\nset<\/span> output<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n
\n
<\/div>\n
\n
\n
\u30bd\u30ec\u30ce\u30a4\u30c9\u306e\u5185\u90e8\u78c1\u5834<\/h3>\n