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/

**デフォルトの名無しさん** · 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
ありがとうございます

**デフォルトの名無しさん** · 2016/04/30(土) 02:50:38.37

Surfの平面内グリッドを消せないものでしょうか

**デフォルトの名無しさん** · 2016/05/02(月) 10:04:30.38

>>319
surf(peaks, 'LineStyle', 'none')
ってこと？

**デフォルトの名無しさん** · 2016/05/04(水) 02:13:43.76

休み明けに試してみます
ありがとうございます

**デフォルトの名無しさん** · 2016/05/13(金) 04:38:18.45

マトラボのダウンロード番って一つのOSしか選べないのん？
当方、Macintosh、Linuxマシンで運用する予定なんやけど

**デフォルトの名無しさん** · 2016/05/13(金) 09:14:23.29

>>322
わいはいつもWin/Mac/Linの最新版ISOファイル落としてきとるがの

**デフォルトの名無しさん** · 2016/05/13(金) 12:37:40.55

>>323
えっ、それって割れちゃうやろ？
正規版でやで？

**デフォルトの名無しさん** · 2016/05/13(金) 20:37:42.06

matlabのhome版のDL購入(？)のクレジット決済ってVプリカでも可能ですか？

**デフォルトの名無しさん** · 2016/05/14(土) 18:46:12.74

>>325
325自身ですが、使えなかったらamazonで使えばいいかって感じで買ったら普通に買えました

**323** · 2016/05/19(木) 16:35:18.19

>>324
複数マシンにインストールしたいとき用に、ISOが公式に用意されてんの。
毎回ダウンロードだと帯域を食うでしょ。

**デフォルトの名無しさん** · 2016/05/23(月) 11:14:31.10

基本的な質問で申し訳ないんですが、
配列A=zeros(5,2)みたいな配列があったとして、
2から4行目をB=[1,1]にするうまい方法ってありませんか？

スカラーの場合みたいにA（2:4,:）=Bのように使える関数みたいなものはありませんか？
A（2:4,:）=repmat(B,3,1) とかやるしかないんでしょうか？

**デフォルトの名無しさん** · 2016/05/23(月) 11:16:11.12

>>328
B=[1,1]はただの例で他の適当な1×2の配列(B[2,3]等)でもできる方法が知りたいです。

**デフォルトの名無しさん** · 2016/05/30(月) 19:23:25.28

F(ｚ,w)を定義して、fsolveでF＝０を解きたいです
ここで最初wは初期値を与えて、zについてまず値を求めます
次に２回目以降はw=z、前回のzの値に更新してふたたびzについてwを更新したFを解きたいです
for文で10回繰り返すとして、wの値をfor文の中で更新しても関数を呼びたすたびに初期値のwになってしまいます
どのようにすればwおよびFを更新できますか？
ちなみに関数定義の中にwの値を与えず、実際のワークスペースでw=...と与えても入力引数が不足してますとエラーになってしまいます

**デフォルトの名無しさん** · 2016/05/30(月) 19:35:57.23

wをローカル変数に汁

**デフォルトの名無しさん** · 2016/05/30(月) 20:36:10.69

>>331
特に変数をglobalで宣言してないのですが、どのように変更するのでしょうか

**デフォルトの名無しさん** · 2016/06/04(土) 15:53:10.19

matlab homeだと出力したファイルに電子透かしっていうのがついてるらしいんだけど、これはどんな感じなの？

**デフォルトの名無しさん** · 2016/06/04(土) 19:01:28.16

waterprint

**デフォルトの名無しさん** · 2016/06/18(土) 03:34:19.81

contorf(x,y,z,20) すると等高線が20本引けるかと思うですが、等高線を消す方法ってありますか？

**323** · 2016/06/18(土) 09:21:22.78

>>335
contourf(x, y, z, 10, 'LineStyle', 'none')

**デフォルトの名無しさん** · 2016/06/23(木) 23:07:27.67

simulinkで、
2次元（x(t),y(t)）のデータをグラフに表示させて、
時間経過(t)とともに変化させたいのですが、
そういうブロックあるんでしょうか？

**デフォルトの名無しさん** · 2016/06/23(木) 23:38:11.33

先生に聞いた方がはやいよ

**デフォルトの名無しさん** · 2016/06/25(土) 15:22:53.64

figureからパラコマの動画を作ってどのように推移していくか
見ようとしています。そこで
figure
for i=1:360;
scatter(x(:,:,i),y(:,:,i),1,z(:,:,i),'filled');
hold on
.....
end
みたいな感じで作ってみたのですが
操作自体はあってるのですが
データ量がすごく多く大変時間がかかるようでした…。
元のデータは640×480というデータ点数に
それぞれxとyとzの座標データ(mm）がありまして、それが360枚分ございます。
scatterじゃなくてこっちのほうが軽い、や
for文でやるべきではないなど全くわからない状態ですので
どなたかご教授していただけるとありがたいです。よろしくお願いします。

**デフォルトの名無しさん** · 2016/06/27(月) 02:36:51.03

自己解決できました。

**デフォルトの名無しさん** · 2016/07/23(土) 23:23:08.11

fsolveを使って２次の非線形連立方程式を解こうとした所、「Equation solved at initial point.」
という結果が出力されました。
最初は本当に初期値がそのまま解になる事もあるかと思い色々と初期値を変えて再計算をした所、
あらゆる初期値で同じ結果が出力されるという結果になってしまい困っております。
この「Equation solved at initial point.」というのは方程式にどのような問題がある時に発生するものなのでしょうか？
私は初め方程式の値が小さすぎるために初期値の段階で０に収束してしまいソルバーが解に到達したと判断している…と考えたのですが、
これだとべらぼうに大きい初期値を用いた場合にも同じ結果になる事との説明が合わないですよね。
（ただ、方程式中の変数は三角関数の中に入っているので初期値の値に関わらず式の形的に０に収束しているという事もあり得るとも考えています。
すみません、この辺は私の理解不足で支離滅裂です）。
また、残差は仮に方程式をf(x,y)+a=0 とするとaの値がほぼそのまま出力され、アルゴリズムは信頼領域法を用いています。

**デフォルトの名無しさん** · 2016/07/24(日) 09:11:21.10

長いな。

**デフォルトの名無しさん** · 2016/07/24(日) 11:31:35.64

>あらゆる初期値で同じ結果

わろた

**デフォルトの名無しさん** · 2016/07/24(日) 11:33:54.85

>>341
そこまで判ってて方程式を晒そうとは思わないのかな

**デフォルトの名無しさん** · 2016/07/24(日) 11:55:17.48

>>341
大きいとか小さいとかじゃなくて最初から微分値が0なんじゃねーの
小学校からやり直せ

**デフォルトの名無しさん** · 2016/09/11(日) 17:45:42.13

メインの.mファイルとクラス定義を1ディレクトリにまとめて実行することってできる？
他人にプログラム使ってもらうときに、検索pathがどうのこうのと説明するのが面倒なんだけど…

**デフォルトの名無しさん** · 2016/09/12(月) 16:21:29.48

そういう香具師と付き合っても疲れるだけだよ

**デフォルトの名無しさん** · 2016/09/13(火) 12:45:54.10

イコライザが動かん

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

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