gnuplot を使おう。 その3 (955ﾚｽ)
上下前次1-新
1(4): 2011/08/25(木) 12:15:50.61 AAS
グラフ描画ソフト gnuplot のスレ。
gnuplot 公式 Web ページ
外部ﾘﾝｸ:www.gnuplot.info
前スレ
gnuplot を使おう。 その2
2chｽﾚ:unix
2(1): 2011/08/25(木) 12:17:02.26 AAS
■過去スレ
gnuplot を使おう。
2chｽﾚ:unix
■役立ちそうなページ
gnuplot tips (not so Frequently Asked Questions)
外部ﾘﾝｸ:t16web.lanl.gov
gnuplot (Takeno Lab) 日本語まぬある
外部ﾘﾝｸ:takeno.iee.niit.ac.jp
GNUPLOTの第一歩
外部ﾘﾝｸ:lagendra.s.kanazawa-u.ac.jp
gnuplot tips
外部ﾘﾝｸ[html]:wwwnucl.ph.tsukuba.ac.jp
3: 2011/08/25(木) 14:17:32.56 AAS
>1 乙
4: 2011/08/25(木) 19:18:05.90 AAS
>1 乙
5: 2011/08/26(金) 16:37:27.77 AAS
>>1 外部ﾘﾝｸ:gnuplot.sourceforge.net
いまは，こっちのほうが公式ホームでなかったですか。
外部ﾘﾝｸ:www.gnuplot.infoはミラーになったのでは？
6: 2011/08/26(金) 16:39:05.94 AAS
ここも忘れてはいけないと思いますが。
外部ﾘﾝｸ[html]:takeno.iee.niit.ac.jp
7: 2011/08/26(金) 17:13:37.59 AAS
>>2
gnuplot tips
wwwnucl.ph.tsukuba.ac.jp/~inakura/gnuplot/gnuplot.html
はリンク切れ。
8: 2011/08/26(金) 18:00:43.10 AAS
>>1
おつ
9(1): 2011/08/30(火) 13:24:32.31 AAS
うに板に即死判定あったっけ？
10(1): 2011/08/30(火) 14:02:35.76 AAS
Unix板の
データ解析ツールoctaveを語ろう
2chｽﾚ:unix
981で即死しました。
11(1): 2011/08/30(火) 14:22:20.82 AAS
>>9
10未満。
12: 2011/08/30(火) 15:52:19.28 AAS
>>10
いや、それは即死じゃない
980越えて24時間レスがなかったら落ちるのは全板共通
即死ってのは立ったばかりのスレが引っかかるやつ
>>11
ありがと
じゃもう大丈夫か
13: 2011/09/12(月) 17:26:00.33 AAS
Release plans: Version 4.4.4 and beyond
もうすぐ，4.4.4がでるようです。
4.4の次は4.6でなくて，Version 5になるのかな。
外部ﾘﾝｸ[php]:sourceforge.net
14(3): 2011/09/20(火) 11:45:58.77 AAS
前スレの972で
>wgnuplotで起動して最初にplotするときだけ時間が異常にかかります
ですがgnuplot\binary\etc\fontsにあるfonts.confを編集すればTEMPで指定してあるところ以外にも自由におけるようです
外部ﾘﾝｸ:swordfish.s53.xrea.com gimpですが同様の方法で変更できました
15(1): 2011/09/20(火) 12:33:57.94 AAS
>>14
<cachedir>WINDOWSTEMPDIR_FONTCONFIG_CACHE</cachedir>
<cachedir>~/.fontconfig</cachedir>
のことですか？
16: 2011/09/21(水) 12:49:41.94 AAS
>>15 です。>>14 の
外部ﾘﾝｸ:swordfish.s53.xrea.com
を見てなかったm(__)m。
これに書いてあるわけですね。
17: 2011/09/21(水) 13:37:37.09 AAS
>>14 前スレの972の現象は，cvs版のバイナリ（Win版）ではおきないようです。
外部ﾘﾝｸ:www.tatsuromatsuoka.com
18(1): 2011/10/07(金) 01:02:24.78 AAS
Fortranのバイナリ形式で出力したデータをgnuplotでプロットしたいのですが、
どう書けば良いかわからないので教えてください。
以下のように出力したデータ
（円柱座標系、R方向,Z方向に0.02mごと、theta方向に1度ことに出力）を、
k=0で指定されるR-Z断面について、Bz成分の強度のグラフをpm3dで描きたいです。
write(10) intA, intB, intC, intD ! 4byte＋単精度整数×4＋4byte＝24byte
write(10) realA, realB, realC, realD ! 4byte＋倍度浮動小数点×4＋4byte＝40byte
do k = 0,359 ! ← 円柱座標系のtheta方向
write(10) ((Bx(i,j,k), By(i,j,k), Bz(i,j,k), i=0,100), j=0,100)
end do
これをプロットしようと、以下の様に途中まで書きましたが、
formatやusingをどのように書けば良いのか教えてください。
splot 'fort.10' binary \
skip=64 \
format="%*int%？？？%*int" \
endian=big \
array=(101,101) \
scan=xy \
dx=0.02 \
dy=0.02 \
using ？ with pm3d title'B_{x}'
gnuplotのバージョンは4.4、Fortranのコンパイラはifortを使っています。
19: 18 2011/10/07(金) 01:19:15.89 AAS
追記
Bx、By、Bzはいずれも倍精度実数です。
20(1): 2011/10/18(火) 21:28:53.26 AAS
windowsのgnuplotのwxtでdashを使いたいんだけど無理なのかな？
アンチエイリアス機能が便利でそのままコピペして使いたい
他の方法でdash扱ってると面倒すぎる
21(1): 2011/10/19(水) 12:40:40.80 AAS
>>20
set term wxt dashed
test
22: 2011/10/20(木) 23:25:29.46 AAS
>>21
まさにすべて解決した気分。この先どうなるかわからないけど
非常に有難い
23: 2011/11/11(金) 19:11:51.83 AAS
データのうち、一部分だけ別のpoint typeでプロットしたいのですができますか？
全ての点をlineでつなぎぐので別々にプロットするのはなしでお願いします。
24: 2011/11/11(金) 19:51:25.96 AAS
すいません自己解決しました。
25(2): 2011/11/20(日) 02:23:59.37 AAS
Windows7でgnuplot4.4.4を使っているのですが、EPSを作成する際にラベルに日本語を使用した場合、
数字などの半角文字の文字間隔が広くなってしまい、見栄えが悪く困っています。
例えば
外部ﾘﾝｸ[html]:www.obihiro.ac.jp
ここに解説されている「日本語を含むEPSファイルの保存」のとおりにEPSを作った場合、
画像ﾘﾝｸ
のような出力が欲しいのですが、私の環境では
画像ﾘﾝｸ
のように出力されてしまうわけです。
どなたか解決方法を教えていただけないでしょうか。
26(2): 2011/11/20(日) 06:22:39.52 AAS
手元にgnuplot-4.4.4がないので，cvs版（4.5）でやってみました。
フォントを"GothicBBB-Medium-83pv-RKSJ-H" で指定して
shift-jisでやったところお望みに近いグラフができました。
ちなみにWindows XPです。
27(2): 2011/11/20(日) 06:30:44.35 AAS
"GothicBBB-Medium-UniJIS-UTF8-H"でUTF-8で保存でもうまくいきました。
ghostscriptのバージョンは9.00です。
28(2): 2011/11/20(日) 07:13:11.15 AAS
gnuplot-4.4.4をダウンロードしてやってみました。
私の環境ではEUCでは，postscriptファイルが壊れて見れませんでした。
"GothicBBB-Medium-83pv-RKSJ-H" shift-jis
"GothicBBB-Medium-UniJIS-UTF8-H" UTF-8
でうまくいきました。ghostscriptのバージョンなどにもよるかもしれませんね。
29(1): 25 2011/11/20(日) 17:23:03.60 AAS
>>26-28
今まではGhostscript9.02を使っていたのですが、9.00を使ったところ問題なく出力されました。
ありがとうございました。
30: 2011/11/26(土) 11:31:26.78 AAS
>>25-29
外部ﾘﾝｸ[html]:takeno.iee.niit.ac.jp
に詳しい解説がなされています。
31: 2011/12/03(土) 20:25:23.89 AAS
こんなサイトがあるようです。
なかなか凝ったプロットがかけるようです
外部ﾘﾝｸ[html]:www.phyast.pitt.edu
32: 2011/12/03(土) 22:52:51.17 AAS
うお、すげえ
33: 2011/12/05(月) 15:16:55.97 AAS
Windowsで特定のアプリケーションのCPU使用率と測定開始時間から終了時間までを測定しGnuplotを使いグラフにしたいのですがCPU使用率を計るにはどうしたらよろしいでしょうか?
34: 2011/12/05(月) 15:32:24.45 AAS
gnuplot関係なさ杉ﾜﾗﾀｗパフォーマンスカウンターでできるよ
35: 2011/12/27(火) 20:09:13.44 AAS
gnuplot -p scriptfile.gp
なんて感じで、スクリプトファイルを読み込ませて起動したとします。
この時、scriptfileg.gp の中で gnuplot が読み込んだスクリプトファイルの名前を取得する事ってできるでしょうか？
36(1): 2011/12/28(水) 02:55:50.93 AAS
俺は label とかに TeX 使いまくてたから、今だに plus 使ってるよ。。。
どうしたら良いのか。。。
37: 2011/12/28(水) 12:42:42.20 AAS
>>36
epslatexとかどう？
gnuplot走らせるところから、出力されたTeXファイルにヘッダ加えてコンパイル、目的の形式に変換
までスクリプトで自動化すればそう面倒ではないけど
38: 2012/01/02(月) 03:58:17.20 AAS
いや、４系の高機能に惹かれるのだが、これまでの資産のせいで乗り移れないってこと。
両方入れれば良いのだが、FreeBSDのポーツだと排他的だし。。。
39: 2012/01/06(金) 02:20:08.50 AAS
何故pyxplotの話にならないのだろう(ぼそ)
40(1): 2012/01/17(火) 08:20:50.04 AAS
グラフをプロットして出力するとデフォルトだと横長の画像が出力されるけど
これを4:3の大きさに指定したり、グラフの枠を指定して(できればピクセル単位で)出力ってできる？
x軸y軸で4:3を比率調整で作っても、それがx軸y1軸y2軸になるとy2軸側の文字が入って
画像としては4:3だけどグラフが4:3じゃなくなってしまう
41: 2012/01/18(水) 18:59:32.18 AAS
出力された画像を4:3に加工するんじゃ駄目なの？
42: 2012/01/18(水) 20:09:00.96 AAS
>>40
set term wxt size 600,800
とか。
詳しくは使用する terminal のヘルプを見てくれ。
43: 2012/01/21(土) 04:41:34.75 AAS
gnuplot 4.6-rc1がでました。
Windows，Cygwin，DJGPPのバイナリも配布されていま。
Windows 版は，今回からインストーラー形式になっています。
また，日本語環境でインストールはデフォールトは日本語環境となります。
初心者には，インストールがずいぶん簡単になるようです。
外部ﾘﾝｸ:sourceforge.net
44: 2012/01/23(月) 05:20:15.28 AAS
理論上+∞に発散して閾値を超えると-∞から始まる曲線を描写しようと考えてるんだけど
実際のデータは無限大に行くわけでもなく、普段0.001程度の値が1まで行って-1まで
一気に落ちるって感じになってる。
こういうデータを無限に発散しているように見せるテクってある？
+∞から-∞の点を繋ぐ線を消したい
45(1): 2012/01/23(月) 10:21:37.76 AAS
2本に分けて描く
46(1): 2012/01/23(月) 13:38:40.04 AAS
(abs(1/x)<c)?(1/x):(0/0)
定数cとsampleの値を適当に調節する
47(1): 2012/01/25(水) 17:52:13.49 AAS
>>45
無限に発散する点と普通に0を通過していく点があるから
2本に分けるにも結構難しかった
もっと効率よくできないものか
>>46
計算ソフトで出力されたデータをgnuplotに読み込ませて描写しているから
関数表現だと無理だと思う
それをプロットに応用できればいいのだが…
計算能力の問題で＞普段0.001程度の値が1まで行って-1まで一気に落ちる
という結果になってる。
あくまで例で綺麗に1から-1まで一瞬に移動してくれるわけでもなく
理論上無限だけど、その点は効果が薄いために
他の発散は1から-1まで下がるのに対して0.05から-0.05に下がる点があったり
発散の度合いにも差があって、発散が弱いと普通に0を通っているのか発散しているのか見分けが難しい
48: 2012/01/25(水) 18:23:29.11 AAS
>>47
>計算ソフトで出力されたデータをgnuplotに読み込ませて描写しているから
> 関数表現だと無理だと思う
Help ternaryをみると次のようなものがあるから，三項演算子とusingで
可能だと思います。
plot 'file' using 1:( $4<0 ? 1/0 : ($2+$3)/2 )
49(2): 2012/01/31(火) 14:14:17.57 AAS
sin 関数のグラフで x 軸のメモリにπを入れる方法ってない？
50(3): 2012/01/31(火) 18:13:06.19 AAS
>>49
gnuplot のturial.pdf のFigure 5 にlatex terminal での例がある。
set terminal latex
set output "eg4.tex"
set format y "$%g$"
set format x "$%.2f$"
set title 'This is $\sin(x)$'
set xlabel "This is the $x$ axis"
set ylabel "$\\sin(x)$"
unset key
set xtics ('$-\pi$' -pi,\
'$-\frac{\pi}{2}$' -pi/2,\
"0" 0,\
'$\frac{\pi}{2}$' pi/2,\
'$\pi$' pi)
plot [-pi:pi] [-1:1] sin(x)
Latex terminalが使えないならpostscript enhanced でやることが考えられる。
51: 2012/01/31(火) 18:15:48.73 AAS
>>50
turial.pdf ==> tutorial.pdf
52(1): 49 2012/02/01(水) 04:48:42.27 AAS
>>50
ありがとう！
よく，emf ファイルで出力するからこの方法で大丈夫か
心配だったけど emf も enhance 付ければ出力できるね！
53: 2012/02/01(水) 11:24:59.88 AAS
>>50 >>52
emfで enhancedでやってみた。
set terminal emf enhanced
set output "sinpi.emf"
set title 'This is sin(x)'
set xlabel "This is the x axis"
set ylabel "sin(x)"
unset key
set xtics ('-{/Symbol p}' -pi, \
'-{/Symbol p} / 2' -pi/2, \
'0' 0 , \
'{/Symbol p} / 2' pi/2, \
'{/Symbol p}' pi)
plot [-pi:pi] [-1:1] sin(x)
54(1): 2012/02/10(金) 21:40:03.29 AAS
"Skipping data file with no valid points"のエラーってなんなんですか?
どんなコードもこれが出てきてしまってグラフィックができないです。
55: 2012/02/10(金) 23:25:38.50 AAS
英語で書かれているとおりだよ。
無効なデータをスキップしましたってこと。
データのフォーマットがおかしいんじゃないの？
56: 2012/02/12(日) 08:06:35.93 AAS
>>54
実行できないコードの一つでもさらせば，解答がつくかも
57(1): 2012/02/15(水) 01:35:35.43 AAS
x データが 48〜120，218〜289 という具合に
データの塊が2つにわかれているんだけど，
これを省略線(波線2つ)を入れて省略するみたいに
表現できないかな？
58(1): 2012/02/15(水) 07:11:17.68 AAS
>>57 データではないけどそういう例はある
外部ﾘﾝｸ[html]:www.phyast.pitt.edu
59: 2012/02/15(水) 15:07:57.76 AAS
>>58
ありがとう！
簡単にできるわけじゃないのねｗ．
とりあえず，じっくり読み考えながらやってみるわ．
60: 2012/03/12(月) 15:18:16.63 AAS
gnuplot version 4.6 released!
61(1): 2012/03/12(月) 19:26:06.94 AAS
↑
あららら・・・
インストール
変わっちゃったね？
62(1): 2012/03/13(火) 08:06:35.29 AAS
gnuplot 4.6.0がでてます。
外部ﾘﾝｸ:sourceforge.net
63: 2012/03/13(火) 08:09:50.61 AAS
>>62 うああ，かぶったm(__)m。
>>61 Windows版ですね。インストーラーがついてますよね。
でも，従来型のzipもあります。
64: 2012/04/09(月) 22:25:18.79 AAS
set yrange [0:2e+9]
set format y "%s%c"
としたとき、Ｙ目盛りが
期待： 500M, 1G, 1.5G, 2G
実際： 500.000000M, 1.000000G, 1.500000G, 2.000000G
となります。
ソースを見ると、gprintfの%sでは、stdioのprintfの%fが使われているっぽい。
%gなら余計な０は表示されないのに、%fが使われている積極的な理由ってあるのかな？
とりあえずソースいじれば何とかなりそうだけど、
自分だけで使ってるわけじゃないので、それはやりたくないんだよな。
65: 2012/04/28(土) 04:15:16.94 AAS
4.4以降でterminal windowsとかでemfを吐くと
それをパワポ上でグループ解除すると枠線らしきものが
かなりのでかさで現れて酷いことになる
wxtだと確かに綺麗だけど
ビットマップ扱いで拡大するとボヤけたりepsに変換すると汚くなる
ボヤけないようにemfで保存すると枠線が・・・となる
term windowsだと破線の太さをgnuplot上で太くできないから
仕方なくwxt使ってるんだけど
wxtで出力された見た目のままベクター形式で出力できないのだろうか
66: 2012/04/28(土) 18:53:34.37 AAS
> wxtで出力された見た目のままベクター形式で出力できないのだろうか
pdfcairo
powerpointで読み込む方法があるのか知らんけど
67(1): 2012/05/03(木) 18:01:13.39 AAS
ylabelを左じゃなくて、軸の上に表示したい
68: 2012/05/03(木) 18:37:55.70 AAS
>>67
ylabel ではなく，label をグラフの外の座標に
設定すれば良いんじゃないだろうか？
69(1): 2012/05/10(木) 16:20:51.55 AAS
gnuplot> set pm3d map
gnuplot> splot [:2*pi][:2*pi] sin(x)*cos(y)
ってやるとy軸の範囲が0~7になって上に空白が入っちゃうのってどうやって直したらいい？
70(1): 2012/05/10(木) 23:20:52.00 AAS
>>69
下限を設定したらうまくいったよ
71: 2012/05/11(金) 06:36:31.14 AAS
>>70
おお！ありがとうございます！
72: 2012/05/13(日) 12:54:54.09 AAS
ver4.6でグラフをemfにしてファイル出力してるんだけど
powerpoint上でグループ解除するとグラフには関係ない謎の透明で巨大な枠が
出現して、作業効率が下がるんだけどこれを解決できないだろうか
ver4.2ではこんなことはなかったんだが
仕様でも変わった・・・？
73(2): 2012/05/18(金) 02:30:06.78 AAS
↓みたいなfor文において
plot for [2:3] "hoge.dat" using 1:i
ここでusingにabsなどの処理をしたいときはどうやればいいですか？
using 1:abs(i) としたら案の定駄目でした
74: [age] 2012/05/28(月) 22:16:09.31 AAS
gnuplot> q=1.60e-19
gnuplot> print q**16
1.84467440737095e-301
gnuplot> print q**20
0.0
上記のように、q=1.60e-19の累乗が大きい値になると、計算結果が0.0になってしまいます。
どうすれば、正確な値がでてきますか？
(次に続く)
75: 2012/05/29(火) 05:41:57.98 AAS
オーバーフロー？
76: 2012/05/29(火) 14:40:59.42 AAS
この場合アンダーフローじゃね？
77: 2012/05/29(火) 23:45:55.64 AAS
そっか、マイナスついていたか。
確かにアンダーフローだな。
78: 2012/05/30(水) 00:27:58.73 AAS
gnuplotの問題ってわけじゃないね。
何事にも限りがあるってことだ
doubleが64ビットより多い処理系でビルドすればましにはなるかな。
いずれにしろ限りはあるが
79: 2012/06/02(土) 15:38:37.78 AAS
set xtics (pi)にした時、横軸に3.14，6.28，…じゃなくて
π，2π，…ってするにはどうしたらいいの？
80: 2012/06/03(日) 05:14:39.77 AAS
set xtics (\
"-2{/Symbol p}" -2*pi, \
"-3{/Symbol p} / 2" -3*pi/2, \
"-{/Symbol p}" -pi, \
"-{/Symbol p} / 2" -pi/2, \
"0" 0, \
"{/Symbol p} / 2" pi/2, \
"{/Symbol p}" pi, \
"3{/Symbol p} / 2" 3*pi/2, \
"2{/Symbol p}" 2*pi \
)
81(1): 2012/06/04(月) 12:18:28.10 AAS
>>73
外部ﾘﾝｸ[html]:takeno.iee.niit.ac.jp
によると，
plot for [i=2:3] 'hoge.dat' using 1:(abs($i))
すべきのようです。
82(1): 2012/06/15(金) 13:18:50.34 AAS
1つのグラフに対して上と下のx軸で別の単位を使うことはできますか?
可視光線の透過スペクトルで下に波長(nm)，上に波数(/cm)を表示したいです．
下軸の波長を100程度の等間隔にして，上に波数(=10^7/波長)の目盛りを付けられませんか
83: 2012/06/15(金) 13:20:14.74 AAS
set x2label
set x2tics nomirror
plot hoge axes x2y1
84: 82 2012/06/16(土) 13:21:48.76 AAS
こういうのです
外部ﾘﾝｸ:christianyates.com
85: 73 2012/06/22(金) 01:40:41.27 AAS
>>81
4.6 patchlevel 0で試したところ (abs($i)) は駄目ですが、 (abs(column(i))) でいけました。
ありがとうございました。
86(5): 2012/07/03(火) 19:49:49.52 AAS
負の数をログスケールにしたいんですが...
set logscale y
set ytics format "-%g"
plot 'hoge' using 1:(-$2)
こんな感じでいけますが，x11などの場合で，
いったん普通にプロットしたあと，キー入力でノーマルスケールと
ログスケールを切り替えたいです。
bind で何とかならないかとも思ったのですが，
plot文を変える必要があるため，特定の用途ならともかく，
汎用的にはやるのはなかなか難しそうです。
理想は
set logscale y -10
とかやると，負数のログスケールになってくれたりするとありがたいんですが...
ちなみにソースを見るとfabsを取っているようで，-10でも10でも同じ結果になります。
87: 2012/07/22(日) 10:00:47.11 AAS
保守
88: 86 2012/08/18(土) 21:27:55.95 AAS
>>86です。
その後ぱったりとレスが途絶えました。
僕は地雷を踏んでしまったのでしょうか？？
89: 2012/08/18(土) 21:49:31.44 AAS
地雷ではなく，単に高度なことを要求して，
誰も答えられないだけだと思われる．
90: 2012/08/20(月) 17:30:11.22 AAS
>>86
ソースいじれ。
もしくは、開発チームに参加しろ。
91: 2012/08/21(火) 00:38:27.78 AAS
定番の、言い出しっぺの法則
92: 2012/08/23(木) 17:43:50.75 AAS
>>86
外部ﾘﾝｸ[html]:takeno.iee.niit.ac.jp
負の数に対数軸を適用したい、ということと、 キー入力でグラフの線形軸と対数軸を切り替えたい (x11)、 という 2 つの質問 (多分) がありました (86)。
93: 86 2012/08/26(日) 16:38:04.67 AAS
>>86です。
竹野先生ありがとうございます。
ちょっと言葉足らずでしたが，やりたいことは，x11ウィンドウで
負数のデータをリニアスケールでプロットした後，キー入力で
（符号を変えた）ログスケールとリニアスケールを交互に切り替えたい，
というものでした。
正数のデータであればｌキーやＬキーで切り替えられますが，
Ｙデータが全部負数の場合はエラーになったりします。
一部が負数の場合は正数部分のみが表示されます。
ちなみに，負数のデータというのは，たとえばあるデバイスを流れる
電流のようなものです。電流方向の定義というか電流計を接続する
向きによって同じものが正になったり負になったりします。
しかもバイポーラトランジスタのGummel plotのようにログスケールで
あらわしたいこともよくあります。
gnuplotソース改変も試みましたが，自分のスキルではちょっと無理そうです。
他には，リニアスケールとログスケールを別ウィンドウでプロットし，
gnuplot自身でなく別のツールで切り替えるというようなこともアイデアと
しては考えています。
94: 2012/08/26(日) 23:51:59.85 AAS
昔、aplacというspiceのパチもんみたいな回路シミュレータでやったな。
spiceでもグラフの切替はできるが、インタラクティブには出来なかった記憶が。
そもそもF77だったし。
95(1): 2012/08/30(木) 21:27:43.82 AA×
96: 86 2012/09/01(土) 23:36:40.90 AAS
>>86です。
>>95さん
いい感じです。変数を使ってbindにplot分を割り当てるのは盲点でした。
自分の目的としてはこれで何とかいけそうです。
ありがとうございました。
97(1): 質問 1/4 2012/09/11(火) 07:49:45.81 AAS
4.6 をインストールしたのだが、上手く動かない。多分4.6固有の問題じゃない気がする。
環境は、Windows7 home 64bit
今までは、4.0 を使っていて、.pltファイルをwgnuplot.exeに関係づけて、、それは正常に動いていた。
4.6は、ディフォルトで .plt に関係付けられているので今までどうりのファイルでいいのだろうと思い、試して見た。
98: 質問 2/4 2012/09/11(火) NY:AN:NY.AN AAS
取りあえず weight06_test.plt として以下のようなファイルを作った。
#cd "D:/USER/TOMURA/Data/体重/"
set term postscript
set title "Body Weight Change (2012/09/01-2013/08/31)"
set xlabel "Date [month/day]"
set xdata time
set timefmt "%Y/%m/%d"
set format x "%m/%d"
set ylabel "Weight [kg]"
set yrange [70:75]
set grid xtics ytics back lw 2
set style line 1 lt 1 lw 1 pt 7 ps 0.5
set output 'weight06_test.ps'
plot 'weight06.txt' using 1:2 with linesp notitle ls 1
これをダブルクリックしても weight06_test.ps は、出来上がらない。 4.0 では、出来上がっていたはず。
99: 質問 3/4 2012/09/11(火) 07:50:48.72 AAS
試しに同じフォルダに作ったバッチファイルに
"C:\Program Files (x86)\gnuplot\bin\wgnuplot.exe"
と書いて実行して weight06_test.plt の中身をコピペで与えると weight_test.ps が出来上がる。
#cd のコメントを削って生かしてやっても .plt ファイルのダブルクリックでは、正常に動かないが、
普通に wgnuplot を起動して cd からの weight06_test.plt の中身をコピペしてやると正常に動く。
(weight06_test.ps が出来上がる)
仕方ないので今は
"C:\Program Files (x86)\gnuplot\bin\wgnuplot.exe" weight06_test.plt
と書いたバッチファイルを作ってそれを実行している。
100: 質問 4/4 2012/09/11(火) 07:51:19.32 AAS
何とか今まで通りに .plt ファイルをダブルクリックしただけで済ませたいのだが、そういう使い方が間違っているのだろうか？
昔、plot してから、 set output あたりをを弄ってやると正常に出力されるようなこと見た覚えがあるのだが、自力で見つけられない。
そういうわけで教えて下さい。(Windows の質問はひょっとするとスレ違いでしょうか？)
101: 2012/09/14(金) 00:52:29.33 AAS
gnuplot で bash の alias のように
長ったらしいコマンドを短く定義したいのですが，
そういうことは gnuplot でできないでしょうか？
特に出力先を変える時の設定が面倒なので…．
102(1): 2012/09/15(土) 21:53:35.94 AAS
今試せないけどマクロ置換は使えないの？
103: 2012/09/16(日) 00:58:09.69 AAS
>>102
まさにその機能です！
ありがとうございます．
gnuplot で関数やらサブルーチンやらエイリアスで
ググっても思った内容がなかなか出て来なかったもので…．
マクロ置換って言うんですね．
助かりました．
104: 2012/09/19(水) 17:45:28.52 AAS
>>97 -100
set term post
set output 'testps.ps'
test
でgnuplot 4.7で試したが，ダブルクリックでファイルができたよ。
105(1): 2012/09/28(金) 07:48:53.82 AAS
いまgnuplot使ってデータ点と関数のフィッティングやってるんですけど、フィッティングの際に各データテント関数の値との相対誤差とか出せませんか？
教えていt抱けると非常に助かります
106(1): 2012/09/29(土) 01:07:45.00 AAS
>>105
日本語でおｋ
107: 2012/10/01(月) 09:34:08.84 AAS
>>106
わかりにくい日本語で済みませんでした。
今(x_n,y_n) n=1〜1000のデータ点を用意して、とある関数（たとえばy=a*x+b）とfitコマンドを使ってgnuplotでフィッティングしています。
その結果、データ点と会うように係数a,bが調整された関数を得ました。
この関数にx_nを代入した値と、おおもとの(x_n,y_n)との間の相対誤差を知りたいのですが、いちいち計算せずにgnuplot上でどうにかできないかという質問です。
108(1): 2012/10/01(月) 11:21:40.72 AAS
x_n:1列目，y_n:2列目のデータが並んでいると仮定して書くけど、
f(x)=a*x+b
fit f(x) なんちゃら using 1:2 via a,b
ここまではできたということだよね？
「相対誤差を知りたい」って相対誤差をプロットしたいってことかな？
相対誤差の定義もいろいろと有ると思うけど，こんなんでどうでしょう？
plot なんちゃら using 1:(($2-f($1))/$2)
109: 2012/10/04(木) 06:45:28.63 AAS
>>108
ありがとうございます！
110(1): 2012/10/09(火) 10:23:10.05 AAS
set formatについて質問です。
例えば、x軸に80 90 100という値があるとして、
その表示を8.0x10^1 9.0x10^1 10x10^1とすることは出来ますか？
仮数部と指数部があるとしたら、指数部を固定したいということになります。
よくある方法で、set format "%2.1l x 10^{%L}"とすると
8.0x10^1 9.0x10^1 1.0x10^2となってしまいます。
111(2): 2012/10/09(火) 22:11:34.14 AAS
>>110
10x10^1 って指数表記にする意味ってあるの？
もし，関数ではなくデータのプロットで
表示している数値の底上げをしたいのであれば，
軸の脇に "x 10" とか書いて，数値に対して
using ($1/10):$2 とかすればいいのでは？
俺のレベルでアドバイスできるのはそれくらい．
112(1): 111 2012/10/09(火) 22:16:12.49 AAS
訂正
× using ($1/10):$2
○ using ($1/10):2
113: 2012/10/09(火) 22:27:48.31 AAS
>>111-112
レスありがとうございます。
残念ながらデータのプロットではなく関数です・・・
一応、関数に10の累乗をかけて整数にしておいて、
目盛りやラベルで誤魔化そうとも思ったのですが、
あまりかっこ良くないなぁと。
114(1): 2012/10/23(火) 19:45:17.38 AAS
ファイルに保存された3次元のデータを描画したいと考えています.
これを普通に描画すると,すべての点が赤で描画されます.
これをzの値で色を変えるにはどうすればよいのでしょうか?
カラーマップを使うと,点ではなく面でベタ塗りされてしまうのです.
私がやりたいのは,こうではなくて,点の色を変えたいのです.
以上,お願いします.
115(1): 2012/10/24(水) 12:16:32.49 AAS
>>114 これでどう？
外部ﾘﾝｸ[html]:gnuplot.sourceforge.net
116: 2012/10/24(水) 13:31:58.71 AAS
>>115
ありがとうございます.
これで何とかなりそうです.
117(1): 2012/10/31(水) 16:16:28.54 AAS
手元に10000行1列の各数字が0 ~ 20範囲のデータがあります。
横軸に 0~5 5~10 10~15 15~20の各範囲、縦軸にそれぞれの範囲内にある数の個数の4本の棒グラフを作成したいのですが、gnuplotで可能でしょうか?
それとも何らかの方法でカウント数をソートしてから読み込むしかないのでしょうか？
118(1): 2012/10/31(水) 21:38:57.02 AAS
それってひとことで言うとヒストグラムのことか？
直接ヒストグラムを描くのは無理だけど，
やりようはあるからあとはググれ。
119: 2012/10/31(水) 21:41:40.79 AAS
>>118
そうです、ヒストグラムです。
やりようがわからないのです。方針だけでいいので教えて頂けませんか。
120: 2012/10/31(水) 22:32:14.34 AAS
つ Google検索 「gnuplot ヒストグラム」
俺だったら直接描く方法がないと分かった時点でカウントするプログラム書いてるけど
121: 2012/11/01(木) 00:57:06.70 AAS
gnuplot って基本グラフ描画ツールなんだよね。
統計とか計算はおまけみたいなもので…。
122: 117 2012/11/01(木) 14:25:34.39 AAS
ぐぐったんですけど、結局良くわからないですね…
gnupotは結構隠しコマンド的なのがあるからできるかもなどと聞いたのですが。
とりあえずカウントするプログラム別に作ることにします。
ありがとうございました。
123: 2012/11/02(金) 23:49:11.73 AAS
ヒストグラムを生成しようとして以下の様に入力しました。
n=20 #number of intervals
max=90. #max value
min=-90. #min value
width=(max-min)/n
hist(x,width)=width*floor(x/width)+width/2.0
set term postscript eps enhanced color
set output "histogram.eps"
set xrange [min:max]
set yrange [0:]
set mxtics 6
set boxwidth width*0.9
set style fill solid 0.5
set tics out nomirror
plot "dec.dat" u (hist($1,width)):(1.0) smooth freq w histeps
この様にhistepsで出力すると綺麗なヒストグラムが得られたのですが、最後の行を
plot "dec.dat" u (hist($1,width)):(1.0):(0.5) with yerrorbars
とすると高さ1のところにプロットが数個現れるだけで上手くエラーバー付きのプロットができません。
どなたか原因が分かる方いませんでしょうか。
124: 2012/11/08(木) 00:08:31.56 AAS
ヒストグラムにエラーバーとは何がしたいかわかんないけど、
うちの環境では普通にエラーバー出るよ。
バージョン古かったりしない？
125(2): 2012/11/08(木) 18:09:34.09 AAS
ｘ軸の名前の位置を調整したいのですがうまく出来ません。例えばx方向に1、y方向に2だけずらしたい時には
set xlabel "x-Name" 1.0,2.0 font "Helvetica,20"
とすれば良いと聞いたのですがエラーが出てしまいます。
使用しているバージョンはgnuplot4.7です。何か方法があったら教えて下さい。
126(1): 2012/11/08(木) 18:28:29.96 AAS
>125
使っているOSや エラーの内容を書かないのは 何かの策略ですか
set xlabel "x-Name" offset 1.0, 2.0 font "..."
と offset を書けばいいのでは
127: 2012/11/08(木) 23:35:53.26 AAS
>>126
お返事ありがとうございます、お陰でできました！
すみません、うっかりしていて書き忘れてしまいました。
Mac OS XでX11を使っています。
offset と書かく必要があったとは。調べたら>>125の表記をしているところが多かったもので…
versionによる違いなのでしょうか。
どちらにしろ無事できたのでよかったです、大変助かりました。
上下前次1-新書関写板覧索設栞歴
あと 828 ﾚｽあります
ｽﾚ情報 赤ﾚｽ抽出 画像ﾚｽ抽出 歴の未読ｽﾚ
ぬこの手 ぬこTOP 0.049s