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

[表示 : 全て最新50 1-99 101- 201- 301- 2chのread.cgiへ]
Update time : 03/17 14:42 / Filesize : 146 KB / Number-of Response : 333
[このスレッドの書き込みを削除する]
[＋板最近立ったスレ＆熱いスレ一覧 : ＋板最近立ったスレ／記者別一覧] [類似スレッド一覧]


↑キャッシュ検索、類似スレ動作を修正しました、ご迷惑をお掛けしました

1 名前：名無しさん＠お腹いっぱい。 mailto:sage [2009/10/28(水) 22:13:01 ]: 知る人ぞ知る強力なデータ解析ツール octave ...

unixを科学研究で使っている香具師には非常に心強い味方のはずだ。

matlab互換を主張しながら、微妙にバージョンアップを続ける
octaveについて語ろう。

そのPart 2です

前スレ

データ解析ツールoctaveを語ろう
pc12.2ch.net/test/read.cgi/unix/1086016412/
128 名前：octave for windows メモの筆者 mailto:sage [2011/06/25(土) 06:29:16.81 ]: >>124
CUDAについては，OctaveはGPLですので，
jwe氏が言うことは仕方のないことなのでしょう。
OctaveがGPLから他のライセンスに変えるか,GPLが変わるか
CUDAがかわるかどちらかでないといけないのでしょうが，
まあ実現はしないのだと思います。

ところで，MySQLは，
www-jp.mysql.com/about/legal/licensing/faq.html
にあるように問題がないように思います。
だからOctave-forgeのdatabaseパッケージでMySQLはサポートされている
のだと思っています。
129 名前：octave for windows メモの筆者 mailto:sage [2011/06/25(土) 06:35:35.07 ]: >>128
ただし，再配布せず自分でやる分には問題がないと思います。
以前GotoBLASが，GPL非互換だったとき自分でソースからビルドして，
octaveのビルド時に組み込んで使用することはできました。

かなり，億劫なことになりますが，
scinesur.wordpress.com/2009/11/12/cuda-on-octave-3-x-in-linux-64-bits/
なこのページもあることですし，自力でやることはできるのかなと思います。
130 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/06/25(土) 08:57:30.61 ]: >>126
mingw版もcygwin版も(どちらも3.2.4)

octave:5> pkg install database-1.0.4.tar.gz

で入れようとしても、内部のmakeでこけますね

tar.gzを直接開いて./configureしてみると
sqlite3でwarningが出るし
結果がこれで

octave-forge is configured with
octave: octave (version 3.2.4)
mkoctfile: mkoctfile for Octave 4
swig: pre-generated
postgres: postgresql/libpq-fe.h not found
mysql: mysql/mysql.h not found
sqlite3: yes
odbc: yes

postgresqlやmysqlはライブラリが必要みたいだし
makeもsqliteの段階でこけるので、unix系なら動くのかもしれませんが
Windowsでは無理っぽいですね
131 名前：octave for windows メモの筆者 mailto:sage [2011/06/25(土) 12:23:55.27 ]: >>130

そうですね。3.0.xのころdatabaseパッケージを四苦八苦していれましたが，
postgres　と　mysqlはあきらめてsqlite3とodbcだけでやったことがあります。

CygwinならたしかpostgresqlがOctaveを意識せずにあったような気がしますが，
それを使ってpkg installをしたことがないので，なんともいえません。
Cygwinの方がUnixに近いのでライブラリをかき集めれば何とかなるかも
知れません。
132 名前：octave for windows メモの筆者 mailto:sage [2011/06/29(水) 02:59:55.51 ]: MinGWの体験版のOctave-3.2.4バイナリを作成し，配布してます。
www.tatsuromatsuoka.com/octave/Eng/Win/

あくまでも体験版で，
Octave-forgeがはいっていませんし，バグもあると思います。
インストーラーなる気の利いたものもなく，zipファイルを解凍して
適当なフォルダにおきます。
ビルドに必要なツールも一応いれたつもりなので，自分でビルドしてみたい
方はどうぞ。ただし，あまり丁寧な説明はしていません。
133 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/07/04(月) 06:29:50.28 ]: どなたか賢い方、助けて下さい…

