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

1 名前：デフォルトの名無しさん [2006/11/10(金) 11:28:22 ] MATLAB　日本語版 www.cybernet.co.jp/matlab/ MATLAB 本家 www.mathworks.com/ Yahoo MATLAB dir.yahoo.co.jp/Science/Mathematics/Software/MATLAB/ MATLAB入門 feedback.mech.fukui-u.ac.jp/matlab/matlab0.htm MATLAB チュートリアル www-isys.eei.metro-u.ac.jp/~sho/pub/ttmatlab/Index.html 218 名前：デフォルトの名無しさん mailto:sage [2007/10/02(火) 15:25:07 ] 3次元プロットで、（例えばsurfなどで）値がaxisで指定した範囲外にあってもプロットされてしまうのですが、 axisで指定した範囲内のみプロットするようにするにはどうしたらよいのでしょうか？ 219 名前：デフォルトの名無しさん mailto:sage [2007/10/02(火) 16:06:06 ] 関数を定義する場合必ず別ファイルにしないといけないのですか？ 同一ファイル内で定義することはできないものでしょうか？ 220 名前：デフォルトの名無しさん mailto:sage [2007/10/02(火) 18:45:41 ] >>218 aってのを3次元の行列だとして surf(a(xmin:xmax,ymin:ymax,zmin:zmax)) じゃだめ？ 意味取り違えてたらごめんね。 >>219 pc11.2ch.net/test/read.cgi/tech/1163125702/82- 同じスレなので出来れば読んでね。 221 名前：デフォルトの名無しさん mailto:sage [2007/10/04(木) 03:19:13 ] matlabってリーもとデスクトップじゃつかえないのだね 222 名前：デフォルトの名無しさん [2007/10/11(木) 14:00:14 ] bata版のNAG ToolBox for MATLAB 無料評価版 www.nag-j.co.jp/ToolBox.htm 223 名前：デフォルトの名無しさん mailto:sage [2007/10/11(木) 14:23:25 ] pythonで大部分の数値計算できるけどmatlabってやっぱりまだ需要あるの？ 224 名前：デフォルトの名無しさん mailto:sage [2007/10/11(木) 16:29:25 ] >>223 C で python のできることは全部できるけど python って需要あんの？ 225 名前：デフォルトの名無しさん mailto:sage [2007/10/12(金) 00:19:45 ] 俺紙とエンピツで全部やってるけどコンピュータって需要(ry 226 名前：デフォルトの名無しさん mailto:sage [2007/10/12(金) 05:57:48 ] さすがにMATLABのスレでそりゃないだろ。 手作業でやったら一生かかっても終わらん計算してるだろうに。 227 名前：デフォルトの名無しさん [2007/10/23(火) 14:25:52 ] 質問なのですが、下記のプログラムを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 228 名前：デフォルトの名無しさん mailto:sage [2007/10/23(火) 18:32:33 ] それは質問じゃなく変換依頼 229 名前：デフォルトの名無しさん mailto:sage [2007/10/25(木) 03:43:18 ] タブ区切りのデータなのですが、一行目が見出しになっています。 一行目をのぞいて簡単に読み込む方法はないでしょうか？ 230 名前：デフォルトの名無しさん mailto:sage [2007/10/25(木) 07:53:38 ] fgetl 231 名前：デフォルトの名無しさん mailto:sage [2007/10/25(木) 16:19:50 ] たぶん激しく初心者な質問でスマソなんですが、サイズがL,M,Nな3次元変数を3次元プロットする方法知りませんか？ 232 名前：デフォルトの名無しさん mailto:sage [2007/10/25(木) 22:13:34 ] 強度イメージから一定範囲の強度値だけを抽出する 方法を教えてください 233 名前：デフォルトの名無しさん [2007/10/25(木) 22:39:50 ] 田中真紀子のﾎﾟﾛﾘ画像ありますか? 234 名前：デフォルトの名無しさん [2007/10/26(金) 09:37:43 ] オナニーをすると身長が伸びないというのは本当ですか？ 235 名前：231 mailto:sage [2007/10/28(日) 16:40:35 ] >>231 自己解決しますた 236 名前：デフォルトの名無しさん [2007/11/01(木) 13:39:31 ] 質問させてください。 specgramの逆の操作がしたいです。 X=specgram(x,128,8000); で、 スペクトログラムXから音xを合成したいのですが、どうしたらいいですか？？ 237 名前：デフォルトの名無しさん mailto:sage [2007/11/02(金) 00:57:30 ] synthesizeとかいう命令があったような気がする。 matlabじゃなくてoctaveの命令だったかも。 238 名前：デフォルトの名無しさん [2007/11/02(金) 06:26:24 ] Core2duoに変えたんだけど、６４BITに変えたら 計算すごく速くなる？ ベンチとかないかなぁ 239 名前：デフォルトの名無しさん [2007/11/02(金) 19:57:08 ] 入力したサウンドデータを，何の音なのか (例えば，パトカーのサイレン，人の叫び声等) 分類するソフトを個人的に作成しています． とりあえずmatlabでスペクトル表示を行い， 様々なサウンドデータのスペクトルを見たのですが， 分類の条件を何にしようか考えています． 「パトカーのサイレンは○○Hzで特に高いピーク値を得られる」 などの特徴を分類する際の条件にできるかなと思ったのですが， この他に，音を分類する上で有効な条件となるようなものを ご存知の方がいましたらご教授お願いします． 240 名前：デフォルトの名無しさん mailto:sage [2007/11/02(金) 20:58:50 ] スレ違い 241 名前：デフォルトの名無しさん mailto:sage [2007/11/02(金) 21:45:16 ] 自分の目でいういろ比較して確かめなさい 242 名前：デフォルトの名無しさん [2007/11/04(日) 12:40:22 ] はじめまして、matlab初心者です。 早速ですが質問させて下さい。 任意形状の２次元物体をメッシュ切りした上に、 コンターグラフの様な物を描かせたいのですが、 どなたか出来る方法をご存じないでしょうか？ 宜しくお願いします。 243 名前：デフォルトの名無しさん [2007/11/04(日) 19:30:34 ] matlabのメインウインドに打ち込むだけでsimulinkの開始ボタンを押したのと同じ処理になるコマンドってありますか？ 244 名前：デフォルトの名無しさん [2007/11/05(月) 20:40:09 ] >>243 SIM 245 名前：デフォルトの名無しさん [2007/11/07(水) 10:02:25 ] 例えば４ｘ４のマトリックスを２ｘ２にして、各要素は元々のヨウ素を加算してくれる関数はないですか？ ループとsum関数を使えばすぐに出来るんですがね 246 名前：デフォルトの名無しさん [2007/11/07(水) 10:05:37 ] >>211 Image Processing Tool Boxだったかにあるよ。座標変換してからinterp2に流し込めばすぐにできるけどね。 247 名前：デフォルトの名無しさん mailto:sage [2007/11/07(水) 22:42:27 ] >>245 じゃぁ、ループとsum関数を使えば？ 248 名前：デフォルトの名無しさん mailto:sage [2007/11/08(木) 01:47:43 ] >>247 matlab ではできるだけループを使わないのが重要 速度がダンチ 249 名前：デフォルトの名無しさん mailto:sage [2007/11/08(木) 05:18:32 ] 書き方が曖昧でわからんがエスパーの俺ならこうかな。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); 250 名前：デフォルトの名無しさん mailto:sage [2007/11/08(木) 05:19:06 ] :と,一カ所ミスった。 251 名前：デフォルトの名無しさん [2007/11/09(金) 15:59:58 ] ??? エラー: ==> anfisedit SCRIPT anfisを関数として実行してみてください. これを解決するにはどうすればよい　のでしょうか？ 一週間前ぐらいにMATLAB始めた者です。 252 名前：デフォルトの名無しさん mailto:sage [2007/11/09(金) 20:08:08 ] >>251 SCRIPT anfisを関数として実行してみてください. 253 名前：デフォルトの名無しさん mailto:sage [2007/11/11(日) 22:32:56 ] MATLABで式をニュートンラプソン法でとくにはどうしたらいいでしょうか？ 分からなさすぎてノイローゼになりそうです。 254 名前：デフォルトの名無しさん mailto:sage [2007/11/12(月) 02:23:50 ] まぁスレ違いだわな 255 名前：デフォルトの名無しさん [2007/11/12(月) 07:39:22 ] GUI上でxの変域をかえたいんですけどどうすれば良いですか？ 256 名前：デフォルトの名無しさん mailto:sage [2007/11/13(火) 21:23:48 ] GUI上で数値を入力させて、 set(h, 'xlim', [x1 x2]); じゃダメなの？ 257 名前：デフォルトの名無しさん [2007/11/16(金) 21:18:38 ] FORTRANで書いたものを、 MATLABに変換したいのですが、手書き以外で 何かいい方法あったら教えてください。 258 名前：デフォルトの名無しさん [2007/11/16(金) 22:16:21 ] fftでパルスの周波数帯域を表示させたいんですが、 周波数軸をどうやって設定すればよいかわからないのです。 パルスはガウシアンパルスでUWBなので５ギガ位まで表示させるつもり何ですが、うまくいかないのです。 どなたかお願いします。 259 名前：デフォルトの名無しさん mailto:sage [2007/11/17(土) 02:50:17 ] >>257 グーグルの使い方くらい知らないと 今後いろいろ苦労すると思うから覚えた方がいいよ。 260 名前：デフォルトの名無しさん mailto:sage [2007/11/17(土) 18:16:24 ] ファイルの読み込みで指定した列だけ読む方法ってありましたっけ。 A = load("hoge.txt"); B = A(:,4); みたいに、全部読んで目的列以外は捨てればいいんですが、 ファイルのサイズが大きいと、 全部読むのに時間がかかるので、効率化したいのです。 261 名前：１８ [2007/11/17(土) 23:13:55 ] matlabでシミュレーションを５回もするとメモリが５００Ｍから１Ｇくらいまで増えてしまうんですが、 これってなぜなんでしょうか？ 262 名前：デフォルトの名無しさん mailto:sage [2007/11/18(日) 00:14:17 ] >>260 dlmreadのマニュアル読むと良いよ。 263 名前：デフォルトの名無しさん mailto:sage [2007/11/18(日) 13:49:03 ] 対数データのフィッティングの方法を教えてください。 通常のデータであれば、最小二乗法ですが、対数データの場合はどのようにやればよいのでしょうか？ 264 名前：デフォルトの名無しさん mailto:sage [2007/11/18(日) 16:04:56 ] 対数の底で冪をとって最小二乗法 265 名前：beta-infinite [2007/11/20(火) 13:16:10 ] 「MATLABの極」というブログを書いています。良ければ、ご参考に beta-infinite.blogspot.com/ 266 名前：デフォルトの名無しさん [2007/11/25(日) 19:51:53 ] 著作権違反ですよｗ 267 名前：デフォルトの名無しさん [2007/11/26(月) 18:14:01 ] MATLABを使ってＢＰ法のニューラルネットワーク作らないといけないんですが 誰か請負いませんか？ 価格は応相談、ＭＡＴＬＡＢ使えてＢＰ法分かる人ならそんなに難しい作業じゃないと思います 仕様とかその辺は請け負ってくれる人がいたらメールなりなんなりで連絡します 268 名前：デフォルトの名無しさん [2007/11/26(月) 18:32:59 ] 今、手一杯だからなー。 月８０万だったら考えるけどね。 それかNNツールキット買えばいいんじゃないの？ 269 名前：デフォルトの名無しさん mailto:sage [2007/11/26(月) 18:46:20 ] ていうか自分で作れって感じ 層やニューロン数を完全に一般化した物でなければ速攻で作れるでしょ もちろんfor文抜きでな 270 名前：267 [2007/11/26(月) 19:12:50 ] 料金は ３時間未満：一律12000円 ３時間以上：時給4000円 ７時間以上：一律28000円 このあたりでよろ おれ自身じゃ作れないので金で解決する 271 名前：デフォルトの名無しさん mailto:sage [2007/11/26(月) 19:19:22 ] 安すぎ 272 名前：デフォルトの名無しさん mailto:sage [2007/11/26(月) 19:29:20 ] 委託研究でMATLABでソフト作ったときは、３０万もらったな 273 名前：デフォルトの名無しさん mailto:sage [2007/11/26(月) 20:31:27 ] 作成物の複雑さとか教えてもらえずに 7時間以上28000円は安いな。 274 名前：267 [2007/11/27(火) 00:00:36 ] ここでは話せないが詳しい仕様を聞いてくれるやつはレスくれ 秘密厳守の上その時点で断ってくれてもＯＫだから 275 名前：デフォルトの名無しさん mailto:sage [2007/11/27(火) 00:06:40 ] こんなところで依頼しようとしている時点でダーティなイメージがついてて嫌だ 全然予算がないようだし学生の卒論か？ 276 名前：267 [2007/11/27(火) 01:14:15 ] >>275 まさにそれ 誰か頼むよ〜 277 名前：デフォルトの名無しさん [2007/11/27(火) 01:19:35 ] 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と書いてあったのでここに質問しました＞＜ スレ違いでしたらすみません。。 278 名前：デフォルトの名無しさん mailto:sage [2007/11/27(火) 06:40:06 ] 別に matlab じゃないとできないわけじゃない。プログラミング言語に依存してたまるか 279 名前：請負人 mailto:i100lao@yahoo.co.jp [2007/11/27(火) 11:53:45 ] 請負内容を見て判断したいです。詳細内容などをメールをください 280 名前：請負人 mailto:i100lao@yahoo.co.jp [2007/11/27(火) 11:59:51 ] >>267 >>270 請負内容を見て判断したいです。詳細内容などをメールをください。 秘密厳守を守ります。 281 名前：デフォルトの名無しさん mailto:sage [2007/11/27(火) 18:59:28 ] 手元にはバックプロパゲーション法の自作ソースくらいあるんだが バカ学生に使いこなせるとは思えん 282 名前：デフォルトの名無しさん mailto:sage [2007/11/28(水) 15:25:32 ] matlab にファイル名を変えるコマンドってありましたっけ？ 同時にファイル内の関数名もあわせて変えてくれるとベスト 283 名前：デフォルトの名無しさん mailto:sage [2007/11/29(木) 23:23:21 ] ttp://www.cybernet.co.jp/matlab/support/techkwdb/search_view.php?id=1303 このページの h=plot(0,0,'.b','MarkerSize',30,'EraseMode','xor'); を使ってアニメーションを作る方法で、 同時に2つの色の違うボールを描画するにはどうすればいいのですか？ 284 名前：デフォルトの名無しさん [2007/11/29(木) 23:34:51 ] matlabで書いたグラフで、線と線に囲まれた部分の 面積を求めることは可能でしょうか。 ファイルをアップしましたので、よろしかったら見てみてください。 www.csync.net/service/file/index.html (ファイル名はmatlabです) 画像の灰色の部分の面積を求めたいのですが、 積分のやり方もわからないのです(:_:) あとMファイル（txtになってしまってます）のコピーが二つありますが、 mensekiの方は画像の元になるグラフが書けます。 danmenの方は、元のプログラムですが、 xの値を０〜０．５まで変えると赤と緑の線が一周するように できているのですが、この二つの線と、黄色の線が作る面積を求めたいんです。 どんな方法でもいいのでどなたか教えてくださいm(__)m 285 名前：beta-infinite [2007/11/30(金) 00:29:49 ] >>283 「MATLABの極」にソースをアップロードしました。いかがでしょうか beta-infinite.blogspot.com/ 286 名前：デフォルトの名無しさん [2007/11/30(金) 03:26:05 ] 教えてください。 MAT LABで書いたFigureのX軸Y軸をX→Y、Y→Xと逆にしたいのですが どうすればいいんでしょうか？お願いします。 287 名前：283 mailto:sage [2007/11/30(金) 10:53:34 ] >>285 ありがとうございます。 今までhold onとか使っていたんですが、こんなシンプルにいけたんですね。 ですが、シミュレーションしたものをアニメーションとして表示したいので この方法だとどうしても描画が遅くてシミュレーションの結果通りに描画できません。 そこで>>283の方法は高速だということなので、 >>283のイレースモードxorを使う方法でのやり方を知りたいのです。 >>283の方法以外でも高速に描画できる方法をご存知の方はお願いします。 具体的にしたいことは、 ・2つボールの0.001秒毎の位置が2つの配列に入っている。 ・それをある時間dt秒毎に描画したい。 ・1回の描画毎にpause(0.001*dt);などして時間に正確にアニメーション化したい。 自分のやり方や>>285さんのやり方だと dtを大きくすればするほどアニメーションがゆっくりになってしまいます 288 名前：283 mailto:sage [2007/11/30(金) 10:56:39 ] >>287の修正 ×　pause(0.001*dt); ○　pause(dt); 289 名前：MATLABの極 [2007/11/30(金) 12:02:44 ] >>288 PAUSE(n) は、n秒間停止して、その後実行を続けます 当然ですが、pause(dt)のdtを大きくすれば、 アニメーションがゆっくりになってしまいます 290 名前：MATLABの極 [2007/11/30(金) 13:29:41 ] >>286 view(90,90)という関数を使えば、いいと思います。 beta-infinite.blogspot.com/2007/11/matlabfigurexyxyyx.html に参考 291 名前：283 mailto:sage [2007/11/30(金) 13:44:06 ] >>289 ですが、dtを大きくするとその分荒いサンプリングをすることになるので 描画時間がdt秒に比べて十分小さければ同じようなアニメになるはずですよね コマ送りになるかどうかの違いはあれど。 でもやはり描画時間はそこまで小さくするのは無理なのでしょうか。 MATLABではこんなことはできないのでしょうか 292 名前：MATLABの極 [2007/11/30(金) 14:04:16 ] >>291 現状の１コマの処理時間は、 = matlabの演算時間 + 描画の時間 + pause(dt)の時間 matlabの演算時間 + 描画の時間をなくする方法は、getframe,movie,movie2avi関数で 使うといいと思います。 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 に参考 293 名前：デフォルトの名無しさん [2007/11/30(金) 15:21:36 ] beta-infinite.blogspot.comってパクリしかないじゃんかｗ 294 名前：MATLABの極 [2007/11/30(金) 15:45:10 ] >>293 MATLABの修業中です。覚え書きのブログだと思っていただければ、 295 名前：283 mailto:sage [2007/11/30(金) 16:30:08 ] >>292 やっぱgetframe,movie,movieしかないですかね ありがとうございました 296 名前：デフォルトの名無しさん [2007/11/30(金) 17:13:21 ] 著作権違反だので、通報しました 297 名前：デフォルトの名無しさん mailto:sage [2007/11/30(金) 17:57:24 ] FFTを使って、スペクトルを計算するのって、どうすればいいの？ 298 名前：デフォルトの名無しさん mailto:sage [2007/11/30(金) 21:40:02 ] だので 299 名前：デフォルトの名無しさん mailto:sage [2007/11/30(金) 22:31:07 ] >>297 help fftとコマンドを打って出てきた説明を理解できなければ 回線を切って死ぬか教授に頭を下げて教えを請いなさい 300 名前：デフォルトの名無しさん mailto:sage [2007/12/01(土) 04:36:12 ] >>297 beta-infinite.blogspot.comだけは見ないように。 301 名前：デフォルトの名無しさん mailto:sage [2007/12/01(土) 12:59:02 ] fftはtwo-sideの出力を出すので、help fftだけでは話は完結しない 302 名前：デフォルトの名無しさん [2007/12/01(土) 16:02:42 ] 教えてください。 プログラムで図を描かせたときにX,Y,Zラベルの刻みのラベルを 変更したいのですが、関数などプログラムで書かせることができるでしょうか？ 少し伝わりにくい文章ですが、よろしくお願いします。 303 名前：デフォルトの名無しさん mailto:sage [2007/12/01(土) 16:26:01 ] なんか楽に関数名リネームする方法ねーのかよ 304 名前：デフォルトの名無しさん mailto:sage [2007/12/01(土) 17:18:44 ] >>303 関数名リネームって>>282が言ってるようなことか? もしそうならファイル名と違ったらファイル名の方が優先なだけだから コード中の関数名は適当でOK。 しかしそういう時は普通呼ぶ側のコードも変えないといけないから どっちにしろ一括置換とかする程度の手間をはぶくのは無理なわけで コードの頭の関数名変えるくらいの手間は取るに足らないと思うが。 305 名前：デフォルトの名無しさん mailto:sage [2007/12/01(土) 19:52:21 ] >>297 ここ見るといいよ www.mathworks.com/support/tech-notes/1700/1702.shtml 306 名前：デフォルトの名無しさん [2007/12/02(日) 18:39:33 ] >>302 できる。 307 名前：デフォルトの名無しさん mailto:sage [2007/12/04(火) 06:02:50 ] 関数の入力にベルトルの集合としての行列をいれるとき皆は NxD 派？ DxN 派？where N is the number of vectors and D is the dimension of each vector. どっちのほうが親切なんだろう？toolbox によって（人によって）違うんだよねぇ。 308 名前：デフォルトの名無しさん mailto:sage [2007/12/04(火) 09:55:13 ] シミュレート板で回答を得られなかったので。 simulinkである変数がある値になったときシミュレーションを終了したいのですが どうすればできますか？ 309 名前：デフォルトの名無しさん mailto:sage [2007/12/05(水) 06:53:25 ] ムービーを作成する際に、figureウインドウを出さずに作成することはできないでしょうか？ figureウインドウを作成してgetframeで取り込む場合、スクリーンキャプチャをしているようで、 他の作業をすると、その画像がムービーの中に入ってしまいます。 そのためムービー作成中は何も作業ができません。 310 名前：デフォルトの名無しさん [2007/12/05(水) 16:25:53 ] DICOMファイルから画像読み込み、JPG形式で書き込みしようとしてるのですが、 X=dicomread('--.dcm'); [Y,map] = gray2ind(X); imwrite(Y,map,'--.jpg','jpg'); とすると、colormapを使用したimshowの表示画像と異なってしまいます。 よろしくお願いします。 311 名前：デフォルトの名無しさん [2007/12/05(水) 18:06:21 ] 正の実数からなる対角行列の、各行と列について、 要素の大きさが対角成分の要素の大きさの0.5から2倍の範囲なら１で置換、 それ以外なら０で置換する、 ってforループを使わずにできるでしょうか？ 312 名前：デフォルトの名無しさん mailto:sage [2007/12/05(水) 23:19:50 ] 質問ばっかで誰も答えてくれない罠ｗｗｗｗｗｗｗｗｗ 313 名前：デフォルトの名無しさん mailto:sage [2007/12/06(木) 00:59:41 ] >>311 対角行列なら対角成分以外すでに全部０やん 314 名前：デフォルトの名無しさん mailto:sage [2007/12/06(木) 02:21:48 ] これぞまさに質問スレ 315 名前：デフォルトの名無しさん mailto:sage [2007/12/06(木) 08:48:34 ] vistaは2007aで正式サポートですが、2006aで動かした場合、どこかに既知の問題のようなリストはないでしょうか？ 316 名前：デフォルトの名無しさん mailto:sage [2007/12/07(金) 01:04:38 ] >>307 僕はMatlabで音系のことをやってるけど、 wavwriteとかのデフォルトが行＝サンプル／列＝チャンネルなので、 ベクトルの集合はいつも縦長ベクトルを横にいくつか並べるというかんじ。 DxNってことになるのかな？ 317 名前：デフォルトの名無しさん mailto:sage [2007/12/07(金) 12:48:35 ] 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」 と、こんな感じなのですが,どなたかサクっと作ってくださいませんでしょうか・・・ 318 名前：デフォルトの名無しさん mailto:sage [2007/12/12(水) 05:19:53 ] 秀丸用の強調定義ファイルはありませんか？

---

---

[ 続きを読む ] / [ 携帯版 ]

read.cgi ver5.27 [feat.BBS2 +1.6] / e.0.2 (02/09/03) / eucaly.net products.
担当:undef