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

MATLABプログラミング ..

168:デフォルトの名無しさん
07/07/23 09:53:55
「何を」for文を使わずに'A.*B'や'A./B'の様にの様に書きたいのかを書いたら?
それは本当にlispのmap関数のようなものがないとできないのか?

169:デフォルトの名無しさん
07/07/23 10:17:18
>>168
おまえはまず興味があるのなら、Yes か No で答えてあげたらどうだ？
疑問文に疑問文で答えろと学校で教えているのか？ってやつだ

170:デフォルトの名無しさん
07/07/23 12:05:52
さいきんなぜだか知らんがこういうやつ多いよ

171:デフォルトの名無しさん
07/07/23 12:06:51
質問するなら何がやりたくて何が分からないのか全部最初に書く。
小出しにして質問者に質問させるような質問に親切に答えてくれる人は少ない。
>>169はそういう親切な人みたいだから、君が答えてあげたら?

172:デフォルトの名無しさん
07/07/23 12:11:23
>>171
あんたが質問の意味も理解できない馬鹿なだけだろ

173:デフォルトの名無しさん
07/07/23 12:13:49
mapでやるようなことをmatlabでは違う方法でやるのが普通。
>>172はちゃかしてるだけで結局答えることも
どう聞くべきかを教えてやることもしないんだろ?

174:デフォルトの名無しさん
07/07/23 13:10:04
普通じゃない方法はやるなと主張するわけだな？

175:デフォルトの名無しさん
07/07/23 14:57:43
>>174
そんなこと一言も言ってないのになんで勝手に決めるんだ?
でお前は質問に答えられもしないし
どう聞けばより答が得られやすいかを教えてあげるわけでもないんだろ?

176:初心者
07/07/23 21:36:26
MATLABのツールボックスを使わずに、
TSPの10都市問題をエネルギー関数を用いて解かせるプログラムを組みたいのですが、
全くわかりません。。
参考にできるサイトなどありましたらお教えいただけませんでしょうか。

177:デフォルトの名無しさん
07/07/23 21:41:48
pythonならmapなかった？

178:177
07/07/23 21:44:44
ちなみに
python matlab compatible
で検索すると幸せになれるかも

179:デフォルトの名無しさん
07/07/24 00:09:10
>>176
検索したら巡回セールスマン問題の解法アルゴリズムは
たくさん出てくると思うんだけど。

それをプログラムにしていく過程でまた具体的な質問が出てきたら
ここで質問したらいいと思うが、その質問だと全部ゼロから書いてくれ
って言ってるようなもんだぞ。

コツとしてはMATLABは忘れてpseudo codeでいいから求める手順を
まずしっかり書き出してみることだな。

>>177,178
matlabとmapで検索しまくって出てきた答?

180:デフォルトの名無しさん
07/07/25 02:09:20
>>176
traveling salesman annealingで検索したらそのまんまなのが見つかるな。
さすがに宿題でこのコード出すわけにはいかないだろうけど(笑)。

181:初心者
07/07/25 22:41:54
デフォルトの名無しさん、丁寧なご回答ありがとうございます。
参考にさせていただきます。

182:デフォルトの名無しさん
07/08/03 18:20:14
3つの列ベクトルx(i),y(j),z(k)をそれぞれx軸、y軸、z軸に設定してグラフを書きたいです。数式によって高さを与えるz平面のプロットは可能ですが実験で得られたデータを高さに与えて平面を書き出すことは可能ですか？

183:デフォルトの名無しさん
07/08/03 22:12:48
高さがあるのに平面とは…

184:デフォルトの名無しさん
07/08/05 04:35:53
mov2aviでPowerPointで再生できる動画を作成したいのですが、何を指定すればよいのでしょうか？

185:デフォルトの名無しさん
07/08/06 12:46:32
>>184
help mov2avi
何が出来て何が出来なかったのか言ってくれないと判らない

186:デフォルトの名無しさん
07/08/10 21:15:59
画像に、フィルターをかけるとき、フィルタのパラメータを調整したものを複数生成して
それらを一覧で表示したいのですが、
cat関数で横一列に繋げると、枚数が増えたとき表示のバランスが取れません。
何枚生成するかは、そのときそのときで違う場合で、よい方法はないでしょうか？

187:デフォルトの名無しさん
07/08/10 22:16:01
3次元のグラフというか立体図面を書きたい。
x,yの平面に高さzで、高さを面で現したような。
n行3列の行列があって(例えばdat)、行毎に[x,y,z]のデータの組が順不同に入っている。
例えば dat = [ 1 2 0.5; 1 3 0.6; 3 8 1.5; 1 1.5 0.6];
plot3(dat(:,1),dat(:,2),dat(:,3),'.')　では上手く行かないし、surfaceを使えば良いのか、検討が付かない。
何か手掛かりだけでも、誰か教えてください。

188:デフォルトの名無しさん
07/08/11 02:38:41
>>186
で、そのフィルターが1個のときと1000個のときで
具体的にどう表示して欲しいんだ?

>>187
よく分からんけどsurfは? とりあえずhelp graph3dしろ。

189:186
07/08/11 18:33:27
>>188
えっと、ある程度表示がよくなればまったく問題ないのですが
画像の右上に、パラメーターの数値を書き込んで
それらの画像をおおよそウィンドウ内部に
おおよそ正方形くらいになるように並べられらばと思っています。

190:デフォルトの名無しさん
07/08/12 00:37:31
1個のときの表示、2個のときの表示、・・・と全部具体的にイメージできてる?
それが出きてればそれをプログラムに書けばいいだけ。文章から推測するに、
君はやろうとしていることを具体的に細分化して理解しようとしていない。
「おおまかにこんな感じ」ではプログラムにならない。
コンピュータはそんなにかしこくないんです。

