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

1 名前：名無しさん＠お腹いっぱい。 [04/06/01 00:13] 知る人ぞ知る強力なデータ解析ツール octave ... unixを科学研究で使っている香具師には非常に心強い味方のはずだ。 matlab互換を主張しながら、微妙にバージョンアップを続ける octaveについて語ろう。 166 名前：名無しさん＠お腹いっぱい。 mailto:sage [2005/05/02(月) 18:59:44 ] FreeBSDではもうoctaveは使えません 167 名前：名無しさん＠お腹いっぱい。 mailto:sage [2005/05/02(月) 19:32:52 ] >>166 (´･ω･`)ｼｮﾎﾞｰﾝ 168 名前：名無しさん＠お腹いっぱい。 [2005/05/03(火) 20:28:50 ] >>166 マジ？ 169 名前：名無しさん＠お腹いっぱい。 mailto:sage [2005/05/08(日) 07:09:16 ] だめじゃん 170 名前：名無しさん＠お腹いっぱい。 mailto:sage [2005/05/08(日) 18:22:22 ] 焦って/usr/ports/distfilesを見たけどoctaveのファイルは残ってなかった 誰かFreeBSD5.4R用のi686バイナリパッケージくれよ 171 名前：名無しさん＠お腹いっぱい。 [2005/05/08(日) 19:39:29 ] こういう時こそportdowngradeの出番では 172 名前：名無しさん＠お腹いっぱい。 mailto:sage [2005/05/10(火) 12:26:18 ] FreeBSDのportsってたまに Octave が Broken になりますよね． またしばらく待っていると復活しますよ． 173 名前：名無しさん＠お腹いっぱい。 mailto:sage [2005/05/10(火) 13:21:18 ] >171 joke かと思ったら本当にあった... sysutils/portdowngrade 174 名前：名無しさん＠お腹いっぱい。 mailto:sage [2005/05/13(金) 00:07:36 ] mahoﾀｿがその内直してくれるんじゃね 175 名前：名無しさん＠お腹いっぱい。 mailto:sage [2005/05/30(月) 22:52:46 ] なおった？ 176 名前：名無しさん＠お腹いっぱい。 mailto:sage [2005/06/01(水) 10:16:48 ] 少なくとも一昨日の段階ではまだ直ってなかったけど... mahoタソたのむ 177 名前：名無しさん＠お腹いっぱい。 mailto:sage [2005/06/01(水) 10:51:19 ] 5.x以降なら直ってるみたい。 4.xはBrokenだけど。 178 名前：名無しさん＠お腹いっぱい。 [2005/06/01(水) 12:51:40 ] 漏れも使っているので、メンテ頼む mahoさん 179 名前：名無しさん＠お腹いっぱい。 mailto:sage [2005/06/01(水) 14:04:37 ] matlab だと、 x{1} = 1; x{2} = 2; x{1} = []; で、x の最初の cell を消去できますが、octave ではどうするんでしょうか? 同様に structure array の要素の消去の仕方も、[] を代入では駄目でした。 180 名前：179 mailto:sage [2005/06/01(水) 14:11:47 ] すみません。嘘付きました。 cell array の要素の消去は、[] を代入でできるんです。 しかし、structure array の要素の消去が分かりません。 181 名前：179 mailto:sage [2005/06/01(水) 14:36:59 ] 解決しました。 web では見つけられませんでしたが、よく考えれば分かりました。 for i=1:10; x(i).f=i; end とし、 x(7) の削除は、 x = x([1:6 8:10]); ですね。 182 名前：名無しさん＠お腹いっぱい。 [2005/06/03(金) 07:45:41 ] 今、FreeBSD で portsから octave を入れたのだけど 微妙に良くなっていっているよね octave メンテや改良ってまだどんどん行われているの? 183 名前：名無しさん＠お腹いっぱい。 [2005/06/03(金) 09:17:17 ] FreeBSD 5.3でportsからoctaveを入れようとしたけれど、 ===> octave-2.1.63 is marked as broken: Does not compile on FreeBSD >=5.x って出て、インスコできなかったよ。 これってFreeBSD 5.x以上では使えないってことなの？ 184 名前：名無しさん＠お腹いっぱい。 mailto:sage [2005/06/03(金) 10:01:00 ] >>183 portsを新しくしてください。 185 名前：183 [2005/06/03(金) 14:38:31 ] >>184 portsは新しくしたはずなんだが…ｵｶｼｲﾅ 186 名前：名無しさん＠お腹いっぱい。 [2005/06/07(火) 19:33:20 ] ttp://adlib.rsch.tuis.ac.jp/~akira/unix/octave/update.html に2.1.Xから多次元配列がサポートされたって書いてあるけど，具体的にどう使えるのでしょうか？ ChangeLogとかNEWSとかにも multi とか　Arrayとか検索しても引っかからないし， zeros(3,3,3)とかも使えないし…． どなたか教えてください． 187 名前：名無しさん＠お腹いっぱい。 mailto:sage [2005/06/07(火) 20:35:35 ] octaveじゃないけど math/sdpaも壊れてるような 188 名前：名無しさん＠お腹いっぱい。 [2005/06/09(木) 00:35:11 ] >>186 バージョンは？ 2.1.50ではだめだけど、2.1.60以降ではzeros(3,3,3)でいける。 189 名前：名無しさん＠お腹いっぱい。 mailto:sage [2005/06/09(木) 22:07:53 ] >>188 2.1.71にしたら出来ました！ありがとうございました！ 190 名前：名無しさん＠お腹いっぱい。 [2005/06/21(火) 11:04:05 ] imageコマンドで起動される外部画像ビューアーが、インストールされた計算機によって、 imagemagicだったり、xvだったりするのですが、この違いはなぜ生じるのでしょうか。 また、明示的にどちらかを使いたい場合、どうすればよいでしょうか。 191 名前：名無しさん＠お腹いっぱい。 mailto:sage [2005/07/06(水) 23:18:05 ] printf('kita---(.v.)-----!'); sleep(10); これが１０秒後まで表示されないんですよね どうしようもないですか？ 192 名前：名無しさん＠お腹いっぱい。 mailto:sage [2005/07/16(土) 22:05:40 ] >>191 printf('kita---(.v.)-----!'); fflush(stdout);sleep(10); 193 名前：名無しさん＠お腹いっぱい。 [2005/07/27(水) 01:39:55 ] ちょっと板違いなんですが、他に内容に沿ったスレ見つけられなかったんで失礼します。 今VisualC++で卒業研究のプログラム作ってるんですが、VC++で使える数値計算ライブラリってないですか？ 欲しいのは特異値分解だけなんですが、そこら辺に転がってるソース改変して使っても計算精度が悪くて使えないんです。 Octave C++もcygwinを通してのgccでしか使えないっぽいし、どうにかVC++で特異値分解をさせるライブラリのようなものお願いします。 194 名前：名無しさん＠お腹いっぱい。 mailto:sage [2005/07/27(水) 08:15:02 ] >>193 ttp://vision.kuee.kyoto-u.ac.jp/~nob/doc/liboctave/liboctave.html これじゃだめ? 195 名前：名無しさん＠お腹いっぱい。 mailto:sage [2005/07/27(水) 14:59:39 ] (.v.)ｵﾚ 196 名前：193 mailto:sage [2005/07/27(水) 17:12:02 ] >>194 若干板違いだったのに、わざわざ有用な情報をありがとうございます。 結構色々調べたつもりだったんですけど、Octave C++をVC++で使えるとは知らなかったです。 Octaveを使うのあきらめてclapack使おうかと思っててたんですが、こっちの方が良さそうですね。 では、卒業研究頑張ってきますﾉｼ 197 名前：193 mailto:sage [2005/07/28(木) 02:03:21 ] なんかうまくいかねーorz 198 名前：名無しさん＠お腹いっぱい。 [2005/10/08(土) 07:15:57 ] Octave で補間や補外ってできますか? 199 名前：名無しさん＠お腹いっぱい。 [2005/10/21(金) 23:36:59 ] octave-2.9.3 に UMFPACK をリンクさせたいんですけど、どーやっても configure が libumfpack.a を見つけてくれません。 誰か助けてー。 200 名前：名無しさん＠お腹いっぱい。 [2005/10/25(火) 23:34:33 ] octaveの可視化機能は、どこまでMatlabに近いですか？ www.cybernet.co.jp/matlab/support/helpdesk/r13/toolbox/matlab/creating_plots/chvolv19.shtml こんな画像を出したいのですが、octave-forgeを入れていても、関数が無いようです。 201 名前：199 [2005/10/26(水) 00:15:58 ] >199 自己レス。 UMFPACKv4.4 の共有ライブラリを自分で作って CVS 版の octave で configure したら認識しました。 でも make がコンパイルエラーで止まる…今度はデバッグか…orz 202 名前：199 [2005/10/27(木) 06:28:25 ] >201 さらに自己レス。解決しますた。静的ライブラリでも大丈夫だった。 足りないライブラリを持ってきて Makefile をやっつけハックでいけた。 2048 * 2048 三重対角行列の LU 分解 - Pentium III Coppermine 1GHz です。 octave:1> a=eye(2048).*rand(2048)+diag(ones(2047,1),1)+diag(ones(2047,1),1)'; octave:2> t=time();[L U]=lu(a);time()-t ans = 46.189 octave:1> b=sparse(eye(2048).*rand(2048)+diag(ones(2047,1),1)+diag(ones(2047,1),1)'); octave:2> t=time();[L U]=splu(b);time()-t ans = 0.30429 (*´д｀*) 203 名前：名無しさん＠お腹いっぱい。 [2005/12/01(木) 11:28:50 ] portsからoctaveをインストールしようと思ったんだけど、 atlasのコンパイルに無茶苦茶時間かかりますね。 CPU毎に最適化したatlasのpackageってどこかに落ちてない? 204 名前：203 mailto:sage [2005/12/01(木) 11:31:15 ] >>203 ここのスレ的にはいちいち言わなくても問題ないと思いますが、 使っているOSはFreeBSD6.0Rです。 ここの板にいると、ついついみんなFreeBSDを使っていると思ってしまう。 205 名前：名無しさん＠お腹いっぱい。 mailto:sage [2005/12/01(木) 22:07:16 ] FreeBSD 6.0Stable, pentium-m 1.3GHzでAtlasのチューニングに36時間かかった パッケージは保存したよ デュアルブートでBSD使ってるからもうあんなのやりたくねえ 206 名前：名無しさん＠お腹いっぱい。 mailto:sage [2005/12/02(金) 00:03:18 ] 今日はじめて使いました。 自分の作ったFFTのプログラムがあまりに遅かったので…。 ところで、起動時のフォルダって変更できるんですか？ 微妙に不便 207 名前：名無しさん＠お腹いっぱい。 [2005/12/03(土) 12:19:24 ] ？/.octaverc　に好きなコマンドを突っ込んでおけば起動時に実行される スタートアップファイルには他にもいろいろあって順番が決まっている 208 名前：名無しさん＠お腹いっぱい。 mailto:sage [2005/12/03(土) 20:42:54 ] >>207 なるほど。 ありがとうございました。 209 名前：名無しさん＠お腹いっぱい。 mailto:sage [2005/12/07(水) 14:24:02 ] しかし、Opteron dual な鯖で、昨夜の8時からatlasのインストールしているのに まだ終わらない。 210 名前：名無しさん＠お腹いっぱい。 [2005/12/26(月) 17:28:06 ] liboctaveでodeをする方法だれかわかりませんか？ 微分方程式がどけない、、、、、 211 名前：名無しさん＠お腹いっぱい。 mailto:sage [2005/12/28(水) 21:59:59 ] 非線形だからMatlabとかという理屈がよくわからん そんな非線形解なんてニュートン法でいいじゃん ニュートン法自体はすげー簡単で単純 学部の頃はQuickBASICで書いてたぞ 奥村晴彦のコンピュータアルゴリズム事典にも載ってるし ニュートン法でポテンシャルのパラメータを出して インパクトパラメータを計算してた 212 名前：名無しさん＠お腹いっぱい。 mailto:sage [2005/12/29(木) 18:19:06 ] 非線形っていうけど、Matlabって線形計算用なんだよね もともとは... 213 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/01/02(月) 07:04:26 ] FreeBSDでoctaveを使いたいのだけど、atlasのコンパイルに時間がかかりすぎます。 それでpackagesを探しているのですが、FreeBSD6.0 or 6-stable用の atlasのpackagesが見付かりません。 どこかにないでしょうか? 214 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/01/02(月) 10:22:57 ] ports にこう書いてあるから、自分で頑張る方が良いんじゃないかな... NO_PACKAGE= Optimizes for the local machine, so a package is not generally useful. Set FORCE_PACKAGE to override 215 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/01/02(月) 13:35:50 ] >>213 オレ、年末にoctaveのportsをmakeしたけど PentiumIII900MHzで30時間かかったよ 216 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/01/04(水) 10:13:40 ] LMItoolboxみたいのありますか？ 217 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/01/04(水) 22:50:27 ] ＞２１３ ATLASは実際にCPUにいろいろ計算させて、その結果を見て最適化していくんだよ。 だから、お前の希望するFreeBSD6.0 or 6-stable用、かつ、お前の使ってるCPU用のバイナリじゃなきゃだめじゃないの？ 218 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/01/04(水) 22:55:46 ] まあatlasが必要ない用途も多いな。 219 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/01/05(木) 02:28:44 ] >>217 Atlas使うといっても何時間もかかるような数値計算をするわけでも ないので、それほどカリカリに最適化してもらわなくてもいいんだけど... 220 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/01/05(木) 02:43:06 ] atlasってそんなにインストールに時間掛かるの？ openoffice位か？ 221 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/01/05(木) 06:06:19 ] 演算の最適化部分を省けばたいしたことないと思うけどね。 簡易インストールみたいなのないの? 222 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/01/05(木) 22:51:42 ] >221 最適化をあきらめる選択肢も必要だよね チューニングに時間がかかりすぎるし ATLASを使わないで、BLASだけでつかえたらいいのにね 223 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/01/05(木) 22:54:08 ] ところで、どうしてFreeBSDのPortsにあるATLASには-ffast-mathオプションが付いてるんだろう はずさないとだめなんじゃないの？ 224 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/01/05(木) 23:59:39 ] カリカリに最適化してくれるから俺はports ひいてはFreeBSDを使っているわけだが 225 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/01/06(金) 00:03:13 ] >222 そういう make option の盛り込み案を maho たんに提案しる 226 名前：名無しさん＠お腹いっぱい。 [2006/01/11(水) 17:41:06 BE:81633582- ] 質問なのですが、save のファイル名の引数に文字列を保持した変数を 与えることはできるのでしょうか？ a=[1,2; 3,4]; filename = "test.dat"; save filename a; とやると,"test.dat"ではなく"filename"というファイルに出力されてしまいます。 変数filenameの中身をみてくれる方法はないのでしょうか？ 227 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/01/13(金) 23:14:55 ] >226 ■俺も初心者。 とりあえず、save関数を使ってやってみてだめだった。つーか、マニュアルを見るとこうある。 たぶんだめなんじゃないのかな？代わりにfopenコマンドを見つけたから試してみたらどうかな？ １．fopenでファイルオープン　２．fputで書き込み　３．fcloseでファイルクローズ perlのファイルハンドルみたいだし、""無かったら変数も展開してくれそうな雰囲気はあるよな。 俺も日頃からクソつまらない作業ばかりの社会人で試す前に寝てしまった。暇ができたらやってみるけど。 ■saveコマンド 15.1.3 Simple File I/O save options file v1 v2 . . . Command Save the named variables v1, v2, . . . in the file file. The special filename '-' can beused to write the output to your terminal. ■15.2.1 Opening and Closing Files [fid, msg] = fopen (name, mode, arch) Built-in Function fid list = fopen ("all") Built-in Function file = fopen (fid) Built-in Function For example, 　　myfile = fopen ("splat.dat", "r", "ieee-le"); fclose (fid) Built-in Function ■15.2.2 Simple Output fputs (fid, string) Built-in Function Write a string to a file with no formatting. ■これを参照した。 GNU Octave　A high-level interactive language for numerical computations Edition 3 for Octave version 2.1.x 228 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/01/14(土) 00:44:12 ] どうしてもsaveにこだわるのなら eval(sprintf("save %s a",filename)) で可能。 こんな方法を使うメリットがあるかは疑問だが。 229 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/01/14(土) 01:29:15 ] >228 おおう なるほど 230 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/01/15(日) 00:02:40 ] MatlabではM関数をエデットすればそれが次ぎの実行で反映されるのですが Octaveでは一度読み込んだ関数はそのMファイルが書き変わっても反映されない様です それでしょうがなく関数を書き換えたらexitしてoctaveを立ち上げ直しているのですが 関数を初期化する方法は無いでしょうか？ 231 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/01/15(日) 08:34:26 ] clear all っぽい でも何かエラーが出る 232 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/01/15(日) 13:00:18 ] x="output.txt" y=[1,2,3,4] save("x","y"); ってなかったけ 233 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/01/15(日) 20:07:34 ] >>230 マニュアルには更新日時を見て 再読込するよな事が書いてあるけど、 実際の動作は違うのか？ help -i reload >>231 clear allは普通にできないか？ >>232 できない。 234 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/01/21(土) 18:48:32 ] Octaveで偏微分方程式は解けますか？ 235 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/02/07(火) 13:45:49 ] octaveとscilabってどう違いますか？ 初心者にお勧めはあります？ 236 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/03/05(日) 16:18:36 ] cygwinでscilabをコンパイルできた人居る？ cygwinのコンソールでscilabを使いたいよ 237 名前：名無しさん＠お腹いっぱい。 [2006/03/12(日) 21:25:19 ] Windows用 Octave 統合環境 Octave Workshop www.math.mcgill.ca/loisel/octave-workshop/ インストールしたらちゃんと動いた！ でもグラフは gnuplot for win エディタはメニューの File で open すれば立ち上がる これの凄い所は全部 MinGW で作ってある所！ Cygwin 無しなので、UNIXになじみの無い職場でも使えるかも 238 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/03/13(月) 16:40:29 ] そういえば、Octave Workshopって、FreeBSDのportsに入っていないね。 239 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/03/16(木) 14:32:49 ] floating ｐｏｉｎｔの形式はｉｅｅｅじゃなきゃだめなのね。 故に、-mfancy-math-387 -ffast-ｍａｔｈなんかつけるとこけてしまう。 240 名前：名無しさん＠お腹いっぱい。 mailto:age [2006/03/17(金) 23:32:02 ] FreeBSDのatlasのportsのfiles/{blas.h,atlas.h}って +#ifdef __cplusplus extern "C" { +#endif --略-- +#ifdef __cplusplus } +#endif なpatchをあてないとc言語の時に困るので 誰かsend-prしてください m(_ _)m 241 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/03/17(金) 23:52:31 ] cygwinでoctaveとatlasをインストールしようとしてたやついる？ gcc3.4でoctaveをコンパイルすると遅くなるって言ってたけど gcc4.0とか普通にコンパイルできるぞ？ 今からこれでoctaveをコンパイルしてみる 242 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/03/18(土) 06:12:02 ] >>241 ａｔｌａｓしらんが、ｏｃｔａｖｅは出来そうだな。 liboctave/lo-ieee.ｃｃかな？にCygwin_ｉｎｉｔとかいう関数がある。 243 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/03/20(月) 21:09:03 ] SciLabは4.0が出たみたいだな 使いにくいから使わないけど 244 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/03/21(火) 01:11:46 ] octaveの方がいろいろルーズでいいな 245 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/03/22(水) 12:30:18 ] >>244 確かにな mファイルをコマンドとして使えないとか 関数mファイルを明示的に読み込まないといけないとか、、etc.. Octaveみたいに操作方法もMatLab互換にすればみんな使うのに 246 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/04/03(月) 00:21:50 ] 俺はLMI tool(MATLAB非互換)が使いたいからSciLabも使うけどね。 慣れればSciLabもいいよ。 247 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/04/05(水) 23:06:34 ] octave-forgeのfillは頂点が4つまでしか描けないの？ 有限要素法の研究をoctaveでやらにゃならんのに 8頂点の要素がかけないのはつらい・・・。 情報をください。ｴﾛい人！！ 248 名前：名無しさん＠お腹いっぱい。 [2006/05/14(日) 07:20:32 ] >241 ATLASはgcc-3.3.3でmakeしました。3.4でもいいはずなのですが ファイルの一部がmakeできなかった？（よく覚えていない） octave以外はgcc-3.2でmake (3.4はbroken octaveになるそう。3.3.3は例外処理が遅いため遅い。 cygwin1.dllの一月のverupより遅さが緩和されているような気がしますが) gcc-3.2は普通にはgetできないので www.nuce.nagoya-u.ac.jp/e8/Matsuoka/OctCpp.html を参考にもってくる。 249 名前：名無しさん＠お腹いっぱい。 [2006/05/14(日) 07:23:26 ] >248 自己フォロー 　>octave以外はgcc-3.2でmake ごめん。 octaveはgcc-3.2でmake の間違い 250 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/05/20(土) 12:12:17 ] Octave Workshop (>>237) を使っています。 起動直後に私のホームフォルダに自動的に移動しておいて欲しいのですが、 どのファイルにコマンドを書き込めばいいのでしょうか。 あちらこちらにある設定ファイルらしきものに書き込んでみたのですが まったく変化がありませんでした。ヒントをもらえると嬉しいです。 251 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/05/20(土) 14:06:11 ] >250 自己解決しました。 普通にリンクファイルでフォルダを設定すればＯＫでした。 失礼しました。 252 名前：名無しさん＠お腹いっぱい。 [2006/06/25(日) 16:16:16 ] >>40のleastsqr.m を動かそうとしているんだけど、なかなかうまくいかない。 なにをやったか書くので、なにがいけないのか、分かる人がいたら教えてください。 x=[1 2 3 4 5]; y=[1 3 6 7 8]; function ret=F(x,p) ret = p(1)*x*x + p(2)*x + p(3); endfunction pin=[1 1 1]; [f,p]=leasqr(x,y,pin,"F") といった形で実行させると、 error: input(x)/output(y) data must have same number of rows といったエラーが出る。leastsqr.mのソースを見ると、 % check data vectors- same length? m=length(y); n=length(pin); p=pin;[m1,m2]=size(x); if m1〜=m ,error('input(x)/output(y) data must have same number of rows ') ,end; といった箇所があって、m=length(y)と[m1,m2]=size(x)のm1を比較している。 つまり、上の例だとmが5となり、m1が1となるので、両者がマッチしない。 だから、エラーになったということのようだった。 このあたり、データの与え方が間違っているのかと思って転置行列 x' で 実行させてみたりしたけれど、なかなかうまくいかない。Octave のバージョンは、 GNU Octave, version 2.1.53 (powerpc-apple-darwin7.5.0). GNU Octave, version 2.1.69 (i386-pc-linux-gnu). 両者とも同じ結果だった。 253 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/06/26(月) 10:54:43 ] >>252 そのメーリングリストで流れてたleasqr.mはバギーで動きません 元々のleasqr.mはMatLAB R12ごろに使われていたものですが フリーの非線形最小二乗関数のleasqr.mでは儲ける事が出来ないので 独自にlsqcurvefit.mを作りOptimization Toolboxに入れて売ってます leasqr.mの方は現在 Octave-forgeに収録されており、 最新版ならバグ無く動く様です(私の所では) 必要なのはleasqr.mとdfdp.mだけなので、octave-forgeをインストールしなくても このソースだけ拾ってくれば使えます。 このleasqr.mはmatlabでもそのまま使えますので optimization toolboxは買えないけどmatlabは使いたい人は使ってみてはどうでしょう lsqcurvefit.mと違って共分散行列(対角項が各説明変数の分散になる)と、 R^2つまり決定係数が得られて便利です。 254 名前：ゆーき [2006/06/27(火) 19:52:37 ] すみません。。。 matlabハイパー初心者なのですが。。。 csapeで曲線を描きたくて調べたのですが、どうしてもcsape(X,Y,CONDS,VALCONDS) のVALCONDSに何を入れたらいいのか分かりません。。。 特異値とか特異条件とかが関係してくるっぽいんですが。 どなたか教えて下さい。 255 名前：名無しさん＠お腹いっぱい。 [2006/06/27(火) 20:31:38 ] 最近，Octaveを使い始めたんですが，図を複数作る方法がわかりません． subplotみたいに１つのグラフにいくつかの波形をプロットするではなく， １つのグラフに１つの波形をプロットしたものをいくつか作る方法を教え てください． GNU Octave2.1.73を使っています．figureを使ってもグラフが表示されま せんでした． 256 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/06/27(火) 20:52:22 ] >>253 octave-forge から持ってきて、leasqrdemo.m をにらめっこしながら、 ようやく動かすことができました。 x=[1; 2; 3; 4; 5]; y=[1; 3; 6; 7; 8]; function ret=F(x,p) ret = p(1)*x.*x + p(2)*x + p(3); endfunction pin=[1 1 1]; [f,p]=leasqr(x,y,pin,"F") まず、xとyを列ベクトル標記に変え、さらにx*xの計算のところを x.*xに変えました。これで、動かし方を理解できました。 どうもありがとうございます。 257 名前：名無しさん＠お腹いっぱい。 [2006/07/03(月) 16:42:51 ] 初めて書き込ませていただきます。 構造体に関する質問なのですが、 例えば octave:1> Kozo(1).name = "ichiro" octave:2> Kozo(2).name = "jiro" octave:3> Kozo(1).height = 170 octave:4> Kozo(2).height = 175 というふうな構造体作ったのですが 変数aにichiroの身長を代入しようと a=Kozo(1).heightとしても scalar cannot be indexed with . とエラーが出てアクセスできません。 どなたか参照する方法を教えて頂けませんでしょうか？ octave:5> Kozo(1) Kozo(1) = { height = 170 name = ichiro } 258 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/07/04(火) 09:27:40 ] >>257 octave-2.1.71:1> Kozo(1).name = "ichiro"; octave-2.1.71:2> Kozo(2).name = "jiro"; octave-2.1.71:3> Kozo(1).height = 170; octave-2.1.71:4> Kozo(2).height = 175; octave-2.1.71:5> a=Kozo(1).height a = 170 octave-2.1.71:6> b=Kozo(1).name b = ichiro octave-2.1.71:7> fieldnames(Kozo) ans = { [1,1] = height [2,1] = name } 何の問題もありませんが？ 259 名前：名無しさん＠お腹いっぱい。 [2006/07/04(火) 13:43:14 ] >>258 別のマシーンでやってみたところ問題ありませんでした。 ということで文法上の問題でなく インストール上の問題のようです。ありがとうございました｡ 260 名前：名無しさん＠お腹いっぱい。 [2006/07/16(日) 02:39:09 ] forge入れて、変わったなーって実感したいんですけど 何かサンプルプログラムありませんか？ 261 名前：名無しさん＠お腹いっぱい。 [2006/07/16(日) 03:22:13 ] octaveで得たグラフとかって論文に載せてもいいの？ 信頼性りかの問題で 262 名前：名無しさん＠お腹いっぱい。 [2006/07/16(日) 03:38:54 ] り、のとこ、と、の間違いでした。 すいません。 263 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/07/16(日) 04:41:49 ] >>262 octaveってグラフはgnuplotで描いていると思うけど、 gnuplotは学術論文のグラフ描画によく使われているよ。 そういった意味では問題は少ないと思うな。 もしもoctaveの計算精度の問題を言っているのであれば、 それは自分で精度を検証してから使うべきだろう。 もちろん、octave自身けっこういろいろなところで使われているので 信頼性はそれほど低くはないと思うけど、 あまり使われていない関数などはバグがけっこう入っていそうだ。 264 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/07/16(日) 09:54:11 ] むしろ >261 は何だったら載せていいと思っているのか問い詰めたい 265 名前：名無しさん＠お腹いっぱい。 [2006/07/16(日) 14:56:30 ] FFTした後の振幅って、その中のピークで規格化するのが一般的なんですか？ 実際のパワースペクトル密度を得たいときはどうすればいいんですか？ 勉強不足ですいません。よろしくおねがいします。 266 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/07/16(日) 21:05:49 ] ニューラルネットの例題プログラム作ってみたけど 思ったより計算速度は速かった。 ベクトル演算が鍵なのかな。

---

---

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

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