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

こんばんわ
早速ですが質問させてください
function のなかで function を定義するのはどうすれば良いのでしょうか。

hoge.m

function hoge
foo
end

function foo
end

>>83
成程、関数の終了時にendが必要なのですね。
早速のご回答、ありがとうございました。

なくても動いたと思うが。

リモートデスクトップ環境ではmatlabは実行できないのでしょうか？？？？

おはようございます
早速ですが，1*pの行ベクトルxがあるときに，
xを縦にn行並べたn*p行列はfor文を使わないで生成できますでしょうか．
よろしくお願いします．

教えてください
Optimization Toolboxのfminuncで
関数y=func(x,z)をxについて最小化したいんですが，
どうすればよいでしょうか．zはこちらで与えたい定数です．
ヘルプには
X = fminunc(FUN,X0,OPTIONS,P1,P2,...) は、
問題に関連したパラメータ P1,P2,... を、直接、関数
FUN に渡します。
と書いてあったんで
x=fminunc(@func,x0,[],z)
というように，zを与えようとしたのですが，うまく行きません．
ご教示いただきますと幸いです．
よろしくお願いいたします．

>>87
repmat

>>87
ones(n, 1) * x;

複素成分の符号を逆（+なら-、-なら+）にする関数ってありますか？

間違えました。複素数の虚数部でした。

conj

最近は共役って言葉を数学で教えなくなったんだな

円周率は大体３、という教育に期待するほうが無駄

わけも分からず３．１４と教えられる方が問題あると思うよ
現代の小学生は正方形の中に
円をかいて４より小さいって自分で調べて３をみつけだしているからね

ご自分のときはどうだったのかな？
私のときは方眼紙に缶とかを利用して円を描き、
その面積を計算したりしたよ。

円を描いて32等分。
並べ替えて平行四辺形みたいな形状をつくり、縦ｘ横。
更に班内で平均。

構造体を与えると固有のパスコードが返されるような関数ってありますか？

例えば
a.b = 1;
a.c = 'test';
なときに
x = somefunc(a);
で
x = 'ZSGADRERASVFDSAD'
みたいなのが返ってくるような。

入力する構造体の構造と値が一致した場合のみ戻り値xも一致して、
そうでない場合は全然違う値が出てくる、というようなものを探して
います。何か良い方法を知っている方がいましたら教えてください。

MATLABやScilabで式長の制限はあるんでしょうか？
Scilabでかなり長い式を定義しようとしたら，
invalid factorと言われて困っております．．

>>99
そのものずばりなものは見たことないよ。
fieldnames()とかを使ってフィールド名を出してきたり、struct2cellで値を出してきたりして、あとはMD5なんかに突っ込むしかないんじゃないかなー。

初心者のしつもんです
simulinkで
y(0)=0
y'=1-xy
の微分方程式を作るにはどうやるのか教えてください
やってみたけどわかりませんでした
よろしくお願いします

すみませんスレ間違いました
出直してきます

ソースにIsim、placeってあるんですがどういう意味ですか？

あげ

isocolorsで3次元上の物体の表面に色をつけることができますが、同様に、alpha値を指定することはできないでしょうか？

一からMATLAB/simulinkを勉強しようと思っているのですが、
どういったところから始めるべきでしょうか？
解説書を見ながら、何かのデータをいじくっていけばいいのでしょうか？

clearで変数をクリアするとき、指定した変数以外をクリアする方法が知りたいです。
-regexpで正規表現を使うってのはわかるんですが
a ab abcって変数があるとき
clear -regexp [^\<ab\>]
とするとabcだけがクリアされちゃいます。

すいません変数abだけクリアしたいのです。

またまたすいません。
変数ab以外です。

１０８ですけど誰か答えられる人いないんですか?

clear allですべて削除した後、必要なものだけ再作成。

正規表現の話だろ？もはや matlab も関係ないよ。正規表現を勉強しろ、とだけ。

>>108
まあそれだとabcにマッチするよな当然。
そしてbはclearされないけどxはクリアされると。

リモートデスクトップでmatlabは使えないのでしょうか？

割れものだと使えません。

自宅で試行錯誤するだけならFreeMATでもいいんじゃないかな
ttp://freemat.sourceforge.net/wiki/index.php/Main_Page

場合によっては octave-forge から持ってくるにしても
関数群が貧弱であれだけど

>>116
そうなの？

初心者です。質問させてください。
imhistでプロットしたヒストグラムの数値を
計算に使うにはどうしたらいいんでしょうか？
どこかに格納されてたりするんでしょうか。
お願いします。

ヘルプをみれば返り値として受け取れると書いてある

ありがとうございます〜。

3Dボリュームの可視化で、表面の色を液体金属のようにしたいのですが、
どのようなパラメータを設定すればよいのでしょうか？

gzipで圧縮されているデータを解凍してから表示しているのですが、
解凍作業をmatlab上から行えないでしょうか？

!でコマンドラインよべなかったっけ？

>>124
dosとかunix

matlabで3次元グラフを表示させてその表示した三次元グラフの表面積って計算可能ですか？

表面積の計算は自分でやってください。

関数外のグローバルな変数を関数内で用いるには、変数として代入するしかないのでしょうか？

function Y = f(X,Za,Zb,･･･)

if X=1
　Y=sum(Za)
elseif X=2
　Y=sum(Zb)
elseif...

というように必要なグローバル変数Z*が非常に多いので、代入以外に方法があればと考えています。
Z*を一つの変数にまとめることも考えましたが、すでにZ*それぞれが３次元変数なので、出来れば避けたいです。

>>128
matlab global でぐぐるくらいしようよ。一発だよ。

x -y
y x
の逆行列を求めるにはどうやるの？

くだらない質問でごめんなさい。
色々調べてみてもイマイチ理解が出来ていないんですが･･･。
matlabでローパスフィルタを作製しようと思ったときに参考にできるサイト
ってありますかね？

>>130
inv([x -y ; y x])

>>131
ローパスフィルタを作ること自体はfilter関数に何を入れてやればいいかだけ。
信号処理論の勉強をすべし。

hann関数を使用すると
未定義の関数、またはメソッド 'hann' (タイプ'double' の入力引数)。
と表示されます
w = hann(1024);
と使っているのですがうまく動きません
だれかご教授を

Signal Processing Toolboxが入ってないんじゃない?

多分それだと思います
すぐに確認できないんで明日確認して報告します

確認したらSignal Processing Toolboxがありませんでしたので
練習ついでにhanning関数作りました。
ありがとうございました。

hannning関数をつくるなんてすごいなぁ

その辺の関数ならOctaveのがそのまま使えるんじゃないの?

2次元データを平滑化する関数はありませんか？

forループで時間を進めて
plotで表示しているのですが、最後のplotしか表示されません。
途中経過も表示したいのですがどのようにしたらよいのか教えてください

figure
hold on
for i=1:n
plot(なんとか)
end
hold off

こんな感じ？
forの前にhold onしてforのあとにhold offかな？

あー、もしかしてこっちかな？
figure
for i=1:n
plot(なんとか)
drawnow
end

drawnowですぐ結果表示だからね。
質問をはっきりして欲しい。

>>141
ソースも書かずになんで原因が分かると思ったんだ?

>>142,143 はひまでしょうがないから
たまたまあてずっぽうで答えてくれただけだからな。

>>143でうまくいきました。
どうもありがとう。

Linux版で動画を作成したいのですが、movie2aviは非圧縮しか対応していません。
圧縮フォーマットは何でもよいのですが、何かよい方法はないでしょうか？

>>146
mplayer とか ffmpeg とかいろいろあるじゃん。

教授から、作成したmファイルがコンパイルできないで困っていると相談を受けました
ソースを解読しながら修正することになったのですが、如何がコンパイルできない、どの関数が実行できないなどの情報はどこかに集まっていないのでしょうか？

教授のくせに頭悪いな

せめてどんなエラーが出るか書くとかしてね。

ふーむ m ファイルをコンパイルか・・・

「移動平均」ならぬ「移動最大値」（or「移動最小値」）
を計算できるような方法はございますでしょうか？
（自分でforとmax使って作ると遅いので・・・）

>>149
どうもすいません。確かに前提のお話が足りていませんね。

>>150
言いたいことは分かりますが、ただ教授も神様じゃないんです。

一応教授から聞いた症状ですが
・使用したコンパイラはmcc
・環境はmatlab6.5でコンパイルできた
・matlab7.1にしたらコンパイルできなくなっていた（これは共同研究者の環境ら

しいです

と、ここまでは聞いた話です。忘れてください。
それで、自分で試した限りでは
・matlab6.5ではコンパイルできず
・matlab7.0ではexeファイルを生成できるけど、実行できない
となっていました。

プログラムは用意したmp3(またはwav)ファイルを読み込んで、その波形から特徴

のあるところを探すものらしいですが、400行ちょっとあって、どこが原因なのか

探すのが一苦労です。
今後もmファイルのコンパイルにはお世話になる可能性が高いので、何が障害にな

るか傾向なんかがあるとやりやすいのですが…

たったの400行で何言ってんだ……

おまえソースコードもってんだろ？エラーメッセージぐらいなんでかかんの？

なんで教授が２ｃｈできくわけ？
教授なら外国の掲示板で聞くと思うが

いつ教授が2chできいたんだ？
それに教授だからとなぜ2chで聞くのがおかしいんだ？
教授は崇高な存在だとでも思っているのか？

スカートの中を手鏡で覗く教授に似てるなと思っただけ

Final
Fantasy
Tactics

以前のhelpだと
inv は、LINPACK のサブルーチン ZGEDI と ZGEFA を使っています。
とありますが、
最近のは
invは、LAPACKルーチンを使って、逆行列を計算します。
行列    関数
実数    DLANGE, DGETRF, DGECON, DGETRI
複素数  ZLANGE, ZGETRF, ZGECON, ZGETRI
となっています。具体的には、何が違うのでしょうか？

対称スペクトルをifftしても虚部の時間信号で出てくるのはなぜでしょうか？

Matlabのアカデミックパックの価格って誰かご存じないでしょうか？

>>162
基本システムは，去年購入した価格で１０万円強だったが，価格改定されてるかも．
サイバネットシステムのサイトからメールで資料請求しなされ．
すぐにカタログと価格表を送ってくれるよ．

うむ。基本パックが11万円くらい。各ツールボックスが4万5000円くらい。

MATLABにLISPのmap関数の様なものはありますでしょうか?

>>165
なんでそんなものが必要だと思うのかを書かないの?

>>166

for文を使わずに，
'A.*B'や'A./B'の様に書きたいからです．

「何を」for文を使わずに'A.*B'や'A./B'の様にの様に書きたいのかを書いたら?
それは本当にlispのmap関数のようなものがないとできないのか?

>>168
おまえはまず興味があるのなら、Yes か No で答えてあげたらどうだ？
疑問文に疑問文で答えろと学校で教えているのか？ってやつだ

さいきんなぜだか知らんがこういうやつ多いよ

質問するなら何がやりたくて何が分からないのか全部最初に書く。
小出しにして質問者に質問させるような質問に親切に答えてくれる人は少ない。
>>169はそういう親切な人みたいだから、君が答えてあげたら?

>>171
あんたが質問の意味も理解できない馬鹿なだけだろ

mapでやるようなことをmatlabでは違う方法でやるのが普通。
>>172はちゃかしてるだけで結局答えることも
どう聞くべきかを教えてやることもしないんだろ?

普通じゃない方法はやるなと主張するわけだな？

>>174
そんなこと一言も言ってないのになんで勝手に決めるんだ?
でお前は質問に答えられもしないし
どう聞けばより答が得られやすいかを教えてあげるわけでもないんだろ?

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

pythonならmapなかった？

ちなみに
python matlab compatible
で検索すると幸せになれるかも

>>176
検索したら巡回セールスマン問題の解法アルゴリズムは
たくさん出てくると思うんだけど。

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

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

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

>>176
traveling salesman annealingで検索したらそのまんまなのが見つかるな。
さすがに宿題でこのコード出すわけにはいかないだろうけど(笑)。

デフォルトの名無しさん、丁寧なご回答ありがとうございます。
参考にさせていただきます。

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

高さがあるのに平面とは…

mov2aviでPowerPointで再生できる動画を作成したいのですが、何を指定すればよいのでしょうか？

>>184
help mov2avi
何が出来て何が出来なかったのか言ってくれないと判らない

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

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を使えば良いのか、検討が付かない。
何か手掛かりだけでも、誰か教えてください。

>>186
で、そのフィルターが1個のときと1000個のときで
具体的にどう表示して欲しいんだ?

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

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

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

すごい初歩的なことなのですが少数点以下をもつ数の累乗根を出すような関数やコマンドは何でしょうか？

ヘルプは当然見たんだろうけど、どうやって調べた?

>>192
helpでコマンド一覧みたいなのを見たのですが整数の累乗根を出すものしか見つけられなかったので質問させてもらいました

>>193
うちでは整数のみのそのコマンドが見つからないです。

等間隔で行を抽出するにはどうすればいいですか？
例えば、1行目、11行目、21行目…という風にです。

1:10:end
これじゃダメ？

matlabで箱ヒゲ図を描きたいです。
ヒゲは、errorbarを使えば何とかなるのですが、
箱を描く方法が思いつきません。
どなたか、ご教授してください。

http://www.sci.kagoshima-u.ac.jp/~itls/Japanese/chapter3/topic/baseball/box.html
こういう図が描きたいのです。

株価のローソク足に似てるなーと思って「MATLAB ローソク足」でググったら一番上にそれらしきものがでてきたよ
Financial Time Series Toolboxがいるらしい

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

>>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？

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

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

>>196>>199
できました、ありがとうございます。

MATLABには行列内部の最小の値を取る関数はあるけど
最小の値のある場所のインデックスを返す関数ってないよね

同じ関数

help その関数
で、わかりますよ＾＾

ygridの位置を指定するにはどうすればいいのでしょうか？

もっと詳しく書いて欲しいところだが、
set(gca, 'YTick', [ ... ]);
とかでどう？

matlabのグラフで，プロットした波形にマウスを当てると，
その時点での縦軸と横軸の値が出ますが，
それを，コマンドプロンプトでやりたいんですが，
関数かなにかありますか？

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

>>207
つ[gnuplot]

>>208
ご教授ありがとうございます．

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

ということは，
そういったことができるmatlabコマンドはないんですね・・・

画像処理でアフィン変換できる関数が見当たらないのですが気のせい？
結構基礎的な関数の気がするんだけど

そのくらい自分でやれ、ってことじゃないの？

まとめサイトはありませんか？

v=1*50000 vector & A=50000*50000 matrix

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

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

>>202
ヘルプ見なさい

>203,204,215
helpみました
吊ってきます…

Axes上に複数のLineオブジェクトを作成しました。
このLineをマウス操作で移動したり、リサイズしたりしたいと思ってます…。

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

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

3次元プロットで、（例えばsurfなどで）値がaxisで指定した範囲外にあってもプロットされてしまうのですが、
axisで指定した範囲内のみプロットするようにするにはどうしたらよいのでしょうか？

関数を定義する場合必ず別ファイルにしないといけないのですか？
同一ファイル内で定義することはできないものでしょうか？

>>218
aってのを3次元の行列だとして
surf(a(xmin:xmax,ymin:ymax,zmin:zmax))
じゃだめ？
意味取り違えてたらごめんね。

>>219
http://pc11.2ch.net/test/read.cgi/tech/1163125702/82-
同じスレなので出来れば読んでね。

matlabってリーもとデスクトップじゃつかえないのだね

bata版のNAG ToolBox for MATLAB
無料評価版

http://www.nag-j.co.jp/ToolBox.htm

pythonで大部分の数値計算できるけどmatlabってやっぱりまだ需要あるの？

>>223
C で python のできることは全部できるけど python って需要あんの？

俺紙とエンピツで全部やってるけどコンピュータって需要(ry

さすがにMATLABのスレでそりゃないだろ。
手作業でやったら一生かかっても終わらん計算してるだろうに。

質問なのですが、下記のプログラムを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

それは質問じゃなく変換依頼

タブ区切りのデータなのですが、一行目が見出しになっています。
一行目をのぞいて簡単に読み込む方法はないでしょうか？

fgetl

たぶん激しく初心者な質問でスマソなんですが、サイズがL,M,Nな3次元変数を3次元プロットする方法知りませんか？

強度イメージから一定範囲の強度値だけを抽出する
方法を教えてください

田中真紀子のﾎﾟﾛﾘ画像ありますか?

オナニーをすると身長が伸びないというのは本当ですか？

>>231
自己解決しますた

質問させてください。
specgramの逆の操作がしたいです。
X=specgram(x,128,8000);
で、
スペクトログラムXから音xを合成したいのですが、どうしたらいいですか？？

synthesizeとかいう命令があったような気がする。
matlabじゃなくてoctaveの命令だったかも。

Core2duoに変えたんだけど、６４BITに変えたら
計算すごく速くなる？
ベンチとかないかなぁ

入力したサウンドデータを，何の音なのか
(例えば，パトカーのサイレン，人の叫び声等)
分類するソフトを個人的に作成しています．

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

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

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

スレ違い

自分の目でいういろ比較して確かめなさい

はじめまして、matlab初心者です。
早速ですが質問させて下さい。

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

宜しくお願いします。

matlabのメインウインドに打ち込むだけでsimulinkの開始ボタンを押したのと同じ処理になるコマンドってありますか？

>>243
SIM

例えば４ｘ４のマトリックスを２ｘ２にして、各要素は元々のヨウ素を加算してくれる関数はないですか？
ループとsum関数を使えばすぐに出来るんですがね

>>211
Image Processing Tool Boxだったかにあるよ。座標変換してからinterp2に流し込めばすぐにできるけどね。

>>245
じゃぁ、ループとsum関数を使えば？

>>247
matlab ではできるだけループを使わないのが重要
速度がダンチ

書き方が曖昧でわからんがエスパーの俺ならこうかな。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);

