gnuplot を使おう。　その3

1 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/08/25(木) 12:15:50.61 .net] グラフ描画ソフト gnuplot のスレ。 gnuplot 公式 Web ページ www.gnuplot.info/ 前スレ gnuplot を使おう。　その2 hibari.2ch.net/test/read.cgi/unix/1164719425/ 45 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/01/23(月) 10:21:37.76 .net] 2本に分けて描く 46 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/01/23(月) 13:38:40.04 .net] (abs(1/x)<c)?(1/x):(0/0) 定数cとsampleの値を適当に調節する 47 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/01/25(水) 17:52:13.49 .net] >>45 無限に発散する点と普通に0を通過していく点があるから 2本に分けるにも結構難しかった もっと効率よくできないものか >>46 計算ソフトで出力されたデータをgnuplotに読み込ませて描写しているから 関数表現だと無理だと思う それをプロットに応用できればいいのだが… 計算能力の問題で＞普段0.001程度の値が1まで行って-1まで一気に落ちる という結果になってる。 あくまで例で綺麗に1から-1まで一瞬に移動してくれるわけでもなく 理論上無限だけど、その点は効果が薄いために 他の発散は1から-1まで下がるのに対して0.05から-0.05に下がる点があったり 発散の度合いにも差があって、発散が弱いと普通に0を通っているのか発散しているのか見分けが難しい 48 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/01/25(水) 18:23:29.11 .net] >>47 >計算ソフトで出力されたデータをgnuplotに読み込ませて描写しているから > 関数表現だと無理だと思う Help ternaryをみると次のようなものがあるから，三項演算子とusingで 可能だと思います。 plot 'file' using 1:( $4<0 ? 1/0 : ($2+$3)/2 ) 49 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/01/31(火) 14:14:17.57 .net] sin 関数のグラフで x 軸のメモリにπを入れる方法ってない？ 50 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/01/31(火) 18:13:06.19 .net] >>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 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/01/31(火) 18:15:48.73 .net] >>50 turial.pdf ==> tutorial.pdf 52 名前：49 mailto:sage [2012/02/01(水) 04:48:42.27 .net] >>50 ありがとう！ よく，emf ファイルで出力するからこの方法で大丈夫か 心配だったけど emf も enhance 付ければ出力できるね！ 53 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/02/01(水) 11:24:59.88 .net] >>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 名前：名無しさん＠お腹いっぱい。 [2012/02/10(金) 21:40:03.29 .net] "Skipping data file with no valid points"のエラーってなんなんですか? どんなコードもこれが出てきてしまってグラフィックができないです。 55 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/02/10(金) 23:25:38.50 .net] 英語で書かれているとおりだよ。 無効なデータをスキップしましたってこと。 データのフォーマットがおかしいんじゃないの？ 56 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/02/12(日) 08:06:35.93 .net] >>54 実行できないコードの一つでもさらせば，解答がつくかも 57 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/02/15(水) 01:35:35.43 .net] x データが 48〜120，218〜289 という具合に データの塊が2つにわかれているんだけど， これを省略線(波線2つ)を入れて省略するみたいに 表現できないかな？ 58 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/02/15(水) 07:11:17.68 .net] >>57 データではないけどそういう例はある www.phyast.pitt.edu/~zov1/gnuplot/html/broken.html 59 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/02/15(水) 15:07:57.76 .net] >>58 ありがとう！ 簡単にできるわけじゃないのねｗ． とりあえず，じっくり読み考えながらやってみるわ． 60 名前：名無しさん＠お腹いっぱい。 [2012/03/12(月) 15:18:16.63 .net] gnuplot version 4.6 released! 61 名前：名無しさん＠お腹いっぱい。 [2012/03/12(月) 19:26:06.94 .net] ↑ あららら・・・ インストール 変わっちゃったね？ 62 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/03/13(火) 08:06:35.29 .net] gnuplot 4.6.0がでてます。 sourceforge.net/projects/gnuplot/files/gnuplot/4.6.0/ 63 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/03/13(火) 08:09:50.61 .net] >>62 うああ，かぶったm(__)m。 >>61 Windows版ですね。インストーラーがついてますよね。 でも，従来型のzipもあります。 64 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/04/09(月) 22:25:18.79 .net] 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 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/04/28(土) 04:15:16.94 .net] 4.4以降でterminal windowsとかでemfを吐くと それをパワポ上でグループ解除すると枠線らしきものが かなりのでかさで現れて酷いことになる wxtだと確かに綺麗だけど ビットマップ扱いで拡大するとボヤけたりepsに変換すると汚くなる ボヤけないようにemfで保存すると枠線が・・・となる term windowsだと破線の太さをgnuplot上で太くできないから 仕方なくwxt使ってるんだけど wxtで出力された見た目のままベクター形式で出力できないのだろうか 66 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/04/28(土) 18:53:34.37 .net] > wxtで出力された見た目のままベクター形式で出力できないのだろうか pdfcairo powerpointで読み込む方法があるのか知らんけど 67 名前：名無しさん＠お腹いっぱい。 [2012/05/03(木) 18:01:13.39 .net] ylabelを左じゃなくて、軸の上に表示したい 68 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/05/03(木) 18:37:55.70 .net] >>67 ylabel ではなく，label をグラフの外の座標に 設定すれば良いんじゃないだろうか？ 69 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/05/10(木) 16:20:51.55 .net] gnuplot> set pm3d map gnuplot> splot [:2*pi][:2*pi] sin(x)*cos(y) ってやるとy軸の範囲が0~7になって上に空白が入っちゃうのってどうやって直したらいい？ 70 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/05/10(木) 23:20:52.00 .net] >>69 下限を設定したらうまくいったよ 71 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/05/11(金) 06:36:31.14 .net] >>70 おお！ありがとうございます！ 72 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/05/13(日) 12:54:54.09 .net] ver4.6でグラフをemfにしてファイル出力してるんだけど powerpoint上でグループ解除するとグラフには関係ない謎の透明で巨大な枠が 出現して、作業効率が下がるんだけどこれを解決できないだろうか ver4.2ではこんなことはなかったんだが 仕様でも変わった・・・？ 73 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/05/18(金) 02:30:06.78 .net] ↓みたいなfor文において plot for [2:3] "hoge.dat" using 1:i ここでusingにabsなどの処理をしたいときはどうやればいいですか？ using 1:abs(i) としたら案の定駄目でした 74 名前：名無しさん＠お腹いっぱい。 mailto:age [2012/05/28(月) 22:16:09.31 .net] 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 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/05/29(火) 05:41:57.98 .net] オーバーフロー？ 76 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/05/29(火) 14:40:59.42 .net] この場合アンダーフローじゃね？ 77 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/05/29(火) 23:45:55.64 .net] そっか、マイナスついていたか。 確かにアンダーフローだな。 78 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/05/30(水) 00:27:58.73 .net] gnuplotの問題ってわけじゃないね。 何事にも限りがあるってことだ doubleが64ビットより多い処理系でビルドすればましにはなるかな。 いずれにしろ限りはあるが 79 名前：名無しさん＠お腹いっぱい。 [2012/06/02(土) 15:38:37.78 .net] set xtics (pi)にした時、横軸に3.14，6.28，…じゃなくて π，2π，…ってするにはどうしたらいいの？ 80 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/06/03(日) 05:14:39.77 .net] 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 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/06/04(月) 12:18:28.10 .net] >>73 takeno.iee.niit.ac.jp/~shige/unix/gnuplot/gpinfo-2012.html#info-20120603 によると， plot for [i=2:3] 'hoge.dat' using 1:(abs($i)) すべきのようです。 82 名前：名無しさん＠お腹いっぱい。 [2012/06/15(金) 13:18:50.34 .net] 1つのグラフに対して上と下のx軸で別の単位を使うことはできますか? 可視光線の透過スペクトルで下に波長(nm)，上に波数(/cm)を表示したいです． 下軸の波長を100程度の等間隔にして，上に波数(=10^7/波長)の目盛りを付けられませんか 83 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/06/15(金) 13:20:14.74 .net] set x2label set x2tics nomirror plot hoge axes x2y1 84 名前：82 [2012/06/16(土) 13:21:48.76 .net] こういうのです christianyates.com/category/subject/flex 85 名前：73 mailto:sage [2012/06/22(金) 01:40:41.27 .net] >>81 4.6 patchlevel 0で試したところ (abs($i)) は駄目ですが、 (abs(column(i))) でいけました。 ありがとうございました。 86 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/07/03(火) 19:49:49.52 .net] 負の数をログスケールにしたいんですが... 　　set logscale y 　　set ytics format "-%g" 　　plot 'hoge' using 1:(-$2) こんな感じでいけますが，x11などの場合で， いったん普通にプロットしたあと，キー入力でノーマルスケールと ログスケールを切り替えたいです。 bind で何とかならないかとも思ったのですが， plot文を変える必要があるため，特定の用途ならともかく， 汎用的にはやるのはなかなか難しそうです。 理想は 　　set logscale y -10 とかやると，負数のログスケールになってくれたりするとありがたいんですが... ちなみにソースを見るとfabsを取っているようで，-10でも10でも同じ結果になります。 87 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/07/22(日) 10:00:47.11 .net] 保守 88 名前：86 mailto:sage [2012/08/18(土) 21:27:55.95 .net] >>86です。 その後ぱったりとレスが途絶えました。 僕は地雷を踏んでしまったのでしょうか？？ 89 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/08/18(土) 21:49:31.44 .net] 地雷ではなく，単に高度なことを要求して， 誰も答えられないだけだと思われる． 90 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/08/20(月) 17:30:11.22 .net] >>86 ソースいじれ。 もしくは、開発チームに参加しろ。 91 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/08/21(火) 00:38:27.78 .net] 定番の、言い出しっぺの法則 92 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/08/23(木) 17:43:50.75 .net] >>86 ttp://takeno.iee.niit.ac.jp/~shige/unix/gnuplot/gpinfo-2012.html#info-20120823 負の数に対数軸を適用したい、ということと、 キー入力でグラフの線形軸と対数軸を切り替えたい (x11)、 という 2 つの質問 (多分) がありました (86)。 93 名前：86 mailto:sage [2012/08/26(日) 16:38:04.67 .net] >>86です。 竹野先生ありがとうございます。 ちょっと言葉足らずでしたが，やりたいことは，x11ウィンドウで 負数のデータをリニアスケールでプロットした後，キー入力で （符号を変えた）ログスケールとリニアスケールを交互に切り替えたい， というものでした。 正数のデータであればｌキーやＬキーで切り替えられますが， Ｙデータが全部負数の場合はエラーになったりします。 一部が負数の場合は正数部分のみが表示されます。 ちなみに，負数のデータというのは，たとえばあるデバイスを流れる 電流のようなものです。電流方向の定義というか電流計を接続する 向きによって同じものが正になったり負になったりします。 しかもバイポーラトランジスタのGummel plotのようにログスケールで あらわしたいこともよくあります。 gnuplotソース改変も試みましたが，自分のスキルではちょっと無理そうです。 他には，リニアスケールとログスケールを別ウィンドウでプロットし， gnuplot自身でなく別のツールで切り替えるというようなこともアイデアと しては考えています。 94 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/08/26(日) 23:51:59.85 .net] 昔、aplacというspiceのパチもんみたいな回路シミュレータでやったな。 spiceでもグラフの切替はできるが、インタラクティブには出来なかった記憶が。 そもそもF77だったし。 95 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/08/30(木) 21:27:43.82 .net] file='file1.txt' plot file u 1:2 log=0 bind 'ctrl-l' '\ 　　if (log==1) \ 　　　　unset logscale y; \ 　　　　plot file u1:2; \ 　　　　log=0; \ 　　else \ 　　　　set logscale y ; \ 　　　　plot file u 1:(abs($2)); \ 　　　　log=1' 96 名前：86 mailto:sage [2012/09/01(土) 23:36:40.90 .net] >>86です。 >>95さん いい感じです。変数を使ってbindにplot分を割り当てるのは盲点でした。 自分の目的としてはこれで何とかいけそうです。 ありがとうございました。 97 名前：質問 1/4 mailto:sage [2012/09/11(火) 07:49:45.81 .net] 4.6 をインストールしたのだが、上手く動かない。多分4.6固有の問題じゃない気がする。 　環境は、Windows7 home 64bit 　今までは、4.0 を使っていて、.pltファイルをwgnuplot.exeに関係づけて、、それは正常に動いていた。 　4.6は、ディフォルトで .plt に関係付けられているので今までどうりのファイルでいいのだろうと思い、試して見た。 98 名前：質問 2/4 mailto:sage [2012/09/11(火) NY:AN:NY.AN .net] 　取りあえず 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 mailto:sage [2012/09/11(火) 07:50:48.72 .net] 　試しに同じフォルダに作ったバッチファイルに "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 mailto:sage [2012/09/11(火) 07:51:19.32 .net] 　何とか今まで通りに .plt ファイルをダブルクリックしただけで済ませたいのだが、そういう使い方が間違っているのだろうか？ 　昔、plot してから、 set output あたりをを弄ってやると正常に出力されるようなこと見た覚えがあるのだが、自力で見つけられない。 そういうわけで教えて下さい。(Windows の質問はひょっとするとスレ違いでしょうか？) 101 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/09/14(金) 00:52:29.33 .net] gnuplot で bash の alias のように 長ったらしいコマンドを短く定義したいのですが， そういうことは gnuplot でできないでしょうか？ 特に出力先を変える時の設定が面倒なので…． 102 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/09/15(土) 21:53:35.94 .net] 今試せないけどマクロ置換は使えないの？ 103 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/09/16(日) 00:58:09.69 .net] >>102 まさにその機能です！ ありがとうございます． gnuplot で関数やらサブルーチンやらエイリアスで ググっても思った内容がなかなか出て来なかったもので…． マクロ置換って言うんですね． 助かりました． 104 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/09/19(水) 17:45:28.52 .net] >>97 -100 set term post set output 'testps.ps' test でgnuplot 4.7で試したが，ダブルクリックでファイルができたよ。 105 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/09/28(金) 07:48:53.82 .net] いまgnuplot使ってデータ点と関数のフィッティングやってるんですけど、フィッティングの際に各データテント関数の値との相対誤差とか出せませんか？ 教えていt抱けると非常に助かります 106 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/09/29(土) 01:07:45.00 .net] >>105 日本語でおｋ 107 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/10/01(月) 09:34:08.84 .net] >>106 わかりにくい日本語で済みませんでした。 今(x_n,y_n) n=1〜1000のデータ点を用意して、とある関数（たとえばy=a*x+b）とfitコマンドを使ってgnuplotでフィッティングしています。 その結果、データ点と会うように係数a,bが調整された関数を得ました。 この関数にx_nを代入した値と、おおもとの(x_n,y_n)との間の相対誤差を知りたいのですが、いちいち計算せずにgnuplot上でどうにかできないかという質問です。 108 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/10/01(月) 11:21:40.72 .net] 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 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/10/04(木) 06:45:28.63 .net] >>108 ありがとうございます！ 110 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/10/09(火) 10:23:10.05 .net] 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 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/10/09(火) 22:11:34.14 .net] >>110 10x10^1 って指数表記にする意味ってあるの？ もし，関数ではなくデータのプロットで 表示している数値の底上げをしたいのであれば， 軸の脇に "x 10" とか書いて，数値に対して using ($1/10):$2 とかすればいいのでは？ 俺のレベルでアドバイスできるのはそれくらい． 112 名前：111 mailto:sage [2012/10/09(火) 22:16:12.49 .net] 訂正 × using ($1/10):$2 ○ using ($1/10):2 113 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/10/09(火) 22:27:48.31 .net] >>111-112 レスありがとうございます。 残念ながらデータのプロットではなく関数です・・・ 一応、関数に10の累乗をかけて整数にしておいて、 目盛りやラベルで誤魔化そうとも思ったのですが、 あまりかっこ良くないなぁと。 114 名前：名無しさん＠お腹いっぱい。 [2012/10/23(火) 19:45:17.38 .net] ファイルに保存された3次元のデータを描画したいと考えています. これを普通に描画すると,すべての点が赤で描画されます. これをzの値で色を変えるにはどうすればよいのでしょうか? カラーマップを使うと,点ではなく面でベタ塗りされてしまうのです. 私がやりたいのは,こうではなくて,点の色を変えたいのです. 以上,お願いします. 115 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/10/24(水) 12:16:32.49 .net] >>114 これでどう？ gnuplot.sourceforge.net/demo_4.6/rgb_variable.html 116 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/10/24(水) 13:31:58.71 .net] >>115 ありがとうございます. これで何とかなりそうです. 117 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/10/31(水) 16:16:28.54 .net] 手元に10000行1列の各数字が0 ~ 20範囲のデータがあります。 横軸に 0~5 5~10 10~15 15~20の各範囲、縦軸にそれぞれの範囲内にある数の個数の4本の棒グラフを作成したいのですが、gnuplotで可能でしょうか? それとも何らかの方法でカウント数をソートしてから読み込むしかないのでしょうか？ 118 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/10/31(水) 21:38:57.02 .net] それってひとことで言うとヒストグラムのことか？ 直接ヒストグラムを描くのは無理だけど， やりようはあるからあとはググれ。 119 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/10/31(水) 21:41:40.79 .net] >>118 そうです、ヒストグラムです。 やりようがわからないのです。方針だけでいいので教えて頂けませんか。 120 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/10/31(水) 22:32:14.34 .net] つ Google検索 「gnuplot ヒストグラム」 俺だったら直接描く方法がないと分かった時点でカウントするプログラム書いてるけど 121 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/11/01(木) 00:57:06.70 .net] gnuplot って基本グラフ描画ツールなんだよね。 統計とか計算はおまけみたいなもので…。 122 名前：117 mailto:sage [2012/11/01(木) 14:25:34.39 .net] ぐぐったんですけど、結局良くわからないですね… gnupotは結構隠しコマンド的なのがあるからできるかもなどと聞いたのですが。 とりあえずカウントするプログラム別に作ることにします。 ありがとうございました。 123 名前：名無しさん＠お腹いっぱい。 [2012/11/02(金) 23:49:11.73 .net] ヒストグラムを生成しようとして以下の様に入力しました。 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 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/11/08(木) 00:08:31.56 .net] ヒストグラムにエラーバーとは何がしたいかわかんないけど、 うちの環境では普通にエラーバー出るよ。 バージョン古かったりしない？ 125 名前：名無しさん＠お腹いっぱい。 [2012/11/08(木) 18:09:34.09 .net] ｘ軸の名前の位置を調整したいのですがうまく出来ません。例えばx方向に1、y方向に2だけずらしたい時には set xlabel "x-Name" 1.0,2.0 font "Helvetica,20" とすれば良いと聞いたのですがエラーが出てしまいます。 使用しているバージョンはgnuplot4.7です。何か方法があったら教えて下さい。 126 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/11/08(木) 18:28:29.96 .net] >125 使っているOSや エラーの内容を書かないのは 何かの策略ですか set xlabel "x-Name" offset 1.0, 2.0 font "..." と offset を書けばいいのでは 127 名前：名無しさん＠お腹いっぱい。 [2012/11/08(木) 23:35:53.26 .net] >>126 お返事ありがとうございます、お陰でできました！ すみません、うっかりしていて書き忘れてしまいました。 Mac OS XでX11を使っています。 offset と書かく必要があったとは。調べたら>>125の表記をしているところが多かったもので… versionによる違いなのでしょうか。 どちらにしろ無事できたのでよかったです、大変助かりました。 128 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/11/17(土) 09:29:01.99 .net] WindowsでXming起動してからputtyでLinuxにsshしてXを飛ばしてるんですが、 gnuplotのX画面をAlt+PrintScreenでキャプチャしてPowerpointへコピペするのが 癖になってしまいました。ベクトル画像じゃなくなるけど、サクッとコピペできるのはすごい楽。 ただ唯一の不満がgnuplotのX画面に枠やタイトルバーが付いている点で、 今でJTrimに毎回削ってます（Powerpoint2010のトリミングは使いづらいので）。 もっと効率良いやりかたないでしょうか。ちょっとスレチかもしれないですが、 自分みたいな使い方してる人って結構いそうな気がするので、お尋ねしました。 129 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/11/17(土) 15:14:42.20 .net] >>128 素直に,gnuplotにpngで保存させろ 130 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/11/17(土) 22:03:48.30 .net] wxtならcopy to clipboardボタンがあるけど、機能するかどうかはXサーバ次第かな 131 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/11/17(土) 23:37:51.53 .net] >>129 それが一手間で面倒で・・・ 普段は set term x11 で作業して、よしこれでOKと思ったら set term png してset output "hoge.png"で保存してWindowsに転送して・・・ 132 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/11/18(日) 22:03:44.23 .net] 作業ディレクトリを Samba ディレクトリにしてしまえばいい 133 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/12/08(土) 20:38:07.16 .net] リリースノートのファイル見てたら，インタラクティブなターミナルとしては 前はX11がオススメだったけど，今はwxtとかqtがナイスだよ，って書いてある。 自分はwxtとかqtとか言われてもナニソレなんだけど，X11にくらべて何が良くなるの？ ちょっとぐぐったら表示がキレイになることはわかったけど，インタラクティブな辺で いいこと（たとえばレンジの指定がインタラクティブにできるとか）があれば， がんばってインストールしてみようかなと思う。 134 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/12/08(土) 22:03:31.36 .net] どんな化石みたいな環境使ってるのか知らないが、たいがい人に聞くより自分で試した方が早いぞ > レンジの指定がインタラクティブにできる マウスホイールを使ってできる でも最近のバージョンなら多分wxt以外でもできると思う 135 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/12/09(日) 12:54:55.24 .net] 最近の Linux ディストリビューションの gnuplot が 採用しているのは wxt なのでは？ その状態で x11 にすると萎える。 136 名前：133 mailto:sage [2012/12/09(日) 17:25:46.98 .net] 環境は仕事で使ってるSoralisの古いやつ。 ライブラリとかも古いままでインストールが大変だからまず聞いてみた。 やる価値あるんなら頑張ってみようかと。 137 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/12/09(日) 20:29:01.42 .net] 近くにWindowsPCくらいあるだろうから、そこで試してみればいいんじゃね？ zip版ならインストールする必要もないし 138 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/12/10(月) 17:32:54.07 .net] Solarisとか、いつまで死んだ子の歳を数えてるんだ… 139 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/12/23(日) 15:32:06.86 .net] 複雑な事しない実験系なら N88BASICとか普通に現役だからなぁ・・・ Solarisみたいな”新しい”環境はまだまだ消えないだろう 140 名前：名無しさん＠お腹いっぱい。 [2013/01/10(木) 19:01:23.74 .net] 区間[0:2*pi]でsin(x), 区間[-2*pi:0]でcos(x), 区間[-2*pi:2*pi]でtan(x)を描画するには, どうすればいいですか？ カンマ区切りやreplotでは範囲が指定できないので困っています. 141 名前：名無しさん＠お腹いっぱい。 mailto:sage [2013/01/10(木) 19:36:57.11 .net] こういうことかな？ ＞plot ( 0<x && x<2*pi ? sin(x) : 1/0) w l,\ ＞ ( -2*pi<x && x<0 ? cos(x) : 1/0) w l,\ ＞ ( -2*pi<x && x<2*pi ? tan(x) : 1/0) w l 142 名前：名無しさん＠お腹いっぱい。 mailto:sage [2013/01/10(木) 19:58:19.98 .net] >>141 できました. 三項演算子とは盲点でした. ありがとうございます. ところで, 1/0はなんですか? 143 名前：名無しさん＠お腹いっぱい。 mailto:sage [2013/01/10(木) 20:47:16.29 .net] help嫁 144 名前：名無しさん＠お腹いっぱい。 mailto:sage [2013/01/10(木) 20:55:40.24 .net] pm3d mapで描画したカラーマップ上に、 境界線をわかりやすくするためにx軸に垂直な直線(例 x=1)を入れたいです 1. noheadのarrowで書こうとしたが、カラーマップの下に入って見えない 2. 以下の通り媒介変数表示を使って書こうとしたが、表示されない set parametric splot 1,u,0　　　　←3次元上には書けてる set pm3d replot　　　　　　 ←同上 set pm3d map replot　　　　　　 ←グラフ上に何も表示されない pm3d map上にx=1等の直線を描画するために 何かいい方法はありませんでしょうか、知恵をおかしくださいｍｍ 145 名前：名無しさん＠お腹いっぱい。 mailto:sage [2013/01/10(木) 22:10:41.01 .net] 描画の順番を入れ替える.

---

---

[ 続きを読む ] / [ 携帯版 ]

read.cgi ver5.27 [feat.BBS2 +1.6] / e.0.2 (02/09/03) / eucaly.net products.
担当:undef