Benjaminさんの
Octave-3.2.4_i686-pc-mingw32_gcc-4.4.0_setup.exe
で、コマンドラインからのpath指定が効かないんですが、なんとかなりませんでしょうか？

> octave.exe --path c:\
> octave.exe --path c:\\
> octave.exe --path c:/
> octave.exe --path c://
> octave.exe --path /c/

いずれもダメでした。
Linuxなら特に問題ないのですが…
134 名前：octave for windows メモの筆者 mailto:sage [2011/07/04(月) 15:21:59.89 ]: >>133 Cygwin版ならちゃんと動くんでバグかもしれません。
135 名前：octave for windows メモの筆者 mailto:sage [2011/07/04(月) 15:28:18.58 ]: >>133

>octave --eval "cd c:\\" --persist
で代用できませんか？
136 名前：octave for windows メモの筆者 mailto:sage [2011/07/04(月) 15:38:46.46 ]: >>134, 135　すみません勘違いです。

マニュアルによれば
--path path
-p path
Add path to the head of the search path for function files.
The value of path specified on the command line will override any
value of OCTAVE_PATH found in the environment, but not any commands in
the system or user startup files that set the internal load path
through one of the path functions.

となっていて，Cygwin版でもMinGW版でも OCTAVE_PATH には，なにも反映され
ない。Linuxでの期待通りの動きとはどんなものですか？

Linuxもってないのでわからないので。
137 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/07/04(月) 23:18:23.00 ]: >>134-136
試して下さって、大変ありがとうございます！！

Linux版なら、
$ octave -q --path /home/user/work
octave:1> path
.
/home/user/work　（←このフォルダが存在する場合のみ）
/usr/lib/octave/3.2.4/site/oct/i686-pc-linux-gnu
/usr/lib/octave/site/oct/api-v37/i686-pc-linux-gnu
...
となります（恐らくご想像通りでしょうが…）。

ちなみに、
$ OCTAVE_PATH=/home/user/work octave -q
でも同義です。

こちらでも迂回策を探してはいるのですが、環境変数も効きませんし、なかなか…
138 名前：octave for windows メモの筆者 mailto:sage [2011/07/05(火) 01:47:59.89 ]: >>137
octave --eval "addpath C:\\" --persist
ではいかがでしょうか。
139 名前：octave for windows メモの筆者 mailto:sage [2011/07/05(火) 01:56:50.28 ]: >>137 蛇足です。Cygwin版のoctaveでは，おっしゃられている動作をする
ようです。MinGWならではのバグのようです。

職場のパソコンで自分でビルドした3.4.2（MinGW）でも，動作しなければバグ
報告をします。
140 名前：octave for windows メモの筆者 mailto:sage [2011/07/05(火) 06:45:34.07 ]: >>139 職場のパソコンでoctave 3.4.2と3.2.4(Benjamin版)で試してみました。

octave --path c:\
＊＊＊＊＊＊＊＊＊＊＊
3.2.4
.
C:\Programs\Octave\3.2.4_gcc-4.4.0\share\octave\packages\zenity-0.5.7
　（以下略）
＊＊＊＊＊＊＊＊＊＊＊＊＊＊
3.4.2
.
C:
C:\octave\3.4.2\lib\octave\3.4.2\site\oct\i686-pc-mingw32
:
(以下略)
となるのですが
octave --path C:\octave

Octave's search path contains the following directories:

.
C:\octave
C:\octave\3.4.2\lib\octave\3.4.2\site\oct\i686-pc-mingw32
:
　（以下略）
＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊
141 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/07/05(火) 07:58:49.52 ]: >>138
ありがとうございます、現在そちらの方向で検討中です。

>>140
え？なぜ3.4.2はそんな残念な動作に…
3.2.4（Benjamin版）ではC:\Octaveを指定しても駄目なようですね

