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

**デフォルトの名無しさん** · 2013/04/02(火) 22:18:12.30

>>117
epsで吐き出すなら、吐き出した先のソフトで色がついてる背景を消去したら
いいだけじゃない

**デフォルトの名無しさん** · 2013/04/15(月) 02:11:41.52

[eigv eign]=eig(Hamiltonian);

などのように行列を対角化したときに
maxやminのように最大か最小を見つけるのではなく
n番目に小さいeignの要素と、その場所を得る関数はありますか？

**デフォルトの名無しさん** · 2013/04/15(月) 09:57:26.38

[Y,I] = sort(X); つかったら？

**デフォルトの名無しさん** · 2013/04/17(水) 15:29:16.26

>>119
texなんですが、どうやってやれば良いのでしょうか？

**デフォルトの名無しさん** · 2013/04/24(水) 14:45:13.87

バンドパスフィルタにはゼロ位相のものと因果律を
満たすものとがあると思うのですが、ネットで落ちていたんですが
以下のはどちらにあたるか教えて下さい。
よろしくお願いします。

**デフォルトの名無しさん** · 2013/04/24(水) 14:49:17.62

function [bandy]=bandp_fft(y,fs,F1,F2)
dt=1/fs;
[k,l]=size(y);
if k==1
n=l;
else
n=k;
end;
ffty=fft(y);
%------------------------------------ Low-pass(f=<F1)
index1=round(F1* (dt*n)+0.5);
index2=n-(index1-2);
Lowpass=ones(size(ffty));
Lowpass(index1+1:index2-1)=zeros(size(Lowpass(index1+1:index2-1)));
lowy=real(ifft(Lowpass.*ffty));
%figure1=figure(1);plot(lowy);
%------------------------------------ High-pass(f>F2)
index1=round(F2*(dt*n)+0.5);
index2=n-(index1-2);
Highpass=zeros(size(ffty));
Highpass(index1+1:index2-1)=ones(size(Highpass(index1+1:index2-1)));
highy=real(ifft(Highpass.*ffty));
%figure2=figure(2);plot(highy);
%------------------------------------ Band-pass (F2<f=<F1)
bandy=real(ifft(Lowpass.*Highpass.*ffty));

**デフォルトの名無しさん** · 2013/04/30(火) 16:36:38.36

MATLABでAM信号を書きたいのですがエラーが出てしまいます。ソースコードは

fc = 5e4;
ac = 1;
wc = 2 * pi * fc;
fm = 5e3;
am = 1;
wm = 2 * pi * fm;
tt = 2e-3;
n = 10000;
ts = tt/n;
t = 0 : ts : tt;
s = ac * (1 + am * sin(wm * t)) * sin(wc * t);
plot(t(1:n/2), s(1:n/2));

時間波形0～2msを10000サンプリングのうちの時間波形0～1msを表示したいです。
情報信号sin(ωc*t),ωc=2*π*5000
搬送波sin(ωm*t),ωm =2*π*50000

内部行列の次元は一致しなければなりません。とエラーが出ます。どなたか知恵を貸してください。

**デフォルトの名無しさん** · 2013/04/30(火) 17:11:39.09

最後のmとc逆でした

**デフォルトの名無しさん** · 2013/04/30(火) 17:49:07.66

ちんぽ噛んで死ね

**デフォルトの名無しさん** · 2013/05/02(木) 00:12:26.05

rcosineとhanningって一緒なの？
説明見ると、よくわからんが
式が一緒っぽいみたいなのと
plotした波形が一緒なんだが
これを説明できる文章が思いつかン
ネットで検索したり、図書館で本読みまくったが、初歩的なことなのか、ぜんぜんのっとらん。ハニングはわかるんだが

**デフォルトの名無しさん** · 2013/05/09(木) 19:26:40.27

あるシミュレーションで100万点ほどサンプルして
計算処理したいのですが、要素が多すぎますってエラーで
計算できません
初歩的なことですみませんが
回避方法教えてください

**デフォルトの名無しさん** · 2013/05/17(金) 17:46:16.03

>>129
解決しました

**デフォルトの名無しさん** · 2013/05/17(金) 22:58:43.48

一挙に読み込むのでなく
読み込みと処理とを適当に分割しました
とかいうんじゃねーだろな

**デフォルトの名無しさん** · 2013/05/19(日) 17:26:54.36

x-5cosx+3sinx-40=0
をニュートン法で求めたいのですが、matlabの数値計算の本を読んでもさっぱりで、どの様に考えたらよろしいでしょうか？
よろしくお願いします。

**デフォルトの名無しさん** · 2013/05/19(日) 19:01:17.44

すいません、自己解決しました。
簡単だった…

**デフォルトの名無しさん** · 2013/05/20(月) 18:43:09.02

http://www.wolframalpha.com/input/?i=Plot[x+-+5+Cos[x]+%2B+3+Sin[x]+-+40%2C+{x%2C+30%2C+50}]

グラフ見て概略把握して

http://www.wolframalpha.com/input/?i=FindRoot[x+-+5+Cos[x]+%2B+3+Sin[x]+-+40+%3D%3D+0%2C+{x%2C+34}]
http://www.wolframalpha.com/input/?i=FindRoot[x+-+5+Cos[x]+%2B+3+Sin[x]+-+40+%3D%3D+0%2C+{x%2C+38}]
http://www.wolframalpha.com/input/?i=FindRoot[x+-+5+Cos[x]+%2B+3+Sin[x]+-+40+%3D%3D+0%2C+{x%2C+42}]

34.1895
38.9164
42.2712

Mathematica

**デフォルトの名無しさん** · 2013/05/21(火) 04:39:15.71

>>131
for１回ごとに１００万この要素を持つ行列をそれぞれ計算させてそれを複数回まわしました
要素が多すぎと出ますが一応計算はできるみたいなので

**デフォルトの名無しさん** · 2013/05/21(火) 18:16:36.83

質問させていただきます。

現在、ある画像ファイルに対し、for文を用い少しずつ違う処理を加えた後、
各データをbmpとしてimwriteで保存しようとしているのですが、
ファイル名の末尾に
filename_1.bmp
filename_2.bmp
filename_3.bmp
...
などと、名前を適宜変えて保存するにはどうすればよいのでしょうか。

for p = 1 : 100
imwrite(data, 'filename_%d.bmp',p ,'bmp' );
end

のような感じで実行できると思っていましたが、どうやら駄目なようなので・・・。
よろしくお願いします。

**デフォルトの名無しさん** · 2013/05/21(火) 23:35:48.25

imwrite(data, ['filename_' num2str(p) '.bmp'] )

**デフォルトの名無しさん** · 2013/05/22(水) 13:02:07.37

imwrite(data, sprintf('filename_%d.bmp', p), 'bmp');
でもいける。

**デフォルトの名無しさん** · 2013/05/24(金) 00:46:13.71

>>137
>>138
ありがとうございます。

**デフォルトの名無しさん** · 2013/06/03(月) 14:42:31.05

質問です
ubuntuでlmstartからサーバ上のmatlabを起動していますが
matファイルやmファイルをダブルクリックで開こうと思ってもmatlab上でロードできません

matファイル等のプロパティからmatlabと関連付けしようとしましたが、一覧にmatlabが表示されず選択出来ませんでした

matファイル等とmatlabを関連付ける方法はありますか？

**デフォルトの名無しさん** · 2013/06/22(土) 13:18:08.13

質問です。よろしくお願いします。
オイラーの運動方程式をルンゲクッタ法で解きたいのですが、どのようにプログラムを組んだら良いのかが分からないという状態です。
よろしければ教えて頂けたら嬉しいです。

