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

1 名前：名無しさん＠お腹いっぱい。 [04/06/01 00:13] 知る人ぞ知る強力なデータ解析ツール octave ... unixを科学研究で使っている香具師には非常に心強い味方のはずだ。 matlab互換を主張しながら、微妙にバージョンアップを続ける octaveについて語ろう。 605 名前：名無しさん＠お腹いっぱい。 mailto:sage [2008/08/16(土) 10:59:45 ] >>596 なんとなくわかった事の覚え書き（チラシの裏みたいなもの） 1) そのnet2comを入れて、ノート内蔵モデムのCOM3をlocalhost:8003にしてStartする 2) octave のコマンドラインから次の順番でコマンドを入れる octave:1> more off sock1 = 1252 octave:2> sock1=socket(AF_INET,SOCK_STREAM,0) octave:3> server_info=struct('addr','localhost','port',8003) server_info = { addr = localhost port = 8003 } octave:4> rc=connect(sock1,server_info) rc = 0 octave:5> [msg_s,len_s]=recv(sock1,255); octave:6> msg='ATI0'; octave:7> msg(5)=13; octave:8> rc=send(sock1,msg) rc = 5 octave:9> [msg_s,len_s]=recv(sock1,20) msg_s = 13 10 53 54 48 48 48 13 10 79 75 13 10 len_s = 13 octave:10> num2str(msg_s,'%c') ans = 56000 octave:11> disconnect(sock1) rc = 0 606 名前：octave for windowsメモの筆者 [2008/08/18(月) 12:47:23 ] >>605 なるほど勘違いしていました。すねわち，forgeにsocket機能がある からNet2comを使用して，socket機能で通信しようというわけですね。 自分の書いたレスが恥ずかしい。　まっ，いっか勉強になったし。ありがとう ございます。 607 名前：名無しさん＠お腹いっぱい。 mailto:sage [2008/08/18(月) 19:13:18 ] >>606 ここの一番下の方にOctaveでのSocketの使用例があります ttp://www.nabble.com/Low-level-sockets-td4253125.html （しかしドキュメント少なすぎ） 昔、Parallel Octaveをやってる人がいましたが あれは独自実装のSocket経由で2台以上のPCのOctaveを接続して 1台をlistenで計算サーバーにさせて、もう一台からsendでデータを送って、 recvで計算結果を得るといみたいな事をやってるようでした。 サーバーが10個ぐらいあれば、データ転送のオーバーヘッドがあるにせよ 10個のパラメータを1度に試せるわけで、特定的な問題なら 一台で解く場合よりも速くなるのでしょう。 608 名前：octave for windowsメモの筆者 [2008/08/19(火) 15:18:25 ] ttp://www.nabble.com/release-3.0.2-to19037149.html のようにmaintainers向けのOctave-3.0.2がリリースされました。 まもなく，公式ページにもソースが載るでしょう。 また，新しいforgeも準備されつつあります。 これも，ここ１，２週間のうちに出ると思います。 609 名前：名無しさん＠お腹いっぱい。 [2008/08/19(火) 22:42:30 ] >>599,600 有難うございます。 分かりました。確かにawkや c + gnuplot は便利ですよね。 最後に一つだけ質問： octave だとdataから簡単にx,yzのデータから３D plot ができます。 ところが、gnuplotを用いるとデータの途中に適宜空行を入れないと いけないのですが皆さんどのような工夫をされているのでしょうか? 1.メインプログラムでデータ出力時に空行を適宜入れる。 2．gnuplotで読み込む前に空行を入れる。 3．他の技がある--> ご教示ください。 610 名前：名無しさん＠お腹いっぱい。 [2008/08/20(水) 16:53:10 ] >>603 です。　>>604様へ まだ、うまくいきませんが・・・ 丁寧な御回答ありがとうございました。 現在は　octave workshop を用いて、 付属のメモ帳を使って頑張っています。 611 名前：octave for windowsメモの筆者 [2008/08/20(水) 18:51:01 ] >>610 octave workshopはちょっとバグ多いので。 OctaveのGUIのフロントエンドを利用するならQtOctaveがいいですよ。 ttp://wiki.octave.org/wiki.pl?OctaveForWindows から飛んでください。確か表形式の行列入力あった気がします。 612 名前：octave for windowsメモの筆者 [2008/08/20(水) 18:55:41 ] >>609 すみません。この板の範囲ではないので，gnuplotには，いい掲示板があるのでそこいって聞いてください。 Gnuplot Q&A 掲示板 ttp://ayapin.film.s.dendai.ac.jp/cgi-bin/trees.cgi 613 名前：octave for windowsメモの筆者 [2008/08/22(金) 16:12:12 ] octave 3.0.2のソースが正式にリリースされました。 新しいForgeもまもなくでるので， 次のバイナリリリースは3.0.2+新forgeとなるでしょう。 ttp://www.gnu.org/software/octave/download.html 614 名前：octave for windowsメモの筆者 [2008/08/25(月) 08:48:37 ] 新しいOctave-Forgeがリリースされました。 ttp://www-old.cae.wisc.edu/pipermail/help-octave/2008-August/010651.html そのうち各種バイナリも整備されるでしょう。 615 名前：名無しさん＠お腹いっぱい。 mailto:sage [2008/08/25(月) 20:50:52 ] いつも乙です >> 中の人 616 名前：名無しさん＠お腹いっぱい。 [2008/08/27(水) 00:50:18 ] FreeBSDのportsがまだ新しくなってないな。 3.0.1のまま 617 名前：octave for windowsメモの筆者 [2008/08/27(水) 12:28:38 ] >>615 こちらこそ乙です。 618 名前：octave for windowsメモの筆者 [2008/08/27(水) 12:35:06 ] >>616 もうちょっと待ってくださいや。Forgeが出たばっかしだし。 3.0.xのシリーズの本体はminor changeなので問題ないけど， Forgeは新しいのがでると新しいパッケージがでて対応大変なんだよ。 Maintainerの苦労もわかってやって下せー。 おねげーしますだ。 619 名前：octave for windowsメモの筆者 [2008/08/28(木) 09:20:12 ] うーん。cygwinでForgeやってるけど前回できたoctgprがうまくビルドできん。 前のバージョンならうまくいくんだけど。 620 名前：octave for windowsメモの筆者 [2008/08/29(金) 11:18:04 ] >>619 自己レスです。Jaroslav Hajek氏（最近はこの方が，3.0.xの主たる コードをかいてます。JWEさんは，3.1ブランチのコードに専念されてます）が ミスがあったので，SVNにあげたと連絡ありで解決しました。 Cygwin gcc-dw2 版3.0.2+forge週明けにはリリースします。 621 名前：octave for windowsメモの筆者 [2008/09/01(月) 17:41:14 ] Octave 3.0.2 + forge 20080831 on cygwin (by gcc-dw2)をリリースしました。 ttp://www.tatsuromatsuoka.com/octave/Eng/Cygwin/ です。 622 名前：名無しさん＠お腹いっぱい。 mailto:sage [2008/09/02(火) 02:51:50 ] >621 あんたみたいになるにはc言語を勉強するの？ 623 名前：octave for windowsメモの筆者 [2008/09/02(火) 09:15:17 ] >>602 Cもそうですが，それよりもUNIXの開発環境の知識が問われます。 でも，一番肝心なのは外国人とメーリングリストでやり取りする度胸ですね。 624 名前：octave for windowsメモの筆者 [2008/09/02(火) 09:17:27 ] >>623 >>602は>>622の間違いです。すいません。 625 名前：名無しさん＠お腹いっぱい。 mailto:sage [2008/09/04(木) 11:03:54 ] 初心者ですが、質問させてください。 インパルス表示をさせようと、 プロット時に"^"を入力したのですが、 線ではなくなぜか△のマークが表示されます。 「それまでのプログラムに問題があるのかな。」 と、思い打ち直しましたが結局同じ症状でした。 違う簡単な式などをプロットしても同じ状況です。 この状況は打開できるでしょうか。 また、ほかにインパルス表示をさせる方法がありますでしょうか。 よろしくお願いいたします。 ちなみにOctaveのバージョンは3.0.1になります。 626 名前：octave for windowsメモの筆者 [2008/09/04(木) 18:27:27 ] Octave 3.0.2 + forge 20080831 on cygwin (by gcc-dw2)の パッケージミスの修正とforgeのftp-1.0.1をビルドしたので これをあわせたパッケージと差分をアップロードしてます。 ttp://www.tatsuromatsuoka.com/octave/Eng/Cygwin/ 627 名前：octave for windowsメモの筆者 [2008/09/04(木) 18:32:28 ] >>625 すみません。MSVC 版ならとGraphic backendがgnuplotかjhandleか 他のものならば，OSの種類とかgnuplotのバージョンとか明記しないと解答はつかないと思います。 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されるでしょう。

---

---

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

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