Microsoft PowerPoint - 14Gnuplot.ppt

Size: px

Start display at page:

Download "Microsoft PowerPoint - 14Gnuplot.ppt"

いとはあかさか
5 years ago
Views:

1 Gnuplot との連携 Gnuplot による関数の描画 Gnuplot によるデータのプロット Gnuplot による 3-D グラフの描画 Gnuplot による 3-D データのプロット今日のポイントグラフ描画ソフト gnuplot を体験してみよう Gnuplot とはグラフ描画専用のフリーウェア Windows 版の名称は wgnuplot インストールはフォルダごとコピーするだけ専用のコマンドウインドウによりおもにコマンドで操作日本語のオンラインヘルプや Web での解説ページも充実 Gnuplot の起動と終了 gnuplot の起動は [ スタート ] [ すべてのフロクラム (P)>] [gnuplot] を < クリック > [ スタート ] [ すべてのフロクラム (P)>] [gnuplot] を < 右クリック > して [ 送る (N)] より [ テスクトッフ ( ショートカットを作成 )] を選びショートカットアイコンを作成ショートカットアイコンを < タフルクリック > すると gnuplot が起動左上の [File] [Eit] または右上の [] X ホタン 3 により終了フォントが小さい場合 gnuplot のウィントウを < 右クリック > し 3つめの [ Choose Font Font... ] を選択しフォントウィントウを開くフォント名 (F): で [Terminal] を選択しサイス (S): で [4] を選んでおく [ OK ] ホタンを押してフォントウィントウを閉じる再び gnuplot のウィントウを < 右クリック > し一番下の [Update wgnuplot.ini] を選んで設定を保存 4

2 練習関数グラフを描いてみよう plot **+3, **3-75* プロンプトべき乗 **+3 **3-75* グラフ描画のための Tips() X 軸ラベル Y 軸ラベルは set label "" set ylabel "y" replot 凡例 (Key) を変えるには Aes メニュー [Replot] ホタン plot **+3 title "y=^+3", > **3-75* title "y=^3-75" プロットの線を線と点に変えるには plot **+3 with linespoints, > **3-75* with linespoints 矢印キーで前のコマンドを呼び出し編集できる Plot メニュー改行して続ける記号 6 グラフ描画のための Tips() X 軸範囲 Y 軸範囲の指定は set range [-5:5] set yrange [-:] replot 目盛線をつけるには set grid 対数グラフにするには set logscale y set コマンドを解除するには unset grid Aes メニュー Chart メニュー Grid on Aes メニュー Logscale yz 練習対数グラフを描いてみよう plot /(+*) set logscale y set range [.:] set yrange [.:] set label "" set ylabel "y" set grid replot Lorentz 関数 Chart メニュー Grid off 7 8 y... /(+*)..

3 sinc () = (sin π/π) sinc () の対数グラフリサージュ図形 sinc() = sin(pi*)/(pi*) plot sinc()** set logscale y set range [.:] set yrange [.:] set label "" set ylabel "y" set grid replot y... 関数定義 sinc()** set parametric plot cos(5*t), sin(7*t) set trange [-pi:pi] set grid unset key.5 replot 凡例なし (No Key) -.5 媒介変数 t 練習 3 データのプロット plot 'Z: nyumon q4.csv' with lines 凡例を書き換える場合 plot 'Z: nyumon q4.csv' > title "y=^+3" with lines さらに線を線と点にかえる plot 'Z: nyumon q4.csv' > title "y=^+3" with linespoints 複数データのプロット plot 'Z: nyumon q4_s.csv' with lines, > 'Z: nyumon q4_s.csv' using :3 with lines 凡例を書き換える場合列目と 3 列目でプロット plot 'Z: nyumon q4_s.csv' > title "^+3" with lines, > 'Z: nyumon q4_s.csv' using :3 > title "^3-75" with lines