**デフォルトの名無しさん** · 2013/06/24(月) 10:55:53.56

質問です。imagescでデータをカラープロット
するときにy軸を対数スケールにしてなおかつy軸は
下から上に増加させたいと思います。自分でやるときに
axis('xy')を付加すると画像が消えます。どなたか
対処法を教えて下さい。よろしくお願いします。

**デフォルトの名無しさん** · 2013/07/04(木) NY:AN:NY.AN

パノラマ画像を作りたいのですが、どうしたらいいのでしょうか？？
stitchingという作業について詳しく知りたいです。

**デフォルトの名無しさん** · 2013/07/19(金) NY:AN:NY.AN

質問です。
文字列1つと数字データを2つ含む構造体を
エクセルに出力したいのですが、どうすれば良いでしょうか。

例；
data=struct('str','','prob',0,'sum',0);
このdataという構造体に例えば100個のデータが入っているとします。
data(1).str=string, data(1).prob=0.9,data(1).sum=5
data(2).str=string2, ...
...
data(100).str=string100,data(100).prob=0.1,data(100).sum=4
のような形です。

このdataという構造体をエクセルに出力するにはどうすれば良いでしょうか。

よろしくお願いします。

**デフォルトの名無しさん** · 2013/07/20(土) NY:AN:NY.AN

質問です。
畳み込み和を計算するユーザ定義関数を
function y = con(x)

y=[];
N=length(x);

for m=0:N-1;
tmp=0;
for m=0:N-1;
tmp=tmp+x(m)*h(n-m);
end
end
上記の様に書いてみたのですが、これであっているのでしょうか？
なにか別の書き方ありましたら教えて欲しいです。

**デフォルトの名無しさん** · 2013/07/20(土) NY:AN:NY.AN

もしよろしければどなたかお答えください。

matlabである時系列のパワースペクトルを白色化（whitening）して
それをまた時間領域に戻したいのですが、ネットをさがしてもでてきません。
どなたかやり方を教えて下さい。
よろしくお願い致します。

**デフォルトの名無しさん** · 2013/07/20(土) NY:AN:NY.AN

>>145
conv関数使って結果を確かめな

**デフォルトの名無しさん** · 2013/08/02(金) NY:AN:NY.AN

imagescで表示した画像中で、マウスを動かした軌跡の座標を取得して保存したいのですが、何か方法はありませんか？
function myBDCallback(src,eventdata)
set(src,'WindowButtonMotionFcn',@myBMCallback);
end
function myBMCallback(src,eventdata)
C = get(gca,'CurrentPoint');
C = round(C);
title(gca, ['(X,Y) = (',num2str(C(1,1)),',',num2str(C(1,2)),')']);
end
function myBUCallback(src,eventdata)
set(src,'WindowButtonMotionFcn','');
end
という風に関数を作ってみたのですが、座標をマトリクスに入れて吐き出すことが出来ません。
よろしくお願いします。

**デフォルトの名無しさん** · 2013/08/28(水) NY:AN:NY.AN

Matlab32bitをWindows8 64bitにインストールできるのでしょうか？

**デフォルトの名無しさん** · 2013/08/29(木) NY:AN:NY.AN

俺は7だから知らないけどとりあえずインストールしてみたら？

**デフォルトの名無しさん** · 2013/10/03(木) 16:31:08.11

複数画像を１つのfigure内に表示する事はできますか？
例えば、I=imread('***.bmp');と読み込んで
それをJ=rgb2gray(I);と処理しなおした後に
左右に並べて比較したいのですが

**デフォルトの名無しさん** · 2013/10/03(木) 22:40:54.71

subplot

**デフォルトの名無しさん** · 2013/10/03(木) 23:14:57.43

>>152
どうもありがとうございます
subplot(2,1,1),imshow(I); subplot(2,1,2),imshow(J);
これで並べる事はできましたが画像そのものが表示されるわけではなく
ヒストグラムのような座標軸が出てきました
画像とそのグレースケール化画像を並べるにはどう修正すればよいのでしょう？

**デフォルトの名無しさん** · 2013/10/04(金) 00:59:47.87

カルマンフィルターが導入されたのはいつのバージョンからだっけ？

**デフォルトの名無しさん** · 2013/10/04(金) 09:45:59.13

>>151
imshowpair(A,B,'montage')

**デフォルトの名無しさん** · 2013/10/16(水) 16:43:31.57

MATLABで、スクリプトでディレクトリ内のファイル名を変更していきたいんですけど、
ファイル名を変更する関数とかコマンドってありますか？。
教えてくださいm(__)m。

**デフォルトの名無しさん** · 2013/10/17(木) 21:39:30.01

simulinkでシミュレーションの固定ステップ0.01秒にしてるのに
simulink上のプログラムembeddedファイルでなぜか0.005秒ごとに
ループが回るんだけど、なんでろ？

**デフォルトの名無しさん** · 2013/10/26(土) 10:50:19.47

>>156
環境によって変わってくるけど、自分はsystemを使ってた。Unixだと
system(sprintf('mv %s %s', filename1, filename2));
みたいに。

**デフォルトの名無しさん** · 2013/10/28(月) 09:58:57.21

>>158
どうもありがとうございます！

**デフォルトの名無しさん** · 2013/11/03(日) 02:29:31.41

>>155
いやimshow([A,B])だろ

**デフォルトの名無しさん** · 2013/11/05(火) 23:45:57.61

画像サイズが同じならそれでもいいね

**デフォルトの名無しさん** · 2013/11/16(土) 07:17:50.91

matlab上でのswitch case文とif else文は
分岐の数が同じ場合どちらが速度が速いですか？

**デフォルトの名無しさん** · 2013/11/24(日) 13:24:30.09

画像解析で色を検出して、例えば青系等で薄い青は低く、濃くなってくと高くなるっていうグラフを作りたいんですが、どうやって作っていけばいいんですかね
ガウス分布で表して、ポテンシャル場にする感じです

日本語下手で申し訳ない

**デフォルトの名無しさん** · 2013/11/24(日) 19:23:37.72

>>163
色の濃淡はrgb色空間向いてないからluvとかにもっていってあとは指定色をしきい値処理して分布みたらいい気がする
ただやりたいことがあまり良くわかんないから参考にならないのかもしれない

**デフォルトの名無しさん** · 2013/11/25(月) 22:43:43.87

S-Functionについての質問です．

InputRealPtrsType u0 = ssGetInputPortRealSignalPtrs(S,0);
real_T *y0 = ssGetOutputPortRealSignal(S,0);

配列u0を配列y0にコピーする関数を別のcファイルに定義しようとしています．
プロトタイプ宣言と関数の引数はどのように渡せばいいのでしょうか？

プロトタイプ宣言を
extern void ORIGINALFCN(real_T u0,real_T y0[]);

関数呼び出しを
ORIGINALFCN(*u0[0],y0);

としているのですが”error C2109: subscript requires array or pointer type”
となり困ってます　御助力いただけないでしょうか

**デフォルトの名無しさん** · 2013/11/27(水) 10:46:04.28

165です
自己解決しましたー

**デフォルトの名無しさん** · 2013/11/28(木) 21:39:25.53

ある座標(x1,y1)を通る直線を座標を中心に回転行列を用いて回転させたいんですが、
どうすればいいのか教えてください

**デフォルトの名無しさん** · 2013/11/28(木) 23:48:23.88

(x1,y1)が中心となるように座標系（直線の式）を平行移動する
任意の角度で座標系を回転させる