こちらの職場PCにはBenjaminさんのメモにしたがってビルドした
3.2.4がありますので、そちらも試してみます。
142 名前：octave for windows メモの筆者 mailto:sage [2011/07/05(火) 08:11:41.89 ]: >>141
140　は長くて，2chにおこられたので，後から続きを書こうとしてわすれていました。
どーもドライブのルートディレクトリの場合は失敗して，C:\とならず，C:
となるのですが，C:\octaveのような場合は大丈夫です。

一応，development branch　でも試すために，現在ビルド中です。　
143 名前：octave for windows メモの筆者 mailto:sage [2011/07/05(火) 09:51:31.50 ]: >>142 予想通りDevelopment Branchのソースからビルドしても，
octave --path c:\
のようなルートディレクトリの指定はバグってますね。
時間があるときにbug trackerに登録します。
（直るかどうかわかりませんが。）
144 名前：octave for windows メモの筆者 mailto:sage [2011/07/05(火) 10:18:47.77 ]: >>143 bug tracker に登録しておきました。
https://savannah.gnu.org/bugs/?33714
145 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/07/05(火) 23:05:49.70 ]: >>143
大変感謝です…。

ちなみに、Benjaminさんメモに従ってビルドした3.2.4ですが、
MSYS環境ではpath指定が通りました。

コマンドプロンプトからはlibgcc_s….dllが足りないとか言われて
実行できなかったので、未検証です。
（MSYS環境以下を探しても見つからなかったので諦めてましたが、
よく考えたらMinGW環境以下にありそうですね…忘れてました）
146 名前：octave for windows メモの筆者 mailto:sage [2011/07/06(水) 01:46:53.64 ]: >>145
Msys 環境下で指定がうまくいくということは，パスのハンドリングが
octave --path でUnix likeになっているからかもしれません。
なお，libgcc_s….dllなどたらない，dllは，Benjaminさんのビルドのbin
ディレクトリからコピーして，自分でビルドしたoctave.exe,
octave-3.2.4.exeのあるディレクトリにペーストすればいいです。

または，バッチファイルなどでBenjaminさんのビルドのbinディレクトリに
パスを通すという手もあります。
147 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/07/06(水) 18:47:15.59 ]: >>146
コマンドプロンプト環境下で試してみました。
dllはOctaveのビルドに使用したc:\mingwのbinにパスを通して
対処してます。

octave --path c:\msys -> OK
octave --path c:\ -> OK（ちゃんとc:\と表示される）
…あれ？

Benjaminさんメモに従うとコンパイルは自動で流れるので、
できるバイナリはインストーラで配布されているものと
同じはず、だと思うんですが…？

なお、本来の目的はOctave起動後にaddpathすることで達して
おりますので、これ以上追いかける意味は無くなってしまったのですが…。
ご報告まで。
148 名前：octave for windows メモの筆者 mailto:sage [2011/07/06(水) 21:34:11.72 ]: >>147
ご報告ありがとうございました。
回避法があるとはいえ，不思議ですね。
149 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/08/03(水) 16:34:09.93 ]: 質問いいでしょうか？
画像をfft2して、その後、元の画像に戻るかiff2でチェックしたいのですが、どうもうまくいきません・・・
プログラムのどこが悪いのでしょうか？

Iin=imread('BOAT.bmp');% 画像の入力
figure,imagesc(Iin,[0 255]);% 画像表示
colormap(gray);
axis image;

%FFT変換
Iin_d=double(Iin);
O_fft=fft2(Iin_d);% FFT
O_fft=fftshift(O_fft);
O_fft=(abs(O_fft));
O_fft=10*log10(O_fft);% FFT結果のlog表示
figure,imagesc(O_fft,[0 MAX]);% 画像表示
colormap(gray);
axis image;

%逆FFT
O_fft2=(abs(O_fft2));
O_fft2=10*log10(O_fft2);
O_fft2=fftshift(O_fft);
O_fft2=real(ifft2(O_fft2));
figure,imagesc(O_fft2,[0 MAX]);% 画像表示
colormap(gray);
axis image;

