texファイルにgnuplotコマンドを埋め込む

はじめに

LaTeXの文書にグラフやデータのプロットを載せるとき，どのような方法を使っていますか？ MatlabやPythonのMatplotlibを使ったり，グラフ用紙に手書するなどして，最終的に画像化したものを\includegraphics{}するのがおそらく一般的な方法かと思います。

ところで，こんな経験はありませんか？

• 画像化したあとで，ミスに気付きやり直し
• 画像をどこに保存したか分からなくなった
• 画像のファイル名を指定するのが面倒
• 何のグラフか分からなくなった

これらを解決し得る一つの方法として，gnuplottexパッケージとその使い方を紹介したいと思います。

gnuplottexパッケージの何が嬉しいか

gnuplottexパッケージを使うと，texファイルの中にgnuplotのコマンドを書くことができます。タイプセットと同時にgnuplotが自動的に呼び出され，きれいなグラフを挿入してくれます。

gnuplotのコマンドを一元管理できる

本文の編集中に，シームレスにグラフの編集ができます。gnuplotファイルを用意する必要がありません。何のグラフだったか忘れても，近くの本文を読めば思い出せます。

作業が減る

自動化されているので，グラフを画像化して保存し，文書に取り込むという一連の作業もなくなります。ミスを見つけた時も，簡単に修正することができます。

グラフに数式を挿入できる

軸や凡例に数式を書くことができます。複雑な関数を書くのも，単位を書くのも簡単です。

環境

• Ubuntu 22.04 (<- Windows 11 + WSL2)
  • TeX Live 2023
  • gnuplot 5.4 patchlevel 2

にて確認。

必要な設定

gnuplotのターミナルには様々な種類があり，gnuplottexパッケージでも複数に対応しています。ここでは，tikzというターミナルを指定することにします^[1]。

gnuplottexパッケージ自体はTeX Liveに含まれていますが，ターミナルにtikzを指定するにはgnuplot-lua-tikzパッケージを準備する必要があります。これが今回やるべき唯一の設定です。

gnuplot-lua-tikzパッケージの準備

まずはスタイルファイルを置く場所を探します。ターミナルで

$ kpsewhich -var-value=TEXMFLOCAL

を実行します。TeX Liveのインストール先がデフォルトのままなら，/usr/local/texlive/texmf-local/が表示されると思います。全ユーザーに適用したくなければ，上でTEXMFLOCALの代わりにTEXMFHOMEとします。これらの場所は，TeX Liveを新しいバージョンにしても変更を受けず，ユーザーが手動でスタイルファイルを置くのに適しています。

次に，表示されたディレクトリに移動し，さらにその下のtex/latex/まで移動します。もしなければ作ってください^[2]。その下にroot権限でgnuplot-lua-tikzというディレクトリを作り，移動します^[3]。その後，root権限でgnuplotを起動します。

$ sudo gnuplot

gnuplotが起動出来たら，次のコマンドを実行します。

set terminal tikz createstyle
exit

これにより，スタイルファイルを含む4つのファイルが生成されます。

$ ls
gnuplot-lua-tikz-common.tex  gnuplot-lua-tikz.sty  gnuplot-lua-tikz.tex  t-gnuplot-lua-tikz.tex

最後に，配置したスタイルファイルを認識してもらうために，次のコマンドでls-Rファイルを更新します。

$ sudo mktexlsr

$ kpsewhich gnuplot-lua-tikz.sty
/usr/local/texlive/texmf-local/tex/latex/gnuplot-lua-tikz/gnuplot-lua-tikz.sty

これで設定は完了です。

使い方

実際に使ってみましょう。読み込むべきパッケージなど，詳しくはサンプルコードを見てください。gnuplot環境の中にgnuplotのコマンドを書くことができます。オプションとしてterminal=tikzを指定します。そのあとのterminaloptionsは必須ではありませんが，グラフのサイズやフォントを変更することができます。

サンプルコード1

適当なディレクトリに次のファイルを置いてください。エンジンはLuaLaTeXを使うことを前提に書いています。

