データ解析ツール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/ 294 名前：octave for windows メモの筆者 mailto:sage [2012/01/16(月) 06:47:46.80 .net] Octave-3.6.0がリリースされました。 octave.1599824.n4.nabble.com/Octave-3-6-0-uploaded-to-ftp-gnu-org-td4297674.html 295 名前：octave for windows メモの筆者 mailto:sage [2012/01/16(月) 12:17:54.38 .net] >>287 よくコードを読むと t=0:dt:t; がfor loopの中に入っているためn=2からは 区間ではなくなってますね。 t=0:dt:t をloopの外に出すと正しく動くようです。 以上 296 名前：octave for windows メモの筆者 mailto:sage [2012/01/26(木) 06:54:39.70 .net] Octave team で開発中のOctave-GUIの現在のスナップショット www.geocities.co.jp/tmoctwin/files/OctaveGUI0.8.5.html Windows上のスナップショットですが，Ubuntuでもほぼ同じイメージでした。 297 名前：octave for windows メモの筆者 mailto:sage [2012/01/26(木) 13:07:09.86 .net] >>213 bode線図でこけるのは，control packageのせいではではなく， fltk graphic_toolkit が対数軸でのマウスオペレーションがこける ためのようです。 https://savannah.gnu.org/bugs/?35292 298 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/01/27(金) 09:11:43.11 .net] >>296 Scilabも5.4.0-alphaからコマンドウインドウの左右に 補助ウインドウが出るようになりました ただしScilabは左が、変数＆履歴ウィンドウ、右がファイルブラウザ まあでもド素人にMatlab系のプログラミング技術を教えるだけなら FreeMATでいいかもしれないです インストールが簡単だし、dos画面じゃないし、 グラフもgnuplotとかではなくてネイティブだし 299 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/01/29(日) 11:13:02.52 .net] >>298 FreeMATって、コントロールパックとかあるんでしたっけ？ 伝達関数とか、それのインパルス応答とか、ボード線図かくやつ。 300 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/01/29(日) 14:56:20.91 .net] >>299 無いですね tfもzpもstepもbodeも 制御をやる目的ならscilabの方がいいと思います scilabはmatlabから入ると流儀がわからないかもしれませんが 慣れればどってことないです octaveは限界が高いですけど 初心者がやるにはハードルも高いと思います dos画面だけってのは初心者は面くらうでしょう それで電卓代わりに使う目的ならfreematで十分という話です 301 名前：octave for windows メモの筆者 mailto:sage [2012/01/30(月) 07:48:33.24 .net] >>300 OctaveはもともとUnix由来のソフトですし，Windowsの開発陣のリソースは実に 薄い。UnixではTerminalでのソフトはまだ一般的ですし，UnixのTerminalは コマンドプロンプトよりは使いやすい。 #Windows 2000からはコマンドプロンプトが正しい名称であって #ドスプロンプトではないです。まあ，それのほうがつうじますが。 302 名前：octave for windows メモの筆者 mailto:sage [2012/01/30(月) 07:54:27.43 .net] >>301 ちょっと舌足らずした。 とはいうもののUnixでもGUIのご時世なので，Octave本体もやっと 真剣にGUIインターフェースをつくりはじめています。 ちとお待ちください。 303 名前：octave for windows メモの筆者 mailto:sage [2012/01/31(火) 12:28:20.15 .net] Octave3.6.0_gcc4.6.2 mingw binaries for testing octave.1599824.n4.nabble.com/Octave3-6-0-gcc4-6-2-mingw-binaries-for-testing-td4338795.html インストールの仕方などは以下を参照してください。 www.octave.org/wiki/index.php?title=Octave_for_Windows#Octave-3.6.0-mingw_.2B_octaveforge_pkgs 304 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/01/31(火) 23:58:19.85 .net] 遅ればせながら、3.6.0のNEWSを読んでみた。 ぱっと目についた新機能は以下。 １．バイナリ演算子のブロードキャスト 　　ex.) rv=(1:5)'、M=magic(5)の時、repmat(rv,1,5).*Mがrv.*Mと書ける ２．プロファイラの実装 　　profile, profshow, profexplore ３．parfor（for文の並列処理）のサポート 　　但し、内部動作はforのまま parfor、一瞬期待したのに…！ 305 名前：octave for windows メモの筆者 mailto:sage [2012/02/01(水) 06:16:26.76 .net] >>303 Release 3.6.0 still not "official" octave.1599824.n4.nabble.com/Release-3-6-0-still-not-quot-official-quot-td4345497.html で3.6.1がでるようです。qhullのルーチンがごく最近改定されたのでちょうど いいかもです。Windows版のビルドもATLASのバージョンを変えてテストが 終わりつつあるのでこれもいいかもです。 parforですか。たしかにあるといいと思いますが今の開発陣のマンパワー だとどうかなという気がします。そこらへんに強い人が開発陣に参戦される ことを期待したいのですが。 306 名前：octave for windows メモの筆者 mailto:sage [2012/02/02(木) 10:06:13.73 .net] Octave for GNU Linux: Binary Octave packages for GNU Linux www.octave.org/wiki/index.php?title=Octave_for_GNU_Linux:_Binary_Octave_packages_for_GNU_Linux 307 名前：octave for windows メモの筆者 mailto:sage [2012/02/08(水) 07:05:13.92 .net] Octave 3.6.1-rc0 release candidate available for ftp octave.1599824.n4.nabble.com/Octave-3-6-1-rc0-release-candidate-available-for-ftp-td4365618.html 308 名前：名無しさん＠お腹いっぱい。 [2012/02/10(金) 15:43:20.25 .net] octaveってデフォルトで上三角の連立方程式 U x = y を高速に解けないの？ ちょっと使ってみようと思ってググったら2003年の記事出てきたんだけど、 www.math.meiji.ac.jp/~mk/labo/text/a-remark-lu.pdf これ未だにfixされてない？ 309 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/02/10(金) 23:24:48.60 .net] >>308 もしかして、試してすらいない？ うちの環境だとこんな感じ。 Octave 3.2.4, Linux 2.6.32, Core2Duo 1.8GHz octave:1> lutest(5000); P の掛け算の計算時間 Elapsed time is 0.0580444 seconds. L による割り算の計算時間 Elapsed time is 0.0310717 seconds. U による割り算の計算時間 Elapsed time is 0.0300882 seconds. leftdivu() による割り算の計算時間 Elapsed time is 0.843658 seconds. 精度 ans = 1.1294e-06 ans = 1.1197e-06 ※実際は100回計算してる。上の時間は平均値 310 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/02/11(土) 17:30:46.00 .net] >>309 試してというか、インストールすらまだしてないｗ ちょっとoctave使ってみようかな、ありがとう 311 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/02/14(火) 21:52:48.99 .net] 久々です 少し訳の分からないエラーが出たので質問したいのですが… 今サクラエディタ上で次のようなプログラムを作ってました function dx = f(x,t) dx = 1/(1+t^2); end t = linspace(0,10,101); x = lsode("f",1,t); 因みにプログラムの内容は微分方程式 X' = 1/(1+t^2)をlsodeを用いて解くというものです コレをコンパイラしたところ以下のようなエラーが出てしまいます warning: function name `f' does not agree with function file name (filename) parse error near line 5 of file (filename) syntax error >>> t = linspace(0,10,101); error: `t' undefined near line 2 column 12 error: called from: error: (filename) こんな感じのエラーでした しかし窓でプログラムを制作した場合は何故かエラーなしでグラフを表示することができました 何でこんなエラーが出たのでしょうか？ ^ 312 名前：octave for windows メモの筆者 mailto:sage [2012/02/15(水) 07:42:59.15 .net] >>311 今サクラエディタ上で次のようなプログラムを作ってました というのは，以下の内容を含むmファイルを作成して（例えばfoo.m） function dx = f(x,t) dx = 1/(1+t^2); end t = linspace(0,10,101); #**********: Octave のプロンプトからfooと打ったということでしょうか。 Octaveのmファイルには，スクリプトmファイルと関数mファイルがある。 Octaveはfunctionで始まるmファイルは関数mファイルと考える。 上の例は関数定義のあとにlsodeがあり，明らかに関数mファイルではない。 したがって上のスクリプトは間違いである。 warning: function name `f' does not agree with function file name (filename) は，関数mファイルは関数名とファイル名が一致しないといけないのでたwarning. 313 名前：octave for windows メモの筆者 mailto:sage [2012/02/15(水) 07:43:43.77 .net] >>311 >>312 上のスクリプトをスクリプトmファイルと認識させたいならfunctionキーワードの 前に命令を入れる。例えば 1; function dx = f(x,t) dx = 1/(1+t^2); end 1;は無意味な命令だがそのおかげで，octaveはこのファイルをスクリプトm ファイルと認識できる。 314 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/02/15(水) 20:58:11.93 .net] >>313 解決しました！ ありがとうございます！ 315 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/02/21(火) 22:33:37.23 .net] 流石に課題丸投げは行けないと思うが、悩みに悩んだ結果積んでしまったので教えてください ttp://loda.jp/subject_2ch/?id=762 これが課題の内容(pdf形式)です 只今小問２に関するプログラムを制作してるのですが、どうしても等高線グラフができません 一応原因はある程度把握してますが(ip(k,n)に関して、kが整数ではないからエラーがでる)それに変わるプログラムが分からず… 以下に小問２に関するプログラムを載せますので、どこがどう違うのかを指摘してくれれば有難いです 1; function dx = f(x,t,n1,k1) dx = -k1*x^n1; endfunction t = (0:1:300); x0 = 10; for k = 10^-9 : 0.1 : 1 for n = 1 : 5 x = lsode(@(x)f(x,t,n,k),x0,t); ip(k,n) = (x(21) - 8)^2 + (x(41) - 6)^2 + (x(61) - 5)^2 + (x(121) - 3)^2 + (x(181) - 2)^2 + (x(301) - 1)^2; endfor endfor contourf(pi) 316 名前：octave for windows メモの筆者 mailto:sage [2012/02/23(木) 11:42:18.48 .net] >>315 ttp://loda.jp/subject_2ch/?id=762 がダウンロードできないので，もとの問題についてはわかりませんが， ip(k,n) = (x(21) - 8)^2 + (x(41) - 6)^2 + (x(61) - 5)^2 + (x(121) - 3)^2 + (x(181) - 2)^2 + (x(301) - 1)^2; はkが正の整数かロジカルでないといけないので，こういう書き方は許されません。 317 名前：octave for windows メモの筆者 mailto:sage [2012/02/23(木) 11:52:06.35 .net] >>316 (ip(k,n)に関して、kが整数ではないからエラーがでる)って書いてありましたね。 すみません。 318 名前：octave for windows メモの筆者 mailto:sage [2012/02/23(木) 12:13:15.01 .net] >>315 1; function dx = f(x,t,n1,k1) dx = -k1*x^n1; endfunction t = (0:1:300); x0 = 10; k=1e-9; for m = 1:10 for n = 1 : 5 x = lsode(@(x)f(x,t,n,k),x0,t); ip(m,n) = (x(21) - 8)^2 + (x(41) - 6)^2 + (x(61) - 5)^2 + (x(121) - 3)^2 + (x(181) - 2)^2 + (x(301) - 1)^2; k=k+0.1; endfor endfor contourf(ip) は一応等高線はでるが，目指しているものと同じかわかりません。 319 名前：octave for windows メモの筆者 mailto:sage [2012/02/23(木) 13:27:38.33 .net] >>318 contourf(ip) は contourf(1:5, 1e-9:0.1:1, ip) か contourf(1e-9:0.1:1, 1:5,ip') にすればよいかも。 320 名前：octave for windows メモの筆者 mailto:sage [2012/02/23(木) 17:57:12.75 .net] >315-319 要は配列の概念とcontoturfの使い方がわかっていないのかなと思います。 ヘルプを読めばいいのですが，英語ですしね。無理もないかもしれません。 pdfをダウンロードできました。 もう少しいじる必要があると思いますが，基本的にはこれでいいのかな。 あとは自力で頑張ってください。 321 名前：octave for windows メモの筆者 mailto:sage [2012/02/23(木) 18:04:53.50 .net] >315-319 要は配列の概念とcontoturfの使い方がわかっていないのかなと思います。 ヘルプを読めばいいのですが，英語ですしね。無理もないかもしれません。 pdfをダウンロードできました。 もう少しいじる必要があると思いますが，基本的にはこれでいいのかな。 課題ですので，あとは自力で頑張ってください。 322 名前：octave for windows メモの筆者 mailto:sage [2012/02/28(火) 15:10:56.60 .net] Octave 3.6.1 がリリースされました。 octave.1599824.n4.nabble.com/GNU-Octave-3-6-1-Released-td4425272.html 3.6.0はバグのため公式版とならなかったため，3.6.1が公式版となります。 323 名前：octave for windows メモの筆者 mailto:sage [2012/02/28(火) 15:18:33.76 .net] まだ，公式には発表されていませんが，Michael氏がMSVC版のoctaveの復活に意欲 を見せておられます。さて，どうなるでしょうか。 324 名前：octave for windows メモの筆者 mailto:sage [2012/02/28(火) 15:18:51.21 .net] まだ，公式には発表されていませんが，Michael氏がMSVC版のoctaveの復活に意欲 を見せておられます。さて，どうなるでしょうか。 325 名前：octave for windows メモの筆者 mailto:sage [2012/03/03(土) 16:10:02.07 .net] MSVC版のoctave-3.6.1がでました。 octave.1599824.n4.nabble.com/Octave-3-6-1-for-Windows-VS2008-VS2010-available-td4440168.html 326 名前：名無しさん＠お腹いっぱい。 [2012/03/03(土) 23:33:37.71 .net] いきなりの質問ですみません。 ２つほど質問があります。 １．２段のワードの形式で論文を書いているのですが、 どうしても、ワードの大きさに縮小してしまうと、メモリどころか 線すらはっきり映りません。これをどう解決すればよいのでしょうか？ ２．グラフを描くときに、色ではなく破線や点線はどのようなコードを 打てばよいのでしょうか？ なにとぞよろしくお願いします。特に１番目の質問が深刻です。 327 名前：名無しさん＠お腹いっぱい。 [2012/03/03(土) 23:42:00.21 .net] ２番目は解決しました。すみません。 328 名前：octave for windows メモの筆者 mailto:sage [2012/03/04(日) 17:42:22.14 .net] >>326 ？？　プロットした図が縮小すると見えなくなるということでしょうか。 ワードだからWindowsですよね。octave-3.2.4でしょうか。 gnuplotの画面をコピーしてワードに貼り付けているということでしょうか。 私はoctave の図をそのまま使わないので良くわかりませんが， print コマンドでemfか何かでセーブした図を使ったほうがいいかもしれません。 329 名前：名無しさん＠お腹いっぱい。 [2012/03/05(月) 20:49:38.55 .net] octave-3.6.0で 1.1 ans = 1.10000000000000 single(1.1) ans = 1.10000002384186 となります。これってどうなんですか？ どなたか教えて下さい。 330 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/03/05(月) 22:52:08.58 .net] >>329 問題ないです 331 名前：octave for windows メモの筆者 mailto:sage [2012/03/07(水) 08:09:54.51 .net] Windows版の書き込みですみませんm(__)m。 MinGW版 の 3.6.1がでました。 sourceforge.net/projects/octave/files/Octave%20Windows%20binaries/Octave%203.6.1%20for%20Windows%20MinGW%20installer/ 332 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/03/07(水) 14:24:18.76 .net] >>329 補足するとsingle(単精度)の有効桁は10進に治すと約6-7桁なので それ以上の桁数の数値についての保証はなく　 >>330 の方がおっしゃるように，単精度の範囲では問題ないことになります。 333 名前：名無しさん＠お腹いっぱい。 [2012/03/29(木) 13:44:37.57 .net] octave 3.6.1-mingwで例えば text(0.5,0.5,"{\pi}","interpreter","tex") とやってもTeXでπで書いてくれません。 対処法を教えていただけませんか。 お願いします。 334 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/03/30(金) 08:04:27.54 .net] >>333 graphics_toolkit gnuplot でなんとかなりませんか？ octave 3.6.1-mingwはgraphics_toolkitはfltkがデフォで，texはまだfltkでは 対応していないと思うので。 335 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/06/20(水) 19:21:31.18 .net] 最近はJuliaというのが出てるらしいです julialang.org/ 336 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/06/21(木) 00:38:00.45 .net] >> 335 へー…。 そういえば、OctaveにもJIT Compilerを実装するプロジェクトが、 今年のGSoCに登録されてました。 これは、期待してもいいんでしょうか… jit-octave.blogspot.jp/ 337 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/06/21(木) 00:54:19.14 .net] SourceForgeに登録されてるWindowsバイナリの3.2.4、3.4.3、3.6.1 （MinGW版、いずれもクリーンインストール）を、現在XPで使用中 なのですが、3.2.4と比較して3.4.3、3.6.1の起動がやたらと遅くなって おります（30secくらい？CPU100%）。 どなたか、回避策をご存知の方とかいらっしゃいませんでしょうか？ 338 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/06/23(土) 01:27:08.49 .net] >>337 全く同じ事で悩んでいたけど、 windows7/64bit+Corei7+16GBRAMに買い換えたら起動が速くなったよ 339 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/06/23(土) 13:28:04.11 .net] >>338 おぉ、やはり他の環境でも再現することがあるのですね。 ぐぐっても何も出てこないので、困ってました。 時間ができたら、起動シーケンスを追っかけてみることにします。Thx! 340 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/06/25(月) 08:31:17.57 .net] >>337 >>338 MinGW版はXPではやたらと起動が遅くなるようです。 MSVC版だとそれほど遅くならないようです。 341 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/06/25(月) 22:54:11.86 .net] >>340 そうなんですか…情報感謝。 ただ、MSVC版は開発環境が揃ってないんで、自作octが コンパイルできないのがネックです。 コンパイラだけでなく、dllとかもぼろぼろ抜けてるような… 342 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/06/28(木) 15:54:52.26 .net] Scilabの5.4.0のbeta-1が出ました www.scilab.org/products/scilab/download/5.4.0-beta-1 343 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/07/03(火) 09:10:21.12 .net] >>341 >コンパイラだけでなく、dllとかもぼろぼろ抜けてるような… それはこまったことですね。 344 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/07/15(日) 19:18:59.66 .net] 337の件、追加でご報告。 どうも、時間が掛かってるのはoctaveのmain関数が呼ばれる前で、 main関数に入ってしまえばコンソールにアクセスできるまで 1秒もかかってません。 dllのロードが重いのかとrebase.exeとか使ってみましたが、効果なし。 私にはここから先のアプローチが思いつかないのですが、どなたか ヒントお願いします… 345 名前：名無しさん＠お腹いっぱい。 [2012/07/25(水) 12:26:52.32 .net] octave-3.6.2 mingw版ですが、peaksなどなんでも良いのですがグラフを書いてもちゃんと描いてくれません。refreshすると描いてくれます。以前にもこんなことありましたが対処法はありますか。 346 名前：名無しさん＠お腹いっぱい。 [2012/07/27(金) 09:14:20.16 .net] 345の件 graphics_toolkit("gnuplot") putenv("GNUTERM","wxt") でとりあえず希望どおりに動きました。 347 名前：名無しさん＠お腹いっぱい。 [2012/07/30(月) 20:02:25.49 .net] 345の続き putenv("GNUTERM","wxt")の状態でpeaksなどとしてグラフを描くと、グラフは ちゃんと描いてくれますがwindowsのターミナルが大きくなってしまいます。 どなたか教えてくれませんか。 最近、書き込みが少なくて寂しいです。 348 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/07/31(火) 19:53:06.42 .net] 答えてあげたいけど、Linux版を使ってるしグラフも余り書かないから答えてあげられない・・・ 349 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/08/17(金) 21:15:22.40 .net] ＞　男性は聯合ニュースの取材に対し、 暴行を受けている間も警察は積極的に止めようとせず、 男らは韓国人を差別する暴言を吐いたという。 警察はチョンの見方ですから有り得ません。 日本人にどうしても有事を惹起せしめたい欧米軍産複合体の意図じゃん 350 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/08/31(金) 22:01:03.02 .net] Congratulations, Max Brister! ttp://jit-octave.blogspot.jp/ Octave用のJIT Compilerに、とりあえず一通りの機能が揃ったらしい！ 速度も期待を裏切らず、例えば下のコードが、34sec->0.1secになるとか。 result = 0; for ii = 0:1e-7:1 　　result += inc * (1/3 * ii * ii); end Octave本家にマージもされたらしく、次回のリリースにexperimentalとして含まれるという話。 ソースはこちら。 ttp://inversethought.com/hg/octave-max 誰か試してみません？ 351 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/09/04(火) 20:33:58.95 .net] Windows(MSYS)版を使ってますが、ファイル選択ダイアログを表示してパスを 取得できる関数、あるいはそういう機能が使えるようになるパッケージって ないもんでしょうか？ 352 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/09/05(水) 13:37:47.23 .net] >>351 使ったことないのでよくはわからないが，uigetfile なるものがあるらしい。 ttp://www.gnu.org/software/octave/doc/interpreter/I_002fO-Dialogs.html#I_002fO-Dialogs 353 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/09/05(水) 20:50:05.79 .net] おぉ、望みどおりのことができました。 GTKなのか見た目はショボいですが、機能的には十分です。 354 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/09/08(土) 14:29:07.59 .net] 画像データを読み込んで行列にしたりそれを表示したりできますが、 この行列データに対して直線や円を描画する関数とかないですか？ 355 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/09/08(土) 18:34:23.91 .net] >>354 それは表示した状態に書き込むのか， それとも画像の実体(行列)に書き込むのか？ 356 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/09/08(土) 19:05:27.78 .net] 行列に対してです。 一応Bresenhamのアルゴリズムとか調べてみましたけど、既に実装 されているものがあればそっちを使いたいなぁと。 357 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/09/11(火) 01:08:11.79 .net] そういえばScilab-5.4.0がβ3まで来ましたが xcosが不安定で使い物にならないですね 5.3.3は安定なんですけど、AFFICH_mがバグって使えない（実行中に数値を表示しない） 358 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/09/19(水) 13:27:37.39 .net] Scilabが5.4.0のβ3が不安定なのを認めてひっこめましたね 原因が特定できたらしいので リリース時には急に落ちる事はなくなりそうです 359 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/10/20(土) 00:59:15.48 .net] MatlabのGUIを自宅でやりたいんだけどOctaveでできますか？ 360 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/10/21(日) 13:47:03.03 .net] できません 361 名前：名無しさん＠お腹いっぱい。 [2012/10/21(日) 18:48:24.67 .net] >>360 ありがとうございます では、他のフリーソフトで行うことはできませんか？ どうしても自宅でやりたいんですが・・・ 362 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/10/22(月) 22:36:45.68 .net] Javaででも書けば？ 363 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/10/22(月) 23:31:44.32 .net] >>361 SIMULINKをやりたいってことですかい？ 364 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/11/10(土) 01:21:39.57 .net] imcropがoctaveにはないけどトリミングはどうすればいいの？ 365 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/11/12(月) 10:39:31.83 .net] >>364 最新のoctave-forgeのimage パッケージにはimcropがあるようです。 octave.sourceforge.net/image/function/imcrop.html 366 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/11/15(木) 03:12:08.14 .net] >>365 ありがとうございます 367 名前：名無しさん＠お腹いっぱい。 [2012/12/08(土) 02:35:03.97 .net] leasqrを使って、２変数関数のフィッティングをやりたいんですが、 どう書けばいいですか？ octave.sourceforge.net/optim/function/leasqr.html によると、 x = vector or matrix of independent variables. と、独立変数がmatrixでもいいように書いてあるので、 ２変数でも使えるんじゃないかと思っています。 試しに、xを2次元で取って、それぞれ-1〜1の整数を代入してやると、 推定対象の関数f(x,p)の中で、x=(-1,0,1,-1,0,1)として認識されています。 その先、関数内で2つの独立変数を式にどう表せばいいのかが分かりません。 368 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/12/08(土) 03:00:32.70 .net] 前のスレから 256 名前：名無しさん＠お腹いっぱい。[sage] 投稿日：2006/06/27(火) 20:52:22 >>253 octave-forge から持ってきて、leasqrdemo.m をにらめっこしながら、 ようやく動かすことができました。 x=[1; 2; 3; 4; 5]; y=[1; 3; 6; 7; 8]; function ret=F(x,p) ret = p(1)*x.*x + p(2)*x + p(3); endfunction pin=[1 1 1]; [f,p]=leasqr(x,y,pin,"F") まず、xとyを列ベクトル標記に変え、さらにx*xの計算のところを x.*xに変えました。これで、動かし方を理解できました。 どうもありがとうございます。 369 名前：名無しさん＠お腹いっぱい。 [2012/12/08(土) 15:59:35.11 .net] [f,p]=leasqr(x,y,pin,"F")と、「x*x→x.*x」ですね。 ありがとうございました。 試してみます（月曜にならないとためせないんですが）。 370 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/12/08(土) 18:04:46.43 .net] Octave-3.6.2-vs2010-setup って、ライブラリもインストールできるのはいいんですが、 pkg rebuild -auto ってやってもロード出来ない（再起動しても、autoLoadされていない） のはしょうがないんですかね？ (pkg load all って打てば大丈夫だから、スタートアップに入れておけばいいんだけどさ） 371 名前：名無しさん＠お腹いっぱい。 [2012/12/09(日) 01:54:34.33 .net] >>368 >>369 ・・・って、思ったら「前のスレから256」は、独立変数xが１つの場合ですよね？ 知りたいのは、例えば独立変数xが２種類(x1とx2とします)あって、 (x1,x2)の組に対して観測値yが１つあり、 ret = p(1)*x1 + p(2)*x2 + p(3) みたいな場合に、x1とx2のそれぞれの取る値をmatrix xに入れた上で、 ret = F(x,p) の式で、matrix xの中のx1に相当する値とx2に相当する値を どう表現すればいいかってことなんですが・・・。 具体的には、画像座標(x,y)について、各整数画像座標位置に対して与えられる とある値（例えば標高値）を関数でフィッティングしたいんです。 372 名前：名無しさん＠お腹いっぱい。 [2012/12/14(金) 23:06:25.18 .net] 参考になりますかどうか。 あと、複素数で表すことも出来ます。 function z = f1(x,p) x0=p(1); y0=p(2); z = exp( -( (x(:,1)-x0).^2 + (x(:,2)-y0).^2 )/2 ); endfunction % ns = xpoints*ypoints; % z(ns,1), x(ns,1) = x of data, x(ns,2) = y of data % p(1, parameters) x0 = 0.1; y0 = 0.2; xt = [-4:0.5:4]; yt = [-4:0.5:4]; [xm, ym] = meshgrid(xt, yt); zm = exp(- ((xm-x0).^2+(ym-y0).^2)/2 ); plot3 (xm, ym, zm, "linestyle","none","marker","+"); input("type enter","s"); x = [xm(:),ym(:)]; z = zm(:); pin = [0.3, 0.4]; [f,p,cvg,iter] = leasqr(x,z,pin,@f1); p 373 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/12/15(土) 09:55:39.53 .net] ちょっと板違いかもしれませんが… WindowsでC++からOctaveを呼ぶ場合、cygwinかMSYSのgccでしか使えないんでしょうか？ VC++で通常のWindowsアプリから使える方法はあるでしょうか？ 374 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/12/16(日) 09:41:45.46 .net] ある。人並みの検索能力がある人間ならば。 375 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/12/19(水) 16:37:50.18 .net] >>373 wiki.octave.org/Octave_for_Windows#Using_the_Visual_C.2B.2B_compiler_with_Octave 376 名前：名無しさん＠お腹いっぱい。 mailto:sage [2013/01/27(日) 21:13:59.07 .net] www.lifeness.jp/525 377 名前：名無しさん＠お腹いっぱい。 mailto:sage [2013/01/27(日) 21:15:14.50 .net] www.lifeness.jp/714 378 名前：名無しさん＠お腹いっぱい。 [2013/01/28(月) 19:37:43.92 .net] octaveを使おうとインストールしたのですが、エラーが出てgnuplotが出ません。 error: feval: function `unimplemented' not found と出ます。どうすればいいのでしょうか。 379 名前：名無しさん＠お腹いっぱい。 mailto:sage [2013/01/29(火) 16:35:29.22 .net] >>378 Octaveのバージョンとプラットフォーム（WIndows,Mac, Ubuntu）などを明示してください 380 名前：名無しさん＠お腹いっぱい。 mailto:sage [2013/01/30(水) 00:04:49.12 .net] >>379 バージョンは3.6.1で、windowsです。 381 名前：名無しさん＠お腹いっぱい。 mailto:sage [2013/01/30(水) 14:14:09.52 .net] >>380 MSVC版ですか？　MinGW版ですか？ 名前にスペースのあるフォルダ　（C:\Program Files 等）にインストールして いませんか？ 382 名前：名無しさん＠お腹いっぱい。 mailto:sage [2013/01/30(水) 16:45:20.23 .net] >>381 申し訳ないのですが、自己解決しました。 ありがとうございました。 383 名前：名無しさん＠お腹いっぱい。 mailto:sage [2013/01/31(木) 11:42:30.56 .net] >>382 どうやって解決したの? 384 名前：名無しさん＠お腹いっぱい。 mailto:sage [2013/01/31(木) 11:44:01.01 .net] まず服を脱ぎました。 385 名前：名無しさん＠お腹いっぱい。 mailto:sage [2013/02/25(月) 23:38:58.12 .net] Windows VS2010版Octave 3.6.2のGUIで 実行するたびにコマンドラインが右にずれるのはなぜでしょうか 日本語ロケールがいけないのかな ちなみにスタート直後にmoreがかかってるので more off した方がいいです 386 名前：名無しさん＠お腹いっぱい。 mailto:sage [2013/03/04(月) 17:30:44.74 .net] >>385 試しに system("chcp 437") としてみてください。 387 名前：名無しさん＠お腹いっぱい。 mailto:sage [2013/03/05(火) 07:01:42.33 .net] >>386 カーソル位置はずれるけど、 shift-jis(932) よりはましですね 英語版Windowsならまともに動くのかもね 388 名前：名無しさん＠お腹いっぱい。 mailto:sage [2013/03/06(水) 13:58:23.32 .net] >>387 >カーソル位置はずれるけど、 Fileメニューのsetteingでterminalタブでfontを Lucida Consoleに選んて,chcp 437をかけると少しましになります。 389 名前：名無しさん＠お腹いっぱい。 mailto:sage [2013/03/06(水) 14:43:53.59 .net] Octave-3.7.2+/Windows [VS2010] available ttp://octave.1599824.n4.nabble.com/Octave-3-7-2-Windows-VS2010-available-td4650194.html 390 名前：名無しさん＠お腹いっぱい。 mailto:sage [2013/03/06(水) 20:47:35.66 .net] >>388 本当だ、すげ！ 使いものにならないと思ってたのが見違えた 391 名前：名無しさん＠お腹いっぱい。 [2013/03/07(木) 04:56:26.85 .net] ねえ　作った関数って保存できないの？ 392 名前：名無しさん＠お腹いっぱい。 mailto:sage [2013/03/07(木) 09:16:37.60 .net] >>389 開発版だとGUIがデフォのようですね。 コマンドライン版Octave-cliになってます。 結構安定に動くのでこちらを使おうかな。 393 名前：名無しさん＠お腹いっぱい。 mailto:sage [2013/03/07(木) 09:56:50.69 .net] >>391 エディタで関数を書いてfoo.mのような名前で保存すればいい。 394 名前：名無しさん＠お腹いっぱい。 mailto:sage [2013/03/07(木) 16:36:46.86 .net] OctaveのGUI環境にはWOctaveというのもあります sourceforge.net/projects/woctave/?source=directory 使い方はインストーラーでインストールした後に メニューバーのEdit>Preferenceで、OctavePathを入力するだけ OctavePathはOctaveをインストールしたbinディレクトリです

---

---

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

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