:と,一カ所ミスった。

??? エラー: ==> anfisedit
SCRIPT anfisを関数として実行してみてください.

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

>>251
SCRIPT anfisを関数として実行してみてください.

MATLABで式をニュートンラプソン法でとくにはどうしたらいいでしょうか？
分からなさすぎてノイローゼになりそうです。

まぁスレ違いだわな

GUI上でxの変域をかえたいんですけどどうすれば良いですか？

GUI上で数値を入力させて、
set(h, 'xlim', [x1 x2]);
じゃダメなの？

FORTRANで書いたものを、
MATLABに変換したいのですが、手書き以外で
何かいい方法あったら教えてください。

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

>>257
グーグルの使い方くらい知らないと
今後いろいろ苦労すると思うから覚えた方がいいよ。

ファイルの読み込みで指定した列だけ読む方法ってありましたっけ。

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

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

>>260
dlmreadのマニュアル読むと良いよ。

対数データのフィッティングの方法を教えてください。
通常のデータであれば、最小二乗法ですが、対数データの場合はどのようにやればよいのでしょうか？

対数の底で冪をとって最小二乗法

「MATLABの極」というブログを書いています。良ければ、ご参考に
http://beta-infinite.blogspot.com/

著作権違反ですよｗ

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

今、手一杯だからなー。
月８０万だったら考えるけどね。
それかNNツールキット買えばいいんじゃないの？

ていうか自分で作れって感じ
層やニューロン数を完全に一般化した物でなければ速攻で作れるでしょ
もちろんfor文抜きでな

料金は
３時間未満：一律12000円
３時間以上：時給4000円
７時間以上：一律28000円
このあたりでよろ

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

安すぎ

委託研究でMATLABでソフト作ったときは、３０万もらったな

作成物の複雑さとか教えてもらえずに
7時間以上28000円は安いな。

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

こんなところで依頼しようとしている時点でダーティなイメージがついてて嫌だ
全然予算がないようだし学生の卒論か？

>>275
まさにそれ
誰か頼むよ〜

http://www.mes.p.dendai.ac.jp/index.php?plugin=attach&refer=%A5%B3%A5%F3%A5%D4%A5%E5%A1%BC%A5%BF%A5%D7%A5%ED%A5%B0%A5%E9%A5%DF%A5%F3%A5%B0%AD%B6&openfile=CP7.pdf
↑この最後のページにある問題を教えてください。
参考の為に問題文だけ書くと・・・

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

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

別に matlab じゃないとできないわけじゃない。プログラミング言語に依存してたまるか

請負内容を見て判断したいです。詳細内容などをメールをください

>>267
>>270
請負内容を見て判断したいです。詳細内容などをメールをください。
秘密厳守を守ります。

手元にはバックプロパゲーション法の自作ソースくらいあるんだが
バカ学生に使いこなせるとは思えん

matlab にファイル名を変えるコマンドってありましたっけ？
同時にファイル内の関数名もあわせて変えてくれるとベスト

ttp://www.cybernet.co.jp/matlab/support/techkwdb/search_view.php?id=1303

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

matlabで書いたグラフで、線と線に囲まれた部分の
面積を求めることは可能でしょうか。
ファイルをアップしましたので、よろしかったら見てみてください。
http://www.csync.net/service/file/index.html
(ファイル名はmatlabです)

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

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

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

>>283
「MATLABの極」にソースをアップロードしました。いかがでしょうか
http://beta-infinite.blogspot.com/

教えてください。
MAT LABで書いたFigureのX軸Y軸をX→Y、Y→Xと逆にしたいのですが
どうすればいいんでしょうか？お願いします。

>>285
ありがとうございます。
今までhold onとか使っていたんですが、こんなシンプルにいけたんですね。

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

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

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

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

>>287の修正

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

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

>>286
view(90,90)という関数を使えば、いいと思います。
http://beta-infinite.blogspot.com/2007/11/matlabfigurexyxyyx.html
に参考

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

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

http://beta-infinite.blogspot.com/search/label/%E3%82%A2%E3%83%8B%E3%83%A1%E3%83%BC%E3%82%B7%E3%83%A7%E3%83%B3
に参考

beta-infinite.blogspot.comってパクリしかないじゃんかｗ

>>293
MATLABの修業中です。覚え書きのブログだと思っていただければ、

>>292
やっぱgetframe,movie,movieしかないですかね
ありがとうございました

著作権違反だので、通報しました

FFTを使って、スペクトルを計算するのって、どうすればいいの？

だので

>>297
help fftとコマンドを打って出てきた説明を理解できなければ
回線を切って死ぬか教授に頭を下げて教えを請いなさい

>>297
beta-infinite.blogspot.comだけは見ないように。

fftはtwo-sideの出力を出すので、help fftだけでは話は完結しない

教えてください。
プログラムで図を描かせたときにX,Y,Zラベルの刻みのラベルを
変更したいのですが、関数などプログラムで書かせることができるでしょうか？
少し伝わりにくい文章ですが、よろしくお願いします。

なんか楽に関数名リネームする方法ねーのかよ

>>303
関数名リネームって>>282が言ってるようなことか?
もしそうならファイル名と違ったらファイル名の方が優先なだけだから
コード中の関数名は適当でOK。

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

>>297
ここ見るといいよ
http://www.mathworks.com/support/tech-notes/1700/1702.shtml

>>302
できる。

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

シミュレート板で回答を得られなかったので。

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

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

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

正の実数からなる対角行列の、各行と列について、
要素の大きさが対角成分の要素の大きさの0.5から2倍の範囲なら１で置換、
それ以外なら０で置換する、
ってforループを使わずにできるでしょうか？

質問ばっかで誰も答えてくれない罠ｗｗｗｗｗｗｗｗｗ

>>311
対角行列なら対角成分以外すでに全部０やん

これぞまさに質問スレ

vistaは2007aで正式サポートですが、2006aで動かした場合、どこかに既知の問題のようなリストはないでしょうか？

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

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」

と、こんな感じなのですが,どなたかサクっと作ってくださいませんでしょうか・・・

秀丸用の強調定義ファイルはありませんか？

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 に行っている。
そして最後の行は、他の行と読みこむ数が違う。

fgetlで１行まるまる読んで、strtokで切り出していくしか無いんじゃないの？Ｃでも普通そうだろ？

>>319
それを一発で君の思い通りに読んでくれる
ファイル読み込み関連の関数が貧弱じゃないやつって何？
煽ってるんじゃなくて単に興味があるから知りたい。
ExcelとかRとかできるんかな?

excelじゃ無理だね。

figureウインドウを表示させずにムービーを作る方法を教えてください

>>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だとフォーマットの指定に不足があるのでしょうか？？

>>324
textread は一行ずつ読み込むんだろ？
最初の行は先頭以降何もないから０をが入っていくと。
[data1,data2,data3,data4,data5,data6,] = textread('ファイル名.txt', '%f %f %f %f %f %f')
として、奇数行と、偶数行に後でわけるとかになるんじゃね？

どうやらそのようだ、と思っていたんですが
以前同様のものを読み込んだ際に行ごとでなかったんですよね

たとえばこんなデータ

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（スペース）を読み込んでるので
行をまたいでも効果があることになると思うんです。

1 1 0.00000000E+00 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.58726901E-12 0.55467926E-06

あああすいません、変な書き込みしてしまいました
1 1 0.00000000E+00 0.00000000E+00 0.00000000E+00
0.00000000E+00 0.58726901E-12 0.55467926E-06

どうやらうまく書き込めないのですが
１番目、２番目の１の下にはスペースが存在していて
３番目の0.00000000E+00 と下の行の0.00000000E+00 は
列がそろっています。

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

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

>>329
textreadにこだわらなくてもいいんじゃないの？

matlabでfor文をまわす時
ループが終わるごとに変数をclearした方がよかったりします？

不振な挙動。。

そんなことないと思うが。
それにしても、Matlabバージョンとか書かないと。

たぶんそれはいくつかの変数を間違って使ってるかロジック間違えてるだけだ。
ちゃんと全部読み直せ。

波形処理をやっていてカーブフィッティングをしています。
fit関数を使ってフィッティングをしているのですが、
フィッティング後の曲線のパラメータが読み取れません。
fit1 =

一般ﾓﾃﾞﾙ Exp1:
fit1(x) = a*exp(b*x)
係数(95% 信頼区間付き):
a = -10 (-10.17, -9.831)
b = -0.01371 (-0.01403, -0.01338)
こんな感じでフィッティングしたあと、bの値を使いたいのですが、
どうすれば読み取れますか？
初心者の質問ですみません。

2値化した画像をカラー画像にするにはどうしたらいいですか

それはもう情報がかなりないからとても難しいあるね

>>337
カラーは白黒でよいのですがそれでも無理ですか?

>>338
それは、フォーマット変換したいと言うこと?
それとも、2値のメッシュをグレイにしたいと言うこと?

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

行列の積で初歩の初歩的な質問だと思うのですが誰か教えてください！
10000x10の行列Aがあるとします。
A*A'だとout of memoryと出力されるのですが
A'*Aだとちゃんと計算結果が出力されます。
考えてもイマイチ理由が浮かんで来なかったのでどなたか教えてください〜！
単純な質問で失礼しました。

前者は 10000x10000
後者は 10x10
の出力になるだろう？明らかすぐる

これはひどい

（・A・）＊（・A・）'

ホントだｗｗｗ頭悪くてすまそｗｗ
ありです。

plotで描ける線の種類を増やせないでしょうか？
6種類のデータを1枚の図に、白黒で書く、ということで困っています。

>>346
解決しないかも知れないけどマーカと併用するのはダメなんですか
線種で区別できる程に線が絡み合ってないのならそれでいけるのではないでしょうか

>>339
フォーマット変換のほうだと思います。

>>348
MATLABでなくていいなら、画像ビューワでも画像エディタでもいくらでもやりようがありそうだが。
何を、何故、MATLABでやりたいのか全然判らんのだよ。

>>349
MATLABで画像処理の勉強をしろといわれ、つい2値化の画像を白黒だけの
カラー画像に変えてやろうと思った。

まずおめーは、白黒だけのカラー画像とかわけわかんねーことをいうな
カラー画像っつったらＲＧＢとか３チャネルあるやつを普通は指すんだよ。

あほす
R == G == B のことだろ

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

ステップ応答の試験をしてMATLABで伝達関数を出したいんですけど、
最小二乗でフィットさせるときの方法とかいうのが分からないので教えてください。

テキストファイルの行数を事前に知る関数はないでしょうか？

>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; % 緑に代入

…でいいはず（うろ覚え

>>353
どういう関数系にフィットさせるの？

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

安価ミス
×　>>357
○　>>356

あ、あと私のMATLABにはSystem Identification Toolboxが入っていないのでこれを使わない方法を教えて欲しいです。

そこまでわかっているのなら、最小自乗でフィットするだけじゃん

ttp://www.wakayama-u.ac.jp/~miw/ApIp/01/2.html

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

fminsearch関数使えば？Simplex方による繰り返し計算だけど。

数学的には間違いだけれども、
対数変換して直線フィットすると、ある値が出る。
それを初期値に>>362にしたがって、解を探す。

>>360
ごめんなさい
その最小二乗でだけっていうところが分からないでの聞きにきたんですよね・・・

polyfit使えば一発。ただ原点に返って正規方程式を逆行列で解くのもありだな。

>>364
MATLAB以前の問題なのか…

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

>>365
ployfitって単なる多項式しか出来ないんですよね

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

だから単純なたとえばy = a*x^2 + …とか>>360に書いてあることはわかるんだが、
伝達関数の場合のプロセスが分からないっていうことなのだが教えてもらえないのか？
なにを調べたらいいか分からないからここに来たんだが

>>368
あなたはどういうデータを持っていますか？
理論的には、どういう操作をしますか
どういう式にあわせて、係数を推定したいですか？

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

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

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

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

370==372
なのか？

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

というよりイライラしていたら言葉が悪くなるんだよ

問題をあなたが理解していないのでは、と思わざるを得ない。
なぜなら、あなたが何をしたいのかが、伝わらないからである。

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

何をしたいか：伝達関数を出したい
教えて欲しいこと：伝達関数の求め方

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

>>376
それで伝わっていると思うのなら、
あんたの理解が甘い。

先生に教えてもらえ。

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

polarでMarkerSizeを変更する方法を教えてください

ふつうのグラフで軸をスケーリングするには
AXIS([XMIN XMAX YMIN YMAX])
とすればできますが、
subplotを使ったグラフの軸をスケーリングするには
どうしたらいのでしょうか？

>>380
同じ

>>381
同じでできるんですか？
AXIS([XMIN XMAX YMIN YMAX])としてもできなかったので
聞きにきたんですが。もう一度試してみます。

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

>>382
subplot(211)
plot(x,y)
axis([xmin xmax ymin ymax])
%
subplor(212)
plot(x,y)
axis([xmin xmax ymin ymax])

>>382
subplot(2,1,1);
plot(x, y);
set(gca, 'xlim', [xmin xmax]);
set(gca, 'ylim', [ymin ymax]);

でもいけるのでは？

>>381>>383-384
できました。ありがとうございました！
自分が馬鹿でした。
axis設定してからplotしてました。
plotしてからaxisなんですね。

>>384
set(gca, 'xlim', [xmin xmax], 'ylim', [ymin ymax]);
で桶では？

逆ラドン変換iradon()を使うとデフォルトでRampフィルタがかかってしまうのですが
フィルタ無しの方法はありませんか？
海外のMATLABのヘルプでは値でNo filteringにできるそうなんですがバージョンが
古いのかバージョン7.2では出来ません。

