gnuplot 2 3次元処理工学 データファイルのプロット まずデータファイルの在るディレクトリに移動 gnuplot> cd '………' gnuplot> plot 'test.txt' gnuplot> plot 'test.txt' using 1:3 gnuplot> plot 'test.txt' with lines gnuplot> plot 'test.txt' with points gnuplot> plot 'test.txt' with linespoints with points w p と省略できる with lines w p with linespoints w lp データファイルの書式(1) # yの値だけ 2.2 2.0 1.9 1.5 1.3 # で始まる行はコメント (無視される) 一行に一つだけ数値がある場合、 その数値が y の値とみなされ、 x の値は 0, 1, 2, … が用いられる。 gnuplot> plot 'test.txt' w lp データファイルの書式(2) # x, 0.0 0.5 0.8 1.5 1.9 y 2.2 2.0 1.9 1.5 1.3 一行に2つの数値がある場合、 1列目を x、2列目を y に用いる。 gnuplot> plot 'test.txt' w lp データファイルの書式(3) # x, 0.0 0.5 0.8 1.5 1.9 y1, 2.2 2.0 1.9 1.5 1.3 y2 0.4 0.8 1.0 1.2 1.8 using 1:2 1列目を x、2列目を y に用いる gnuplot> plot 'test.txt' using 1:2 w lp, 'test.txt' using 1:3 w lp fitting データに式をあてはめるのに fit コマンドが使える。 たとえば、y = ax + b という直線で fit するなら gnuplot> f(x) = a*x + b gnuplot> fit f(x) 'test.txt' via a,b (結果が画面に出力される) gnuplot> plot f(x), 'test.txt' 常に fit が成功するとは限らない。 収束しない場合、パラメタ a, b などの初期値を与 えるとうまくいくこともある。 fitting (例) damp.txt にあるデータを t / f (t ) ae sin(t ) という関数で fit してみよ。 gnuplot> set dummy t とすると、2次元プロットの横軸を t にできる。 動的粘弾性測定 (t ) 0 cost (t ) 0 cost 動的粘弾性測定(2) (t ) 0 cost (t ) 0 cost 0 cos cost 0 sin sin t 0 G() cost G() sin t 0 G( ) cos 0 0 G( ) sin 0 動的粘弾性測定(3) (t ) 0 cost Re (t ) it (t ) 0e (t ) 0 cost Re (t ) (t ) 0e i (t ) 0 i it e 0e 0 (t ) G () (t ) 0 i G ( ) e G( ) iG( ) 0 動的粘弾性測定(4) (t ) ReG () (t ) 0 ReG () exp(it ) * * 0 ReG() iG()cost i sin t 0 G() cost G() sin t 演習 stress.txt を (t) の測定データ( t の単 位は秒、 の単位は Pa )とみなし、 (t ) 0 cost で fit し、 0 , を求める。 これから、この周波数での 貯蔵弾性率G'() 、損失弾性率G"() を求めよ。ひずみ振幅は 0 0.1 とせよ。 パラメータプロット gnuplot> set parametric gnuplot> set parametric gnuplot> plot cos(t), sin(3*t) gnuplot> splot sin(v)*(3+cos(u)), cos(v)*(3+cos(u)),sin(u)
© Copyright 2025 ExpyDoc