どなたかお願いします
150 名前：octave for windows メモの筆者 mailto:sage [2011/08/08(月) 09:37:41.35 ]: >>149
画像処理は詳しくありませんが，演算が逆演算になっていないような気がします。
151 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/08/08(月) 14:03:14.87 ]: abs()をつけてる時点で虚数次元が飛んでますからね
152 名前：octave for windows メモの筆者 mailto:sage [2011/08/08(月) 16:59:04.63 ]: >>149 >>150
自信なし：一応画像で確認した。

Iin_d=double(Iin);
O_fft=fft2(Iin_d);% FFT
O_fft2=real(ifft2(O_fft)); %ifft
figure,imagesc(Iin,[0 255]);

fftした結果を演算されていますが，もとのイメージのフーリエ空間での
画像に演算をされておられるので情報が変わってしまっています。
したがって，ifftで戻るはずは原理的にないような気がします
153 名前：octave for windows メモの筆者 mailto:sage [2011/08/08(月) 17:55:52.63 ]: >>152
訂正

Iin_d=double(Iin);
O_fft=fft2(Iin_d);% FFT
O_fft2=real(ifft2(O_fft)); %ifft
figure,imagesc(O_fft2,[0 255]); %ここを訂正
154 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/08/09(火) 15:43:46.03 ]: むしろこんな感じ

sigma=40;
for XX=1:128,
for YY=1:128;
ZZ(XX,YY)=63*exp(-1*(XX-64)^2/(2*sigma^2))*exp(-1*(YY-64)^2/(2*sigma^2));
end;
end;
ZZ1=fft(ZZ);
ZZ2=ifft(ZZ1);
image(ZZ)
colorbar
pause(2);
image(abs(ZZ1))
colorbar
pause(2);
image(abs(ZZ2));
colorbar
pause(2);

要点は、imageで表示するとき以外はabsとかを使わない
155 名前：octave for windows メモの筆者 mailto:sage [2011/09/02(金) 10:37:46.27 ]: Blas としてOpenBLASかなり使えます。
https://github.com/xianyi/OpenBLAS

なお，Cygwin 上にmingw64のクロスコンパイラが入りました。
mingw64は32 bitもサポートしていますので，windows 用には，
それを使います。
（setup.exeでインストールしてください。）

makeは以下のようにしてください。
make CC=i686-pc-mingw32-gcc FC=i686-pc-mingw32-gfortran OSNAME=WINNT TARGET=CORE2
TARGETは各自のCPUにあわせてください。
できた
libopenblas_(cpu id)-r0.1alpha2.2.dll
octaveのあるbinフォルダにコピペして，
octave-3.2.4 mingwなら blas.dllとリネームすればいいです。
（octave-3.4.2 ならlibblas.dllに）
156 名前：octave for windows メモの筆者 mailto:sage [2011/09/05(月) 10:46:31.60 ]: Octave-3.4.2-mingw + octaveforge pkgs が出ました。
old.nabble.com/Octave-3.4.2-mingw-%2B-octaveforge-pkgs-to32394771.html
157 名前：octave for windows メモの筆者 mailto:sage [2011/09/07(水) 17:24:14.86 ]: Gnu Octave windows版のOctave-forge packageがアップデートされました。
old.nabble.com/Re%3A-Octave-3.4.2-mingw-%2B-octaveforge-pkgs-p32413743.html
158 名前：octave for windows メモの筆者 mailto:sage [2011/09/15(木) 08:01:26.94 ]: Gnu Octave windows版のOctave-forge packageがアップデートされました。
octave.1599824.n4.nabble.com/Updated-Octave3-4-2-Mingw-for-testing-was-Mingw-Octave-3-4-2-blank-gnuplot-printouts-td3814047.html
159 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/09/22(木) 00:01:20.14 ]: memory exhaustエラーが、少し大きめのデータで出てしまうのが、
改善されるといいなぁ。
160 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/09/22(木) 09:47:51.52 ]: >>159
64bit版ならそのエラーが簡単には出ないって話ですけどね
ただし正式なバイナリーがリリースされてない
161 名前：octave for windows メモの筆者 mailto:sage [2011/09/22(木) 14:14:55.62 ]: >>159 >>160

ttp://www.gnu.org/software/octave/doc/interpreter/Compiling-Octave-with-64_002dbit-Indexing.html

