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 307 名前：デフォルトの名無しさん mailto:sage [2007/12/04(火) 06:02:50 ] 関数の入力にベルトルの集合としての行列をいれるとき皆は NxD 派？ DxN 派？where N is the number of vectors and D is the dimension of each vector. どっちのほうが親切なんだろう？toolbox によって（人によって）違うんだよねぇ。 308 名前：デフォルトの名無しさん mailto:sage [2007/12/04(火) 09:55:13 ] シミュレート板で回答を得られなかったので。 simulinkである変数がある値になったときシミュレーションを終了したいのですが どうすればできますか？ 309 名前：デフォルトの名無しさん mailto:sage [2007/12/05(水) 06:53:25 ] ムービーを作成する際に、figureウインドウを出さずに作成することはできないでしょうか？ figureウインドウを作成してgetframeで取り込む場合、スクリーンキャプチャをしているようで、 他の作業をすると、その画像がムービーの中に入ってしまいます。 そのためムービー作成中は何も作業ができません。 310 名前：デフォルトの名無しさん [2007/12/05(水) 16:25:53 ] DICOMファイルから画像読み込み、JPG形式で書き込みしようとしてるのですが、 X=dicomread('--.dcm'); [Y,map] = gray2ind(X); imwrite(Y,map,'--.jpg','jpg'); とすると、colormapを使用したimshowの表示画像と異なってしまいます。 よろしくお願いします。 311 名前：デフォルトの名無しさん [2007/12/05(水) 18:06:21 ] 正の実数からなる対角行列の、各行と列について、 要素の大きさが対角成分の要素の大きさの0.5から2倍の範囲なら１で置換、 それ以外なら０で置換する、 ってforループを使わずにできるでしょうか？ 312 名前：デフォルトの名無しさん mailto:sage [2007/12/05(水) 23:19:50 ] 質問ばっかで誰も答えてくれない罠ｗｗｗｗｗｗｗｗｗ 313 名前：デフォルトの名無しさん mailto:sage [2007/12/06(木) 00:59:41 ] >>311 対角行列なら対角成分以外すでに全部０やん 314 名前：デフォルトの名無しさん mailto:sage [2007/12/06(木) 02:21:48 ] これぞまさに質問スレ 315 名前：デフォルトの名無しさん mailto:sage [2007/12/06(木) 08:48:34 ] vistaは2007aで正式サポートですが、2006aで動かした場合、どこかに既知の問題のようなリストはないでしょうか？ 316 名前：デフォルトの名無しさん mailto:sage [2007/12/07(金) 01:04:38 ] >>307 僕はMatlabで音系のことをやってるけど、 wavwriteとかのデフォルトが行＝サンプル／列＝チャンネルなので、 ベクトルの集合はいつも縦長ベクトルを横にいくつか並べるというかんじ。 DxNってことになるのかな？ 317 名前：デフォルトの名無しさん mailto:sage [2007/12/07(金) 12:48:35 ] u1=[ [ { arctan( ( 1-k ) * cot( u1*v*w ) -b*tan(u1*v*w) ) } /( 2*a) ] + ( q*π ) ] / [ { arctan( (1-k) * cot(v*w) -b*tan(v*w) ) } / (2*a) ] u2=[ [ { arctan( ( 1-k ) * cot( u2*v*w ) -b*tan(u2*v*w) ) } /( 2*a) ] + ( r*π ) ] / [ { arctan( (1-k) * cot(v*w) -b*tan(v*w) ) } / (2*a) ] a=sqrt(k)*(y/x-x/y) b=k*(y/x)^2-(x/y)^2 matlab流の書き方じゃなくてすみません。 上記の連立非線型方程式からxとwを数値解として求めるプログラムをつくりたいのですが・・・ x,w 以外の値はすべて既知で、 u1,u2 は 1<u1<u2<=3 である任意の実数 q,rは 0または1 それ以外は任意の正の実数 分母は必ず正の値(になるはず) 最適解を求める際に、分子のarctan(・・・)の項が負の値になってしまう時にはqまたはrを「1」にしてπを加える それ以外の場合にはq,rは「0」 と、こんな感じなのですが,どなたかサクっと作ってくださいませんでしょうか・・・ 318 名前：デフォルトの名無しさん mailto:sage [2007/12/12(水) 05:19:53 ] 秀丸用の強調定義ファイルはありませんか？ 319 名前：デフォルトの名無しさん mailto:sage [2007/12/12(水) 10:23:37 ] matlab ってファイル読み込み関連の関数が貧弱で困った。 一行にいくつ数があるかわからないときにどうやって読み込んだらいいんだ？ ちなみにこんなファイルフォーマット 1 1:8392 2:8384.23 4:9832.82 2 2:8329 3:8389.34 5:9328.82 3 4:9838 先頭の 1,2,3 は最終的には捨てる行数 1:8392 は要素[1]に8392という値ということ。そして最初の行は 3 を飛ばして 4 に行っている。 そして最後の行は、他の行と読みこむ数が違う。 320 名前：デフォルトの名無しさん [2007/12/12(水) 10:36:23 ] fgetlで１行まるまる読んで、strtokで切り出していくしか無いんじゃないの？Ｃでも普通そうだろ？ 321 名前：デフォルトの名無しさん mailto:sage [2007/12/13(木) 04:09:16 ] >>319 それを一発で君の思い通りに読んでくれる ファイル読み込み関連の関数が貧弱じゃないやつって何？ 煽ってるんじゃなくて単に興味があるから知りたい。 ExcelとかRとかできるんかな? 322 名前：デフォルトの名無しさん mailto:sage [2007/12/13(木) 12:21:50 ] excelじゃ無理だね。 323 名前：デフォルトの名無しさん mailto:sage [2007/12/13(木) 13:11:34 ] figureウインドウを表示させずにムービーを作る方法を教えてください 324 名前：デフォルトの名無しさん [2007/12/13(木) 13:25:52 ] >>319と似た質問なんですが フリーフォーマットで他のソフトで出力されたデータを取り込みたいのですが textreadがうまくいかず… 1 1.00000E+00 2.00000E+00 3.00000E+00 4.00000E+00 -0.40763E-13 -0.22240E-04 2 5.00000E+00 6.00000E+00 7.00000E+00 8.00000E+00 -0.84319E-13 -0.22422E-04 と並んでるデータがありまして、このデータの６つ並んだうちの 後ろ二つを取り込んで行ベクトルにしたいんです。 で、textreadを用いて [head1,data1,data2,data3,data4,data5,data6,] = textread('ファイル名.txt', '%u %f %f %f %f %f %f') とやってみたものの、うまくいきません。 どうやら１の後ろに何かがあるとmatlabが認識し（ホントは改行だけのはず） data1〜6 (%fの指定部分)に0をいれとります。 そして次の行に移行し、そこではきちんとデータを拾えてます。 周波数を表す　１　と　次の行の　1.00000E+00　を取り込むとき %u %fだとフォーマットの指定に不足があるのでしょうか？？ 325 名前：デフォルトの名無しさん mailto:sage [2007/12/13(木) 14:21:09 ] >>324 textread は一行ずつ読み込むんだろ？ 最初の行は先頭以降何もないから０をが入っていくと。 [data1,data2,data3,data4,data5,data6,] = textread('ファイル名.txt', '%f %f %f %f %f %f') として、奇数行と、偶数行に後でわけるとかになるんじゃね？ 326 名前：デフォルトの名無しさん [2007/12/13(木) 15:38:43 ] どうやらそのようだ、と思っていたんですが 以前同様のものを読み込んだ際に行ごとでなかったんですよね たとえばこんなデータ 1 1 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.58726901E-12 0.55467926E-06 を読み込むとき、 [head1,head2,data1,data2,data3,space1,space1,data4,data5,data6] = textread('testATV.txt','%u %u %f %f %f %c %c %f %f %f'); で、data5,data6に 0.58726901E-12 0.55467926E-06 を格納できたんです。 これって、五つ目の%f（浮動小数点）のあとに%c（スペース）を読み込んでるので 行をまたいでも効果があることになると思うんです。 327 名前：デフォルトの名無しさん [2007/12/13(木) 15:39:32 ] 1 1 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.58726901E-12 0.55467926E-06 328 名前：デフォルトの名無しさん [2007/12/13(木) 15:41:18 ] あああすいません、変な書き込みしてしまいました 1 1 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.58726901E-12 0.55467926E-06 329 名前：デフォルトの名無しさん [2007/12/13(木) 15:42:56 ] どうやらうまく書き込めないのですが １番目、２番目の１の下にはスペースが存在していて ３番目の0.00000000E+00 と下の行の0.00000000E+00 は 列がそろっています。 330 名前：ks [2007/12/16(日) 14:41:45 ] ＝＝＝＝＝＝＝＝ 教えてください。 ＝＝＝＝＝＝＝＝ １次遅れ＋無駄時間要素を含む制御対象の限界感度法による PIDパラメータをmatlab上で、どう求めるんですかね〜 ＝＝＝＝＝＝＝＝ 教えてください。 ＝＝＝＝＝＝＝＝ 331 名前：デフォルトの名無しさん [2007/12/21(金) 01:33:00 ] >>329 textreadにこだわらなくてもいいんじゃないの？ 332 名前：デフォルトの名無しさん [2007/12/27(木) 03:52:31 ] matlabでfor文をまわす時 ループが終わるごとに変数をclearした方がよかったりします？ 不振な挙動。。 333 名前：デフォルトの名無しさん mailto:sage [2007/12/27(木) 09:04:36 ] そんなことないと思うが。 それにしても、Matlabバージョンとか書かないと。 334 名前：デフォルトの名無しさん mailto:sage [2007/12/28(金) 00:55:49 ] たぶんそれはいくつかの変数を間違って使ってるかロジック間違えてるだけだ。 ちゃんと全部読み直せ。 335 名前：デフォルトの名無しさん [2008/01/07(月) 02:46:00 ] 波形処理をやっていてカーブフィッティングをしています。 fit関数を使ってフィッティングをしているのですが、 フィッティング後の曲線のパラメータが読み取れません。 fit1 = 一般ﾓﾃﾞﾙ Exp1: fit1(x) = a*exp(b*x) 係数(95% 信頼区間付き): a = -10 (-10.17, -9.831) b = -0.01371 (-0.01403, -0.01338) こんな感じでフィッティングしたあと、bの値を使いたいのですが、 どうすれば読み取れますか？ 初心者の質問ですみません。 336 名前：デフォルトの名無しさん [2008/01/08(火) 00:30:09 ] 2値化した画像をカラー画像にするにはどうしたらいいですか 337 名前：デフォルトの名無しさん mailto:sage [2008/01/08(火) 00:51:57 ] それはもう情報がかなりないからとても難しいあるね 338 名前：336 [2008/01/10(木) 14:39:29 ] >>337 カラーは白黒でよいのですがそれでも無理ですか? 339 名前：デフォルトの名無しさん mailto:sage [2008/01/10(木) 16:15:22 ] >>338 それは、フォーマット変換したいと言うこと? それとも、2値のメッシュをグレイにしたいと言うこと? 340 名前：デフォルトの名無しさん [2008/01/10(木) 22:09:11 ] 2次元座標から、各座標の点をラインで連結してプロットさせ それを画像として処理したいのですが、現在は plot(point(:,1),point(:,2));または line(point(:,1),point(:,2)); pointは100*2行列、xy座標の点100個としています。 これだとプロット結果を保存しても画像のサイズが変更されてしまいます。 点と点を直線で結んでいく関数などありましたら教えていただけないでしょうか？ 341 名前：デフォルトの名無しさん mailto:sage [2008/01/11(金) 04:51:05 ] 行列の積で初歩の初歩的な質問だと思うのですが誰か教えてください！ 10000x10の行列Aがあるとします。 A*A'だとout of memoryと出力されるのですが A'*Aだとちゃんと計算結果が出力されます。 考えてもイマイチ理由が浮かんで来なかったのでどなたか教えてください〜！ 単純な質問で失礼しました。 342 名前：デフォルトの名無しさん mailto:sage [2008/01/11(金) 07:34:46 ] 前者は 10000x10000 後者は 10x10 の出力になるだろう？明らかすぐる 343 名前：デフォルトの名無しさん mailto:sage [2008/01/11(金) 08:40:10 ] これはひどい 344 名前：デフォルトの名無しさん [2008/01/11(金) 09:11:51 ] （・A・）＊（・A・）' 345 名前：デフォルトの名無しさん mailto:sage [2008/01/11(金) 15:54:33 ] ホントだｗｗｗ頭悪くてすまそｗｗ ありです。 346 名前：デフォルトの名無しさん [2008/01/12(土) 16:46:00 ] plotで描ける線の種類を増やせないでしょうか？ 6種類のデータを1枚の図に、白黒で書く、ということで困っています。 347 名前：デフォルトの名無しさん mailto:sage [2008/01/12(土) 21:13:17 ] >>346 解決しないかも知れないけどマーカと併用するのはダメなんですか 線種で区別できる程に線が絡み合ってないのならそれでいけるのではないでしょうか 348 名前：デフォルトの名無しさん mailto:sage [2008/01/12(土) 22:26:34 ] >>339 フォーマット変換のほうだと思います。 349 名前：デフォルトの名無しさん mailto:sage [2008/01/15(火) 10:21:20 ] >>348 MATLABでなくていいなら、画像ビューワでも画像エディタでもいくらでもやりようがありそうだが。 何を、何故、MATLABでやりたいのか全然判らんのだよ。 350 名前：デフォルトの名無しさん mailto:sage [2008/01/15(火) 22:46:03 ] >>349 MATLABで画像処理の勉強をしろといわれ、つい2値化の画像を白黒だけの カラー画像に変えてやろうと思った。 351 名前：デフォルトの名無しさん mailto:sage [2008/01/16(水) 00:07:40 ] まずおめーは、白黒だけのカラー画像とかわけわかんねーことをいうな カラー画像っつったらＲＧＢとか３チャネルあるやつを普通は指すんだよ。 352 名前：デフォルトの名無しさん mailto:sage [2008/01/16(水) 00:15:23 ] あほす R == G == B のことだろ それより2値化の画像って何だ？ 2値の画像だろ 353 名前：デフォルトの名無しさん mailto:sage [2008/01/17(木) 06:03:13 ] ステップ応答の試験をしてMATLABで伝達関数を出したいんですけど、 最小二乗でフィットさせるときの方法とかいうのが分からないので教えてください。 354 名前：デフォルトの名無しさん [2008/01/17(木) 12:20:28 ] テキストファイルの行数を事前に知る関数はないでしょうか？ 355 名前：デフォルトの名無しさん [2008/01/17(木) 17:09:55 ] >336 エスパー回答すると、im2bw関数とかで2値化してしまった画像は、0が黒、1が白になります。 そしてどちらかの値しか取れない。 これをグレースケールの画像に変換したい場合は、0〜255の値を扱えるように拡張して、その後に掛け算をしてあげると良いです。 具体的には BW; %　２値化された画像 I = uint8(BW)*255; % グレースケール化された画像 な、感じです。 実際は、Image Processing Toolbox内のim2uint8とかいう関数が同じ事をしているよ。 RGB化したい場合はさらに [w,h] = size(I); CI = zeros(w,h,3); % カラーイメージを入れるための準備 CI(:,:,1) = I; % 赤に代入 CI(:,:,2) = I; % 青に代入 CI(:,:,3) = I; % 緑に代入 …でいいはず（うろ覚え 356 名前：デフォルトの名無しさん mailto:sage [2008/01/18(金) 00:32:14 ] >>353 どういう関数系にフィットさせるの？ 357 名前：デフォルトの名無しさん mailto:sage [2008/01/18(金) 06:14:54 ] >>357 よくある１次遅れ、2次遅れ、むだ時間とかいうパラメータを使った低次元伝達関数です たとえば G(s) = Kωn^2 / (s^2 + 2ζωns + ωn^2) * e^(-τs) のパラメータを出したいんです 358 名前：デフォルトの名無しさん mailto:sage [2008/01/18(金) 06:15:39 ] 安価ミス ×　>>357 ○　>>356 359 名前：デフォルトの名無しさん mailto:sage [2008/01/18(金) 06:31:51 ] あ、あと私のMATLABにはSystem Identification Toolboxが入っていないのでこれを使わない方法を教えて欲しいです。 360 名前：デフォルトの名無しさん [2008/01/18(金) 08:58:30 ] そこまでわかっているのなら、最小自乗でフィットするだけじゃん ttp://www.wakayama-u.ac.jp/~miw/ApIp/01/2.html 361 名前：デフォルトの名無しさん mailto:sage [2008/01/18(金) 10:53:24 ] 指数関数的に減衰していくデータにフィットさせたいのですが、このように非線形の場合どのようにやればよいのでしょうか？ 対数プロットでは直線になるので、何か最小自乗法のような物は使えないでしょうか？ 362 名前：デフォルトの名無しさん [2008/01/18(金) 17:47:29 ] fminsearch関数使えば？Simplex方による繰り返し計算だけど。 363 名前：デフォルトの名無しさん mailto:sage [2008/01/18(金) 19:14:34 ] 数学的には間違いだけれども、 対数変換して直線フィットすると、ある値が出る。 それを初期値に>>362にしたがって、解を探す。 364 名前：デフォルトの名無しさん mailto:sage [2008/01/19(土) 12:02:15 ] >>360 ごめんなさい その最小二乗でだけっていうところが分からないでの聞きにきたんですよね・・・ 365 名前：デフォルトの名無しさん [2008/01/19(土) 12:47:42 ] polyfit使えば一発。ただ原点に返って正規方程式を逆行列で解くのもありだな。 366 名前：デフォルトの名無しさん [2008/01/19(土) 14:03:31 ] >>364 MATLAB以前の問題なのか… >>365 polyfitって、フィットさせる式に制限があるからね。 367 名前：デフォルトの名無しさん [2008/01/19(土) 21:23:54 ] >>365 ployfitって単なる多項式しか出来ないんですよね >>366 MATLAB以前の問題かもしれないんですけど、 やるとしたらデジタルで出てきたデータをz変換してそれを行列使って最小二乗してやるということになるのでしょうか？ 368 名前：デフォルトの名無しさん [2008/01/20(日) 15:18:25 ] だから単純なたとえばy = a*x^2 + …とか>>360に書いてあることはわかるんだが、 伝達関数の場合のプロセスが分からないっていうことなのだが教えてもらえないのか？ なにを調べたらいいか分からないからここに来たんだが 369 名前：デフォルトの名無しさん [2008/01/20(日) 18:50:24 ] >>368 あなたはどういうデータを持っていますか？ 理論的には、どういう操作をしますか どういう式にあわせて、係数を推定したいですか？ 370 名前：デフォルトの名無しさん mailto:sage [2008/01/20(日) 19:08:25 ] わかった具体的に質問しよう まず>>357のような伝達関数で出力されたデータを>>360のような行列表記するとどうなるのか そもそもこのやり方以外に方法があるのならどういったものがあるのか ということをよければ教えてくれ 371 名前：デフォルトの名無しさん [2008/01/20(日) 20:11:17 ] 実空間で解いたらいいんじゃないの？ ２次の遅れ系なら解は２つの指数関数の和だろ？ それに入力パルス（今の場合e^(-τs) なのかな？）との叩き込み積分が掛かっているだけ。 それをコスト関数にしてfminsearchにかければいい。 行列表記ができるのはpolyfit関数なんかが使える線形な式のみじゃないの？ 372 名前：デフォルトの名無しさん [2008/01/21(月) 13:58:02 ] 実空間でfminsearch使ったら出来そうだけど、ζの大きさによって式が変わるから、 せっかくz空間で来まった式で表現できるG(z)がもったいないなぁと思ったんですが、そんなものですか？ あと、e^(-τ*s)っていうのはむだ時間を考えたときのです。いるかいらないかは実験データを見てから決めようと思うんですけど。 今functionの使い方を勉強してる段階ですが、fminsearchでやってみます。 373 名前：デフォルトの名無しさん mailto:sage [2008/01/21(月) 23:06:45 ] 370==372 なのか？ なにか教えてもらえると思うとすぐに下手に出るのか。 374 名前：372 [2008/01/21(月) 23:24:44 ] というよりイライラしていたら言葉が悪くなるんだよ 375 名前：デフォルトの名無しさん [2008/01/22(火) 21:13:48 ] 問題をあなたが理解していないのでは、と思わざるを得ない。 なぜなら、あなたが何をしたいのかが、伝わらないからである。 伝達関数の求め方はいろいろあるんだよ。 376 名前：デフォルトの名無しさん [2008/01/22(火) 23:54:33 ] 何をしたいか：伝達関数を出したい 教えて欲しいこと：伝達関数の求め方 求め方として最小二乗かなと思って>>353と書き込んだんだがどこが伝わらなかったんでしょうか？ 377 名前：デフォルトの名無しさん [2008/01/23(水) 09:50:12 ] >>376 それで伝わっていると思うのなら、 あんたの理解が甘い。 先生に教えてもらえ。 378 名前：デフォルトの名無しさん [2008/01/23(水) 13:29:48 ] linux(openSuSE 10.3)でmatlabを使おうとするとjavaのAssertionに引っかかって 死にます。java1.5だとそうなるらしいのですが、-nojvmというオプションを渡すと GUIなしで起動するという話もみました。でも_JAVA_OPTIONという環境変数に-nojvmを 渡しても、Matlabが起動しません。(matlabのactivationはまだできてない)。どなたか同様の問題を乗り越えた方はいませんか？ 379 名前：デフォルトの名無しさん mailto:sage [2008/01/23(水) 18:24:49 ] polarでMarkerSizeを変更する方法を教えてください 380 名前：デフォルトの名無しさん mailto:sage [2008/01/23(水) 20:08:24 ] ふつうのグラフで軸をスケーリングするには AXIS([XMIN XMAX YMIN YMAX]) とすればできますが、 subplotを使ったグラフの軸をスケーリングするには どうしたらいのでしょうか？ 381 名前：デフォルトの名無しさん mailto:sage [2008/01/23(水) 21:17:09 ] >>380 同じ 382 名前：デフォルトの名無しさん mailto:sage [2008/01/23(水) 21:51:14 ] >>381 同じでできるんですか？ AXIS([XMIN XMAX YMIN YMAX])としてもできなかったので 聞きにきたんですが。もう一度試してみます。 もしこれでできない場合、他に方法はないでしょうか？ 383 名前：デフォルトの名無しさん [2008/01/23(水) 23:35:12 ] >>382 subplot(211) plot(x,y) axis([xmin xmax ymin ymax]) % subplor(212) plot(x,y) axis([xmin xmax ymin ymax]) 384 名前：デフォルトの名無しさん mailto:sage [2008/01/24(木) 09:38:12 ] >>382 subplot(2,1,1); plot(x, y); set(gca, 'xlim', [xmin xmax]); set(gca, 'ylim', [ymin ymax]); でもいけるのでは？ 385 名前：デフォルトの名無しさん mailto:sage [2008/01/24(木) 13:55:38 ] >>381>>383-384 できました。ありがとうございました！ 自分が馬鹿でした。 axis設定してからplotしてました。 plotしてからaxisなんですね。 386 名前：デフォルトの名無しさん [2008/01/27(日) 13:14:09 ] >>384 set(gca, 'xlim', [xmin xmax], 'ylim', [ymin ymax]); で桶では？ 387 名前：デフォルトの名無しさん [2008/01/27(日) 18:20:12 ] 逆ラドン変換iradon()を使うとデフォルトでRampフィルタがかかってしまうのですが フィルタ無しの方法はありませんか？ 海外のMATLABのヘルプでは値でNo filteringにできるそうなんですがバージョンが 古いのかバージョン7.2では出来ません。 388 名前：デフォルトの名無しさん [2008/01/27(日) 20:46:38 ] 逆投影なんかツールボックスを使わなくても数行でかけるのに(´･ω･`) 389 名前：387 [2008/01/28(月) 04:32:33 ] なんとかヒントでもよければ教えてもらえないでしょうか？ 専門外のことなのでどうにもこうにもツールに頼ってしまってorz 390 名前：デフォルトの名無しさん mailto:sage [2008/01/28(月) 09:50:41 ] edit その関数 で、中身見てフィルタ部分カットしたのを書くのはどうですか？ 391 名前：デフォルトの名無しさん [2008/01/28(月) 12:49:16 ] まあ、簡単に考えると配列の方をirotateで回転させ、そこに加算していけば２行で書けるんだけどね。 それで角度でループさせるの。 iradonの中を触るんだったら、Hというのがフィルター関数なので、この部分を取り去るか、1にしてしまうか。 勉強のために１から作りたいんだったら、Ｃ言語で出ている教科書を買って、ネチネチと移植してもいい。 おれは院の先輩の薦めで「SPECT画像再構成の基礎」（医療科学社）という本の中に出ているフーリエ変換法のプログラムを移植したよ。 392 名前：デフォルトの名無しさん mailto:sage [2008/01/29(火) 00:08:17 ] グラフの大きさを変えるにはどうすればいいのでしょうか？ いつもマウスでウィンドウのサイズを変えて 縦横比を見やすく調整しているんですか コマンドでできないでしょうか？ 393 名前：デフォルトの名無しさん mailto:sage [2008/01/29(火) 02:07:23 ] position 394 名前：デフォルトの名無しさん mailto:sage [2008/01/29(火) 12:54:18 ] set(2,'position',[30 246 1120 420] って感じですね？ できました！ ありがとうございます 395 名前：デフォルトの名無しさん mailto:sage [2008/01/29(火) 15:41:21 ] figure(1); print -deps 〜.eps のようにグラフをepsファイルに出力してるんですが その出力したepsファイルのサイズ（縦横比）を変えることはできないでしょうか？ set 'position'でできるかと思ったらできませんでした。 396 名前：デフォルトの名無しさん [2008/01/29(火) 16:20:06 ] 自由落下した物体が床で跳ね返る時、速度の符号が 変わりますよね。その時ゼロクロッシングとかいう エラーでちゃいます。ゼロクロッシングの検知を 無効にするとうまくシミュレーションできません。 何でですか？ 397 名前：デフォルトの名無しさん mailto:sage [2008/01/29(火) 20:45:50 ] 483 名前：名無しさん＠５周年[] 投稿日：2008/01/29(火) 16:04:56 ID: 自由落下した物体が床で跳ね返る時、速度の符号が 変わりますよね。その時ゼロクロッシングとかいう エラーでちゃいます。ゼロクロッシングの検知を 無効にするとうまくシミュレーションできません。 何でですか？ 398 名前：デフォルトの名無しさん [2008/01/30(水) 16:50:48 ] set(gcf,'PaperPositionMode','auto')かな？ 399 名前：395 mailto:sage [2008/01/30(水) 20:04:51 ] >>398 printに-zbufferとか-openglオプションってのを つけてみてもできなくて困っていたのですが ずばりそれでできました！ありがとうございます 400 名前：デフォルトの名無しさん [2008/02/05(火) 13:58:03 ] MEXをVisualC++2008ExpressEditionで 開発するときの設定の仕方ってどこか にありませんでしょうか？m(_ _)m 401 名前：デフォルトの名無しさん [2008/02/05(火) 17:04:41 ] x=[0:1000]に対応するf(x)が手元にあって 各xについてx-x/10<=x <= x+x/10の範囲でｆの移動平均をとって平滑化をしたいんですが ループを使わない方法はないでしょうか？ 402 名前：デフォルトの名無しさん mailto:sage [2008/02/06(水) 21:10:33 ] 100*100*100のxyz空間を用意して,各点に対して同じ100*100*100のAという行列の値を カラーマップに対応する色で描画するにはどのようにすればいいのでしょうか？ X = 1:100;　Y = 1:100;　Z = 1:100; [XX,YY,ZZ] = meshgrid(X,Y,Z); によってxyz空間とそれぞれの点に値を持つ100*100*100のA行列は用意できている状態です。 例えばカラーマップの範囲が0(青)〜10(赤)でありA(50,50,50)=10の場合 出力されるfigureの(50,50,50)の部分を赤に色付けしたいです 403 名前：デフォルトの名無しさん [2008/02/08(金) 03:53:49 ] Control System Toolboxを使われている方に質問です。 ゼロ次ホールドの連続時間伝達関数である、 H = (1-exp(-sT))/s ;　　T:サンプリング周期[sec] という表現をmatlab上で再現するにはどのようにしたらよいのでしょうか？ s = tf('s'); を使って直接表現しようと試みたのですが、上手くいきませんでした。 連続時間システムの制御器の設計及びコンピュータへの実装の際に 使用したいと考えています。 よろしくお願いします。 404 名前：デフォルトの名無しさん [2008/02/08(金) 12:02:16 ] 構造体にセットしたデータをfwriteとかでバイナリファイルに書き込むのはどうしたら良いのでしょうか？ 405 名前：デフォルトの名無しさん [2008/02/12(火) 21:18:32 ] 素人質問ですが教えていただければ幸いです。 10個の自動計測されたデータが列方向に並んでいるcsvファイルを csvreadで一気に読み込んだ後に(例えばAとする)、 10個の列要素に小分け(A1〜A10)しようと考えています。 各要素の1行目は説明なので2行目から最終行までを選択して取り込むには どうすればよいのでしょうか? 思いつきでA1=A(:,2:inf)とかしようとしたらエラーがでてしまいました。 多分while文とかでNaNがでるまで読み込み続けるとかすると思ったのですが 力尽きてしまいました。 やさしい人。お願い。 406 名前：デフォルトの名無しさん [2008/02/12(火) 21:48:36 ] infじゃなくってendじゃないの？ 407 名前：デフォルトの名無しさん mailto:sage [2008/02/12(火) 22:15:32 ] >>405 よく読まずに書いてしまうが、infではなくendを使ってみたらどう

---

---

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

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