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

1 名前：名無しさん＠お腹いっぱい。 [04/06/01 00:13] 知る人ぞ知る強力なデータ解析ツール octave ... unixを科学研究で使っている香具師には非常に心強い味方のはずだ。 matlab互換を主張しながら、微妙にバージョンアップを続ける 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 ] ニューラルネットの例題プログラム作ってみたけど 思ったより計算速度は速かった。 ベクトル演算が鍵なのかな。 267 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/07/16(日) 21:08:49 ] >>266 そうだろうね。 同じ演算でも行列使った時とループ使った時では格段に計算時間が違う。 268 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/07/22(土) 00:57:02 ] >>261 Excelよりは信頼性あるだろ 269 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/07/22(土) 23:03:41 ] そういや、あまり知られてないけどExcelってバグ多いよね。 一度、共同研究先が用意したデータを処理しようとしてExcel使ったけど、 マニュアルの記載にも間違いがあったりして、まともに 計算結果がでなくても困ったことがある。 そういった意味ではOctaveの方がましだと思った。 270 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/07/22(土) 23:16:57 ] 特に財務計算の関数はアレなんだっけ。 金の問題だけに検算しないとｗ 271 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/07/23(日) 00:12:51 ] 実験データを excel でいじろうとする奴を見ると 殴りたくなるよなっ! 272 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/07/23(日) 00:21:08 ] 現状のExcelじゃ扱える行数が少なくてどうにもならん。 次ので増えると聞いたが。 273 名前：名無しさん＠お腹いっぱい。 [2006/07/25(火) 01:33:11 ] ＞＞271 そうなの！？ 結構普通じゃない！？ 274 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/07/25(火) 11:20:44 ] >273 ガッ 275 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/07/27(木) 12:06:56 ] >>271 俺の周りはどんなデータもFORTRANでいじろうとする奴ばかりだな Octave使ってる俺だけ浮いてる orz 276 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/08/07(月) 17:07:19 ] すごくFAQっぽくて恐縮なのですが、質問させてください。 cygwinのコンソール上でoctaveを使おうとしています。 （正確にはoctaveのC++向けライブラリを使いたい） それで、cygwinのインストール時にg++関連（3.4.4）とoctave関連（2.1.73）をすべて導入し、 gnuplotのWindows版バイナリ（v4.0）にもパスを通して、 \usr\share\octave\2.1.73\m\startup\octavercと \usr\share\octave\site\m\startup\octaverc（違いがよく分からないので両方）に gnuplot_binary = "pgnuplot";を追加しました。 この状態でコンソールからoctaveと入力しても、「何もおきません」。 何か根本的に勘違いしているでしょうか？ 277 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/08/10(木) 01:33:29 ] octaveはscilabとかMatlabみたいに G=1/(1+s) って記述は出来ないんすか？ G=tf([0 1],[1 1]) ってわかりにくい。 278 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/08/10(木) 08:02:00 ] octave:1> s = 1 s = 1 octave:2> G = 1/(1+s) G = 0.50000 279 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/08/10(木) 10:19:02 ] function F(t) F=sin(t); endfunction t=0:0.1:2*pi; y=F(t) この最後のところで error: value on right hand side of assignment is undefined になるのは何故ですか？ 280 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/08/10(木) 10:20:17 ] >>277 そういう事がしたければSCILAB使えばいいじゃん。 281 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/08/10(木) 12:12:32 ] >>278 あとあと、ｔｆ関数とか使って、ボード線図が書きたいという趣旨ぢゃないすか？ 282 名前：名無しさん＠お腹いっぱい。 [2006/08/10(木) 12:13:40 ] >>266 ベクトル演算が速くなるように書くっていうのはCray1の昔からですねｗ 283 名前：278 mailto:sage [2006/08/10(木) 21:54:29 ] 趣旨とか言われてもシラネ 284 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/08/12(土) 17:35:48 ] >>283 G=tf([0 1],[1 1]) って書くと bode(G) って出来るけど G=1/(1+s) って書いても bode(G) って出来ないってことでしょ？ scilabなら s=poly(0,"s") G=1/(1+s) bode(G) って書ける。 285 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/08/14(月) 11:29:50 ] >>278 ﾜﾛｽ 286 名前：名無しさん＠お腹いっぱい。 [2006/08/18(金) 14:34:02 ] liboctave使っているのですが、liboctaveの説明十分なマニュアル・チュートリアル的な ものは存在するのですか？ 探してみましたが、liboctave wikiぐらいしか見つからなくて困っています・・・ 287 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/08/19(土) 16:23:09 ] マニュアルあってもいいけど Doxygenで作ったAPIリファレンスで事足りる 288 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/08/26(土) 19:57:32 ] FreeBSDの倍精度が64bitじゃなくて53bitで扱われてるって書いてある本があったけど・・・ どうしたらいいんだ 289 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/08/27(日) 05:50:44 ] 何が53bitなの? IEEE規格で計算しているから64bitのはず。 290 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/08/27(日) 10:58:54 ] >>288 仮数部の有効桁じゃないのか? 291 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/08/27(日) 13:18:12 ] 仮数部52ビット、符号部1ビット、指数部11ビット。 FreeBSDと限定する意味がわからんなあ、その本。 292 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/08/27(日) 14:23:30 ] >>291 「正規化方式で実質的に 53bit」って話は聞いたことがある 293 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/08/27(日) 16:57:17 ] intel のFPUは内部レジスタは80bitの精度を持っているけど、 IEEE規格だと64bitだからね。 もっともコンパイルオプションを指定すれば80bitで計算できたはずだけど... 294 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/08/27(日) 23:39:22 ] IEEEにも拡張精度は定義されてるけど 295 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/08/28(月) 02:05:31 ] 「IEEE」って… まあ意味分かるけど。 ttp://ja.wikipedia.org/wiki/%E6%B5%AE%E5%8B%95%E5%B0%8F%E6%95%B0%E7%82%B9%E6%95%B0 ttp://docs.sun.com/source/806-4847/ncg_goldberg.html 296 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/08/28(月) 07:08:00 ] man fpsetprec 297 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/08/28(月) 22:53:00 ] これか www.jp.freebsd.org/QandA/news/news-2003-10-16.html FreeBSD QandA 2003年10月16日 更新分 管理番号 587 (更新) 更新履歴 Q. gccで "long double" を使用して数値計算を行ないましたが、同じ CPU での Linux の結果と比べて計算精度が悪く、計算結果が異なります。 A. FreeBSD の標準カーネルでは浮動小数点の仮数部の精度が (IEEE754 と異なり) 53bit に設定されているために起きる現象です。 　 Intel x86 系の FPU の仮数部の精度は 24/53/64bit に設定可能です。 FreeBSD では仮数部の精度が 53bit に設定されているために、long double を使用しても仮数部の精度が 53bit (つまり、double と同じ) しかありません。 1. 関数 fpsetprec()/fpgetprec() を使用して、FPU の仮数部の精度を 64bit に変更する事が出来ます。fpsetprec(3) を御覧下さい。 2. GNU Multiple Precision Arithmetic Library (GMP)、PARI 等 FreeBSD 標 準以外の数値演算ライブラリを利用してもよいでしょう。 298 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/08/29(火) 22:29:53 ] カーネルの設定で64bitにするにはどうしたらいいの 299 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/08/30(水) 10:51:41 ] hackしろ 300 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/08/30(水) 11:22:55 ] ついでに見てやったぞ。 i386なら、/sys/i386/include/npx.hの__INITIAL_NPXCW__ amd64なら/sys/amd64/include/fpu.hの__INITAIL_FPUCW__ を適当に変えろ。 301 名前：２８８，２９７ mailto:sage [2006/08/30(水) 20:48:29 ] ソースを見て怪しい部分をさがせばいいのか デーモン君みたいだな よくわかんないけどがんばるよ 302 名前：名無しさん＠お腹いっぱい。 [2006/09/11(月) 14:40:45 ] octaveって(GNU Octave, version 2.1.73 (i686-pc-cygwin). A=[ 1 2 3 4 ] は A=[\ [1,2];\ [3,4]\ ] って、CRLFをエスケープしなきゃダメ？ 303 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/09/12(火) 14:43:02 ] やってみればわかる 304 名前：名無しさん＠お腹いっぱい。 [2006/09/12(火) 16:53:55 ] >>303 matlabだと A=[ 1 2 3 4 ] でいいのに、octaveだと\をつけなきゃならんのは、片手落ちだとおもう 305 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/09/13(水) 00:32:48 ] それ、 octaveじゃなくてターミナルの問題じゃね？ 306 名前：名無しさん＠お腹いっぱい。 [2006/09/13(水) 10:45:57 ] >>305 あ、ちょっと言葉足らずだった。ごめん A=[ 1 2 3 4 ] って、g.mに書いて、 source（"g.m") ってすると、ダメなんですよ。 って書いていて気がついた。 そか。改行をCRLFじゃなくてLFにしなきゃだめなんだ→やってみた→うまくいった。 307 名前：名無しさん＠お腹いっぱい。 [2006/10/30(月) 13:21:51 ] 卒研でoctave使って信号解析することになってるんだが。 音声ファイルの読み込みとかがいまいちよくわからん。 MATLABの本なんか見るとwavreadとかwavwriteとかsoundとかそういうのがあるっぽいんだけど octaveではそういうのできるのかい？ 308 名前：名無しさん＠お腹いっぱい。 [2006/10/30(月) 17:41:07 ] >>307 scilabならあるらしい？ 309 名前：名無しさん＠お腹いっぱい。 [2006/10/30(月) 17:42:16 ] >>307 nis-ei.eng.hokudai.ac.jp/~aoki/book01.html >(4) octave_programの中に，wavread.m，wavwrite.m，sound.m，test01.m，guitar_2048.wavをコピーする． 310 名前：307 [2006/10/31(火) 11:31:52 ] >>309 うぉすげえ。ｻﾝｸｽ。 っていうかこの本読んだんだけど wavplay.exeについては書いてなかった気がするんだよなぁ。 311 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/11/04(土) 11:58:22 ] CygwinでfftwとATLAS込みのOctaveを2.1.73をソースからビルドしたけど あんまり速度が速くならなくて悲しい 312 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/11/05(日) 16:53:22 ] Intel compiler + Intel MKLでOctaveビルドしてるいませんか？ 313 名前：名無しさん＠お腹いっぱい。 [2006/11/10(金) 02:01:59 ] SCILABでニュートン法を使ってtan(x)=C を解くプログラムを作りたいんですけど、 だれか、スクリプトを教えてくれません？ 314 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/11/10(金) 08:19:44 ] >>313 1)なんでわざわざニュートン法をプログラムするのか不明。解きたいならfsolveで。 2)その辺のCのプログラムでニュートン法のところ探してきて、まぁそれをそのまま書けばほぼOK 4)scilabのことはscilabのスレで 315 名前：名無しさん＠お腹いっぱい。 [2006/11/10(金) 21:03:30 ] >>313 宿題は自分でやりましょう 316 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/11/10(金) 22:11:47 ] (´-`).｡ｏO(　3)はどうしたんだろう…) 317 名前：名無しさん＠お腹いっぱい。 [2006/11/17(金) 01:25:59 ] かなり時間のかかるプログラムを走らせる際に，進行状況を見るためにループ文の中で if mod(n,50)==0 disp(n) ,end という感じで表示させたいのですが，サブ関数のループ中でやろうとすると，計算中には 表示をせずに，計算が終わったあとにまとめて全部表示する形になってしまいます 上手くサブ関数中にリアルタイムで表示させるような方法はありますでしょうか? 318 名前：317 mailto:sage [2006/11/17(金) 01:28:31 ] 最後の文のサブ関数中には計算中にの間違いです･･･ 319 名前：317 mailto:sage [2006/11/21(火) 23:49:38 ] ttp://dl.cybernet.co.jp/matlab/support/manual/r14/toolbox/optim/tutori34.shtml にあるように options = optimset('Display', 'iter'); を使うことで解決しました． 320 名前：名無しさん＠お腹いっぱい。 [2006/11/22(水) 01:43:34 ] workshop使ってみたんだけど、これエラーメッセージが文字化けして使えんね 俺だけ? ちなみにOSはXPsp2 321 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/11/22(水) 10:05:34 ] scilab 4.0でint16の挙動がWindows版とLinux版(どちらもバイナリ版をインス トールしました)で異なっているんですが、どちらが正しいんでしょうか。 Linux版 ------------------------------------------------------------------ -->int16([32767, 32768, 32769]) ans = ! 32767 -32768 -32768 ! ------------------------------------------------------------------ Windows版 ------------------------------------------------------------------ -->int16([32767, 32768, 32769]) ans = ! 32767 -32768 -32767 ! ------------------------------------------------------------------ 322 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/11/22(水) 12:44:56 ] 一応Scilab関連スレを張っとく SCILAB&SCICOSについて教えてください science4.2ch.net/test/read.cgi/sim/980480122/ Scilabを触ってみるｽﾚ science4.2ch.net/test/read.cgi/kikai/1068735555/ SCILABな science4.2ch.net/test/read.cgi/denki/1142648802/1-100 323 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/11/22(水) 13:52:21 ] >>322 誘導どうもです。そちらで聞いてみます。 324 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/11/24(金) 00:16:46 ] Octave（ライブラリ？）ってC言語（C++除く）から呼び出して使えますか？ 325 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/11/24(金) 06:29:38 ] liboctaveはC言語からは無理じゃないでしょうか 326 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/11/25(土) 00:55:44 ] そっか残念。 327 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/11/27(月) 08:44:46 ] ？ 328 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/12/12(火) 17:50:45 ] octave はガロア体上での計算できますか？ GF(2^16) くらいでいろいろやりたいんですが。 329 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/12/12(火) 18:43:49 ] とりあえず手元の Debian にパッケージ入れてみたら・・ gf() なんていう関数が見あたらんと・・ 330 名前：名無しさん＠お腹いっぱい。 mailto:sage [2006/12/18(月) 00:28:52 ] 質問です。 二重積分の計算をしようとして、以下のようなoctave スクリプトを書きました。 function y=a(x1) y = exp(-x1); endfunction function y=b(x2) y = quad("a",0,x2); endfunction quad("b",0,4) すると、以下のようなエラーが出ました。 error: quad: invalid recursive call error: evaluating assignment expression near line 8, column 5 error: called from `b' error: quad: evaluation of user-supplied function failed Octave のバージョンは 2.1.53 です。 なぜこのようなエラーとなるのか、教えていただけませんか。 331 名前：名無しさん＠お腹いっぱい。 [2006/12/30(土) 19:34:23 ] ループ中のグラフのアニメーションを、キー入力で途中で止める機能を付けたいので、ループ中に x=kbhit(); として、xの中身によりbreakするようにしたいのですが、上記の文を書くと、 error: value on right hand side of assignment is undefined error: evaluating assignment expression near line 4, column 2 error: near line 4 of file `kbhit.m' といったエラーが出ます。左辺のx=を書かずに、kbhit();だけだとエーラが出ません。 また、上記の文だけを書いたMスクリプトでも同様のエラーが出ます。 ちなみに、GNU Octave, version 2.9.9 (x86_64-redhat-linux-gnu)です。 アドバイスください！

---

---

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

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