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

1 名前：名無しさん＠お腹いっぱい。 [04/06/01 00:13] 知る人ぞ知る強力なデータ解析ツール octave ... unixを科学研究で使っている香具師には非常に心強い味方のはずだ。 matlab互換を主張しながら、微妙にバージョンアップを続ける octaveについて語ろう。 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)です。 アドバイスください！ 332 名前：名無しさん＠お腹いっぱい。 [2007/01/06(土) 15:27:37 ] アゲ☆アゲ☆ 333 名前：名無しさん＠お腹いっぱい。 mailto:sage [2007/01/06(土) 23:38:42 ] サゲ☆サゲ☆ 334 名前：名無しさん＠お腹いっぱい。 [2007/01/07(日) 20:28:55 ] アゲ♂アゲ♂everynight 335 名前：名無しさん＠お腹いっぱい。 mailto:sage [2007/01/07(日) 21:47:52 ] サゲ♂サゲ♂everynight 336 名前：名無しさん＠お腹いっぱい。 [2007/01/07(日) 22:08:07 ] 上げてしまう？ なかなかbuildがうまくいかない、Mac OSX版octave/octave-forge HPCで、2007/1月版が出てました。(gcc 4.3/gfortranも） hpc.sourceforge.net/ が、依然として、画像系が駄目。Octave.mのベンチマークすらも途中でこける。 Linuxでdual bootか、と転びそうな日々です。 337 名前：名無しさん＠お腹いっぱい。 mailto:sage [2007/01/09(火) 14:08:24 ] 答えを入れたファイルＸ個。これの平均をとってプロットしたいんだがどうすればよいのですか？ 338 名前：名無しさん＠お腹いっぱい。 mailto:sage [2007/01/12(金) 00:09:50 ] >>337 > xmean = mean(x)*ones(0,length(x)) ; > plot(xmean,'-') こんな感じかな。 339 名前：名無しさん＠お腹いっぱい。 mailto:sage [2007/01/13(土) 18:01:35 ] cygwin版のoctaveが起動せずに終了してしまうんだけど みんなもそうですか？ 340 名前：339 mailto:sage [2007/01/13(土) 18:16:40 ] パスに/usr/lib/lapakが入ってないだけだった すまん 341 名前：名無しさん＠お腹いっぱい。 mailto:sage [2007/01/27(土) 04:12:29 ] すんません、ちょっと遊んでると memory exhausted とか言われてしまうのですが、 どのくらいのサイズのリストまで作れるんでしょうか。 342 名前：名無しさん＠お腹いっぱい。 mailto:age [2007/04/12(木) 09:26:58 ] Octaveでシリアル通信する為の備忘録 (無保証です) octaveはcygwinに付属のもの、テスト用シリアルは内蔵モデム(COM3) octave:1> f=fopen('/dev/ttyS2','r+'); octave:2> fcntl(f,F_SETFL,O_NONBLOCK) ans = 0 octave:3> fprintf(f,'AT\r') octave:4> id1=fscanf(f,'%s') id1 = OK octave:5> id1=fscanf(f,'%s') id1 = octave:6> fclear(f) octave:7> fprintf(f,'ATI\r') octave:8> id1=fscanf(f,'%s') id1 = 56000OK octave:9> fclose(f); なんか動いてるようです。fcntlはoctaveはデフォルトではI/Oをブロックしてるので それを解除する為のおまじないです。baudrateなどの通信設定はsystemを使って windows の mode で設定するようです(modeについては cmd で help mode して下さい) 343 名前：名無しさん＠お腹いっぱい。 [2007/04/20(金) 09:23:28 ] FreeMATというのがあるみたいなんですが 使い心地はどうでしょうか？ ttp://journal.mycom.co.jp/news/2007/01/15/361.html 344 名前：名無しさん＠お腹いっぱい。 [2007/04/20(金) 10:22:51 ] 沖縄県の方へ（命に関わる注意事項です） 沖縄県での選挙ですが、どうか民主党だけは避けてください。県民の生命に関わる可能性があります。 民主党の最大の公約は一国二制度（※）ですが、一度「一国二制度　沖縄　三千万」で検索をお願いします。 この際、民主党のＨＰで調べても良いです。以下の注釈↓と矛盾することは書いてないはずですから… ※一国二制度 　簡単に言えば沖縄を中国と日本の共有物にし、そこに3000万人の中国人を入植させます。 　（つまり沖縄人口の 96% を中国人にして、実質、沖縄を中国人の居住地とします。） 　さらに「自主」の名の下、沖縄で有事が起きても自衛隊は干渉できません。 　3000万人の中国人が、少数派となった130万人の日本人に何をしても、です。 　そして反日教育を受けた中国人の反日感情の強さは、ほとんどの日本人の理解を超えるものです。 今回の選挙で民主党が勝った場合、「自主」「発展」を連呼しつつ段階的に進めていくことになります。 自主と言っても、自主を認めるのが「住人の96%が中国人となった」後だということに気をつけてください。 発展と言っても、新沖縄の少数派となった「少数民族日本人」の発展ではないことに気をつけてください。 345 名前：名無しさん＠お腹いっぱい。 [2007/10/15(月) 00:41:40 ] MAC OSX やってる人いますか？ scilabはMACでやってる人いるみたいだけど 346 名前：名無しさん＠お腹いっぱい。 [2007/10/15(月) 00:49:31 ] MAC OSX やってる人いますか？ scilabはMACでやってる人いるみたいだけど 347 名前：名無しさん＠お腹いっぱい。 [2007/11/03(土) 10:25:35 ] mac os 10.4 octave 2.9.9でやっているが ここ何日か他人に使わせるためのマニュアル作製の下調べやっているが、 どうにも素人に簡単にやらせるわけにはいかないようだ 出来る限りバイナリを貰って入れると言う方針で考えているのだが 2.9.15とそれに入っているgnuplotをドラックドロップ forgeのimageとioをしょうがないからpkg installで入れる でimreadが動かない。imagemagick6.3.6を入れる。 エラーを吐く。 identify の戻り値がver6.3以降？ではbits ではなくbitになっているので そこを直す。 さらにエラー freadが返すcountがあるべき値の2倍になっている。 糞だ。そびえ立つ糞だ つーかimread130行辺りのこの判定は何のためなんだ？ いやなんかしているのは分かるが、コメントアウトしていい？ というかこんな俺も分からんもんど素人に勧められねー なんでmacバイナリはforge入ってないんだ？ windowsのはあるのに。 というか俺が使っている2.9.9はまとめてバイナリだったのに どうして外すの？なんで 涙目 348 名前：名無しさん＠お腹いっぱい。 mailto:sage [2007/11/04(日) 15:25:57 ] 文体に心当たりが… 349 名前：347 mailto:sage [2007/11/08(木) 12:06:50 ] 心当たられるような事はしていないので人違いだろう。 で、結局あきらめた。 ○ーコしたMatlabでやってもらう事になるだろうな。 俺は努力した。 わざわざMatlab用のコードまで埋め込む事になるが どちらかといえばその方が早そうだ。 やれやれ 350 名前：名無しさん＠お腹いっぱい。 mailto:sage [2007/11/24(土) 13:29:12 ] >>347 Fink版でOctave2.9.14が出た（for Leopard?）。 自分は legend で困ってたんだけど、それも直った。 Forge も相当マージされた模様。 imreadも直ってるかも？ それとは別に質問なんですが、plot の横軸に日付って出せたっけ？ MATLABでは可能らしいというところまでは分かったんだけど・・。 351 名前：名無しさん＠お腹いっぱい。 mailto:sage [2007/11/25(日) 15:21:56 ] wiki.octave.org/wiki.pl?action=browse&diff=1&id=MatlabOctaveCompatibility ここでdatetickを探すとそれっぽい話してる？

---

---

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

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