にしたがってビルドするしかないようですね。
Debian 上での例が載ってます。
162 名前：octave for windows メモの筆者 mailto:sage [2011/09/22(木) 17:00:01.71 ]: >>159 >>160
Scilab なら64bit版があるようです。

ttp://www.scilab.org/products/scilab/download

でも，文法かなり違うからなあ。
163 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/09/26(月) 13:34:26.79 ]: >>162
Scilabの32bitも64bitもメモリー管理としては大差無いようです

Windows7 64bit に4GBのメモリーを積んだ条件で
Octave 3.2.4で

AA=ones(12000,10000);

とやると通りますが

AA=ones(13000,10000);

では error: memory exhausted と出ます。

Scilab 5.3.3 では32bitでも64bitでも

stacksize('max');

を実行した後に

AA=ones(13000,10000);

は通りますが

AA=ones(14000,10000);

はエラーが出ます。
倍精度が16bitだとして(複素数の実部と虚部にそれぞれ8bit倍精度)
2GBぐらいに壁があるようです。
164 名前：octave for windows メモの筆者 mailto:sage [2011/09/26(月) 17:59:23.61 ]: >>163
私もやってみました。
-->AA=ones(13600,10000);
!--error 17
stack size exceeded!

Use stacksize function to increase it.
Memory used for variables: 135007946
中間的なメモリが必要です: 136000002
Total memory available: 268435455

Scilabの場合，スタックサイズの上限に問題があるようですね。
165 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/09/27(火) 09:00:38.67 ]: >>163
倍精度複素数16bitではなくて16Bytesの間違いでした

全部の変数で2GB未満しかアロケート出来ないっぽいです
2^32は42億9496万7296でしたか、つまり4Gですが
CPUは2の補数をとった2Gまでしかアクセスできないのでしょう
166 名前：octave for windows メモの筆者 mailto:sage [2011/09/27(火) 10:26:32.43 ]: >>165
ttp://www.gnu.org/software/octave/doc/interpreter/Compiling-Octave-with-64_002dbit-Indexing.html
をよく読むとArray Indexとなる整数のバイト数を64 bitの体系にあわせて，
ライブラリまで含めて構築しなおす必要があるようです。
（たとえばFortran のライブラリなら（たとえばLAPACK），
interger のデフォルトのバイトを8バイトにする。）
Scilab 64 bitは，そこまで考えてビルドされていなのではと推測します。
167 名前：octave for windows メモの筆者 mailto:sage [2011/09/30(金) 10:42:53.78 ]: Octave 3.4.3-rc0 release candidate available for ftp

ttp://octave.1599824.n4.nabble.com/Octave-3-4-3-rc0-release-candidate-available-for-ftp-td3857294.html

だそうです。
168 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/10/01(土) 22:47:50.57 ]: Octaveの連想配列はどう書けば良いのですか?
データに番号がついていて連続していない状態で、他のデータから連続していないデータを参照しているので
内部で番号をつけ直しています。つけ直しの時に総ナメのような感じになってしまうので連想配列が使えたら
と思いまして質問しました。
169 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/10/02(日) 01:13:05.09 ]: >>168
構造体で、それっぽいことはできたはず。
その用途でしたら、スパース行列でもよさげ。
どっちが速いかは未検討。
170 名前：168 mailto:sage [2011/10/02(日) 09:20:50.75 ]: >>169
a.idみたいな感じで丁度構造体を使っているので構造体を使った連想っぽい処理をもう少し具体的に知りたいです
どんな感じで使えば良いのですか?
171 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/10/02(日) 14:17:52.86 ]: >>169
問題がややあいまいなので、もうちっと具体的にどうぞ。
情報の後出しは、無しということで。
172 名前：168 mailto:sage [2011/10/02(日) 18:24:07.84 ]: どこまで必要な情報なのか・・・書ききれるのか分かりませんが書いてみます
有限要素法のプログラムを勉強のために作っているのですが、データが節点、
要素等番号がついています。NASTRANのデータ形式をGRIDとかCHEXAとか
を使っています。(手を抜いて継続行を一行にしていたりしますが・・・)
これらの番号(ID)は連続した数字ではなく、プログラム内部で参照して番号を
関連付けるしかないのですが、現状総ナメで番号(ID)を内部の連続した配列番号
に関連付けています。変数はelem(:).id、node(:).idがデータの番号で関連付
けした変数はelem(:).no、node(:).noとかにしています。主にelemがnode
を参照するのでnode(:).idを連想配列風の処理で内部の配列番号に変換したいのです
読みにくいと思いますが、構造体を使った連想配列風の処理を教えていただけない
でしょうか？
173 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/10/02(日) 23:48:46.22 ]: >>172
もっとスマートな方法もありそうですが…
期待外れだったら、すいません。