191:デフォルトの名無しさん
07/08/24 10:49:42
すごい初歩的なことなのですが少数点以下をもつ数の累乗根を出すような関数やコマンドは何でしょうか？

192:デフォルトの名無しさん
07/08/24 11:14:08
ヘルプは当然見たんだろうけど、どうやって調べた?

193:デフォルトの名無しさん
07/08/25 03:05:55
>>192
helpでコマンド一覧みたいなのを見たのですが整数の累乗根を出すものしか見つけられなかったので質問させてもらいました

194:デフォルトの名無しさん
07/08/25 04:42:51
>>193
うちでは整数のみのそのコマンドが見つからないです。

195:デフォルトの名無しさん
07/09/12 18:12:42
等間隔で行を抽出するにはどうすればいいですか？
例えば、1行目、11行目、21行目…という風にです。

196:デフォルトの名無しさん
07/09/13 01:12:14
1:10:end
これじゃダメ？

197:デフォルトの名無しさん
07/09/13 18:58:31
matlabで箱ヒゲ図を描きたいです。
ヒゲは、errorbarを使えば何とかなるのですが、
箱を描く方法が思いつきません。
どなたか、ご教授してください。

URLﾘﾝｸ(www.sci.kagoshima-u.ac.jp)
こういう図が描きたいのです。

198:デフォルトの名無しさん
07/09/13 19:20:26
株価のローソク足に似てるなーと思って「MATLAB ローソク足」でググったら一番上にそれらしきものがでてきたよ
Financial Time Series Toolboxがいるらしい

そこでCandle Plotを描くだけならだれかが自作のmファイルを後悔してるかも知れないと思って
「MATLAB Candle Plot」でググったらCandle.mを無性公開しているサイトを発見したぜ

199:デフォルトの名無しさん
07/09/13 19:20:30
>>195
fid = fopen(filename);
while 1
line = fgetl(fid);
　if length(line)==0
break;
end
// 任意の操作をここに。
for ii = 1:10 // 10行を読み飛ばす。
fgetl(fid);
end
end
fclose(fid);

でOK？

200:197
07/09/13 19:57:27
>>198
多謝。それです。
「ローソク足」というものがあるのですね。知りませんでした。
見つけていただいたコードを参考にしたいと思います。

以下、チラ裏。
plot( [x1 x2 x2 x1 x1], [uh uh lh lh uh] ));
で、箱が描けるとは。
matlab暦5ヶ月目で、新鮮な発見。
検索スキルとmatlab力を磨いて出直してきます。

201:195
07/09/15 16:14:20
>>196 >>199
できました、ありがとうございます。

202:デフォルトの名無しさん
07/09/23 16:27:08
MATLABには行列内部の最小の値を取る関数はあるけど
最小の値のある場所のインデックスを返す関数ってないよね

203:デフォルトの名無しさん
07/09/23 16:40:45
同じ関数

204:デフォルトの名無しさん
07/09/23 19:07:22
help その関数
で、わかりますよ＾＾

205:デフォルトの名無しさん
07/09/25 11:19:11
ygridの位置を指定するにはどうすればいいのでしょうか？

206:デフォルトの名無しさん
07/09/26 00:01:34
もっと詳しく書いて欲しいところだが、
set(gca, 'YTick', [ ... ]);
とかでどう？

207:デフォルトの名無しさん
07/09/26 16:05:17
matlabのグラフで，プロットした波形にマウスを当てると，
その時点での縦軸と横軸の値が出ますが，
それを，コマンドプロンプトでやりたいんですが，
関数かなにかありますか？

検索しても見つからなかったので，
やはり自分でプログラムを組むべきなのでしょうか？

208:デフォルトの名無しさん
07/09/27 15:34:32
>>207
つ[gnuplot]

209:207
07/09/27 23:44:56
>>208
ご教授ありがとうございます．

gnuplotは，さわりだけしか使ったことないんで，
どうなるか分かりませんが，やってみます．

210:207
07/09/27 23:47:12
ということは，
そういったことができるmatlabコマンドはないんですね・・・

211:デフォルトの名無しさん
07/09/28 14:52:23
画像処理でアフィン変換できる関数が見当たらないのですが気のせい？
結構基礎的な関数の気がするんだけど

212:デフォルトの名無しさん
07/09/28 22:45:39
そのくらい自分でやれ、ってことじゃないの？

213:デフォルトの名無しさん
07/09/29 11:00:49
まとめサイトはありませんか？

214:デフォルトの名無しさん
07/09/30 15:08:48
v=1*50000 vector & A=50000*50000 matrix

v*AをやろうとするとOut of memoryになっていまいます。

メモリーを追加しないで、効率よく計算するにはどうしたらいいですか？
よいBuilt-in関数とかないっすか？

215:デフォルトの名無しさん
07/09/30 19:10:11
>>202
ヘルプ見なさい

216:デフォルトの名無しさん
07/09/30 22:31:40
>203,204,215
helpみました
吊ってきます…

217:デフォルトの名無しさん
07/10/01 17:27:37
Axes上に複数のLineオブジェクトを作成しました。
このLineをマウス操作で移動したり、リサイズしたりしたいと思ってます…。

Axes上にLineを引くときは、function axes1_ButtonDownFcn( )でいけました。
しかし、function axes1_ButtonDownFcn( )の中にselectmoverisizeを書くと、
Axes自体が移動・リサイズされてしまいました。

そもそも、Lineオブジェクトをselectmoveresizeで操作するのは無理なんでしょうか？