**デフォルトの名無しさん** · 2013/12/01(日) 22:49:24.34

離散データを特殊関数でフィッティングする
やり方を教えて下さい。
よろしくお願いします。

**デフォルトの名無しさん** · 2013/12/01(日) 23:36:03.84

lsqcurvefit

**デフォルトの名無しさん** · 2013/12/04(水) 04:11:20.84

mesh関数のview(2)で
連続的に変化する複数のグラフを作ったのですが
これをgif化するにはどうしたらできますか？

**デフォルトの名無しさん** · 2013/12/04(水) 04:11:53.13

グラフ→マップ図

**デフォルトの名無しさん** · 2013/12/05(木) 17:39:46.79

170さん、ありがとうございます。
lsqcurvefitはネットで使い方が
書いてあるのですが、使い方がよくわかりません。
ちなみに、xdataとydataは同じ大きさのベクトル
なのですが、初期値x0の意味がよくわかりません。
ちなみに、ベッセル関数で離散データをフィッティング
したいです。

**デフォルトの名無しさん** · 2013/12/08(日) 02:15:45.83

>>168
>>167です。遅れましたが、ありがとうございました

**デフォルトの名無しさん** · 2013/12/16(月) 18:16:03.97

離散データの極値を求めたいのですが
ネットで調べるとSymbolic Math Toolbox
が必要とでてきました。
Symbolic Math Toolboxがなくても微分
できますか？

**デフォルトの名無しさん** · 2013/12/17(火) 01:02:35.46

複雑じゃなきゃ関数の微分くらい手計算でやれ
面倒なら数値微分で近似

**デフォルトの名無しさん** · 2013/12/17(火) 01:38:40.35

今時微分積分もスマホのCASでも出来るよ
ネットでPCならSage math cloudがお勧め

**デフォルトの名無しさん** · 2013/12/17(火) 11:41:27.66

>>175
オイラー法でいけんじゃね？しらないけどww

**デフォルトの名無しさん** · 2013/12/18(水) 01:04:37.19

単純に数式計算がしたいならフリーのmaximaでも使えば？

**デフォルトの名無しさん** · 2013/12/19(木) 04:20:30.18

学生版では配列のサイズに制限があるようですが
オプション関連でこの配列制限が影響するオプションはありますか？
SimscapenやSimElectronics，SimPowerSystemsに影響があるか知りたいです

**デフォルトの名無しさん** · 2013/12/20(金) 21:58:08.76

ある行列testからゼロ以上の最小値を得たい場合、
findをどのように使えば良いでしょうか。

[r, c] = find(test == min(test(:)) && test >= 0);

という書き方はできないようで…

**デフォルトの名無しさん** · 2013/12/20(金) 22:21:38.61

min( test(test>=0) )

**デフォルトの名無しさん** · 2013/12/22(日) 02:39:15.24

>>182
大変参考になりました
ありがとうございます

**デフォルトの名無しさん** · 2014/01/26(日) 22:02:35.70

ポップアップメニューの中身を可変にしたいんですけど可能ですか？
ファイルを読み込み、それによってポップアップメニューの中身の項目数が
３個だったり、１０個だったりと変化するのですが、うまくいきません。

あとfigure1からfigure2を呼び出すことはできますか？
figure1がメインの画面、figure2がプロパティの画面になります。
引数付きで呼び出したいのですが、やり方がわかりません。

わかる人、よろしくお願いします。

**デフォルトの名無しさん** · 2014/01/26(日) 23:24:39.34

ゴミゴミゴミゴミゴミゴミゴミゴミゴミゴミゴミゴミゴミゴミゴミゴミゴミゴミゴミゴミゴミゴミゴミｗｗｗｗｗｗｗｗｗｗ

**デフォルトの名無しさん** · 2014/01/27(月) 16:18:59.66

関数定義の方法について非常に困っています。どうかお助けください。

入力変数にx,y(サイズN)を取るN元連立非線形方程式F_i(x,y)=0 (i=1,2,...N)があり、
ワークスペースに測定によって得られた既知の変数xの値x_measuredがあります
この時F_i(x_measured,y)=0は(yについて)fsolveで解くことができるかと思うのですが、
F_i(x_measured,y)=0の上手い定義の仕方が分かりません

最も愚直にはワークスペース上のx_measuredの値を見ながら、
function F = fun(y)
F = [f_1(x_measured,y); f_2(x_measured,y); ... ; f_N(x_measured,y)]
end
と関数定義してスクリプトで"fsolve(@myfun, y0, ... )"とすれば一応解けますが、
これではx_measuredが変わる度にfunctionの定義しなおさなくてはなりません
(実際問題Nは非常に大きく、また入力はx,yだけではないのでほぼ不可能です)

できれば関数の定義は次のように
function F = fun(x,y)
F = [f_1(x,y); f_2(x,y); ... ; f_N(x,y)]
end
としておいて、fsolveに受け渡す際に"fsolve(@myfun(x_measured,y), y0, ... )"
のような書き方ができれば非常に嬉しいのですが、可能でしょうか？

もし可能であれば、関連するドキュメントや簡単な例など紹介していただければ幸いです。

**デフォルトの名無しさん** · 2014/01/27(月) 22:28:34.31

関連するドキュメントと簡単な例
ttp://www.mathworks.co.jp/jp/help/optim/ug/fsolve.html

**デフォルトの名無しさん** · 2014/01/28(火) 00:30:34.29

調べたところ、無名関数を使えば良いっぽいです
何とか解決できそうなので、もう少し頑張ってみます
ttp://www.mathworks.com/matlabcentral/answers/28728-fsolve-function-with-values-read-from-the-workspace

**デフォルトの名無しさん** · 2014/02/01(土) 12:53:34.55

FFTしたwavファイルの周波数をグラフにプロットしているのですが縦軸が五線譜に見えるように
「縦軸を(12√2)*n毎に目盛りを配置」尚且つ「目盛り同士が等間隔」
になるよう表示するやり方はありますか？

またmatlabではなくscilabでプログラムしているので差異があればそこも教えてほしいです

答えていただける方がいましたらもう少し詳細をお伝えしますのでお願いします

**デフォルトの名無しさん** · 2014/02/01(土) 13:30:36.68

matlabのプロットならわかるがscilabは知らん

**デフォルトの名無しさん** · 2014/02/01(土) 13:36:51.58

>>190
分かりました
ではmatlabのままで構いませんので教えていただけないでしょうか？

**デフォルトの名無しさん** · 2014/02/01(土) 16:39:06.25

目盛り幅は基本的に等間隔
目盛りを手動で指定するにはytickをいじる
ttp://sach1o.blog80.fc2.com/blog-entry-58.html

**デフォルトの名無しさん** · 2014/02/08(土) 13:13:23.34

仕事でやむを得ずMATLABを使用しているが、
言語仕様が酷すぎてストレスが溜まる

**デフォルトの名無しさん** · 2014/02/08(土) 14:46:16.67

どういう点がストレス？

**デフォルトの名無しさん** · 2014/02/09(日) 11:07:55.88

数値計算にはかなり広く使いやすい仕様だとおもうけどな

**デフォルトの名無しさん** · 2014/02/12(水) 20:46:40.36

特にセル配列周りは頓珍漢仕様

**デフォルトの名無しさん** · 2014/02/12(水) 21:25:06.09

たしかに仕様は変だけどセル配列ってそんなに多用するかな
それよりも文字列をシングルクォートで囲むことが変態的
行列転置と被るからダブルクォートにしてくれりゃいいのに。どうせ使ってないんだし

