MATLABプログラミング 質問箱 その2
■ このスレッドは過去ログ倉庫に格納されています
0001デフォルトの名無しさん
2006/11/10(金) 11:28:22http://www.cybernet.co.jp/matlab/
MATLAB 本家
http://www.mathworks.com/
Yahoo MATLAB
http://dir.yahoo.co.jp/Science/Mathematics/Software/MATLAB/
MATLAB入門
http://feedback.mech.fukui-u.ac.jp/matlab/matlab0.htm
MATLAB チュートリアル
http://www-isys.eei.metro-u.ac.jp/~sho/pub/ttmatlab/Index.html
0685デフォルトの名無しさん
2008/08/21(木) 10:07:19******
for文を使わないことはできますか？
２次関数(ax^2+bx+c)を入力にして、総あたりでa,b,cのパラメータを推定する、というようなことをやっています。
入力をf(x)として、-5<x<5の範囲を見てるのですが、コードは
error=inf;
for a=-10:1:10
for b=-10:1:10
for c=-10:1:10
for t=-5:0.1:5
e=(f(t)-(a*t*t+b*t+c))^2;
if (error>e);
error=e;
end
end
end
end
end
というような感じになってます。これをfor文を使わないでコロンを使った形式にできるものなのでしょうか？
よろしくお願いします。
0686デフォルトの名無しさん
2008/08/21(木) 10:10:21よろしくお願いします。
0687デフォルトの名無しさん
2008/08/21(木) 10:16:19で、あとは min で argmin うけとる。
0688684
2008/08/21(木) 16:35:55ご指摘いただいた方法でできました。
ありがとうございます。
0689デフォルトの名無しさん
2008/08/21(木) 21:09:03だったらrepmatじゃだめじゃね？
ってか俺も今
for a=1:10
for b=1:5
kotae(a,b)=a+b
end
end
って感じのやってるけど、これfor文じゃなくせるの？
0690デフォルトの名無しさん
2008/08/22(金) 00:13:06「500の最大の再帰呼出の制限達しました。制限を変更するには、・・・」
となって動きませんでした。
制限を2000あたりまで変更しても同じ結果で。3000くらいになるとMatlabが
強制終了してしまいました。
これ何か対処法ないでしょうか。。？？
0691デフォルトの名無しさん
2008/08/22(金) 02:39:48つか、少しは自分でk
0692デフォルトの名無しさん
2008/08/22(金) 02:56:10なくせう。for 文使うとゲロのろいぞ matlab。
0693デフォルトの名無しさん
2008/08/22(金) 09:31:000694デフォルトの名無しさん
2008/08/22(金) 10:26:09kotaeって10x5のマトリックスになるだろ？
0695684
2008/08/22(金) 17:33:09meshgridというのを使うとてっとり早いです。
行列にして掛け算する、ということです。
687さんがおっしゃった通り、repmatを使ってもできます。
というか、repmatを使ってmeshgridが実装されてるようなイメージです。
0696MATLAB初心者
2008/08/23(土) 12:06:04ＭＡＴＬＡＢをはじめたばかりの初心者なのですが、連番をつけるところで悩んでおります。
文字列の配列で0001、0002、..... 、9999のように整数で任意の桁数にしたいのですが
うまい方法はありませんでしょうか？
エクセルでやっていたときは、上記の場合であれば”000”という文字列と整数の値を結合し
左から4文字を切り出して使用しておりました。
ＭＡＴＬＡＢでもエクセルでやっていたような力技でやるしかないのでしょうか？
もしご存知の方がおりましたら、ご教示いただけませんでしょうか？
よろしくお願いいたします。
0697MATLAB初心者
2008/08/23(土) 12:12:11上記のレスで左から切り出していたというのは誤りで、右から切り出しておりました。
0698デフォルトの名無しさん
2008/08/23(土) 14:51:58ttp://dl.cybernet.co.jp/matlab/support/manual/r13/toolbox/matlab/ref/?/matlab/support/manual/r13/toolbox/matlab/ref/sprintf.shtml
0699MATLAB初心者
2008/08/23(土) 15:46:04ありがとうございます。
大変申し訳ございません。私には、このページではsprintfでたとえば”1”を”00001”と定義する方法が理解できませんでした。
整数を5桁で表して、それを文字型にnum2strか何かで文字列に変換するようなことができるということでしょうか？
私の思い付きではstrcatをつかって696に書いたような力技で書く以外に方法がわかりません。
お手数をおかけしますが、ご教示いただいたページのどれを使用すればよいか教えていただけませんでしょうか？
0700デフォルトの名無しさん
2008/08/23(土) 16:25:44printf フォーマット
でぐぐればいろいろ解説があるよ。
Ｃ言語が元だけど自由度がかなり高いから覚えとくとよい
↓とりあえず解決策
str=sprintf('%05d',i)
0701MATLAB初心者
2008/08/23(土) 16:35:09取り合えずというか、そのものずばり私の欲しかった回答を書いていただいて、感謝の限りです。
お忙しい中、どうもありがとうございました。
これからも精進いたします。
0702675
2008/08/26(火) 03:23:35みたいなベクトルがあって、[1 2] のある位置を知るのに、
for ループを使わずにすむ方法ってありますか？
matlab ならありそうな気がするのですが
B = [0 1 0 0 0 0 0 0 1 0 0]
が返ってほしい。
0703デフォルトの名無しさん
2008/08/27(水) 10:05:550704デフォルトの名無しさん
2008/08/28(木) 08:09:50x=find(A==1);
y=find(A==2)-1;
z=intersect(x,y);
B=zeros(1,length(A));
B(z)=1;
もっと簡単なのきぼん
0705デフォルトの名無しさん
2008/08/28(木) 22:26:43((A(1:length(A)-1)+A(2:length(A))).*(A(1:length(A)-1)==1))==3
正の値で且つ後ろ一つは無視する前提
一般性があるかどうかは知らない
0706デフォルトの名無しさん
2008/08/28(木) 22:38:400707デフォルトの名無しさん
2008/08/29(金) 01:15:31レベル低いなあ。
こんなのすらすっと出てこないってことは
よっぽど効率悪いコード書いてるんだろうなお前ら。
0708デフォルトの名無しさん
2008/08/29(金) 01:40:140709デフォルトの名無しさん
2008/08/29(金) 01:46:06スマン
0710デフォルトの名無しさん
2008/08/29(金) 09:07:52勉強になるんだけどな。
0711デフォルトの名無しさん
2008/09/04(木) 13:25:51MATLAB上で画像処理をする際に
画像のある画素(座標)を保存して
さらにその保存した座標のｘ方向成分だけ読み出すには
どういったコマンドを入力すればいいんですか？
0712デフォルトの名無しさん
2008/09/04(木) 13:42:08まず入門サイトでも読んだ方がはやくね？ってかはやい
0713デフォルトの名無しさん
2008/09/04(木) 16:05:56それについて詳しく解説してある本を知っていたら
誰か紹介していただけないでしょうか。
よろしくお願いします
0714デフォルトの名無しさん
2008/09/08(月) 12:56:090715デフォルトの名無しさん
2008/09/08(月) 18:24:34例えば、
Data1 = 1;
Data2 = 2;
・
・
・
Data30 = 30;
と30個の変数を定義したいとき、for文を使って定義することは可能でしょうか？
0716デフォルトの名無しさん
2008/09/08(月) 23:43:57キーワードを英語でlookforするくらいはしようよ。
>>715
help eval
0717デフォルトの名無しさん
2008/09/09(火) 11:36:43matlabのエディットウィンドウで、文字の挿入ってどうやるんですか？
現状では入力してある行に文字を挿入すると後ろの文字が消えてしまい、
全文入力し直しになってしまうのです、これをwordなどみたいに
文字を挿入したら、文が自動で後ろに移動してくれるようにできないですか？
また、同様にスペースキーも文字間のスペースをあけてくれず
文字を消してしまうのですが？どうしたらいいですか？
OSは vinelinux4
matalab v7.0.1 SP1です
0718デフォルトの名無しさん
2008/09/09(火) 22:12:45キーボードのinsert
0719デフォルトの名無しさん
2008/09/10(水) 13:58:00ボタンを押すとシミュレーションが実行するようにしたんですが、
GUIなしでコマンドから直接実行した場合だと１６０秒くらいで終わるのに、
なぜかGUIからだと２６０秒もかかるんですが、
なんでこんなにかかると思いますか？
やってることはまったく変わらないですが・・・。
全然見当がつかん。
プロファイラを実行した状態で GUI 経由でシミュレーションを行って、
時間のかかっている箇所を割り出してから相談してもらった方がよいと思われ。
0721デフォルトの名無しさん
2008/09/11(木) 00:44:20アドバイスありがとう。
実はこういう経験あまりしてなかったんで、
プロファイラとか使ったことなかったんです。
ちょっと調べて割り出してみます。
0722初心者
2008/09/11(木) 01:43:15Cでいうポインタのような概念はMatlabにはないのでしょうか？
0723デフォルトの名無しさん
2008/09/11(木) 01:54:01ないはずだよ
wikipediaにも書いてある
まあ、google codeとかmatlab centralで検索すれば
自作のポインタ関数出て来ても可笑しくないけど
0724デフォルトの名無しさん
2008/09/11(木) 02:15:40デフォはポインタ
渡された変数を編集しようとした時点でコピーが始まる
0725デフォルトの名無しさん
2008/09/11(木) 06:15:09何をしたいのかをはっきりさせて
それをMATLABではどうやるのか、と考える方がいいよ。
他の言語のやり方をそのまま持って来ようとしても
そうは行かない場合は多い。MATLABに限らず。
0726デフォルトの名無しさん
2008/09/14(日) 02:33:41例えば
3,6,4,5,2
4,15,7,5,3
6,4,4,7,3
のような行列があってノイズで15だけ値が大きい場合には、
それを周りの数字に近い値に矯正（たしか平滑？）するにはどうすればいいですか？
自分が知ってるやり方だと、ノイズも含め全ての成分が平均化されてしまいます。
どうにかノイズ以外はそのままの値で、ノイズ部分だけを周りの値に馴染ませたいんです。
0727デフォルトの名無しさん
2008/09/14(日) 02:59:5015を適当に5とかにしたら「馴染んだ」ことになるの?
自分が何をやりたいかをまずは具体的にする。
君以外にそれがわかる人はここにはいないよ。
0728デフォルトの名無しさん
2008/09/14(日) 07:57:320729デフォルトの名無しさん
2008/09/14(日) 12:03:36if ノイズ
平滑
else
なにもしない
0730デフォルトの名無しさん
2008/09/15(月) 07:18:580731デフォルトの名無しさん
2008/09/15(月) 13:38:27新しいやり方を考えればいいだけ。
0732デフォルトの名無しさん
2008/09/17(水) 23:20:14Ｃのプログラム等に記載するための離散時間系のデジタル１次ローパスフィルタを作るためにＭＡＴＬＡＢでフィルタの係数を求めようとしています。
そのためのプログラムとして、サンプリングインターバル＝６ｍｓ、フィルタ周波数＝１０Ｈｚとして、
omega = 2 * pi * 10
T = 1 / omega
sys = tf(1, [T, 1])
[num, den] = tfdata(sys)
[numd, dend] = c2dm(num, den, 0.006)
として求めると、y[n] = 0.6859 * y[n-1] + (1 - 0.6859) * x[n]という式が求められます。
一方で、連続時間系の1次ローパスフィルタの伝達関数（サンプリングインターバル＝6ms、フィルタ周波数＝10Hz） G(s) = 1 / (1+Ts) を手計算でz変換を行い、離散系のローパスフィルタの式を求めると、（T = 1 / omega、omega = 2 * pi * 10）
y[n] = 1 / (1+omega*0.006) * y[n-1] + (1 - 1 / (1+omega*0.006)) * x[n]
となるため、値を代入すると、
y[n] = 0.7263 * y[n-1] + (1 – 0.7263) * x[n]
となり、先ほどＭＡＴＬＡＢ関数から求めた式と乖離が出てきてしまいました。
そこで、有識者の方々にお聞きしたいのは
１．そもそも上記のデジタル１次ローパスフィルタの作り方は間違っているのか？
２．c2dm コマンドの中身がよくわからないが、これは計算誤差が出るようなものなのか？
３．離散系に変換したシステムのボード線図を見て特性を確認したいのだが、その場合
も連続系と同じく
bode(numd, dend)
でよいのか？
長文となってしまい、見にくくて申し訳ありませんがどなたか有識者の方がいらっしゃればご教授ください。
0733デフォルトの名無しさん
2008/09/18(木) 14:54:340734デフォルトの名無しさん
2008/09/18(木) 22:11:28と言いますと・・・？
0735デフォルトの名無しさん
2008/09/18(木) 22:13:54間違いを探したいのなら、
他の数字の組み合わせにしてみて、計算してみるのはどうでしょう？
0736デフォルトの名無しさん
2008/09/18(木) 23:19:11すいません、だいぶ他の数字でもやってみたのですが・・・
基本的な結果としては、低周波（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
どなたかご教授いただけると助かります
0737デフォルトの名無しさん
2008/09/18(木) 23:34:44サンプリングインターバルも変えてみたらどう？
0738デフォルトの名無しさん
2008/09/19(金) 00:13:17MATLAB計算値 手計算値
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)
という見方でよろしいのでしょうか？明らかに異なった波形が出てくるので、離散系の場合はボード線図を違う見方で見なければいけないのかなと予想しているのですが。。。
見方がわかりませんので、もしご存知でしたら教えてください。
0739デフォルトの名無しさん
2008/09/19(金) 00:33:01ttp://www.ikko.k.hosei.ac.jp/joomla1506/matlab-topmenu-33/17-matlab20070307.html
0740デフォルトの名無しさん
2008/09/23(火) 00:35:38stem3でつくったグラフを等高線のグラフにそのまま変えたいのですが、
どなたかわかる方いらっしゃいましたらお願いします。
0741デフォルトの名無しさん
2008/09/23(火) 02:50:240742デフォルトの名無しさん
2008/09/23(火) 11:49:33値がどのように分布しているか等高線の方がわかりやすいと思い
変えたいのですが。
0743デフォルトの名無しさん
2008/09/23(火) 15:11:34matlab 等高線
でぐぐれ
0744デフォルトの名無しさん
2008/09/30(火) 15:46:00ファイルguimain.m内で
handles.a=1;
guidata(hObject, handles);
とhandles.aを設定したとき
別のファイルgui2.m内で
handles.aを編集する方法はあるのでしょうか
function gui2(hObject,handles)
handles.a=2;
guidata(hObject, handles);
ではできませんよね?
0745デフォルトの名無しさん
2008/09/30(火) 16:32:46function out = gui2(handles)
handles.a=2;
out = handles;
として、呼び出し側でguidata(hObject, handles); することで解決しました。
0746t
2008/10/07(火) 19:46:19困っています。
0747デフォルトの名無しさん
2008/10/15(水) 20:23:010748デフォルトの名無しさん
2008/10/18(土) 10:33:41構造体使う必要がないならglobal宣言すれば簡単かと
0749デフォルトの名無しさん
2008/10/19(日) 22:46:09どうやったらよいのでしょうか？
試しに
----------------------
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)
----------------------
としたのですが、うまくいきませんでした。
おそらく行列のプロットを理解していない
と思うのですが・・・。
非常に低レベルな質問で申し訳ありませんが、
よろしくお願いします。
0750デフォルトの名無しさん
2008/10/19(日) 23:52:000751749
2008/10/20(月) 00:01:45早速のご回答、ありがとうございます！
無事、3D plot をすることができました。
ありがとうございます！
0752デフォルトの名無しさん
2008/10/23(木) 17:54:34どうやって、変形して最終の形にするのか
見る方法ありませすか？
0753デフォルトの名無しさん
2008/10/23(木) 23:45:04mathematica？
0754デフォルトの名無しさん
2008/10/23(木) 23:47:52シンボリックマスツールズの
0755いかわ
2008/10/24(金) 15:35:48もしも、''がなく、"2008-09-23 02:30:00.1" というフォーマットならば、sscanf(%d-%d-%d %d:%d:%d); によって時間データを抽出できるのですが。。
0756デフォルトの名無しさん
2008/10/24(金) 20:06:56どうすればいいかわからなくなるのかなこういう人は・・。
0757デフォルトの名無しさん
2008/10/25(土) 04:01:54C で使ったことある人いない？
ttp://www.mathworks.com/access/helpdesk/help/toolbox/compiler/index.html?/access/helpdesk/help/toolbox/compiler/f2-972343.html
これの通りやってみて、C をコンパイルできるようにまでは持って
いったんだけど、実行すると
libmatrixInitialization: unhandled exception
と出る。なにが足りないんでしょうかね？
0758デフォルトの名無しさん
2008/10/25(土) 05:49:280759デフォルトの名無しさん
2008/10/25(土) 10:11:57回答にはなってないかもしれんけど、
おいらは、加えてBuilder for .NETを使ってMファイルを.NETコンポーネント化
してから、.NET系言語から、そのコンポーネント呼び出して使ってるよ。
0760デフォルトの名無しさん
2008/10/27(月) 15:42:45だけでメモリ不足になります
余分な変数消去や仮想メモリ等をいじってみましたが効果はありませんでした
これを解消するには結局メモリ増設しかないのでしょうか？
0761デフォルトの名無しさん
2008/10/27(月) 16:38:49メモリ増設しても無理じゃね?
100000*100000って、数値だけで10Gになるからアクセスすることさえ厄介になってくると思う。
メモリ量的には更に数倍になるわけで、オンメモリで処理するのはどうせ無理だし。
0762760
2008/10/27(月) 18:09:13ご回答ありがとうございます
低レベルな質問ですみませんでした
このような行列はMATLABではそもそも扱えない（普通は扱わない）と考えてよいのでしょうか？
株価のシミュレーションでforを使うより行列を使った方が早いため
あのような行列が処理できると助かったのですが
0763デフォルトの名無しさん
2008/10/28(火) 01:30:51１０Ｇのメモリっていうといわゆるスーパーコンピューターの領域
0764デフォルトの名無しさん
2008/10/30(木) 04:26:55処理の進行状況を表示したいのですが、
改行を抑制する方法はありませんか？
0765デフォルトの名無しさん
2008/10/30(木) 21:20:25複素数の実部虚部それぞれ七桁で表示させたいのですが
zに複素数が入っていて
sprintf('%7f',z)
とすると実部だけが七桁で表示されてしまいます
どうしたらよいでしょうか？
0766デフォルトの名無しさん
2008/10/30(木) 22:41:170767デフォルトの名無しさん
2008/10/30(木) 23:07:490768デフォルトの名無しさん
2008/10/31(金) 00:06:48なんで今どうやってるのか書かないの?
0769デフォルトの名無しさん
2008/10/31(金) 14:35:01ありがとうございます
sprintf('%7f %+7f',real(z) imag(z)i)
のような感じでしょうか？
今試せる環境ではないのでiとかがどうなるのかよく分かりませんが・・・・・
0770デフォルトの名無しさん
2008/10/31(金) 19:10:46１辺の長さが１の正三角形を上下に並べて出来る平行四辺形をプロットして、
それに含まれるかどうかの判定をしたいのですが、判定条件はどうしたらいいでしょうか。
N個の点を行列で与えて、図形に入っていれば赤、入っていなければ黒という風に
表示させたいのです。
よろしくお願いします
0771デフォルトの名無しさん
2008/10/31(金) 19:25:24+の位置が不適切だけどそれを直せばイケル
0772デフォルトの名無しさん
2008/11/01(土) 17:39:35すいませんまだ試せていないんですがこうすると実数値の場合に
3+0i
みたいな表示になっちゃいませんかね？
もっと簡単にできる方法ないですかね
0773デフォルトの名無しさん
2008/11/07(金) 14:26:53例えばデータシートが1000枚あって、各データシートからX,Y,Zの値を取り出して、
座標(X,Y,Z)をplot3でプロットします。これがこのシートの3次元での点（位置）になりますよね。
それを1000枚分やると3次元空間に点が1000点でき、動きの軌道のように繋がって表示されます（そういうデータなので）。
問題は次なのですが、数十点ほど動きの軌道からはみ出してプロットされた点（多分計測時のノイズ）が
存在します。なのでプロット点全体を平滑化したいのですが、三次元プロットした点をフィルタするのは
可能なんでしょうか。
よろしくお願いします
0774デフォルトの名無しさん
2008/11/07(金) 21:04:04成分ごとにやれば？
0775デフォルトの名無しさん
2008/11/08(土) 07:04:44フィルタしてからプロットするんでしょ?
0776デフォルトの名無しさん
2008/11/08(土) 09:03:51mlfMatrixPrint とか
mlfLength とかってどこに行ったの？
昔は matlab.h ってのがあったみたいなんだが、
R2007b だけど、ないんだよね。
matlab 内のビルトイン関数を C から直接呼び出せなくなったの？
0777デフォルトの名無しさん
2008/11/09(日) 14:40:02>>675の解決策ってありますか?
0778デフォルトの名無しさん
2008/11/10(月) 07:24:170779デフォルトの名無しさん
2008/11/10(月) 20:04:54250*350でフレーム数200のAVIの動画をRGBで格納するために
T=zeros(250,350,3,200)
としてみたのですが、メモリが足りませんと出てしまいました。
解決策がわかりません
ちなみにPCはXP Core(TM)2 Duo CPU メモリ2Gを使っています。
よろしくお願いします
0780デフォルトの名無しさん
2008/11/10(月) 22:36:17使わずに 作りたいのですが可能でしょうか。
基礎的な質問でしょうがよろしくお願いします。
0781デフォルトの名無しさん
2008/11/10(月) 23:00:26[A;A;B;B]だったらn*2行m列だと思うんですが。
0782780
2008/11/10(月) 23:08:49[A;A;A;B;B;B;C;C;C]という行列を作りたいのです。
0783デフォルトの名無しさん
2008/11/10(月) 23:58:44メモリが足りないのものはしょうがないのである。
本当に全部格納しなければいけないのか考えてみるのである
0784デフォルトの名無しさん
2008/11/11(火) 00:01:09つまり *3 でいいんだね？
X = [A;B;C] からA, B, C の分解は終わっているものとする。
Y = [repmat(A, 3, 1); repmat(B, 3, 1); repmat(C, 3, 1)]
■ このスレッドは過去ログ倉庫に格納されています