218:デフォルトの名無しさん
07/10/02 15:25:07
3次元プロットで、（例えばsurfなどで）値がaxisで指定した範囲外にあってもプロットされてしまうのですが、
axisで指定した範囲内のみプロットするようにするにはどうしたらよいのでしょうか？

219:デフォルトの名無しさん
07/10/02 16:06:06
関数を定義する場合必ず別ファイルにしないといけないのですか？
同一ファイル内で定義することはできないものでしょうか？

220:デフォルトの名無しさん
07/10/02 18:45:41
>>218
aってのを3次元の行列だとして
surf(a(xmin:xmax,ymin:ymax,zmin:zmax))
じゃだめ？
意味取り違えてたらごめんね。

>>219
ｽﾚﾘﾝｸ(tech板:82-番)
同じスレなので出来れば読んでね。

221:デフォルトの名無しさん
07/10/04 03:19:13
matlabってリーもとデスクトップじゃつかえないのだね

222:デフォルトの名無しさん
07/10/11 14:00:14
bata版のNAG ToolBox for MATLAB
無料評価版

URLﾘﾝｸ(www.nag-j.co.jp)

223:デフォルトの名無しさん
07/10/11 14:23:25
pythonで大部分の数値計算できるけどmatlabってやっぱりまだ需要あるの？

224:デフォルトの名無しさん
07/10/11 16:29:25
>>223
C で python のできることは全部できるけど python って需要あんの？