案１（スパース行列）
node_id = [1 10 100 1000 10000];
node_id_t = sparse(0,0);
node_id_t(node_id) = 1:length(node_id);
% node_id_t(10000) -> 5

案2（連想配列？）
node.id_1 = 1;
node.id_10 = 2;
node.id_100 = 3;
node.id_1000 = 4;
node.id_10000 = 5;
% getfield(node,['id_' num2str(10000)]) -> 5

ちなみに、FEMをOctaveで、ということでしたら、OpenFEMというのも参考になるかもしれません。
Matlab用で、コードがめちゃくちゃ読みづらくて、かなりbuggyで信用なりませんが、
各種アルゴリズムの実装を知るには非常に有用です（個人的に、固有値解析のとこは感動しました）。
174 名前：168 mailto:sage [2011/10/04(火) 20:03:02.40 ]: >>173
お返事有り難う御座います。
処理速度を早くするためにはスパース行列が良さそうですね、
処理手順を多少弄らなければいけなさそうですが少し考えてみます。
案2の方も少し考えてみます。（使えそうな使えなさそうな・・・）

OpenFEMの情報有り難う御座います。少しググってみます。
175 名前：octave for windows メモの筆者 mailto:sage [2011/10/11(火) 10:17:19.15 ]: GNU Octave 3.4.3 Released　
ttp://octave.1599824.n4.nabble.com/GNU-Octave-3-4-3-Released-td3890670.html
176 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/10/11(火) 21:28:31.48 ]: >>175
毎度ありがとうございます。
目玉機能は何かあります？
(simulinkがついに出来たとか :-P )
177 名前：octave for windows メモの筆者 mailto:sage [2011/10/12(水) 16:12:10.66 ]: >>176
残念ながら喜ぶような機能の追加はなく，バグフィックスのようです。
3.4はこれでおしまいで，次は3.6になります。
GUIのインターフェースがつくかも知れませんが，はっきりとしません。
178 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/10/12(水) 20:10:34.63 ]: >>177
あ、それは残念。
純正のOctave Workshopと、ネイティブEXCELリーダーが付けばいいなぁ。
179 名前：octave for windows メモの筆者 mailto:sage [2011/10/13(木) 11:33:00.99 ]: >>ネイティブEXCELリーダー
Octave-forgeのio packageである程度のことはできるようです。
私はつかったことがないのでよくわかりません。
ttp://octave.sourceforge.net/io/overview.html

>>純正のOctave Workshop
Qt BaseのGUIが検討されていてもしかしたら，3.6で実装
（たぶんExperimental）されるやもしれません。
Windows版については，最近活動がストップしているBenjaminさんに
変わる人が出てくればいいのですが。
私のスキル，体調，仕事上の都合でそこまでついていけないのが実情です。
180 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/10/13(木) 13:02:18.29 ]: scilabでは伝達関数を定義するのに

s=poly(0,"s")
sys=1/(1+s)

ってできますが
Octaveでは
sys=tf( [1] , [1 1] );
ですよね。
scilabみたいに記法する方法はないですかね？
181 名前：octave for windows メモの筆者 mailto:sage [2011/10/13(木) 17:08:44.70 ]: 制御系には全く無知なのではずしていたらすみません。
Octave 3.4に対応した
Octave-forgeのcontrol packageでは，
ttp://octave.sourceforge.net/control/function/tf.html