**デフォルトの名無しさん** · 2014/02/13(木) 23:56:01.06

コンパイラでスタンドアロンアプリを作ろうとすると
セル配列を使う場面は結構あるのよ

**デフォルトの名無しさん** · 2014/03/04(火) 16:43:10.80

MATLAB Coderを使ってC言語にする時、複数のfunctionAがfunctionBを呼び出す場合はどうやればいいですか？
教えてください。

例えば、fprintf(1,'Hello¥n');を実行するfunction hello()をmainから呼び出すとしたら
codegen mainとすると
「このテキストは、空でない最上位の式を含んでいます。それはスクリプトのようです。」とエラーが出てしまいます。

**デフォルトの名無しさん** · 2014/03/25(火) 03:00:28.60

質問です。
例えば、
A=cat(3, [1 2 0 ; 0 5 6 ; 0 8 9], [1 1 1 ; 1 1 1 ; 1 1 1], [2 2 2 ; 2 2 2 ; 2 2 2])

という行列Aがあるとして、
index=find(A(:,:,1)==0)

このindexを利用して、A(:,:,1)の0の位置と同じ、A(:,:,2)およびA(:,:,3)の要素をxに置き換えて、
A=
1 2 x
x 5 6
x 8 9

1 1 x
x 1 1
x 1 1

2 2 x
x 2 2
x 2 2

という形にしたいのですが、うまい方法が思いつきません。

ご助言お願いします。

**デフォルトの名無しさん** · 2014/04/01(火) 22:46:42.16

初心者ですが、ＭＡＴＬＡＢをつかってリアプノフ指数を求めることはできるでしょうか？

**デフォルトの名無しさん** · 2014/04/02(水) 04:38:59.48

>>200
亀レスだけど
[row,col] = find(A(:,:,1)==0);
A(row,col,:) = x;

**デフォルトの名無しさん** · 2014/04/14(月) 01:11:48.38

x軸の上下で異なる目盛り(例えば、下：秒、上：分、など)を付ける方法はあるのでしょうか？
もしご存知の方がいらっしゃいましたら、お教え頂けますと大変助かります。
初心者の質問で恐縮ですが、どうぞよろしくお願いします。

**デフォルトの名無しさん** · 2014/04/14(月) 22:43:23.43

matlabって固有値が重解もつときの固有ベクトル正しく求められるの？
Ver5.3くらいの時は求められないことが普通にあったけど

**デフォルトの名無しさん** · 2014/04/28(月) 17:03:08.35

MATLABでSNRを求めないといけないんですが全然進みません。
一応ここ
www.mathworks.com/matlabcentral/fileexchange/35479-stochastic-resonance-in-the-duffing-oscillator-with-matlab/content/runduffsr.m
のプログラムを基にしたいんですが・・・

**デフォルトの名無しさん** · 2014/07/22(火) 00:00:42.43

scilabからの関数のプログラムの移植で困ってます。

元のscilabプログラムは関数の引数に伝達関数を入力している部分がありますが、
matlabではうまくいきません。

必要なのは伝達関数の分母分子のsの係数だけなので、
それを抽出した配列を引数にしようとしてもダメでした。

どうすればよいでしょうか。
よろしくお願いします。

**デフォルトの名無しさん** · 2014/07/23(水) 09:03:30.02

>>206
Matlabでも同じことができるはずだけど、
コードをさらしてくれないと何も言えね。

**デフォルトの名無しさん** · 2014/07/26(土) 12:55:18.49

ここはsimulinkの質問もOKなのですか？

**デフォルトの名無しさん** · 2014/07/28(月) 22:21:24.90

forの中にこのような表記で計算結果を連続させアニメーションを表示させています。
この表示中に右ドラッグで自由に回転させたいのですが、回転しようとすると直ぐに元の位置に戻ってしまいます。
アニメーションを表示させながら回転させられる方法はないでしょうか？
drawlater();
delete();
surf(X,Y,Z);
xset('colormap',jetcolormap(32))
drawnow();
sleep(1);

**デフォルトの名無しさん** · 2014/07/28(月) 22:59:32.75

>>207
function RetVal = calcParam(Gp, cont_model, ref_model)

となっていて、Gpがtf型です。
2つ目、3つ目の引数は整数。

関数を呼び出す際にGpにtf型の値を入れると、
『関数 'calcParam' (タイプ'tf' の入力引数) が未定義です。』
となります。

tf型の引数を与えられるんでしょうか。
それとも別の方法になるのでしょうか。

**デフォルトの名無しさん** · 2014/07/29(火) 03:23:52.76

ＭＡＴＬＡＢで磁界のシミュレーションをしていきたいと考えています。
しかしネットで調べても
「マトラボで磁界のシミュレーションができるよ」
とは書いてあるものの具体的な方法を記載しているサイトが見つかりません

サイトに限らず書籍などでもいいので、参考になるものを知っている方がいらしたら教えてください。
よろしくお願いします。

**206=210** · 2014/07/30(水) 22:46:46.11

functionを使わないことで、引数の件は本質的ではないですが解決(？)しました。

ところで、scilabではrootsで3次方程式を解いていて、
解がすべて小数で求められているのですが、
matlabでsolveコマンドでは解がatanやらcosやらが混ざった式の解として求められ、
この求められた式を計算してもscilabの解と一致しません。

3次方程式はテキストファイル上にscilab、matlabの2つを上下に並べて比較し、
同じであることを確認してあります。

matlabの方程式の変数はsymsで定義してあります。

何か考えられる原因はあるでしょうか。
あと解をすべて小数表示で求める方法はあるでしょうか。

**デフォルトの名無しさん** · 2014/08/04(月) 17:32:53.13

Octave環境において、他ファイルに記述した関数を利用するために
source ./functions.m

とsourceコマンドを用いていたのですが、
MATLAB環境に移行するにはどのように記述すればいいのでしょうか？

よろしくお願いします。

**デフォルトの名無しさん** · 2014/08/22(金) 12:01:14.93

scilab+SIVPで画像処理をしてるけど処理が激遅で困ってます
┌──┐
│ ○　│→┌───┐
└──┘　 └───┘
↑みたいな480×640の画像で、真ん中の円の部分を長方形の画像に展開したいんですが、
1ドットずつ長方形の枠に色情報を移すようなループ処理でやってるせいか超遅い・・・

円の半径は480/2の240として、円の最外周から30ドットくらいまでを長方形の画像に
展開するのに高速で終わるようにするにはどうしたら良いでしょうか。
今は下みたいな処理をしています。

MARUは480*640*3の配列で、SIKAKUは30*1500*3の配列
MARUで指定する座標の色情報を一つずつSIKAKUの配列に入れて行く
r=240　//半径r
while(r>210)　 //半径の最大から最小値まで繰り返す（行置換）
　　a=360　　　　//360°ぐるりと回りながらドット取得
　　j=1　　　　　//長方形展開後の色情報を格納する列番号
　　while(a>0) 　//角度を360から0まで繰り返す（列置換）
　　　　x=round(px+r*cosd(a))　　//○画像の抽出したいドットの横方向座標取得
　　　　y=round(py+r*sind(a))　　//○画像の抽出したいドットの縦方向座標取得
　　　　SIKAKU(i,j,1)=MARU(y,x,1)　　//R・G・Bそれぞれ円の座標の色情報を
　　　　SIKAKU(i,j,2)=MARU(y,x,2)　　//長方形の画像の指定する場所に格納
　　　　SIKAKU(i,j,3)=MARU(y,x,3)
　　　　j=j+1　//長方形展開後の次の色情報を格納する場所の列指定
　　　　a=a-0.24　//0.24°角度を移動　(半径480ドットの円の円周≒1500ドット、360/1500=0.24)
　　end
　　r=r-1　//最外周から1ドットずつ内側に移動
　　i=i+1　//長方形の画像の次の色情報を格納する場所の行指定
end

