MATLABプログラミング　質問箱　その2

1 名前：デフォルトの名無しさん [2006/11/10(金) 11:28:22 ] MATLAB　日本語版 www.cybernet.co.jp/matlab/ MATLAB 本家 www.mathworks.com/ Yahoo MATLAB dir.yahoo.co.jp/Science/Mathematics/Software/MATLAB/ MATLAB入門 feedback.mech.fukui-u.ac.jp/matlab/matlab0.htm MATLAB チュートリアル www-isys.eei.metro-u.ac.jp/~sho/pub/ttmatlab/Index.html 555 名前：デフォルトの名無しさん mailto:sage [2008/05/22(木) 01:41:05 ] その質問は、もう質問の体をした荒らしみたいなもの。 google もまともに使えんのに matlab で何をするつもりなのか、答え給え。 ぼくぁいつもマジレスなんだ。 556 名前：デフォルトの名無しさん mailto:sage [2008/05/22(木) 14:24:57 ] たいしたことじゃないんですが、 a = [0 0 1 1 1 0 0 0 0 1 1 1]; とあったときに、3:4, 9:11 のように連続点を抜き出すのに楽な方法ってないですかね？ b = find(a); として b = [3 4 9 10 11]; の状態から 3:4, 9:11 と２つに分ける方法でも良い。 今はループ回しているんですが、matlab だし、何か楽な方法がないかなぁと。 ちなみに、この 3:4, 9:11 の表現は、 c = [3 4; 9 11]; のようにしているのだけれど、別になんでもいいです。 557 名前：デフォルトの名無しさん mailto:sage [2008/05/22(木) 17:36:53 ] 伝達関数からステップ入力に対する オーバーシュートや静定時間を求める関数ってありますか？ 558 名前：デフォルトの名無しさん [2008/05/22(木) 21:12:07 ] >>556 aを、微分（差分）してみるといいかもしれない 559 名前：デフォルトの名無しさん mailto:sage [2008/05/22(木) 23:42:58 ] >>558 おぉおもしろいアイデアだ。 D = conv(a, [1 -1]); starts = find(D == 1); ends = find(D == -1) -1; ind = [starts(:) ends(:)]; でいけました。うほっ。とんくそ。 560 名前：デフォルトの名無しさん mailto:sage [2008/05/23(金) 22:01:47 ] >>557 stepinfo ただしR2006a以上のControl System Toolbox 561 名前：デフォルトの名無しさん mailto:sage [2008/05/28(水) 20:12:10 ] >>560 ありがとうございます。でもうちのじゃ使えなかったみたいです。 あきらめます。 562 名前：デフォルトの名無しさん mailto:sage [2008/05/29(木) 19:44:51 ] 過去にMATLABにLISPのmap関数の様なものがあるかと質問したものです． その後，mexを作って解決しましたが，MATLABには標準で，arrayfunやcellfunといった同様の関数があることが分かりましたので一応報告しておきます． 563 名前：デフォルトの名無しさん mailto:sage [2008/05/30(金) 03:32:38 ] ちゃんと聞けば教えてあげたのに。 564 名前：デフォルトの名無しさん mailto:sage [2008/05/30(金) 14:18:05 ] すいません．コンソールで値の入力を要求して，その結果を変数に代入するような文はどう書くのでしょうか？ Please imput value of x : と表示させてxの値を入力させたいのです． 565 名前：デフォルトの名無しさん mailto:sage [2008/05/30(金) 14:18:29 ] inputですた 566 名前：デフォルトの名無しさん mailto:sage [2008/05/30(金) 22:32:32 ] すいません．自己解決しました． 567 名前：デフォルトの名無しさん [2008/06/01(日) 21:22:19 ] 図をepsで保存してイラレで読み込むと マイナスが見えなくなるのですが これってどうしたらいいのですか？ 568 名前：デフォルトの名無しさん [2008/06/03(火) 14:59:19 ] feval(fnc,r) でfncに'sin'とかは入れてsin(r)を計算できるのですが、 x^2などのxにrを入れて計算するのはどうすればいいですか。 569 名前：デフォルトの名無しさん mailto:sage [2008/06/03(火) 17:23:50 ] >>568 feval(@(x) x^2, r) 570 名前：デフォルトの名無しさん [2008/06/04(水) 18:58:14 ] matlabで多目的整数計画問題って解けるんですか？ fminimax()の変数が整数版みたいなやつが欲しいんです。 fminimax()の条件をいじって、変数を整数に縛ろうと思ったんですけど、 やり方もわかりませんし・・・ 571 名前：デフォルトの名無しさん [2008/06/04(水) 21:26:37 ] 縦軸の値から横軸の値を表示することはできますか？ たとえば、y=2xのグラフをxを1〜10まで0.1刻みで描いたとします。 そのとき、y=5のときのx値を知りたいわけです。 刻み幅が大きく合致する値がない場合はどうなるのでしょう？ 572 名前：デフォルトの名無しさん [2008/06/04(水) 21:40:04 ] >>571 interp1dをつかったら？ 573 名前：デフォルトの名無しさん mailto:sage [2008/06/04(水) 23:27:43 ] >>571 推定するしかねーだべが。 574 名前：デフォルトの名無しさん mailto:sage [2008/06/05(木) 02:21:50 ] matlabの関数の説明が日本語で書かれているサイトないでしょうか？ できればURL張ってくださいorz 575 名前：デフォルトの名無しさん mailto:sage [2008/06/05(木) 23:35:26 ] matlab の中の人たちって何者？ 特に toolbox とか数学的に難しいものもあるから、 それを実装できる知識を持っているが、 研究者やらずに matlab 開発しているとか良くわからない。 ph.d 持ちの人たちなのかな？ 576 名前：デフォルトの名無しさん mailto:sage [2008/06/06(金) 19:13:20 ] 大丈夫、私も数学的知識はないけど数学関数の実装は(数学屋に書いてもらったアルゴリズムに基づいて)行っている。 577 名前：デフォルトの名無しさん [2008/06/07(土) 17:31:53 ] >>574 「ＭＡＴＬＡＢ　ヘルプ」でぐぐれ 578 名前：デフォルトの名無しさん mailto:sage [2008/06/13(金) 02:40:44 ] >>575 「ＭＡＴＬＡＢ　ヘルプ」でぐぐれ 579 名前：デフォルトの名無しさん mailto:sage [2008/06/13(金) 14:51:55 ] Unix版とWin版で何か、違いはありますか？ 580 名前：デフォルトの名無しさん [2008/06/16(月) 09:27:12 ] 変わんない 581 名前：デフォルトの名無しさん mailto:sage [2008/06/16(月) 20:48:40 ] 「ＭＡＴＬＡＢ　ヘルプ」でぐぐれ 582 名前：デフォルトの名無しさん [2008/06/25(水) 20:53:55 ] 配列のindexing をC/C++みたいに 0 始まりにすることってできますか? 583 名前：デフォルトの名無しさん mailto:sage [2008/06/26(木) 10:32:35 ] 「ＭＡＴＬＡＢ　ヘルプ」でぐぐれ 584 名前：デフォルトの名無しさん [2008/06/26(木) 19:17:32 ] >>583 お前は首くくれ 585 名前：デフォルトの名無しさん mailto:sage [2008/06/26(木) 21:14:55 ] >>582 できないと思う 586 名前：デフォルトの名無しさん mailto:sage [2008/06/26(木) 22:06:29 ] >>584 カレーはククレ 587 名前：デフォルトの名無しさん [2008/06/27(金) 08:35:39 ] 金かけたくないなら、pylabいいよ 588 名前：デフォルトの名無しさん [2008/07/01(火) 04:20:23 ] h=[1 0.8 0.4 0.2]; x=[1 0 0 0 0]; y=convol(h,x) subplot(4,1,1) xa=1:7; plot2d(xa,y,style=-2,rect=[0 0 8 2]) のプログラムを実行すると plot2d: first and second arguments have incompatible dimensions とエラーが出ます。 どう対処すればいいですか？誰か教えて下さい。 589 名前：デフォルトの名無しさん [2008/07/01(火) 06:45:02 ] xaとyの数が違うんじゃないの？ plot2d(xa,y(1:length(x) ... にしたら？ 590 名前：しっとり [2008/07/04(金) 17:42:43 ] 突然すいません。 MATLABでLeonhardの軌跡ってどのように描けますか？教えて下さい。 591 名前：デフォルトの名無しさん [2008/07/04(金) 18:01:39 ] こんにちはいきなりですが質問です。 例えば100行100列のランダムな値をもつ行列から上位30点の要素を抽出するにはどうしたらいいでしょうか？ 誰か教えてください 592 名前：デフォルトの名無しさん mailto:sage [2008/07/04(金) 20:23:06 ] a=rand(100:100); b=create_set(a); c=b(length(b)-29:end) 593 名前：デフォルトの名無しさん [2008/07/05(土) 10:59:13 ] >>592 素早い回答ありがとうございます 早速やってみます 594 名前：デフォルトの名無しさん [2008/07/07(月) 04:26:21 ] 質問です　次のようなアプリを作りたいのですがMATLABを使い始めて日が浅いので見当がつきません、 作り方・コードを教えて下さい　お願いします。 フォルダの中に「ＰＣ」等の２文字が書かれたbmp画像がN枚ある 一枚一枚を任意の場所で左右にカットする（ＰとＣが分かれるように自分で分ける線の位置を決める） それぞれ○○○○○L.bmp　○○○○○R.bmpとリネームして同じフォルダに保存する この作業をN回繰りかえし終了 図：www.dotup.org/uploda/www.dotup.org22131.jpg 595 名前：デフォルトの名無しさん [2008/07/07(月) 18:25:14 ] あんまり甘やかすのは、アレな気がするけどこれでできるはず… 問題は、分割位置が固定されているのか、問題毎に違うのか…。 プログラムで分割位置を判別は簡単じゃないと思うけど % Image Prosessing Toolboxがインストールされていて % 画像ファイルはRGB形式の仮定で行きます % 違う場合は若干書き換えてください filename = { 'image1.png', 'image2.png' }; % ファイル名を書いたリスト l=length(filename); for i=1:1 I=imread( filename{i} ); % 画像の読み込み I=rgb2gray(I); % 画像をグレースケールに変換 [m,n]=size(I); sep=100; % 分割位置を指定 % 画像の保存領域を作成 IR=zeros(m,sep); IL=zeros(m,n-sep); % 画像を転写 IR = I(:,1:sep); IL = I(:,sep:n); % 書き出し imwrite( IR, strcat([filename{i},'R.bmp']), 'bmp' ); imwrite( IL, strcat([filename{i},'L.bmp']), 'bmp' ); end 596 名前：594 [2008/07/08(火) 00:29:39 ] >>595 ありがとうございます >>問題は、分割位置が固定されているのか、問題毎に違うのか…。 画像ごとに違う位置で分割したいので分割位置を一枚一枚その度に自分で決めれるようにしたいのです 説明が分かりにくかったかもしれません、↓実際のファイル内はこのような感じです www.dotup.org/uploda/www.dotup.org22827.jpg これをそれぞれ「退」と「院」に分割して保存したいのです しかし画像によって分ける位置が違うので手動で分割位置を決めれれば…と思っています >>% Image Prosessing Toolboxがインストールされていて >>% 画像ファイルはRGB形式の仮定で行きます Image Prosessing Toolboxについて簡単に調べてみましたが 多分インストールされていると思います、MATLABのタイプ（？）は[R2007b]です。 画像ファイルはモノクロビットマップです >>filename = { 'image1.png', 'image2.png' }; % ファイル名を書いたリスト 画像が1000枚近くあって名前もバラバラなのですが・・・全ての名前を書き出さずにできないでしょうか cut.mと名付けこのまま実行してみたところ、このようなエラーが出ました。一応載せておきます ??? エラー ==> imread at 315 ファイル "image1.png" は存在しません エラー ==> cut at 5 I=imread( filename{i} ); % 画像の読み込み 597 名前：デフォルトの名無しさん mailto:sage [2008/07/08(火) 07:31:40 ] コード中で.pngとファイル名を指定しておきながら.bmp画像しか用意してなかったんだろ そのまま実行するという暴挙に出てエラーが出ましたとか言ってしまう神経が以下略だぜ 598 名前：596 [2008/07/08(火) 13:26:40 ] >>597 その通りですね…冷静に考えればなんと当たり前な… image1.pngとimage2.pngを用意して実行したところ、今度はこのようなエラーが出てしまいました ??? エラー ==> rgb2gray>parse_inputs at 82 MAP must be a m x 3 array. エラー ==> rgb2gray at 35 X = parse_inputs(varargin{:}); エラー ==> cut at 6 I=rgb2gray(I); % 画像をグレースケールに変換 599 名前：595 [2008/07/09(水) 16:16:09 ] えっと、答える私も私なんだけど少しだけ叱らせて。 ここは初心者に対して伝統的に厳しいスレではないけれど、それを考えてもあなたのレベルはちょっと酷い。 私も大した事は言えた立場ではないけれどMatlabが可哀想なのでせめてもう少し勉強してあげて。 一応回答するけど、画像を一括指定したい場合は、知恵を絞りましょう。 指定したディレクトリにある画像をすべて取得するならこんな感じ % 指定したディレクトリにあるファイルのファイル名をリストで取得 imgdir = 'sampleimages'; % 画像がおいてあるフォルダの名前 files = dir(imgdir); % 指定したフォルダのファイル名を取得する filename = {}; % ファイル名を書いたリスト % ディレクトリのファイルからファイル名を取得 for i=1:length(files) if exist(imdir+files(i).name) filename{i}=imdir+files(i).name; end end バグがおこっている行はRGB画像をグレースケールに変換するものなので画像のフォーマットの問題。 該当行を I=im2uint8(I); % 画像をバイト型からグレースケールに変換 とかに差し替えてやればよいと思う。 あとは画像の分割位置を手動で決定するにはGUI画面を作らないといけないね。 ここからが初心者には大変なんだろうけど、面倒くさいな… 600 名前：デフォルトの名無しさん mailto:sage [2008/07/10(木) 13:22:44 ] アルゴリズムを考えることが レポートの課題なり、研究のポイントなんだから、 そこまでここで教えてもらおうなんて思わないように。 601 名前：今風若者気質 mailto:sage [2008/07/10(木) 13:56:22 ] >>600 何言ってるんですか、レポートにしろ卒研にしろ、抜ける手は抜くのが一番じゃないですか。 602 名前：デフォルトの名無しさん mailto:sage [2008/07/10(木) 14:20:33 ] 勉強するために大学行くんだろ。 就職のためなら専門学校いけよ。 603 名前：デフォルトの名無しさん mailto:sage [2008/07/10(木) 14:26:22 ] 頭は使わないと劣化していくよ 教授を見ていれば分かるだろ 604 名前：今風若者気質 mailto:sage [2008/07/10(木) 17:38:02 ] >>602 何言ってるんですか、専門学校の就職率なんてカスですよ。 605 名前：595 [2008/07/11(金) 12:51:43 ] uproda11.2ch-library.com/src/11100794.zip.shtml あぷろだ初めて使うんだけどこれでいいのかな？ pass:1111 削除キー:1111 使い方 ・samplesフォルダに全ての画像を移す ・imageselect.mをmatlabから起動 ・開始ボタンを押した後、画像の分割位置をクリックで選択 見た目に良く分からないかもしれないけどちゃんと分割された画像が保存されています。 手抜き仕事なのでバグあったら頑張れ 606 名前：デフォルトの名無しさん mailto:sage [2008/07/11(金) 21:26:56 ] >>605 おまえさんのアホさに感動を覚える 607 名前：605 [2008/07/14(月) 13:46:22 ] あっはっは、出来が酷いのは手抜き中の手抜きだから勘弁してやってください。まぁ、赤の他人に対しては良くやったよね… ってそういう意味じゃないか… 一応やっている分野の近い相手なので、何となく情が働きまして… matlabのGUIプログラミングはVBなどと比べても決して難しい物ではないけれど、あのレベルで習得させてやらせるのは辛いしね 608 名前：デフォルトの名無しさん mailto:sage [2008/07/15(火) 08:28:19 ] GUIにせずとも、マウス入力で区域指定できるのではないかと 609 名前：デフォルトの名無しさん [2008/07/15(火) 23:33:57 ] filter関数がＣ言語でかけねぇ。 誰か作ってないかなぁ 610 名前：デフォルトの名無しさん [2008/07/16(水) 00:42:37 ] どんなフィルター？平滑化なんかの２次元フィルターってこと？ 611 名前：デフォルトの名無しさん mailto:sage [2008/07/16(水) 02:46:42 ] matlab の filter 関数でしょ。 612 名前：デフォルトの名無しさん [2008/07/16(水) 07:39:55 ] そうそう。matlabのbutter関数で係数を求めておいて、 それをつかってfilter(b,a,signal)って使うことのできる関数を Cで作りたいわけ 613 名前：デフォルトの名無しさん mailto:sage [2008/07/16(水) 10:06:30 ] 簡単じゃん。 614 名前：デフォルトの名無しさん mailto:sage [2008/07/24(木) 05:49:41 ] matlab での代表的なコーディング規約ってないの？ built-in 関数みてると大文字はおろか _ すら使ってないのがほとんどな気がするんだけど、そういうもん？ ちょっと、これじゃやりにくい。 615 名前：デフォルトの名無しさん [2008/07/24(木) 06:52:25 ] 昔のFORTRANの流れを汲んでるから、ひょっとすると関数名は６文字までだったりしてね。 個人的には _ を使った関数名は自作関数というぐあいに区別しているけど。 616 名前：デフォルトの名無しさん [2008/07/24(木) 18:29:52 ] MATLAB(R2008a)で関数自体の引数の数を得る方法はありますか? 無名関数を引数としてその無名関数自体の引数の数を返す関数を作りたいのですが。 また，関数自体はまだオブジェクト化されていないのでしょうか？ その様なクラスも見当たらないですし。 617 名前：デフォルトの名無しさん [2008/07/24(木) 21:09:47 ] 関数自体ってどういう意味？引数の数はnarginだろ普通 618 名前：デフォルトの名無しさん [2008/07/24(木) 21:44:37 ] >>617 >引数の数はnarginだろ普通 それはある関数呼び出しに関する引数の数だろ。 phpspot.org/blog/archives/2008/04/javascript_65.html これで分かるかな? 619 名前：デフォルトの名無しさん [2008/07/24(木) 23:09:42 ] 無理じゃね？ matlab は数値演算はいいけど、そういう補助機能が弱いよな。 >関数自体はまだオブジェクト化されていないのでしょうか？ そんなことする気あるの？なんかで言ってた？ 620 名前：デフォルトの名無しさん [2008/07/25(金) 02:27:23 ] というか、MATLABをプログラム言語としてみるとこれ以上ダサい文法の言語は正直無いと言えると思う。 最初に行列を格好良く計算できる機構と入れたのは良いけど、やっぱり後で条件分岐やループが必要だと入れたあたりから破綻が始まっている そこからは無計画と言うしかない文法の拡張が行われている 比較的最初に気づくものでも 関数作りたかったら別にmファイルを作ってそこに書かないといけないとか 関数は最後にendで括っても括らなくても良いとか 何故か決して使いやすいとは言えない関数名がC言語風とか 行列計算の書式しか無かった頃だけ直感的で利点があった、何故かインデックスが1から始まるとか 正直、簡潔な書式の行列計算＋強力なplot機能＋インタプリタですぐ結果が見える使いやすさなのだと思う 本当は文法関係は言語の進化に合わせて再設計が必要な時に来ているのだと思う。もう無理だけど…。 >>614 mlint関数で文法チェックとかどうでしょう？ 621 名前：デフォルトの名無しさん mailto:sage [2008/07/26(土) 17:50:46 ] スライドバーのコールバックで処理するとき、 ドラッグが終わってバーの位置が確定した時だけ処理を行うことはできないでしょうか? バーが変化するごとにいちいち変化してると重いですorz 622 名前：デフォルトの名無しさん mailto:sage [2008/07/29(火) 12:20:46 ] >>621 アルゴリズムの問題だと思う。 ちょっと考え直したら？ 623 名前：MATLAB初心者 mailto:sage [2008/07/30(水) 00:02:55 ] MATLABを使って 下記のようなニュートン法でα(ステップサイズ)を 適応的に求めるプログラムを作りたいのですが、 一向に分かりません。FORTRAN77を今まで使っていましたが 私のいる会社では、今年からMATLABになってしまったのです。 申し訳ありませんがこの掲示板をごらんになっておられる皆様 ご教授御願いいたします。 x(k+1)=x(k)−αinv(H(x(k)))g(x(k)) g(x(k))はf(x(k))の1次微分 H(x(k))はf(x(k))の2次微分 624 名前：デフォルトの名無しさん mailto:sage [2008/07/30(水) 02:42:36 ] 宿題ぐらい自分でやろうぜ つーか、ほとんどmatlab関係なくね？ 625 名前：デフォルトの名無しさん [2008/07/30(水) 03:43:34 ] FORTRANならわかるの？なのに matlab だとわからないの？ 626 名前：デフォルトの名無しさん mailto:sage [2008/07/30(水) 12:54:20 ] アルゴリズムが分からないのなら 教科書を読め 627 名前：デフォルトの名無しさん mailto:sage [2008/07/30(水) 23:23:53 ] matlab の関数名、変数名は lower_case のように小文字と _ を使っている人が多いようだけど、 なにか理由があるのかね？ なんとなくそういうスタイルのテキスト多いのかな？ datatool.com/downloads/matlab_style_guidelines.pdf の人は変数名には lowerCamel を押しているようなんだけど。 ＃関数名は _ もなしの lowercase 押してるけど、これはないよな。よみづらい。 628 名前：デフォルトの名無しさん mailto:sage [2008/07/30(水) 23:45:26 ] matlab を起動すると勝手に IME が日本語入力モードになります。 matlab で日本語入力モードになる意味がわからないので、 止めたいのですが、何か手法はあるのでしょうか？ 629 名前：デフォルトの名無しさん mailto:sage [2008/07/31(木) 11:09:24 ] ここのスレ、性格悪い低脳多いな。 630 名前：デフォルトの名無しさん mailto:sage [2008/07/31(木) 13:49:10 ] >>627 空白をあらわすのにlower caseを使う人がアメリカ流の人には多いと思う 631 名前：デフォルトの名無しさん [2008/07/31(木) 23:33:08 ] matlabのmファイルからwavをsound関数を使って再生しているのだけど、 途中で止めることができない（；＿；）。 途中で止めることはできないのでしょうか？ コードは以下 [y, Fs] = wavread('test.wav'); sound(y, Fs); 632 名前：デフォルトの名無しさん mailto:sage [2008/08/01(金) 05:27:56 ] >>630 lower_case と lowerCamel どっちがいいかなぁ・・・ lower_case は古い流行な気がしないでもないんですよね。 UNIX というか GNU C コーディングスタイル的な。 633 名前：デフォルトの名無しさん mailto:sage [2008/08/02(土) 07:11:06 ] >>631 最初から止めたい位置がわかっているのなら、 マトリックスの部分指定。 とにかく止めたいのなら、CTRL+C 634 名前：デフォルトの名無しさん [2008/08/02(土) 17:41:28 ] >>632 JavaScriptなんかの場合は最初は自分で書いたライブラリは_で区切るなんて風習もあったけど そのうちライブラリを自分で配る人達が出てきてそれが普及してからはlowerCamelとJavaScriptの言語の標準に従う様になった。 最初は遊びで書いた関数も、そのうち人に配る可能性が無いわけでもないし、MATLAB標準に黙って従うのが一応正しい姿勢となるのでは？ 635 名前：デフォルトの名無しさん mailto:sage [2008/08/02(土) 17:43:38 ] MATLABの大文字小文字の区別って、完全に可能なんだっけ？ 636 名前：デフォルトの名無しさん mailto:sage [2008/08/03(日) 05:07:36 ] >>634 でも、matlab 標準というと、例えば bankrupt company が bkrptcmpny みたいなかんじじゃない？ これはひどいでしょ？実際は、ここまでひどくはないけど、 _ もなにもないし、まぁひどい。 >>635 ファイル名 == 関数名、だから windows では大文字小文字が区別つかない。 けど、pca という関数がある時に、Pca を作るようなセンスの人はいないだろぉと思うし、 このために、lowerCamel は絶対ダメ！っていうのはもったいない気がする。 637 名前：デフォルトの名無しさん mailto:sage [2008/08/03(日) 05:18:56 ] >>634,635 で、 toolbox 作っている人をみると、関数名は、 lower_case みたいに書いている人が多いので、そちらが標準かと思ったら、 ttp://www.mathworks.com/access/helpdesk/help/techdoc/matlab_prog/f4-39683.html を見ると関数名が makeParabola みたいに lowerCamel になってたりするんだよね。 でも、やっぱりファイル名の問題で、lower_case にしようかと思ったとすると、 変数名も lower_case にしたくなってくるわけだけど、 ttp://www.mathworks.com/access/helpdesk/help/techdoc/ref/persistent.html とかみると変数名は lowerCamel だったりするんだよね。 今、リンクないけど、lower_case な関数名、変数名を使っているマニュアル例もあった。 というわけで、matlab 標準がどっちつかずなもんで、むきーという状況。 638 名前：デフォルトの名無しさん mailto:sage [2008/08/03(日) 05:29:17 ] なんだか長くなって申し訳ないけど、 C も標準は、略称使いの _ もなしの lowercase と思っていいいと思うけど、 OpenCV とか OpenGL は lowerCamel を使っている。 個人的には lowerCamel を使いたいのだけれど、 ファイル名 == 関数名の所で悩んでいる。 関数名は lower_case のような形で書いて、変数名は lowerCamel ってアリなセンスなのかな？ ところで、private 関数的なネスト関数の名前は最初に _ でも付けようか、 と思ったら matlab では変数、関数名の最初を _ にはできないのね。がっかり。 639 名前：デフォルトの名無しさん mailto:sage [2008/08/05(火) 05:03:28 ] お前プログラム書くの遅いだろ 640 名前：デフォルトの名無しさん [2008/08/05(火) 13:46:17 ] >というわけで、matlab 標準がどっちつかずなもんで、むきーという状況。 で、個人的意見なのですが、基本的にはMATLABのユーザーは物理系の人で、プログラミングの専門家やほかの言語に詳しい人は少ないのでそういうリテラシは少し低めなのではないでしょうか？ OpenCV とか OpenGLとかはC++の文化を意識してああいう風になっているんでしょうねぇ オブジェクト指向を基本的には出来ないとはいえ、一応matlabの関数名の文化は、Cからの影響なのを考えるとlowerCamelの記法は可なのではないでしょうか？ > ところで、private 関数的なネスト関数 Lispでは「高階関数」と呼んでいるのでこれに合わせると良いのでは？ 641 名前：デフォルトの名無しさん mailto:sage [2008/08/05(火) 21:59:29 ] >>640 >lowerCamelの記法は可なのではないでしょうか？ やっぱり問題は関数名なんだよね >Lispでは「高階関数」と呼んでいるので 高階関数ってラムダ関数ちゅうか関数ポインタを受け取ったり 返したりって関数でしょ？fminsearch みたいな。 ネスト関数ってのはこっちね。 ttp://www.mathworks.com/access/helpdesk/help/techdoc/matlab_prog/f4-39683.html 642 名前：デフォルトの名無しさん mailto:sage [2008/08/06(水) 11:40:35 ] 速く走らせるには工夫が要るけど それほど大きなパッケージを作るわけでもなかったら、 書法は気にしなくてもいいのでは、って思う。 おれはそのレベルのユーザー。 643 名前：デフォルトの名無しさん mailto:sage [2008/08/07(木) 02:02:31 ] z = (2*cos(x.)+1)(2*cos(y.)+1) のグラフを書きたいのですが、上記の式で エラー: 一致していないか、予期せぬ括弧か大括弧です。 とでます。x,yの値は x=-3:0.1:3; y=-3:0.1:3; このように設定しております。 どなたか教えてください。 644 名前：デフォルトの名無しさん mailto:sage [2008/08/07(木) 02:35:46 ] すみません。.の意味を取り違えてました。 645 名前：デフォルトの名無しさん [2008/08/07(木) 12:05:48 ] 質問があります エクセルデータを読んでグラフにしたいのですが、ファイルネームが ABC-01.xls ABC-02.xls . . . ABC-09.xls と続くので、forを使って for I=1:9 data=xlsread('ABC-0'num2str(I)'.xls'); hold on plot(data) end としたいんです。 でも2行目にどうやっても忠告がでてしまいます。 文法うんぬん以前に考え方自体が間違ってるのでしょうか？ 646 名前：デフォルトの名無しさん mailto:sage [2008/08/07(木) 12:28:00 ] こんなシンプルなのに考え方もなんも無いって。単純に文法の問題。 「こう書けるだろう」と勘で書いたところがあるでしょ? そういうところが本当に正しいかどうかチェックしてけばいい。 help strings 647 名前：デフォルトの名無しさん [2008/08/07(木) 12:55:00 ] やっぱりそうですよね。 プログラム自体が始めてから数週間なのでグダグダです。 ありがとうございました 648 名前：デフォルトの名無しさん mailto:sage [2008/08/07(木) 15:07:13 ] アイディアとしてはいい線いっているので、 あと少しだけ、知れば、サクサク作業できるようになると思う。 649 名前：デフォルトの名無しさん mailto:sage [2008/08/07(木) 20:58:54 ] すいません、ぐぐっても情報が少なすぎたので質問させてください。 JMatLinkというJAVAとMatlabをつなげるツールをインストールしたいのですが、 www.held-mueller.de/JMatLink/install.html に書いてある通りにしても上手くいきません。 上のリンク先に書いてあるように, c:\WINDOWSにJMatLink.dll をコピーし, c:\JMatLinkでDos窓を開き、 %java TestGui と打ち込んだのですが, java.lang.NoClassDefFoundError というエラーが出てテストウインドウが開きません. OSはWindows XP, matlabのバージョンは7.3.0、 JAVAはjdk1.6.0_05を使用しています. どなたか、JMatLinkを使用している方がいましたら、 何が悪いのか教えて頂けないでしょうか。 650 名前：デフォルトの名無しさん [2008/08/07(木) 23:14:34 ] Matlab 2007bと、Viual Studio 2005の質問です。 既存のC++で作られたライブラリの関数をラップして （要は、mdlOutputs内で既存の関数をコールする）、 C++ MEX S-functionを作ろうとしています。 Visual C++ 2005上で作成／ビルドして、*.mexw32を作っているのですが、 simulink上でこれを実行しようとすると、 「無効なMEX-ファイル*****; 指定されたモジュールが見つかりません。」 というメッセージが出ます。 サイバネットが出している、C-MEX-S-ファンクションチュートリアル （www.cybernet.co.jp/matlab/library/library/file/T/TT034.pdf?id=TT034） の、「5.1 既存Cコードの利用」に、 1. C MEX S-Function 内で外部関数のextern 宣言を行う。 2. コンパイル時に外部ファイルもコンパイルする とあったので、1.のみ行って、なんとかビルドはできたのですが、まだ同じ エラーメッセージが出ます。 やはり2.もしないといけないようなのですが、VC++で外部関数も一緒にコンパイルする 方法が、わかりません。 上記のマニュアルに書いてある方法は、Matlabコマンドラインからの すごく単純な構成のもののコンパイルなので、今やろうとしている、 たくさんのヘッダファイルやらdllやらをリンクしているモジュールの コンパイルは、とても難しそうです。 651 名前：デフォルトの名無しさん mailto:sage [2008/08/07(木) 23:25:33 ] >>647 >>646に書いてあるように help strings すれば 必要なことが書いてるよ。 MATLABはちょっとクセがあるから他の言語を知ってる人が いろいろひっかかるポイントがあると思うけど ヘルプはけっこうちゃんとしてると思うから どんどんヘルプ読んでいけばすぐ覚えるよ。 652 名前：デフォルトの名無しさん [2008/08/07(木) 23:34:38 ] >>616 ezplotとかだと関数じたいを引数として渡してるよね。 無名関数が可能なのかは知らない（たぶん出来ない）。 653 名前：デフォルトの名無しさん mailto:sage [2008/08/07(木) 23:50:55 ] help function_handle 654 名前：デフォルトの名無しさん mailto:sage [2008/08/08(金) 01:53:14 ] 無名関数は知ってるけど、 その引数の数を調べられないか？ってのが大事なんじゃろ。 で、matlabの引数って function hoge(a,b,c) と書いてあっても hoge(a,b,c,d) と呼べるっちゃよべるし、呼び出しの時点になるまでわからないのでは富田 655 名前：デフォルトの名無しさん mailto:sage [2008/08/08(金) 03:24:53 ] matlab.exe のオプションってどう見るの？ --help でも matlab 起動するし。 656 名前：デフォルトの名無しさん mailto:sage [2008/08/08(金) 06:02:13 ] >>654 >>652に答えてるのにお前が知ってたらなんなんだよ。。 だいたい呼び出しの時点の引数の数を呼び出す前に知る方法があるわけないし。 何言ってるのかわからんが。 >>616で書いてるのは定義にある引数の数だろ。だったらnargin。 657 名前：デフォルトの名無しさん mailto:sage [2008/08/15(金) 05:55:31 ] [X,Y,Z] = peaks(30); surfc(X,Y,Z) で出る図を png で保存すると、 真っ黒画像が作られるのですが、 みなさんはきちんと保存できたりしますか？ ぐぐってもだれも surf がおかしいとか言っていないようで。 658 名前：デフォルトの名無しさん mailto:sage [2008/08/15(金) 06:18:21 ] GUIでエディットテキストを30個作り、これにあらかじめ定義しておいた15*2サイズの行列Aの 各成分を一つずつ表示しようとしてます。edit1=>A(1,1), edit2=>A(1,2), edit3=>A(2,1)という 具合にです。 set(handles.edit1,'string',num2str(A(1,1))); などとSet文を30個書くのは効率が悪いんで、簡潔に書く方法を探しています。とりあえずevalを 使えばいいのかと思い、 for i=1:30 　if rem(i,2)==1 　　eval(['set(handles.edit' num2str(i) ',''string'','... 　　　'num2str(A(' num2str(ceil(i/2)) ',1)));']); 　else 　　eval(['set(handles.edit' num2str(i) ',''string'','... 　　　'num2str(A(' num2str(2/i) ',2)));']); 　end end としてみましたが、「evalのインデックスはpositive integerかlogicalでないといけない」みたいな エラーが出てしまいます。何が悪いんでしょうか。あるいは、eval以外にもっと良い方法があれば 教えてもらえると助かります。 659 名前：デフォルトの名無しさん mailto:sage [2008/08/15(金) 14:37:19 ] 別ディレクトリの関数を addpath せずに実行する方法ってないでしょうか？ その関数を実行して、中で一連の addpath をさせたいのです。 660 名前：デフォルトの名無しさん mailto:sage [2008/08/15(金) 14:41:36 ] >>658 とりあえず sprintf のほうが見やすそう。 661 名前：デフォルトの名無しさん mailto:sage [2008/08/15(金) 15:55:16 ] >>657 オプションを設定 662 名前：デフォルトの名無しさん mailto:sage [2008/08/15(金) 15:56:25 ] >>658 エラーの内容そのまま。 １行づつ、確認しろ。 663 名前：658 mailto:sage [2008/08/15(金) 16:53:11 ] レスくれた人達ありがとう。 うおお、下から３行目の num2str(2/i) は num2str(i/2) とすべきでしたね。evalの性質の理解不足 とか以前の問題でしたorz　あと、sprintfを提案してくれた人もありがとう。これからいろいろ試してみます。 664 名前：デフォルトの名無しさん mailto:sage [2008/08/16(土) 03:08:34 ] >>661 オプションって surf のオプション？どれ？ 665 名前：658 mailto:sage [2008/08/16(土) 05:21:46 ] 昨日はどうも。もう一つだけ質問させて下さい。 >>658とは逆のケースで、30個のエディットテキストにユーザーが入力した値を、15*2サイズの 行列Aとして格納することを考えています。 guideでエディットテキストを配置すると、対応するm-ファイルにCallback関数が30個自動的に 生成されるので、それぞれの下部に各行列成分への格納文を書けばいいわけですが、やはり非効率的 なので、もっと良い方法を探しています。 この場合、これら30個の自動生成された関数宣言文を削除し、evalなりfevalなりを使って、30個の 関数を変数を使ってまとめて定義するような文を手で打ってもいいんでしょうか？　確かこの手の、 「guideを経由しないでm-ファイルに直打ち」はNGのケースがある、と何かの本で見たことがある んですが。 あるいは、このような場合での常套手段が他にありましたら教えて下さい。 666 名前：デフォルトの名無しさん mailto:sage [2008/08/17(日) 12:10:20 ] >>659 スクリプトでいいならrunでできるけど 関数を呼ぶにははそれと組合せて function handleを使うとかかな。 それかベタだけど wd=cd(dir); func; cd(wd); とか。 667 名前：デフォルトの名無しさん mailto:sage [2008/08/17(日) 12:23:29 ] >>666 addpath をしたいだけなので、 スクリプトにして、run で呼べばよさそうですね。ありがとう。 668 名前：デフォルトの名無しさん [2008/08/17(日) 19:18:43 ] ランダムな行列の中から局所的な最大値とその座標を求めたいのですが、 何かいい方法がありますでしょうか？ max関数を２回使うことで、大域的な最大値は求まるのですが。。 局所的な最大値を求めるアイデアがどうしても出てきません。。 お願いします。 669 名前：デフォルトの名無しさん mailto:sage [2008/08/17(日) 23:01:34 ] >>668 大域的なmax値をNaNで置換して、複数回走らせる あるいは 元データを微分する 670 名前：デフォルトの名無しさん mailto:sage [2008/08/18(月) 07:09:50 ] >>668 局所的な最大値って極大値でいいのか? いずれにしても「ランダムな行列」がどういうものであるかと 何を「局所的な最大値」と定義するかが決まらないと 答えられないのでは? 671 名前：668 [2008/08/18(月) 11:29:38 ] レスありがとうございます。 具体的に書くべきでした。 例えば、 1 4 6 8 4 6 4 5 1 9 0 7 4 0 9 4 1 6 4 7 5 9 0 2 とあった場合、３列目の９と４列目の２つの９、６列目の７が最大値（＝極大値） と考えたときです。 お願いします。 672 名前：デフォルトの名無しさん mailto:sage [2008/08/18(月) 22:10:41 ] 座標が整数値でない画像データを可視化したいんですが、どうやったらいいですか？ データの形式としては(x, y)が座標、(R, G, B)が与えたい色として x y R G B 0.1 -1.1 10 20 30 0.5 0.3 112 30 11 -2.4 3.2 3 15 205 のようにn*5（nは行数）のベクトルデータとなっています。 これを穴があってもいいから通常の画像のように可視化したいんです。 点ごとにプロットする色を変えたらいいのかと思ったんですがうまくいかなくて 673 名前：デフォルトの名無しさん mailto:sage [2008/08/18(月) 22:58:01 ] >>672 座標を例えば 1 - 128 に正規化して、I = zeros(128); 後は for ループで x = data(n,1); y = data(n,2); I(x,y,1) = data(n,3); I(x,y,2) = data(n,4); I(x,y,3) = data(n,5); とでも格納しなおして imshow するとかそういう単純な話だろ？ >>671 matlab central れ。 674 名前：デフォルトの名無しさん [2008/08/19(火) 04:39:11 ] ある関数Ｆの最適化問題を解きたいんだけど、Ｆがゼロになった瞬間に 最適化を止めて次のことをさせたいような場合、どうしたらいいですか？ 最適化にはＮｅｌｄｅｒ−Ｍｅａｄを使いたいので、Ｆｍｉｎｕｎｃより Ｆｍｉｎｓｅａｒｃｈを使いたいです。 Ｆの値は定義上０より小さくなることはないので、一度０がでてからさらに パラメター空間でサーチするのは単なる時間の無題です。 Ｆｍｉｎｓｅａｒｃｈのコードがわかれば外にＩＦ節をつけれてＦの値が ０になったところで止めさせられるものの、どうしていいかわかりません。 675 名前：デフォルトの名無しさん mailto:sage [2008/08/19(火) 07:27:11 ] A = [1 1 2 1 3 4 5 1 1 2 3 6] みたいなベクトルがあって、[1 2] のある位置を知るのに、 for ループを使わずにすむ方法ってありますか？ matlab ならありそうな気がするのですが B = [0 1 0 0 0 0 0 0 1 0 0] が返ってほしい。 676 名前：672 mailto:sage [2008/08/19(火) 20:22:10 ] >>673 即レスありがとうございます ちなみに正規化ってどうやったらいいんでしょうか xの列でソート→代入、yの列でソート→代入じゃ効率悪すぎですよね 手早くやれる方法はありますか？ 677 名前：デフォルトの名無しさん [2008/08/19(火) 21:11:55 ] １から１００の整数をランダムに並び替えるプログラムを教えてください！ 678 名前：デフォルトの名無しさん mailto:sage [2008/08/19(火) 21:27:19 ] >>677 宿題は自分でやれ 679 名前：デフォルトの名無しさん mailto:sage [2008/08/20(水) 00:42:24 ] >>677 randperm >>676 ソート？なぜにソートが・・・・ low = [1 1]; high = [128 128]; [N, d] = size(data); % d == 5 mini = min(data(:,1:2), [], 1); maxi = max(data(:,1:2), [], 1); data(:, 1:2) = data(:,1:2) - repmat(mini, N, 1); % move min to 0 data(:, 1:2) = data(:,1:2) ./ repmat(maxi-mini, N, 1); % [0, 1] data(:, 1:2) = data(:, 1:2) .* repmat(high-low, N, 1); % [0, high-low] data(:, 1:2) = data(:, 1:2) + repmat(low, N, 1); % [low, high] data(:, 1:2) = fix(data(:, 1:2)); % int 680 名前：デフォルトの名無しさん mailto:sage [2008/08/20(水) 02:50:49 ] すみません、どなたか>>665が分かる方いらっしゃいますか？ 本質的には、いかにしてCallBack関数をまとめて設定するか、ということだと思います。 681 名前：668 [2008/08/20(水) 14:53:04 ] >>673 matlab centralにあるということでしょうか？ 探してみましたが、わかりませんでした。。 もしよければどれか教えてもらえないでしょうか？ お願いします。 682 名前：デフォルトの名無しさん mailto:sage [2008/08/20(水) 15:05:44 ] えー、なんで？極大値は英語にしたら extreme か extrema だよ？ extrema.m 683 名前：668 [2008/08/20(水) 23:06:13 ] >>682 完全に見落としてました。。 これでいけそうです。 ありがとうございました！ 684 名前：デフォルトの名無しさん [2008/08/21(木) 10:04:17 ] for文を使わないことはできますか？ ２次関数(ax^2+bx+c)を入力にして、総あたりでa,b,cのパラメータを推定する、というようなことをやっています。 入力をf(x)として、-5<x<5の範囲を見てるのですが、コードは error=inf; for a=-10:1:10 for b=-10:1:10 for c=-10:1:10 for t=-5:0.1:5 e=(f(t)-(a*t*t+b*t+c))^2; if (error>e); error=e; end end end end end というような感じになってます。これをfor文を使わないでコロンを使った形式にできるものなのでしょうか？ よろしくお願いします。 685 名前：デフォルトの名無しさん [2008/08/21(木) 10:07:19 ] すみません、間違えてタブを使ってしまったので投稿しなおします ****** for文を使わないことはできますか？ ２次関数(ax^2+bx+c)を入力にして、総あたりでa,b,cのパラメータを推定する、というようなことをやっています。 入力をf(x)として、-5<x<5の範囲を見てるのですが、コードは error=inf; for a=-10:1:10 for b=-10:1:10 for c=-10:1:10 for t=-5:0.1:5 e=(f(t)-(a*t*t+b*t+c))^2; if (error>e); error=e; end end end end end というような感じになってます。これをfor文を使わないでコロンを使った形式にできるものなのでしょうか？ よろしくお願いします。 686 名前：デフォルトの名無しさん [2008/08/21(木) 10:10:21 ] スペースもダメみたいですね…みにくくてすみません。 よろしくお願いします。 687 名前：デフォルトの名無しさん mailto:sage [2008/08/21(木) 10:16:19 ] repmat を駆使して４次元配列にまでもっていけばきっとできる。 で、あとは min で argmin うけとる。 688 名前：684 [2008/08/21(木) 16:35:55 ] >>687 ご指摘いただいた方法でできました。 ありがとうございます。 689 名前：デフォルトの名無しさん [2008/08/21(木) 21:09:03 ] >>684って行列の中身が一緒じゃないんでしょ？ だったらrepmatじゃだめじゃね？ ってか俺も今 for a=1:10 for b=1:5 kotae(a,b)=a+b end end って感じのやってるけど、これfor文じゃなくせるの？ 690 名前：デフォルトの名無しさん [2008/08/22(金) 00:13:06 ] matlab centralで拾ったファイルを動かそうと思ったんですが、 「500の最大の再帰呼出の制限達しました。制限を変更するには、・・・」 となって動きませんでした。 制限を2000あたりまで変更しても同じ結果で。3000くらいになるとMatlabが 強制終了してしまいました。 これ何か対処法ないでしょうか。。？？ 691 名前：デフォルトの名無しさん mailto:sage [2008/08/22(金) 02:39:48 ] mファイルの作成者に聞くべきじゃない？ つか、少しは自分でk 692 名前：デフォルトの名無しさん mailto:sage [2008/08/22(金) 02:56:10 ] >>689 なくせう。for 文使うとゲロのろいぞ matlab。 693 名前：デフォルトの名無しさん [2008/08/22(金) 09:31:00 ] コロンでできるの？ 694 名前：デフォルトの名無しさん mailto:sage [2008/08/22(金) 10:26:09 ] できる。 kotaeって10x5のマトリックスになるだろ？ 695 名前：684 [2008/08/22(金) 17:33:09 ] >>689 meshgridというのを使うとてっとり早いです。 行列にして掛け算する、ということです。 687さんがおっしゃった通り、repmatを使ってもできます。 というか、repmatを使ってmeshgridが実装されてるようなイメージです。 696 名前：MATLAB初心者 [2008/08/23(土) 12:06:04 ] はじめまして。 ＭＡＴＬＡＢをはじめたばかりの初心者なのですが、連番をつけるところで悩んでおります。 文字列の配列で0001、0002、.....　、9999のように整数で任意の桁数にしたいのですが うまい方法はありませんでしょうか？ エクセルでやっていたときは、上記の場合であれば”000”という文字列と整数の値を結合し 左から4文字を切り出して使用しておりました。 ＭＡＴＬＡＢでもエクセルでやっていたような力技でやるしかないのでしょうか？ もしご存知の方がおりましたら、ご教示いただけませんでしょうか？ よろしくお願いいたします。 697 名前：MATLAB初心者 [2008/08/23(土) 12:12:11 ] 申し訳ございません。 上記のレスで左から切り出していたというのは誤りで、右から切り出しておりました。 698 名前：デフォルトの名無しさん mailto:sage [2008/08/23(土) 14:51:58 ] >>697 ttp://dl.cybernet.co.jp/matlab/support/manual/r13/toolbox/matlab/ref/?/matlab/support/manual/r13/toolbox/matlab/ref/sprintf.shtml 699 名前：MATLAB初心者 [2008/08/23(土) 15:46:04 ] >>698 ありがとうございます。 大変申し訳ございません。私には、このページではsprintfでたとえば”1”を”00001”と定義する方法が理解できませんでした。 整数を5桁で表して、それを文字型にnum2strか何かで文字列に変換するようなことができるということでしょうか？ 私の思い付きではstrcatをつかって696に書いたような力技で書く以外に方法がわかりません。 お手数をおかけしますが、ご教示いただいたページのどれを使用すればよいか教えていただけませんでしょうか？ 700 名前：デフォルトの名無しさん mailto:sage [2008/08/23(土) 16:25:44 ] >>699 printf フォーマット でぐぐればいろいろ解説があるよ。 Ｃ言語が元だけど自由度がかなり高いから覚えとくとよい ↓とりあえず解決策 str=sprintf('%05d',i) 701 名前：MATLAB初心者 [2008/08/23(土) 16:35:09 ] >>700 取り合えずというか、そのものずばり私の欲しかった回答を書いていただいて、感謝の限りです。 お忙しい中、どうもありがとうございました。 これからも精進いたします。 702 名前：675 mailto:sage [2008/08/26(火) 03:23:35 ] A = [1 1 2 1 3 4 5 1 1 2 3 6] みたいなベクトルがあって、[1 2] のある位置を知るのに、 for ループを使わずにすむ方法ってありますか？ matlab ならありそうな気がするのですが B = [0 1 0 0 0 0 0 0 1 0 0] が返ってほしい。 703 名前：デフォルトの名無しさん [2008/08/27(水) 10:05:55 ] findじゃだめ？ 704 名前：デフォルトの名無しさん mailto:sage [2008/08/28(木) 08:09:50 ] 俺の力じゃ5行が限界だった x=find(A==1); y=find(A==2)-1; z=intersect(x,y); B=zeros(1,length(A)); B(z)=1; もっと簡単なのきぼん 705 名前：デフォルトの名無しさん mailto:sage [2008/08/28(木) 22:26:43 ] 面白そうなので参加 ((A(1:length(A)-1)+A(2:length(A))).*(A(1:length(A)-1)==1))==3 正の値で且つ後ろ一つは無視する前提 一般性があるかどうかは知らない 706 名前：デフォルトの名無しさん mailto:sage [2008/08/28(木) 22:38:40 ] size(A)も使ってみよう 707 名前：デフォルトの名無しさん mailto:sage [2008/08/29(金) 01:15:31 ] A(1:end-1)==1&A(2:end)==2 レベル低いなあ。 こんなのすらすっと出てこないってことは よっぽど効率悪いコード書いてるんだろうなお前ら。 708 名前：デフォルトの名無しさん mailto:sage [2008/08/29(金) 01:40:14 ] 煽るのは構わんがそれ全部ゼロになるぞ 709 名前：デフォルトの名無しさん mailto:sage [2008/08/29(金) 01:46:06 ] って１ズラしてるか スマン 710 名前：デフォルトの名無しさん mailto:sage [2008/08/29(金) 09:07:52 ] MATLABの関数自体のコードを読んでみると 勉強になるんだけどな。 711 名前：デフォルトの名無しさん [2008/09/04(木) 13:25:51 ] 質問です。 MATLAB上で画像処理をする際に 画像のある画素(座標)を保存して さらにその保存した座標のｘ方向成分だけ読み出すには どういったコマンドを入力すればいいんですか？ 712 名前：デフォルトの名無しさん mailto:sage [2008/09/04(木) 13:42:08 ] >>711 まず入門サイトでも読んだ方がはやくね？ってかはやい 713 名前：デフォルトの名無しさん [2008/09/04(木) 16:05:56 ] 2008からオブジェクト指向が充実したみたいで、 それについて詳しく解説してある本を知っていたら 誰か紹介していただけないでしょうか。 よろしくお願いします 714 名前：デフォルトの名無しさん [2008/09/08(月) 12:56:09 ] 誰か、MATLAB(2007b)で３次元の楕円球体の描き方教えてください。 715 名前：デフォルトの名無しさん [2008/09/08(月) 18:24:34 ] 定義したい変数名を、文字列を使って定義することは可能なのでしょうか？ 例えば、 Data1 = 1; Data2 = 2; ・ ・ ・ Data30 = 30; と30個の変数を定義したいとき、for文を使って定義することは可能でしょうか？ 716 名前：デフォルトの名無しさん mailto:sage [2008/09/08(月) 23:43:57 ] >>714 キーワードを英語でlookforするくらいはしようよ。 >>715 help eval 717 名前：デフォルトの名無しさん [2008/09/09(火) 11:36:43 ] はじめまして、今matlabの勉強をしている者です。 matlabのエディットウィンドウで、文字の挿入ってどうやるんですか？ 現状では入力してある行に文字を挿入すると後ろの文字が消えてしまい、 全文入力し直しになってしまうのです、これをwordなどみたいに 文字を挿入したら、文が自動で後ろに移動してくれるようにできないですか？ また、同様にスペースキーも文字間のスペースをあけてくれず 文字を消してしまうのですが？どうしたらいいですか？ OSは vinelinux4 matalab v7.0.1 SP1です 718 名前：デフォルトの名無しさん mailto:sage [2008/09/09(火) 22:12:45 ] MATLABっていうか キーボードのinsert 719 名前：デフォルトの名無しさん [2008/09/10(水) 13:58:00 ] GUIを作成して、 ボタンを押すとシミュレーションが実行するようにしたんですが、 GUIなしでコマンドから直接実行した場合だと１６０秒くらいで終わるのに、 なぜかGUIからだと２６０秒もかかるんですが、 なんでこんなにかかると思いますか？ やってることはまったく変わらないですが・・・。 720 名前：133.43.185.64 mailto:sage [2008/09/10(水) 21:20:24 ] >>719 全然見当がつかん。 プロファイラを実行した状態で GUI 経由でシミュレーションを行って、 時間のかかっている箇所を割り出してから相談してもらった方がよいと思われ。 721 名前：デフォルトの名無しさん [2008/09/11(木) 00:44:20 ] >>720 アドバイスありがとう。 実はこういう経験あまりしてなかったんで、 プロファイラとか使ったことなかったんです。 ちょっと調べて割り出してみます。 722 名前：初心者 [2008/09/11(木) 01:43:15 ] 非常に低レベルなことで申し訳ないのですが、 Cでいうポインタのような概念はMatlabにはないのでしょうか？ 723 名前：デフォルトの名無しさん mailto:sage [2008/09/11(木) 01:54:01 ] >>722 ないはずだよ wikipediaにも書いてある まあ、google codeとかmatlab centralで検索すれば 自作のポインタ関数出て来ても可笑しくないけど 724 名前：デフォルトの名無しさん mailto:sage [2008/09/11(木) 02:15:40 ] >>722 デフォはポインタ 渡された変数を編集しようとした時点でコピーが始まる 725 名前：デフォルトの名無しさん mailto:sage [2008/09/11(木) 06:15:09 ] >>722 何をしたいのかをはっきりさせて それをMATLABではどうやるのか、と考える方がいいよ。 他の言語のやり方をそのまま持って来ようとしても そうは行かない場合は多い。MATLABに限らず。 726 名前：デフォルトの名無しさん [2008/09/14(日) 02:33:41 ] 質問です。 例えば 3,6,4,5,2 4,15,7,5,3 6,4,4,7,3 のような行列があってノイズで15だけ値が大きい場合には、 それを周りの数字に近い値に矯正（たしか平滑？）するにはどうすればいいですか？ 自分が知ってるやり方だと、ノイズも含め全ての成分が平均化されてしまいます。 どうにかノイズ以外はそのままの値で、ノイズ部分だけを周りの値に馴染ませたいんです。 727 名前：デフォルトの名無しさん mailto:sage [2008/09/14(日) 02:59:50 ] 10はノイズ? 13はノイズ? 15を適当に5とかにしたら「馴染んだ」ことになるの? 自分が何をやりたいかをまずは具体的にする。 君以外にそれがわかる人はここにはいないよ。 728 名前：デフォルトの名無しさん mailto:sage [2008/09/14(日) 07:57:32 ] 実験データの改竄捏造でもしたいのかなｗ 729 名前：デフォルトの名無しさん mailto:sage [2008/09/14(日) 12:03:36 ] >>726 if ノイズ    平滑 else    なにもしない 730 名前：デフォルトの名無しさん mailto:sage [2008/09/15(月) 07:18:58 ] 　ああ　 731 名前：デフォルトの名無しさん mailto:sage [2008/09/15(月) 13:38:27 ] >>726 新しいやり方を考えればいいだけ。 732 名前：デフォルトの名無しさん [2008/09/17(水) 23:20:14 ] 質問です。 Ｃのプログラム等に記載するための離散時間系のデジタル１次ローパスフィルタを作るためにＭＡＴＬＡＢでフィルタの係数を求めようとしています。 そのためのプログラムとして、サンプリングインターバル＝６ｍｓ、フィルタ周波数＝１０Ｈｚとして、 　　　omega = 2 * pi * 10 T = 1 / omega sys = tf(1, [T, 1]) [num, den] = tfdata(sys) [numd, dend] = c2dm(num, den, 0.006) として求めると、y[n] = 0.6859 * y[n-1] + (1 - 0.6859) * x[n]という式が求められます。 一方で、連続時間系の1次ローパスフィルタの伝達関数（サンプリングインターバル＝6ms、フィルタ周波数＝10Hz）　G(s) = 1 / (1+Ts) を手計算でz変換を行い、離散系のローパスフィルタの式を求めると、（T = 1 / omega、omega = 2 * pi * 10） 　　　y[n] = 1 / (1+omega*0.006) * y[n-1] + (1 - 1 / (1+omega*0.006)) * x[n] となるため、値を代入すると、 　　　y[n] = 0.7263 * y[n-1] + (1 – 0.7263) * x[n] となり、先ほどＭＡＴＬＡＢ関数から求めた式と乖離が出てきてしまいました。 そこで、有識者の方々にお聞きしたいのは 　　　１．そもそも上記のデジタル１次ローパスフィルタの作り方は間違っているのか？ 　　　２．c2dm コマンドの中身がよくわからないが、これは計算誤差が出るようなものなのか？ 　　　３．離散系に変換したシステムのボード線図を見て特性を確認したいのだが、その場合 　　　　　も連続系と同じく 　　　　　　　　　　bode(numd, dend) 　　　　　でよいのか？ 長文となってしまい、見にくくて申し訳ありませんがどなたか有識者の方がいらっしゃればご教授ください。 733 名前：デフォルトの名無しさん mailto:sage [2008/09/18(木) 14:54:34 ] 簡単な計算ミスのような気もするのですが、、、 734 名前：デフォルトの名無しさん [2008/09/18(木) 22:11:28 ] >>733 と言いますと・・・？ 735 名前：デフォルトの名無しさん mailto:sage [2008/09/18(木) 22:13:54 ] 周波数とインターバル（周期）とで混乱していない？ 間違いを探したいのなら、 他の数字の組み合わせにしてみて、計算してみるのはどうでしょう？ 736 名前：デフォルトの名無しさん [2008/09/18(木) 23:19:11 ] >>735 すいません、だいぶ他の数字でもやってみたのですが・・・ 基本的な結果としては、低周波（0.01Hz、1Hz）では値の乖離はほとんど見られないのですが、高周波になればなるほどどんどん値が離れていくということになっています・・・ なので、アナログからデジタルに変換するときの計算誤差なのではと思った次第なのですが・・・ 　　　　　　　　　　　MATLAB計算値　　　　　　　　　手計算値 　　　0.01Hz　　　　　　0.9996　　　　　　　　　　　　0.9996 0.5Hz 0.9813　　　　　　　　　　　 0.9815 5Hz 0.8282 0.8414 10Hz 0.6859 0.7260 どなたかご教授いただけると助かります 737 名前：デフォルトの名無しさん mailto:sage [2008/09/18(木) 23:34:44 ] >>736 サンプリングインターバルも変えてみたらどう？ 738 名前：デフォルトの名無しさん [2008/09/19(金) 00:13:17 ] 変えてみた結果は以下の通りになりました。 　　　　　　　　　　　　　MATLAB計算値　　　　　　　　手計算値 　0.5Hz　　12msec　　　　　　0.9630　　　　　　　　　　0.9637 2Hz 24msec 0.7396 0.7683 0.01Hz 480msec 0.9703 0.9707 10Hz 480msec 0.0000 0.0032 という感じで、サンプリングインターバルが大きくなる、もしくは周波数が大きくなると乖離がどんどん大きくなるといった傾向が見られます。 そこで離散系に直したときのボード線図を見て、連続系のそれとの特性の違いを見たいのですが、離散系に直したシステム 　　　　　[numd, dend] = c2dm(num, den, 0.006) 　　　　　bode(numd, dend) という見方でよろしいのでしょうか？明らかに異なった波形が出てくるので、離散系の場合はボード線図を違う見方で見なければいけないのかなと予想しているのですが。。。 見方がわかりませんので、もしご存知でしたら教えてください。 739 名前：デフォルトの名無しさん mailto:sage [2008/09/19(金) 00:33:01 ] ttp://www.ikko.k.hosei.ac.jp/joomla1506/matlab-topmenu-33/17-matlab20070307.html 740 名前：デフォルトの名無しさん mailto:sage [2008/09/23(火) 00:35:38 ] matlabをはじめたばかりなんですが質問です。 stem3でつくったグラフを等高線のグラフにそのまま変えたいのですが、 どなたかわかる方いらっしゃいましたらお願いします。 741 名前：デフォルトの名無しさん mailto:sage [2008/09/23(火) 02:50:24 ] ステムを輪切りにして上から見て、意味あるの？ 742 名前：デフォルトの名無しさん mailto:sage [2008/09/23(火) 11:49:33 ] >741 値がどのように分布しているか等高線の方がわかりやすいと思い 変えたいのですが。 743 名前：デフォルトの名無しさん mailto:sage [2008/09/23(火) 15:11:34 ] >>742 matlab 等高線 でぐぐれ 744 名前：デフォルトの名無しさん [2008/09/30(火) 15:46:00 ] GUIのプログラムについての質問です ファイルguimain.m内で handles.a=1; guidata(hObject, handles); とhandles.aを設定したとき 別のファイルgui2.m内で handles.aを編集する方法はあるのでしょうか function gui2(hObject,handles) handles.a=2; guidata(hObject, handles); ではできませんよね? 745 名前：デフォルトの名無しさん [2008/09/30(火) 16:32:46 ] 効率はわかりませんが function out = gui2(handles) handles.a=2; out = handles; として、呼び出し側でguidata(hObject, handles); することで解決しました。 746 名前：t [2008/10/07(火) 19:46:19 ] ベルヌーイの螺旋を描ける方いらっしゃいませんか？ 困っています。 747 名前：デフォルトの名無しさん mailto:sage [2008/10/15(水) 20:23:01 ] 描けるよ 748 名前：デフォルトの名無しさん mailto:sage [2008/10/18(土) 10:33:41 ] >>745 構造体使う必要がないならglobal宣言すれば簡単かと 749 名前：デフォルトの名無しさん [2008/10/19(日) 22:46:09 ] z = f(x, y) を、3次元プロットしたいのですが、 どうやったらよいのでしょうか？ 試しに ---------------------- z = [ 0 0.1 0.3 0.5 0.3 0; 0 0.3 0.5 0.3 0.3 0; 0 0.5 1.0 0.7 0.5 0; 0 0.3 0.5 0.3 0.3 0; 0 0.1 0.3 0.5 0.3 0 ]; [x, y] = size(z); surfc(x, y, z) ---------------------- としたのですが、うまくいきませんでした。 おそらく行列のプロットを理解していない と思うのですが・・・。 非常に低レベルな質問で申し訳ありませんが、 よろしくお願いします。 750 名前：デフォルトの名無しさん mailto:sage [2008/10/19(日) 23:52:00 ] mesh(z)でどうだ 751 名前：749 [2008/10/20(月) 00:01:45 ] >>750 様 早速のご回答、ありがとうございます！ 無事、3D plot をすることができました。 ありがとうございます！ 752 名前：デフォルトの名無しさん mailto:sage [2008/10/23(木) 17:54:34 ] MATLABで式の簡単にする変形を一ステップずつやって どうやって、変形して最終の形にするのか 見る方法ありませすか？ 753 名前：デフォルトの名無しさん mailto:sage [2008/10/23(木) 23:45:04 ] 意味がわからない。 mathematica？ 754 名前：デフォルトの名無しさん mailto:sage [2008/10/23(木) 23:47:52 ] シンプリファイとかやるときの途中の式がしりたい。 シンボリックマスツールズの 755 名前：いかわ [2008/10/24(金) 15:35:48 ] '"2008-09-23 02:30:00.1"'　というフォーマットで書かれたストリングのセルから時間データを抽出するよいアイディアをどなたかご存知でないでしょうか？ もしも、''がなく、"2008-09-23 02:30:00.1"　というフォーマットならば、sscanf(%d-%d-%d %d:%d:%d);　によって時間データを抽出できるのですが。。 756 名前：デフォルトの名無しさん mailto:sage [2008/10/24(金) 20:06:56 ] 日付の区切りがスラッシュにかわったりしただけでも どうすればいいかわからなくなるのかなこういう人は・・。 757 名前：デフォルトの名無しさん mailto:sage [2008/10/25(土) 04:01:54 ] だれか Matlab compiler で matllab スクリプトコンパイルして、 C で使ったことある人いない？ ttp://www.mathworks.com/access/helpdesk/help/toolbox/compiler/index.html?/access/helpdesk/help/toolbox/compiler/f2-972343.html これの通りやってみて、C をコンパイルできるようにまでは持って いったんだけど、実行すると libmatrixInitialization: unhandled exception と出る。なにが足りないんでしょうかね？ 758 名前：デフォルトの名無しさん mailto:sage [2008/10/25(土) 05:49:28 ] 愛 759 名前：デフォルトの名無しさん mailto:sage [2008/10/25(土) 10:11:57 ] >>757 回答にはなってないかもしれんけど、 おいらは、加えてBuilder for .NETを使ってMファイルを.NETコンポーネント化 してから、.NET系言語から、そのコンポーネント呼び出して使ってるよ。 760 名前：デフォルトの名無しさん mailto:sage [2008/10/27(月) 15:42:45 ] A=zeros(100000,100000); だけでメモリ不足になります 余分な変数消去や仮想メモリ等をいじってみましたが効果はありませんでした これを解消するには結局メモリ増設しかないのでしょうか？ 761 名前：デフォルトの名無しさん mailto:sage [2008/10/27(月) 16:38:49 ] >>760 メモリ増設しても無理じゃね? 100000*100000って、数値だけで10Gになるからアクセスすることさえ厄介になってくると思う。 メモリ量的には更に数倍になるわけで、オンメモリで処理するのはどうせ無理だし。 762 名前：760 mailto:sage [2008/10/27(月) 18:09:13 ] >>761さん ご回答ありがとうございます 低レベルな質問ですみませんでした このような行列はMATLABではそもそも扱えない（普通は扱わない）と考えてよいのでしょうか？ 株価のシミュレーションでforを使うより行列を使った方が早いため あのような行列が処理できると助かったのですが 763 名前：デフォルトの名無しさん mailto:sage [2008/10/28(火) 01:30:51 ] 「普通のＰＣ」では扱えない　というのが正解かと。 １０Ｇのメモリっていうといわゆるスーパーコンピューターの領域 764 名前：デフォルトの名無しさん [2008/10/30(木) 04:26:55 ] たとえば*を増やしていって、 処理の進行状況を表示したいのですが、 改行を抑制する方法はありませんか？ 765 名前：デフォルトの名無しさん mailto:sage [2008/10/30(木) 21:20:25 ] 質問です 複素数の実部虚部それぞれ七桁で表示させたいのですが zに複素数が入っていて sprintf('%7f',z) とすると実部だけが七桁で表示されてしまいます どうしたらよいでしょうか？ 766 名前：デフォルトの名無しさん mailto:sage [2008/10/30(木) 22:41:17 ] 関数realとimagを使用して実部と虚部を取り出せばいいじゃない 767 名前：デフォルトの名無しさん [2008/10/30(木) 23:07:49 ] サポート受けられないということは、ここの住人はコピーや労か 768 名前：デフォルトの名無しさん [2008/10/31(金) 00:06:48 ] >>764 なんで今どうやってるのか書かないの? 769 名前：デフォルトの名無しさん mailto:sage [2008/10/31(金) 14:35:01 ] >>766 ありがとうございます sprintf('%7f %+7f',real(z) imag(z)i) のような感じでしょうか？ 今試せる環境ではないのでiとかがどうなるのかよく分かりませんが・・・・・ 770 名前：デフォルトの名無しさん mailto:sage [2008/10/31(金) 19:10:46 ] 教えてください。 １辺の長さが１の正三角形を上下に並べて出来る平行四辺形をプロットして、 それに含まれるかどうかの判定をしたいのですが、判定条件はどうしたらいいでしょうか。 N個の点を行列で与えて、図形に入っていれば赤、入っていなければ黒という風に 表示させたいのです。 よろしくお願いします 771 名前：デフォルトの名無しさん mailto:sage [2008/10/31(金) 19:25:24 ] >>769 +の位置が不適切だけどそれを直せばイケル 772 名前：デフォルトの名無しさん mailto:sage [2008/11/01(土) 17:39:35 ] >>771 すいませんまだ試せていないんですがこうすると実数値の場合に 3+0i みたいな表示になっちゃいませんかね？ もっと簡単にできる方法ないですかね 773 名前：デフォルトの名無しさん [2008/11/07(金) 14:26:53 ] 質問があります 例えばデータシートが1000枚あって、各データシートからX,Y,Zの値を取り出して、 座標(X,Y,Z)をplot3でプロットします。これがこのシートの3次元での点（位置）になりますよね。 それを1000枚分やると3次元空間に点が1000点でき、動きの軌道のように繋がって表示されます（そういうデータなので）。 問題は次なのですが、数十点ほど動きの軌道からはみ出してプロットされた点（多分計測時のノイズ）が 存在します。なのでプロット点全体を平滑化したいのですが、三次元プロットした点をフィルタするのは 可能なんでしょうか。 よろしくお願いします 774 名前：デフォルトの名無しさん [2008/11/07(金) 21:04:04 ] >>773 成分ごとにやれば？ 775 名前：デフォルトの名無しさん mailto:sage [2008/11/08(土) 07:04:44 ] 結局の所作業としてはプロットした点をフィルタせずに フィルタしてからプロットするんでしょ? 776 名前：デフォルトの名無しさん mailto:sage [2008/11/08(土) 09:03:51 ] Matlab C ライブラリの関数 mlfMatrixPrint　とか mlfLength とかってどこに行ったの？ 昔は matlab.h ってのがあったみたいなんだが、 R2007b だけど、ないんだよね。 matlab 内のビルトイン関数を C から直接呼び出せなくなったの？ 777 名前：デフォルトの名無しさん mailto:sage [2008/11/09(日) 14:40:02 ] 私は>>675じゃないんですけど >>675の解決策ってありますか? 778 名前：デフォルトの名無しさん mailto:sage [2008/11/10(月) 07:24:17 ] スレを675以降もネットリと見てみよう 779 名前：デフォルトの名無しさん [2008/11/10(月) 20:04:54 ] こんばんは。最近MATLABで動画像処理を始めたものですが、 250*350でフレーム数200のAVIの動画をRGBで格納するために T=zeros(250,350,3,200) としてみたのですが、メモリが足りませんと出てしまいました。 解決策がわかりません ちなみにPCはXP　Core(TM)2 Duo CPU　メモリ2Gを使っています。 よろしくお願いします 780 名前：デフォルトの名無しさん mailto:sage [2008/11/10(月) 22:36:17 ] [A;B]という行列(n行m列)から [A;A;B;B]という行列(n^2行m列)をfor文を 使わずに 作りたいのですが可能でしょうか。 基礎的な質問でしょうがよろしくお願いします。 781 名前：デフォルトの名無しさん mailto:sage [2008/11/10(月) 23:00:26 ] >>780 [A;A;B;B]だったらn*2行m列だと思うんですが。 782 名前：780 mailto:sage [2008/11/10(月) 23:08:49 ] 例が分かりづらくて済みません。別の例ですが[A;B;C]という行列だったら [A;A;A;B;B;B;C;C;C]という行列を作りたいのです。 783 名前：デフォルトの名無しさん mailto:sage [2008/11/10(月) 23:58:44 ] >>779 メモリが足りないのものはしょうがないのである。 本当に全部格納しなければいけないのか考えてみるのである 784 名前：デフォルトの名無しさん mailto:sage [2008/11/11(火) 00:01:09 ] >>782 つまり *3 でいいんだね？ X = [A;B;C] からA, B, C の分解は終わっているものとする。 Y = [repmat(A, 3, 1); repmat(B, 3, 1); repmat(C, 3, 1)] 785 名前：780 mailto:sage [2008/11/11(火) 00:57:46 ] >>784 ありがとうございます。確かに直接記述すればそうなるんですが、 行列の行数が任意の数nであってもそのような行列が作れるか 、つまり各行ベクトルをn個並べてn^2行の行列を作れるかということが 分からないんです。 786 名前：デフォルトの名無しさん mailto:sage [2008/11/11(火) 06:44:37 ] matlab Mex で X(1,:) ってどうやんの？ 787 名前：デフォルトの名無しさん mailto:sage [2008/11/11(火) 10:22:08 ] >>785 repmat(x, 1, length(x)) 788 名前：780 mailto:sage [2008/11/11(火) 13:01:35 ] >>787 おかげで解決しましたーありがとうございます。 789 名前：デフォルトの名無しさん [2008/11/11(火) 18:11:41 ] [id,check] = textread('01.txt','%s%s'); と読み込んだ後 id 　とすると中身があるのは確認できるのですが fprintf('filename','%s',id); というふうにして書き込もうとするとFunction is not defined for 'cell' inputs. というエラーがでてidをtxtに書き込むことができません。 textreadで読み込んだidをtxtに書き込むにはどのようにすればよいでしょうか？ よろしくお願いします。 790 名前：デフォルトの名無しさん mailto:sage [2008/11/11(火) 21:31:48 ] >>789 www.excite.co.jp/world/ 791 名前：789 [2008/11/12(水) 16:00:23 ] >>789エラーの内容は分かったんですけれど セルをtxtに書き込めるようにする方法が調べてみてもよく分かりません もしよかったら教えてもらえませんか？ よろしくお願いします。 792 名前：シン [2008/11/14(金) 12:57:46 ] MATRABのプログラムで、ｘ１＝０、ｘ２＝１．２ ｘ３＝２．１ ｘ４＝３．５ ｙ１＝０．３ ｙ２＝１．２ ｙ３＝２．４ ｙ４＝２．７ の４点のラグランジュ補間をとりグラフを出力したいのですが教えてください 。できればfor文をつかった形式でお願いします。あと範囲は０＜ｘ＜１０です。 793 名前：デフォルトの名無しさん mailto:sage [2008/11/14(金) 18:27:37 ] レベル低過ぎ 説明書読め 794 名前：デフォルトの名無しさん mailto:sage [2008/11/14(金) 21:16:11 ] MATLABわかる人 pc11.2ch.net/test/read.cgi/tech/1226576323/l50 スレからのマルチポスト 795 名前：デフォルトの名無しさん [2008/11/15(土) 17:26:18 ] MATLABのプログラムで、xy平面上の点P(i)=(x(i),y(i)), i=1,2,3,4,5について ５個の点P(i)を適当に設定し、これらの点を順に線分で結んでできる多角形の図を作れ。ただし各点P(i)は*印で示すこと 最後の点P(5)を最初の点P(1)に結んでみよ。 のプログラムを作りたいんですけど、誰かわかる方教えてください。 796 名前：デフォルトの名無しさん mailto:sage [2008/11/16(日) 12:12:10 ] >>792 plot([1:10], lagrange(x, y)); >>795 line(x, y, '*'); 宿題は宿題スレに頼めよ 797 名前：デフォルトの名無しさん [2008/11/16(日) 13:49:26 ] 宿題スレのURL教えて！ 798 名前：770 mailto:sage [2008/11/17(月) 07:40:28 ] ｼｮﾎﾞｰﾝ(´・ω・｀) 799 名前：デフォルトの名無しさん mailto:sage [2008/11/17(月) 11:33:58 ] >>770 まはやMATLABですらねぇな 三角形の中に点があるかどうかは時計回り（または反時計回り）に 各辺を回ったときに常に辺から左側（または右側）に点があるかどうかでわかる N角形〜閉曲線の場合でも一緒 800 名前：デフォルトの名無しさん [2008/11/18(火) 19:50:56 ] MATLABで2の5乗根の値を出したいのですが何かいい関数がありませんか？ よろしくお願いします 801 名前：デフォルトの名無しさん mailto:sage [2008/11/18(火) 22:09:27 ] 2^(1/5) 802 名前：800 mailto:sage [2008/11/19(水) 15:32:16 ] ありがとうございます。 803 名前：デフォルトの名無しさん [2008/11/20(木) 14:09:32 ] 質問です。 MATLAB5.0でフォルダ内のファイル全てに同じ処理を施すようにしたいのですが、 どのようにしたらいいでしょうか？ サンプルプログラムを書いていただけたらうれしいです。 804 名前：デフォルトの名無しさん mailto:sage [2008/11/20(木) 14:58:26 ] 自分でかけよ 805 名前：デフォルトの名無しさん [2008/11/20(木) 17:45:43 ] こんばんは。 MATLABで近似曲線を描きたい場合、データがXとYの２軸のグラフでプロットされてるなら、 近似曲線用のコマンドがあるのでできます。しかし今回はXYZの三軸の要素があるデータをplot3コマンドで グラフ化したものに近似曲線をひきたいと考えていて、方法をずっと考えていたのですが全く思い浮かびません。 どなたかアイデアがありませんか？ 806 名前： ◆DcKrXqTCB2 mailto:sage [2008/11/22(土) 04:46:22 ] >>803 name=dir ------- gnuplotを、Windows 2000のOctave Ver3.1.50 Winows で使いたい、どうすれば？ 807 名前：デフォルトの名無しさん [2008/11/22(土) 13:26:48 ] 工学系の学生です。質問させてください。 常微分方程式の境界値問題用のbvp4cソルバや、偏微分方程式のpdeソルバは様々な微分方程式で汎用的に使えるものなのでしょうか。 実験系の研究室にも関わらず、実験結果を説明するための数値計算をやらなくてはいけないことになってしまったのですが、 自分でコードを作る技術がなく、既存のソルバで手早く計算できないものかと思っています。 情けない質問ですがよろしくお願いします。 これらのソルバを使う時に収束性や初期値の設定等にコツが必要でしたら、そのコツについても教えていただけませんしょうか。 「この本を読め！」でもうれしいです。 勉強用にこれ↓を買おうと思ってます。もっと良い本ってありますかね。 ttp://www.amazon.co.jp/MATLAB%E6%95%B0%E5%80%A4%E8%A7%A3%E6%9E%90-G-J-Borse/dp/427402380X 808 名前：ｐ [2008/11/23(日) 09:07:11 ] 上の本、高いですね−。 質問です。 plotで線をつなげずにヒストグラムのように描くにはどうすればいいですか？ （データは600行ほどのテキストファイルで、それぞれの行に数値が入っています。 それをdlmreadで読み込んでplotしたら横の軸は600までで思い通りなのですが それぞれのデータ点がつながってplotされてしまいました） 809 名前：デフォルトの名無しさん mailto:sage [2008/11/23(日) 12:19:59 ] hist()があるよ 810 名前：デフォルトの名無しさん mailto:sage [2008/11/24(月) 12:21:35 ] アッコにおまかせでやっていたけど、 古物商の免許を取って古本屋で、売れ残りかつ専門的な本を安値で買って アマゾンで売るという商売で月にだいたい１０万円、多いときで３０万円 の商売している人がいるらしいよ。 >>807の紹介した本もそういう人が提供している可能性大だね。 ではではノシ 811 名前：デフォルトの名無しさん mailto:sage [2008/11/24(月) 16:42:29 ] アッコにおまかせで微分方程式 まで読んだ。 812 名前：ｐ [2008/11/24(月) 19:25:16 ] 809さん、hist()をうまく使えばいいんですね！ありがとうございます。 813 名前：age [2008/11/26(水) 00:56:33 ] MATLAB R2008b studentバージョン を買いました。 mccを実行しようとしたところ、mccは用意されてなかったみたいです。 studentバージョンだからですか？ 814 名前：デフォルトの名無しさん mailto:sage [2008/11/28(金) 14:04:59 ] matlabの正式バージョンってマジでいくらするんですか？ 815 名前：デフォルトの名無しさん mailto:sage [2008/11/29(土) 17:04:20 ] >>814 MATLABだけ？それとも、Simulinkも？ まぁ、国産の軽〜高級車くらい見ておけばいいんでない？ 816 名前：デフォルトの名無しさん mailto:sage [2008/11/29(土) 22:52:25 ] 随分幅のある例を出してきたな。10倍以上あるじゃないか。 817 名前：デフォルトの名無しさん mailto:sage [2008/12/01(月) 19:50:28 ] ライセンスは細かく分かれているから、組み合わせでかなり変わるしな。 818 名前：デフォルトの名無しさん mailto:sage [2008/12/03(水) 13:38:54 ] 今日、MATLAB EXPOやってるな 聞いてくれば？ 819 名前：デフォルトの名無しさん mailto:sage [2008/12/05(金) 22:31:49 ] ある計算で、メモリ不足に陥っています。 512ＭＢ×４のメモリから、２Ｇ×２のメモリに増設しました。 zerosでどれだけのサイズを確保できるかを試したんですが、増設しても変わりありませんでした・・・。 10000×10000のゼロ行列でメモリ不足に陥ります。 サイズ増やせると期待したのですが、こんなもんなんでしょうか。 問題の計算のほうも改善できませんでした。どこかおかしいところがあればアドバイスお願いします。 バージョンはR2006b、ＯＳは32bitＸＰです。 820 名前：デフォルトの名無しさん mailto:sage [2008/12/05(金) 23:17:47 ] 10000x10000x8bytes=80,000,000 100Mそこそこだから確保は出来るはずだが、他にも巨大な変数があったら不明 メモリが断片化してる可能性もあるのでpackコマンドを使ってみる それでも足りなければ↓に従い64bitOSを使うなどの対処を。 www.cybernet.co.jp/matlab/support/techkwdb/search_view.php?id=1628 821 名前：デフォルトの名無しさん mailto:sage [2008/12/05(金) 23:47:44 ] 覆面算をエクセルのマクロで解きたいのですが教えてください Ｂ、Ｃ、Ｄは0を除外 　　　　Ｃ　Ａ　Ｂ 　+）　 Ｂ　Ａ　Ｃ ---------------------- 　　Ｄ　Ａ　Ｂ　Ａ 822 名前：デフォルトの名無しさん mailto:sage [2008/12/06(土) 00:31:47 ] >>820 800M だべ？ >>819 他のソフトも入れて (matlab 自体もいれて)、 32bit CPU && Windows なら 2GB の壁があるから、 そんなもんじゃねーがの。 823 名前：820 mailto:sage [2008/12/06(土) 14:35:28 ] >>882 おおう 単純な計算ミスorz 仮想メモリもふやして一応matlabが2G使えるはずだけど800Mなら他の変数しだいか。 やっぱ64bitでやるとかコマンドラインモードでやるとかが必要かもねえ 824 名前：デフォルトの名無しさん mailto:sage [2008/12/06(土) 15:21:17 ] >>820>>823 有難う御座います。 clearやって不要な変数などは一切消去して試した結果です。 因みにタスクマネージャーでmatlabのメモリ使用量みたら700Ｍ程度でした。 ウィンドウとか最小化すると200Ｍまで減ったりしますが、こんなもんでしょうか。。 やっぱり64ビット環境かぁ。。。 その場合、ゼロ行列ではどれだけのサイズが確保可能になるんでしょうか？ 825 名前：デフォルトの名無しさん mailto:sage [2008/12/06(土) 20:36:07 ] 変数名を文字列に変換することって出来ないのでしょうか？ 例えば、行列Aのプロット結果のtitleに「plot of A」とか書きたいんだけど、数が多くて手入力では限界が・・・ 826 名前：デフォルトの名無しさん [2008/12/09(火) 00:28:04 ] その変数をシングルクォートではさめば文字列になるよ。 827 名前：デフォルトの名無しさん mailto:sage [2008/12/09(火) 03:06:40 ] >>825 文字列から変数っつーか、式に eval で変換する、ほうが楽かも 828 名前：825 mailto:sage [2008/12/09(火) 21:52:58 ] >>826>>827 アドバイスありがとうございます。しかし私の説明が不十分でした。 私がやりたいことは、例えば、 func = [1 2 3] plot(func) で[1 2 3]がプロットされますが、そのタイトルに'func'を含めたいのです。 title(func) ではfuncの中身、つまり[1 2 3]がタイトルに表示されてしまいます。 title('func') だと、funcとは表示されますが、変数名自信(funcの部分)が変わる場合には自動で対応できません。 そこで、別変数(例えばfunc_name)に変数名を格納できれば title(func_name) とすることで、自動でタイトルを付けることが出来ると思ったのです。 そういう関数とかないかなと探しているのですが、なんだかなさそうな気配です・・・ 829 名前：デフォルトの名無しさん [2008/12/09(火) 22:41:47 ] title('func')だと自動対応できない、というが その例で変数名が変わったとき plot(func)は対応できるの? 君自身、自分が何をやりたいかが 分かってないのではないか。 830 名前：デフォルトの名無しさん mailto:sage [2008/12/09(火) 23:42:36 ] >>828 だから、変数名を格納して、eval で展開すれば、というとるやんけ 831 名前：デフォルトの名無しさん [2008/12/10(水) 03:25:58 ] eval使わなくてもMATLABには>825の目的の機能があるよ。 825が何やりたいかちゃんと説明できてないから予想だけど。 832 名前：デフォルトの名無しさん mailto:sage [2008/12/10(水) 06:55:04 ] >>831 その予想を書けよ。言われる前に書けよ。 833 名前：デフォルトの名無しさん [2008/12/10(水) 23:04:31 ] まともに質問できない人には質問になってないよって 教えてあげるのが親切ってもんだよ。 834 名前：825 mailto:sage [2008/12/11(木) 00:09:46 ] 不評を買ってしまったようで・・・ 実はプログラミング自体にあまり慣れていないもので、上手く伝えることが出来ないみたいです。 アドバイスを元にもうちょっと試行錯誤してみようと思います。 >>829 そこは泥臭い手法を採ってまして・・・ あまりに下手なので説明は遠慮させていただきます。 泥臭くて良いのなら引数で渡せばいいじゃんとか言われるとそれまでなのですが、 機能として有りそうにも拘わらず見つけることが出来なかったため質問した次第です。 >>830 ありがとうございます。調べて試してみます。 >>831 説明が上手くなくて申し訳ないです。 835 名前：デフォルトの名無しさん mailto:sage [2008/12/11(木) 01:48:56 ] >>834 >機能として有りそうにも拘わらず見つけることが出来なかったため質問した次第です。 そんな機能必要ないから存在しないだけだから 836 名前：デフォルトの名無しさん [2008/12/11(木) 01:57:53 ] プログラムに慣れてないとかいう問題ではたぶんなくて 自分がやってることとやりたいことを全部書いてないだけだと思うよ。 >825 > 変数名を文字列に変換する をしてくれるvarnameという関数が仮に存在したとしても あなたの質問の例では意味がない。 func = [1 2 3]; plot(func) title('func') を func = [1 2 3]; plot(func) title(varname(func)) と書けても「手入力」の節約になってない。 plotとtitleの部分を関数にしたいということなら分かるけど。 help inputname 837 名前：デフォルトの名無しさん mailto:sage [2008/12/11(木) 03:15:16 ] そう。だから文字列配列を作って eval で展開するのさ。 funcnames = { 'func1', 'func2', 'func3' }; for i = 1:length(funcs) 　eval(sprintf('plot(%s)', funcnames{i])); 　title(funcnames{i]); end >>831 eval 使わないとなるとどうやるんだ？参考までに教えてくれ。 838 名前：デフォルトの名無しさん [2008/12/11(木) 07:16:07 ] なんで>>837はえらそうなのか。 839 名前：デフォルトの名無しさん [2008/12/11(木) 17:53:58 ] > eval 使わないとなるとどうやるんだ？ >>836に書いてあるんだけど。 evalもいいけどデバッグしにくくなるし よほどの理由がある時以外は使わない方がいいよ。 840 名前：デフォルトの名無しさん mailto:sage [2008/12/12(金) 02:41:03 ] >>839 書いてないよ >varnameという関数が仮に存在したとして 841 名前：デフォルトの名無しさん [2008/12/12(金) 02:59:31 ] > help inputname 842 名前：デフォルトの名無しさん mailto:sage [2008/12/12(金) 07:44:47 ] >>841 は？inputname は全然ちゃうやろ。help みたんかいな 843 名前：デフォルトの名無しさん [2008/12/12(金) 07:57:51 ] 頭悪いのにえらそうな人ってやっかいだね。 >>828で func1 = [1 2 3]; func2 = [4 5 6]; plot(func1) title('func1') plot(func2) title('func2') とすると数が多くなったときにtitleを いちいち書くのが手間だっていうんだから plotfunc(func1) plotfunc(func2) function plotfunc(f) plot(f) title(inputname(1)) end でしょ。 844 名前：デフォルトの名無しさん [2008/12/12(金) 08:01:52 ] あ、テストせずに書いたけど plotfuncの中にfigure追加しないとだめですね。 もとの質問者がsubplotとか使ってるなら そのへんは適当に。 845 名前：デフォルトの名無しさん [2008/12/15(月) 18:02:33 ] m個のデータから正規方程式を使い A=Cexp(Bcosφ)を知りたい時どうすればいいですか？ Bは定数で、未知のCとφが知りたいんです。 未知定数が2つで積っていうのがよくわかりません 846 名前：デフォルトの名無しさん [2008/12/15(月) 18:20:25 ] 韓国に攻撃されてるぞｗｗｗｗ７ jfk.2ch.net/test/read.cgi/news4vip/1229319433/ 世界フィギュア選手権で浅田真央が韓国を破り金賞に輝く ↓ 韓国人が腹いせに2ちゃんねるを攻撃 ↓ vipper今夜19:00 反撃凸　←いまここ★ 目標 www.dcinside.com/ Dcinside（ディーシーインサイド）は、韓国の電子掲示板サイトである。 デジタルカメラの情報を提供する専門サイトとして有名になりつつ、だんだん人が集まり韓国最大のコミュニティーサイトとなった。 韓国の2ちゃんねるともいえるサイト。 2004年に2ちゃんねるのサーバをDoS攻撃したことで有名。 847 名前：デフォルトの名無しさん mailto:sage [2008/12/17(水) 18:08:04 ] ｍ質点の串団子モデルの剛性マトリクスをfor文を使って次のように書いたのですが、 for文を使わないで書く方法はありますでしょうか？ for文をなるだけ使わないのがMATLABっと以前のレスで見たのですが、 何分プログラミング自体初めてでして。。 row_m = input('Enter m = '); k0 = input('Enter 0 or 1 ='); switch k0 case 0 k1 = input('Enter k = '); kk = k1*ones(1,row_m); otherwise kk = zeros(1,row_m)+input('Enter kmatrix = '); end K = kk(:,1); for N = 2:row_m K0 = [kk(:,N),-kk(:,N);-kk(:,N),kk(:,N)]; K1 = [K,zeros(N-1,1);zeros(1,N-1),zeros(1)]; K2 = [zeros(N-2),zeros(N-2,2);zeros(2,N-2),K0]; K = K1+K2; end ここはこうしたほうが良いなどなんでも結構ですのでご助言を頂けると幸いです。 848 名前：デフォルトの名無しさん mailto:sage [2008/12/18(木) 00:06:06 ] for N = 2:row_mなくして Nの代わりにまんま2:row_mを使え 849 名前：デフォルトの名無しさん [2008/12/19(金) 03:19:52 ] すいませんがよろしくお願いします。 waveファイルをインポートするにはどうすればいいのでしょうか？ helpをみてもできません。 [y,Fs,Bits] = wavread('onsei1.wav'); と入力したら、 ??? エラー ==> wavread at 67 ファイルを開けません とでます。 誰か教えてください！ 850 名前：デフォルトの名無しさん mailto:sage [2008/12/19(金) 04:42:52 ] plotで図示したグラフとcontourで図示したグラフを重ねて表示する際、 X軸Y軸が逆になってしまいます。 contour(y,x,z)とすれば軸はあうのですが、いまいち納得いきません。 どうしてこうなるのでしょうか？ 851 名前：デフォルトの名無しさん [2008/12/20(土) 07:17:39 ] 画像をimresizeでbicubicで補間しながら拡大・縮小したときの動作確認がしたいのですが どうしたらいいでしょうか？ もしくはbicubicで補間しつつ拡大縮小を少数精度の倍率を与えたときに正確に動作が 保証できるプログラムを組んだことのある方いらっしゃいませんか？ ソースをアップして欲しいのです・・・。修士論文が間に合わない・・。 852 名前：デフォルトの名無しさん mailto:sage [2008/12/22(月) 12:44:18 ] Y = find([0 1 0 0 1]) Y 　2 5 の逆ってどうやるんでしたっけ？ 853 名前：デフォルトの名無しさん mailto:sage [2008/12/22(月) 12:49:22 ] >>849 ファイルがそこにないんじゃない？ >>850 meshgrid を使ったと仮定して、 meshgrid の使い方間違ったんじゃない？ 854 名前：デフォルトの名無しさん [2008/12/24(水) 00:27:15 ] >>852 X(Y)=1 855 名前：デフォルトの名無しさん [2008/12/25(木) 13:57:23 ] プログラム初心者です。 MATLABを使えるようになりたいんですが 基本的な使い方とか解説してあるいいサイトとかあったら教えていただけないでしょうか。 856 名前：デフォルトの名無しさん mailto:sage [2008/12/25(木) 15:36:00 ] サイトではないですが、「MATLAB+Scilabﾌﾟﾛｸﾞﾗﾐﾝｸﾞ事典」という書籍にはお世話になりました。MATLABに似たscilabについても記載されていますが。 後はMATLABのhelpでしょうか。 857 名前：855 mailto:sage [2008/12/26(金) 11:56:00 ] >>856 ありがとうございます。 探してみようと思います。 858 名前：デフォルトの名無しさん mailto:sage [2008/12/28(日) 22:52:26 ] Student Editionって以前は自分名義のクレジットカード+HPへの登録+学生証のコピーの郵送が必要って Q&Aに書いてあったと思うのですが、今は生協で買えるんですねー その場合は上記の郵送とか面倒な作業は必要ないんでしょうか。 859 名前：デフォルトの名無しさん mailto:sage [2008/12/29(月) 23:57:12 ] >>858 生協で学生証提示でオケ 860 名前：デフォルトの名無しさん mailto:sage [2009/01/04(日) 02:36:50 ] >>859 返事遅くなってすいません。 煩雑な手続きは無くなったんですね。ちょっと興味有ったので買ってみることにします。 ありがとうございました。 861 名前：デフォルトの名無しさん mailto:sage [2009/01/06(火) 00:42:52 ] scilabってコマンドはマトラボとおなじなんでしょうか？ 家でもやろうかと思ってるんですが 862 名前：デフォルトの名無しさん [2009/01/06(火) 00:54:09 ] 言い直すと学校でマトラボを使っているんですが 家でもやりたいのでscilabをだいようできるでしょうか？ 863 名前：デフォルトの名無しさん mailto:sage [2009/01/06(火) 01:35:49 ] www.ipc.yamanashi.ac.jp/bul/final01/ohki/matlab.html >Scilab の文法は MATLAB と良く似てはいるが， 細部には違いがある． >Octave は， 基本的な言語仕様の点で MATLAB 互換となるように設計されており， >Octave で作成したプログラムの多くは修正なしでそのまま MATLAB で実行可能である >(MATLABの方が機能が豊富なので，逆は真でない)． 864 名前：デフォルトの名無しさん mailto:sage [2009/01/06(火) 03:08:02 ] >>863 わかりやすくありがとうございます 865 名前：デフォルトの名無しさん mailto:sage [2009/01/06(火) 13:51:44 ] freematは? 866 名前：デフォルトの名無しさん mailto:sage [2009/01/06(火) 15:34:35 ] octaveを語ろう pc11.2ch.net/test/read.cgi/unix/1086016412/ Scilabを触ってみるｽﾚ science6.2ch.net/test/read.cgi/kikai/1068735555/ SCILAB&SCICOSについて教えてください science6.2ch.net/test/read.cgi/sim/980480122/ どれも過疎化しているが… freematスレは見当たらなかった。 867 名前：MATRAB_Mathematica [2009/01/08(木) 21:47:19 ] 下記をMathematica上で使うにはどうすればいいですか？ Seq = {-0.2266091476053614,0;[0.2155407618197651,-1.2670686037583443],0; [-4.2551584226048398,0.5047757263881194],2; [0.2331599353469357,0.0447459687134724],-2; [6.6244572505007815,-18.389000853969371],4; [-0.0567684937266291,0.1443950619899142],-4;-5.5119344180654508,5}; ScaleS = -0.5985483742581210; ScaleD = -1.6707087396895259; Family = 'Symlet'; 868 名前：デフォルトの名無しさん mailto:sage [2009/01/13(火) 16:23:57 ] matlabでH.264動画を表示させるにはどうしたらよいのでしょうか。 「mmreader」で読み込んでも走査線が分かれて表示されます。 インターレース動画は奇数偶数をいじって解除・表示できたのですが、 H.264やmp4コーデックでエンコードされた動画は上1/3が黒く表示され、 走査線がバラバラに表示されます。 H.264動画の走査線をどのように分解したらよいのか詳しくわかりません。 詳しい方がいらっしゃいましたらよろしくお願いします。 869 名前：デフォルトの名無しさん mailto:sage [2009/01/16(金) 16:05:57 ] プログラミング初心者です。画像処理の質問です。 ある画像から2値化を行い，1つの物体だけを取り出そうと考えています。 bwlabelを使ってラベル画像を作った後に，ラベル付けされた各物体の特徴量として面積と円形度を求め，ある範囲に入っている物体1つを取り出したいのですが，よくわかりません。 調べたところbwselectがあったのですが，これは画像をクリックする操作があるため出来れば避けたいと思っています。 わかる方がいましたらアドバイスを頂けないでしょうか？ よろしくお願いします 870 名前：デフォルトの名無しさん mailto:sage [2009/01/26(月) 18:35:45 ] MATLABを使っていて分からないことがあったので、画像処理について質問させてください。 0〜1の間の乱数を128×128の行列に入れて、これをXと置いた後、 imshow(X,128)で濃淡レベル128のグレースケールで表示します。 これで、白〜黒の画像が表示されるのですが、これを白〜赤のように色をつけて出力したいです。 colormap(Red);imshow(X,128) でできるかと思ったのですが、上手くいきません。どうすればよいか教えていただけないでしょうか。 871 名前：デフォルトの名無しさん mailto:sage [2009/01/27(火) 10:49:11 ] >>870 MATLABよく知らんから逆に聞くけど、colormap(Red)するとどう出力されるの? もしかして、黒-灰-白が黒-赤-白になるとか? colormapの仕様をよく調べてみるといいんでない? 872 名前：デフォルトの名無しさん mailto:sage [2009/01/27(火) 13:32:39 ] Redがサポートされてないってオチだろ Hotとかにしてみ 873 名前：デフォルトの名無しさん [2009/01/27(火) 21:20:12 ] 2次元画像処理で、ハイパスフィルタを書いたのですが このコードで正しいのか自信がもてません。 何方か、分かる方がいらっしゃいましたらご指摘願います function I=highpass2(I,THV) % ハイパスフィルタ function I=hp2(I,THV) I=double(I); I2=fft2(I); I3=abs(I2); [m,n]=size(I); for k=1:m; for l=1:n; if( I3(k,l)<THV ) I2(k,l)=0; end end end I = uint8(real(ifft2(I2))); end I = hp2(I,THV*1000); end 874 名前：873 [2009/02/03(火) 20:31:01 ] スレ汚し失礼いたしました。 どうも上のコードで良かったみたいですので、高速化ver貼り付けておきます。 function I=highpass2(I,THV) % ハイパスフィルタ % THVは0〜255の間でお願いします(0〜10推奨) function I=hp2(I,THV) I=double(I); I2=fft2(I); [m,n]=size(I); I2=I2.*double( im2bw( abs(I2)./(m*n*255), min(THV./255,1.0) ) ); I = uint8(real(ifft2(I2))); end [m,n,o]=size(I); for i=1:o I(:,:,1) = hp2(I(:,:,1),THV); end end 875 名前：デフォルトの名無しさん [2009/02/05(木) 12:28:10 ] Linux用Matlabについての質問です。 32 bit版のmatlabがあるのですが、64bitカーネルのLinuxで動きますか？ 876 名前：デフォルトの名無しさん mailto:sage [2009/02/05(木) 13:13:25 ] つかってるLinuxによるとしかいえない。 877 名前：デフォルトの名無しさん [2009/02/10(火) 07:48:30 ] plotでデータ数5000の系列を、同じグラフに5本書きます。 それぞれの系列をカラーでわけたくない場合はマーカで区別すると思うのですが、 データ数が5000だとどのマーカも潰れて見えません。 データ数100個飛ばしでマーカをつけるなど、そのような飛び飛びでマーカをつける方法はありますか？ 878 名前：デフォルトの名無しさん [2009/02/10(火) 11:47:31 ] カラーバーにおいて、特定の値の色を指定する方法を 教えて頂けないでしょうか。 今、カラーバーの範囲を[a, b] (a < 0 < b) にして、 0 の色を白に指定したいのですが、苦労しています。 よろしくお願いします。 879 名前：デフォルトの名無しさん [2009/02/14(土) 07:56:25 ] www.osaka-kyoiku.ac.jp/~ashino/topics/wavehist.html? このページの真ん中当たりにコッホ曲線のウェーブレット変換が書かれてあるのですが、 この上側の図に書かれてあるコッホ曲線はよく見かける三角形をたくさん合わせたものと全く違いますが、何を表しているのでしょうか？ 下の図は3次元的な図になっていますが、それぞれ横、縦、手前の軸は何を表しているのでしょうか？ 880 名前：デフォルトの名無しさん [2009/02/22(日) 20:44:57 ] eval関数について質問です． 変数を文字列化してテキストボックスとしてfigureに挿入しようと 思っているのですが，以下のようにやってもうまくいきません． eval(['text(20,70,',num2str(k),',''FontSize'',16,''FontName'',''Century'');']); 基本的に([])内の'と,の使い方についてよく理解していません． '〜',〜,'〜',〜,…　のような順番で結合したい要素を配置して， ''内に更に'〜'のような表現がある場合''〜''に書き換える決まり？が あるんじゃないかと勝手に解釈して上記のコマンド書きました． 非常にややこしい文になりましたが，そこらへんの事も含めて ご教授いただきたいです.よろしくお願いします． 881 名前：デフォルトの名無しさん mailto:sage [2009/02/23(月) 22:13:59 ] >>880 俺の脳内MATLABで実行してみたら、こんな感じだと俺の妄想どおりになった eval(['text(20,70,''',num2str(k),''',''FontSize'',16,''FontName'',''Century'');']); >>''内に更に'〜'のような表現がある場合''〜''に書き換える そうそうそんな感じ dl.cybernet.co.jp/matlab/support/manual/r13/toolbox/matlab/ref/?/matlab/support/manual/r13/toolbox/matlab/ref/eval.shtml 882 名前：880 [2009/02/24(火) 00:16:41 ] >>881 仰るとおりのコマンドでうまくいきました． リンク先も一応見た事あるのですが，evalで長いコマンド括る方法が無かった物で… おそらくevalで一括するよりも，先に文字列作るなりしたほうが良いのかも知れませんね． ありがとうございました． 883 名前：デフォルトの名無しさん mailto:sage [2009/02/25(水) 13:58:40 ] plot([0.1:0.1:0.9 0.91:0.01:0.99], [0.1:0.1:0.9 0.91:0.01:0.99].^2) こんな意味のプロットを plot(1:18, [0.1:0.1:0.9 0.91:0.01:0.99].^2) のような表示でプロットしたい（x 軸の点ごとの間隔を等間隔にしたい）のですが、 上のようにプロットすると、x 軸のラベルがもちろん、1,2,3,4 のようになってしまいます。 このラベルは　[0.1:0.1:0.9 0.91:0.01:0.99]　のままにしたいのですが、 どうようにすれば実現できるでしょうか？よろしくお願いします。 884 名前：デフォルトの名無しさん mailto:sage [2009/03/01(日) 09:09:55 ] >>883 plot(1:18, [0.1:0.1:0.9 0.91:0.01:0.99].^2) でプロットしてからXTickLabelを変更すればいいんじゃない？ 885 名前：デフォルトの名無しさん mailto:sage [2009/03/04(水) 20:50:40 ] そんな姑息な真似はいやなんです。 886 名前：デフォルトの名無しさん mailto:sage [2009/03/04(水) 22:08:33 ] そうかい、すまんね。 887 名前：デフォルトの名無しさん [2009/03/06(金) 14:46:30 ] 既出かもしれませんが、 たとえば、(z-0.99)/(z-0.58)みたいな伝達関数があって、 任意のzdをｚに代入したい場合、つまり、ｚｄ＝0.45＋0.53ｊをｚにいれるには どうしたらできるでしょうか？ それとも、そもそもMatLabではできないですか？ 888 名前：デフォルトの名無しさん mailto:sage [2009/03/07(土) 01:41:31 ] 伝達関数、と限定した言い方をするってことは こんな↓単なる台数方程式的な話ではないんですよね。 z=0.45+0.53*j x=(z-0.99)/(z-0.58) 889 名前：デフォルトの名無しさん mailto:sage [2009/03/07(土) 10:47:17 ] なんだか分野を書かない質問が多いよね。 信号処理の質問を経済学の人が見てもしょうがないのにさ。 ちなみに filter だったと思う。 signal processing toolbox に一度ひととおり目をとおしておくとよいでしょう。 890 名前：デフォルトの名無しさん mailto:sage [2009/03/11(水) 02:58:33 ] >>888 方法が見つからなかったので、結局、代数方程式を新たに作って解決したのですが、 もし、tfで作った関数に当てはめることができれば、そういった手間を省けると 思って質問しました。 >>889 すいません。分野は機械工学のデジタル信号処理です。離散時間を扱ってます。 あと、 そういったツールボックスがあるのですね。勉強になりました。 891 名前：デフォルトの名無しさん [2009/03/17(火) 22:46:36 ] どうもすいません。matlabのコードを.Net向けにコンパイルを行ってその結果を受け取っているのですが function [a,b] = func( c ) の様に引数を2つ以上返すとき、どうすれば２つとも受け取れるのでしょうか？ というより文法的に何か可能は方法があるのでしょうか？ 892 名前：デフォルトの名無しさん mailto:sage [2009/03/17(火) 22:57:06 ] >>891 何の言語を使っているのか知らんが、その言語スレで聞いたら? まぁ、構造体を作ってって処だろうけど。 893 名前：abc [2009/03/18(水) 16:27:20 ] 行列の列ごとに、0で無い値の中での平均を求める方法を教えてください。 894 名前：デフォルトの名無しさん mailto:sage [2009/03/18(水) 23:53:24 ] for i = 1:size(A,2) mu(i) = mean(A(A(:,i) ~= 0, i)); end いや、 mu = sum(A, 1) ./ sum(A ~= 0, 1) こっちのほうがスマートかな 895 名前：デフォルトの名無しさん mailto:sage [2009/03/19(木) 08:37:26 ] meanの引数に直接行列をぶち込むのと結果が変わらなくね? 896 名前：デフォルトの名無しさん mailto:sage [2009/03/19(木) 11:15:21 ] mean なら0もカウントしちゃうぜしょ。 897 名前：デフォルトの名無しさん [2009/03/19(木) 22:33:19 ] figureで軸を指定した関数で表示するにはどうしたらいいでしょうか？ logプロットの様な感じでlog以外の指定した関数で表示したいです。 898 名前：デフォルトの名無しさん mailto:sage [2009/03/24(火) 06:48:41 ] 基本的な0次、1次のBessel関数でMATLAB組み込みより速いのってありますか?

---

---

[ 新着レスの取得/表示 (agate) ] / [ 携帯版 ]

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