4 練習 4 3-D グラフの描画 R(,y)=sqrt(*+y*y) gauss(,y)=ep(-r(,y)**) splot gauss(,y) set range [-:] set yrange [-:] set isosample 4,4 set label "" set ylabel "y" replot gauss(,y) 関数定義 default のでは足りない y 3 練習 5 等高線の描画 R(,y)=sqrt(*+y*y) gauss(,y)=ep(-r(,y)**) splot gauss(,y) set range [-:] set yrange [-:] set contour base unset surface set view, set cntrparam levels set label "" set ylabel "y" replot y D/ 等高線グラフのための Tips Milk drop 陰線処理 set hidden3d カラーマップ表示 set pm3d カラーパレットの指定 set palette gray set palette color gauss(,y) y Milkdrop(,y) = ep( ep( -R(,y) ) ( ep( cos( R(,y) )** ) + 8*sin( R(,y) )** + *sin( *R(,y) )**8 ) ) Milkdrop(,y) y

5 3-D プロット用のデータ構造 Ecelとは異なり -Dプロット用データを縦に並べることで3-Dプロット用データとしている -Dプロット用データの区切りには < 空行 > を用いる 3-Dプロットデータをさらに複数重ねるときには行分の < 空行 > を境界とする練習 6 3-D データのプロットまず f (, y)=ep( ( + y )), -, y のgnuplot 用データを作成し gauss.dat に保存するつぎに gnuplotにより gauss.dat の3-D プロットを描く 7 8 プログラム例 (gaussg.c) #include <stdio.h> #include <math.h> int main(void) { int i, j, n=; double, y, z; for (j = -n; j <= n; j++) { y =.*j; for (i = -n; i < n; i++) { =.*i; z = ep(-(*+y*y)); printf("%.f, %.f, %.3f n",,y,z); printf(" n"); return ; 実行手順 Z: nyumon>cl gaussg.c... Z: nyumon>gaussg > gauss.dat このあと gnuplot を起動し splot 'Z: nyumon gauss.dat' > with lines 9

6 スキルアップタイム :(gaussg.c) スキルアップタイムの出力例. 最初に f (r) = ep( r ) cos(πr) 絶対値は fabs( ) "gauss.dat" のデータを出力したあと回改行し続けて f (r) = ep( r ) + を出力するプログラムを作成し gauss.dat に保存せよただし r = +y,, y.. つぎに gnuplotにより gauss.dat の3-D プロットを描けヘルプその他 Help メニュー : コマンドの詳しい使い方などコマンドは誤解されない範囲で短く省略可 plot p, splot sp, replot rep with w, lines l, title t, using u range r, label l, logscale logs palette pal, hidden3d hid などデモは C: Gnuplot demo all.dem 3 faceg.dat の表示 set view 8, 真下から見る unset surface メッシュを非表示 unset tics 軸も非表示 unset ytics 正方形にセット unset ztics グレイスケールで表示 set size square set pm3d gnuplot 用にデータを並べ替えた faceg.dat を3D 表示 set palette gray splot 'Z: nyumon faceg.dat' 4

7 face.dat faceg.dat の変換 face.dat faceg.dat の変換 #include <stdio.h> face4g.c #include <string.h> int main(void) { int i; char *p,s[5]; while (fgets(s,5,stdin)!=null) { p=strtok(s," "); if (p && (*p!=' ')) printf("%s n", p); for (i=;i<8;i++) { p=strtok(null," "); if (p && (*p!=' ')) printf("%s n", p); printf(" n"); return ; 5 Z: nyumon>cl face4g.c Z: nyumon>face4g < face.dat > faceg.dat 6

Microsoft Word - gnuplot

Microsoft Word - gnuplot GNUPLOT の使い方 I. 初期設定 GNUPLOT を最初に起動させたときの Window の文字は小さいので使い難いそこで文字フォントのサイズを設定します 1.GNUPLOT を起動させます ( 右のような Window が起動します ) 2. 白い領域のどこでも構わないのでポインタを移動してマウスの右ボタンをクリックします ( 右のようにメニューが起動します ) 3. Choose