この後さらに複数の長方形に展開した画像を結合して長い円筒にカメラを突っ込んだような画像から、
円筒を唐竹割りしたような展開画像を作りたいけど良い方法有りますでしょうか。

**デフォルトの名無しさん** · 2014/08/22(金) 12:25:56.14

>>213
MATLABなら関数が置いてあるパスを登録する

OCTAVEも同様に.octavercに
パスの一覧登録しておく方法がある

**214** · 2014/08/22(金) 16:30:56.24

まだループが有るけど大分早くなったので途中経過・・・
後は動画から連続する画像取得とその展開、展開画像の結合か・・・先は長い(´･ω･｀)

i=1
r=240　//半径r
a=0:(360/1500):360
while(r>210)
　　MARU_R=DR(round(py+r*sind(a)),round(px+r*cosd(a)))
　　MARU_B=DG(round(py+r*sind(a)),round(px+r*cosd(a)))
　　MARU_G=DB(round(py+r*sind(a)),round(px+r*cosd(a)))

　　MARU_R=diag(MARU_R)'
　　MARU_B=diag(MARU_B)'
　　MARU_G=diag(MARU_G)'

　　SIKAKU(i,ANS_yoko,1)=sum(MARU_R,'r')
　　SIKAKU(i,ANS_yoko,2)=sum(MARU_B,'r')
　　SIKAKU(i,ANS_yoko,3)=sum(MARU_G,'r')

　　r=r-1
　　i=i+1
end

**デフォルトの名無しさん** · 2014/08/22(金) 23:53:25.04

quiver関数を使っているのに矢印がプロットされません
ためしにscale=0にして自動スケーリングを使用しないよう設定したら
矢印が表示されました。しかしすごい大きな矢印とちいさい矢印が複数ありました
これは大きいのと小さいのが入り乱れてるせいでスケールすると表示されなくなってしまうものなのでしょうか？
ネットで調べてもいまいちquiverを掘り下げて説明しているところがなかったので
こちらで質問させていただきます。

**デフォルトの名無しさん** · 2014/10/11(土) 19:15:16.44

textscanを使ってテキストファイル(.txt)のデータを読み取った後、読み取ったテキストファイルを削除したいと考えています
しかし、textscan後にdeleteを使ったり、手動で削除しようとすると「MATLABによってファイルは開かれているため、操作を完了できません。」
と表示され、削除できません。
現状では一度MATLABを再起動してからテキストファイルの削除を行っているのですが、再起動せずにMATLABからテキストファイルへの読み込みを終了させるコマンドはないでしょうか

**デフォルトの名無しさん** · 2014/10/12(日) 13:17:28.68

>>218
fclose()やってないからじゃないの？

**デフォルトの名無しさん** · 2014/10/12(日) 15:13:51.74

>>219
できました！！
ありがとう！！

**デフォルトの名無しさん** · 2014/10/20(月) 22:10:45.56

>>146
filter(lpc(x,3),1,x)

**デフォルトの名無しさん** · 2014/10/26(日) 12:30:23.03

シリアル日付値をx軸にすると、xlimで指定した範囲ぴったりにならないことがあります。
どうすればいいのか、どなたかご教示いただければ幸いです。
（例）
t = datenum(2014,10,26,12,30,0:1/100:60-1/100);
x = rand(1,6000);
plot(t,x);
xlim([t(1) t(100)]);
datetick('x','MM:SS','keeplimits');

**デフォルトの名無しさん** · 2014/10/31(金) 18:24:26.17

matlabのmaxは最初に見つけられたインデックスを返すのですが、最後に
見つけられたインデックスをかえすようにしたいのですが、方法が分かりません。
よかったら教えて下さい。

**デフォルトの名無しさん** · 2014/11/12(水) 22:37:23.56

maxの引数になっている変数をA(end:-1：1)ってなのが一番楽じゃないか
後ろから数えたインデックスになってるけど全体サイズ分かってれば前からの位置もすぐ分かるし
一週間以上前のだからもう見てなさそうだけど

**デフォルトの名無しさん** · 2014/11/13(木) 00:28:03.60

fliplrとかflipudとか使ってもなんとかなりそう

**デフォルトの名無しさん** · 2015/01/25(日) 18:58:32.20

FDAtool使ったFIRフィルタの設計って、設定したフィルタ特性や量子化パラメータでは
どんなに頑張っても実現出来ないときに、「出来ませんでした><」って教えてくれないの？
自分で特性グラフみて判断するしかない？？

**デフォルトの名無しさん** · 2015/03/06(金) 23:20:44.66

matlabで振動にFFTかけてるんですけど、
周波数帯域毎のパワーを時間軸に沿って表示する、というのを出来る関数は無いんですか？

**デフォルトの名無しさん** · 2015/03/06(金) 23:28:35.80

ウェーブレット使えばいいじゃん

**デフォルトの名無しさん** · 2015/03/08(日) 19:40:16.17

それか短時間フーリエ変換かな？
短時間フーリエ変換なら
specrogram
って関数がＳignal Processing Toolboxにある

関数なくても短時間フーリエ変換程度ならなら
自作するかネットに転がってるだろ？

**デフォルトの名無しさん** · 2015/03/09(月) 20:17:54.85

ありがとうございます。
挙げて頂いたワードを元に調べた所、非定常振動の解析にはウェーブレットが向いているようなので、ウェーブレットでやってみようと思います。

**デフォルトの名無しさん** · 2015/03/27(金) 22:24:22.78

MATLABって過去のバージョンとか使えないんですか？
今からだと2015しか使えない？

**デフォルトの名無しさん** · 2015/03/28(土) 16:44:57.12

自分は過去バージョンダウンロードできるけど、
以前から保守契約してるからかもしれない。

**デフォルトの名無しさん** · 2015/04/08(水) 22:48:39.72

すんません、y軸のメモリを10の二乗等で表示したいですが、どうしても"1e2"となってしまいます。
どうしても"10^2"のように表示させたいのですが、設定方法をご教示頂けますか？

**デフォルトの名無しさん** · 2015/04/09(木) 23:36:13.70

自分はset(gca, 'xtick', ～)とset(gca, 'xticklabel', ～)使ってる。

**デフォルトの名無しさん** · 2015/04/10(金) 02:57:47.83

>>234
返信ありがとうございます。

＞自分はset(gca, 'xtick', ～)とset(gca, 'xticklabel', ～)使ってる。

set(gca, 'xtick', ～)の～の中身を教えて頂けますか？そこが分かればおそらく解決です。

**デフォルトの名無しさん** · 2015/04/12(日) 18:58:23.99

>>235
みんなそういう使い方しませんし、そのようなことはマトラボじゃできません。
私の知る限りではgnuplot ならばできたはずですよ。
マトラボはあくまでも数値計算に特化したソフトウェアなので、グラフの調整などは不得意です。

**デフォルトの名無しさん** · 2015/04/12(日) 19:54:49.79

>>235
たとえば、
set(gca, 'xtick', [1 10 100 1000]);
set(gca, 'xticklabel', {'10^0', '10^1', '10^2', '10^3'});
とやればOK。

**デフォルトの名無しさん** · 2015/04/13(月) 00:10:27.46

>>237
アドバイスありがとうございます。
しかしそのやり方だと文字通り"10^2"のように表されるんですね。
10の上付き表示はできないものなのでしょうか？

