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

1 名前：名無しさん＠お腹いっぱい。 [04/06/01 00:13] 知る人ぞ知る強力なデータ解析ツール octave ... unixを科学研究で使っている香具師には非常に心強い味方のはずだ。 matlab互換を主張しながら、微妙にバージョンアップを続ける octaveについて語ろう。 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 ] ニューラルネットの例題プログラム作ってみたけど 思ったより計算速度は速かった。 ベクトル演算が鍵なのかな。 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ならあるらしい？

---

---

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

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