極座標プロットのあれこれ データを2次元極座標でプロットする． 通常の2次元データはデカルト座標系の(X,Y)の点の組として与えられていま す．これを極座標で表示します．データファイルの書き方は，普通の2次元データと同じです．極座標表示では， (X,Y)の組が，(θ,r)として取り扱われます．つまり，X=角度，Y=動径長さ， です．角度の単位は，デフォルトではradianで，0〜2πでプロットは一周 します．度(degree)を使いたいときは， set angles degrees と します． 次の様な2次元データがあるとします．これはデカルト座標系でプロッ トしたもので，Xは0〜180度までの角度，Yは-1〜1で変化する量です． これを極座標表示します． gnuplot> set polar dummy variable is t for curves gnuplot> set a

ガウス分布の外形もわからない「「interpolate」なし」の場合と比べて、パラメータの数値が大きいほど、滑らかになっていることがわかります。 「interpolate 5,5」とすることで、たったに対するデータでも、ガウシ分布の外形を理解することができます。 ただし、「interpolate 10,10」「interpolate 15,15」では、ガウス分布の中心に本来はないはずの十字の線があるように見えます。 つまり、「interpolate 」オプションを利用した補間時の問題点といえるかもしれません。大きくしすぎないようにすべきですね。 gnuplot によるカラーマップ表示の最適値 本稿では terminalで「gif」を選択し、gifファイルのサイズ「600×480」ピクセルに指定しています。 その際に、上記ののような平面の分割数がでは足りません。 そこで、元データでに分割に対

Ruby Gnuplot is a pipe-based interface to the ever popular gnuplot package. Through this interface almost every capability of gnuplot is useable from within Ruby. The goals of this package are as follows (in order of my priorities): Provide a general object oriented interface for scientific plotting that completely hides the details of the underlying implementation. ie. The fact that this uses gnu

Gnuplot入門 第0.08版

次 前 GNUPLOTの第一歩 GNUPLOTで簡易アニメーションを作る gnuplotのコマンド(plotやsplotなど)をファイルに書いておいて、それを gnuplotに読ませることで実行することができます。コマンドの書かれた ファイルをスクリプトといいます。 スクリプトsample1をカレントディレ クトリに準備しておいて % gnuplot sample1 とすることで、data1のグラフが表示されます。pause -1と書くことで、 リターンキーが押されるのを待たせています。リターンキーを押してみる と次のグラフ data2 が表示されます。 もし、pause -1を書かれていないとすると、data1とdata2のグラフを リターンを待たずに連続して表示し、その後終了します。その例はスク リプトsample2です。 % gnuplot sample2 pause -1の代わりに

3項演算子 (A ? B : C) の使い方． 3項演算子は，C言語のそれと同様に働きます，と言うとあまりに不親切な ので，もう少し詳しく説明しましょう．3項演算子とは A ? B : C のように書 かれたもので，Aは条件式，BとCは実行文です．これは，Aが真ならBを実行し， そうでなければCを実行せよ，という意味を表します．X>0の範囲ではexp(-x)， Xが負の領域ではexp(4*x)となる関数をプロットしてみましょう． gnuplot> set xrange [-5:5] gnuplot> plot x>0 ? exp(-x) : exp(4*x) 条件式を入れ子にすれば，もう少し複雑な関数を定義することも可能です． Xの絶対値の範囲に応じてY値が離散的に変化する関数を作ってみましょう． gnuplot> f(x)= (abs(x)<1) ? 100 : ( (abs(x)<2

[Rubyforge Project page] [ChangeLog] [Authors] [License] History and Background Gnuplot is a program that has a rich language for the generation of plots. It has a unique place in academia as it was one of the first freely available programs for plot generation. I started using gnuplot over 10 years ago while pursuing my Master's degree in Physics and have been using it actively ever since. Versi

京都大学固体量子物性研究室 米澤進吾のホームページです。

シェルスクリプトによるデータ処理 （一部工事中、順次拡張予定） シェルスクリプトの基本 シェルスクリプトは、言語というよりはコマンド群のようなもので、 CやFortranと異なりコンパイルせずに使うことができる。 ここでは、大概のUNIXシステムで使うことのできるBシェルを扱う。 シェルスクリプトの一番簡単な使い方は、ただコマンドを羅列することである。 #!/bin/sh echo "This is B-shscript." ls cd .. 上記のファイルを例えばtest.shとする。実行するには、ファイルに実行許可を与えてからファイル名を打つ。 % chmod +x test.sh % test.sh (環境によっては、% ./test.sh) echo文が実行され、"This is B-shscript."という出力がなされ、ディレクトリ内のファイル一覧が出され、最後にディレクトリが

作業のバッチ化 ここでは，数値計算編の最後に保存したプロットファイル "output.plt" を使って，「その後」の仕事を楽に する方法を概説します．このファイルを覗いてみると，作業中に設定した部分 (軸の範囲やプロットするファイル名)と，その他のデフォルトで設定されてい る部分が混在しているのが分かると思います．ちょっと見づらいので，デフォ ルトで設定されている部分を消して，自分で設定したものだけを残してみましょ う．但し，termianlとplotの部分は長くなりすぎるので，継続行記号(backslash)を 使って折り返しています． #!/usr/local/bin/gnuplot -persist # set terminal postscript landscape noenhanced monochrome \ # dashed defaultplex "Helvetica

wgnuplotでシンボリックリンクやショートカットは読み込めない。 gnuplot> load 'symlink.plt' Cannot open load file 'symlink.plt' ハードリンクなら読み込める。 gnuplot> load 'hardlink.plt' グラフを指定した範囲だけプロットする 例えば sin(x)を5<x<10 の範囲でプロットするには、三項演算子を用いて f(x)=(x>5 && X<10) ? sinx(x) : 1/0 plot [0:20] f(x) とします。つまり指定した範囲以外を無限大（1/0）に飛ばしてしまいます。 ステップ関数 三項演算子を用います。 f(x)= x>0?1:0 動画（GIFアニメーション）を作る まずgnuplotで表示させるためのデータを作成しておく。ここではサンプルとして拡散方程式を使う。拡散方程式を解く