**デフォルトの名無しさん** · 2015/04/14(火) 01:52:47.90

>>238
だーかーらーマトラボはそんな変な設定しても結局あんたがしたいような表示はできません。
前も言ったけど、マトラボは計算に特化したソフトウェアだからね？
大体、set(gca,xrick,[........])でいちいち中身書くなんてアホでしょ？
Gnuplot はそんなとこしなくても%Lみたく書けば、OKだよ。
ここにいるアホオタクは知らないだけのオタクだから、気をつけな～（＾ω＾）/

**デフォルトの名無しさん** · 2015/04/18(土) 13:18:46.48

>>238
自分の環境（Matlab 2015a）では上付き文字で表示されるので問題ないのですが、
2014バージョンあたりでグラフィックス周りが刷新されたのでそのせいかもしれません。
set(gca, 'TickLabelInterpreter', 'tex');
が使えれば、TeX文法でティックラベルが表記できます。

**デフォルトの名無しさん** · 2015/04/18(土) 13:19:27.60

ちなみに、TeX文法つかいたくないときは
set(gca, 'TickLabelInterpreter', 'none');
です。

**デフォルトの名無しさん** · 2015/04/18(土) 13:19:48.22

ちなみにTeX文法をつかいたくないときは
set(gca, 'TickLabelInterpreter', 'tex');
です。

**デフォルトの名無しさん** · 2015/04/18(土) 18:09:50.51

二重ごめんなさい。

**デフォルトの名無しさん** · 2015/04/26(日) 13:18:58.54

他にmatalabスレが見つからなかったので教えて頂けないでしょうか？
個人ライセンス版の購入を考えているのですが、商用利用不可である以外は
通常のmatlabライセンスと機能等に変わりは無いのでしょうか？
またライセンスサーバーは特定のpcに紐付けされるかたちになるのでしょうか？

**デフォルトの名無しさん** · 2015/04/26(日) 15:54:04.25

おまえの知りたい事は全てMathWorksのホムペに書いてある
超初歩的な疑問位少しは自分で探す努力しろ
1問目yes
2問目yes

**デフォルトの名無しさん** · 2015/04/27(月) 02:56:23.14

レスありがとう御座います！
何だかんだで親切すぎるw
mathworksって振り込みクレカしか駄目なんですね。。。
さすがにVプリカはダメだろうし、クレカ入力は抵抗あるので一時見送りますorz

**デフォルトの名無しさん** · 2015/05/01(金) 12:40:31.09

ＪＮＢのワンタイムデビットで支払えないのかな。

**デフォルトの名無しさん** · 2015/06/04(木) 00:03:13.97

test

**デフォルトの名無しさん** · 2015/06/04(木) 00:06:36.97

実験データから伝達関数を求める方法を教えてください。
具体的には、ハンマリング試験した加振力と加速度の時系列データがあって、
それらを用いて伝達関数(ボード線図)を描きたいです。
そして、その伝達関数を2次の振動系に近似し、変数を同定する方法も教えていただきたいです。

また、その伝達関数の実部、虚部に分解したコクアド線図も描きたいです。

どなたかよろしくお願いします。

**デフォルトの名無しさん** · 2015/06/05(金) 03:25:11.16

>>249
matlab以前の問題もするけど・・・
一応聞くけどコードに落とせる程度に手計算は終わってるんだよね？

**249** · 2015/06/05(金) 23:03:31.24

>>250
手計算とは具体的にどういったものでしょうか。
実験データ(エクセル、テキスト)の時系列データから
そのまま関数を打つだけで伝達関数が出たりするのではないかと思っています。
探したりしていろいろ試しているのですがうまくいきません。

**249** · 2015/06/09(火) 21:42:07.28

解決しました

**デフォルトの名無しさん** · 2015/06/10(水) 01:40:41.68

>>252
氏ね
何をどうしたら解決したのかくらい書けや
ゴキが掲示板荒らしたとしか考えられん
二度と書き込むな、クズ

**デフォルトの名無しさん** · 2015/06/10(水) 03:27:36.11

小野測器とか見れば参考になるよって解決したのか
ダブルタップはしっかり省いてね

**デフォルトの名無しさん** · 2015/06/11(木) 01:14:15.04

ただ本人が解決したから、もう答えなくてもいいですよって意思表示しただけなのに
きれるとか短気すぎんだろ

**デフォルトの名無しさん** · 2015/06/12(金) 08:17:52.04

>>255
短気すぎないわ。
そんな奴ただのスレ汚しのカスだから。
生きてる意味ないよ。端的にゴミゴキと言い表したいね。
ただしこの話題で盛り上がる気はないので以後、スルーで

**デフォルトの名無しさん** · 2015/06/12(金) 21:56:58.78

かすだなぁ

**デフォルトの名無しさん** · 2015/07/05(日) 03:51:35.78

すんません。

Plotyyで２つの異なるy軸使ってグラフを書いたんだけど、y軸の色がそれぞれ曲線の色になってます。
軸はすべて黒にしたいのですが、なにかいい方法ありますか

**デフォルトの名無しさん** · 2015/07/08(水) 19:48:27.68

>>258
plotyyってそんなふうになるんだっけ・・・？
面倒だけどset(gca,～)だったかで、直接軸の色を指定しちゃったら駄目かな

**デフォルトの名無しさん** · 2015/07/09(木) 06:32:17.79

>>259
ありがとうございます。
ちょっとそれでやってみますね。

ちなみにoctaveでも確認してみましたが、同様にこんな現象起きます

**デフォルトの名無しさん** · 2015/07/10(金) 03:04:35.50

>>260
グラフまわりで困ったら取り合えずset使っとけば良いよ
何でも出来るから慣れとくと役に立つかも

**デフォルトの名無しさん** · 2015/07/12(日) 14:51:51.14

表示しているfigureを新しいwindowで開きつつコピーする方法は無いでしょうか？

**デフォルトの名無しさん** · 2015/08/24(月) 21:58:13.27

Windows環境で作成したMATLABのスクリプトをlinux環境で動作させたいです。

Win環境で動作するexeは生成できましたが、同様にlinuxで動作するバイナリファイルは作成できますか？

**デフォルトの名無しさん** · 2015/08/31(月) 17:19:06.69

すんません、複数のコメントアウト行と列の数が途中から異なる数値データを読み込む方法ってどうすればよいですか？
具体的には以下のようなデータを読み込みたいです。

== this is comment line
2.56 2.67 2.34 3.45
2.56 2.67 2.34 3.45
2.56 2.67 2.34 3.45
2.56 2.67 2.34 3.45
2.56 2.67 2.34 3.45
== this is comment line
2.56 2.67 2.34 3.45 2.56 2.67 2.34 3.45
2.56 2.67 2.34 3.45 2.56 2.67 2.34 3.45
2.56 2.67 2.34 3.45 2.56 2.67 2.34 3.45
2.56 2.67 2.34 3.45 2.56 2.67 2.34 3.45
2.56 2.67 2.34 3.45 2.56 2.67 2.34 3.45
2.56 2.67 2.34 3.45 2.56 2.67 2.34 3.45

textscan を使買ってみてもなかなかうまくいきません。

Fid = fopen('filename', 'r');
N = 8;
While ~feof(did)
C = textscan(Fid, '%f %f %f %f %f %f %f %f', N, commentStyle, '==');
End

として, c{1}(一列目の数値データを見たい)をみても 2.56 が11行出てきません。
どこを改善すればよいでしょうか？