225:デフォルトの名無しさん
07/10/12 00:19:45
俺紙とエンピツで全部やってるけどコンピュータって需要(ry

226:デフォルトの名無しさん
07/10/12 05:57:48
さすがにMATLABのスレでそりゃないだろ。
手作業でやったら一生かかっても終わらん計算してるだろうに。

227:デフォルトの名無しさん
07/10/23 14:25:52
質問なのですが、下記のプログラムをsimulinkでブロック図にて表すにはどうすればよろしいでしょうか？

clear allsyms z n p T
syms z n p T
h=0.1;
A=[exp(-h) int(exp(-T),0,h);-0.5 1];
b=[0;0.5];
c=[exp(-p) int(exp(-T),0,p)];
G=c/(z*eye(2)-A)*b;
R=z/(z-1);
Y=R*G;
y=iztrans(Y);
t=1;
for m=0:9,
v(1,m+1)=subs(y,n,m);
for l=0:9,
w(1,l+t)=subs(v(1,m+1),p,l*1e-2);
end
t=t+10;
end
figure
plot(w)
grid on

228:デフォルトの名無しさん
07/10/23 18:32:33
それは質問じゃなく変換依頼

229:デフォルトの名無しさん
07/10/25 03:43:18
タブ区切りのデータなのですが、一行目が見出しになっています。
一行目をのぞいて簡単に読み込む方法はないでしょうか？

230:デフォルトの名無しさん
07/10/25 07:53:38
fgetl

231:デフォルトの名無しさん
07/10/25 16:19:50
たぶん激しく初心者な質問でスマソなんですが、サイズがL,M,Nな3次元変数を3次元プロットする方法知りませんか？

232:デフォルトの名無しさん
07/10/25 22:13:34
強度イメージから一定範囲の強度値だけを抽出する
方法を教えてください

233:デフォルトの名無しさん
07/10/25 22:39:50
田中真紀子のﾎﾟﾛﾘ画像ありますか?

234:デフォルトの名無しさん
07/10/26 09:37:43
オナニーをすると身長が伸びないというのは本当ですか？

235:231
07/10/28 16:40:35
>>231
自己解決しますた

236:デフォルトの名無しさん
07/11/01 13:39:31
質問させてください。
specgramの逆の操作がしたいです。
X=specgram(x,128,8000);
で、
スペクトログラムXから音xを合成したいのですが、どうしたらいいですか？？

237:デフォルトの名無しさん
07/11/02 00:57:30
synthesizeとかいう命令があったような気がする。
matlabじゃなくてoctaveの命令だったかも。

238:デフォルトの名無しさん
07/11/02 06:26:24
Core2duoに変えたんだけど、６４BITに変えたら
計算すごく速くなる？
ベンチとかないかなぁ

239:デフォルトの名無しさん
07/11/02 19:57:08
入力したサウンドデータを，何の音なのか
(例えば，パトカーのサイレン，人の叫び声等)
分類するソフトを個人的に作成しています．

とりあえずmatlabでスペクトル表示を行い，
様々なサウンドデータのスペクトルを見たのですが，
分類の条件を何にしようか考えています．

「パトカーのサイレンは○○Hzで特に高いピーク値を得られる」
などの特徴を分類する際の条件にできるかなと思ったのですが，

この他に，音を分類する上で有効な条件となるようなものを
ご存知の方がいましたらご教授お願いします．

240:デフォルトの名無しさん
07/11/02 20:58:50
スレ違い

241:デフォルトの名無しさん
07/11/02 21:45:16
自分の目でいういろ比較して確かめなさい

242:デフォルトの名無しさん
07/11/04 12:40:22
はじめまして、matlab初心者です。
早速ですが質問させて下さい。

任意形状の２次元物体をメッシュ切りした上に、
コンターグラフの様な物を描かせたいのですが、
どなたか出来る方法をご存じないでしょうか？

宜しくお願いします。

243:デフォルトの名無しさん
07/11/04 19:30:34
matlabのメインウインドに打ち込むだけでsimulinkの開始ボタンを押したのと同じ処理になるコマンドってありますか？

244:デフォルトの名無しさん
07/11/05 20:40:09
>>243
SIM

245:デフォルトの名無しさん
07/11/07 10:02:25
例えば４ｘ４のマトリックスを２ｘ２にして、各要素は元々のヨウ素を加算してくれる関数はないですか？
ループとsum関数を使えばすぐに出来るんですがね

246:デフォルトの名無しさん
07/11/07 10:05:37
>>211
Image Processing Tool Boxだったかにあるよ。座標変換してからinterp2に流し込めばすぐにできるけどね。

247:デフォルトの名無しさん
07/11/07 22:42:27
>>245
じゃぁ、ループとsum関数を使えば？

248:デフォルトの名無しさん
07/11/08 01:47:43
>>247
matlab ではできるだけループを使わないのが重要
速度がダンチ

249:デフォルトの名無しさん
07/11/08 05:18:32
書き方が曖昧でわからんがエスパーの俺ならこうかな。Aのサイズは偶数として。
B=A(1:2:end,1:2,end)+A(2:2:end,1:2:end)+A(1:2:end,2:2:end)+A(2:2:end,2:2:end);

250:デフォルトの名無しさん
07/11/08 05:19:06
:と,一カ所ミスった。

251:デフォルトの名無しさん
07/11/09 15:59:58
??? エラー: ==> anfisedit
SCRIPT anfisを関数として実行してみてください.

これを解決するにはどうすればよい　のでしょうか？
一週間前ぐらいにMATLAB始めた者です。

252:デフォルトの名無しさん
07/11/09 20:08:08
>>251
SCRIPT anfisを関数として実行してみてください.

253:デフォルトの名無しさん
07/11/11 22:32:56
MATLABで式をニュートンラプソン法でとくにはどうしたらいいでしょうか？
分からなさすぎてノイローゼになりそうです。

254:デフォルトの名無しさん
07/11/12 02:23:50
まぁスレ違いだわな

255:デフォルトの名無しさん
07/11/12 07:39:22
GUI上でxの変域をかえたいんですけどどうすれば良いですか？

256:デフォルトの名無しさん
07/11/13 21:23:48
GUI上で数値を入力させて、
set(h, 'xlim', [x1 x2]);
じゃダメなの？

257:デフォルトの名無しさん
07/11/16 21:18:38
FORTRANで書いたものを、
MATLABに変換したいのですが、手書き以外で
何かいい方法あったら教えてください。

258:デフォルトの名無しさん
07/11/16 22:16:21
fftでパルスの周波数帯域を表示させたいんですが、
周波数軸をどうやって設定すればよいかわからないのです。
パルスはガウシアンパルスでUWBなので５ギガ位まで表示させるつもり何ですが、うまくいかないのです。
どなたかお願いします。

259:デフォルトの名無しさん
07/11/17 02:50:17
>>257
グーグルの使い方くらい知らないと
今後いろいろ苦労すると思うから覚えた方がいいよ。

260:デフォルトの名無しさん
07/11/17 18:16:24
ファイルの読み込みで指定した列だけ読む方法ってありましたっけ。

A = load("hoge.txt");
B = A(:,4);
みたいに、全部読んで目的列以外は捨てればいいんですが、
ファイルのサイズが大きいと、
全部読むのに時間がかかるので、効率化したいのです。

261:１８
07/11/17 23:13:55
matlabでシミュレーションを５回もするとメモリが５００Ｍから１Ｇくらいまで増えてしまうんですが、
これってなぜなんでしょうか？

262:デフォルトの名無しさん
07/11/18 00:14:17
>>260
dlmreadのマニュアル読むと良いよ。

263:デフォルトの名無しさん
07/11/18 13:49:03
対数データのフィッティングの方法を教えてください。
通常のデータであれば、最小二乗法ですが、対数データの場合はどのようにやればよいのでしょうか？

264:デフォルトの名無しさん
07/11/18 16:04:56
対数の底で冪をとって最小二乗法

265:beta-infinite
07/11/20 13:16:10
「MATLABの極」というブログを書いています。良ければ、ご参考に
URLﾘﾝｸ(beta-infinite.blogspot.com)

266:デフォルトの名無しさん
07/11/25 19:51:53
著作権違反ですよｗ

267:デフォルトの名無しさん
07/11/26 18:14:01
MATLABを使ってＢＰ法のニューラルネットワーク作らないといけないんですが
誰か請負いませんか？
価格は応相談、ＭＡＴＬＡＢ使えてＢＰ法分かる人ならそんなに難しい作業じゃないと思います
仕様とかその辺は請け負ってくれる人がいたらメールなりなんなりで連絡します

268:デフォルトの名無しさん
07/11/26 18:32:59
今、手一杯だからなー。
月８０万だったら考えるけどね。
それかNNツールキット買えばいいんじゃないの？

269:デフォルトの名無しさん
07/11/26 18:46:20
ていうか自分で作れって感じ
層やニューロン数を完全に一般化した物でなければ速攻で作れるでしょ
もちろんfor文抜きでな

270:267
07/11/26 19:12:50
料金は
３時間未満：一律12000円
３時間以上：時給4000円
７時間以上：一律28000円
このあたりでよろ

おれ自身じゃ作れないので金で解決する

271:デフォルトの名無しさん
07/11/26 19:19:22
安すぎ

272:デフォルトの名無しさん
07/11/26 19:29:20
委託研究でMATLABでソフト作ったときは、３０万もらったな

273:デフォルトの名無しさん
07/11/26 20:31:27
作成物の複雑さとか教えてもらえずに
7時間以上28000円は安いな。

274:267
07/11/27 00:00:36
ここでは話せないが詳しい仕様を聞いてくれるやつはレスくれ
秘密厳守の上その時点で断ってくれてもＯＫだから

275:デフォルトの名無しさん
07/11/27 00:06:40
こんなところで依頼しようとしている時点でダーティなイメージがついてて嫌だ
全然予算がないようだし学生の卒論か？

276:267
07/11/27 01:14:15
>>275
まさにそれ
誰か頼むよ～

277:デフォルトの名無しさん
07/11/27 01:19:35
URLﾘﾝｸ(www.mes.p.dendai.ac.jp)
↑この最後のページにある問題を教えてください。
参考の為に問題文だけ書くと・・・

１）40名のクラスがあった時、すべての人の誕生日が異なる確立をモンテカルロ法で計算しなさい。
２）その結果を解析的に求めた結果と比較・考察しなさい。
よろしくお願いします。

モンテカルロ法が使えるプログラムはMATLABと書いてあったのでここに質問しました＞＜
スレ違いでしたらすみません。。

278:デフォルトの名無しさん
07/11/27 06:40:06
別に matlab じゃないとできないわけじゃない。プログラミング言語に依存してたまるか

279:請負人
07/11/27 11:53:45
請負内容を見て判断したいです。詳細内容などをメールをください

280:請負人
07/11/27 11:59:51
>>267
>>270
請負内容を見て判断したいです。詳細内容などをメールをください。
秘密厳守を守ります。

281:デフォルトの名無しさん
07/11/27 18:59:28
手元にはバックプロパゲーション法の自作ソースくらいあるんだが
バカ学生に使いこなせるとは思えん

282:デフォルトの名無しさん
07/11/28 15:25:32
matlab にファイル名を変えるコマンドってありましたっけ？
同時にファイル内の関数名もあわせて変えてくれるとベスト

283:デフォルトの名無しさん
07/11/29 23:23:21
URLﾘﾝｸ(www.cybernet.co.jp)

このページの
h=plot(0,0,'.b','MarkerSize',30,'EraseMode','xor');
を使ってアニメーションを作る方法で、
同時に2つの色の違うボールを描画するにはどうすればいいのですか？

284:デフォルトの名無しさん
07/11/29 23:34:51
matlabで書いたグラフで、線と線に囲まれた部分の
面積を求めることは可能でしょうか。
ファイルをアップしましたので、よろしかったら見てみてください。
URLﾘﾝｸ(www.csync.net)
(ファイル名はmatlabです)

画像の灰色の部分の面積を求めたいのですが、
積分のやり方もわからないのです(:_:)

あとMファイル（txtになってしまってます）のコピーが二つありますが、
mensekiの方は画像の元になるグラフが書けます。
danmenの方は、元のプログラムですが、
xの値を０～０．５まで変えると赤と緑の線が一周するように
できているのですが、この二つの線と、黄色の線が作る面積を求めたいんです。

どんな方法でもいいのでどなたか教えてくださいm(__)m

285:beta-infinite
07/11/30 00:29:49
>>283
「MATLABの極」にソースをアップロードしました。いかがでしょうか
URLﾘﾝｸ(beta-infinite.blogspot.com)

286:デフォルトの名無しさん
07/11/30 03:26:05
教えてください。
MAT LABで書いたFigureのX軸Y軸をX→Y、Y→Xと逆にしたいのですが
どうすればいいんでしょうか？お願いします。

287:283
07/11/30 10:53:34
>>285
ありがとうございます。
今までhold onとか使っていたんですが、こんなシンプルにいけたんですね。

ですが、シミュレーションしたものをアニメーションとして表示したいので
この方法だとどうしても描画が遅くてシミュレーションの結果通りに描画できません。
そこで>>283の方法は高速だということなので、
>>283のイレースモードxorを使う方法でのやり方を知りたいのです。

>>283の方法以外でも高速に描画できる方法をご存知の方はお願いします。

具体的にしたいことは、
・2つボールの0.001秒毎の位置が2つの配列に入っている。
・それをある時間dt秒毎に描画したい。
・1回の描画毎にpause(0.001*dt);などして時間に正確にアニメーション化したい。

自分のやり方や>>285さんのやり方だと
dtを大きくすればするほどアニメーションがゆっくりになってしまいます

288:283
07/11/30 10:56:39
>>287の修正

×　pause(0.001*dt);
○　pause(dt);

289:MATLABの極
07/11/30 12:02:44
>>288
PAUSE(n) は、n秒間停止して、その後実行を続けます
当然ですが、pause(dt)のdtを大きくすれば、
アニメーションがゆっくりになってしまいます

290:MATLABの極
07/11/30 13:29:41
>>286
view(90,90)という関数を使えば、いいと思います。
URLﾘﾝｸ(beta-infinite.blogspot.com)
に参考

291:283
07/11/30 13:44:06
>>289
ですが、dtを大きくするとその分荒いサンプリングをすることになるので
描画時間がdt秒に比べて十分小さければ同じようなアニメになるはずですよね
コマ送りになるかどうかの違いはあれど。
でもやはり描画時間はそこまで小さくするのは無理なのでしょうか。
MATLABではこんなことはできないのでしょうか

292:MATLABの極
07/11/30 14:04:16
>>291
現状の１コマの処理時間は、
= matlabの演算時間 + 描画の時間 + pause(dt)の時間
matlabの演算時間 + 描画の時間をなくする方法は、getframe,movie,movie2avi関数で
使うといいと思います。

URLﾘﾝｸ(beta-infinite.blogspot.com)
に参考

293:デフォルトの名無しさん
07/11/30 15:21:36
beta-infinite.blogspot.comってパクリしかないじゃんかｗ

294:MATLABの極
07/11/30 15:45:10
>>293
MATLABの修業中です。覚え書きのブログだと思っていただければ、

295:283
07/11/30 16:30:08
>>292
やっぱgetframe,movie,movieしかないですかね
ありがとうございました

296:デフォルトの名無しさん
07/11/30 17:13:21
著作権違反だので、通報しました

297:デフォルトの名無しさん
07/11/30 17:57:24
FFTを使って、スペクトルを計算するのって、どうすればいいの？

298:デフォルトの名無しさん
07/11/30 21:40:02
だので

299:デフォルトの名無しさん
07/11/30 22:31:07
>>297
help fftとコマンドを打って出てきた説明を理解できなければ
回線を切って死ぬか教授に頭を下げて教えを請いなさい

300:デフォルトの名無しさん
07/12/01 04:36:12
>>297
beta-infinite.blogspot.comだけは見ないように。

301:デフォルトの名無しさん
07/12/01 12:59:02
fftはtwo-sideの出力を出すので、help fftだけでは話は完結しない

302:デフォルトの名無しさん
07/12/01 16:02:42
教えてください。
プログラムで図を描かせたときにX,Y,Zラベルの刻みのラベルを
変更したいのですが、関数などプログラムで書かせることができるでしょうか？
少し伝わりにくい文章ですが、よろしくお願いします。

303:デフォルトの名無しさん
07/12/01 16:26:01
なんか楽に関数名リネームする方法ねーのかよ

304:デフォルトの名無しさん
07/12/01 17:18:44
>>303
関数名リネームって>>282が言ってるようなことか?
もしそうならファイル名と違ったらファイル名の方が優先なだけだから
コード中の関数名は適当でOK。

しかしそういう時は普通呼ぶ側のコードも変えないといけないから
どっちにしろ一括置換とかする程度の手間をはぶくのは無理なわけで
コードの頭の関数名変えるくらいの手間は取るに足らないと思うが。

305:デフォルトの名無しさん
07/12/01 19:52:21
>>297
ここ見るといいよ
URLﾘﾝｸ(www.mathworks.com)

306:デフォルトの名無しさん
07/12/02 18:39:33
>>302
できる。

307:デフォルトの名無しさん
07/12/04 06:02:50
関数の入力にベルトルの集合としての行列をいれるとき皆は
NxD 派？ DxN 派？where N is the number of vectors and D is the dimension of each vector.
どっちのほうが親切なんだろう？toolbox によって（人によって）違うんだよねぇ。

308:デフォルトの名無しさん
07/12/04 09:55:13
シミュレート板で回答を得られなかったので。

simulinkである変数がある値になったときシミュレーションを終了したいのですが
どうすればできますか？

309:デフォルトの名無しさん
07/12/05 06:53:25
ムービーを作成する際に、figureウインドウを出さずに作成することはできないでしょうか？
figureウインドウを作成してgetframeで取り込む場合、スクリーンキャプチャをしているようで、
他の作業をすると、その画像がムービーの中に入ってしまいます。
そのためムービー作成中は何も作業ができません。

310:デフォルトの名無しさん
07/12/05 16:25:53
DICOMファイルから画像読み込み、JPG形式で書き込みしようとしてるのですが、
X=dicomread('--.dcm');
[Y,map] = gray2ind(X);
imwrite(Y,map,'--.jpg','jpg');
とすると、colormapを使用したimshowの表示画像と異なってしまいます。
よろしくお願いします。

311:デフォルトの名無しさん
07/12/05 18:06:21
正の実数からなる対角行列の、各行と列について、
要素の大きさが対角成分の要素の大きさの0.5から2倍の範囲なら１で置換、
それ以外なら０で置換する、
ってforループを使わずにできるでしょうか？

312:デフォルトの名無しさん
07/12/05 23:19:50
質問ばっかで誰も答えてくれない罠ｗｗｗｗｗｗｗｗｗ

313:デフォルトの名無しさん
07/12/06 00:59:41
>>311
対角行列なら対角成分以外すでに全部０やん

314:デフォルトの名無しさん
07/12/06 02:21:48
これぞまさに質問スレ

315:デフォルトの名無しさん
07/12/06 08:48:34
vistaは2007aで正式サポートですが、2006aで動かした場合、どこかに既知の問題のようなリストはないでしょうか？

316:デフォルトの名無しさん
07/12/07 01:04:38
>>307
僕はMatlabで音系のことをやってるけど、
wavwriteとかのデフォルトが行＝サンプル／列＝チャンネルなので、
ベクトルの集合はいつも縦長ベクトルを横にいくつか並べるというかんじ。
DxNってことになるのかな？

317:デフォルトの名無しさん
07/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:デフォルトの名無しさん
07/12/12 05:19:53
秀丸用の強調定義ファイルはありませんか？

319:デフォルトの名無しさん
07/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:デフォルトの名無しさん
07/12/12 10:36:23
fgetlで１行まるまる読んで、strtokで切り出していくしか無いんじゃないの？Ｃでも普通そうだろ？

321:デフォルトの名無しさん
07/12/13 04:09:16
>>319
それを一発で君の思い通りに読んでくれる
ファイル読み込み関連の関数が貧弱じゃないやつって何？
煽ってるんじゃなくて単に興味があるから知りたい。
ExcelとかRとかできるんかな?

322:デフォルトの名無しさん
07/12/13 12:21:50
excelじゃ無理だね。

323:デフォルトの名無しさん
07/12/13 13:11:34
figureウインドウを表示させずにムービーを作る方法を教えてください

324:デフォルトの名無しさん
07/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:デフォルトの名無しさん
07/12/13 14:21:09
>>324
textread は一行ずつ読み込むんだろ？
最初の行は先頭以降何もないから０をが入っていくと。
[data1,data2,data3,data4,data5,data6,] = textread('ファイル名.txt', '%f %f %f %f %f %f')
として、奇数行と、偶数行に後でわけるとかになるんじゃね？

326:デフォルトの名無しさん
07/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:デフォルトの名無しさん
07/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:デフォルトの名無しさん
07/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:デフォルトの名無しさん
07/12/13 15:42:56
どうやらうまく書き込めないのですが
１番目、２番目の１の下にはスペースが存在していて
３番目の0.00000000E+00 と下の行の0.00000000E+00 は
列がそろっています。

330:ks
07/12/16 14:41:45
＝＝＝＝＝＝＝＝
教えてください。
＝＝＝＝＝＝＝＝
１次遅れ＋無駄時間要素を含む制御対象の限界感度法による
PIDパラメータをmatlab上で、どう求めるんですかね～

＝＝＝＝＝＝＝＝
教えてください。
＝＝＝＝＝＝＝＝

331:デフォルトの名無しさん
07/12/21 01:33:00
>>329
textreadにこだわらなくてもいいんじゃないの？

332:デフォルトの名無しさん
07/12/27 03:52:31
matlabでfor文をまわす時
ループが終わるごとに変数をclearした方がよかったりします？

不振な挙動。。

333:デフォルトの名無しさん
07/12/27 09:04:36
そんなことないと思うが。
それにしても、Matlabバージョンとか書かないと。

334:デフォルトの名無しさん
07/12/28 00:55:49
たぶんそれはいくつかの変数を間違って使ってるかロジック間違えてるだけだ。
ちゃんと全部読み直せ。

335:デフォルトの名無しさん
08/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:デフォルトの名無しさん
08/01/08 00:30:09
2値化した画像をカラー画像にするにはどうしたらいいですか

337:デフォルトの名無しさん
08/01/08 00:51:57
それはもう情報がかなりないからとても難しいあるね

338:336
08/01/10 14:39:29
>>337
カラーは白黒でよいのですがそれでも無理ですか?

339:デフォルトの名無しさん
08/01/10 16:15:22
>>338
それは、フォーマット変換したいと言うこと?
それとも、2値のメッシュをグレイにしたいと言うこと?

340:デフォルトの名無しさん
08/01/10 22:09:11
2次元座標から、各座標の点をラインで連結してプロットさせ
それを画像として処理したいのですが、現在は
plot(point(:,1),point(:,2));または
line(point(:,1),point(:,2));
pointは100*2行列、xy座標の点100個としています。
これだとプロット結果を保存しても画像のサイズが変更されてしまいます。
点と点を直線で結んでいく関数などありましたら教えていただけないでしょうか？

341:デフォルトの名無しさん
08/01/11 04:51:05
行列の積で初歩の初歩的な質問だと思うのですが誰か教えてください！
10000x10の行列Aがあるとします。
A*A'だとout of memoryと出力されるのですが
A'*Aだとちゃんと計算結果が出力されます。
考えてもイマイチ理由が浮かんで来なかったのでどなたか教えてください～！
単純な質問で失礼しました。

342:デフォルトの名無しさん
08/01/11 07:34:46
前者は 10000x10000
後者は 10x10
の出力になるだろう？明らかすぐる

343:デフォルトの名無しさん
08/01/11 08:40:10
これはひどい

344:デフォルトの名無しさん
08/01/11 09:11:51
（・A・）＊（・A・）'

345:デフォルトの名無しさん
08/01/11 15:54:33
ホントだｗｗｗ頭悪くてすまそｗｗ
ありです。

346:デフォルトの名無しさん
08/01/12 16:46:00
plotで描ける線の種類を増やせないでしょうか？
6種類のデータを1枚の図に、白黒で書く、ということで困っています。

347:デフォルトの名無しさん
08/01/12 21:13:17
>>346
解決しないかも知れないけどマーカと併用するのはダメなんですか
線種で区別できる程に線が絡み合ってないのならそれでいけるのではないでしょうか

348:デフォルトの名無しさん
08/01/12 22:26:34
>>339
フォーマット変換のほうだと思います。

349:デフォルトの名無しさん
08/01/15 10:21:20
>>348
MATLABでなくていいなら、画像ビューワでも画像エディタでもいくらでもやりようがありそうだが。
何を、何故、MATLABでやりたいのか全然判らんのだよ。

350:デフォルトの名無しさん
08/01/15 22:46:03
>>349
MATLABで画像処理の勉強をしろといわれ、つい2値化の画像を白黒だけの
カラー画像に変えてやろうと思った。

351:デフォルトの名無しさん
08/01/16 00:07:40
まずおめーは、白黒だけのカラー画像とかわけわかんねーことをいうな
カラー画像っつったらＲＧＢとか３チャネルあるやつを普通は指すんだよ。

352:デフォルトの名無しさん
08/01/16 00:15:23
あほす
R　==　G　==　B　のことだろ

それより2値化の画像って何だ？
2値の画像だろ

353:デフォルトの名無しさん
08/01/17 06:03:13
ステップ応答の試験をしてMATLABで伝達関数を出したいんですけど、
最小二乗でフィットさせるときの方法とかいうのが分からないので教えてください。

354:デフォルトの名無しさん
08/01/17 12:20:28
テキストファイルの行数を事前に知る関数はないでしょうか？

355:デフォルトの名無しさん
08/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:デフォルトの名無しさん
08/01/18 00:32:14
>>353
どういう関数系にフィットさせるの？

357:デフォルトの名無しさん
08/01/18 06:14:54
>>357
よくある１次遅れ、2次遅れ、むだ時間とかいうパラメータを使った低次元伝達関数です
たとえば
G(s) = Kωn^2 / (s^2 + 2ζωns + ωn^2) * e^(-τs)
のパラメータを出したいんです

358:デフォルトの名無しさん
08/01/18 06:15:39
安価ミス
×　>>357
○　>>356

359:デフォルトの名無しさん
08/01/18 06:31:51
あ、あと私のMATLABにはSystem Identification Toolboxが入っていないのでこれを使わない方法を教えて欲しいです。

360:デフォルトの名無しさん
08/01/18 08:58:30
そこまでわかっているのなら、最小自乗でフィットするだけじゃん

URLﾘﾝｸ(www.wakayama-u.ac.jp)

361:デフォルトの名無しさん
08/01/18 10:53:24
指数関数的に減衰していくデータにフィットさせたいのですが、このように非線形の場合どのようにやればよいのでしょうか？
対数プロットでは直線になるので、何か最小自乗法のような物は使えないでしょうか？

362:デフォルトの名無しさん
08/01/18 17:47:29
fminsearch関数使えば？Simplex方による繰り返し計算だけど。

363:デフォルトの名無しさん
08/01/18 19:14:34
数学的には間違いだけれども、
対数変換して直線フィットすると、ある値が出る。
それを初期値に>>362にしたがって、解を探す。

364:デフォルトの名無しさん
08/01/19 12:02:15
>>360
ごめんなさい
その最小二乗でだけっていうところが分からないでの聞きにきたんですよね・・・

365:デフォルトの名無しさん
08/01/19 12:47:42
polyfit使えば一発。ただ原点に返って正規方程式を逆行列で解くのもありだな。

366:デフォルトの名無しさん
08/01/19 14:03:31
>>364
MATLAB以前の問題なのか…

>>365
polyfitって、フィットさせる式に制限があるからね。

367:デフォルトの名無しさん
08/01/19 21:23:54
>>365
ployfitって単なる多項式しか出来ないんですよね

>>366
MATLAB以前の問題かもしれないんですけど、
やるとしたらデジタルで出てきたデータをz変換してそれを行列使って最小二乗してやるということになるのでしょうか？

368:デフォルトの名無しさん
08/01/20 15:18:25
だから単純なたとえばy = a*x^2 + …とか>>360に書いてあることはわかるんだが、
伝達関数の場合のプロセスが分からないっていうことなのだが教えてもらえないのか？
なにを調べたらいいか分からないからここに来たんだが

369:デフォルトの名無しさん
08/01/20 18:50:24
>>368
あなたはどういうデータを持っていますか？
理論的には、どういう操作をしますか
どういう式にあわせて、係数を推定したいですか？

370:デフォルトの名無しさん
08/01/20 19:08:25
わかった具体的に質問しよう
まず>>357のような伝達関数で出力されたデータを>>360のような行列表記するとどうなるのか
そもそもこのやり方以外に方法があるのならどういったものがあるのか
ということをよければ教えてくれ

371:デフォルトの名無しさん
08/01/20 20:11:17
実空間で解いたらいいんじゃないの？
２次の遅れ系なら解は２つの指数関数の和だろ？
それに入力パルス（今の場合e^(-τs) なのかな？）との叩き込み積分が掛かっているだけ。
それをコスト関数にしてfminsearchにかければいい。
行列表記ができるのはpolyfit関数なんかが使える線形な式のみじゃないの？

372:デフォルトの名無しさん
08/01/21 13:58:02
実空間でfminsearch使ったら出来そうだけど、ζの大きさによって式が変わるから、
せっかくz空間で来まった式で表現できるG(z)がもったいないなぁと思ったんですが、そんなものですか？
あと、e^(-τ*s)っていうのはむだ時間を考えたときのです。いるかいらないかは実験データを見てから決めようと思うんですけど。

今functionの使い方を勉強してる段階ですが、fminsearchでやってみます。

373:デフォルトの名無しさん
08/01/21 23:06:45
370==372
なのか？

なにか教えてもらえると思うとすぐに下手に出るのか。

374:372
08/01/21 23:24:44
というよりイライラしていたら言葉が悪くなるんだよ

375:デフォルトの名無しさん
08/01/22 21:13:48
問題をあなたが理解していないのでは、と思わざるを得ない。
なぜなら、あなたが何をしたいのかが、伝わらないからである。

伝達関数の求め方はいろいろあるんだよ。

376:デフォルトの名無しさん
08/01/22 23:54:33
何をしたいか：伝達関数を出したい
教えて欲しいこと：伝達関数の求め方

求め方として最小二乗かなと思って>>353と書き込んだんだがどこが伝わらなかったんでしょうか？

377:デフォルトの名無しさん
08/01/23 09:50:12
>>376
それで伝わっていると思うのなら、
あんたの理解が甘い。

先生に教えてもらえ。

378:デフォルトの名無しさん
08/01/23 13:29:48
linux(openSuSE 10.3)でmatlabを使おうとするとjavaのAssertionに引っかかって
死にます。java1.5だとそうなるらしいのですが、-nojvmというオプションを渡すと
GUIなしで起動するという話もみました。でも_JAVA_OPTIONという環境変数に-nojvmを
渡しても、Matlabが起動しません。(matlabのactivationはまだできてない)。どなたか同様の問題を乗り越えた方はいませんか？

379:デフォルトの名無しさん
08/01/23 18:24:49
polarでMarkerSizeを変更する方法を教えてください

380:デフォルトの名無しさん
08/01/23 20:08:24
ふつうのグラフで軸をスケーリングするには
AXIS([XMIN XMAX YMIN YMAX])
とすればできますが、
subplotを使ったグラフの軸をスケーリングするには
どうしたらいのでしょうか？

381:デフォルトの名無しさん
08/01/23 21:17:09
>>380
同じ

382:デフォルトの名無しさん
08/01/23 21:51:14
>>381
同じでできるんですか？
AXIS([XMIN XMAX YMIN YMAX])としてもできなかったので
聞きにきたんですが。もう一度試してみます。

もしこれでできない場合、他に方法はないでしょうか？

次ページ