Exampleに
octave:1> s = tf ("s");
octave:2> G = 1/(s+1)

Transfer function "G" from input "u1" to output ...
1
y1: -----
s + 1
とあります。
Windowsの場合
ttp://www.tatsuromatsuoka.com/octave/jpn/OctaveWinMemo.html#110906
こちらをご覧ください。
182 名前：octave for windows メモの筆者 mailto:sage [2011/10/13(木) 17:10:32.55 ]: >>181
1
y1: -----
s + 1
が崩れて表示されていますね。すみませんm(__)m。
183 名前：octave for windows メモの筆者 mailto:sage [2011/10/13(木) 17:18:17.15 ]: >>181
何度もすみません。
***********＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊
Windowsの場合
ttp://www.tatsuromatsuoka.com/octave/jpn/OctaveWinMemo.html#110906
こちらをご覧ください。
***********＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊
は
WindowsでOctave-3.4.2でcontrol packageなどのOctave-forgeを使いたい場合
ということです。舌足らずな説明ですみません。
184 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/10/13(木) 23:33:20.68 ]: 早速ありがとうございます。
使ってみます。

ginputのバッティングなんかも解消できているといいのですが。
185 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/10/14(金) 15:57:12.92 ]: >>180
OctaveのはMatlab互換なので馴れてしまえば後々で便利です

例えば化学の一次反応を
(d/dt)x=k(x0-x)
と置きましょう。xは反応量、x0は初期量、kは反応速度定数(毎秒)です。
式を展開すると
(1/k)(d/dt)x=x0-x
x0=(1/k)(d/dt)x+x
これをラプラス変換すると
X0=(s/k)X+X=(s/k+1)X
X0=((s+k)/k)X
伝達関数G(s)は(伝達関数はOutput/Input)
G(s)=X/X0=k/(s+k)
これを普通にラプラス逆変換すると
x/x0=k exp(-kt)
ですが、Octaveで計算すると、k=0.02 として
Octave:4> k=0.02; G=tf([k] [1,k],0);
Octave:5> [y,t]=step(G);plot(t,y);
であっさり描けてしまいます
時定数は1/k=1/0.02=50秒で、図でも50秒で1-exp(-1)=63%にまで上がってるはずです。
186 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/10/15(土) 13:48:46.19 ]: >>185
>>180のヒトは

伝達関数G(s)は(伝達関数はOutput/Input)
G(s)=X/X0=k/(s+k)
で

Octave:4> k=0.02; G=tf([k] [1,k],0);
を
Octave:4> k=0.02; G=k/(s+k)
で書ければいいなぁ　

って言っているんだと思います。

Octave3.4だと出来るのか。
187 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/10/15(土) 14:33:56.65 ]: >>186
それはわかってるつもりなのですが

伝達関数は計算してから入れればいいという考え派なので

しかし制御の本がScilabばっかりなのは
伝達関数の計算なんかめんどい＆scicos(5.2以降ならxcos)で制御を勉強したい
という人が多いのでしょうね
188 名前：octave for windows メモの筆者 mailto:sage [2011/10/17(月) 07:33:56.72 ]: >>186
>Octave3.4だと出来るのか。
正確にいうと最近のControl package でできるようになり，
それにはoctaveが3.4以上である必要があるということのようです。

ところで3.4.2 for winodws用のOctave-forge Packageが更新されました。
　wiki.octave.org/wiki.pl?OctaveForWindows
を参考にしてください。Control package もバージョンが
あがっているようです。
189 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/10/17(月) 20:24:04.79 ]: >>188
すいません。
octave.sourceforge.net/
のWindows installer (Octave and Octave-Forge)
の方は更新されないんですか？
190 名前：octave for windows メモの筆者 mailto:sage [2011/10/18(火) 08:54:52.33 ]: >>189
される見込みは，残念ながら今のところありません。
MaintainerのBenjamin Lindnerさんが5月ごろから，svnの作業ツリーやML
などから，姿を消しておられます。

