データ解析ツールoctaveを語ろう Part 2

1 名前：名無しさん＠お腹いっぱい。 mailto:sage [2009/10/28(水) 22:13:01 .net] 知る人ぞ知る強力なデータ解析ツール octave ... unixを科学研究で使っている香具師には非常に心強い味方のはずだ。 matlab互換を主張しながら、微妙にバージョンアップを続ける octaveについて語ろう。 そのPart 2です 前スレ データ解析ツールoctaveを語ろう pc12.2ch.net/test/read.cgi/unix/1086016412/ 23 名前：名無しさん＠お腹いっぱい。 mailto:sage [2009/12/16(水) 22:13:45 .net] matlabのkmeans関数ってoctaveには（少なくとも同じ名前では）無いようですが、 別名で同じ機能の関数ってあるのでしょうか？ それともこれは自前で作成するしかないですか？ 24 名前：octave for windows メモの筆者 mailto:sage [2009/12/21(月) 18:11:04 .net] >>23 私はmatlabユーザーではないのでよく分かりませんが， そもそもkmeans関数はどのような機能の関数なのでしょうか？ また，mファイルで書かれている関数で，ソースが手に入るのならば， octaveで動かせる可能性があります。 mexでもソースが手に入れば，octaveでmexを作成しなおせば， 動く可能性があります。 25 名前：名無しさん＠お腹いっぱい。 mailto:sage [2010/01/06(水) 09:48:14 .net] >>19 >>21 dynare の version が 4.1.0にアップデートされました。 Octave 3.2.3 で動きます。 26 名前：名無しさん＠お腹いっぱい。 [2010/01/15(金) 01:17:52 .net] 23 さん Octave で書かれた K-mean プログラムについて． 横から失礼します. k-means のプログラムを Octave で書きましたので，参考になれば と思いお伝えします． ・133.5.18.167/rinkou/octave/okmeans.html 　すべて Octave で書いたもの． ・133.5.18.167/rinkou/octave/kmeans.html 　OpenCV の kmeans 関数を呼び出すようにしているもの 一応、動作確認しています。バグ残っているかも知れません． （すぐには対応できないと思います）．ご参考になるかと思いお伝えしました． 27 名前：：octave for windows メモの筆者 mailto:sage [2010/01/19(火) 18:54:33 .net] Jaroslav Hajek　さんによるとoctave 3.2以上であればfsolveを使って 非線形カーブフィットができるようです。 例 b0 = 3; a0 = 0.2; x = 0:.5:5; noise = 1e-5 * sin (100*x); y = exp (-a0*x) + b0 + noise; c_opt = [a0, b0]; tol = 1e-5; [c, fval, info, output] = fsolve (@(c) (exp(-c(1)*x) + c(2) - y), [0, 0]); Nabble の　アーカイブ ttp://old.nabble.com/Re:-i-just-want-to-use-nonlinear-curve-fitting-function-,-what-should--i-do--p27206660.html 28 名前：：octave for windows メモの筆者 mailto:sage [2010/01/20(水) 08:52:38 .net] octave 3.2.4のRelease Candidateがでてます。 もうすぐ，octave 3.2.4がでると思います。 おそらく，3.2シリーズは3.2.4で終わりで，次はoctave 3.4がでると思います。 ビルドシステムの大幅な変更があります (libtoolを使ったものになります。)。 一般の方にはあまり関係がないと思いますが，自分でビルドされる方は maintainersのMLをご覧になって情報を収集されることをお勧めします。 思いますが， 29 名前：octave for windows メモの筆者 mailto:sage [2010/01/24(日) 12:37:21 .net] octave 3.2.4のソースがでました。（unofficialだそうです。） ttp://old.nabble.com/3.2.4-release-to27272214.html をごらんください。 30 名前：octave for windows メモの筆者 mailto:sage [2010/01/30(土) 09:07:05 .net] Octave公式Webにはまだのっていないようですが，octave-3.2.4 のソースがアップロードされています。 ttp://old.nabble.com/Octave-Version-3.2.4-Released-to27353121.html#a27353121 おそらく，3.2シリーズでは最後のリリースとなると思います。 3.2.3で問題となったバグがかなりフィックスされているようです。 各種バイナリもそのうち整備されると思います。　 31 名前：：octave for windows メモの筆者 mailto:sage [2010/02/03(水) 17:40:31 .net] Cygwin 版のOctave 3.2.4のバイナリが配布されています。 32 名前：：octave for windows メモの筆者 mailto:sage [2010/02/18(木) 13:58:35 .net] Octave のMLにGNU ProjectのリーダーのRichard Stallman氏が登場して， Mathworks社のライセンス（Matlabが関係すると思います）についての話題が 提供され,現在，精力的に議論が行われています。 法律的なこともかかわっているので正直言って理解できないことが多いのですが ライセンスについては，敏感な立場の方がおられると思いNabbleのアーカイブ のURLを書いておきます。 ttp://old.nabble.com/Mathworks-hosted-GPL'd-software-to27550496.html をご覧ください。 33 名前：：：octave for windows [2010/02/23(火) 12:35:16 .net] 計算のたびに以下のように octave-3.2.3.exe:16:C:\Octave\3.2.3_gcc-4.4.0\bin が表示され，非常に見難いのですが 非表示にする方法はありませんか？ > x/y ans = 0.50000 octave-3.2.3.exe:14:C:\Octave\3.2.3_gcc-4.4.0\bin > x^y ans = 16 octave-3.2.3.exe:15:C:\Octave\3.2.3_gcc-4.4.0\bin > x**y ans = 16 octave-3.2.3.exe:16:C:\Octave\3.2.3_gcc-4.4.0\bin > sin(pi/2) ans = 1 octave-3.2.3.exe:17:C:\Octave\3.2.3_gcc-4.4.0\bin 34 名前：：octave for windows メモの筆者 mailto:sage [2010/02/23(火) 18:09:54 .net] >>33 ttp://old.nabble.com/Re:-Really-Long-Command-Prompt-p27222063.html 35 名前：octave for windows メモの筆者 mailto:sage [2010/03/26(金) 08:55:03 .net] Windows 版Octave-3.2.4がでました ttp://sourceforge.net/projects/octave/files/ からどうぞ。 36 名前：octave for windows メモの筆者 mailto:sage [2010/03/26(金) 18:01:51 .net] >>35 2点不具合を発見しました。 ttp://old.nabble.com/Re:-Octave-3.2.4-mingw32-available-p28037642.html ttp://old.nabble.com/flicking-problem-again-Octave-3.2.4-mingw32-td28038688.html 対処はできそうです。 37 名前：octave for windows メモの筆者 mailto:sage [2010/03/26(金) 19:22:22 .net] たびたびすみません。<m(__)m> >>36　のについては ttp://old.nabble.com/flicking-problem-again-Octave-3.2.4-mingw32-td28038688.html なんか何もしなくても，うまくいきだしました。^^; 38 名前：名無しさん＠お腹いっぱい。 mailto:age [2010/04/06(火) 19:30:51 .net] mingwの3.2.4を入れてみたけど wxtターミナルで、plotしたあとに xlabel('hoge'); drawnow; とやっても反映されないのね plotコマンドを再実行すると反映されるので xlabel, ylabelで定義してから plotする流儀に変わったのかな 39 名前：名無しさん＠お腹いっぱい。 mailto:sage [2010/04/06(火) 20:05:32 .net] mingwの3.2.4を入れて気が付いたのが、binフォルダのblas.dll 前のスレ(2chsearch.jp/unix/1086016412/?first=1)の969-970で ttp://machnum.net/base/index.php/articles-computational-physics/37-computational-physics/47-assembler-code-and-mingw ここを見ながらmingw上でgotoblasをコンパイルしてblasのdllを作りました。 またcygwin上でもblasのdllを作りました。 そのblas.dllを、3.2.4デフォルトのblas.dllと入れ替えて動くか試すと、なんと！動きました。 環境変数に GOTO_NUM_THREADS があり、2がセットしてあるのを確認した上で n=2000; A=randn(n); B=randn(n);tic; C=A*B; t=toc, MFLOPS=2*n^3/t*1e-6 これでベンチをとってみると Octave 3.2.4 default blas.dll t = 4.2031 MFLOPS = 3806.7 GotoBLAS mingw t = 2.1719 MFLOPS = 7366.9 GotoBLAS cygwin t = 1.2344 MFLOPS = 1.2962e+004 PCは2.1GHzのCore2Duoで理論演算能力は2.1*4*2＝16.8GFLOPSなので、cygwin版のGotoBlasと入れ替えれば近い線までいきそう。 mingwの場合はマルチコア未対応なのでしょう。 最近発売された6CoreのCore-i7を使えば、3.3*4*6=79.2GFLOPS近くの性能が出るかもしれません。 40 名前：名無しさん＠お腹いっぱい。 mailto:sage [2010/04/08(木) 23:17:13 .net] >38 3.2.4はgnuplotとの連携のとこバグってると思う x = [0, 0.2, 0.4, 0.6, 0.8, 1.0, 1.2, 1.4, 1.6 ,2.0] plot( x, x.*x ) とかもっともシンプルに打っても動きがヘン 41 名前：octave for windows メモの筆者 mailto:sage [2010/04/09(金) 15:38:24 .net] >>38,>>40 たぶん，oct2matというoctave-forgeのパッケージがぶつかっているせいだと 思います。 その他，Octave 本家のMLで話題になったことの対処法を拙ページに書いて おりますのでよろしければ参考にしてください。 ttp://www.tatsuromatsuoka.com/octave/jpn/OctaveWinMemo.html#100403 42 名前：octave for windows メモの筆者 mailto:sage [2010/04/09(金) 15:56:46 .net] >>39 GotoBLAS cygwin t = 1.2344 MFLOPS = 1.2962e+004 これはすごいですね。残念ながらいまだに HT Pentium なので恩恵にはあずか れませんが，導入したらいつかやってみたいですね。 >>mingwの場合はマルチコア未対応なのでしょう。 はい，そうです。SSE3でもなかったような気がします。 ALTASの場合は，cygwinコンパイラでdllを作ると動きません。 また，cygwinで配布しているコンパイラでを使うとcygwinでしか使えない atlasになります。したがって，gcc-4.3.4あたりを自力でcygwin上でビルドし それでスタティックライブラリを作成し，そのスタティックライブラリを mingwのgccでリンクしてdllを作るという技がいります。 pthread-win32も改良されて，HTの場合は速くなりませんが，mingwで動く multi-threaded ATLASの作成には個人的には成功しています。 43 名前：octave for windows メモの筆者 mailto:sage [2010/04/09(金) 15:58:21 .net] >>42 からの続きです 使うマシンがないので公開はしていませんが， mingwでatlasを作成するためだけにCygwin上でビルドした， 特別なGCCとどうやってスタティックライブラリから dllを作るかは書いたほうがいいですかね。 あまり暇がないので，GCCの公開ぐらいならやりますが。 また，Benjmaminさんはどうもうまくmingw上のAtlasをハンドリングできて いないようです。 彼は本当に高速な計算にはLinuxをお使いのようなので， そこまで期待するのはちょっときついかもしれません。 今回octave-forgeのサポートがかなり増えていますので，相当なご尽力を なさっていますのであまり多くを望むのは難しいでしょう。 44 名前：40 mailto:sage [2010/04/09(金) 18:02:17 .net] >41 おお、サンクス。まともにプロットするようになった。 45 名前：名無しさん＠お腹いっぱい。 mailto:sage [2010/04/10(土) 09:44:04 .net] そういえば、wxtの場合なのか現バージョンのoctaveの場合なのかわかりませんが 補助目盛(minor tics)がデフォルトでは表示されません リニアプロットでは気が付かないかもしれませんがloglogだと気付くでしょう そんな場合は set (gca, 'xminortick', 'on', 'yminortick', 'on') とやればいいので、octavercにでも書き込んでおけば問題ないと思ってたのですが warning: mark_as_command is obsolete and will be removed from a future version of Octave という警告が出ますな。正式な手続きがわからん 46 名前：octave for windows メモの筆者 mailto:sage [2010/04/10(土) 13:33:31 .net] >>45 >補助目盛(minor tics)がデフォルトでは表示されません たしかにそうですね。 x=1:100; loglog(x,x) で確認しました。 mingw版でもcygwin版でも（ともに3.2.4）ともにご指摘のようななので そういう仕様なのかもしれませんね。 一度本家MLで聞いてみます。 > warning: mark_as_command is obsolete and will be removed from a future version of Octave たしかにこれ気になっているんですが，octave --norc でstartupファイル を無視するとでないので Octave\3.2.4_gcc-4.4.0\share\octave\3.2.4\m\startup\octavercで 使っている関数がどこかで使っているのでしょう。 パット見た限りわかりませんので気が向いたら，原因探りますが まあ，やらないかな。 いまのとこ実害がないし，mark_as_commandが なくなったらwarningでなくてerrorになるからそのときはBenjaminさん も対応されるでしょう。 47 名前：octave for windows メモの筆者 mailto:sage [2010/04/10(土) 13:56:50 .net] >>45 >>46 >Octave\3.2.4_gcc-4.4.0\share\octave\3.2.4\m\startup\octavercで 使っている関数がどこかで使っているのでしょう。 のではないようですね。 octave --norcでたちあげて pkg load jhandles やると同じメッセージがでます。一度でたら二度目はでないので octave-forgeのどれかをloadするときでるのでしょう。 ちなみにjhandleはデフォールトではloadされない仕様なのでこれは原因 ではないでしょう。 対処しようとするとおそらくpkg周りをいじらないといけないでしょうね。 いまは，我慢して無視するのが現実的かと思います。 どこかの 48 名前：名無しさん＠お腹いっぱい。 mailto:sage [2010/04/10(土) 15:35:16 .net] 警告の原因は、これ（↓）じゃないんですか。了解です set (gca, 'xminortick', 'on', 'yminortick', 'on') 49 名前：octave for windows メモの筆者 mailto:sage [2010/04/11(日) 12:18:20 .net] >>48 はい set (gca, 'xminortick', 'on', 'yminortick', 'on') は警告の原因ではないようです。 mark_as_commandのコードをみましたが，3.2ではdeprecatedというフォルダに はいっていて実質何もしていません。mark_as_commandが実行されたときに エラーになるとコードが動かなくなるのを防ぐために存在しているようです。 というわけで，多分，octave-forgeのパッケージが将来的に対応すべき 問題なのでしょう。現時点ではやはり無視してつかうのでいいでしょう。 50 名前：octave for windows メモの筆者 mailto:sage [2010/04/17(土) 17:31:46 .net] >>45 >>46 >>47 >>48 >>49 warning: mark_as_command is obsolete and will be removed from a future version of Octave の消し方わかりました。 \Octave\3.2.4_gcc-4.4.0\share\octave\3.2.4\m\startup のoctaverc で pkg ("load", "auto"); という文があるので warning ("off","Octave:deprecated-function"); pkg ("load", "auto"); warning ("on","Octave:deprecated-function"); のようにすれば消すことができます。 51 名前：名無しさん＠お腹いっぱい。 mailto:sage [2010/04/24(土) 09:29:28 .net] >>50 了解です それとは関係ありませんが 前スレは、ここで安定して見れるっぽいです ttp://ranran2.net/app/2ch/unix/1086016412 52 名前：octave for windows メモの筆者 mailto:sage [2010/04/26(月) 12:23:21 .net] >>51 おお！　ありがとうございます<m(__)m> 53 名前：名無しさん＠お腹いっぱい。 [2010/04/26(月) 19:12:25 .net] CSVファイルを読み込むためのdlmread関数を使うために pkg load io とすると package io is not installed と表示されるのですが、 ioパッケージって標準では入ってない物なのですか？ 54 名前：octave for windows メモの筆者 mailto:sage [2010/04/27(火) 18:05:28 .net] >>53 octaveのバージョンはなんでしょう。また，platformはなんでしょう。 3.2.4ならdlmreadは，io パッケージでなくoctave本体に入って入ってますが。 また，octave-forgeのパッケージが何がはいるかはdistributorにもよるので 具体的な情報がないと誰もこたえられないと思います。 55 名前：octave for windows メモの筆者 mailto:sage [2010/04/27(火) 18:12:08 .net] >>53 書き忘れましたが，csvならばcsvreadも使えます。dmlreadで セパレータ指定してもいいと思います。 56 名前：名無しさん＠お腹いっぱい。 mailto:sage [2010/05/19(水) 19:54:09 .net] 今さっき気がつきましたが pkg load jhandles でGNUTERMに何が設定してあっても、jhandleでグラフ描画するんですね wxtはなんかグラフがはみ出すし、オーバーヘッドが大きいような気がしてました 57 名前：octave for windows メモの筆者 mailto:sage [2010/05/20(木) 19:23:11 .net] >>56 そうですね。3Dとかでmeshを細かく刻むとjhandlesじゃないと遅いですよね。 ただ，ginputをやろうとするとエラーがでたりします^^;。 3.4がでてFLTK backendが，まともになればjhandleもいらないかなとは 思います。 ただ，3.4リリースのための開発版3.3はMinGWではまだ不安定です。 gnulibでportabilityをあげようとしての結果なのだと思っています。 3.4のリリースはいつごろになるのかな。 58 名前：名無しさん＠お腹いっぱい。 mailto:sage [2010/05/20(木) 22:44:31 .net] >>57 使ってみたらprintがだめっぽいですね gif, png, jpgはだめ MSVC版ではいんちきくさいながらも出来たのに それとAuto Redrawしないので表示が消えちゃいます drawnowすればいいだけですけど 実装があんまりよくないです でもまあ、表示速度は段違いなので 経過を表示させる場合はいいかもです 59 名前：octave for windows メモの筆者 mailto:sage [2010/05/21(金) 11:30:57 .net] >>58 そうですね。まあ，MSVC版ではoctaveのメンテナーとjhandleのメンテナーが 同じ人でしたからね。 MinGWでSun Javaを扱うのが難しいなどで時間もかかったようですし， FLTK Backendがまともになれば，jhandlesの意味も薄れるので力がそれ 実装がMSVC版ほど力が入っていないのはある程度しょうがないと思います。 pipeがwindowsでは遅いのでgnuplotが遅いのはあるていど仕方ないですね。 octaveからだとgnuplotはwindows NativeよりCygwinの方が大分速いですね。 私は，本当につかうグラフはデータ書き出してgnuplotでやります。 octaveからだとgnuplotの機能をフルには引き出せないので octaveのグラフィックスはあくまでも経過表示用にしか使っていません。 60 名前：名無しさん＠お腹いっぱい。 mailto:sage [2010/05/21(金) 13:57:11 .net] >>59 今出来るかどうか知りませんが、昔はこんな事してました gset term postscript color gset output 'hogehoge.ps' plot (x1,y1) gset term X11 単にprint文を知らなかっただけですけど 61 名前：octave for windows メモの筆者 mailto:sage [2010/05/21(金) 16:34:10 .net] >>60 今はできません。それらの関数は3.2からはすべて削除されました。 ただ，octave-forgeにgnuplotパッケージが入ったのでこれをつかうと いう手もあるかもしれません。 MinGW-3.2.4には，gnuplotパッケージ入っているのでそれを使うという手も ありそうです。 ttp://octave.sourceforge.net/functions_by_package.php のgnuplotのreferenceはあまり詳しい記述がないので， これでどんな関数があるか調べてhelp g_plotなどどする必要があります。 3.4ではgnuplotはまだbackendとして残るとおもいますが，将来的には backendはfltkのみなり，gnuplotはoctave-forgeになる流れではと思います。 62 名前：octave for windows メモの筆者 mailto:sage [2010/05/26(水) 19:40:02 .net] 今日，MinGW上でgcc-4.5.0でdevelopment branchのoctaveビルドしました。 先月ビルドしたときは，極めて不安定ですぐseg-faultを起こしていましたが 今日は安定して動いています。 fltk　バックエンドも軸がちゃんと表示できるようになっていましたし， 連続プロットもgnuplotに比べかなり高速です。 ただしマウスで3Dぐりグリができなさそうなのが残念ですが。 しかし，MinGWはdevelopment branchのなかで一番対応が遅れていたのですが これでなんとかなるめどがついた気がします。 あまり遠くない時期にoctave3.4が出てくるような気がします。 63 名前：名無しさん＠お腹いっぱい。 mailto:sage [2010/06/26(土) 11:09:20 .net] Core2DuoのCULV(低電圧)版のSU2300(1.2GHz, 1MBキャッシュ)を積んだノートを買ったので テストしてみました(OSはXP Home) ちなみにCore2以降のCPUに載ってるSSE3は 64bit倍精度浮動小数点演算を1クロックに4個実行できるSIMDユニットで このSSE3が有効な場合の理論FLOPSをを計算すると 　1.2G(Hz)*4(FLOPS/Clock)*2(Core)＝9.6GFLOPS となります。AtomなどにもSSE3のファンクションは載ってますが それを実現するハードウェアはないため、速くはなりません。 GotoBLAS2は新たに1.13をダウンロードしてCygwin1.7上でmakeします makeオプションは www.db.is.kyushu-u.ac.jp/rinkou/cygwin/cblas.html ここのをそのまま使います makeが終わると、exports/にlibgoto2.dllが出来上がってますので これをMingw版Octave 3.2.4のbinのblas.dllと置き換えればOKです その前に環境変数GOTO_NUM_THREADSに2をセットしておき cygwinのbinから、cygwin1.dll、cyggcc_s-1.dll、cyggfortran-3.dllをoctaveのbinにコピーしておきます ベンチマークはいつものを使います n=2000; A=randn(n); B=randn(n);tic; C=A*B; t=toc, MFLOPS=2*n^3/t*1e-6 デフォルト(たぶんatlas.dllを呼んでいる)のblas.dllだと 2.07 GFLOPS GotoBLAS2だと 5.72 GFLOPS で倍以上に速度が上がりました GOTO_NUM_THREADSに1をセットすると 3.24 GFLOPS なのでデュアルコア自体は効いてるようです 64 名前：octave for windows メモの筆者 mailto:sage [2010/07/26(月) 18:27:24 .net] >>62 最近のdevelopmentブランチをビルドしてみました(MinGW，gcc)。 fltkでepsを作成できるようになっており，またマウスによる3dの回転も できるようになっています。 fltkもだいぶ実用に耐えるようになってきた感じです。 なにしろ高速なのがいいです。やっぱりOpenGLをつかっているからでしょう。 pkg installにもバグがあったようですが，それも解消されているし， developmentブランチへの乗り換えを画策しています。 65 名前：名無しさん＠お腹いっぱい。 mailto:sage [2010/08/14(土) 01:19:11 .net] 数学セミナーの9月号で数学を発展させるコンピュータソフトという特集が組まれております www.nippyo.co.jp/magazine/maga_susemi.html そこでは、Rは取り上げられてますが MATLAB/Octaveはさっぱりですけど、 その特集とは別に 高速化・最適化のためのBLAS入門／藤澤克樹　50 と中央大学でSDPAの研究をしている藤澤先生が MATLABとBLASを使った計算について説明してます やってる事は >>63 とほぼ同じ（行列の１辺が2000ではなくて4000）スクリプトで 4コアCPUを2個使って、スレッド数を1個から8個まで増やすとどうなるか という事について書かれてます。 66 名前：名無しさん＠お腹いっぱい。 mailto:sage [2010/09/14(火) 21:48:37 .net] 既知の問題であれば申し訳ないのですが、ちょっと困っています。 以下のようにして、 R = 10.0; L = 0.1e-3; Rc = 100e3; Cc = 0.1e-9; Ca = 0.01e-9; Z = tf([Cc*Rc 1],[Ca*Cc*Rc Cc+Ca 0]); G = tf([1],[L R]); sys = sysmult(G,Z); sysout(sys,"tf"); とすると、本来numはsの一次式になるはずがそうなりません。 Cc、Caを数桁大きくすると正しく計算できるようです。 お手数ですが、この問題の回避策があれば教えていただけないでしょうか？ 67 名前：66 mailto:sage [2010/09/14(火) 21:52:39 .net] すみません、書き忘れました。 Win版 3.2.4です。 68 名前：名無しさん＠お腹いっぱい。 [2010/10/12(火) 03:41:06 .net] 保守 69 名前：octave for windows メモの筆者 mailto:sage [2010/10/16(土) 17:54:10 .net] >>51 現在　前スレ見れるところ ttp://logsoku.com/thread/pc8.2ch.net/unix/1086016412/ 70 名前：名無しさん＠お腹いっぱい。 [2010/11/13(土) 10:00:19 .net] OctaveやFreeMATではMATLABのstemにあたるマッチ棒のようなプロットはどうやるんですか？ 71 名前：名無しさん＠お腹いっぱい。 [2010/11/14(日) 00:32:04 .net] cvytvftvyghgbvyvybuhjkkvghghvgvgfvgfcjnhjgjhhuihuyfgvcg 72 名前：octave for windows メモの筆者 mailto:sage [2010/11/15(月) 16:46:45 .net] >>70 octave manual 15.1.1 Two-Dimensional Plotsより Function File: h = stem (x, y, linespec) ― Function File: h = stem (..., "filled") 　： x = 1:10; y = ones (1, length (x))*2.*x; stem (x, y, "b"); もっと詳しくはマニュアルを 73 名前：octave for windows メモの筆者 mailto:sage [2010/11/24(水) 08:02:34 .net] 前スレでも，>>63 にもあるように，ちょくちょく話題になるGotoBLAS2 に ついてライセンスが変わったとoctave　ML 情報が入っています。 octave.1599824.n4.nabble.com/GotoBLAS2-license-td3055645.html 参考までに 74 名前：名無しさん＠お腹いっぱい。 mailto:sage [2010/11/24(水) 09:50:17 .net] >>73 >GotoBLAS2 has been released by the Texas Advanced Computing Center >as open source software under the BSD license. 素晴らしい。 従来は独特のライセンスのせいでユーザが自分で組み込む必要がありましたが、 今後は公開ソフトウェアに最初から組み込めるわけですね。 GotoBLAS2-1.13で打ち止めへ　(2010-04-09 13:29:50) blog.goo.ne.jp/nakatamaho/e/61d0cfbe1222364a3ef70505dce4f463 後藤和茂さんがヘッドハントされてからバージョンアップが止まってますが、 これも改善するかもしれません。 75 名前：octave for windows メモの筆者 mailto:sage [2010/11/24(水) 19:04:09 .net] >>73 Goto BLASの開発は，終了したようです。 blog.goo.ne.jp/nakatamaho/e/61d0cfbe1222364a3ef70505dce4f463 残念ですが。 76 名前：octave for windows メモの筆者 mailto:sage [2010/11/24(水) 19:06:02 .net] >>75 すみません。 74をみていませんでした<m(__)m>。 77 名前：octave for windows メモの筆者 mailto:sage [2010/11/25(木) 14:49:21 .net] >>74 >後藤和茂さんがヘッドハントされてからバージョンアップが止まってますが、 >これも改善するかもしれません。 そうですね。 天才ハッカー（いい意味でのハッカーです。）後藤さんのあとを 続くような方が現れてプロジェクトとして 立ち上がってくれればいいのですが。 78 名前：名無しさん＠お腹いっぱい。 mailto:sage [2010/12/21(火) 14:14:37 .net] scilabとの違いが良くわからない 79 名前：名無しさん＠お腹いっぱい。 mailto:sage [2010/12/21(火) 15:21:18 .net] >>78 Scilabの使いにくい点 ・実行 　MATLAB/Octaveは、拡張子.mを省いたmファイルをコマンドラインで打ち込めばいい 　Scilabは、コマンドラインでexec('hoge.sce')とやらないといけない ・関数定義 　MATLAB/Octaveは、パスの通ってるディレクトリかカレントディレクトリにあれば 　　勝手に読み込んでくれる。 　Scilabは、明示的にgetf('hoge.sci')と読み込まないといけない。 　※Scilab5では途中からgetfではなくexec('hoge.sci')になった ・コメント行 　MATLABは%、Octaveは%と# 　Scilabは// ・その他、ScilbはMATLABから見て方言が多いし、Octave Forgeのような豊富な関数を使えない。 利点 ・グラフがOctaveよりはきれい ・Serial-0.2.tgzという非オフィシャルのパッケージを導入すると 　RS-232Cのコントロールが可能 例えば、COM3のモデムを操作する場合 fp=openserial(3,"9600,n,8,1") writeserial(fp,'ATI0'+ascii(13)) xpause(200000) readserial(fp) xpause(200000) writeserial(fp,'ATI1'+ascii(13)) xpause(200000) readserial(fp) closeserial(fp) こんな感じ。xpauseはパラメータのμ秒数だけ待つ命令 80 名前：名無しさん＠お腹いっぱい。 mailto:sage [2010/12/22(水) 13:33:15 .net] データを全部テキストデータで保存すると ファイル増えてしかたない sqliteあたりにまとめて保存できないでしょうか 81 名前：名無しさん＠お腹いっぱい。 mailto:sage [2010/12/23(木) 22:20:35 .net] >>80 saveコマンドがあります aaa=[1:1000]; save hoge.sav aaa とやれば変数aaaがhoge.savにテキスト形式で保存されます 読み出す場合は bbb=load('hoge.sav') とやればいいです この場合テキスト形式ですけど save -binary hoge.bin aaa とやればバイナリーで保存できますが、テキストよりもサイズが小さくなるかどうかはケースバイケースです 読み出しは ccc=load('hoge.bin') とやれば、勝手に保存形式を判断してくれます。 またMATLABのsaveコマンドで作られたファイルとの互換性もあるようです 82 名前：octave for windows メモの筆者 mailto:sage [2010/12/24(金) 17:39:00 .net] >>79 ・グラフがOctaveよりはきれい Ver. 3.3.54+ のFLTK Backendは,グラッフィクス，かなりよいですよ。 でも，自分でビルドしないと使えないのとMinGWではまだバグがある。 Ver 3.4はどれくらいかかるかな。 JWEさんのエンジンかかってきたから，数ヶ月後にはでるとおもうのですが。 83 名前：octave for windows メモの筆者 mailto:sage [2010/12/30(木) 05:32:50 .net] >>80 help save `-zip', `-z'とか`-hdf5'を使うと小さくなると思うのですが。 84 名前：octave for windows メモの筆者 mailto:sage [2010/12/31(金) 06:32:59 .net] 東京電機大学の松田先生の「Octaveの精義」がでました。 85 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/01/02(日) 08:27:31 .net] ttp://www.amazon.co.jp/dp/4877832319/ 86 名前：octave for windows メモの筆者 mailto:sage [2011/01/19(水) 18:39:37 .net] 開発版も様々なプラットフォームでのビルドがかなり安定してきています。 私は，MinGWのビルドしていますが状況はかなり改善しています。 3.4リリースが確実に近づいています。 新しい fltk Backend は,すばらしいものです。 できるだけ早く皆様が見る日がくればいいなと思います。 Debianなどでは，自力でビルドできる方がいらっしゃると思います。 ビルドできる環境におられるかたは 是非自力で開発版をビルドされて， その実力をみていただきたいと思います。 87 名前：octave for windows メモの筆者 mailto:sage [2011/01/24(月) 07:26:47 .net] Octave 3.4のリリースの準備としてOctave 3.3.90がリリース（ソースのみ） されました。 alpha.gnu.org/gnu/octave/ 多くの方にビルドいただき，3.4のリリースにご協力いただきたいと思います。 88 名前：octave for windows メモの筆者 mailto:sage [2011/01/28(金) 08:30:06 .net] Octave 3.3.90のMinGWでのビルドを行ってバイナリパッケージを作成してみました。 www.tatsuromatsuoka.com/octave/Eng/Win/ 89 名前：octave for windows メモの筆者 mailto:sage [2011/02/02(水) 09:18:42 .net] Octave 3.3.91のMinGWでのビルドを行ってバイナリパッケージを作成してみました。 www.tatsuromatsuoka.com/octave/Eng/Win/ いくつかのバグがフィックスされています。 90 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/02/08(火) 08:12:40 .net] >>89 テストパッケージ、ありがたく使わせてもらっています。 octaveと直接関係ないところで申し訳ありませんが。。。 パッケージに入っているckwなのですが、ckw.cfgの以下行が Ckw*chdir: d:\usr\Tatsu 自分の環境に存在しないディレクトリだとckwが 立ち上がらないようです。 91 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/02/08(火) 10:59:39 .net] >>89 入れてみましたけど octave:1> backend("fltk") error: `backend' undefined near line 1 column 1 となりますね。 plotするとgnuplotで描画されます。 92 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/02/08(火) 12:32:24 .net] >>91 graphics_toolkitに変わったみたいですよ（たしか） 93 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/02/08(火) 14:27:09 .net] >>92 確認しました octave:1> graphics_toolkit('fltk') octave:2> xx=[0:100]; octave:3> yy=sin(xx./5); octave:4> plot(xx,yy) これで表示されました 94 名前：octave for windows メモの筆者 mailto:sage [2011/02/09(水) 09:54:37 .net] >>90 ckw のご指摘ありがとうございました。 95 名前：octave for windows メモの筆者 mailto:sage [2011/02/09(水) 10:00:50 .net] Octave-3.4.0のソースがリリースされました。 www.gnu.org/software/octave/download.html 各ディストリビユーターによるバイナリのリリースはこれからだと思います。 なお，Windows版はMaintanerのBenjamin Lindnerさんの準備状況から 少し時間がかかると思います。 混乱を避けるため，3.4.0のtesting binaryを私の方から出すことは しないつもりです。 直前のスナップショットの3.3.92を出すかもしれません。 （たぶん中身はほとんど変わらない） 96 名前：octave for windows メモの筆者 mailto:sage [2011/02/10(木) 07:16:38 .net] GUI Octave　なるものがあるそうです。 https://sites.google.com/site/guioctave/ 97 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/02/10(木) 12:50:16 .net] >>96 GNUplotを使ってたり残念感がありますね ckwとgnuplotをwxt拡張で使ってる方がまだましというか どうせならScilabやR言語のグラフ機能を乗っ取ったとかの方が良かった 98 名前：octave for windows メモの筆者 mailto:sage [2011/02/27(日) 16:18:39.27 .net] >>97 まあ，3.4になってfltkになればよくなると思いますが。Windowsでは， まだしばらく3.4のリリースかかりそう。 99 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/03/27(日) 03:29:38.85 .net] 3.4をビルドしてみたけど、fltkでついにui系の機能が実装され始めたみたいだね。 今のとこuigetfile、uiputfile、uigetdir、uimenuだけだけど、今まで手つかずだっただけに素晴らしい！ ところで、fltkって画像ファイルの出力ってできないんだろうか？ graphics_toolkit fltk; plot(1:10); print('test.eps','-deps'); とかやっても、空っぽのファイルしか出来ないんだけど…。 他の形式でも同様。 …うん、printはgnuplot backend使えって話はごもっともなんだけど、アニメーション作成用に もうちょっと高速な出力手段を模索中。 誰か情報持ってない？ 100 名前：名無しさん＠お腹いっぱい。 [2011/04/15(金) 10:51:34.50 .net] すみません，2ch初心者なんですが…． Octaveのプログラミングについて質問させていただける板ってどこにありますか？ 「Octave」「octave」で検索かけてもこの板しか見つかりませんでした． エラーの返ってきちゃう関数をどうすれば使えるか調べたいんですけど(汗) 101 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/04/18(月) 22:25:16.10 .net] >>100 MATLABのスレを探してみてください OctaveはMATLABのクローンなので あと、エラーの返ってくる関数を調べる方法は知りませんが エラー処理はtry, catch でやります 例 A=[1,2] B=[1;2] try C=A+B catch %lasterr disp('error!!') D=A+B' end AとBはベクトルの向きが違うので、足すとエラーが出るのですが エラーの出そうな処理をtryの後に書き、エラーが出た後の処理をcatchの後に書きます コメントアウトしているlasterrは直前のエラーが何かを表示します 102 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/05/08(日) 12:20:29.28 .net] octave for windowsはビルドが難航してるのかな？ 103 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/06/13(月) 00:08:14.52 .net] octave から gnuplot を使うとき、最近のバージョンのoctave だと gset や __gnuplot_set とかが使えないんですが set term tgif みたいな処理はoctave からはどのように すれば良いんでしょうか？ 104 名前：octave for windows メモの筆者 [2011/06/16(木) 03:40:35.38 .net] >>103 最近バージョンでは，gnuplotのコマンドを直接たたくことはできません。 plotしたあと drawnow のdebug オプションを使って >> drawnow ("tgif", "foo.tgif", false, "debug.gp") で >> system('gnuplot debug.gp'); でfoo.tgifを作成する。 ぐらいでしょうか 105 名前：octave for windows メモの筆者 mailto:sage [2011/06/16(木) 03:56:40.97 .net] >>104 すみません。やってみましたが，うまくいきませんでした。<m(__)m> 106 名前：octave for windows メモの筆者 mailto:sage [2011/06/16(木) 04:07:12.45 .net] >>105 うまくいかなかったのはWindows 版のgnuplot のbug のせいのようです。 windowsでは system('gnuplot debug.gp'); system ('type debug.gp | gnuplot')でうまくいきそうです。 tgifをもっていないのでfoo.tgifをテキストエディタで確認しただけですが。 postscriptでうまくいったので，tgifでもおそらくいくと思います。 107 名前：octave for windows メモの筆者 mailto:sage [2011/06/16(木) 06:58:04.50 .net] >>103 たびたびすみません。ケアレスミスがあるので104-106は無視してください。 plotしたあと drawnow のdebug オプションを使って >> drawnow ("tgif", "foo.obj", false, "debug.gp"); でgnuplotのスクリプトdebug.gpを出力し それをsystemコマンドでgnuplotでgif出力（foo.obj）させます。 >> system('gnuplot debug.gp'); debug.gp　のデータはバイナリであり，windows 上ではgnuplotバクにより 正しくプロットされません。 108 名前：octave for windows メモの筆者 mailto:sage [2011/06/16(木) 07:02:12.97 .net] 107の続き windowsでは >> system ('type debug.gp | gnuplot'); としてください。 109 名前：octave for windows メモの筆者 mailto:sage [2011/06/16(木) 07:33:34.95 .net] 108　の続き tgif はもっていないので，fig ファイルを生成するようにして Cygwin 上Xfigで でfoo.figを作成したところうまくいきました。 　>> drawnow ("fig", "foo.fig", false, "debug.gp"); 　>> system ('type debug.gp | gnuplot'); 　できたfoo.fig をXfigでよみこんでプロットを確認しました。 tgifでもうまくいくと思います。 110 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/06/18(土) 10:54:38.14 .net] octave for windows メモの筆者さん >>103です。 いろいろお答えいただきましてありがとうございました。 tgif への出力はできるようになりました。 私がoctaveを使いこなしていないからかも 知れませんが 個人的には gnuplot の細かい設定が octave上からできる方が便利なのですが。 111 名前：octave for windows メモの筆者 mailto:sage [2011/06/18(土) 17:08:46.76 .net] >>110 Octave は3.2(3.2では試験的，3.4から本格的)から Open GL ベースのグラフィックスを扱うようになり， 整合性の観点から低レベルでgnuplotを直接触ることはやめになりました。 その代わり，octave-forgeにgnuplot パッケージがあり， gnuplotを直接いじりたい方はそちらを使うという思想のようです。 112 名前：octave for windows メモの筆者 mailto:sage [2011/06/18(土) 17:15:21.41 .net] >>111 gnuplotとパッケージのReferenceはこちらから， octave.sourceforge.net/gnuplot/overview.html 説明がわかりにくいのでつかったことはありません。 113 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/06/18(土) 20:09:09.69 .net] >>112 意外と良いかも 114 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/06/20(月) 11:49:46.70 .net] cygwin x上で以下のようなエラーが出てgnuplotを使うことができません．どなたか解決策が分かる方いましたら教えて下さい． octave:1> x = [0.0, 0.2, 0.4] x = 0.00000 0.20000 0.40000 octave:2> plot(x, sin(x)) error: Invalid call to strcat. Correct usage is: -- Function File: strcat (S1, S2, ...) error: called from: error: /usr/share/octave/3.4.0/m/help/print_usage.m at line 87, column 5 error: /usr/share/octave/3.4.0/m/strings/strcat.m at line 92, column 5 error: /usr/share/octave/3.4.0/m/plot/__gnuplot_get_var__.m at line 113, column 11 error: /usr/share/octave/3.4.0/m/plot/__go_draw_axes__.m at line 40, column 18 error: /usr/share/octave/3.4.0/m/plot/__go_draw_figure__.m at line 167, column 19 error: /usr/share/octave/3.4.0/m/plot/__gnuplot_drawnow__.m at line 86, column 5 115 名前：octave for windows メモの筆者 mailto:sage [2011/06/21(火) 00:54:07.79 .net] >>114 どうもCygwin1.7.9-1とOctave+gnuplotの相性が悪いようで Cygwin を1.7.8-1のダウングレードすると動きます。 116 名前：名無しさん＠お腹いっぱい。 [2011/06/21(火) 15:01:05.47 .net] 114です． cygwinのバージョンを1.7.8-1に直したところグラフ表示できました．ありがとうございました！！ 117 名前：名無しさん＠お腹いっぱい。 [2011/06/21(火) 21:24:25.60 .net] 最近octaveを始めたんだけど octaveのサポート状況の良いディストリってどれ？ ubuntu のoctave はバージョン古い。 FreeBSDのportsのoctaveなんかは最新のおっかけは凄いけど portsからのインストールは遅い。 118 名前：octave for windows メモの筆者 mailto:sage [2011/06/22(水) 03:04:37.49 .net] >>117 jweさんがDebianだからのぞいてみたけど，まだ3.2.4のままですね。 octaveの本家のMLをみるとLinuxの人は割りと自分でビルドしているようです。 119 名前：octave for windows メモの筆者 mailto:sage [2011/06/22(水) 06:18:06.49 .net] >>117 >>118 Unix はよくわかりませんがぐぐったら以下がでてきました。 octave Download - Linux Packages Search (RPM, DEB, XZ, TXZ, TGZ) pkgs.org/package/octave 3.4.0のようです。 120 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/06/24(金) 03:56:24.07 .net] FreeBSDでは、現在 3.4.1 です。 なんか凄い気合入って追っかけているような感じ。 121 名前：octave for windows メモの筆者 mailto:sage [2011/06/24(金) 04:04:41.54 .net] >>120 FreeBSDのOctaveのメンテナーの方を知っています。 （日本人です。） その方は多倍長BLASの作者でもあります。 最近連絡を取りました。 ところでoctave 3.4.1には，mkoctfileにバグが見つかり, 3.4.2の準備段階です。（現在RC2がでている）。 FreeBSDは，3.4.2がでたらすぐ対応するようです。 122 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/06/24(金) 23:22:10.88 .net] しばらくはFreeBSDだな。 Octave のためだけに。 123 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/06/25(土) 00:23:16.92 .net] VMwareにFreeBSD入れておくか

---

---

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

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