**デフォルトの名無しさん** · 2015/09/04(金) 07:13:32.01

質問です。

X = [0.1:0.1:1.0]
Y = [a(1) ...a(i)] i=10
Z = [b(1) ...b(I)]
の行列があった時、X-Y平面上にZの高さ情報を色で表した図って書けますか？
イメージ的には plot3 で書いたものをxy平面に射影したもので、zの高さ情報を色で表したものといった感じです。

**デフォルトの名無しさん** · 2015/09/04(金) 07:39:55.95

>>265
pcolor
Zは行列になってないとダメかもしれない

**デフォルトの名無しさん** · 2015/10/02(金) 00:17:11.01

conterfで描いた図に特定の値の等高線の太さだけを変えることって、できますか？

**デフォルトの名無しさん** · 2015/10/04(日) 21:18:03.30

>>267に関して誰かいませんか？？？

**デフォルトの名無しさん** · 2015/10/05(月) 21:19:42.11

>>267
しょうがないマジレスしてやるか
matlab 特定の等高線レベルの強調表示
でググってみな

hold on、hold off 使って線の太さと
等高線描画の範囲変えて図を重ねればいい

**デフォルトの名無しさん** · 2015/10/09(金) 02:20:54.07

>>269
全然反映されなちんです。。。
どうしまひょ？

**デフォルトの名無しさん** · 2015/10/09(金) 07:40:05.99

Mathworksのサイトにサンプルあるんだ
から、先ずそれが動くか確認
それが問題無いなら自分のコードが
間違ってんだろ？

**デフォルトの名無しさん** · 2015/10/09(金) 22:47:17.87

>>271
マトラボのページはいつも peak関数を使うじゃないですか…
それが本当に煩雑にしてるなぁと思うんですよね。
私の場合それぞれ独立したデータ x, y, z でコンター作ってるんです。

で、z = z0 となるときの等高線レベルのみを強調したいと考えておるのです。

**デフォルトの名無しさん** · 2015/10/10(土) 06:45:27.91

情報後出し萎えた
もう勝手にやってくれ

**デフォルトの名無しさん** · 2015/10/10(土) 10:09:02.17

>>273
そんなこと仰らずに助けてください。
お願いいたします

**デフォルトの名無しさん** · 2015/10/10(土) 15:52:42.74

conterfって使ったこと無いんだけど、そのz0の座標だけ抜き出して重ね書き＆太さ調整できんの？
どういうデータのグラフを描こうとしてるのかわからんけど、等高線で示してokな程度のものなら
適当に散布図とかでプロットの色・サイズを弄ったりしてもいいんじゃない？

**デフォルトの名無しさん** · 2015/10/10(土) 16:29:03.30

先ずconterfなんて無えから…

**デフォルトの名無しさん** · 2015/10/10(土) 16:51:10.71

あーcontourfか orz
普段はimagescばっかだなぁ

**デフォルトの名無しさん** · 2015/10/22(木) 00:34:00.05

x = 0.1:0.1:1.0
y = [2 4 6 7 9 13 26 25 67 89]
みたいなデータがあって
y >= 10 なら破線、10以下なら実線てな具合でプロットすることできる関数ってあったっけか？
条件分岐とplotを使って書くと、y = 9と13の間で空白ができちまってなんともダサいのよね(9 < y < 13 までプロットされず空白になる)

**デフォルトの名無しさん** · 2015/10/23(金) 03:57:18.87

条件ダブらせればいいんでないの？
それか10の値内装してマーカー有り実線グラフの上からマーカー無し白点線とか

**デフォルトの名無しさん** · 2015/10/28(水) 23:12:16.86

互換なソフトがあるのになぜ有料のものを使うのですか？

**デフォルトの名無しさん** · 2015/10/30(金) 03:05:26.91

matlab使えば分かるよ
ってか言葉に互換性が有るだけで使い心地は全然違う

**デフォルトの名無しさん** · 2015/10/30(金) 08:01:42.18

>>281
どのくらい違うの？
今octave使ってるんだけど

**デフォルトの名無しさん** · 2015/10/30(金) 18:11:03.37

Pythonのnumpy使ってるけど、行列計算はmatlabが圧倒的に楽に書けるね

**デフォルトの名無しさん** · 2015/10/31(土) 10:00:16.08

その楽ってのがわからんのよねー

マトラボはなんだがオプションが多すぎてなに追加したらいいのかわからんからなぁ

**デフォルトの名無しさん** · 2015/10/31(土) 15:24:38.19

取りあえずデフォのmatlabだけで大体問題ないと思う
信号処理ならsignal~
toolbox追加するとか、後追いで対応できるしね
というか、やっぱ本当に必要な関数は自分で作っちゃたほうが安心な気がするけどな
homeライセンスならoffice感覚な値段だし試ししてみたら？

**デフォルトの名無しさん** · 2015/10/31(土) 21:29:01.17

>>285
そうしてみようかしら。。。

octave使いなんだけど、やはりマトラボさんの方が絶対いい！って言える？

ごめんね、しつこくて。
(背中押して欲しい)

**デフォルトの名無しさん** · 2015/11/01(日) 04:12:05.31

浮気しろ＾＾

**デフォルトの名無しさん** · 2015/11/01(日) 17:31:34.94

octaveだとデフォルトでは使えないfuctionが結構あるけど、
追加しようと思えばできるから個人的にはoctaveでも問題ないけどなぁ

**デフォルトの名無しさん** · 2015/11/04(水) 22:57:37.38

MatlabはIDEが良くなってきたのとドキュメントがしっかりしてるところが魅力。
Matlabに10年ほど課金してきたけど、もうドキュメント読まなくても覚えちゃってるし、IDEなくてもEmacsでいいし、そろそろOctaveでいいやってなってる。

**デフォルトの名無しさん** · 2015/11/05(木) 02:48:42.34

まぁぶっちゃけoctave自分から選ぶような層はあえてmatlab選ぶ必要なんて無いよねぇ
cでもpythonでも金のかからない適当な方法選んでばっちり解析出来るイメージ
自分は仕事でどっぷりすぎてmatlab以外の選択肢は考えなかった口だわ；；

**デフォルトの名無しさん** · 2015/11/05(木) 07:00:23.05

matlab 使う仕事って学術系？企業だと使わないイメージ

**デフォルトの名無しさん** · 2015/11/05(木) 07:32:06.44

トヨタとか日産なんかの自動車の
研究開発で使ってるよ

mathworksの創始者もmatlabは最初
数学科とかには見向きもされず
工学系の学生に受けが良かったって
語ってなかったっけ？

**デフォルトの名無しさん** · 2015/11/06(金) 03:02:03.67

>>288-290
なるほど～
octave使えてればmatlabしなくても。。。ってことなのかな？
それともoctave使う層の人間はmatlabにない関数やらは自分で作っちゃうよね？ってこと？
そしたら最悪描画に限ったらgnuplotでもよくないかい？

**デフォルトの名無しさん** · 2015/11/06(金) 21:19:40.51

どっちかって言うと後者かな、自分で作っちゃうイメージ

気のせいかもしれないけど前にoctave触ったとき、メモリ馬鹿食いするプログラム運用するときや数千点プロットする時は
matlabの方が軽いかなーって印象があったような気がする
大きいデータをいっぺんに扱わにゃならんとか、マルチコアフル回転させたいとかならmatlabがいいのかも（octaveあんまり詳しくないです・・・）

**デフォルトの名無しさん** · 2015/11/07(土) 11:13:07.82

