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

1 名前：名無しさん＠お腹いっぱい。 mailto:sage [2009/10/28(水) 22:13:01 ] 知る人ぞ知る強力なデータ解析ツール octave ... unixを科学研究で使っている香具師には非常に心強い味方のはずだ。 matlab互換を主張しながら、微妙にバージョンアップを続ける octaveについて語ろう。 そのPart 2です 前スレ データ解析ツールoctaveを語ろう pc12.2ch.net/test/read.cgi/unix/1086016412/ 200 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/12/02(金) 10:02:30.70 ] FreeMATが4.0から4.1に上がりました freemat.sourceforge.net/ Octave互換モードも出来たらしいです FreeMATとOctaveを統合してはどうかという話もありますがはたして 201 名前：octave for windows メモの筆者 mailto:sageg [2011/12/03(土) 10:38:58.08 ] >>200 情報ありがとうございます。確かにOctave Integrationとありますね。 OctaveのMLでは，なにも話しは出ていませんが，どうなんでしょうね。 Octave MLでは3.6のリリースに向けて議論がなされています。 GUIも開発されつつあります。 202 名前：octave for windows メモの筆者 mailto:sageg [2011/12/03(土) 20:10:07.47 ] >>200 Octave のMLで情報を流してみました octave.1599824.n4.nabble.com/FYI-FreeMat-4-1-includes-Octavee-Integration-td4151738.html レスポンスがどうなるかですね。 203 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/12/04(日) 10:21:08.63 ] FreeMatって G=1/(1+0.1*s) bode(G) みたいな事はできるんですかね？ Octaveは最新版でやっと出来るようになったみたいですが。 204 名前：：octave for windows メモの筆者 mailto:sage [2011/12/05(月) 07:09:35.32 ] control-2.2.2 released だそうです。 octave.1599824.n4.nabble.com/control-2-2-2-released-in-package-forum-please-upload-td4157808.html octave.svn.sourceforge.net/viewvc/octave/trunk/octave-forge/main/control/NEWS?revision=9260&view=markup 205 名前：：octave for windows メモの筆者 mailto:sage [2011/12/08(木) 07:03:41.15 ] control-2.2.3 released in package forum - please upload octave.1599824.n4.nabble.com/control-2-2-3-released-in-package-forum-please-upload-td4169495.html [2] octave.svn.sourceforge.net/viewvc/octave/trunk/octave-forge/main/control/NEWS?revision=9303&view=markup 206 名前：octave for windows メモの筆者 mailto:sage [2011/12/09(金) 12:18:21.65 ] NizanさんのWindows版Octave-3.4.3のバイナリがOctaveForgeにアップロード されました。 https://sourceforge.net/projects/octave/files/Octave_Windows%20-%20MinGW/Octave%203.4.3%20for%20Windows%20MinGW%20Installer/ 207 名前：名無しさん＠お腹いっぱい。 [2011/12/10(土) 02:50:59.58 ] >>206 おぉお。 ついに！　 ありがとうございます 208 名前：octave for windows メモの筆者 mailto:sageg [2011/12/10(土) 17:54:33.57 ] >>206 >>207 今回のバイナリーには，インストーラーがありません。そのため上記のページに 詳しい解説とインストール方法が記載されています。英語なので，拙Webページに自分の 改定を少し加えた翻訳の解説とインストール方法を記載しましたので，ご利用 ください。 www.tatsuromatsuoka.com/octave/jpn/OctaveWinMemo.html#111210 209 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/12/11(日) 23:57:09.22 ] cygwin(install all)いれるときに一緒に入れたoctave使いたい。 初期設定ファイルいじらずに、初めて起動した。 startxからoctave起動して、簡単な計算とかはできるんだけど グラフがまったく出力されない。↓こんなのでる。 gnuplotにパス通してないから？(´・ω・｀) octave:1> x=[0,0.2,0.4] x = 0.00000 0.20000 0.40000 octave:2> plot(x,sin(x)) 78 [main] octave-3.4.2 2720 exception::handle: Exception: STATUS_ACCESS_VIOLATION 1712 [main] octave-3.4.2 2720 open_stackdumpfile: Dumping stack trace to octave-3.4.2.exe.stackdump error: popen2: popen2 (child): unable to start process -- Bad address error: called from: error: /usr/share/octave/3.4.2/m/plot/__gnuplot_open_stream__.m at line 30, column 44 error: /usr/share/octave/3.4.2/m/plot/__gnuplot_drawnow__.m at line 72, column 19 error: Invalid call to strcat. Correct usage is: -- Function File: strcat (S1, S2, ...) error: called from: error: /usr/share/octave/3.4.2/m/help/print_usage.m at line 87, column 5 error: /usr/share/octave/3.4.2/m/strings/strcat.m at line 92, column 5 error: /usr/share/octave/3.4.2/m/plot/__gnuplot_get_var__.m at line 113, column 11 error: /usr/share/octave/3.4.2/m/plot/__go_draw_axes__.m at line 40, column 18 error: /usr/share/octave/3.4.2/m/plot/__go_draw_figure__.m at line 167, column 19 error: /usr/share/octave/3.4.2/m/plot/__gnuplot_drawnow__.m at line 86, column 5 210 名前：209 mailto:sage [2011/12/12(月) 00:01:14.96 ] 過去レスにサンプルまでいっしょのレスがありました。驚きです。 cygwinのダウングレードで解決されていましたが、 ダウングレードはするつもりがないので相性がよくなるまで待ちます。 211 名前：octave for windows メモの筆者 mailto:sageg [2011/12/12(月) 05:14:49.53 ] >>210 CygwinのMLを見ているとOctaveとの相性の悪い1.7.9に代わって1.7.10が テストリリースされているようです。まだ，試していませんが問題が解決して いることを願います。 212 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/12/13(火) 21:31:23.85 ] 　 213 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/12/13(火) 21:31:51.59 ] >>208 すいません s=tf("s") bode(s) ってやるとボード線図は出るのですが、出てきたボード線図のグラフィックウインドウ上でマウスをドラッグしたり ホイールをまわしたりすると、落ちちゃいます。 僕の環境（WinXP32bit)だけかもしれませんが。 普通にplotだと、出力されたグラフをマウスでいじっても大丈夫なんですが。 www.42ch.net/UploaderSmall/source/1323769033.png 214 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/12/15(木) 09:13:34.34 ] Mingw正式リリースの3.2.4を使っると巨大データでこけるので困っていました (Memory Exhaustedのエラーを吐く) 特にWindows7 64bitだとすごくこけやすかったので (理論的には2GBまで大丈夫なはずなのに500MBぐらいのメモリー使用でこける) 環境を変えるのはあまり好きでは無かったのですが 3.4.3を入れてみたら全くこけなくなりました、 このバージョンなら2GBまでこけないかも 助かりました OpenBLASは ttp://www.tatsuromatsuoka.com/octave/Eng/Win/ のを使わせていただきました 自分でmakeするとエラーが出たので 215 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/12/15(木) 19:28:05.14 ] >>214 MemoryExhaustedは3.2.4でも悩まされてるけど、 3.4.3は、さらに起きやすくなった。実際、同じスクリプトが動かなくなった。 216 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/12/16(金) 09:38:28.61 ] >>215 こちらはメモリーを4GB積んでます それとWin7 32bit では試してないので WinXP 32bit で 3.2.4 Win7 32bit で 3.2.4 Win7 64bit で 3.2.4 WinXP 32bit で 3.4.3 Win7 32bit で 3.4.3 Win7 64bit で 3.4.3 この６つの組合せとメモリー搭載量で起きやすい、起きにくいがあるのかもしれません よっぽどMigw64bit でmakeしようと思ったけど、色々挫折した(笑 3.4.3を試そうと思ったのはOctaveのMailingリストで、3.2.4から3.4に上げれば Memory Exhaustedが出ないよという書き込みがあったから 217 名前：名無しさん＠お腹いっぱい。 [2011/12/18(日) 02:33:58.07 ] iphone用octaveでないかなー 218 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/12/18(日) 11:03:15.98 ] >>217 そこでLinuX zaurusですよ 219 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/12/19(月) 10:36:54.47 ] Scilab 5.4.1のalpha版が出ました www.scilab.org/products/scilab/download/5.4.0-alpha-1 5.3.3は使っていてbugだらけだったので、少しはましになるかどうか (xcos の AFFICH_mでリアルタイムで数値が表示されないのは 　仕様じゃなくてbugだったのか orz) 220 名前：octave for windows メモの筆者 mailto:sage [2011/12/20(火) 03:33:44.62 ] 最近，HDDが壊れたThink Pad X40にUSB HDDを利用してUbuntu 11.10をいれて触っています。 Ubuntuのソフトウェアセンターのoctaveは3.2.4なので，3.4.3をビルドしてみました。 依存ライブラリが揃っており，Synapticパッケージマネージャーで簡単にインストールできました。 Windowsでのビルドの経験が役にたって少々の試行錯誤で無事ビルドできました。 なかなか，いい感じです。しかし，ビルドツールがUnbuntu上ではWindows上より はるかに速いのでおどろきました。 また，3.6.0のプレリリースである3.5.90もビルドできました。 221 名前：octave for windows メモの筆者 mailto:sage [2011/12/20(火) 03:37:19.84 ] 最近，HDDが壊れたThink Pad X40にUSB HDDを利用してUbuntu 11.10をいれて触っています。 Ubuntuのソフトウェアセンターのoctaveは3.2.4なので，3.4.3をビルドしてみました。 Windowsでのビルドの経験が役にたって少々の試行錯誤で無事ビルドできました。 なかなか，いい感じです。 しかし，ビルドツールがUnbuntu上ではWindows上よりはるかに速いのでおどろきました。 また，3.5.90のビルドもできました 222 名前：octave for windows メモの筆者 mailto:sage [2011/12/21(水) 07:20:59.27 ] 最近，HDDが壊れたThink Pad X40にUSB HDDを利用してUbuntu 11.10をいれて 触っています。Ubuntuのソフトウェアセンターのoctaveは3.2.4なので， 3.4.3をビルドしてみました。 Windowsでのビルドの経験が役にたって少々の試行錯誤で無事ビルドできました。 なかなか，いい感じです。しかし，ビルドツールがUnbuntu上ではWindows上より はるかに速いのでおどろきました。 223 名前：octave for windows メモの筆者 mailto:sage [2011/12/21(水) 07:22:37.20 ] 最近，HDDが壊れたThink Pad X40にUSB HDDを利用してUbuntu 11.10をいれて 触っています。Ubuntuのソフトウェアセンターのoctaveは3.2.4なので， 3.4.3をビルドしてみました。 Windowsでのビルドの経験が役にたって少々の試行錯誤で無事ビルドできました。 なかなか，いい感じです。しかし，ビルドツールがUnbuntu上ではWindows上より はるかに速いのでおどろきました。 224 名前：octave for windows メモの筆者 mailto:sage [2011/12/23(金) 10:37:59.12 ] ＞＞220　-223 ２ｃｈブラウザの特性を知らず同じ書き込みを何度もしてすみません。m(__)m 225 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/12/23(金) 10:49:14.91 ] >>214 >OpenBLASは >ttp://www.tatsuromatsuoka.com/octave/Eng/Win/ >のを使わせていただきました >自分でmakeするとエラーが出たので MinGW上ではふつうにビルドしただけではエラーになるようです。 同時にできる静的ライブラリから，次の例のようにしてdllを作成しました。 i686-pc-mingw32-gcc -shared-libgcc -shared \ -Wl,--enable-auto-import \ -Wl,--enable-runtime-pseudo-reloc-v2 \ -Wl,--allow-multiple-definition \ -Wl,--out-implib=libblas.dll.a -o libblas.dll \ -Wl,--output-def=libblas.exp \ -Wl,--whole-archive libopenblas_nehalemp-r0.1alpha2.4.lib \ -Wl,--no-whole-archive -lgfortran こうしてできたものを ttp://www.tatsuromatsuoka.com/octave/Eng/Win/ にアップしました。 226 名前：octave for windows メモの筆者 mailto:sage [2011/12/23(金) 10:51:16.90 ] >>225 i686-pc-mingw32-gcc -shared-libgcc -shared \ で始まる行が一行ずつ開いているのは詰めてください。m(__)m 227 名前：octave for windows メモの筆者 mailto:sage [2011/12/23(金) 12:16:03.65 ] >>255 たびたびすみません。 >MinGW上ではふつうにビルドしただけではエラーになるようです。 Cygwin上でふつうにビルドしただけではMinGWでエラーになるようです。 で，CygwinのMinGW64の32bit版GCCで静的ライブラリからリンクしなおして， dllを作成しています。 228 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/12/23(金) 13:13:13.43 ] >>227 どうもありがとうございました dllのビルドがおかしかったんだ 229 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/12/23(金) 13:22:53.36 ] >>224 あ、おかえりなさいｗ 普通は自動追尾するはずなんですが、 板欄→板一覧の更新 スレ欄→スレ一覧の更新 を手動でやればいいんでしたっけか。 自分も適当にやったので、本当はどうすればいいのか、実はよくわかってないですｗ というわけなんですが、 >>213 の状態になるのは僕の環境だけなんでしょうかねぇ？ 230 名前：octave for windows メモの筆者 mailto:sage [2011/12/23(金) 16:04:18.95 ] >>213 >>229 再現確認しました。 現在配布のバイナリは，share\octave\site\m\startup\octaverc でfltkがgraphics_toolkitに指定してあります。 graphics_toolkit gnuplot でやると落ちないようです。 gnuplotのmultiplotを使っているので，マウスぐりぐりは 変なことになりますが，落ちることはないようです。 fltkで落ちるのは，バグかなと思います。 ちょっと調べてみます。 231 名前：octave for windows メモの筆者 mailto:sage [2011/12/23(金) 16:04:31.43 ] >>213 >>229 再現確認しました。 現在配布のバイナリは，share\octave\site\m\startup\octaverc でfltkがgraphics_toolkitに指定してあります。 graphics_toolkit gnuplot でやると落ちないようです。 gnuplotのmultiplotを使っているので，マウスぐりぐりは 変なことになりますが，落ちることはないようです。 fltkで落ちるのは，バグかなと思います。 ちょっと調べてみます。 232 名前：octave for windows メモの筆者 mailto:sage [2011/12/23(金) 16:06:37.91 ] また，２重投稿になってるm(__)m 233 名前：octave for windows メモの筆者 mailto:sage [2011/12/23(金) 16:16:09.81 ] >>230 subplotが怪しいのかなあとおもいましたが，単にsubplotしただけではfltkでの グラフ落ちませんでした。Controlパッケージのせいなのかな？ ちょっと見当つかないですが，余力があればBug Trackerにあげてみます。 それで解決するかはわかりませんが。 234 名前：octave for windows メモの筆者 mailto:sage [2011/12/23(金) 16:32:41.70 ] >>233 bug trackerにあげときました。 savannah.gnu.org/bugs/?35159 235 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/12/23(金) 18:01:28.68 ] >>234 ありがとうございます。 お世話かけます。 236 名前：octave for windows メモの筆者 mailto:sage [2011/12/24(土) 12:24:11.51 ] Octave-3.5.91がでました。 ttp://octave.1599824.n4.nabble.com/Octave-3-5-90-snapshot-available-for-ftp-tp4200255p4230198.html 3.6.0へ向けての最後のスナップショットとなると思います。 237 名前：octave for windows メモの筆者 mailto:sage [2011/12/24(土) 18:57:01.99 ] 少し前ですが，Octave のWikiが新しくなりました。　 ttp://www.octave.org/wiki/index.php?title=Main_Page の左の Navigation の index からいろいろなページに飛べます。 238 名前：octave for windows メモの筆者 mailto:sage [2011/12/25(日) 05:39:00.48 ] ttp://sourceforge.net/projects/octave/files/Octave_Windows%20-%20MinGW/Octave%203.4.3%20for%20Windows%20MinGW%20Installer / にインストーラー版 Octave-3.4.3 for windowsがおかれました。 ttp://gnu-octave-repository.2306053.n4.nabble.com/Installer-for-mingw32-Octave-3-4-3-td4225772.html Octave-forgeは含まれないようなので，従来からあるOctave-Forgeを利用される 方は Octave3.4.3_gcc4.5.2_pkgs_20111026.7z の内容をを上書きして， 設定する必要があると思います。 239 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/12/25(日) 22:30:58.01 ] 本当に初歩的な質問ですが、octaveで乱数(メルセンヌツイスタ)を使いたいのですが、その乱数がcファイルなので、どうすればoctaveで使えるのですか？ 因みに目的はメルセンヌツイスタを用いた乱数発生で円を求める奴です 240 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/12/26(月) 10:42:26.06 ] >>239 help rand を見ると、現在のOctaveのrandで用いているのは 擬似Mersenne Twisterのようです 241 名前：octave for windows メモの筆者 mailto:sage [2011/12/26(月) 13:59:09.11 ] >>239 初歩的な質問でなく，かなり高度な質問（バリア）も高い質問です。 www.gnu.org/software/octave/doc/interpreter/Oct_002dFiles.html#Oct_002dFiles を読んでoct ファイルを作る。C++だから，ラッパーを書く必要があるかも。 MATLABでも使うなら,mexがいいかもしれない。 www.gnu.org/software/octave/doc/interpreter/Mex_002dFiles.html#Mex_002dFiles 英語だから苦労するかもしれないけど，勉強だと思って。 242 名前：octave for windows メモの筆者 mailto:sage [2011/12/26(月) 14:00:58.82 ] >>241外した答えになっているような気がする。許してください。 243 名前：239 mailto:sage [2011/12/26(月) 16:06:13.41 ] >>240 ごめんなさい、メルセンヌツイスタ縛りなのでrand関数は不可です あくまでも公式サイトから落としたメルセンヌツイスタプログラムを用いて計算しなければならないので… てか何故PCから見れないんだ、このスレ(´・ω・`) 244 名前：octave for windows メモの筆者 mailto:sage [2011/12/26(月) 17:11:53.04 ] いまは↓ ttp://toro.2ch.net/test/read.cgi/unix/1256735581/ で見れるようです。私もある方に教えてもらいました。 245 名前：octave for windows メモの筆者 mailto:sage [2011/12/26(月) 17:19:19.18 ] >>243 ttp://www.math.sci.hiroshima-u.ac.jp/~m-mat/MT/MT2002/CODES/mt19937ar.c ソースはこれですか。 コンパイルして，実行形式を作る。 そのあと，Octaveから system コマンドをつかって実行して， load コマンドなどでOctaveの変数に取り込むのが手っ取り早いかと。 246 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/12/26(月) 18:33:11.68 ] それです、正にそれです 流石に馬鹿過ぎる質問ですがコンパイルと実行ファイル制作ってどうするんですか？ コマンドプロンプトを使ってやるのですか？ それともoctave上でコンパイルするのですか？ あ、当然ググりました 247 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/12/26(月) 19:02:34.54 ] >>246 課題？ $ gcc -o mt19937ar mt19937ar.c $ ./mt19937ar 1000 outputs of genrand_int32() 1067595299 955945823 477289528 4107218783 4228976476 3344332714 3355579695 227628506 810200273 2591290167 (以下略) 248 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/12/26(月) 23:46:31.06 ] >>247 無事に窓を用いてコンパイルは出来ました(mt19937ar.exe) けどそこからがよく分かりません… system関数？load関数？ 因みにこれは微妙に課題です ゼミでoctaveを使うことになったのでその予行練習課題みたいなものです 249 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/12/27(火) 00:10:48.73 ] >>239 確認ですが、その問題というのは、 0〜1の乱数を２個発生させる。(r1,r2) r1^2+r2^2<　1　だったら　カウンタを１増加させる。 この作業を100万回行う。 カウンタの値/100万×4がπになるはず っていう問題ですかい？ 問題を丸投げするんだから、情報を小出しにせず、 堂々と丁寧に全部書くことをオススメします。 250 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/12/27(火) 00:40:58.98 ] 了解しました 課題は2次元球(円)から10次元球までの体積(面積)をoctaveを用いて解くものです 条件としましては ・乱数を用いること、その際条件としてメルセンヌツイスタを用いる ・半径は1(x1+x2+x3+…+xn≦1を満たす、但しx1,x2,x3,…,xnはメルセンヌツイスタを用いて発生させた乱数とする) ・乱数発生回数は問わない(但し回数は考察対象とする) 因みに背景としましては、自分卒研にてベイズ推計を用いての逆問題について研究していますので、その際用いられるMCMC法でメルセンヌツイスタ乱数およびoctaveが必要なので、先生から予行練習課題としてこの課題が出されたのです 251 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/12/27(火) 00:46:21.14 ] 訂正です 乱数の範囲は[0,1]です そして体積はV=(octaveで求めた値)×2^n(n:次元)で求めるものとします 一応octaveの基本はある程度理解をしたのですが、高度なのはさっぱり分からないのでここで聞きました 252 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/12/27(火) 00:49:10.19 ] 書き忘れた…orz >>249で書いてるのと一緒です オレの馬鹿、連レスサーセンです 253 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/12/27(火) 01:24:12.38 ] >>252 多分、方法は２つあって、 １）　メルセンヌツイスタのアルゴリズムをOctaveで書き下ろして使う ２）　dirk.eddelbuettel.com/code/octave-mt.html　を読む のどっちかじゃないすか？ ＃そもそもなんで、この問題をわざわざoctaveでやるのかの意味がよくわからんです。 　octaveでどんな問題でも解けるように　というトレーニングなら１）の方法でしょーね。 254 名前：octave for windows メモの筆者 mailto:sage [2011/12/27(火) 06:47:11.09 ] >>253 番さんの ２）の方法を使うためには，>>241のoctファイルを作成する方法をOctaveの マニュアルでよんでください。 randmt.cc　をダウンロードして mkoctfile コマンドをつかって，Octファイル作ればいいのですが， なんかエラーがでてしまいます。 octave:2> mkoctfile randmt.cc randmt.cc:42:17: warning: inline function 'uint32 randomMT()' used but never defined randmt.o: In function `randomMTunit': E:\usr\Tatsu\test\octave/randmt.cc:48: undefined reference to `randomMT' randmt.o: In function `Z11FrandmtseedRK17octave_value_listi': E:\usr\Tatsu\test\octave/randmt.cc:105: undefined reference to `seedMT' randmt.o: In function `randomMTunit': E:\usr\Tatsu\test\octave/randmt.cc:48: undefined reference to `randomMT' E:\usr\Tatsu\test\octave/randmt.cc:48: undefined reference to `randomMT' E:\usr\Tatsu\test\octave/randmt.cc:48: undefined reference to `randomMT' E:\usr\Tatsu\test\octave/randmt.cc:48: undefined reference to `randomMT' E:\usr\Tatsu\test\octave/randmt.cc:48: undefined reference to `randomMT' randmt.o:E:\usr\Tatsu\test\octave/randmt.cc:48: more undefined references to `randomMT' follow collect2: ld returned 1 exit status うーむ。 255 名前：octave for windows メモの筆者 mailto:sage [2011/12/27(火) 06:57:19.79 ] >>254　ページをよく読んでいなかったorz。 Download code.からoctave-mt.tar.gzをダウンロードする。 適当なフォルダに解凍する。 たとえば，E:\usr\Tatsu\test\octave\octave-mt に解凍したら， octave:3> cd E:\usr\Tatsu\test\octave\octave-mt octave:4> system('make') のようにすると gcc -Wall -O6 -c -o cokus.o cokus.c mkoctfile randmt.cc cokus.o randmt.cc:42:17: warning: inline function 'uint32 randomMT()' used but never defined ln -sf randmt.oct randmtseed.oct ln -sf randmt.oct randmtn.oct ans = 0 とoctファイルができる。 helpを読むには， octave:5> help randmt とする。 256 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/12/27(火) 12:06:42.27 ] やってみました octave:1> cd "C:/Documents and Settings/DateMasa/デスクトップ/MT" octave:2> system('make') mkoctfile randmt.cc cokus.o randmt.cc:42: warning: inline function `uint32 randomMT()' used but never define d c:\MinGW\bin\..\lib\gcc\mingw32\3.4.5\..\..\..\..\mingw32\bin\ld.exe: unrecogniz ed option '--enable-runtime-pseudo-reloc-v2' c:\MinGW\bin\..\lib\gcc\mingw32\3.4.5\..\..\..\..\mingw32\bin\ld.exe: use the -- help option for usage information collect2: ld returned 1 exit status make: *** [randmt.oct] Error 1 ans = 2 こんな感じのエラーが出てしまいました… 確認で一応 help randmt　を実行しましたが「そんな関数なんてねーよ」と吐き出してきました 後マニュアルに書いてあるとおりmkoctfileでccファイルをoctファイルにしてみようと試みたのですが これもエラーが出てしまいました 257 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/12/27(火) 12:23:14.32 ] 因みにmkoctfileでのエラーは次のとおりです octave:6> mkoctfile randmt.cc g++: randmt.cc: No such file or directory g++: no input files g++: randmt.o: No such file or directory 258 名前：octave for windows メモの筆者 mailto:sage [2011/12/27(火) 12:50:27.95 ] >>256 >>257 パスにスペースがあるとだめですよ。 だから cd "C:/Documents and Settings/DateMasa/デスクトップ/MT" ではだめです。 259 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/12/27(火) 13:54:13.73 ] アドバイスを参考にスペースが入らないよう解凍位置を変更してみました octave:3> cd C:/MT-file octave:4> system('make') mkoctfile randmt.cc cokus.o randmt.cc:42: warning: inline function `uint32 randomMT()' used but never define d c:\MinGW\bin\..\lib\gcc\mingw32\3.4.5\..\..\..\..\mingw32\bin\ld.exe: unrecogniz ed option '--enable-runtime-pseudo-reloc-v2' c:\MinGW\bin\..\lib\gcc\mingw32\3.4.5\..\..\..\..\mingw32\bin\ld.exe: use the -- help option for usage information collect2: ld returned 1 exit status make: *** [randmt.oct] Error 1 ans = 2 これでhelp randmtをやっても「だからねーって言ってるだろ」とエラーが出てしまいました どこが違うのか… 260 名前：octave for windows メモの筆者 mailto:sage [2011/12/27(火) 14:08:52.94 ] >>259 c:\MinGW\bin\..\lib\gcc\mingw32\3.4.5\..\..\..\..\mingw32\bin\ld.exe: unrecogniz ed option '--enable-runtime-pseudo-reloc-v2' c:\MinGW\bin\..\lib\gcc\mingw32\3.4.5\..\..\..\..\mingw32\bin\ld.exe: use the -- help option for usage information collect2: ld returned 1 exit status make: *** [randmt.oct] Error 1 ans = 2 ちがったGCCのバージョンが認識されてますね。（3.4.5） Octave 3.4.3ではGCC-4.5.2がパッケージされていて， それが認識されないといけないのですが。 ttp://octave.1599824.n4.nabble.com/Re-Fw-Re-msys-from-octave-td4234096.html が参考になるかもです。 ちなみに'--enable-runtime-pseudo-reloc-v2’はGCC-3.4.5にはないオプションです。 261 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/12/27(火) 15:44:46.34 ] minGWを新しいものに変え、今までの事をやりました 無事randmtが起動できました そして課題もまとめることができました 本当にありがとうございます！ 262 名前：octave for windows メモの筆者 mailto:sage [2011/12/28(水) 06:35:47.21 ] >>261 よかったですね。卒研頑張ってください。 263 名前：octave for windows メモの筆者 mailto:sage [2011/12/28(水) 06:58:12.76 ] >>213 >>234 savannah.gnu.org/bugs/?35159 Control PackageとFLTKの組み合わせでBode線図にグラフに問題があるのは， Windows だけでないようです。 3.6.0で改善すればよいのですが。 3.6のプレリリースの3.5.90や現在developmentソースでControl packgeが 組み込めない状態なのでちょっと時間がかかると思います。 しばらくは，gnuplot を graphics_toolkit にしてお使いください。 264 名前：名無しさん＠お腹いっぱい。 [2011/12/28(水) 22:51:43.25 ] 本当に初歩的な質問で申し訳ございません。 今octave3.4.3をインストールして使おうと思っています。 octaveはmatlabと非常に相関性があるということで家ではこのoctaveを使おうと思っています。 そこでいくつか質問があります。 １．octaveを起動させ、matlabコードを張り付けて実行したいのですが 張り付けすることができません。 いちいち、タイプしなければならないんでしょうか？ ２．matlabでは、memo帳（functionやscript）が出てきてそこにコードを打って 実際のmatlab上で実行させるのですが、このoctaveはメモ帳のようなものはなく いちいち打たなければならないのでしょうか？ 調べてみると、matlabと同じようにfunctionがあるみたいですが、 これはfunctionが別のところになければ実行できないのですが、 どうすればよいのでしょうか？ よろしくお願いいたします。 265 名前：octave for windows メモの筆者 mailto:sage [2011/12/29(木) 06:23:58.44 ] >>264 文章からして、windows上のことと思います。Octaveは、マルチプラットフォーム なソフトなので、OS(Windows, Mac, Linix)などを記載してください。 私の書く文章にわからない点があると思いますが適宜Googleなどで 調べてください。 １．Windows上では、Octaveはコマンドプロンプト上で動作しています。 octaveの実行ウィンドウをマウスで右クリックして、プロパティを開き、 オプションタブで簡易編集モードがアクティブになっていることを確認 してください。octave　のプロンプトに張り付けたい文章をクリップボード にコピーしてoctaveのプロンプト上で右クリックをすれば、張り付ける ことができます。 これは、Octaveの機能というよりコマンドプロンプトの機能ですので 「コマンドプロンプト　張り付け　コピー　簡易編集モード」などで ググってください。NotePad++を既定のテキストエディタにするには、 Octaveをダウンロードされたページの説明（英語）をよく読んでください。 設定をすれば、octave のプロンプトからeditとうてば起動できます。 ２．将来的には、octaveもMATLABのようなGUIを備えたソフトになるべく 作業がされていますが、現在はGUIソフトではありません。 テキストエディタは自分で準備してください。 Notepad++がいいのではと思います。 ＞調べてみると、matlabと同じようにfunctionがあるみたいですが、 これはfunctionが別のところになければ実行できないのですが、 どうすればよいのでしょうか？ なにが言いたいのかよくわかりません。もう少し実例をあげてわかりやすく 質問してください。 266 名前：octave for windows メモの筆者 mailto:sage [2011/12/29(木) 06:28:52.00 ] >>265 いかん文章の推敲があまいなあ。 NotePad++を既定のテキストエディタにするには、 Octaveをダウンロードされたページの説明（英語）をよく読んでください。 設定をすれば、octave のプロンプトからeditとうてば起動できます。 は ２．の部分の説明でNotepad++ を既定のエディタにするために何を読めばいいか をかきました。 なお、拙ページに説明の翻訳をあげていますので、よければ参考にしてください。 ttp://www.tatsuromatsuoka.com/octave/jpn/OctaveWinMemo.html#111210 267 名前：octave for windows メモの筆者 mailto:sage [2011/12/29(木) 06:34:02.08 ] >>264 バージョンは、3.2.4ですが、私のoctaveを使った講義の資料を紹介した部分を あげていますので参考になれば幸いです。 ttp://www.tatsuromatsuoka.com/octave/jpn/OctaveWinMemo.html#110613 268 名前：名無しさん＠お腹いっぱい。 [2011/12/29(木) 11:36:23.43 ] どういうわけかグラフを一つ書きません。 訳がわかりません。 ## WindowsXP + octave 3.2.4 and 3.4.3 clf for i = 1:4 subplot(4,4,i+12); h=plot(randn(1,5),"r+"); set(gca(),"position",[(i-1)*0.225+0.1, 0.25, 0.2,0.2]) endfor input("3番目が抜けているはず。書き直すと今度は書く。改行で継続","s"); i = 3; subplot(4,4,i+12); h=plot(randn(1,5),"r+"); set(gca(),"position",[(i-1)*0.225+0.1, 0.25, 0.2,0.2]) 269 名前：名無しさん＠お腹いっぱい。 [2011/12/29(木) 13:53:53.94 ] Notepad++（テキストエディタ） の インストール（オプション，推奨） 最新の Notepad++ を notepad-plus-plus.org/ からダウンロードし，インストールします。 <Install folder>\share\octave\site\m\startup\octaverc をエディタで編集する。 次の行のコメントを解除し，Octave の既定のエディタに設定します。 EDITOR('C:\\Program Files\\Notepad++\\notepad++.exe'); 注意：notepad++.exe がインストールされるフォルダが上記でない場合は，設定を適宜変更します これをやったのですが、 editとうっても error: C:\Octave\Octave3.4.3_gcc4.5.2\Octave3.4.3_gcc4.5.2\share\octave\3.4.3\ m\miscellaneous\edit.m at line 212, column 7 と出てきて起動できません。 270 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/12/29(木) 20:44:28.38 ] >>264 エスパー的解釈で補完して答えてあげやう。 どうせwindows版だろうし。 １）C:\Octave　というフォルダがあるだろうから（なければ作る） その中に MyProj とかいうフォルダを作る。 １）Windowsのメモ帖を立ち上げる。 ２）そこに a=1 b=a+5 ってな感じでmatlabっぽいコードを書いて見る。 ３）　名前をつけて保存で　ファイルの種類をテキスト文書(*.txt)じゃなく、全ての文書(*.*)にしてから 1)で作ったc:\Octave\MyProjの中に　mytest.m という名前で保存する。　ファイル名は適当でいいが、 拡張子は必ず .m にすること　（テキスト文書で保存すると、 mytest.m.txt っていうワケワカランファイル名に なってしまうので注意） ４）Octaveを立ち上げる。 ５）　cd c:\octave\myproj と打つ。　（１で作ったフォルダ） ６） pwd と打って、ans=c:\octave\myproj って出るのを確かめる ７) dir と打って、　mytest.m が表示されるのを確かめる ８）　mytest と打つ（　拡張子 .m は つけちゃダメ） ９）　動いた？ 10）　本当は source mytest.m と打つべきなんだけどね（これは拡張子をつける。ややこしいね） １１）　mytest とうつと、　mytest という関数を起動しようとして、カレントディレクトリを探すと 　　 mytest.m があるので　それが起動されるというシカケなのだ。 12) 例えば which mytestってやると、それがどこにあるかわかるだろう 13）　たとえば　which plotってやると実は plot も plot.m を遥か深い漆黒の穴から古（いにしえ）の契約により呼び出している事がわかるだろう 14）　こんなもんでどうでしょか？ 271 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/12/30(金) 00:06:35.73 ] >>263 対処ありがとうございます。 どうも複合技でバグっているっぽいですね。 慣れ親しんだgnuplotで使い続けさせていただきます。 それにしてももう少し起動が早くならないですかねぇ。 scicoslab (scilabじゃなくて）　はとても起動が速いのに。 272 名前：octave for windows メモの筆者 mailto:sage [2011/12/30(金) 07:12:56.31 ] >>271 起動はUbuntuなどのUnixではとても速いです。起動がWindows用に最適化 されていないように感じます。 273 名前：octave for windows メモの筆者 mailto:sage [2011/12/30(金) 08:18:47.62 ] >>269 どうもoctavercおかしい、文字列をの''であって、""でないから、\\でなくて\が パスセパレータですね。でも、おかしい octave:1> edit warning: unrecognized escape sequence `\P' -- converting to `P' warning: unrecognized escape sequence `\N' -- converting to `N' error: system: CreateProcess failed -- can't create child process error: called from: error: C:\Octave\Octave3.4.3_gcc4.5.2\share\octave\3.4.3\m\miscellaneous\edit.m at line 212, column 7 とりあえず、Notepad++をOctaveとは別に起動してそこで編集してみてください。 274 名前：octave for windows メモの筆者 mailto:sage [2011/12/30(金) 08:33:41.29 ] >>271 >>263 Cygwin版octaveで、試してみました。FLTKマウスぐりぐりでも、Octaveハング しません。ただし、マウスぐりぐり自体は少し変ですが。 # Cygwin版では、octave-forgeパッケージはauto loadでないですね。 275 名前：octave for windows メモの筆者 mailto:sage [2011/12/30(金) 15:21:23.37 ] >>273 windows 7 64ビットだとNotePad++は C:\Program Files (x86)\Notepad++ にインストールされる。 （）は、認識されるかわからないので、コマンドプロンプトで dir c:\ /X 古いプログラム互換のための短いパス名を調べて、Program Files (x86) は、PROGRA~2であることがわかる。ついてにパスセパレータを\から /に変更して、私の環境では、 EDITOR('C:/PROGRA~2/Notepad++/notepad++.exe'); としてみた。 これで、octave のコマンドプロンプトから、editで NotePad++ が （少なくとも私の環境では）起動するようになった。 276 名前：名無しさん＠お腹いっぱい。 [2011/12/31(土) 10:25:08.41 ] >>270 できました！ すごい。9の時点で動いたので、これからはnotepadにコードを作って そのたびに5から8の動作をして確かめたいと思います。 本当にありがとうございました：D 277 名前：octave for windows メモの筆者 mailto:sage [2012/01/06(金) 15:14:04.04 ] New dependency libraries and toolkits for MinGW build are uploaded. ttp://octave.1599824.n4.nabble.com/New-dependency-libraries-and-toolkits-for-MinGW-build-are-uploaded-td4267845.html 278 名前：名無しさん＠お腹いっぱい。 [2012/01/06(金) 23:34:58.90 ] matlab風に使いたいのでoctave workshopのinstaller付のwin版を ダウンロードしたのですが、インストールはできたのに起動しません。 原因をご存知の方がいればお教えください。 279 名前：octave for windows メモの筆者 mailto:sage [2012/01/07(土) 12:55:10.79 ] >>278 Octave Workshopのoctaveあまりにも古く，使用するのはあきらめてください。 GUI Octaveの方が現実的だと思います。 ttps://sites.google.com/site/guioctave/ GUI OctaveではOctave-3.2.4ですが，頑張ればOctave-3.4.3も使用できます。 ttp://guioctave79.freeforums.org/guioctave-stuck-with-gnu-octave-3-4-2-t108.html （OpenGLのgraphicsが使用できませんが。） 280 名前：名無しさん＠お腹いっぱい。 [2012/01/07(土) 21:22:43.08 ] 278です。 GUI Octaveをインストールしたのですが Error message: Redirect console error: 3 と表示されてしまいます。 たびたび申し訳ありませんが、原因をご存知の片がいればお教えください。 281 名前：octave for windows メモの筆者 mailto:sage [2012/01/08(日) 07:16:25.00 ] >>281 私は，GUI に興味がないので使用したことがないので，お役に立てません。 ここで回答がつかなければ，GUI Octaveは生きているプロジェクトのようです ので敷居は高いですが，GUI Octaveの英語のフォーラムで尋ねられるとよいか と思います。 282 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/01/08(日) 12:10:39.94 ] >>280 Octave-3.2.4をインストールしていない に１票 283 名前：名無しさん＠お腹いっぱい。 [2012/01/08(日) 12:30:30.35 ] Octave-3.2.4をインストールする前にGUI Octaveをインストールすると、 Error message: Redirect console error: 3 が出る。 これは、GUIOctaveの　File/Options/Octave/GNUOctave Executable Files Directory に パスが正しく設定されていないことでもわかる。 Octave-3.2.4を入れているのに、なぜかここにパスが設定されていないときは、 ここに　 C:\Octave\3.2.4_gcc-4.4.0\bin を設定する。 そもそもOctave-3.2.4を入れていないのであれば、 一度GUIOctaveをアンインストールして、Octave-3.2.4を入れる. これが正しく入っていると、GUIOctaveはインストール時に自動的にOctaveパスをセットしてくれるので 特に設定は要らない。 Octave-3.2.4を入れるときには octave for windows メモの筆者さんの www.tatsuromatsuoka.com/octave/jpn/OctaveWinMemo.html#100403 を、ありがたく、よく読んでからな〜 284 名前：octave for windows メモの筆者 mailto:sage [2012/01/09(月) 13:38:37.75 ] 3.6.0がでそうでなかなか，いろんな話題がでてなかなかでませんが， そういうなか 3.5.92のスナップショットがでました。 ビルドできる人はトライしてみてください。 octave.1599824.n4.nabble.com/Octave-3-5-92-snapshot-available-for-ftp-td4270829.html 285 名前：octave for windows メモの筆者 mailto:sage [2012/01/11(水) 08:09:58.85 ] Octave-3.6.0のrc0がでました。 octave.1599824.n4.nabble.com/Octave-3-6-0-rc0-release-candidate-available-for-ftp-td4283677.html Octave-3.6.0のリリースもまもなくでしょう。 286 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/01/11(水) 13:30:25.43 ] FreeBSD 9.0Rが出たので入れてみたメモ octaveはpackagesから入れると起動時にGLIBCXXなにがしというエラーが出ます portsから入れるとmake途中で同じエラーが出ます 原因はglibcのバージョンの問題で LD_LIBRARY_PATH の先頭に /usr/local/lib/gcc46 を入れてください。 一応正常そうに動きました。 287 名前：名無しさん＠お腹いっぱい。 [2012/01/14(土) 20:02:14.37 ] for i=1:1:n t=0:dt:t; figure(i) subplot(1,3,1) plot(t,xx1(i,:)) grid on title('Displacement-time') xlabel('time') ylabel('displacement') subplot(1,3,2) plot(t,xx2(i,:)) grid on title('Velocity-time') xlabel('time') ylabel('Velocity') subplot(1,3,3) plot(t,xx3(i,:)) grid on title('Acceralation-time') xlabel('time') ylabel('Acceralation') hold on end 288 名前：名無しさん＠お腹いっぱい。 [2012/01/14(土) 20:04:41.86 ] 上のコードは自分で作成しました。 もし、n=5であれば、５つの図を出したいのですが、 上記の場合ですと、１つ目の図はプロットされていますが、 ２つ目以降のグラフは、grid とtitleなどが表示された図が表示されるだけで プロットされてません。 なお、xx1やxx2はちゃんと値が入っているので数字的ミスはまったくないです。 forで回すのは、figureを作成するうえでやめたほうがいいんでしょうか？ よろしくお願いします。 289 名前：octave for windows メモの筆者 mailto:sage [2012/01/14(土) 22:22:53.77 ] >>288 それぞれの図を描くとき drawnow が入っていないのが原因なのでは？ 290 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/01/15(日) 12:05:33.07 ] t=0:dt:t; octaveにかけると、この文でエラーになるけど？ 291 名前：octave for windows メモの筆者 mailto:sage [2012/01/15(日) 14:17:54.69 ] >>290 上のコードはこれ単独では動きません。 動かすには，tとdtの設定のほかにも，xx1とxx2の定義が必要ですので。 できれば，単独で動いて動作確認できるコードをあげていただけると 答えがつきやすいと思います。 ところで t=0:dt:t; はループの外に出すでしょう。ループの中にいれるのは， 間違いでないけど意味が不明です。 292 名前：octave for windows メモの筆者 mailto:sage [2012/01/15(日) 18:49:34.46 ] やってみました。たしかにループでなく描画するとグラフかけますね。 バグかもしれません。 Octaveのバージョン，Windows，Mac，Linuxなどお知らせください。 293 名前：octave for windows メモの筆者 mailto:sage [2012/01/15(日) 18:50:29.52 ] >> おっと, windows版octave 3.4.3で試してみました。 294 名前：octave for windows メモの筆者 mailto:sage [2012/01/16(月) 06:47:46.80 ] 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 ] >>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 ] 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 ] >>213 bode線図でこけるのは，control packageのせいではではなく， fltk graphic_toolkit が対数軸でのマウスオペレーションがこける ためのようです。 https://savannah.gnu.org/bugs/?35292 298 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/01/27(金) 09:11:43.11 ] >>296 Scilabも5.4.0-alphaからコマンドウインドウの左右に 補助ウインドウが出るようになりました ただしScilabは左が、変数＆履歴ウィンドウ、右がファイルブラウザ まあでもド素人にMatlab系のプログラミング技術を教えるだけなら FreeMATでいいかもしれないです インストールが簡単だし、dos画面じゃないし、 グラフもgnuplotとかではなくてネイティブだし 299 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/01/29(日) 11:13:02.52 ] >>298 FreeMATって、コントロールパックとかあるんでしたっけ？ 伝達関数とか、それのインパルス応答とか、ボード線図かくやつ。 300 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/01/29(日) 14:56:20.91 ] >>299 無いですね tfもzpもstepもbodeも 制御をやる目的ならscilabの方がいいと思います scilabはmatlabから入ると流儀がわからないかもしれませんが 慣れればどってことないです octaveは限界が高いですけど 初心者がやるにはハードルも高いと思います dos画面だけってのは初心者は面くらうでしょう それで電卓代わりに使う目的ならfreematで十分という話です 301 名前：octave for windows メモの筆者 mailto:sage [2012/01/30(月) 07:48:33.24 ] >>300 OctaveはもともとUnix由来のソフトですし，Windowsの開発陣のリソースは実に 薄い。UnixではTerminalでのソフトはまだ一般的ですし，UnixのTerminalは コマンドプロンプトよりは使いやすい。 #Windows 2000からはコマンドプロンプトが正しい名称であって #ドスプロンプトではないです。まあ，それのほうがつうじますが。 302 名前：octave for windows メモの筆者 mailto:sage [2012/01/30(月) 07:54:27.43 ] >>301 ちょっと舌足らずした。 とはいうもののUnixでもGUIのご時世なので，Octave本体もやっと 真剣にGUIインターフェースをつくりはじめています。 ちとお待ちください。 303 名前：octave for windows メモの筆者 mailto:sage [2012/01/31(火) 12:28:20.15 ] 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 ] 遅ればせながら、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 ] >>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 ] 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 ] 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 ] 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 ] >>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 ] >>309 試してというか、インストールすらまだしてないｗ ちょっとoctave使ってみようかな、ありがとう 311 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/02/14(火) 21:52:48.99 ] 久々です 少し訳の分からないエラーが出たので質問したいのですが… 今サクラエディタ上で次のようなプログラムを作ってました 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 ] >>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 ] >>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 ] >>313 解決しました！ ありがとうございます！ 315 名前：名無しさん＠お腹いっぱい。 mailto:sage [2012/02/21(火) 22:33:37.23 ] 流石に課題丸投げは行けないと思うが、悩みに悩んだ結果積んでしまったので教えてください 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 ] >>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 ] >>316 (ip(k,n)に関して、kが整数ではないからエラーがでる)って書いてありましたね。 すみません。 318 名前：octave for windows メモの筆者 mailto:sage [2012/02/23(木) 12:13:15.01 ] >>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 ] >>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 ] >315-319 要は配列の概念とcontoturfの使い方がわかっていないのかなと思います。 ヘルプを読めばいいのですが，英語ですしね。無理もないかもしれません。 pdfをダウンロードできました。 もう少しいじる必要があると思いますが，基本的にはこれでいいのかな。 あとは自力で頑張ってください。 321 名前：octave for windows メモの筆者 mailto:sage [2012/02/23(木) 18:04:53.50 ] >315-319 要は配列の概念とcontoturfの使い方がわかっていないのかなと思います。 ヘルプを読めばいいのですが，英語ですしね。無理もないかもしれません。 pdfをダウンロードできました。 もう少しいじる必要があると思いますが，基本的にはこれでいいのかな。 課題ですので，あとは自力で頑張ってください。 322 名前：octave for windows メモの筆者 mailto:sage [2012/02/28(火) 15:10:56.60 ] 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 ] まだ，公式には発表されていませんが，Michael氏がMSVC版のoctaveの復活に意欲 を見せておられます。さて，どうなるでしょうか。 324 名前：octave for windows メモの筆者 mailto:sage [2012/02/28(火) 15:18:51.21 ] まだ，公式には発表されていませんが，Michael氏がMSVC版のoctaveの復活に意欲 を見せておられます。さて，どうなるでしょうか。 325 名前：octave for windows メモの筆者 mailto:sage [2012/03/03(土) 16:10:02.07 ] 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 ] いきなりの質問ですみません。 ２つほど質問があります。 １．２段のワードの形式で論文を書いているのですが、 どうしても、ワードの大きさに縮小してしまうと、メモリどころか 線すらはっきり映りません。これをどう解決すればよいのでしょうか？ ２．グラフを描くときに、色ではなく破線や点線はどのようなコードを 打てばよいのでしょうか？ なにとぞよろしくお願いします。特に１番目の質問が深刻です。 327 名前：名無しさん＠お腹いっぱい。 [2012/03/03(土) 23:42:00.21 ] ２番目は解決しました。すみません。 328 名前：octave for windows メモの筆者 mailto:sage [2012/03/04(日) 17:42:22.14 ] >>326 ？？　プロットした図が縮小すると見えなくなるということでしょうか。 ワードだからWindowsですよね。octave-3.2.4でしょうか。 gnuplotの画面をコピーしてワードに貼り付けているということでしょうか。 私はoctave の図をそのまま使わないので良くわかりませんが， print コマンドでemfか何かでセーブした図を使ったほうがいいかもしれません。 329 名前：名無しさん＠お腹いっぱい。 [2012/03/05(月) 20:49:38.55 ] 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 ] >>329 問題ないです 331 名前：octave for windows メモの筆者 mailto:sage [2012/03/07(水) 08:09:54.51 ] 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 ] >>329 補足するとsingle(単精度)の有効桁は10進に治すと約6-7桁なので それ以上の桁数の数値についての保証はなく　 >>330 の方がおっしゃるように，単精度の範囲では問題ないことになります。

---

---

[ 新着レスの取得/表示 (agate) ] / [ 携帯版 ]

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