逆投影なんかツールボックスを使わなくても数行でかけるのに(´･ω･`)

なんとかヒントでもよければ教えてもらえないでしょうか？
専門外のことなのでどうにもこうにもツールに頼ってしまってorz

edit その関数
で、中身見てフィルタ部分カットしたのを書くのはどうですか？

まあ、簡単に考えると配列の方をirotateで回転させ、そこに加算していけば２行で書けるんだけどね。
それで角度でループさせるの。

iradonの中を触るんだったら、Hというのがフィルター関数なので、この部分を取り去るか、1にしてしまうか。

勉強のために１から作りたいんだったら、Ｃ言語で出ている教科書を買って、ネチネチと移植してもいい。
おれは院の先輩の薦めで「SPECT画像再構成の基礎」（医療科学社）という本の中に出ているフーリエ変換法のプログラムを移植したよ。

グラフの大きさを変えるにはどうすればいいのでしょうか？
いつもマウスでウィンドウのサイズを変えて
縦横比を見やすく調整しているんですか
コマンドでできないでしょうか？

position

set(2,'position',[30 246 1120 420]
って感じですね？
できました！
ありがとうございます

figure(1);
print -deps 〜.eps

のようにグラフをepsファイルに出力してるんですが
その出力したepsファイルのサイズ（縦横比）を変えることはできないでしょうか？

set 'position'でできるかと思ったらできませんでした。

自由落下した物体が床で跳ね返る時、速度の符号が
変わりますよね。その時ゼロクロッシングとかいう
エラーでちゃいます。ゼロクロッシングの検知を
無効にするとうまくシミュレーションできません。
何でですか？

483 名前：名無しさん＠５周年[] 投稿日：2008/01/29(火) 16:04:56 ID:
自由落下した物体が床で跳ね返る時、速度の符号が
変わりますよね。その時ゼロクロッシングとかいう
エラーでちゃいます。ゼロクロッシングの検知を
無効にするとうまくシミュレーションできません。
何でですか？

set(gcf,'PaperPositionMode','auto')かな？

>>398
printに-zbufferとか-openglオプションってのを
つけてみてもできなくて困っていたのですが
ずばりそれでできました！ありがとうございます

MEXをVisualC++2008ExpressEditionで
開発するときの設定の仕方ってどこか
にありませんでしょうか？m(_ _)m

x=[0:1000]に対応するf(x)が手元にあって
各xについてx-x/10<=x <= x+x/10の範囲でｆの移動平均をとって平滑化をしたいんですが
ループを使わない方法はないでしょうか？

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)の部分を赤に色付けしたいです

Control System Toolboxを使われている方に質問です。

ゼロ次ホールドの連続時間伝達関数である、
H = (1-exp(-sT))/s ;　　T:サンプリング周期[sec]
という表現をmatlab上で再現するにはどのようにしたらよいのでしょうか？
s = tf('s');
を使って直接表現しようと試みたのですが、上手くいきませんでした。

連続時間システムの制御器の設計及びコンピュータへの実装の際に
使用したいと考えています。
よろしくお願いします。

構造体にセットしたデータをfwriteとかでバイナリファイルに書き込むのはどうしたら良いのでしょうか？

素人質問ですが教えていただければ幸いです。

10個の自動計測されたデータが列方向に並んでいるcsvファイルを
csvreadで一気に読み込んだ後に(例えばAとする)、
10個の列要素に小分け(A1〜A10)しようと考えています。

各要素の1行目は説明なので2行目から最終行までを選択して取り込むには
どうすればよいのでしょうか?

思いつきでA1=A(:,2:inf)とかしようとしたらエラーがでてしまいました。

多分while文とかでNaNがでるまで読み込み続けるとかすると思ったのですが
力尽きてしまいました。

やさしい人。お願い。

infじゃなくってendじゃないの？

>>405
よく読まずに書いてしまうが、infではなくendを使ってみたらどう

>>２　請けた仕事について、２ｃｈで無料で意見伺う無神経さｗ
というか綱渡りで仕事してるからけつに火がついてしまうのに懲りてないな。
ベンチャー高専隊

A(:,1)=[]
で１列消せるはず

407さん、ありがとうございます！おっしゃる通りで解決しました。

さらに色々触ってみて疑問が沸いたので厚かましく質問させてください。
（長文ですみません、誰も呼んでくれないだろうなぁ・・・）

例えば、
-----------------
time　data1 data2
1 10 100
2 11 99
3 10 102
4 9 98
| | |
-----------------

というデータをcsvファイルから読み込んで、
data1とdata2を1/1000倍してから
csvファイルとかExcelファイルに上と同じ書式でE書き込もうとして
下記のコードを書いてみました。

結果、xlswriteではうまく動くのですが、csvwriteではエラーが出ました。
セル配列を使えなかったようです。

私が自宅で使っているMATLABはバージョンが古く(ver.6.5)、xlswriteが使えません。
できれば古いバージョンでも動くファイル処理関数で対応できないものかと思っています。

(続き）書いてみて失敗したコードです。ご指導お願いします！！
------------
clear all

Alldata = csvread('csvsample.csv', 1, 0);

Time = Alldata(2:end, 1);
Data1 = Alldata(2:end, 2);
Data2 = Alldata(2:end, 3);

Data1 = Data1/1000;
Data2 = Data2/1000;

Dataset = [Time Data1 Data2];
Dataset_cell = num2cell(Dataset)

Label = {'Time', 'Data1', 'Data2'};

Outdata = {Label ; Dataset_cell};

%xlswrite('newxls', Outdata)
csvwrite('newcsv.csv',Outdata)

みなさんMATLABのfigureをWordとかPowerPointにはるのにどうしてますか？

今まではBMPで保存してXnViewの自動トリミングで余白を削ってはってました。
ただ、BMPだとWordでグラフを拡大縮小したりするとフォントがかすれたりして見た目がかなり汚くなります
ベクターフォーマットのemfファイルで保存するのが一番きれいなのですが、グラフの余白を簡単に削る方法が見つからないです
XnViewでもemfファイルを読み込めるんですが、XnViewで自動トリミングしてemfファイルに保存すると、
ファイルサイズが10倍以上になって、しかもビットマップデータみたいに拡大縮小すると汚くなります

グラフの余白を簡単に削れてきれいに貼り付けられるうまい方法とかないですか？

自己解決しました。スレ汚しすみませんでした。
泥臭く下記のように書いてみました。
データ整理するのにxlsreadとかxlswriteって便利だなぁと思いました。
---------
clear all
Alldata = csvread('csvsample.csv', 1, 0);
Time = Alldata(2:end, 1);
Data1 = Alldata(2:end, 2);
Data2 = Alldata(2:end, 3);
Data1 = Data1/1000;
Data2 = Data2/1000;

Fid=fopen('newcsv3.csv', 'wt');
Label = {'Time', 'Data1', 'Data2'};
for i=1:length(Label)
fprintf(Fid, '%s,', Label{i});
end
fprintf(Fid,'\n');

Dataset = [Time, Data1, Data2];
[Rownum, Colnum] = size(Dataset);
for i=1:Rownum
for j=1:Colnum
fprintf(Fid, '%d,', Dataset(i,j));
end
fprintf(Fid, '\n');
end
fclose(Fid);

すいません、皆さんみたいにハッカーになるためにはどんなことが必要ですか？

>>412
グラフの大きさを指定して、拡大縮小せずに使えるようにする

MATLABでドロネー分割をつかって、３角形のパッチを利用した形状モデルを作っているのですが、パッチにテクスチャを貼り付けるにはどうすればいいでしょうか？

3次元の立体形状から有限要素法用のメッシュを作り出す方法を教えてください

テクスチャマッピングの仕方を教えて欲しいです。
わかりやすいサイトとかがあればお願いします。
よろしくお願いします。

要求が高いなぁｗ

まず、www.google.co.jp。日本語以外のサイトも調べよう。

２次元のグレースケールの画像を、Z軸がピクセル値で３次元表示させたいんですけど、どうすればよろしいですかね？

>>420
surf

すいません，
a = 'Hello World'
を
b[1] = 'Hello'
b[2] = 'World'
に分割するコマンドはありませんでしょうか．

doc regexp

>>422
[ str,rem ] = strtok( a,' ' );

これで、最初のスペースまでの単語がstrに入り、残りがremに入る。
remの中をさらに切り分けたければこれを繰り返す。
str系の関数はC言語とほぼ同じだな。

WAVファイルをパワースペクトル表示する方法を教えてください
お願いします

ぐぐれ

ググっても見つからなかったんで質問させてください。

h1 = plot(x1)
hold on
h2 = plot(x2)
h3 = plot(x3) ...

といくつかのプロットを重ねた後、一つだけ(例えばx2のプロット)だけ消したい、または変更したいという場合、どうしたらいいですか？
全部プロットしなおせばいいんですが、ちょっと重い上に何度もやるプロットなんで、全部やり直すというのは避けたいんです。

書き忘れてました。プロットツールは使わず、関数で出来る方法で探しています。

http://www.google.co.jp/search?hl=ja&rls=GGGL%2CGGGL%3A2006-17%2CGGGL%3Aja&q=matlab+line&btnG=%E6%A4%9C%E7%B4%A2&lr=

>>429
ありがとうございます！

グラフの凡例中で改行するにはどうすればよいでしょうか？

plot(X1)
hold on
plot(X2)
plot(X3)
legend({'X1';'test'},'X2','X3')

とタイトルや座標軸中の改行と同じように入力しても無理でした。

100枚の画像あってサイズが縦240横352なんですが
横の左50画素、右46画素をカットして横幅を256画素に、
30×30の位置の色と同じものを上に縦16×横240くっつけて
256×256の画像にしたいのですが、どうすればいいでしょうか？

つ[ImageMagick]
なんでもかんでもMATLABでやりたいならしらね。
つーか、MATLABエキスパートに任せた。

>>423
>>424
ありがとうございました！
6.5なので423での'split'は無理みたいでした．

質問です！
ある画像の一部分（オブジェクト）の画素値の平均をとりたいのですが、
ラベル付けして、オブジェクトの全ピクセルの座標は取ったのですが、
for i=r(min):r(max)
for j=c(min):c(max)
heikin=mean2(Img(i,j,:));
だと、オブジェクトの形が四角ではないのでオブジェクトの外側の画素値まで
計算してしまいます。

何かいい方法がありましたら教えてください。

>>435
オブジェクトがなんなのかわかりませんが、その内側と外側で違う特徴的なものをfindで引っかけてやればいいのでは？
たとえばオブジェクトの外側が-100より必ず小さいのであれば

m = find( Img( : ) < -100 );
heikin = mean( Img( m ));

の２行でできるはず。
Imgの範囲を絞りたいのだったら、

i = r(min):r(max);
j = c(min):c(max);
m = find( Img( i : j ) < -100 );
heikin = mean( Img( m ));

forループは極力使うのは止めましょう。

ご回答ありがとうございます。
バッチリ上手く行きました。

差しでがましいですが
もう一つ質問が…

二値化の画像にラベル付けをしてあるのですが、
ある座標を与えて、その座標を含むオブジェクトのラベル番号、
もしくはそのオブジェクトを構成している全座標を入手したいです。
[r,c] = find(bwselect(Img2,i,j));
これが僕の絞り出した答えでしたが無理でした。。

いいアイディアがありましたら教えてください。

ご回答ありがとうございます。
バッチリ上手く行きました。

差しでがましいですが
もう一つ質問が…

二値化の画像にラベル付けをしてあるのですが、
ある座標を与えて、その座標を含むオブジェクトのラベル番号、
もしくはそのオブジェクトを構成している全座標を入手したいです。
[r,c] = find(bwselect(Img2,i,j));
これが僕の絞り出した答えでしたが無理でした。。

いいアイディアがありましたら教えてください。

>>432
画像の読み込みをまずやってみよう

フルカラーを使ってのムービーを作成する際のコーデックは何がおすすめですか？

無圧縮AVI
後で別のソフトで変換

質問させてください！
統計的に（最小２乗法などを用いて）２次近似モデルを式として作成する関数を探しています。

例えば実験にて２つの入力を調整し、１つの出力が得られるものがあるとします。
そのときの計測条件を[x,y]（２入力、ともにベクトル）、計測値を[z]（１出力、ベクトル）としたときの
計測値Zに対するX,Yの２次近似モデルを作成したいのです。
アウトプットイメージとしては
Z=係数*X^2+係数*Y^2+係数*X*Y+係数*X+係数*Y+係数
のような感じです。

プロの皆さん、どうか教えてください〜

>>442

z=[ z1,z2,......,zn ] '、x=[ x1,x2,......,xn ] ' , y = [ y1,y2,......,yn ] ' とし

z = A*x^2 + B*y^2 + C*x*y + D * x + E*y + F

とすれば

a = [ x.^2 y.^2 x.*y x y ones( length( x ),1) ];
x = [ A B C D E F] ' ;

とおいて

x = a\z;

で解けるはず。未知数が６つなので、データ数が少なくとも６つ以上ないと無理だけどね。
重み係数を使う必要がある場合は、ノルムを最小にするようにしないといけないけどね。

>>437
find関数の中は条件文じゃないと機能しないのでは？

>>443
ありがとうございますー、やってみます！

>重み係数を使う必要がある場合は、ノルムを最小にするようにしないといけないけどね。

この意味が理解できないｗ

>>446
ノルムというのは距離のこと

重み付き最初二乗法はどうやって記述できるの？

なぜmatlabスレで数式をきくんだ？
matlab への翻訳がわからないのなら数式は記述しておくとか、
数式が記述してあるページとかを示しておくべきなんじゃないのか？
数式がわからないというのはスレ違いだろ？

>>449
ウルサイ！

シミュ板の scilab のスレッドが過疎で
レスがいただけないようなので、こちらで質問させてください。

MATLAB で

( x^2 ) + ( x * y ) + ( y^2 ) = 1;

主軸が斜めになっている楕円は
どのように plot2d すればいいのでしょうか？

>>447
いや、そもそも変なこと書いてるのではないか？

あるwaveファイルにQ値、ゲインをある特定の周波数に持たせたい
（EQをかけたい）のですが、どうすればよいのでしょうか？
フィルタを使う（fdatool）ではダメですよね？
宜しくお願いします

>>453
ここを参考にすべし
http://www.musicdsp.org/files/Audio-EQ-Cookbook.txt

http://dl.cybernet.co.jp/matlab/support/manual/r2006/toolbox/matlab/visualize/?/matlab/support/manual/r2006/toolbox/matlab/visualize/chvolv19.shtml
このような3次元の物体の可視化を行っているのですが、物体の表面で金属のような質感を出すにはどのように
設定したらよいのでしょうか？

質問させていただきます。
「MATLABによる制御系設計」という本で勉強しているんですが、
MATLABの制御系コマンドを使えるフリーソフトを探しています。

検索すると無料でScilab,Octave,Freemat,有料でMATXというのがあるようですが、
制御系コマンドを使う時、互換性が最も高いなどお勧めの
ものがあれば教えてください。

有料だけどMATLABっていうのがあるらしいよ

>>457
456
もちろんMATLABは知ってますが・・・・
研究でMATLABは使ってきたのですが、これからは
家でも勉強したいのであのような質問をしたのですが。
拙い質問と感じられたのでしたら申し訳。

Scilabでいいんじゃね

Octaveでいいんじゃね

Freematでいいんじゃね

とりあえず見つけたやつ全部自分で試して見ればいいんじゃね

学生ならMATLAB買え

>462
ごもっともだと思います。
知り合いにもそういう意見をもらったんで、
そうしてみます。とりあえず、Scilabインストールしてみました。
金銭的に厳しいのでフリーのソフトから自分なりに試して
みようと思います。

>463
それもそうですね。
MATLABをインストールできるだけお金が貯まったそうしようと思います。
４月から社会人なので。

意見ありがとうございます。

前の職場でmatlab使ってたけど今のところでmatlab買えない
pyrhon.numpy
と
maxima
どっちがお勧め？

俺的には;
python + scipy + matplotlib

グリッドの幅を変更したいのですがどうやって指定すればいいのでしょう？
grid minor　より細かくしたいでっす。

Unix用Matlabで圧縮ムービーを作る方法はありませんか？
movie2aviだと非圧縮しかサポートしていないようなのですが。

行列の平方根を求める方法を探したら、関数「sqrtm」が使えることが分かった。
例えば、
Y = X * X;
でYが与えられたとき、解Xを求めるには、
X = sqrtm(Y);
とすればいい。


ところで、次の場合、解Xを求めるにはどうすればいいんだろう？
分かる人いる？

Y = X * X'
※ 行列はすべて複素行列で、「'」は複素共役転置を表す

教授から言われて、MATLABで作ったプログラムをコンパイルしています。
コンパイルは普通に通るのですが、エラーが発生する箇所が時々出ます。
エラーが出た箇所や解決法をリスト化して書いてあるサイトなどは無いのでしょうか？
今回は解決したのですが、毎度これを解決するのは難儀なので。

今回はsystem関数で
-----------------------------------------
"プログラムのフルパス" -コマンドラインオプション
-----------------------------------------
としてあったのを
-----------------------------------------
cd "プログラムのあるフォルダ"
プログラム名 -コマンドラインオプション
-----------------------------------------
としたら動きました。

>>470
どうでもいいが、ディレクトリをほいほい移動しちまっていいのか?
その後ファイルアクセスしてたら厄介だぞ。

UNIXのoctaveスレの人たちをここに誘導した方がいいのだろうか

お互いに行き来すればいいんじゃない？
全員誘導に従うとは思えないけど

>>471
ご指摘どうもです。
今回の場合は、外部のプログラムを1回呼び出すだけだったので問題はなかったです。

コンパイル時に動かないのは、プログラムのパスをダブルクォテーション(")で囲っているのにその後にコマンドライン引数あげている事なので…

あと、どうも他にMATLAB関連スレが見当たらないので一人事を書いておきます

・eval関数を使うとコンパイル後実行時にエラーが出る
mファイルをコンパイル時に外部mファイルが必要と判断された場合は、必要な外部mファイルは暗号化されてctfファイルにまとめられるが
eval文中で外部mファイルを
eval('foo(1);');
と呼び出すとコンパイル時に解析されないのでctfファイルに同梱されない
おそらくfevalやinline等の関数も同じ症状が出る

邪魔だったら文句言って下さい…

独り言その２
evalやfeval、inline等の関数を使うと、コンパイル後、プログラムをを実行中にエラーが出ることがある
これはMATLABに同梱されている関数でもその対象となるものがある。
その関数の中でfevalなどが使われていると、やはり同梱されない関数が出てくるため実行エラーがある
Signal Processing Toolboxのaryule 関数でそういう症状があった。

間違いがありました
aryule関数ではなく、pyulear 関数でした。
独り言その３
このようにeval関数関係でエラーが起こる場合、一番有効な方法は
ctfファイルに必要な関数を同梱させるために一度も使わない関数を書いて
function foo()
関数1();
関数2();
関数3();

とコンパイルするmファイルに追記しておくと良いと思います。

ムービーをQuickTime形式で保存する方法を教えてください

>>477
>>441

質問です

MATLABにおいてベクトルの要素が複素数の場合に最大値を求めると、
最大値はなぜそれが最大値だとみなされるのでしょうか？

長さじゃないの？

もうちょっと詳しくお願いできませんか？
3+iと4-3iの場合、どのように大小を決めるかということです

ベクトルの長さだから、sqrt( 3^2 + 1^2 ) と　srqt( 4^2 + (-3)^2 )じゃないの？

ああなるほど、そういうことですか
親切にありがとうございました

matlabを使う前に中学数学を勉強しよう

MATLABはコンパイラがあって、それを使うとC/C++で組んだコードと同じくらいの
速さになる、と聞いたのですが本当ですか？

本当だよ。１００万近くするけどね。研究室のを使っている。

>>486
ありがとうございます。友人からそのことを聞いて、いまC++を勉強してるのが無駄な感じがしてきてたのですが、
やっぱり、MATLABがそれなりにできてコンパイラもある状態なら、あたらしくC++を勉強する必要はない、という
くらいですか？実行速度で、C++で作ったコードには勝てますか？
C++で行列使うのいやなんですよね。。。それでも速度でC++の方が速いなら、やる価値はあるのかと思ってたけど、
そうじゃないならC++とはバイバイしたいス。

C/C++からmatlabを使って高速に動かす、という手もある。
こちらは100円も掛からない。

>>488
MATLABの関数を外部ライブラリみたいにつかうということですか？

uicontrol('style','text','string','TEST')

など、GUIで文字を入力した図をeps形式で保存しました。
これをillustratorで開くと、入力した文字が図形オブジェクトとして認識されてしまいます。
illustrator内で、文字オブジェクトとして認識させる方法はないのでしょうか？

movie2aviを実行すると

エラー: ==> avifile.addframe>ValidateFrame
Frame must be 440 by 439.

というエラーが出ます。ウインドウのサイズが問題なのだと思うのですが、
movie2aviコマンドに適切なウインドウサイズをプログラムで指定できませんか？

2つの変数を数十万通り、方程式に代入して解を求めたいのですが、Matlabのプログラムで可能でしょうか？
可能でしたら、具体的にどのようにすれば求めればいいでしょうか？
（Mathematica、Maple、C言語など、他のソフトも使える環境にあります。Matlabより向いているソフトがあれば教えてください。
ちなみにソフトは苦手で、どれもにわかです。）

数十万だろうが一つだろうが本質的にはかわらんよ。方程式は非線形なの線形なの？

非線形です。
方程式は、今ある状態ではニュートン・ラプソン法等で近似を行う必要があり、計算が膨大になりそうなので
ある程度ｎ、簡略化を行おうと努力してるところです。
一応、線形にできる可能性もないこともないと思います。

知りたいのは、何十万通りの2つの変数を外部データから読み込んで計算するのにMatlabは適しているのか、ということです。

質問です
rootsを使って多項式の解を求めるとき、
f(z) = z^-2 +9 のように逆数になってる場合どうすればいいのでしょうか？
また複素数面上に極座標表示するにはどれを使えばいいでしょうか？

>>494
本当に何十万通りやらなければいけないのか、考えてみることが必要かと。
プログラム自体は何パターンでも計算可能

二つの変数を取る、その問題に最適化されたプログラムを自分の手で書いたほうがいいと思うんだが

画像処理のテンプレートマッチングで、xcorrを利用している人いますか？
xcorrは、相互相関係数を計算する関数だそうですが、
http://dl.cybernet.co.jp/matlab/support/manual/r14/toolbox/signal/?/matlab/support/manual/r14/toolbox/signal/xcorr.shtml
を読んでもよくわからないです。
マッチングする領域の画素データをx、yとすると、xcorr(x,y)はひとつの値になると
思うんですが、行列で出てくるみたいです。
xcorrをどうやって使うか、教えてください。

xcorrに限らず、配列がマトリックスの場合、それぞれの縦ベクトルで別々に計算するというルールがあるの。sumとかmeanとか。
だから、これを一本の縦ベクトルにしてからぶち込んだらいい。
あとは自分で調べてね。
ここ基本だからね。

>>498
ん？なんで cross correlation がひとつの値になるんだ？
語弊を恐れず書くと cross correlation は convolution みたいなもんなんだから、|N1| + |N2| -1 の長さになるだろう？
あと、２次元データなら２次元拡張 xcorr2 がある。

matlab買えないのでRにした
octaveっていま一つ

>>498
相互相関から勉強しよう

とういうか単に相関corrを使いたいだけなんじゃないの？

PC変更したらアクチベーションは2台までらしく
アクチベーションできなくなって
PC移行時は問い合わせろと書いてあったので
米・日両方問い合わせるも反応ナッシングな俺はどうしたらいいですか。
Studentの2007年版デス。

積分が行われる様子を可視化することは可能でしょうか？
球の表面積を面積要素で表したいのですが。

>>504
土日だからじゃないの？
もう返事きたでしょ？

>>506
来ました
Machine IDの変更申請をしろとな

ウインドウを表示せずにアニメーションを作成する方法を教えてください。

MATLABで
A=atan(imag(B)./(real(B));
と位相応答を求めた後、群遅延を求めるにはどうしたらよいでしょうか？
上式を周波数で微分したらよいと思うのですが・・・
やり方がわかりません・・・

すみません。スレ違いだったらスマソ

友達がPCに詳しい自分にmatlabを使えるPCはどれがいいかと聞いてきたんですが、自分はこのソフトについて一切知らないのが現状です。
みなさんがmatlabを快適に使う環境とはどのようなものなのでしょうか？


よろしくお願いします

>>510
matlabを快適に使うとなればメモリは8GBくらい欲しいし、となれば64bitLinuxをWoodcrestXeonで動かせばいいかな。
matlabだけならCPUはDualCore一個でいいんじゃない? ディスクはできればRAIDを組むとして、HPのProliant辺りでどうでしょ。

>>511さん
メモリ４GBに64bitは自分のメインマシンですがそれでも足りないんですねｗ

すみません。言い方を間違えたようです。
とりあえずコスト１０万ぐらいのレベルで考えたときどのような構成にすればよいでしょうか？

>>511さんのレスからするとメモリは多めでCPUはdual coreぐらいでも問題ないと。
ぜんぜん理解していなくて申し訳ないのですがmatlabは大容量のデータをいじることにポイントを絞り、それに伴うCPUがあればいい感じですか？


どうぞよろしくお願いします

そりゃあんた、>510の友達がmatlabをどう使うのか判らんから判らんがな。

>> 486
MATLABにコンパイラがあるのは本当。C++のコードに一度変換してコンパイルしている。
で、間違いのほうなんだけどある程度高速にはなるけどPureなC++でマシなコードを書いたほうが断然早い。
私の場合画像処理ライブラリ関係でお世話になったけど3割程度しか高速化しなかった。
ちなみにMATLAB7.1環境。

>>510
私は研究室のノートPCに入れて使っているけどストレスを感じませんよ。
余程に型落ちなPC(1G以下のCPU)でもない限りOSの推奨スペック程度あれば基本的には大丈夫かと。

スペックはあるに越したことはないけど、必要かどうかは、何の研究をやるかによる
計算量が爆発するようなアルゴリズムが必要な分野でないなら気にしなくて良い。

matlabの超初心者です。
既にあるプログラムを手探り状態で書き換えているのですが、
行き詰ってしまったので教えて下さい。

ボタンが２個とスライダーが１個あるフォーム(diologueform)を作りました。
そのフォームのmファイルは以下のとおりです。

function varargout = dialogueform(varargin)
・・・
省略
・・・
% Get default command line output from handles structure
varargout{1} = handles.output;

% The figure can be deleted now
delete(handles.figure1);

% --- Executes on button press in pushbutton1.
function pushbutton1_Callback(hObject, eventdata, handles)
handles.output = get(hObject,'String');　% 　←　この行は自分で追加
guidata(hObject, handles);
uiresume(handles.figure1);

% --- Executes on button press in pushbutton2.
function pushbutton2_Callback(hObject, eventdata, handles)
handles.output = get(hObject,'String');　% 　←　この行は自分で追加
guidata(hObject, handles);
uiresume(handles.figure1);

（続き）
% --- 以下はスライダーについての処理
function slider1_CreateFcn(hObject, eventdata, handles)
usewhitebg = 1;
if usewhitebg
set(hObject,'BackgroundColor',[.9 .9 .9]);
else
set(hObject,'BackgroundColor',get(0,'defaultUicontrolBackgroundColor'));
end

% --- Executes on slider movement.
function slider1_Callback(hObject, eventdata, handles)
slider1 = get(hObject,'Value');　% 　←　この行は自分で追加

・・・
省略
・・・


そして、他の関数のmファイルで、

result = dialogueform;

とすれば、
resultに押したボタンの名前が保存されます。
しかし、どうすればslider1の値（スライダーの値）を
他の関数のmファイルで参照できるのかが、
わかりません。

わかりづらい、文章で申し訳ございませんが、
とても困っているので教えてください。

まぁ>511は極端な例だけど、matlabで「待つ」ような処理をさせるのであれば2coreは欲しいね。
メモリはあるに越したことはないけど何の計算をするか次第かな。
ってことで、普通にCore2duoのPCを(メモリとディスクをけちらずに)組めばいいんじゃね?
GPUをNVIDIAの8800GT辺りにしておいて、CUDAを使ったライブラリを使える可能性を視野に入れてもいいし。

>>518
MATLABって、まだマルチCPUに完全に対応していなかったような気がするのだが、
どのバージョンから対応したのか教えてくれないか？

いや、matlabの結果を待っている間、他の作業もできるようにってこと。
matlab自身が対応しているかどうかが理由じゃないよ。
まぁ、Pen4でmatlab動かしてる間にデータの整理しようと思ったら結構かったるかったんでね。

>>517
１）slider1 をグローバル変数にするのが一番早いような木がする。
２）荒技だけど、handlesの構造体に変数を勝手に付け加えてそこに書き込むと言う手もある。
handles.slider1 = xxxx;

>>521
早速のご返答ありがとうございます！
複雑さを避けるために２）を試してみます。

>>521
handles.slider1 = get(hObject,'Value');
としたのですが、他の関数のmファイルで、
handlesを参照する方法がわかりません。。。
たびたびすみませんが、ご教授下さい。

自力で解決しました。ありがとうございました。

もう一度書きこみさせてください

位相特性を周波数で微分して、群遅延を求めたいのですが
「周波数で微分する」の部分をどう書けばいいのかわかりません

位相特性は求めることができているのですが・・・
どなたか助けてください
お願いします

位相特性が周波数なのでそれを微分する

この春から、学生はMATLABが18900円で買える件。
アップデートは無しだが、そこそこToolboxも付いてるし使えそうよ。

開いてみるまでPICのスレだと思ってた事はヒミツです

>>527
まだ届いてないけど注文しますた。
社会人になってからも株のシミュレーションで使えるし迷わず購入

以前に購入したMATLAB 2006b 学生版を新しく買う予定のPCに入れようと思っています。
OSがビスタなのですが、インストールできますでしょうか。
サイバネットのホームページを見ても分かりませんでした(探し方が甘いかもしれませんが)。
親切な方、教えてください。

InternetArchiveから2006年のFAQを一通り見てみたけど、vistaに関する記述はないな。

>>530
以下、僕の経験則だけど。
6.x はライセンスマネージャが Vista に対応していないのでインストールできない。
メインバージョンが 7 になって以降は問題なくインストールできてる。
従って、学生版でも R2006b なら問題なくインストールできると思われ。

余談だけど Vista で 6.x を動作させたい場合は
VMware や Virtual PC を wrapper にしてる。

>>531,532さん
ご回答ありがとうございます。
すごく助かりました。
おかげさまで安心してPC購入できます。

ところで、別の研究室でVistaにver.6.5(R13)をインストールしようとしててできない奴がいます。
532さんが書かれた要因だったようですね。

その場合、532さんが書かれている方法をとるか、XP以前のOSが入っているPCにインストールすることになるかと思います。

そこでもう一度質問させていただきたいと思います。
恥ずかしながら、「VMware」や「Virtual PCをwrapper」の内容が分かりません。

その設定を使うことでVISTAでもver6.xを使うことは可能なのでしょうか。
その設定は素人が手を出さないほうが良い内容でしょうか。
もし手を出してよさそうなのであれば、どのあたりを調べれば使い方が習得できますでしょうか。

聞いてばかりですみませんがよろしくお願いします。

>>532
＞VMware や Virtual PC を wrapper にしてる。
速度はあんまり落ちないってこと？

>>533

お役に立てて何より。

簡単に説明します。
VMware や Virtual PC は、Vista 上で別の OS を走らせるソフト、と考えて下さい。
以下、僕の踏んだ手順です。
1. Vista に Virtual PC をインストールする
2. その Virtual PC に Windows98 を導入する
3. その Windows98 に Matlab 6.x をインストール
4. めでたしめでたし。

Virtual PC は MS 謹製でフリー、動作はしっかりしてるけど速度は気持ち遅め。
VMware はこの手のソフトの始祖で歴史があり、様々な OS に対応してます。
各々使い方を知ってて損のないソフトなので、とりあえず触ってみるのがよいかも。
使用法に関して細かく説明しているサイトは山ほどあるので、ググってみて下さい。
ちなみに wrapper というのは、「ラップするもの」の意です。
Vista を 98 の皮でラップして、その環境下で 6.x を動作させるという文脈で
使ってます。good luck!


>>534

おそらく予想どおりの回答で申し訳ないけど、スペック次第。
上述のとおり VMware の方が速度的には勝っている模様だけど、
個人的には Virtual PC でもストレスなく使えてる。
まぁ試してみるのがよいよね。

MATLABとIDLってどっちが実行速度速いですか？

MATLAB学生版、科研費（学振）で買えるそうですよ。
と報告。

>>537
どうやって？
HP見ても学生個人が買わないといけないって書いてあるけど。

学振だったらいいんじゃないのかな？

ぐぐっても探せなかったので質問させてください。
OSのコマンドプロンプトから直接Matlabスクリプトは叩けないのでしょうか？
自分の調べた限りでは、みんなMatlabを立ち上げて、Matlabコマンドプロンプトからの
実行ですた。
CGIとかでWebからMatlabを動かせたらステキだと思ったので。
やっぱりMatlabを立ち上げないとムリなのでしょうか。

>>540
matlab -r foo.m

大学の講義でMATLABを使って、非線形方程式を自分で考えて解けという講義内容だったのですが、MATLABの使い方が全くわからないので誰か助けてください。

いいよ

>> 541
ありがとうございました！！
教えてもらったものを元に、いろいろやってみたら、できました！

>>537
それは問題になると思うよ。

MATLABで作った図をイラレで使いたいのですが
どうすればいいのでしょうか？

>>546
ファイル→エクスポート→ベクタを保持できるファイル形式(epsやemf)で保存
→おもむろにillustratorで開く
コマンドラインから画像を保存する場合はprint命令を使う

これでいいかい

>>547
フォントとか座標軸とか、勝手に書き換えられるのではないか？

matlabでオブジェクト指向プログラミングをしたいのですが,
どなたか勉強に良いサイトを教えていただけませんか？

http://www.mathworks.com/access/helpdesk_r13/help/techdoc/matlab_prog/ch14_o10.html

本家のサイトが意外によかったようです．
自己解決しました．

フーリエ変換のグラフで横軸がタイムフレーム、縦軸が周波数binで値が大きいほど濃いグラフ書きたいんですがどんな関数使えばいいでしょうか？

初歩的な質問で申し訳ないんですが、
matlabを開いた時に開かれるcurrent directoryの
初期フォルダを変更することって可能でしょうか？
わかる方いましたらよろしくお願いします。

matlab のショートカットアイコンを作っといて
そのプロパティの作業フォルダを変えるとか。
もしくは startup.m を作って cd(foo); とか。

ショートカットアイコンってなんですか?

その質問は、もう質問の体をした荒らしみたいなもの。
google もまともに使えんのに matlab で何をするつもりなのか、答え給え。

ぼくぁいつもマジレスなんだ。

たいしたことじゃないんですが、
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];
のようにしているのだけれど、別になんでもいいです。

伝達関数からステップ入力に対する
オーバーシュートや静定時間を求める関数ってありますか？

>>556
aを、微分（差分）してみるといいかもしれない

>>558
おぉおもしろいアイデアだ。
D = conv(a, [1 -1]);
starts = find(D == 1);
ends = find(D == -1) -1;
ind = [starts(:) ends(:)];
でいけました。うほっ。とんくそ。

>>557
stepinfo
ただしR2006a以上のControl System Toolbox

>>560
ありがとうございます。でもうちのじゃ使えなかったみたいです。
あきらめます。

過去にMATLABにLISPのmap関数の様なものがあるかと質問したものです．
その後，mexを作って解決しましたが，MATLABには標準で，arrayfunやcellfunといった同様の関数があることが分かりましたので一応報告しておきます．

ちゃんと聞けば教えてあげたのに。

すいません．コンソールで値の入力を要求して，その結果を変数に代入するような文はどう書くのでしょうか？
Please imput value of x :
と表示させてxの値を入力させたいのです．

inputですた

すいません．自己解決しました．

図をepsで保存してイラレで読み込むと
マイナスが見えなくなるのですが
これってどうしたらいいのですか？

feval(fnc,r)
でfncに'sin'とかは入れてsin(r)を計算できるのですが、
x^2などのxにrを入れて計算するのはどうすればいいですか。

>>568

feval(@(x) x^2, r)

matlabで多目的整数計画問題って解けるんですか？
fminimax()の変数が整数版みたいなやつが欲しいんです。
fminimax()の条件をいじって、変数を整数に縛ろうと思ったんですけど、
やり方もわかりませんし・・・

縦軸の値から横軸の値を表示することはできますか？
たとえば、y=2xのグラフをxを1〜10まで0.1刻みで描いたとします。
そのとき、y=5のときのx値を知りたいわけです。
刻み幅が大きく合致する値がない場合はどうなるのでしょう？

>>571
interp1dをつかったら？

>>571
推定するしかねーだべが。

matlabの関数の説明が日本語で書かれているサイトないでしょうか？
できればURL張ってくださいorz

matlab の中の人たちって何者？
特に toolbox とか数学的に難しいものもあるから、
それを実装できる知識を持っているが、
研究者やらずに matlab 開発しているとか良くわからない。
ph.d 持ちの人たちなのかな？

大丈夫、私も数学的知識はないけど数学関数の実装は(数学屋に書いてもらったアルゴリズムに基づいて)行っている。

>>574
「ＭＡＴＬＡＢ　ヘルプ」でぐぐれ

>>575
「ＭＡＴＬＡＢ　ヘルプ」でぐぐれ

Unix版とWin版で何か、違いはありますか？

変わんない

「ＭＡＴＬＡＢ　ヘルプ」でぐぐれ

配列のindexing をC/C++みたいに 0 始まりにすることってできますか?

「ＭＡＴＬＡＢ　ヘルプ」でぐぐれ

>>583
お前は首くくれ

>>582
できないと思う

>>584
カレーはククレ

金かけたくないなら、pylabいいよ

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
とエラーが出ます。
どう対処すればいいですか？誰か教えて下さい。

xaとyの数が違うんじゃないの？

plot2d(xa,y(1:length(x) ...

にしたら？

突然すいません。
MATLABでLeonhardの軌跡ってどのように描けますか？教えて下さい。

こんにちはいきなりですが質問です。
例えば100行100列のランダムな値をもつ行列から上位30点の要素を抽出するにはどうしたらいいでしょうか？
誰か教えてください

a=rand(100:100);
b=create_set(a);
c=b(length(b)-29:end)

>>592
素早い回答ありがとうございます
早速やってみます

質問です　次のようなアプリを作りたいのですがMATLABを使い始めて日が浅いので見当がつきません、
作り方・コードを教えて下さい　お願いします。

フォルダの中に「ＰＣ」等の２文字が書かれたbmp画像がN枚ある
一枚一枚を任意の場所で左右にカットする（ＰとＣが分かれるように自分で分ける線の位置を決める）
それぞれ○○○○○L.bmp　○○○○○R.bmpとリネームして同じフォルダに保存する
この作業をN回繰りかえし終了
図：http://www.dotup.org/uploda/www.dotup.org22131.jpg

あんまり甘やかすのは、アレな気がするけどこれでできるはず…
問題は、分割位置が固定されているのか、問題毎に違うのか…。
プログラムで分割位置を判別は簡単じゃないと思うけど

% 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

>>595
ありがとうございます

>>問題は、分割位置が固定されているのか、問題毎に違うのか…。
画像ごとに違う位置で分割したいので分割位置を一枚一枚その度に自分で決めれるようにしたいのです
説明が分かりにくかったかもしれません、↓実際のファイル内はこのような感じです
http://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} ); % 画像の読み込み

コード中で.pngとファイル名を指定しておきながら.bmp画像しか用意してなかったんだろ
そのまま実行するという暴挙に出てエラーが出ましたとか言ってしまう神経が以下略だぜ

>>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); % 画像をグレースケールに変換

えっと、答える私も私なんだけど少しだけ叱らせて。
ここは初心者に対して伝統的に厳しいスレではないけれど、それを考えてもあなたのレベルはちょっと酷い。
私も大した事は言えた立場ではないけれど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
何言ってるんですか、レポートにしろ卒研にしろ、抜ける手は抜くのが一番じゃないですか。

勉強するために大学行くんだろ。
就職のためなら専門学校いけよ。

頭は使わないと劣化していくよ
教授を見ていれば分かるだろ

>>602
何言ってるんですか、専門学校の就職率なんてカスですよ。

http://uproda11.2ch-library.com/src/11100794.zip.shtml
あぷろだ初めて使うんだけどこれでいいのかな？
pass:1111
削除キー:1111

使い方
・samplesフォルダに全ての画像を移す
・imageselect.mをmatlabから起動
・開始ボタンを押した後、画像の分割位置をクリックで選択
見た目に良く分からないかもしれないけどちゃんと分割された画像が保存されています。
手抜き仕事なのでバグあったら頑張れ

>>605
おまえさんのアホさに感動を覚える

あっはっは、出来が酷いのは手抜き中の手抜きだから勘弁してやってください。まぁ、赤の他人に対しては良くやったよね…
ってそういう意味じゃないか…
一応やっている分野の近い相手なので、何となく情が働きまして…
matlabのGUIプログラミングはVBなどと比べても決して難しい物ではないけれど、あのレベルで習得させてやらせるのは辛いしね

GUIにせずとも、マウス入力で区域指定できるのではないかと

filter関数がＣ言語でかけねぇ。
誰か作ってないかなぁ

どんなフィルター？平滑化なんかの２次元フィルターってこと？

matlab の filter 関数でしょ。

そうそう。matlabのbutter関数で係数を求めておいて、
それをつかってfilter(b,a,signal)って使うことのできる関数を
Cで作りたいわけ

簡単じゃん。

matlab での代表的なコーディング規約ってないの？
built-in 関数みてると大文字はおろか _ すら使ってないのがほとんどな気がするんだけど、そういうもん？
ちょっと、これじゃやりにくい。

昔のFORTRANの流れを汲んでるから、ひょっとすると関数名は６文字までだったりしてね。
個人的には _ を使った関数名は自作関数というぐあいに区別しているけど。

MATLAB(R2008a)で関数自体の引数の数を得る方法はありますか?

無名関数を引数としてその無名関数自体の引数の数を返す関数を作りたいのですが。
また，関数自体はまだオブジェクト化されていないのでしょうか？
その様なクラスも見当たらないですし。

関数自体ってどういう意味？引数の数はnarginだろ普通

>>617
>引数の数はnarginだろ普通
それはある関数呼び出しに関する引数の数だろ。

http://phpspot.org/blog/archives/2008/04/javascript_65.html
これで分かるかな?

無理じゃね？
matlab は数値演算はいいけど、そういう補助機能が弱いよな。
>関数自体はまだオブジェクト化されていないのでしょうか？
そんなことする気あるの？なんかで言ってた？

というか、MATLABをプログラム言語としてみるとこれ以上ダサい文法の言語は正直無いと言えると思う。
最初に行列を格好良く計算できる機構と入れたのは良いけど、やっぱり後で条件分岐やループが必要だと入れたあたりから破綻が始まっている

そこからは無計画と言うしかない文法の拡張が行われている
比較的最初に気づくものでも
関数作りたかったら別にmファイルを作ってそこに書かないといけないとか
関数は最後にendで括っても括らなくても良いとか
何故か決して使いやすいとは言えない関数名がC言語風とか
行列計算の書式しか無かった頃だけ直感的で利点があった、何故かインデックスが1から始まるとか

正直、簡潔な書式の行列計算＋強力なplot機能＋インタプリタですぐ結果が見える使いやすさなのだと思う
本当は文法関係は言語の進化に合わせて再設計が必要な時に来ているのだと思う。もう無理だけど…。
>>614
mlint関数で文法チェックとかどうでしょう？

スライドバーのコールバックで処理するとき、
ドラッグが終わってバーの位置が確定した時だけ処理を行うことはできないでしょうか?

バーが変化するごとにいちいち変化してると重いですorz

>>621
アルゴリズムの問題だと思う。
ちょっと考え直したら？

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次微分

宿題ぐらい自分でやろうぜ
つーか、ほとんどmatlab関係なくね？

FORTRANならわかるの？なのに matlab だとわからないの？

アルゴリズムが分からないのなら
教科書を読め

matlab の関数名、変数名は lower_case
のように小文字と _ を使っている人が多いようだけど、
なにか理由があるのかね？
なんとなくそういうスタイルのテキスト多いのかな？
http://datatool.com/downloads/matlab_style_guidelines.pdf
の人は変数名には lowerCamel を押しているようなんだけど。
＃関数名は _ もなしの lowercase 押してるけど、これはないよな。よみづらい。

matlab を起動すると勝手に IME が日本語入力モードになります。
matlab で日本語入力モードになる意味がわからないので、
止めたいのですが、何か手法はあるのでしょうか？

ここのスレ、性格悪い低脳多いな。

>>627
空白をあらわすのにlower caseを使う人がアメリカ流の人には多いと思う

matlabのmファイルからwavをsound関数を使って再生しているのだけど、
途中で止めることができない（；＿；）。

途中で止めることはできないのでしょうか？

コードは以下
[y, Fs] = wavread('test.wav');
sound(y, Fs);

>>630
lower_case と lowerCamel どっちがいいかなぁ・・・
lower_case は古い流行な気がしないでもないんですよね。
UNIX というか GNU C コーディングスタイル的な。

>>631
最初から止めたい位置がわかっているのなら、
マトリックスの部分指定。
とにかく止めたいのなら、CTRL+C

>>632
JavaScriptなんかの場合は最初は自分で書いたライブラリは_で区切るなんて風習もあったけど
そのうちライブラリを自分で配る人達が出てきてそれが普及してからはlowerCamelとJavaScriptの言語の標準に従う様になった。

最初は遊びで書いた関数も、そのうち人に配る可能性が無いわけでもないし、MATLAB標準に黙って従うのが一応正しい姿勢となるのでは？

MATLABの大文字小文字の区別って、完全に可能なんだっけ？

>>634
でも、matlab 標準というと、例えば bankrupt company が
bkrptcmpny みたいなかんじじゃない？
これはひどいでしょ？実際は、ここまでひどくはないけど、
_ もなにもないし、まぁひどい。

>>635
ファイル名 == 関数名、だから windows では大文字小文字が区別つかない。
けど、pca という関数がある時に、Pca を作るようなセンスの人はいないだろぉと思うし、
このために、lowerCamel は絶対ダメ！っていうのはもったいない気がする。

>>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 標準がどっちつかずなもんで、むきーという状況。

なんだか長くなって申し訳ないけど、
C も標準は、略称使いの _ もなしの lowercase と思っていいいと思うけど、
OpenCV とか OpenGL は lowerCamel を使っている。

個人的には lowerCamel を使いたいのだけれど、
ファイル名 == 関数名の所で悩んでいる。
関数名は lower_case のような形で書いて、変数名は lowerCamel
ってアリなセンスなのかな？

ところで、private 関数的なネスト関数の名前は最初に _ でも付けようか、
と思ったら matlab では変数、関数名の最初を _ にはできないのね。がっかり。

お前プログラム書くの遅いだろ

>というわけで、matlab 標準がどっちつかずなもんで、むきーという状況。
で、個人的意見なのですが、基本的にはMATLABのユーザーは物理系の人で、プログラミングの専門家やほかの言語に詳しい人は少ないのでそういうリテラシは少し低めなのではないでしょうか？
OpenCV とか OpenGLとかはC++の文化を意識してああいう風になっているんでしょうねぇ
オブジェクト指向を基本的には出来ないとはいえ、一応matlabの関数名の文化は、Cからの影響なのを考えるとlowerCamelの記法は可なのではないでしょうか？

> ところで、private 関数的なネスト関数
Lispでは「高階関数」と呼んでいるのでこれに合わせると良いのでは？

>>640
>lowerCamelの記法は可なのではないでしょうか？
やっぱり問題は関数名なんだよね

>Lispでは「高階関数」と呼んでいるので
高階関数ってラムダ関数ちゅうか関数ポインタを受け取ったり
返したりって関数でしょ？fminsearch みたいな。
ネスト関数ってのはこっちね。
ttp://www.mathworks.com/access/helpdesk/help/techdoc/matlab_prog/f4-39683.html

速く走らせるには工夫が要るけど
それほど大きなパッケージを作るわけでもなかったら、
書法は気にしなくてもいいのでは、って思う。
おれはそのレベルのユーザー。

z = (2*cos(x.)+1)(2*cos(y.)+1)
のグラフを書きたいのですが、上記の式で
エラー: 一致していないか、予期せぬ括弧か大括弧です。
とでます。x,yの値は
x=-3:0.1:3;
y=-3:0.1:3;
このように設定しております。
どなたか教えてください。

すみません。.の意味を取り違えてました。

質問があります


エクセルデータを読んでグラフにしたいのですが、ファイルネームが
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行目にどうやっても忠告がでてしまいます。
文法うんぬん以前に考え方自体が間違ってるのでしょうか？

こんなシンプルなのに考え方もなんも無いって。単純に文法の問題。
「こう書けるだろう」と勘で書いたところがあるでしょ?
そういうところが本当に正しいかどうかチェックしてけばいい。
help strings

やっぱりそうですよね。
プログラム自体が始めてから数週間なのでグダグダです。
ありがとうございました

アイディアとしてはいい線いっているので、
あと少しだけ、知れば、サクサク作業できるようになると思う。

すいません、ぐぐっても情報が少なすぎたので質問させてください。
JMatLinkというJAVAとMatlabをつなげるツールをインストールしたいのですが、
http://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を使用している方がいましたら、
何が悪いのか教えて頂けないでしょうか。

Matlab 2007bと、Viual Studio 2005の質問です。
既存のC++で作られたライブラリの関数をラップして
（要は、mdlOutputs内で既存の関数をコールする）、
C++ MEX S-functionを作ろうとしています。
Visual C++ 2005上で作成／ビルドして、*.mexw32を作っているのですが、
simulink上でこれを実行しようとすると、
「無効なMEX-ファイル*****; 指定されたモジュールが見つかりません。」
というメッセージが出ます。

サイバネットが出している、C-MEX-S-ファンクションチュートリアル
（http://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やらをリンクしているモジュールの
コンパイルは、とても難しそうです。

>>647
>>646に書いてあるように help strings すれば
必要なことが書いてるよ。

MATLABはちょっとクセがあるから他の言語を知ってる人が
いろいろひっかかるポイントがあると思うけど
ヘルプはけっこうちゃんとしてると思うから
どんどんヘルプ読んでいけばすぐ覚えるよ。

>>616

ezplotとかだと関数じたいを引数として渡してるよね。
無名関数が可能なのかは知らない（たぶん出来ない）。

help function_handle

無名関数は知ってるけど、
その引数の数を調べられないか？ってのが大事なんじゃろ。
で、matlabの引数って
function hoge(a,b,c)
と書いてあっても
hoge(a,b,c,d)
と呼べるっちゃよべるし、呼び出しの時点になるまでわからないのでは富田

matlab.exe のオプションってどう見るの？
--help でも matlab 起動するし。

>>654
>>652に答えてるのにお前が知ってたらなんなんだよ。。
だいたい呼び出しの時点の引数の数を呼び出す前に知る方法があるわけないし。
何言ってるのかわからんが。

>>616で書いてるのは定義にある引数の数だろ。だったらnargin。

[X,Y,Z] = peaks(30);
surfc(X,Y,Z)
で出る図を png で保存すると、
真っ黒画像が作られるのですが、
みなさんはきちんと保存できたりしますか？
ぐぐってもだれも surf がおかしいとか言っていないようで。

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以外にもっと良い方法があれば
教えてもらえると助かります。

別ディレクトリの関数を addpath せずに実行する方法ってないでしょうか？
その関数を実行して、中で一連の addpath をさせたいのです。

>>658
とりあえず sprintf のほうが見やすそう。

>>657
オプションを設定

>>658
エラーの内容そのまま。
１行づつ、確認しろ。

レスくれた人達ありがとう。

うおお、下から３行目の num2str(2/i) は num2str(i/2) とすべきでしたね。evalの性質の理解不足
とか以前の問題でしたorz　あと、sprintfを提案してくれた人もありがとう。これからいろいろ試してみます。

>>661
オプションって surf のオプション？どれ？

昨日はどうも。もう一つだけ質問させて下さい。
>>658とは逆のケースで、30個のエディットテキストにユーザーが入力した値を、15*2サイズの
行列Aとして格納することを考えています。

guideでエディットテキストを配置すると、対応するm-ファイルにCallback関数が30個自動的に
生成されるので、それぞれの下部に各行列成分への格納文を書けばいいわけですが、やはり非効率的
なので、もっと良い方法を探しています。

この場合、これら30個の自動生成された関数宣言文を削除し、evalなりfevalなりを使って、30個の
関数を変数を使ってまとめて定義するような文を手で打ってもいいんでしょうか？　確かこの手の、
「guideを経由しないでm-ファイルに直打ち」はNGのケースがある、と何かの本で見たことがある
んですが。

あるいは、このような場合での常套手段が他にありましたら教えて下さい。

>>659
スクリプトでいいならrunでできるけど
関数を呼ぶにははそれと組合せて
function handleを使うとかかな。

それかベタだけど
wd=cd(dir); func; cd(wd);
とか。

>>666
addpath をしたいだけなので、
スクリプトにして、run で呼べばよさそうですね。ありがとう。

ランダムな行列の中から局所的な最大値とその座標を求めたいのですが、
何かいい方法がありますでしょうか？
max関数を２回使うことで、大域的な最大値は求まるのですが。。
局所的な最大値を求めるアイデアがどうしても出てきません。。
お願いします。

>>668
大域的なmax値をNaNで置換して、複数回走らせる
あるいは
元データを微分する

>>668
局所的な最大値って極大値でいいのか?
いずれにしても「ランダムな行列」がどういうものであるかと
何を「局所的な最大値」と定義するかが決まらないと
答えられないのでは?

レスありがとうございます。

具体的に書くべきでした。
例えば、
1 4 6 8 4 6
4 5 1 9 0 7
4 0 9 4 1 6
4 7 5 9 0 2
とあった場合、３列目の９と４列目の２つの９、６列目の７が最大値（＝極大値）
と考えたときです。
お願いします。

座標が整数値でない画像データを可視化したいんですが、どうやったらいいですか？
データの形式としては(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は行数）のベクトルデータとなっています。
これを穴があってもいいから通常の画像のように可視化したいんです。

点ごとにプロットする色を変えたらいいのかと思ったんですがうまくいかなくて

>>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 れ。

ある関数Ｆの最適化問題を解きたいんだけど、Ｆがゼロになった瞬間に
最適化を止めて次のことをさせたいような場合、どうしたらいいですか？
最適化にはＮｅｌｄｅｒ−Ｍｅａｄを使いたいので、Ｆｍｉｎｕｎｃより
Ｆｍｉｎｓｅａｒｃｈを使いたいです。

Ｆの値は定義上０より小さくなることはないので、一度０がでてからさらに
パラメター空間でサーチするのは単なる時間の無題です。

Ｆｍｉｎｓｅａｒｃｈのコードがわかれば外にＩＦ節をつけれてＦの値が
０になったところで止めさせられるものの、どうしていいかわかりません。

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]
が返ってほしい。

>>673
即レスありがとうございます
ちなみに正規化ってどうやったらいいんでしょうか
xの列でソート→代入、yの列でソート→代入じゃ効率悪すぎですよね
手早くやれる方法はありますか？

１から１００の整数をランダムに並び替えるプログラムを教えてください！

>>677
宿題は自分でやれ

>>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

すみません、どなたか>>665が分かる方いらっしゃいますか？
本質的には、いかにしてCallBack関数をまとめて設定するか、ということだと思います。

>>673
matlab centralにあるということでしょうか？
探してみましたが、わかりませんでした。。
もしよければどれか教えてもらえないでしょうか？
お願いします。

えー、なんで？極大値は英語にしたら extreme か extrema だよ？
extrema.m

>>682
完全に見落としてました。。
これでいけそうです。
ありがとうございました！

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文を使わないでコロンを使った形式にできるものなのでしょうか？
よろしくお願いします。

すみません、間違えてタブを使ってしまったので投稿しなおします
******

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文を使わないでコロンを使った形式にできるものなのでしょうか？
よろしくお願いします。

スペースもダメみたいですね…みにくくてすみません。
よろしくお願いします。

repmat を駆使して４次元配列にまでもっていけばきっとできる。
で、あとは min で argmin うけとる。

>>687
ご指摘いただいた方法でできました。
ありがとうございます。

>>684って行列の中身が一緒じゃないんでしょ？
だったらrepmatじゃだめじゃね？
ってか俺も今
for a=1:10
for b=1:5
kotae(a,b)=a+b
end
end
って感じのやってるけど、これfor文じゃなくせるの？

matlab centralで拾ったファイルを動かそうと思ったんですが、
「500の最大の再帰呼出の制限達しました。制限を変更するには、・・・」
となって動きませんでした。
制限を2000あたりまで変更しても同じ結果で。3000くらいになるとMatlabが
強制終了してしまいました。
これ何か対処法ないでしょうか。。？？

mファイルの作成者に聞くべきじゃない？
つか、少しは自分でk

>>689
なくせう。for 文使うとゲロのろいぞ matlab。

コロンでできるの？

できる。
kotaeって10x5のマトリックスになるだろ？

>>689
meshgridというのを使うとてっとり早いです。
行列にして掛け算する、ということです。
687さんがおっしゃった通り、repmatを使ってもできます。
というか、repmatを使ってmeshgridが実装されてるようなイメージです。

はじめまして。
ＭＡＴＬＡＢをはじめたばかりの初心者なのですが、連番をつけるところで悩んでおります。
文字列の配列で0001、0002、.....　、9999のように整数で任意の桁数にしたいのですが
うまい方法はありませんでしょうか？

エクセルでやっていたときは、上記の場合であれば”000”という文字列と整数の値を結合し
左から4文字を切り出して使用しておりました。

ＭＡＴＬＡＢでもエクセルでやっていたような力技でやるしかないのでしょうか？

もしご存知の方がおりましたら、ご教示いただけませんでしょうか？
よろしくお願いいたします。

申し訳ございません。
上記のレスで左から切り出していたというのは誤りで、右から切り出しておりました。

>>697
ttp://dl.cybernet.co.jp/matlab/support/manual/r13/toolbox/matlab/ref/?/matlab/support/manual/r13/toolbox/matlab/ref/sprintf.shtml

>>698

ありがとうございます。

大変申し訳ございません。私には、このページではsprintfでたとえば”1”を”00001”と定義する方法が理解できませんでした。
整数を5桁で表して、それを文字型にnum2strか何かで文字列に変換するようなことができるということでしょうか？

私の思い付きではstrcatをつかって696に書いたような力技で書く以外に方法がわかりません。

お手数をおかけしますが、ご教示いただいたページのどれを使用すればよいか教えていただけませんでしょうか？

>>699
printf フォーマット
でぐぐればいろいろ解説があるよ。
Ｃ言語が元だけど自由度がかなり高いから覚えとくとよい

↓とりあえず解決策
str=sprintf('%05d',i)

>>700

取り合えずというか、そのものずばり私の欲しかった回答を書いていただいて、感謝の限りです。

お忙しい中、どうもありがとうございました。
これからも精進いたします。

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]
が返ってほしい。

findじゃだめ？

俺の力じゃ5行が限界だった
x=find(A==1);
y=find(A==2)-1;
z=intersect(x,y);
B=zeros(1,length(A));
B(z)=1;
もっと簡単なのきぼん

面白そうなので参加

((A(1:length(A)-1)+A(2:length(A))).*(A(1:length(A)-1)==1))==3

正の値で且つ後ろ一つは無視する前提
一般性があるかどうかは知らない

size(A)も使ってみよう

A(1:end-1)==1&A(2:end)==2

レベル低いなあ。
こんなのすらすっと出てこないってことは
よっぽど効率悪いコード書いてるんだろうなお前ら。

煽るのは構わんがそれ全部ゼロになるぞ

って１ズラしてるか
スマン

MATLABの関数自体のコードを読んでみると
勉強になるんだけどな。

質問です。

MATLAB上で画像処理をする際に
画像のある画素(座標)を保存して
さらにその保存した座標のｘ方向成分だけ読み出すには
どういったコマンドを入力すればいいんですか？

>>711
まず入門サイトでも読んだ方がはやくね？ってかはやい

2008からオブジェクト指向が充実したみたいで、
それについて詳しく解説してある本を知っていたら
誰か紹介していただけないでしょうか。
よろしくお願いします

誰か、MATLAB(2007b)で３次元の楕円球体の描き方教えてください。

定義したい変数名を、文字列を使って定義することは可能なのでしょうか？
例えば、
Data1 = 1;
Data2 = 2;
・
・
・
Data30 = 30;
と30個の変数を定義したいとき、for文を使って定義することは可能でしょうか？

>>714
キーワードを英語でlookforするくらいはしようよ。

>>715
help eval

はじめまして、今matlabの勉強をしている者です。

matlabのエディットウィンドウで、文字の挿入ってどうやるんですか？
現状では入力してある行に文字を挿入すると後ろの文字が消えてしまい、
全文入力し直しになってしまうのです、これをwordなどみたいに
文字を挿入したら、文が自動で後ろに移動してくれるようにできないですか？

また、同様にスペースキーも文字間のスペースをあけてくれず
文字を消してしまうのですが？どうしたらいいですか？

OSは vinelinux4
matalab v7.0.1 SP1です

MATLABっていうか
キーボードのinsert

GUIを作成して、
ボタンを押すとシミュレーションが実行するようにしたんですが、
GUIなしでコマンドから直接実行した場合だと１６０秒くらいで終わるのに、
なぜかGUIからだと２６０秒もかかるんですが、
なんでこんなにかかると思いますか？
やってることはまったく変わらないですが・・・。

>>719
全然見当がつかん。
プロファイラを実行した状態で GUI 経由でシミュレーションを行って、
時間のかかっている箇所を割り出してから相談してもらった方がよいと思われ。

>>720
アドバイスありがとう。
実はこういう経験あまりしてなかったんで、
プロファイラとか使ったことなかったんです。
ちょっと調べて割り出してみます。

非常に低レベルなことで申し訳ないのですが、
Cでいうポインタのような概念はMatlabにはないのでしょうか？

>>722
ないはずだよ
wikipediaにも書いてある
まあ、google codeとかmatlab centralで検索すれば
自作のポインタ関数出て来ても可笑しくないけど

>>722
デフォはポインタ
渡された変数を編集しようとした時点でコピーが始まる

>>722
何をしたいのかをはっきりさせて
それをMATLABではどうやるのか、と考える方がいいよ。

他の言語のやり方をそのまま持って来ようとしても
そうは行かない場合は多い。MATLABに限らず。

質問です。

例えば
3,6,4,5,2
4,15,7,5,3
6,4,4,7,3

のような行列があってノイズで15だけ値が大きい場合には、
それを周りの数字に近い値に矯正（たしか平滑？）するにはどうすればいいですか？
自分が知ってるやり方だと、ノイズも含め全ての成分が平均化されてしまいます。
どうにかノイズ以外はそのままの値で、ノイズ部分だけを周りの値に馴染ませたいんです。

10はノイズ? 13はノイズ?
15を適当に5とかにしたら「馴染んだ」ことになるの?

自分が何をやりたいかをまずは具体的にする。
君以外にそれがわかる人はここにはいないよ。

実験データの改竄捏造でもしたいのかなｗ

>>726
if ノイズ
   平滑
else
   なにもしない

　ああ　

>>726
新しいやり方を考えればいいだけ。

質問です。

Ｃのプログラム等に記載するための離散時間系のデジタル１次ローパスフィルタを作るためにＭＡＴＬＡＢでフィルタの係数を求めようとしています。

そのためのプログラムとして、サンプリングインターバル＝６ｍｓ、フィルタ周波数＝１０Ｈｚとして、
　　　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
と言いますと・・・？

周波数とインターバル（周期）とで混乱していない？
間違いを探したいのなら、
他の数字の組み合わせにしてみて、計算してみるのはどうでしょう？

>>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

どなたかご教授いただけると助かります

>>736
サンプリングインターバルも変えてみたらどう？

変えてみた結果は以下の通りになりました。

　　　　　　　　　　　　　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)
という見方でよろしいのでしょうか？明らかに異なった波形が出てくるので、離散系の場合はボード線図を違う見方で見なければいけないのかなと予想しているのですが。。。
見方がわかりませんので、もしご存知でしたら教えてください。

ttp://www.ikko.k.hosei.ac.jp/joomla1506/matlab-topmenu-33/17-matlab20070307.html

matlabをはじめたばかりなんですが質問です。
stem3でつくったグラフを等高線のグラフにそのまま変えたいのですが、
どなたかわかる方いらっしゃいましたらお願いします。

ステムを輪切りにして上から見て、意味あるの？

>741
値がどのように分布しているか等高線の方がわかりやすいと思い
変えたいのですが。

>>742
matlab 等高線
でぐぐれ

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);

ではできませんよね?

効率はわかりませんが

function out = gui2(handles)
handles.a=2;
out = handles;

として、呼び出し側でguidata(hObject, handles); することで解決しました。

ベルヌーイの螺旋を描ける方いらっしゃいませんか？
困っています。

描けるよ

>>745
構造体使う必要がないならglobal宣言すれば簡単かと

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)
----------------------

としたのですが、うまくいきませんでした。
おそらく行列のプロットを理解していない
と思うのですが・・・。
非常に低レベルな質問で申し訳ありませんが、
よろしくお願いします。

mesh(z)でどうだ

>>750 様

早速のご回答、ありがとうございます！
無事、3D plot をすることができました。
ありがとうございます！

MATLABで式の簡単にする変形を一ステップずつやって
どうやって、変形して最終の形にするのか
見る方法ありませすか？

意味がわからない。
mathematica？

シンプリファイとかやるときの途中の式がしりたい。
シンボリックマスツールズの

'"2008-09-23 02:30:00.1"'　というフォーマットで書かれたストリングのセルから時間データを抽出するよいアイディアをどなたかご存知でないでしょうか？

もしも、''がなく、"2008-09-23 02:30:00.1"　というフォーマットならば、sscanf(%d-%d-%d %d:%d:%d);　によって時間データを抽出できるのですが。。

日付の区切りがスラッシュにかわったりしただけでも
どうすればいいかわからなくなるのかなこういう人は・・。

だれか 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
と出る。なにが足りないんでしょうかね？

愛

>>757
回答にはなってないかもしれんけど、
おいらは、加えてBuilder for .NETを使ってMファイルを.NETコンポーネント化
してから、.NET系言語から、そのコンポーネント呼び出して使ってるよ。

A=zeros(100000,100000);
だけでメモリ不足になります

余分な変数消去や仮想メモリ等をいじってみましたが効果はありませんでした
これを解消するには結局メモリ増設しかないのでしょうか？

>>760
メモリ増設しても無理じゃね?
100000*100000って、数値だけで10Gになるからアクセスすることさえ厄介になってくると思う。
メモリ量的には更に数倍になるわけで、オンメモリで処理するのはどうせ無理だし。

>>761さん

ご回答ありがとうございます
低レベルな質問ですみませんでした

このような行列はMATLABではそもそも扱えない（普通は扱わない）と考えてよいのでしょうか？
株価のシミュレーションでforを使うより行列を使った方が早いため
あのような行列が処理できると助かったのですが

「普通のＰＣ」では扱えない　というのが正解かと。
１０Ｇのメモリっていうといわゆるスーパーコンピューターの領域

たとえば*を増やしていって、
処理の進行状況を表示したいのですが、
改行を抑制する方法はありませんか？

質問です
複素数の実部虚部それぞれ七桁で表示させたいのですが
zに複素数が入っていて
sprintf('%7f',z)
とすると実部だけが七桁で表示されてしまいます
どうしたらよいでしょうか？

関数realとimagを使用して実部と虚部を取り出せばいいじゃない

サポート受けられないということは、ここの住人はコピーや労か

>>764
なんで今どうやってるのか書かないの?

>>766
ありがとうございます
sprintf('%7f %+7f',real(z) imag(z)i)
のような感じでしょうか？
今試せる環境ではないのでiとかがどうなるのかよく分かりませんが・・・・・

教えてください。
１辺の長さが１の正三角形を上下に並べて出来る平行四辺形をプロットして、
それに含まれるかどうかの判定をしたいのですが、判定条件はどうしたらいいでしょうか。
N個の点を行列で与えて、図形に入っていれば赤、入っていなければ黒という風に
表示させたいのです。
よろしくお願いします

>>769
+の位置が不適切だけどそれを直せばイケル

>>771
すいませんまだ試せていないんですがこうすると実数値の場合に
3+0i
みたいな表示になっちゃいませんかね？
もっと簡単にできる方法ないですかね

質問があります


例えばデータシートが1000枚あって、各データシートからX,Y,Zの値を取り出して、
座標(X,Y,Z)をplot3でプロットします。これがこのシートの3次元での点（位置）になりますよね。
それを1000枚分やると3次元空間に点が1000点でき、動きの軌道のように繋がって表示されます（そういうデータなので）。

問題は次なのですが、数十点ほど動きの軌道からはみ出してプロットされた点（多分計測時のノイズ）が
存在します。なのでプロット点全体を平滑化したいのですが、三次元プロットした点をフィルタするのは
可能なんでしょうか。

よろしくお願いします

>>773
成分ごとにやれば？

結局の所作業としてはプロットした点をフィルタせずに
フィルタしてからプロットするんでしょ?

Matlab C ライブラリの関数
mlfMatrixPrint　とか
mlfLength とかってどこに行ったの？
昔は matlab.h ってのがあったみたいなんだが、
R2007b だけど、ないんだよね。
matlab 内のビルトイン関数を C から直接呼び出せなくなったの？

私は>>675じゃないんですけど
>>675の解決策ってありますか?

スレを675以降もネットリと見てみよう

こんばんは。最近MATLABで動画像処理を始めたものですが、
250*350でフレーム数200のAVIの動画をRGBで格納するために
T=zeros(250,350,3,200)
としてみたのですが、メモリが足りませんと出てしまいました。
解決策がわかりません

ちなみにPCはXP　Core(TM)2 Duo CPU　メモリ2Gを使っています。

よろしくお願いします

[A;B]という行列(n行m列)から [A;A;B;B]という行列(n^2行m列)をfor文を
使わずに 作りたいのですが可能でしょうか。
基礎的な質問でしょうがよろしくお願いします。

>>780
[A;A;B;B]だったらn*2行m列だと思うんですが。

例が分かりづらくて済みません。別の例ですが[A;B;C]という行列だったら
[A;A;A;B;B;B;C;C;C]という行列を作りたいのです。

>>779
メモリが足りないのものはしょうがないのである。
本当に全部格納しなければいけないのか考えてみるのである

>>782
つまり *3 でいいんだね？
X = [A;B;C] からA, B, C の分解は終わっているものとする。
Y = [repmat(A, 3, 1); repmat(B, 3, 1); repmat(C, 3, 1)]

>>784
ありがとうございます。確かに直接記述すればそうなるんですが、
行列の行数が任意の数nであってもそのような行列が作れるか
、つまり各行ベクトルをn個並べてn^2行の行列を作れるかということが
分からないんです。

matlab Mex で
X(1,:)
ってどうやんの？

>>785
repmat(x, 1, length(x))

>>787
おかげで解決しましたーありがとうございます。

[id,check] = textread('01.txt','%s%s');
と読み込んだ後
id 　とすると中身があるのは確認できるのですが
fprintf('filename','%s',id);
というふうにして書き込もうとするとFunction is not defined for 'cell' inputs.
というエラーがでてidをtxtに書き込むことができません。

textreadで読み込んだidをtxtに書き込むにはどのようにすればよいでしょうか？
よろしくお願いします。

>>789
http://www.excite.co.jp/world/

>>789エラーの内容は分かったんですけれど
セルをtxtに書き込めるようにする方法が調べてみてもよく分かりません
もしよかったら教えてもらえませんか？
よろしくお願いします。

MATRABのプログラムで、ｘ１＝０、ｘ２＝１．２ ｘ３＝２．１ ｘ４＝３．５ ｙ１＝０．３ ｙ２＝１．２ ｙ３＝２．４ ｙ４＝２．７ の４点のラグランジュ補間をとりグラフを出力したいのですが教えてください
。できればfor文をつかった形式でお願いします。あと範囲は０＜ｘ＜１０です。

レベル低過ぎ
説明書読め

MATLABわかる人
http://pc11.2ch.net/test/read.cgi/tech/1226576323/l50
スレからのマルチポスト

MATLABのプログラムで、xy平面上の点P(i)=(x(i),y(i)), i=1,2,3,4,5について
５個の点P(i)を適当に設定し、これらの点を順に線分で結んでできる多角形の図を作れ。ただし各点P(i)は*印で示すこと
最後の点P(5)を最初の点P(1)に結んでみよ。
のプログラムを作りたいんですけど、誰かわかる方教えてください。

>>792
plot([1:10], lagrange(x, y));

>>795
line(x, y, '*');

宿題は宿題スレに頼めよ

宿題スレのURL教えて！

ｼｮﾎﾞｰﾝ(´・ω・｀)

>>770
まはやMATLABですらねぇな

三角形の中に点があるかどうかは時計回り（または反時計回り）に
各辺を回ったときに常に辺から左側（または右側）に点があるかどうかでわかる

N角形〜閉曲線の場合でも一緒

MATLABで2の5乗根の値を出したいのですが何かいい関数がありませんか？
よろしくお願いします

2^(1/5)

ありがとうございます。

質問です。
MATLAB5.0でフォルダ内のファイル全てに同じ処理を施すようにしたいのですが、
どのようにしたらいいでしょうか？

サンプルプログラムを書いていただけたらうれしいです。

自分でかけよ

こんばんは。
MATLABで近似曲線を描きたい場合、データがXとYの２軸のグラフでプロットされてるなら、
近似曲線用のコマンドがあるのでできます。しかし今回はXYZの三軸の要素があるデータをplot3コマンドで
グラフ化したものに近似曲線をひきたいと考えていて、方法をずっと考えていたのですが全く思い浮かびません。

どなたかアイデアがありませんか？

>>803
name=dir
-------
gnuplotを、Windows 2000のOctave Ver3.1.50 Winows で使いたい、どうすれば？

工学系の学生です。質問させてください。
常微分方程式の境界値問題用の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

上の本、高いですね−。

質問です。
plotで線をつなげずにヒストグラムのように描くにはどうすればいいですか？
（データは600行ほどのテキストファイルで、それぞれの行に数値が入っています。
それをdlmreadで読み込んでplotしたら横の軸は600までで思い通りなのですが
それぞれのデータ点がつながってplotされてしまいました）

hist()があるよ

アッコにおまかせでやっていたけど、
古物商の免許を取って古本屋で、売れ残りかつ専門的な本を安値で買って
アマゾンで売るという商売で月にだいたい１０万円、多いときで３０万円
の商売している人がいるらしいよ。
>>807の紹介した本もそういう人が提供している可能性大だね。
ではではノシ

アッコにおまかせで微分方程式 まで読んだ。

809さん、hist()をうまく使えばいいんですね！ありがとうございます。

MATLAB R2008b studentバージョン
を買いました。
mccを実行しようとしたところ、mccは用意されてなかったみたいです。
studentバージョンだからですか？

matlabの正式バージョンってマジでいくらするんですか？

>>814
MATLABだけ？それとも、Simulinkも？
まぁ、国産の軽〜高級車くらい見ておけばいいんでない？

随分幅のある例を出してきたな。10倍以上あるじゃないか。

ライセンスは細かく分かれているから、組み合わせでかなり変わるしな。

今日、MATLAB EXPOやってるな
聞いてくれば？

ある計算で、メモリ不足に陥っています。
512ＭＢ×４のメモリから、２Ｇ×２のメモリに増設しました。
zerosでどれだけのサイズを確保できるかを試したんですが、増設しても変わりありませんでした・・・。
10000×10000のゼロ行列でメモリ不足に陥ります。
サイズ増やせると期待したのですが、こんなもんなんでしょうか。
問題の計算のほうも改善できませんでした。どこかおかしいところがあればアドバイスお願いします。
バージョンはR2006b、ＯＳは32bitＸＰです。

10000x10000x8bytes=80,000,000
100Mそこそこだから確保は出来るはずだが、他にも巨大な変数があったら不明
メモリが断片化してる可能性もあるのでpackコマンドを使ってみる

それでも足りなければ↓に従い64bitOSを使うなどの対処を。
http://www.cybernet.co.jp/matlab/support/techkwdb/search_view.php?id=1628

覆面算をエクセルのマクロで解きたいのですが教えてください
Ｂ、Ｃ、Ｄは0を除外

　　　　Ｃ　Ａ　Ｂ
　+）　 Ｂ　Ａ　Ｃ
----------------------
　　Ｄ　Ａ　Ｂ　Ａ

>>820
800M だべ？
>>819
他のソフトも入れて (matlab 自体もいれて)、
32bit CPU && Windows なら 2GB の壁があるから、
そんなもんじゃねーがの。

>>882
おおう
単純な計算ミスorz
仮想メモリもふやして一応matlabが2G使えるはずだけど800Mなら他の変数しだいか。
やっぱ64bitでやるとかコマンドラインモードでやるとかが必要かもねえ

>>820>>823
有難う御座います。
clearやって不要な変数などは一切消去して試した結果です。
因みにタスクマネージャーでmatlabのメモリ使用量みたら700Ｍ程度でした。
ウィンドウとか最小化すると200Ｍまで減ったりしますが、こんなもんでしょうか。。

やっぱり64ビット環境かぁ。。。
その場合、ゼロ行列ではどれだけのサイズが確保可能になるんでしょうか？

変数名を文字列に変換することって出来ないのでしょうか？
例えば、行列Aのプロット結果のtitleに「plot of A」とか書きたいんだけど、数が多くて手入力では限界が・・・

その変数をシングルクォートではさめば文字列になるよ。

>>825
文字列から変数っつーか、式に eval で変換する、ほうが楽かも

>>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)
とすることで、自動でタイトルを付けることが出来ると思ったのです。
そういう関数とかないかなと探しているのですが、なんだかなさそうな気配です・・・

title('func')だと自動対応できない、というが
その例で変数名が変わったとき
plot(func)は対応できるの?

君自身、自分が何をやりたいかが
分かってないのではないか。

>>828
だから、変数名を格納して、eval で展開すれば、というとるやんけ

eval使わなくてもMATLABには>825の目的の機能があるよ。
825が何やりたいかちゃんと説明できてないから予想だけど。

>>831
その予想を書けよ。言われる前に書けよ。

まともに質問できない人には質問になってないよって
教えてあげるのが親切ってもんだよ。

不評を買ってしまったようで・・・
実はプログラミング自体にあまり慣れていないもので、上手く伝えることが出来ないみたいです。
アドバイスを元にもうちょっと試行錯誤してみようと思います。

>>829
そこは泥臭い手法を採ってまして・・・
あまりに下手なので説明は遠慮させていただきます。
泥臭くて良いのなら引数で渡せばいいじゃんとか言われるとそれまでなのですが、
機能として有りそうにも拘わらず見つけることが出来なかったため質問した次第です。

>>830
ありがとうございます。調べて試してみます。

>>831
説明が上手くなくて申し訳ないです。

>>834
>機能として有りそうにも拘わらず見つけることが出来なかったため質問した次第です。
そんな機能必要ないから存在しないだけだから

プログラムに慣れてないとかいう問題ではたぶんなくて
自分がやってることとやりたいことを全部書いてないだけだと思うよ。

>825
> 変数名を文字列に変換する
をしてくれるvarnameという関数が仮に存在したとしても
あなたの質問の例では意味がない。

func = [1 2 3];
plot(func)
title('func')

を

func = [1 2 3];
plot(func)
title(varname(func))

と書けても「手入力」の節約になってない。
plotとtitleの部分を関数にしたいということなら分かるけど。

help inputname

そう。だから文字列配列を作って eval で展開するのさ。
funcnames = { 'func1', 'func2', 'func3' };
for i = 1:length(funcs)
　eval(sprintf('plot(%s)', funcnames{i]));
　title(funcnames{i]);
end

>>831
eval 使わないとなるとどうやるんだ？参考までに教えてくれ。

なんで>>837はえらそうなのか。

> eval 使わないとなるとどうやるんだ？

>>836に書いてあるんだけど。

evalもいいけどデバッグしにくくなるし
よほどの理由がある時以外は使わない方がいいよ。

>>839
書いてないよ
>varnameという関数が仮に存在したとして

> help inputname

>>841
は？inputname は全然ちゃうやろ。help みたんかいな

頭悪いのにえらそうな人ってやっかいだね。
>>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

でしょ。

あ、テストせずに書いたけど
plotfuncの中にfigure追加しないとだめですね。
もとの質問者がsubplotとか使ってるなら
そのへんは適当に。

m個のデータから正規方程式を使い
A=Cexp(Bcosφ)を知りたい時どうすればいいですか？
Bは定数で、未知のCとφが知りたいんです。
未知定数が2つで積っていうのがよくわかりません

韓国に攻撃されてるぞｗｗｗｗ７
http://jfk.2ch.net/test/read.cgi/news4vip/1229319433/

世界フィギュア選手権で浅田真央が韓国を破り金賞に輝く
↓
韓国人が腹いせに2ちゃんねるを攻撃
↓
vipper今夜19:00 反撃凸　←いまここ★

目標
http://www.dcinside.com/

Dcinside（ディーシーインサイド）は、韓国の電子掲示板サイトである。
デジタルカメラの情報を提供する専門サイトとして有名になりつつ、だんだん人が集まり韓国最大のコミュニティーサイトとなった。 韓国の2ちゃんねるともいえるサイト。
2004年に2ちゃんねるのサーバをDoS攻撃したことで有名。

ｍ質点の串団子モデルの剛性マトリクスを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

ここはこうしたほうが良いなどなんでも結構ですのでご助言を頂けると幸いです。

for N = 2:row_mなくして
Nの代わりにまんま2:row_mを使え

すいませんがよろしくお願いします。
waveファイルをインポートするにはどうすればいいのでしょうか？
helpをみてもできません。
[y,Fs,Bits] = wavread('onsei1.wav');
と入力したら、
??? エラー ==> wavread at 67
ファイルを開けません
とでます。
誰か教えてください！

plotで図示したグラフとcontourで図示したグラフを重ねて表示する際、
X軸Y軸が逆になってしまいます。
contour(y,x,z)とすれば軸はあうのですが、いまいち納得いきません。
どうしてこうなるのでしょうか？

画像をimresizeでbicubicで補間しながら拡大・縮小したときの動作確認がしたいのですが
どうしたらいいでしょうか？

もしくはbicubicで補間しつつ拡大縮小を少数精度の倍率を与えたときに正確に動作が
保証できるプログラムを組んだことのある方いらっしゃいませんか？
ソースをアップして欲しいのです・・・。修士論文が間に合わない・・。

Y = find([0 1 0 0 1])
Y
　2 5
の逆ってどうやるんでしたっけ？

>>849
ファイルがそこにないんじゃない？

>>850
meshgrid を使ったと仮定して、
meshgrid の使い方間違ったんじゃない？

>>852

X(Y)=1

プログラム初心者です。
MATLABを使えるようになりたいんですが
基本的な使い方とか解説してあるいいサイトとかあったら教えていただけないでしょうか。

サイトではないですが、「MATLAB+Scilabﾌﾟﾛｸﾞﾗﾐﾝｸﾞ事典」という書籍にはお世話になりました。MATLABに似たscilabについても記載されていますが。
後はMATLABのhelpでしょうか。

>>856
ありがとうございます。
探してみようと思います。

Student Editionって以前は自分名義のクレジットカード+HPへの登録+学生証のコピーの郵送が必要って
Q&Aに書いてあったと思うのですが、今は生協で買えるんですねー
その場合は上記の郵送とか面倒な作業は必要ないんでしょうか。

>>858
生協で学生証提示でオケ

>>859
返事遅くなってすいません。
煩雑な手続きは無くなったんですね。ちょっと興味有ったので買ってみることにします。
ありがとうございました。

scilabってコマンドはマトラボとおなじなんでしょうか？
家でもやろうかと思ってるんですが

言い直すと学校でマトラボを使っているんですが
家でもやりたいのでscilabをだいようできるでしょうか？

http://www.ipc.yamanashi.ac.jp/bul/final01/ohki/matlab.html
>Scilab の文法は MATLAB と良く似てはいるが， 細部には違いがある．

>Octave は， 基本的な言語仕様の点で MATLAB 互換となるように設計されており，
>Octave で作成したプログラムの多くは修正なしでそのまま MATLAB で実行可能である
>(MATLABの方が機能が豊富なので，逆は真でない)．

>>863
わかりやすくありがとうございます

freematは?

octaveを語ろう
http://pc11.2ch.net/test/read.cgi/unix/1086016412/

Scilabを触ってみるｽﾚ
http://science6.2ch.net/test/read.cgi/kikai/1068735555/

SCILAB&SCICOSについて教えてください
http://science6.2ch.net/test/read.cgi/sim/980480122/

どれも過疎化しているが…
freematスレは見当たらなかった。

下記を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';

matlabでH.264動画を表示させるにはどうしたらよいのでしょうか。
「mmreader」で読み込んでも走査線が分かれて表示されます。
インターレース動画は奇数偶数をいじって解除・表示できたのですが、
H.264やmp4コーデックでエンコードされた動画は上1/3が黒く表示され、
走査線がバラバラに表示されます。
H.264動画の走査線をどのように分解したらよいのか詳しくわかりません。
詳しい方がいらっしゃいましたらよろしくお願いします。

プログラミング初心者です。画像処理の質問です。

ある画像から2値化を行い，1つの物体だけを取り出そうと考えています。
bwlabelを使ってラベル画像を作った後に，ラベル付けされた各物体の特徴量として面積と円形度を求め，ある範囲に入っている物体1つを取り出したいのですが，よくわかりません。

調べたところbwselectがあったのですが，これは画像をクリックする操作があるため出来れば避けたいと思っています。

わかる方がいましたらアドバイスを頂けないでしょうか？
よろしくお願いします

MATLABを使っていて分からないことがあったので、画像処理について質問させてください。

0〜1の間の乱数を128×128の行列に入れて、これをXと置いた後、
imshow(X,128)で濃淡レベル128のグレースケールで表示します。
これで、白〜黒の画像が表示されるのですが、これを白〜赤のように色をつけて出力したいです。
colormap(Red);imshow(X,128)
でできるかと思ったのですが、上手くいきません。どうすればよいか教えていただけないでしょうか。

>>870
MATLABよく知らんから逆に聞くけど、colormap(Red)するとどう出力されるの?
もしかして、黒-灰-白が黒-赤-白になるとか? colormapの仕様をよく調べてみるといいんでない?

Redがサポートされてないってオチだろ
Hotとかにしてみ

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

スレ汚し失礼いたしました。
どうも上のコードで良かったみたいですので、高速化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

Linux用Matlabについての質問です。
32 bit版のmatlabがあるのですが、64bitカーネルのLinuxで動きますか？

つかってるLinuxによるとしかいえない。

plotでデータ数5000の系列を、同じグラフに5本書きます。
それぞれの系列をカラーでわけたくない場合はマーカで区別すると思うのですが、
データ数が5000だとどのマーカも潰れて見えません。

データ数100個飛ばしでマーカをつけるなど、そのような飛び飛びでマーカをつける方法はありますか？

カラーバーにおいて、特定の値の色を指定する方法を
教えて頂けないでしょうか。
今、カラーバーの範囲を[a, b] (a < 0 < b) にして、
0 の色を白に指定したいのですが、苦労しています。
よろしくお願いします。

http://www.osaka-kyoiku.ac.jp/~ashino/topics/wavehist.html?

このページの真ん中当たりにコッホ曲線のウェーブレット変換が書かれてあるのですが、
この上側の図に書かれてあるコッホ曲線はよく見かける三角形をたくさん合わせたものと全く違いますが、何を表しているのでしょうか？
下の図は3次元的な図になっていますが、それぞれ横、縦、手前の軸は何を表しているのでしょうか？

eval関数について質問です．
変数を文字列化してテキストボックスとしてfigureに挿入しようと
思っているのですが，以下のようにやってもうまくいきません．

eval(['text(20,70,',num2str(k),',''FontSize'',16,''FontName'',''Century'');']);

基本的に([])内の'と,の使い方についてよく理解していません．
'〜',〜,'〜',〜,…　のような順番で結合したい要素を配置して，
''内に更に'〜'のような表現がある場合''〜''に書き換える決まり？が
あるんじゃないかと勝手に解釈して上記のコマンド書きました．

非常にややこしい文になりましたが，そこらへんの事も含めて
ご教授いただきたいです.よろしくお願いします．

>>880
俺の脳内MATLABで実行してみたら、こんな感じだと俺の妄想どおりになった
eval(['text(20,70,''',num2str(k),''',''FontSize'',16,''FontName'',''Century'');']);

>>''内に更に'〜'のような表現がある場合''〜''に書き換える
そうそうそんな感じ
http://dl.cybernet.co.jp/matlab/support/manual/r13/toolbox/matlab/ref/?/matlab/support/manual/r13/toolbox/matlab/ref/eval.shtml

>>881
仰るとおりのコマンドでうまくいきました．
リンク先も一応見た事あるのですが，evalで長いコマンド括る方法が無かった物で…
おそらくevalで一括するよりも，先に文字列作るなりしたほうが良いのかも知れませんね．
ありがとうございました．

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]　のままにしたいのですが、
どうようにすれば実現できるでしょうか？よろしくお願いします。

>>883
plot(1:18, [0.1:0.1:0.9 0.91:0.01:0.99].^2)
でプロットしてからXTickLabelを変更すればいいんじゃない？

そんな姑息な真似はいやなんです。

そうかい、すまんね。

既出かもしれませんが、
たとえば、(z-0.99)/(z-0.58)みたいな伝達関数があって、
任意のzdをｚに代入したい場合、つまり、ｚｄ＝0.45＋0.53ｊをｚにいれるには
どうしたらできるでしょうか？
それとも、そもそもMatLabではできないですか？

伝達関数、と限定した言い方をするってことは
こんな↓単なる台数方程式的な話ではないんですよね。
z=0.45+0.53*j
x=(z-0.99)/(z-0.58)

なんだか分野を書かない質問が多いよね。
信号処理の質問を経済学の人が見てもしょうがないのにさ。
ちなみに filter だったと思う。
signal processing toolbox に一度ひととおり目をとおしておくとよいでしょう。

>>888
方法が見つからなかったので、結局、代数方程式を新たに作って解決したのですが、
もし、tfで作った関数に当てはめることができれば、そういった手間を省けると
思って質問しました。

>>889
すいません。分野は機械工学のデジタル信号処理です。離散時間を扱ってます。
あと、
そういったツールボックスがあるのですね。勉強になりました。

どうもすいません。matlabのコードを.Net向けにコンパイルを行ってその結果を受け取っているのですが
function [a,b] = func( c )
の様に引数を2つ以上返すとき、どうすれば２つとも受け取れるのでしょうか？
というより文法的に何か可能は方法があるのでしょうか？

>>891
何の言語を使っているのか知らんが、その言語スレで聞いたら?
まぁ、構造体を作ってって処だろうけど。

行列の列ごとに、0で無い値の中での平均を求める方法を教えてください。

for i = 1:size(A,2)
mu(i) = mean(A(A(:,i) ~= 0, i));
end
いや、
mu = sum(A, 1) ./ sum(A ~= 0, 1)
こっちのほうがスマートかな

meanの引数に直接行列をぶち込むのと結果が変わらなくね?

mean なら0もカウントしちゃうぜしょ。

figureで軸を指定した関数で表示するにはどうしたらいいでしょうか？
logプロットの様な感じでlog以外の指定した関数で表示したいです。

基本的な0次、1次のBessel関数でMATLAB組み込みより速いのってありますか?

３次元において、ある点(x、y、z)のz座標によって点の色を変えたい(例えば正なら赤系、負なら青系)のですが、何か良い方法はありませんでしょうか？

答えていただけたら幸いです。
n×m行列から、任意の数の列をランダムに抜き出すプログラムはありますか？
宜しくお願いします。

あります。

>>901さん
返事ありがとうございます。
具体的にはどんなプログラムですか？

>>902
「具体的」も何も、一行二行で書ける程度の簡単なものです。

このような場合、礼儀として、まず自分が試したことを説明すべきでは？
Matlab をどの程度習得しているのか。何が出来て何が出来なかったのか。
都合よく答えだけ貰おうとする姿勢は物乞いのそれです。

>>903

4月に入ってからmatlabの勉強をし始め、書いてあるプログラムは
理解できるのですが、自分で組むことができません。
サーバーネットのサイトや本を参考にしようと探したのですが、
見つけることができず2chにきました。
ぜひ教えていただけませんか？

>>904

サーバーネット→サイバーネット

なるほど。
個人的に一番役に立つリソースは本家 mathworks のサイトだと思います。
そちらも参照して頑張ってください。

>>679 にもありますが、今回は randperm を使えばよいでしょう。
「任意の数の」列を抜き出す箇所はご自身で組んでみてください。
dat = [1 2 3; 4 5 6];
result = dat(:, randperm(size(dat, 2)));

>>906

ありがとうございます。
randpermを使って組んでみようと思います。

>>900
宿題投げるならこんな場所もあるぞ
http://www21.atwiki.jp/matlaber/
やってることはこのスレと大して変わらんが

MATLAB環境で動くプログラムに関する質問です。

最近MATLAB2008b student version を購入した素人なのですが困っています。
「MATLAB環境であれば特に設定はいらず、非MATLAB環境であればMCRを
インストールすれば使えます。」
のような説明がされているプログラムを起動しようとしても
"Could't find the mclmcrrt79..."
のようなエラーが発生して起動できませんでした。
しょうがないので適切なMCRをインストールしたら起動できるように
なったのですが、どうしてMATLAB環境であるはずなのに
そのままで起動できないのでしょうか？
student versionだからでしょうか？

MATLAB環境用にしか作られていないプログラムも使いたいので、困っています。
m-fileを開くとMATLABは普通に動き出しているのでMATLAB自体は正常だと
思うんですが・・・

どなたかご教授くださりましたもう。
Win/Vista/Home Premium　SP1

そのプログラムとは何ぞや？
「MATLAB環境であれば特に設定はいらず、非MATLAB環境であればMCRを
インストールすれば使えます。」ってそれらしく説明してても、
実は無条件に mclmcrrt79.dll を読みに行ってる駄目プログラムの可能性もあるよね。

質問させて頂きます。

urlread('http://eow.alc.co.jp/love/UTF-8/')

としてwebの情報を取得しようとしているのですが、以前はできたのに突然できなくなってしまいました。
サーバーの方からアクセスを拒否しているのか、それとも私の方の設定が間違っているのかわかりません。
どなかた上記のコマンドを実行してwebの情報が獲得できるか確認していただけないでしょうか？

>>911
アルクの仕様が変わったんだろ。

>>912
ですかね？
そうするとプログラミング初めから書き直しですね・・・

>>911
>>913
ヒント
そのサイトは HTTP_USER_AGENT を見て、
コンテンツの表示を制限している。

>>913
プログラミングは書き直しじゃないと思うぞ。

逆はさておき、MATLABの話題じゃない希ガス。

>>914
ありがと。ということは
java.lang.System.setProperty('http.agent','agent-name')
でUserAgentを変更すればいいんですね。

再度質問なのですが、UserAgent名を変更した後、どのようにUserAgentをデフォルトに戻すことができるのでしょうか？

java.lang.System.setProperty('http.agent','agent-name')
じゃできんかった・・・

MＡTLABで成分が具体的な数字ではなくa,b,c等の文字からなる行列の
行列式をdetという関数で表示させることは出来ますか？
それとも全て具体的な数値出なければ出力してくれないのでしょうか
｜a b c｜
｜d e f｜　　　　←このような感じの行列です
｜1 2 3｜

matlab知らんけど、試してみればいいんじゃね?

シンボリック関数

Mathematica使えばいいじゃない

>>911
そうか、確かに駄目プログラムという可能性もありますね。
実はある種の映像解析ソフトでありまして、ジャーナルNature Methodsに
載った論文で公開してるやつなんですが、
http://www.nature.com/nmeth/journal/v6/n4/extref/nmeth.1310-S12.zip
といきなりリンクはっていいのか知りませんがどなたかご興味があれば
試していただけたらなぁなんて、、
ここでこういうのはちょっと場違いですかね。。

ただ自分の環境（MATLAB student 2008b on Win Vista)だからダメなのかソフトが悪いのか
わからなくて。

ちなみに
MCRをインストールする前：　"Could not find the mclmcrrt79..." というエラーが出て起動できない
MCRをインストールした後：　一応起動はできるが、動画ファイル(avi or wmv形式)を読み込もうとするといきなり終了してしまう
という状況で困っているんです。

まぁ、製作者に聞けって話ですがねぇ（英語でのコンタクトがちょっと面倒ですが‥）

なら直リンするなよ・・・

>>923
まぁ直リンの何がマズいかも知らなさそうな子だから大目に見てあげなよ。

>>922
このアプリケーション、Matlab 入ってる環境でも実行時にエラー出るね。
ソースを見ても MCR のラッパー呼びまくってるから、
MCR 必須ということで終了。

>>924
直リン以後気を付けます。。

MATLAB2008bが必要とだけ説明があったんですが（それ以上のヴァージョンは不明）
結局はMCRが必要だなんて。。
でもってMCR入れたところで、起動はできても動画読み込めないし。。

ちゃんとしてほしいわぁ。

というかわざわざ試して、
そしてソースまで見てただいてありがとうございました。
優しい方もいるんだなぁとちょっと感動。

直リンの何がまずいのか本気で分からん
ボットがDOSでもしかけんの？

>>920
遅くなりましたが有難うございました

時間変化する３次元の座標データがあり、それの時間変化をプロットしたいのですが、どのような手法を使えばいいのでしょうか？？

>>928
どう見せたい（見たい）か、によるんでないの？

４次元データにxyz座標と時間が入っているので、時間に連動して座標変化の軌跡を表示したいです。
figureで時間変化を追ってプロットしていくのに必要な関数などがわかりません。
ご教授お願います。

どう見せたい（見たい）か、によるんでないの？

3次元ディスプレイなんてものが使えるわけじゃないんだろうから。

>>931
こぴられたｗ

>>930
時間変化しない3次元のデータだったらプロットの仕方はわかる？
それがわかれば、あとは一定時間で区切って再描画コマ送りするなり、
それをアニメーションとして保存するなり、お好きなように

Cソースで書くS-Functionで、複数のデータ型を含むバス信号を入力として
受け取りたいのですが、方法が見つかりません。
もしかして、S-Functionではデータ型が混在するバスやベクトル信号は
受け付けられないのでしょうか?
もし受け付ける方法があるなら教えていただけないでしょうか?
よろしくお願いします。

>>933
ポートをわけて渡さないとダメじゃないかなあ

>>934
やはりそうですか・・・
データ型を取得するマクロがポート単位のものしかないですもんね・・・
複数データ型のバスを受け取らなくても良い方法を考えることにします。
ありがとうございました。

>>932
その都度figureを立ち上げてコマ送りのように表示することはできたのですが、もう少しスマートにできないものでしょうか？？
一定時間で区切りコマ送りする手法を教えていただきたいです。

>>936
少しは情報を探したのか？
http://www.cybernet.co.jp/matlab/support/techkwdb/search_view.php?id=1303

差分方程式
x(n+2)+a_1x(n+1)+a_0x(n)=0 ,
x(0)=x_0 , x(1)=x_1 , n=0,1,2,∙∙∙
の階x(k)=0,1,2,∙∙∙,t_fを求めるプログラミングを作成せよ。また、次のデータについて実行せよ。
a_0=0.3 . a_1=-1.5 , x(0)=1 , x(1)=2 , t_f=10

のプログラムが作れません。誰か教えてください！お願いします！！

質問です

みなさんマトラボはいつ頃から勉強はじめました？

ぜひ教えてください

>>939
初めて触れたのは学部４年前期の制御工学の授業

>>939
大学2年から

>>938
久々のレスが宿題かよ！
x(n+1)+a_1x(n)=0, X(0)=X_0
なら作れる？
>>939
おれも940と同じ

>>939
院にいってから

ちょっと質問させてください。

assignin関数でfunction 内の変数をWorkSpaceに置きたいんですけど、
その際に変数名を配列で指定するとエラーがでて出来ません。
何か良い方法はないですかね？

こんな感じで書いてます。

function hoge
str =[ test1 test2 ]
for i=1:2
　　assignin('base',str(1,i), 2); ←ここでエラー
end

よろしくお願いします。

三つアドバイス。
1. エラーメッセージをしっかり読む
2. ブレークポイントを設定して、関数実行時の変数の中を見てみる
3. 文字配列とセル配列の違いと扱い方を体で覚える

function hoge

str = { 'test1' 'test2' };
for k = 1:length(str)
assignin('base', str{ks}, 1);
end

あと、matlab では i は虚数を意味する（i * i = -1）ので、
できるだけカウンタとしての使用は避けるべきなのだそうな。
標準的な関数群でも平気で i がカウンタとして使われているけれども。

ごめん typo があった。

function hoge

str = { 'test1' 'test2' };
for k = 1:length(str)
assignin('base', str{k}, 1);
end

>>945-946
素早く丁寧な返事をありがとうございます。

１と２に関しては実行しました。
エラー文は読んだのですが、不適切な文ですとしか出ず
str(1,i)を　'test1' におきかえると正しく動作したので、str(1.i)が
間違っているということはわかったのですが、そこから進まなくて・・・。
セル配列{}? と文字配列()の違いでエラーが出てたという事ですかね？
str(1,1)とコマンドラインに入れればtest1　が帰ってきたので、
eval(str(1,1))
などを入れてやれば解決するかなと考えたのですが、これでもうまくいかず・・・。

また、明日学校で試させてもらいます。ありがとうございました。

この例題なら自宅でOctave使って試すこともできるぜ

>>948
お勧めされたので早速やってみました。
実は質問外なので省いたのですが、最初の
str = { 'test1' 'test2' };
を取ってくる際にimportdataを用いてエクセルファイルのヘッダを取得しているのですが、
その辺りの関数がOctaveにはないようですね。

ただ、上記の関数がOctave上で正しく動くことは確認できました。ありがとうございます。
明日学校でMATLAB上でも正しく動くか検証してみます。

すごく初歩的な質問ですまんが


%f(x)=x^2-5
N0=100;
TOL=0.0000001;
XL=1;
XR=5;
i=1;
FXL=XL^2-5;

while (i<N0)
XM=(XL+XR)/2;
FXM=XM^2-5;
if (XM==0 | (XR-XL)/2<TOL)
disp(XM);
return;
end
i=i+1;
if (FXL*FXM>0)
XL=XM;
FXL=FXM;
else
XR=XM;
end
end
disp('解が得られない');

このプログラムを関数x^3-7x^2+14x-6で動かしたいんだけど、どこをどう書き換えればいいのかわからん
自分でやるといつもFXM=のところでエラーになる
教授から丸投げされて、もうお手上げ状態なんで誰かアドバイスください

>>950

学校のＰＣからじゃ書き込めないから携帯からになるけど

７行目を
FXL = XL^3-7*XL^2+14*XL-6

１１行目を
FXM = XM^3-7*XM^2+14*XM-6

にすればいいのでは？

>>945,946
遅くなりましたがＭＡＴＬＡＢでも動作確認できました。
丁寧なアドバイスありがとうございます。

俺はmaximaに移行した

質問です。
データにフィルターをかけるとデータの最初の方が０に近い値になってしまうのですが、回避する方法はあるのでしょうか？
それとも、フィルターをかけた上で最初の方を切り取った方がいいのでしょうか。

使っているフィルターは３次のバターワースフィルタです。

すみません、sageわすれましたorz

位相遅れを無くしたいならfiltfiltを使えばいい

>最初の方が０に近い値になってしまう
これが何か問題になるの？

2次のルンゲクッタ法のホイン法と中点法でプログラムを作っています。
xx = [];
i=0;
for t=0:dt:tf
i=i+1;
xx(:,i)=x;
xt = x;
for j=1:2
f = A*x + B*u;
d(:,j) = f*dt;
x = xt + d(:,j)*0.5;
end
x = xt + (d(:,1) + d(:,2))/2;
end
t=0:dt:tf;

どこか違うとこありますか？あったら教えてください。

2次のルンゲクッタ法のホイン法と中点法でプログラムを作っています。
ホイン法
xx = [];
i=0;
for t=0:dt:tf
i=i+1;
xx(:,i)=x;
xt = x;
for j=1:2
f = A*x + B*u;
d(:,j) = f*dt;
x = xt + d(:,j);
end
x = xt + d(:,2);
end
t=0:dt:tf;
中点法
どこか違うとこありますか？あったら教えてください。

2次のルンゲクッタ法のホイン法と中点法でプログラムを作っています。
中点法
xx = [];
i=0;
for t=0:dt:tf
i=i+1;
xx(:,i)=x;
xt = x;
for j=1:2
f = A*x + B*u;
d(:,j) = f*dt;
x = xt + d(:,j)*0.5;
end
x = xt + d(:,2);
end
t=0:dt:tf;

どこか違うとこありますか？あったら教えてください。

2次のルンゲクッタ法のホイン法と中点法でプログラムを作っています。
中点法
xx = [];
i=0;
for t=0:dt:tf
i=i+1;
xx(:,i)=x;
xt = x;
for j=1:2
f = A*x + B*u;
d(:,j) = f*dt;
x = xt + d(:,j)*0.5;
end
x = xt + d(:,2);
end
t=0:dt:tf;

どこか違うとこありますか？あったら教えてください。

なんだこの馬鹿は

Student Editionてamazon.co.jpでも買える？

matlabにはfortrunでいうgoto文のようなものはありますでしょうか？
○○行から××行の間の処理を繰り返したいのですが、
結構な行数なので見辛くなってしまいまして。。
ご存知の方がいましたら教えて頂けると幸いです。

結論：無し
ついでに書くならmatlabはループが苦手

>>964
やっぱりないですよね。。
ご回答有難うございました。

uicotrol関数のプッシュボタンやトグルボタンを
他の形(丸や多角形)にできますかね？
プロパティをみるとpositionによる縦横長さを
弄る以外にはなさそうですが、何か手段はありますでしょうか？

デジタル信号の加速度時系列波形を
２階積分して変位の時系列波形を求めたいのですが、
どのように組めばいいのかがわかりません。
サンプルプログラムが載っていると
順を追えて理解がしやすいのですが、
何か参考となる本や文献があれば教えて頂きたいです。

FFTの勉強中でプログラムを実際に作ってみたんですが
データ数８までは成功するのにどうしてか１６以上になるとできません。
XXには既に並び替えたデータが入っています。
処理された結果がYです。

phase=log2(N);
for m=1:phase
W=exp(-i*2*pi/(2^m));
for n=1:N/2^m
for k=1:2^(m-1)
K=k+2^m*(n-1);
Y(K)=XX(K)+W^(k-1)*XX(K+2^(m-1));
Y(K+2^(m-1))=XX(K)-W^(k-1)*XX(K+2^(m-1));
end
end
XX=Y;
end

質問させてください。
matlabでImageJと同等のリニアスケールのトゥルーカラー画像を表示させたいのですが、うまくいきません。
matlabでは明度がはっきりしない画像となってしまいます。

何か良い方法があれば教えてください。

Matlabで大きな画像（1000x4000）をimagescで表示させ、saveasで保存すると、
元の画像よりも小さくなってしまいます。。

元のピクセルサイズを壊さずに保存できる方法を、
もしご存知の方がいたら教えてくださいｍ（_　_）ｍ

自己解決しました；；

imwrite(X, 'test.png');
のようにして解決しました。

拘束条件（壁）があるときの2自由度ロボットアームのシミュレーションをしたいのですが、
拘束条件をどのように設定したらいいのか分かりません；

漠然とした質問で申し訳ないのですが、教えていただけたら幸いです。

エスパーするとode+eventでなんとかなりそう

画像をR成分、G成分、B成分と段階的に表示するプログラムを作りたいです。

A=imread('画像のアドレス');
B=zeros(size(A,1),size(A,2),3);

for(a=1:1:3)
for(i=1:1:size(A,1))
for(j=1:1:size(A,2))
B(i,j,a)=A(i,j,a);
B(i,j,a)=B(i,j,a)/255;
end
end
imshow(B)
end

このプログラムだとimshowが三回実行されますが、
実行間隔が短いので最後の1回の結果しか確認できません。
それぞれの間に10秒ずつ時間を置きたい場合などはどうしたらよいですか？

自己解決しました。

pause(n)
でn秒プログラムの実行を停止できました。

画像を読み込み、読み込んだ画像情報をファイルとして書き出す際に、
連番をファイル名に付け様と思ったのですが、上手くいきません。

write = imread('lex.jpg');
for i = 1 : 10
imwrite(write, ['file_' i '.jpg'], 'JPEG');
end

使用したのは上記のプログラムですが、imwriteのところでjpegを書き出すためのファイルが開けない、
というエラーがでます。
これを改善するにはどうすればいいのでしょうか。
よろしければ誰かご教授のほどをお願いいたします。

int2str

int2str

質問です。
　A = diag([1,1,1,0,0,0]);
　eigs(A, 2, 'SM')
とすると
警告: Matrix is close to singular or badly scaled.
Results may be inaccurate. RCOND = 0.000000e+000.
が出て、結果が NaN, NaN となります。

どう対処すればよいでしょうか？
なお、実際に使いたい対象はもっとでかい素行列で、
full したくない、という事情があります。

singularでない行列を放り込むしかないんじゃね？

>>980
実際、今は入力の行列に小さな摂動を加えて
singularではないようにして回避していますが、
どうにも筋が悪いように思い、質問した次第です。