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

1 名前：名無しさん＠お腹いっぱい。 [04/06/01 00:13] 知る人ぞ知る強力なデータ解析ツール octave ... unixを科学研究で使っている香具師には非常に心強い味方のはずだ。 matlab互換を主張しながら、微妙にバージョンアップを続ける octaveについて語ろう。 628 名前：625 [2008/09/04(木) 22:46:19 ] >>627 ご丁寧にありがとうございます。 聞く側としての初歩的な所を忘れていました。 現在、OctaveはMSVC版3.0.1を使用。 graphic backendはgnuplotで、バージョンは4.2になります。 ちなみにOSはXPのSP3です。 上記のバージョンのoctaveにて、 どうすればgnuplotで言う所の、 「plot 〜with impulse」 の様な表示が出来るのでしょうか。 以前使っていたバージョン(windows版octave2.173)では、 plot (〜,"^")と入力をしたら出来ていたのですが、 現在のバージョンで同様に入力をすると、 三角マークが列をなしてしまいました。 ご教授、よろしくお願いいたします。 629 名前：octave for windowsメモの筆者 [2008/09/05(金) 09:02:23 ] >>628　ウーム。どうやらバグですな。 help plotで表示すると `^' Set impulses plot style. ってでてくるし，Cygwin版でも同じ症状ですね。 いまから，本家に報告します。 630 名前：octave for windowsメモの筆者 [2008/09/05(金) 09:17:48 ] >>629 ttp://www-old.cae.wisc.edu/pipermail/bug-octave/2008-September/006673.html 本家bug@octave,orgにあげときました。お待ちください。 631 名前：octave for windowsメモの筆者 [2008/09/05(金) 10:31:06 ] >>630 どうもhelp,manualのミスでstemというのを使うらしい。 ttp://www-old.cae.wisc.edu/pipermail/bug-octave/2008-September/006675.html 詳しくは，help stemかマニュアルを見てください。 632 名前：名無しさん＠お腹いっぱい。 [2008/09/05(金) 21:55:26 ] >>631 先ほどstemで実行をしてみた所、 エラーこそ出ますが、なんとか表示されました。 ここの辺りは自分のプログラムと、helpを参照しながら 解消していこうと思います。ありがとうございました＾＾ 633 名前：名無しさん＠お腹いっぱい。 mailto:sage [2008/09/06(土) 18:56:48 ] 良かったね＾＾ 634 名前：octave for windowsメモの筆者 [2008/09/16(火) 12:57:12 ] gcc-4.3.2-1(dw2 例外処理)によるcygwin octave3.0.2-2が Marco Atzeriさんのサイトより手に入るそうです。 私もまだ試していません。 ttp://matzeri.altervista.org/ 635 名前：octave for windowsメモの筆者 [2008/09/16(火) 15:26:46 ] >>634 testしました。動作確認しました。 ttp://www-old.cae.wisc.edu/pipermail/octave-maintainers/2008-September/008486.html ということで私のgcc-3.4.4-dw2によるoctaveの配布は役割を終えました。 まもなくページは封鎖します。 少し寂しい気もしますが，Cygwin上のOctaveはこれですっきりするでしょう。 ただし，gnuplotは4.3(cvs)でないとマウスズームが使えません。 したがって，gnuplotは4.3(cvs)を配布することにしました。 ttp://www.tatsuromatsuoka.com/gnuplot/Eng/cygbin/ 636 名前：octave for windowsメモの筆者 [2008/09/17(水) 16:12:00 ] うーん。cygwin版は3.0.2版は出たのですが， msvc版は3.0.2はMichaelさんはスキップするのかな？ BenjaminさんのMINGW版は，いつになったら出るんだろう。 彼は慎重派だからな。 で, Windows Native版3.0.2は，どうなるんだろう と思っています。 それでmingwの3.0.2のあるんですが，BenjaminさんのMINGW版がでるんならと 出していませんが，興味のあるかたいますか。一人でもいたら，私のページ にあげますが。 637 名前：octave for windowsメモの筆者 [2008/09/18(木) 12:10:15 ] >>634 Octave 3.0.2-2 on cygwinがcygwin setupからでもinstall できるようになりました。 ttp://cygwin.com/ml/cygwin/2008-09/msg00398.html でもforge20080831-2はインストールできないので手動となります。 638 名前：octave for windowsメモの筆者 [2008/09/18(木) 14:37:00 ] >>637 なんかHangしますね。しばらくは手動でインストールしてください。 一応報告したので解決すると思いますので手動インストールができない方は 対応するまでお待ちください。 639 名前：octave for windowsメモの筆者 [2008/09/19(金) 17:53:15 ] >>636 Michaelさんに確認したところ，MSVC版のoctaveは3.0.2はスキップ するそうです。 あまり時間をおかずに Jaroslavさんが3.0.3を出しそうないきおいですからね。 640 名前：octave for windowsメモの筆者 [2008/09/22(月) 08:47:48 ] >>638 一部あちらのミスもありましたが，私の思い違いでした。 experimental versionha expのボタンをアクティブにしてインストール するそうです。 ttp://cygwin.com/ml/cygwin-apps/2008-09/msg00081.html 641 名前：名無しさん＠お腹いっぱい。 [2008/09/24(水) 17:49:01 ] >>640 expボタンを選ばなくても3.0.2-2を選び binところののチェックボタン （□（四角））をチェックするのを忘れないようにすればいいようです。 642 名前：octave for windowsメモの筆者 [2008/09/29(月) 09:22:42 ] maintainers向けに3.0.3のプレリリースがありました。 各maintainerも確認したところです。 ただし，マイナーチェンジなため各バイナリーが出るかは分かりません。 MSVC版は3.0.2をスキップしたため3.0.3の準備をミカエルさんがされている ようです。 643 名前：名無しさん＠お腹いっぱい。 [2008/09/29(月) 11:24:48 ] 初めて立ち寄りました。 場違いですが、matlabにオイラー法を使えるコマンドはないのでしょうか？ 644 名前：octave for windowsメモの筆者 [2008/09/30(火) 18:35:06 ] >>643 なぜまたオイラー法が必要なんでしょう。 MATLABには ODE45などのルンゲクッタ法の発展型の優れた関数があるのに。 以下はOctaveでしか動きませんが， ODEの解法の原理として講義で使ったものです。 関数を引数として受けるやり方は良く分かりませんが。 % Octave script m file % Eulerのテスト 1; %dummy % ****** Definition of a function function dy=f(y,t); dy=-y; end % ******* Main script t0=0; tend=5; N=100; %t0 初期値，tend 範囲の終わり, N　間隔の数（点の数は，N+1） t=linspace(t0,tend,N+1); % tの配列を定義 dt=(tend-t0)/N; % tの刻み幅 y(1)=1; % yの初期値 % for n=1:N y(n+1)=y(n) + f(y(n),t(n))*dt; end % plot(t,exp(-t), "", t, y, "@"); 645 名前：名無しさん＠お腹いっぱい。 mailto:sage [2008/10/02(木) 12:10:19 ] pc11.2ch.net/test/read.cgi/tech/1221144557/123 646 名前：名無しさん＠お腹いっぱい。 [2008/10/02(木) 12:46:32 ] あれなんかリンク貼ってある。 >>645じゃないですけど、リンク先の質問したものです。 インテルコンパイラでのコンパイルはどのように行うのでしょうか？ 647 名前：octave for windowsメモの筆者 [2008/10/02(木) 13:44:52 ] 自分でOctaveをintel complierでビルドする必要があります。 そういう方はLinux userには多いですね。 本家のMLで検索して調べてください。（Nabbleサーチが便利です。） 日本でもやっている方はいるとおもうんですが，ここを見ておられるかと いうとそういう猛者の方は英語のMLを参考にされていると考えたほうが いいとおもいます。 がんばってください。 648 名前：646 [2008/10/02(木) 14:01:01 ] >>647 "intel compiler"で検索してみたら結構な数がでてきたのでちょっと読んでみます。 アドバイスありがとうございます。 649 名前：octave for windowsメモの筆者 [2008/10/02(木) 16:40:37 ] Mingw版の正規版のOctave-3.0.2がリリースされました。 ttp://sourceforge.net/project/showfiles.php?group_id=2888 からダウンロードできます。 gnuplotは4.3cvsリリースをpguplotとして使っているようなので Michaelさんのものと違いマウスズームが効くとおもいます。 週末試してから私のWebには，リリースニュースとしてあげようとおもいます。 650 名前：octave for windowsメモの筆者 [2008/10/06(月) 08:46:06 ] Octave 3.0.3 のソースが公開されました artax.karlin.mff.cuni.cz/~hajej2am/ulozna/octave/ Octave の公式ページには，そのうち公開されるでしょう。基本的にはバグフィックスが中心です。Michael氏はおそらく この3.0.3でMSVC版のUpdateをするのではと踏んでいます。 651 名前：名無しさん＠お腹いっぱい。 [2008/10/07(火) 00:28:21 ] （　ﾞﾟ'ωﾟ'）ﾓﾏﾏﾏ…ﾏｼﾞで！？ love6.2ch.net/test/read.cgi/kao/1198685446/ 652 名前：octave for windowsメモの筆者 [2008/10/10(金) 18:15:52 ] >>650 読みどおりMSVC Ver.3.0.3はまもなく出るようです。 sourceforge.net/mailarchive/forum.php?thread_name=128f38bd0810100131q4b66f14aq1af1db4e03d79c06%40mail.gmail.com&forum_name=octave-dev 653 名前：名無しさん＠お腹いっぱい。 [2008/10/18(土) 00:33:36 ] スレの皆様に質問があります。 MSVC版を愛用しているものです。 私の目的では32bitだとメモリが少なすぎるように感じています。 64bit windowsではMSVC走りますか？ それよりもunixに移行するのが手っ取り早いのでしょうか。 64bitを使う上で気をつけたほうがいいことなどありましたら 教えてください。 (Matlabを買うというのが一番だとは思うのですが) 654 名前：octave for windowsメモの筆者 [2008/10/20(月) 08:14:00 ] octdevスレッドに同様の質問があり，MicahelさんとBenjaminさんが答えてます。 ただし，このスレッドは他のTopicと混じってますので， Status of Windows 64-bit support?　のタイトルの記事を見てください。 自分でコンパイルしろ書いてあります。 LInuxでも64bitでは自分でビルドする必要があると思います。 ただし，自分でビルドする手間は，windowsよりLinuxの方が簡単かとおもいます。 sourceforge.net/mailarchive/message.php?msg_id=128f38bd0810100418i2caf46f5rbb118d8b9031954c%40mail.gmail.com 655 名前：船木康博 [2008/10/20(月) 08:22:56 ] 信濃の国でも伊那方面は、有線でxDSLを最初にやったりさぁ、 無線じゃあ一昔前に妖怪がネット組んでいるのさ。おもしれぇぞ。 656 名前：名無しさん＠お腹いっぱい。 [2008/11/02(日) 19:24:11 ] >>654 ありがとうございます。急に立て込んでしまって、まだ試せていません。 とりあえずお礼まで。 657 名前：octave初心者 [2008/11/08(土) 09:50:11 ] cygwinからoctaveをインストールして、imageパッケージを入れて 画像処理を行いたいのですが、imreadするところで下記エラーが出力されて 画像を読み込むことができません。 （ImageMagickもcygwinからインストールしました） im=imread('im1.bmp'); convert: unable to open image `/cygwin/tmp/oct-22d688.0.ppm': No such file or directory. error: imread: error running ImageMagick's 'convert' error: evaluating if command near line 113, column 5 error: called from `imread' in file `/home/tosiyuki/octave/image-1.0.8/imread.m' error: evaluating assignment expression near line 1, column 5 どなたか解決方法を教えてもらえませんか？ 658 名前：octave for windowsメモの筆者 [2008/11/13(木) 07:59:52 ] Macroさんに教えてもらいました。 Hi Tatsuro, it is not a fault of image. I don't know for wich reason on cygwinthe TMPDIR enviroment is set to C:/CYGWIN/TMP instead of /TMP see src/sysdep.cc call to cygwin32_conv_to_full_win32_path Probably is for old compatibility with non cygwin programs, but it confuse cygwin programs likeconvert of ImageMagick. As workound is possible to setenv ("TMPDIR","/tmp") before im=imread('im1.bmp'); I will try to remove the issue on the next octave cygwin release. 659 名前：659 [2008/11/22(土) 15:24:37 ] はじめまして。 windowsVistaで、スタンドアローンoctaveを使おうとしています。 使いたいのは画像処理のpackage内のedge.mです。 image packageをインストールまでしたのですが、 edge("***.jpg", 'canny', [0 7], 1.6);などしても error: nonmax_supress underined...と言われ、エラー箇所を見ると bw = nonmax_supress(Es, Eo, thresh(1), thresh(2));とありました。 nonmax_supress()は関数のようですが、これが見あたらないようです。 どうすればいいのでしょうか？ (cannyを使わないでprewittを使うと動作が確認できました。BW = edge(I,'prewitt')) 660 名前：659 [2008/11/22(土) 15:25:51 ] UNIX板なのですが、octaveに関する場所がここしかないので、 ここで聞かせていただきました。場所違いでしたらすみません。 661 名前：659 [2008/11/22(土) 15:33:46 ] nonmax_supress()はCanny専用の関数のようです。 なぜ読み込まれていない（または組み込まれていない）のでしょう？ octave.sourceforge.net/doc/f/nonmax_supress.html 662 名前：659 [2008/11/22(土) 15:48:16 ] octave.svn.sourceforge.net/viewvc/octave/trunk/octave-forge/main/image/src/nonmax_supress.cc?revision=HEAD&content-type=text/plain を見るとnonmax_supress()はC++で書かれていて、 それをoctaveで動的に利用できるようです（www.obihiro.ac.jp/~suzukim/masuda/octave/html/octave_71.html）。 C++のソースコードは公開されているので、 それをコピーして手動でoctaveでコンパイルして、 nonmax_supress()を持っておけばいいのだと思うのですが。。。 663 名前：octave for windowsメモの筆者 [2008/11/25(火) 09:34:43 ] >>659 image packageは自分で使わないのでよく分かりません つい最近，octave-3.0.3 (MSVC)がでたのでとりあえずそれを使って 試してください。だめならhelp@octave.orgに投稿します。 sourceforge.net/project/showfiles.php?group_id=2888&package_id=40078 664 名前：名無しさん＠お腹いっぱい。 mailto:sage [2008/11/25(火) 13:47:04 ] >>309 Octave3.01では wavread.m, wavwrite.m は不要 sound.m の中身を一行 system("wavplay \"%s\"", wavefile); とダブルクォーテーションに変えてやる必要がある。 検索してたどり着く人用コメントでした 665 名前：octave for windowsメモの筆者 [2008/11/25(火) 18:47:19 ] >>659 文法違いませんか。 以下のように使うのでは，edgeの第一引数はMxNの グレースケールの行列のはずでファイル名ではないですよ。 help imread，help edge　でhelpを読みましょう。 octave:2> IM=imread('test.jpg'); octave:3> image(IM); octave:4> BM=edge(IM(:,:,1),'Canny'); octave:5> image(BM); 666 名前：octave for windowsメモの筆者 [2008/11/25(火) 19:04:32 ] >>665 すみませんm(__)m。書き方が刺々しいですね。 反省m(__)m。 667 名前：名無しさん＠お腹いっぱい。 mailto:sage [2008/12/06(土) 01:09:30 ] (;'ｰ`)FreeBSDのportsにあるoctaveのコンパイルオプションに-ffast-mathつけるのやめてくれ 668 名前：名無しさん＠お腹いっぱい。 mailto:sage [2008/12/06(土) 10:56:47 ] 理由があるならそれ書いて send-pr しれ 669 名前：octave for windowsメモの筆者 [2008/12/07(日) 14:34:16 ] >>666 確かに変ですね。いつか本家で聞いてみます。 670 名前：octave for windowsメモの筆者 [2008/12/07(日) 14:37:22 ] さっきのは>>667、>>668の間違い。m(__)m 671 名前：octave for windowsメモの筆者 [2008/12/08(月) 13:10:45 ] >>667 >>668 いま，octaveのMLできいていますが，FreeBSDで聞けといわれたらすみませんが 自分でFreeBSDの方で聞いてください。m(__)m 672 名前：名無しさん＠お腹いっぱい。 mailto:sage [2008/12/08(月) 15:27:43 ] あ，668 の人だけど >668 は >667 宛てね (FreeBSD 事情なら >668 ということで) 673 名前：octave for windowsメモの筆者 [2008/12/08(月) 18:35:09 ] Octave　MLより 3.0.4 call for patches (Jaroslav Hajek-2 Dec 08, 2008; 03:49pm) www.nabble.com/3.0.4-call-for-patches-to20890571.html ということでOctave-3.0.4のリリースの準備が始まりました。 674 名前：octave for windowsメモの筆者 [2008/12/09(火) 09:51:01 ] >>667 最初の書き込みがいけなかったため，該当する返事はまだ得られていません。 （-ffast-mathオプションが使えるかという風にとられたため。 違う趣旨の返事しか来てません。再度趣旨を訂正しましたがこれまでのところ 答えはついていません） では， www.nabble.com/-ffast-math-option-at-compling-octave-in-FreeBSD-ports---to20878991.html 675 名前：octave for windowsメモの筆者 [2008/12/10(水) 13:28:53 ] >>668 --- John W. Eaton wrote: > I think you should ask the people who packaged Octave for FreeBSD. > > Based on the description in the GCC manual, I would generally > recommend against using -ffast-math. > > jwe ということで，FreeBSDで聞いてくださいとのこと。m(__)m jwe氏もagainstと書いている通り，おすすめしないといっておられます。 -ffast-mathなんでついてるのかな。 676 名前：名無しさん＠お腹いっぱい。 mailto:sage [2008/12/12(金) 03:34:33 ] (;'ｰ`)つーか、-ffast-math付けると計算をはしょるから意図しない答えになることがある 677 名前：octave for windowsメモの筆者 [2008/12/12(金) 09:58:45 ] >>676 なるほどそうなんですか。私はJaroslav Hajekさんが書かれたことぐらいしか 知らないですが。どちらにしろつけるのは，変ですよね。 678 名前：名無しさん＠お腹いっぱい。 mailto:sage [2008/12/12(金) 10:50:43 ] portsのMakefileには-ffast-mathついてないんだけど、 /etc/make.confに自分で書いてないかな？ 679 名前：octave for windowsメモの筆者 [2008/12/12(金) 14:57:09 ] まあ，私が書き込んだことでoctaveのhelpとmaintainersのMLに残るのでFreeBSD portsを使っている方へのアナウンスにはなっているわけでだし あとはFreeBSD のportsをメンテされてる方がどう考えているかですよね。 ところでmkoctfile -vで出てくるg++でのコンパイルの表示に -ffast-mathついてますか? 基本的にはconfigureで指定したoptionがついてますけどね。 680 名前：名無しさん＠お腹いっぱい。 mailto:sage [2008/12/12(金) 20:02:11 ] (;'ｰ`)うわっ、ごめんよ！ FreeBSDのportsに-ffast-mathついてるのはmath/atlasだったよ 681 名前：octave for windowsの筆者 [2008/12/14(日) 12:32:14 ] >>680 了解。上にあげときます。 682 名前：octave for windowsメモの筆者 [2008/12/15(月) 11:52:18 ] >>680 上に報告しときました。 683 名前：名無しさん＠お腹いっぱい。 mailto:sage [2009/01/08(木) 16:48:40 ] 誰かovtaveで区分求積法を使って積分の計算をやってる例とか貼ってくれない？ 684 名前：octave for windowsメモの筆者 [2009/01/08(木) 19:29:46 ] >>683 区分求積って台形公式でいいんでしょうか。 I = trapz (Y) ; I = trapz (X,Y)； 台形公式法（trapezodial method）による，数値積分をする関数。X,Yはともにベクトルであり，引数がYだけ時はX軸の間隔が1であるとしてyの台形の面積の和をとる。通常は，X,Yの両変数を指定し，和をとる。 685 名前：octave for windowsメモの筆者 [2009/01/08(木) 19:30:25 ] octave> format long octave> x=linspace(0,pi,100+1); octave> I=trapz(x,sin(x)) I = 1.99983550388744 octave> x=linspace(0,pi,200+1); octave> I=trapz(x,sin(x)) I = 1.99995887647922 octave> x=linspace(0,pi,1000+1); octave> I=trapz(x,sin(x)) I = 1.99999835506566 686 名前：octave for windowsメモの筆者 [2009/01/08(木) 19:32:58 ] 精度が必要な時は，適応的自動積分のquadを使ってください。 687 名前：octave for windowsメモの筆者 [2009/01/08(木) 19:36:22 ] [v, ier, nfun, err] = quad (f, a, b, tol, sing) Quadpackという積分パッケージに収録された適応的自動積分の関数。 第1引数は，被積分関数の名前で関数は， y = f (x) の形を持ち，xとyは、スカラーである。第２，第３引数は，積分範囲であり，−無限大や+無限大を片側または両方に指定できる。Octaveで無限大は，infまたはInfであらわす。 第４引数のtol は，ベクトルであり，第一要素は，要求したい精度の絶対値である。第二要素は，相対許容度である。 第５引数のsingは，ベクトルで被積分関数が特異点となる独立変数の値からなり，必要に応じて指定する。 第４，第５引数は，オプション引数なので省略は可能である。 左辺にvは，積分結果の値（関数の返り値）。ierはエラーコードで0は，積分の成功を意味する。nfunは関数の評価回数でと積分の困難さの指標になる。Errは，積分結果の評価誤差である。 688 名前：octave for windowsメモの筆者 [2009/01/08(木) 19:37:13 ] octave> format long octave> function y=f(x);y=1./sqrt(1-x.^2);end; octave> [v, ier, nfun, err] = quad (@f,0,1) v = 1.57079632679489 ier = 0 nfun = 315 err = 5.44135847491134e-011 octave> pi/2 ans = 1.57079632679490 689 名前：名無しさん＠お腹いっぱい。 mailto:sage [2009/01/09(金) 03:56:42 ] 『octave for windowsメモの筆者』さんすごい！素敵です！ 690 名前：名無しさん＠お腹いっぱい。 mailto:sage [2009/01/09(金) 16:31:09 ] >>683ですが、貼ってくれた方ありがとうございます！ 691 名前：octave for windowsメモの筆者 [2009/01/09(金) 16:57:46 ] >>689 >>690 大学の講義で使っている資料の一部を貼りました。 著作権の問題があり公開はしていませんが， 著作権問題がない部分については公開したほうがいいですかね。 ちょっと手間がかかるので要望あればぼちぼちやります。 692 名前：名無しさん＠お腹いっぱい。 [2009/01/16(金) 16:35:33 ] 大学の宿題で「観測値と予測値をプロットするoctaveのプログラムを示しなさい。」 っていう問題が出たんですけど授業を聞いてないボクにはさっぱりです・・・。。 誰か教えてくれませんか？ 693 名前：名無しさん＠お腹いっぱい。 mailto:sage [2009/01/16(金) 16:52:59 ] ノート買えばいいじゃん 694 名前：octave for windowsメモの筆者 [2009/01/16(金) 18:26:40 ] Octave 3.0.4のテストリリースのソースのリリースがありました。 しばらく，すれば3.0.4が出ることになります。 ttp://www.nabble.com/3.0.4-RC1-to21495212.html 695 名前：octave for windowsメモの筆者 [2009/01/16(金) 18:39:20 ] >>692 まあ，>>693さんがいっているように授業受けた人に聞きましょう。 ↓に私の講義資料で著作権に関係のない部分をのせといたからこれでも読んで勉強しなさい。 ttp://www.tatsuromatsuoka.com/octave/jpn/08CompAlgo_02.pdf 696 名前：名無しさん＠お腹いっぱい。 mailto:sage [2009/01/16(金) 22:36:07 ] 予測値ってのはもう用意されてるのかね 予測手法も適切に選択するまでが課題だったら随分大雑把で困難な宿題だ 697 名前：名無しさん＠お腹いっぱい。 [2009/01/17(土) 03:41:12 ] すいません。 勉強します・・・。 698 名前：名無しさん＠お腹いっぱい。 [2009/01/26(月) 02:16:03 ] 質問させてください。 cos(pi*90/180) を実行すると 6.1232e-17 というネイピア数で表示されてしまうのですが、これを0と表示させるにはどうすれば良いのでしょうか？ 699 名前：名無しさん＠お腹いっぱい。 mailto:sage [2009/01/26(月) 08:55:39 ] つ round 700 名前：名無しさん＠お腹いっぱい。 mailto:sage [2009/01/26(月) 10:33:35 ] ところでネイピア数は0.271828...e-1だが？ 701 名前：名無しさん＠お腹いっぱい。 mailto:sage [2009/01/26(月) 10:38:31 ] じゃない、0.271828...e1 だ というか質問の意味がいまわかった。 そのeはネイピア数じゃなくて、浮動小数点数の10のべき乗表記で、 6.1232 × 10 ^ -17 じゃまいか？ 702 名前：名無しさん＠お腹いっぱい。 [2009/01/26(月) 13:03:46 ] >>699 roundで丸めるのではなく、不動小数点の10のベキ乗表記を小数点第何位までにするかという設定（要するに0の極限値の設定方法）があった気がしたので、そのやり方がわからないのです。 >>701 ネイピア数の定義がよくわかってなくてすみませんでした。 おっしゃるとおり6.1232 × 10 ^ -17です。 703 名前：octave for windowsメモの筆者 [2009/01/26(月) 19:31:40 ] octave3.0.4のソースのRC2が出ました。 ttp://www-old.cae.wisc.edu/pipermail/octave-maintainers/2009-January/010310.html octave3.0.4ももうじきです。 一応3.0シリーズのラストと位置づけられているようで， octave 3.2シリーズの準備がoctave 3.1で行われています。 704 名前：octave for windowsメモの筆者 [2009/01/27(火) 08:36:09 ] >>703 しかし，3.0になるときも，次は3.0とでながらかなり細かいバージョンアップを続けた経緯があるのでどうなるかは。分かりませんが。 705 名前：octave for windowsメモの筆者 [2009/01/28(水) 07:54:39 ] cygwin版のoctave-3.0.3がでました。ただし，パッケージングミスがあり, mkoctfileが動きません。/usr/lib/octaveのlibctave.dll.a(oが抜けている) をliboctave.dll.aにすれば動きました。報告をして，Macroさんが確認された ので数日中にはfixされるでしょう。 706 名前：名無しさん＠お腹いっぱい。 [2009/02/01(日) 00:21:14 ] output_precision=16って入力して表示桁数を変えようと思ったら, output_precision is no longer a built-in variable と出ました。 その後output_precition=16となりましたが数値の表示桁数が変わりませんでした。 どうすればよいのでしょうか 707 名前：octave for windowsの筆者 [2009/02/01(日) 18:03:15 ] >>706 output_precision(16)としてください。octaveはver3から殆この手のものは=ではなくて、()で設定します。 708 名前：octave for windowsの筆者 [2009/02/01(日) 18:05:12 ] >>706 output_precision(16)としてください。octaveはver3から殆この手のものは=ではなくて、()で設定します。 709 名前：octave for windowsの筆者 [2009/02/01(日) 18:05:36 ] >>706 output_precision(16)としてください。octaveはver3から殆この手のものは=ではなくて、()で設定します。 710 名前：octave for windowsの筆者 [2009/02/01(日) 18:06:06 ] >>706 output_precision(16)としてください。octaveはver3から殆この手のものは=ではなくて、()で設定します。 711 名前：octave for windowsの筆者 [2009/02/01(日) 18:08:53 ] 携帯から書込みしたら重複投稿してしまいました。ごめんなさい。m(_ _)m 712 名前：名無しさん＠お腹いっぱい。 mailto:sage [2009/02/02(月) 17:07:41 ] octaveでシミュレーションをやっているんですがどうしても下のようなエラーが出てしまいます 先生や先輩に聞いてもなぜかよくわからないというしまつ・・・orz どなたか教えていただけないでしょうか？よろしくお願いします octave:2> test2 error: invalid column index = 2001 error: evaluating binary operator `+' near line 72, column 28 error: evaluating assignment expression near line 72, column 17 error: evaluating for command near line 67, column 1 error: near line 73 of file `/home/soturon/test2.m' ソース→www.uploda.org/uporg1984631.tex.html 713 名前：octave for windowsメモの筆者 [2009/02/03(火) 10:48:16 ] 712 ソースの場所に行ってみましたが見れませんでした^^; 714 名前：名無しさん＠お腹いっぱい。 mailto:sage [2009/02/03(火) 10:54:34 ] めでたしめでたし 715 名前：712 mailto:sage [2009/02/03(火) 11:36:10 ] すいませんでした！ 次こそ大丈夫なはずです ソース→ultraup.net/man/src/u_man2312.zip.html パス：test 716 名前：名無しさん＠お腹いっぱい。 [2009/02/03(火) 12:25:39 ] >>715 Get Files Thank youとでますが，どこに落ちたかわからん。 サイトがサイトだけに，落ちないのかも，octave for windowsメモで Googleれば私のサイトは見つかるのでそこのアドレスに添付で 送ってください。 717 名前：名無しさん＠お腹いっぱい。 [2009/02/03(火) 12:26:31 ] >>715 Get Files Thank youとでますが，どこに落ちたかわからん。 サイトがサイトだけに，落ちないのかも，octave for windowsメモで Googleれば私のサイトは見つかるのでそこのメアドに添付で 送ってください。 718 名前：octave for windowsメモの筆者 [2009/02/03(火) 12:28:23 ] >>716 m(__)m そこのアドレス==>そこに記載のメールアドレス　 719 名前：名無しさん＠お腹いっぱい。 mailto:sage [2009/02/03(火) 13:32:15 ] メール欄にsageと入れよう！ 720 名前：octave for windowsメモの筆者 mailto:sage [2009/02/03(火) 13:40:12 ] >>719 入れましたがどのような意味で？　おじさんはワカラナイ？？ 721 名前：712 mailto:sage [2009/02/03(火) 13:50:48 ] たびたび申し訳ないです＞＜ メールで添付しときました！よろしくです つーかアップローダーはどこがいいんだろ・・・ うｐろだとか有名なとこはすぐ消えちゃうし 722 名前：octave for windowsメモの筆者 mailto:sage [2009/02/03(火) 14:06:50 ] >>721 メールとどきましたこれから見ます。 Uploaderではなく，Yahooのアカウントをお持ちなら，YahooのGeocitiesをフリーで 広告付きだけど50Mbyteのweb spaceを借りれますよ。 私は，自分のドメインと1Gbyteのスペースと持ちたいから 有料の契約してますが，普段使う分には50Mbyteあればいいでしょう。 723 名前：名無しさん＠お腹いっぱい。 mailto:sage [2009/02/03(火) 14:13:21 ] >>722 Geocities？ブリーフケースのこと？？ 724 名前：名無しさん＠お腹いっぱい。 mailto:sage [2009/02/03(火) 14:15:17 ] >>720 メール欄にsageと入れると書き込みをしてもスレッド一覧の最上位にいちいち浮上しません。 sageと入れないと一部の人からくどく思われる場合がありますので、普通はsageと入れます。 ちなみにsageは日本語の「下げ」をローマ字で書いたものです。他言語のsageの意味はありません。 725 名前：octave for windowsメモの筆者 mailto:sage [2009/02/03(火) 14:21:16 ] >>724 m(__)m　ありがとうございます。 726 名前：octave for windowsメモの筆者 mailto:sage [2009/02/03(火) 14:22:03 ] >>724 m(__)m　ありがとうございます。 727 名前：octave for windowsメモの筆者 mailto:sage [2009/02/03(火) 14:36:30 ] >>712 error: invalid column index = 2001ですね。 66行目　b(:,1)=[0;0];によって octave:3> size(b) ans 2 1 で1行2列の行列になってしまうわけです。 66行目　b=zeros(2,2001); とすれば大丈夫です。 しかし，For Loopの塊ですね。 octaveでは，遅くなりますね。 漸化式likeなので スパース行列を使えば高速化できるかも知れませんね。 ttp://www.nuce.nagoya-u.ac.jp/e8/Matsuoka/07OctaveNum/07OctaveNum.html 偏微分方程式なのでまともにやったら2重ループなのですが，遅くなるので， 漸化式を係数行列にしてスパースを使って高速化してます。 古い資料なので，octaveのバージョンが古いですが，基本は同じから 参考にしてください。 でもいまのままでも，答えは少し待てば出るので，まあ， 現状でいいかもしれませんがこういうテクニックも あることもしっておくといいですよ。 728 名前：octave for windowsメモの筆者 mailto:sage [2009/02/03(火) 14:41:57 ] >>723 ブリーフケースではありません。Webをおけるスペースです。 まあ，Yahooにいってみてください。 webがおけるのでファイルもおけて， <html> <a href='test2.m'>test2.m</a> </html> とでもして，index.htmlをつくり同じフォルダーに test2.mをいれとけばいいです。

---

---

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

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