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

**デフォルトの名無しさん** · 2012/09/03(月) 17:35:08.97

MATLAB　日本語版
http://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

前スレ
MATLABプログラミング　質問箱　その3
http://toro.2ch.net/test/read.cgi/tech/1254528108/

**デフォルトの名無しさん** · 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