matlabにおいて既定の設定ですとdoubleで保存されます。
そのため変数が例えば a=2.63890000000... というふうに
小数点に無駄な0が多く含まれておりますmatファイルで保存を行う場合、容量が大きくなってしまいます。
出来れば無駄な0をさくじょしたいのですがどのように行えばよいでしょうか
よろしくお願いします

**デフォルトの名無しさん** · 2015/11/07(土) 13:14:30.83

ハローワークでも、自動車系で、
matlabの募集をたまに見る

でも徐々に、R, matlabの仕事は、
Pythonへ置き換わっていくかも

**デフォルトの名無しさん** · 2015/11/07(土) 15:07:03.07

仕事始めてから触りだした口だけど手計算のノリというか関数電卓の延長で使えるのが好きかな

**デフォルトの名無しさん** · 2015/11/07(土) 15:15:14.61

画像処理も最新研究のmファイルが落とせるニッチな世界

**デフォルトの名無しさん** · 2015/11/07(土) 16:07:24.48

後方互換性捨てるpythonなんて
怖くて仕事で使えねえよw

数年前のmatlabのparforはcやfortranの
openmpと比べて並列効果が小さかっ
たんだけど今はどうなんだろ？

**デフォルトの名無しさん** · 2015/11/08(日) 00:31:39.80

どの程度性能上がったかはしらんが単体だと4コアまでというクソ仕様

**デフォルトの名無しさん** · 2015/11/11(水) 13:34:26.77

乱数使ったジャンケンゲームどうやればいいんだ？

**デフォルトの名無しさん** · 2015/11/21(土) 15:33:21.26

rps関数とかかなー

**デフォルトの名無しさん** · 2015/11/24(火) 13:54:40.03

matlabでCAN通信を使いたいんだけどRTI CAN blocksetが見当たらないんだけどどうすればいいんですか？
dSPACEでファームウェアの更新でDS1103CANをアップデートしてdSPACEでは対応してるみたいなんですがちょっとこれからどうすればいいか分からないです

**デフォルトの名無しさん** · 2015/11/25(水) 17:51:08.97

>>286
遅レスなんで見てないかもだが
非線形関数を最適化する
fminsearch とか明らかに挙動が変わる
octaveは最適解探すの直ぐ諦める感じ

このレベルの関数自作できるのなら
octaveでも構わないんだろうけど

**デフォルトの名無しさん** · 2015/11/25(水) 23:48:39.89

lsqcurvefitで非線形曲線の近似曲線を求めようとしてるんですが、よい曲線形に収束しません。
どうすればいいのでしょう？

**デフォルトの名無しさん** · 2015/12/17(木) 17:16:26.68

pitchを増やす

**デフォルトの名無しさん** · 2015/12/29(火) 18:20:33.35

Matlab大好きな人がOctave は重すぎるって行ってけど
今はi686に最適化されてるから使ってる使っているCPUに最適化したら
そこまで重くなさそうなんだけどどうなんだろう

**デフォルトの名無しさん** · 2016/02/11(木) 00:49:42.98

質問失礼します。

※下付き文字を便宜上[]で示します。
　例えば水はH[2]O、二酸化炭素はCO[2]と表しています。

生物の個体数の変動を表したモデルとして、以下のようなロジスティック写像がある。
x[n+1]=ax[n](1-x[n])

x[n]は正規化されたn年目の生物の個体数、x[n+1]は次年の個体数。aはその生物の出生率として与えられる。
出生率aがもし3以下なら、その生物の個体数は一定に収束する。
しかし、aが3以上になると周期的な挙動を引き起こし、さらにaがある臨界点を超えると、非常に複雑な変動を観察することができる。

このロジスティック写像を視覚的に理解するため、MATLABを利用して以下に述べる様々なグラフィックを描させてみよ。

(1)時系列
横軸に時間(年次)n、縦軸に個体数x[n]をとり、ロジスティック写像の時間的な変動を観察せよ。
a=3とする。

(2)リターンマップ
横軸にx[n]、縦軸にx[n+1]をとり、ロジスティック写像におけるx[n]とx[n+1]の関係をプロットせよ。
a=3とする。

(3)分岐図
横軸に出生率a、縦軸にx[n]をとり。aの値によってx[1]、x[2]、…、x[n]がどんな値をとるか観察せよ。

よろしくお願いします。

**デフォルトの名無しさん** · 2016/02/11(木) 08:29:13.07

何も難しいことない。
指示されたとおりにグラフを描くだけ。

**デフォルトの名無しさん** · 2016/02/11(木) 15:45:36.50

spectrogram()関数について質問です
出力引数無しの場合スペクトログラムがプロットされますが、
出力引数［s,f,t,p］として、関数を使ったあと出力引数無しの場合と同じようなフィギュアをプロットしたいのですが何関数を使えばいいのでしょうか？

**デフォルトの名無しさん** · 2016/02/12(金) 13:48:38.84

以下のような非線形の3階常微分方程式がある。

dx/dt =-y-z
dy/dt =x+ay
dz/dt =bx-(c-x)z
(a=0.344 , b=0.4 , c=4.5)

これをレスラー方程式という。

このレスラー方程式をMATLABによる数値微分によって解き、x(t),y(t),z(t)の挙動を次の2通りの方法でコンピュータ画面上に描いてみよ。
なお、常微分方程式の数値積分法についてはルンゲッタ法を参考にするとよい。

(1)
x(t),y(t),z(t)の時間波形(横軸:時間(t),縦軸:x or y or zとした図)を描け。
また数値積分の時間の刻みΔtを何通りか変えて、Δtの値により波形が異なる様子を観察せよ。

(2)
x,y,zを3軸にとったこの方程式の解軌道の3次元プロットを作成し(1)と同様に解軌道に及ぼす影響について観察せよ。

正直何を言っているのかすらわからずお手上げ状態です。
よろしくお願いします。

**デフォルトの名無しさん** · 2016/02/12(金) 16:14:11.68

宿題やってもらいたいなら何でも質問スレにでも行けば？

**デフォルトの名無しさん** · 2016/02/25(木) 02:56:33.46

Neuralnetwork Toolboxを利用しようと思いサンプルを試していました。
しかし、パターン認識や時系列、クラスタリングと全てでサンプルデータを動かしてみることができませんでした。
どのエラーも、
「関数または変数 'network' が未定義です。」
と表示されます。
正直、サンプルを動かそうとしているのにエラーを吐かれて困っています。
VERコマンドでインストールされていることは確認しました。
R2013aとR2015aの両方で試しましたが同じエラーで実行できませんでした。
どうすればエラーを回避できるようになるでしょうか？

**デフォルトの名無しさん** · 2016/03/12(土) 22:51:46.82

PDEtoolboxの日本語解説サイトはここか？

**デフォルトの名無しさん** · 2016/03/13(日) 19:30:06.85

vec = [1 2 3]
この行ベクトルを指定した数だけ行の次元を増やして同じ値を入れるにはどうしたらいいですか？
2だったら
vec2 = [1 2 3;1 2 3]
3だったら
vec3 = [1 2 3;1 2 3;1 2 3]
みたいなのを
vec2 = ?(2,vec)
vec3 = ?(3,vec)
のように行いたいです

**デフォルトの名無しさん** · 2016/03/13(日) 21:52:47.20

>>315
まずrankを変えて行列を定義して、roopで順に入れ込んでいく。

**デフォルトの名無しさん** · 2016/03/13(日) 22:58:26.61

>>315
vec2 = repmat(vec,2,1)
vec3 = repmat(vec,3,1)

**デフォルトの名無しさん** · 2016/03/13(日) 23:11:42.92

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

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

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