私は，ビルドしかできないので3.4.2を自分のWebにのせたら，
Nitzan Araziさんが7zipによるアーカイブと octave-forge のパッケージを
(7 zip)作ってくださいました。
191 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/10/19(水) 20:34:16.86 ]: lsode関数で微分方程式を解いて時間経過を見るとき
変数xが負の値をとる場合にはゼロになるように条件を加えたいのですが
そういったことは出来るのでしょうか。
192 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/10/20(木) 12:35:32.94 ]: >>191 マニュアルを読む限りではそのようなオプションはないと思います。
数学的な工夫でできるかもしれませんが。すみませんがわかりません。
193 名前：名無しさん＠お腹いっぱい。 [2011/10/21(金) 23:07:54.32 ]: >>188
Octave3.4.2_gcc4.5.2\share\octave\packages\control-2.1.54
の中身を見ているのですが、
tf2ss とか　sysout とか、その変の関数がゴッソリと無くなっているのは
なんか理由があるんでしょうかね？

s=tf("s") で　G=1/s , H = 1/(s^2+1) , GH=G*H
っていう記述が出来るようになったのは有りがたい反面、使い慣れた関数が
無くなっているのはちょっと寂しい。
(matlabにはあるのに）
194 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/10/24(月) 20:16:09.57 ]: >>193 本家のoctave-dev ML (Octave-forgeのメーリングリスト）で聞くしか
ないでしょう。たしかMaintainer が最近変わったようなので。
195 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/10/24(月) 21:27:34.49 ]: >>193
tf2ssはssでイケルっぽい
196 名前：octave for windows メモの筆者 mailto:sage [2011/10/27(木) 07:58:32.39 ]: Mingw Octave-3.4.3 binaries for testing on windows - Updated
octave.1599824.n4.nabble.com/Mingw-Octave-3-4-3-binaries-for-testing-on-windows-Updated-td3942303.html
197 名前：名無しさん＠お腹いっぱい。 [2011/10/29(土) 15:09:02.44 ]: ↑のやつですが、Win7 x64でcommunications-1.1.0がmakeでエラーを吐いてインストールできません。
communicationsを使いたいので、どなたか情報を持ってる方お願いします

pkg install signal image optim miscellaneous struct specfun -auto -forge -global　（これはOK）
pkg install communications -auto -forge -global
ov-galois.o:C:\Users\hoge\AppData\Local\Temp\oct-2\communications\src/ov-galois.cc:556: undefined reference to `H5Gcreate'
ov-galois.o:C:\Users\hoge\AppData\Local\Temp\oct-2\communications\src/ov-galois.cc:567: undefined reference to `H5Screate_simple'
(snip)
collect2: ld returned 1 exit status
make: *** [gf.oct] Error 1
'make' returned the following error: make: Entering directory `/tmp/oct-2/communications/src'
(snip)
make: Leaving directory `/tmp/oct-2/communications/src'
error: called from `pkg>configure_make' in file C:\Octave\Octave3.4.3_gcc4.5.2\share\octave\3.4.3\m\pkg\pkg.m near line
1339, column 9
error: called from:
error: C:\Octave\Octave3.4.3_gcc4.5.2\share\octave\3.4.3\m\pkg\pkg.m at line 790, column 5
error: C:\Octave\Octave3.4.3_gcc4.5.2\share\octave\3.4.3\m\pkg\pkg.m at line 361, column 9
198 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/10/30(日) 05:18:24.99 ]: >>197
回線が細いので確認できませんが，
Octave3.4.3_gcc4.5.2_20111025.7zにくわえて
Octave3.4.3_gcc4.5.2_pkgs_20111026.7zもいれて
pkg rebuild -auto
でいけませんか。
不要なものはpkg uninstallで消す。
Octave3.4.3_gcc4.5.2_pkgs_20111026.7zにすでにoctave-forgeは入っていて，
それを再認識させるには pkg installでなく，pkg rebuildを使います。
199 名前：名無しさん＠お腹いっぱい。 mailto:sage [2011/10/31(月) 01:21:57.38 ]: >>199
すみません、最初にpkg rebuild -autoしたのにどうしてもダメだと思ってましたが
octaveの再起動が必要なことに今の今まで気づいてませんでした；；
どうもお騒がせしました。
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