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

1 名前：名無しさん＠お腹いっぱい。 [04/06/01 00:13] 知る人ぞ知る強力なデータ解析ツール octave ... unixを科学研究で使っている香具師には非常に心強い味方のはずだ。 matlab互換を主張しながら、微妙にバージョンアップを続ける octaveについて語ろう。 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をいれとけばいいです。 729 名前：名無しさん＠お腹いっぱい。 mailto:sage [2009/02/03(火) 14:48:34 ] 日本のyahooの利用規約では公開目的ではなく単にファイルを置くだけなのは禁止されているはず。 国内で利用するんだったらwww.axfc.net/uploader/が良いかもしれないね。 流れてしまう速度はそんなに速くないはずだから。 730 名前：octave for windowsメモの筆者 mailto:sage [2009/02/03(火) 14:48:32 ] >>727 おっとこの講義資料の最後の回のものです。 731 名前：名無しさん＠お腹いっぱい。 mailto:sage [2009/02/03(火) 14:55:26 ] ほんとループが長い！まだ回ってる・・・ >>728,>>729 つまりヤッホーのほうは利用しないで>>729のほうがいいってことですかね？？ なにはともあれありがとうございます >>730 是非拝見させていただきます 732 名前：octave for windowsメモの筆者 mailto:sage [2009/02/03(火) 15:07:17 ] >>729 あっとそうでしたね。 <html> <body> <P>数値計算ソフトのソースの例(octave) <a href='test2.m'>test2.m</a> </body> </html> とでもしとかなきゃだめでしたね。 こうすれば，一応公開目的なのでおけますよね。 ただし，パスワードがかけられませんからね。 まあ，見られてもいい（つまり公開）であればいいのですよね。 733 名前：名無しさん＠お腹いっぱい。 mailto:sage [2009/02/04(水) 15:14:49 ] グラフを印刷する時、軸のラベルの文字の大きさを変えるのはどうするんですか？ print "hoge.eps" -deps -color -F:20 ってやると、目盛の数字は大きくなるけど、ラベルは小さい文字のままになりま す。 734 名前：octave for windowsメモの筆者 mailto:sage [2009/02/04(水) 16:20:56 ] >>733 xlabel('{\fontsize{25} x}'); 735 名前：名無しさん＠お腹いっぱい。 mailto:sage [2009/02/04(水) 16:31:14 ] 使ってるプロットツールによると思うんだけど 文字はそのままでグラフを小さく出力すると印刷時には相対的に文字が大きくなる 例えば俺はgnuplot使ってるから __gnuplot_set__ size 0.4 0.4; とかやって全体的に文字サイズを変えてしまうのはどうだろう 736 名前：octave for windowsメモの筆者 mailto:sage [2009/02/04(水) 16:40:56 ] __gnuplot_set__ size 0.4 0.4; この設定は3.0からは使えません。 ttp://www.gnu.org/software/octave/NEWS-3.The following plot commands are now considered obsolete and will be removed from a future version of Octave: __gnuplot_set__ __gnuplot_show__ __gnuplot_plot__ __gnuplot_splot__ __gnuplot_replot__ Additionally, these functions no longer have any effect on plots created with the Matlab-style plot commands (plot, line, mesh, semilogx, etc.).html　より The following plot commands are now considered obsolete and will be removed from a future version of Octave: __gnuplot_set__ __gnuplot_show__ __gnuplot_plot__ __gnuplot_splot__ __gnuplot_replot__ Additionally, these functions no longer have any effect on plots created with the Matlab-style plot commands (plot, line, mesh, semilogx, etc.). 737 名前：octave for windowsメモの筆者 mailto:sage [2009/02/04(水) 16:49:14 ] >>736 すいません。m(__)m。また，確認ミス。 とろろで「この設定は3.0からは使えません。」は使えませんはうそで， 正しくは，plot関数などによるplotには影響せず，__gnuplot_set__などは あくまでも，別のルーチンでgnuplotのコマンドをダイレクトに使用する描画 にのみ影響する。 が正しいです。 といいのは，octaveの次のバージョンから独自のgraphicsエンジンを 使うという構想があり，このようのgnuplotの直接命令は，plotなどの 関数に影響を与えないように変更されています。 738 名前：名無しさん＠お腹いっぱい。 mailto:sage [2009/02/04(水) 17:27:26 ] 今プログラムを作ってるんですが下のようなエラーがどうしても消えません どこがおかしいのでしょうか？ error: invalid column index = 0 error: evaluating binary operator `*' near line 35, column 14 error: evaluating binary operator `+' near line 35, column 12 error: evaluating assignment expression near line 35, column 4 error: evaluating for command near line 33, column 2 error: evaluating for command near line 29, column 1 error: near line 42 of file `/home/tse550u/soturon/ndynamicsg.m' ソース：www1.axfc.net/uploader/He/so/190986 739 名前：名無しさん＠お腹いっぱい。 mailto:sage [2009/02/04(水) 17:29:18 ] すいません、>>738のやつパスはｄｗです 740 名前：名無しさん＠お腹いっぱい。 mailto:sage [2009/02/04(水) 17:35:28 ] 新しいOctave入ってるPCで見てみたらマジだ 昔gset sizeでやってた頃に「これは使えなくなるから__gnuplot_set__ でやれ」 ってOctaveに怒られた時からいつかこうなるとは思ってたけど 741 名前：octave for windowsメモの筆者 mailto:sage [2009/02/04(水) 18:00:14 ] >>738 error: invalid column index = 0 が怪しい。行列やベクトルは最初に現れたsizeで初期化されますから， sizeを指定したいときは， A=zeros(N,M); などであらかじめ指定してください。 742 名前：octave for windowsメモの筆者 mailto:sage [2009/02/04(水) 18:06:08 ] >>740 今は，JhandleというMATLABそっくりな図を描くグラフィックバックエンドも ありますしね。MATLAB互換路線になっていますから，まあ仕方ないでしょう。 私は，データ吐き出して, gnuplotで直接図描きますから関係あまりないのですが^^;。 ところで，jhandleは，epsはファイル形式は，postscriptのみだったと 思います。 emfなどはpstoeditをつかって変換するしかないでしょう。 743 名前：octave for windowsメモの筆者 mailto:sage [2009/02/05(木) 09:44:33 ] >>741 error: invalid column index = 0 おっと　octaveはindexが1から始まるからでたエラーですねm(__)m。 744 名前：名無しさん＠お腹いっぱい。 mailto:sage [2009/02/05(木) 18:07:53 ] >>734-737 ありがとうございました。 自前のグラフィックスですかあ。私はgnuplotが結構好きなんで、微妙だなあ。 なんかOctave全体が重くなるのもいやだし。 745 名前：名無しさん＠お腹いっぱい。 mailto:sage [2009/02/06(金) 10:43:06 ] 筆者さん、Octaveの授業もやってるんですね。わたしも筆者さんと同じ分野の人 間です。授業でOctave使おうとしたことがあるんですが、他の教員から「そんな ブラックボックス使わせても教育にならん」と反対されました。かといって、現 行の授業でやってるようにC言語だと、病的な文法の瑣末なところで引っ掛る学生 がいっぱいいて、シミュレーションそのものを体感させるところまで行かないん ですよね。 "プログラミング演習"の例題をOctaveで解かせるのは、研究室の学生にさせたこ とがあります。"プログラミング演習"は過去におそらく3種類出版されているかと 思いますが、わたしは80 年代末に出たBasic版が一番のお気に入りでした。基本 算法の解説が一番詳しかったように記憶しています。現在はFORTRAN版もBasic版 も入手困難なようですね。FORTRAN版は手元にあるんですが、BASIC版はどっか行っ ちゃいましたorz。今出まわってる改訂版は、うーん、解説は適当だし、Cと FORTRAN両方のコード載せるのも場所の無駄だし、おまけに今時FORTRAN77だし。 746 名前：名無しさん＠お腹いっぱい。 mailto:sage [2009/02/06(金) 11:12:37 ] いきなりOctaveやらせると感動がないじゃん 「スゲーこの言語スゲー!!俺こんなに手抜いちゃって良いの?」 っていう感じで 747 名前：名無しさん＠お腹いっぱい。 mailto:sage [2009/02/06(金) 13:12:13 ] FORTRANでもいいんじゃね 748 名前：名無しさん＠お腹いっぱい。 mailto:sage [2009/02/06(金) 13:26:18 ] Matlab/Octaveの醍醐味は初級プログラミングのレベルではなくて 普通にループを回すと何時間も何日もかかるループ処理が 行列演算で１発で出るという所なので Excelで済む程度のデータなら必要無いでしょうね 逆に例えばモンテカルロ計算なんかだと、 一般的なプログラムではループ10万回ぐらい回して そのループの中で乱数を発生させて、関数に入れて結果を見る事になりますが Matlab/Octaveでは最初に、10万個の乱数配列を作り 行列演算化した関数に入れて1発で結果が出るみたいな感じなので 発想を根本的に変えないといけません 749 名前：octave for windowsメモの筆者 mailto:sage [2009/02/06(金) 18:05:40 ] >>748 そうですね。いかに，行列で問題を記述するかにかかってて， 偏微分の差分法の陰解法では良くガウスザイデル法とかつかいますが， Octaveなら行列で方程式書いて左除算で済ませる。 係数行列は疎行列だからsparseを使って高速化するとか 発想が違いますよね。 まあ，ブラックボックスっいう方いますが，硬い微分方程式なんぞ ルンゲクッタで解こうとすると大変で，lsode使うとあっという間に 解けますよね。 適応的自動積分でも同じだとおもいます。 シンプソン公式を自分でCで書いて使うような時代ではないと思ってます。 ブラップボックスでもいいから，数値計算ができることの意味を分からせる のが大事だとおもいます。 750 名前：octave for windowsメモの筆者 mailto:sage [2009/02/10(火) 09:42:13 ] 開発版octave 3.1.52がでました。 ttp://www.nabble.com/Octave-3.1.52-available-for-ftp-to21923205.html#a21923205 Octaveの次のMajor Version Upはver 3.2であり，octave-3.1.xxは開発版です。 その3.1.52がでました。 gnuplotは一応defaultのgraphic backendでopenglによるgraphic backend を機能させるためにはいろいろライブラリィが必要なようです。 詳しくは以下を ttp://www.nabble.com/Octave-3.1.52-available-for-ftp-to21923205.html#a21923205

---

---

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

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