\documentclass{ltjsarticle}   % お好きなドキュメントクラスでよいと思います。
\usepackage{luatexja}         % LuaLaTeXで日本語を使うのに必要です。
\usepackage{shellesc}         % shell escapeするのに必要です。
\usepackage{gnuplottex}       % 今回の主人公。
\usepackage{tikz}             % gnuplot-lua-tikzパッケージはtikzを使う前提なので。
\usepackage{gnuplot-lua-tikz} % 手動で配置したgnuplot-lua-tikzパッケージです。
\usepackage{siunitx}          % 本文中やグラフの軸に単位を書くために使います。
\usepackage{float}            % 図を出力する位置を自分で指定するため。

\begin{document}
$\sin x$と，その近似式をプロットしてみましょう。
\begin{figure}[H]
  \centering
  \begin{gnuplot}[terminal=tikz, terminaloptions={plotsize 9cm, 6cm}]
    set samples 1000
    set xrange [-5:5]
    set yrange [-3:3]
    set grid
    set xlabel "$x$"
    set ylabel "$y$"
    plot \
    sin(x) title "$\\sin x$", \
    x title "$x$", \
    x - x**3 / 6 title "$x - x^3/3!$", \
    x - x**3 / 6 + x**5 / 120 title "$x - x^3/3! + x^5/5!$"
  \end{gnuplot}
\end{figure}
\end{document}

実行は次のコマンドで行います。-shell-escapeオプションを忘れないでください。

$ lualatex -shell-escape example1.tex

下の画像のようなpdfが得られたでしょうか？

サンプルコード2

こちらは参考までに。

\documentclass{ltjsarticle}   % お好きなドキュメントクラスでよいと思います。
\usepackage{luatexja}         % LuaLaTeXで日本語を使うのに必要です。
\usepackage{shellesc}         % shell escapeするのに必要です。
\usepackage{gnuplottex}       % 今回の主人公。
\usepackage{tikz}             % gnuplot-lua-tikzパッケージはtikzを使う前提なので。
\usepackage{gnuplot-lua-tikz} % 手動で配置したgnuplot-lua-tikzパッケージです。
\usepackage{siunitx}          % 本文中やグラフの軸に単位を書くために使います。
\usepackage{float}            % 図を出力する位置を自分で指定するため。

\begin{document}
ここに本文を書きます。ここに本文を書きます。ここに本文を書きます。ここに本文を書きます。
ここに本文を書きます。ここに本文を書きます。ここに本文を書きます。ここに本文を書きます。

実験データをプロットしてみましょう。図\ref{fig:example2}のようになりました。
\begin{figure}[H]
  \centering
  \begin{gnuplot}[terminal=tikz, terminaloptions={plotsize 9cm, 6cm}]
    set samples 1000
    set log x
    set xtics 10
    set ytics 30
    set mytics 3
    set xrange [10:1e6]
    set yrange [-120:150]
    set grid
    set xlabel "周波数$f \\, \\unit{\\per\\hertz}$"
    set ylabel "位相差$\\theta \\, \\unit{\\per\\deg}$"
    
    R = 2206
    C = 105e-9
    L = 889.8e-6
    r = 22.2
    plot \
    (180 / pi) * atan(R * (2 * pi * x * L / (r**2 + (2 * pi * x * L)**2) - 2 * pi * x * C) / (r * R / (r**2 + (2 * pi * x * L)**2) + 1)) title "理論曲線（$r$を考慮）", \
    (180 / pi) * atan(R * (1 / (2 * pi * x * L) - 2 * pi * x * C)) title "理論曲線（$r$を無視）", \
    'exp1.dat' using 1:4:5 with yerrorbars title "測定値"
  \end{gnuplot}
  \caption{微分回路に正弦波を入力したときの電圧の位相差。縦軸に位相差を取り，横軸に周波数の対数を取った。理論曲線は，インダクタンスの内部抵抗$r$を考慮した場合と無視した場合の2つを書き込んだ。}
  \label{fig:example2}
\end{figure}
数式を埋め込むのも簡単です。
\end{document}

是非レポートなどで使ってみてください！

参考

• TeX Wiki：日本語のLaTeX関連の情報が集まっています。
• gnuplottexの公式ドキュメント：ここに書いてある情報がおそらく一番正しい。

他にも参考になる記事を2つ挙げておきます。1つ目の記事では，TikZのexternalライブラリを使ってコンパイルにかかる時間を短縮する方法について説明されています。

https://qiita.com/satl/items/0c11c8808b43f806ee21

https://qiita.com/kou-JP/items/309dfd3729f602cd0042