mathematicaによるシミュレーション教えてください
■ このスレッドは過去ログ倉庫に格納されています
0001名無しさん＠１周年
NGNG求めたのですが，この関数の例えば0.1秒刻みの値を
出力させるにはどうしたらよいのでしょうか？
できれば，Exelに出力させたいのですが。
ご存知の方，教えていただけませんか。
0153名無しさん＠１周年
NGNGわかりやすいのがあんまり無いな。
http://florida.mes.titech.ac.jp/ga-j.html
0154名無しさん＠１周年
NGNG0155152
NGNGありがとうございました。
さっそく勉強してみます。
0156あぼーん
NGNG0157名無しさん＠１周年
NGNG0158名無しさん＠１周年
NGNGこれができないものかと悩んでいます。。。
0159名無しさん＠１周年
NGNGじゃだめ？
0160158
NGNG一般にｎ重となるとどうしていいかわからなくて。。。。
∫…∫ｄｘ１ｄｘ２・・・・ｄｘｎ
っていうのをmathmaticaを使って計算させてみたいんです。
0161名無しさん＠１周年
NGNG具体的にどういう計算をしたいか良くわからんが、
Ver4.1以降ならBoole関数使うと上手く行くかも。
0162名無しさん＠１周年
NGNGつまり、引数が可変ということ？
だとすると、こんなのはどうでしょう（頭で考えただけで未検証）
１）積分変数とその区間の組を List で作る
２）そのリストに被積分関数を Prepend
３）Apply で List を Integrate
0163162
NGNGすいません。
大人しくリーマン和かモンテカルロを構成するのかなぁ。。。
0164名無しさん＠１周年
NGNG0165名無しさん＠１周年
NGNG時間分割タイプの式をちらっとしか見たことありませんが
ちょっとしたソフトウエアプロジェクト並になると思います。
ググったら、２～３の米国大学でシラバスのトピックで
見かけた内容ですね。
領域が変なタイプの場合の対応として
積分領域を含む単純タイプの新領域 +
領域内外区別用の Heviside x 被積分関数
にしてみたり
有限分割で式を組み立てて、 Limit もしくは
分割数大で数値的な収束を見る
C とかでライブラリがありそうな気がしますので、
そのソースを解読し、 Mathematica ならではの
リスト処理をするのが早いと思います。
0167名無しさん＠１周年
NGNG0168名無しさん＠１周年
NGNG0169名無しさん＠１周年
NGNG多分、任意のn個のdataを自動的にファイルの名前にaaa_###.txt
としてセーブしたいんですよね？
めちゃめんどうです。
0170169
NGNGに比べてかなり面倒です。
data[[n]]としてわりあてて
Table[Export[data[[n]].txt,data[[n]],"Table"],{n,1,100}]
として発生しておき、実行するしかないでしょう。
やったことないのでちゃんと動くか分かりませんが…
0171名無しさん＠３周年
NGNG各層3層でのさくっと書いてくださる方，お願いします．
0172BP
NGNGバックプロパゲーションならネットでも簡単にひっかかるだろ
検索してみな
0173名無しさん＠３周年
NGNGhttp://www.catkk.co.jp/mathcd/index.html
これ使うと勉強効率あがるかしら？
0174名無しさん＠３周年
NGNG詳しく書かれてるＨＰとかありませんかね？
wolframのＨＰだけだと、よくわかりません、、、
0175名無しさん＠３周年
NGNG0176200
NGNGついでにage
0178超初心者
NGNGImplicitPlotはどうやったら使えるようになるのでしょうか?
ヘルプにある以下のプログラムを打ち込んでも
下の行の式をそのまま返してくるのですが･･･
<< Graphics`ImplicitPlot`
ImplicitPlot[x^2 + 2 y^2 == 3, {x, -2, 2}]
きちんと各行shift+enterで入力決定しております。
よろしくお願いします。
0179名無しさん＠３周年
NGNGPackageがうまく組み込めていないかもと思われ。
Input menu->Get File Pathで、<<のあとにパスを入れてミソ
きれいな楕円だね？
0180超初心者
NGNGお返事ありがとうございます。
申し訳ありませんが、具体的なプログラムを書いていただけますでしょうか？
<<C:\Program Files\Wolfram Research\Mathematica\5.0\SystemFiles\Graphics\Packages\FileGraphics.m Graphics`ImplicitPlot`
ImplicitPlot[x^2 + 2 y^2 == 3, {x, -2, 2}]
ではダメでした。
0181名無しさん＠３周年
NGNGResearch\\Mathematica\\3.0\\AddOns\\StandardPackages\\Graphics\\ImplicitPlot.m"
で、shift+enter
ImplicitPlot[x^2 + 2 y^2 == 3, {x, -2, 2}]
で、また、shift+enter
0183名無しさん＠３周年
NGNG解きたい変数の数が何十個とあります。
仕方ないので
vars = {v1,v2,…全部の変数を入力してる…,v30}
Solve[ eqns, vars ]
とやってるんですが、これではあんまり大変なので
なにか上手い方法はないでしょうか。
0184名無しさん＠３周年
NGNGvar=Table[v[[i]],{i,1,30,1}];
にした少しは楽になるのでは？
>>182
WindowsでMathematica使うとよくそんなことあるの？（パッケージがうまくよめないとか）
Mac使ってますが、一度もそんなことないなぁ
0185名無しさん＠３周年
NGNGSolveじゃなくて行列使って計算するのじゃだめなのかな？
0186名無しさん＠３周年
NGNG0187名無しさん＠３周年
NGNG横レスですが、Mathematicaは操作の途中で生成した、
オブジェクト変数名を起動時からずっと記憶しますよね。
途中で重複した変数名が残ってしまって悪さをしたのでは？
それら途中履歴をクリアする、Clear文っていうのもあります。
>186
Mahtematica関係の本は入門書からとんでもない高度な話題のまで
いろいろありますが、自分のオススメは、
牧野書店 Mathematica 数値数式プログラミング
ISBN4-7952-0113-7 C3055
0188名無しさん＠３周年
NGNG何が原因なんでしょ。
0189名無しさん＠３周年
NGNG1+1を１＋１っていれているからじゃないか？
0191名無しさん＠３周年
NGNGたとえば ケインジアン交差図と投資関数から導かれるIS曲線
Y=C(Y-T)+I(r)+G = a+b(Y-T) + c-dr + G
⇔Y=(1/1-b){a+c+G-bT-dr}
流動性選好理論から導かれるLM曲線
Ld(r.Y)=eY-fr≡Ld(M/P)
⇔r=(e/f) - (1/f)(M/P)
これらを (r.p.y)空間上で交差させて
この交点から導かれる総需要曲線
Y=(z/1-b)[a+c+G-bT-{d/f-fb)}(M/P)}]
z=f/[f+de/(1-b)] をその接点を切り取る形で表示して
これに 政府がマネーサプライを増加させたときや
租税を増加させたときどのような変化がおこるかわかりやすく図示
またピグー効果が起こった時 C(Y-T,M/P) = a+b(Y-T) +ν(M/P)
と消費関数がなるわけですが
このとき 総需要関数がどのように変化するかもわかりやすく
示したいんです
でもmatematica で 数理的に説明するよりも
自分で勝手に三次元グラフを作ってしまったほうが手っ取り早いんですね
Y=(z/1-b)[a+c+G-bT-{d/f-fb)}(M/P)}] この式から見て Mの増加や
C(Y-{T+kY}) = a+b(Y-{T+kY}) となる累進課税時に 曲線の傾きがどうなるかは簡単にわかりますが
実際にグラフで視覚的に体験して見たいのです
Mの増加は LM曲線を左方にシフトさせますが
それが需要曲線の変動とどうつながるか
三次元グラフを使って立体的に示したいんです
新着レス 2004/06/25(金) 21:57
88 名前： デフォルトの名無しさん 投稿日： 04/06/25 21:57
例えば コブダグラス型の生産関数です
z = 60 x^0.6 y^0.4
g = Plot3D[z, {x, 0, 10}, {y, 0, 10}, ViewPoint -> {-1, -3, 0.7}]
ContourPlot[z, {x, 0, 10}, {y, 0, 10}, PlotPoints -> 30]
0192名無しさん＠３周年
NGNGr = 4 - 0.5y (is曲線)
g = Plot3D[r, {x, 0, 10}, {y, 0, 10}, ViewPoint -> {1.5, 0.5, 0.7}]
ContourPlot[r, {x, 0, 10}, {y, 0, 10}, PlotPoints -> 30]
164 名前： １３２人目の素数さん 投稿日： 04/06/25 11:35
s = 10y - p^2 (lm曲線 1/pではエラーなので-p^2にした)
g = Plot3D[r, {p, 0, 10}, {y, 0, 10}, ViewPoint -> {1.5, 1, 0.7}]
ContourPlot[r, {p, 0, 10}, {y, 0, 10}, PlotPoints -> 30]
t = Show[r, s]
とすると
Show::"gcomb": "\!\(Show[\(\(2, \
y\)\)]\)のグラフィックオブジェクトの合成処理においてエラーが発生しました．"
と
なってしまう
0193名無しさん＠３周年
NGNGhttp://www.geocities.jp/toritsumacro2004/
0194名無しさん＠３周年
NGNG変数とパラメータの区別ってどうしたらできますか．
例えば，x_^y_ にマッチさせて，
x と y のどれが変数なのかによって，
返す結果を変えたいということなんですけど．
0195名無しさん＠３周年
NGNG変数とパラメータの区別ってどうしたらできますか．
例えば，x_^y_ にマッチさせて，
x と y のどれが変数なのかによって，
返す結果を変えたいということなんですけど．
0196名無しさん＠３周年
NGNG0197名無しさん＠３周年
NGNG0198名無しさん＠３周年
NGNG0199名無しさん＠３周年
NGNGなんか、面白そうな機能がついてるな。
http://www.wolfram.com/products/mathematica/newin51/equationtrekker.html
0200名無しさん＠３周年
NGNGﾋﾟｭ.ｰ （ ・３・） ( ＾＾ ） ＜これからも僕たちを応援して下さいね（＾＾）。
＝〔~∪￣￣￣∪￣￣〕
＝ ◎――――――◎ 山崎渉&ぼるじょあ
0201名無しさん＠３周年
05/02/04 00:23:35ImplicitPlot3D[x^2 + y^2 + z^2, {x, -3, 3}, {y, -3, 3}, {z, -3, 3},
Contours -> {1.5}]
この様な式が本に載っていたのでそのまま書いているのですが、
Options[ImplicitPlot]のタグImplicitPlotはプロテクトされています。
と出てしまいグラフ化できません。
どう言ったことが原因なのでしょうか？
0202ぼるじょあ ◆yBEncckFOU
05/02/20 09:25:04∧＿∧ ∧＿∧
ﾋﾟｭ.ｰ （ ・３・） ( ＾＾ ） ＜これからも僕たちを応援して下さいね（＾＾）。
＝〔~∪￣￣￣∪￣￣〕
＝ ◎――――――◎ 山崎渉&ぼるじょあ
0203名無しさん＠３周年
05/02/21 01:05:24良く見かけるが・・・
0204山.崎 渉
05/02/22 19:37:24━―━―━―━―━―━―━―━―━[JR山崎駅（＾＾）]━―━―━―━―━―━―━―━―━―
∧＿∧
ﾋﾟｭ.ｰ ( ＾＾ ） ＜これからも僕を応援して下さいね（＾＾）。
＝〔~∪￣￣〕
＝ ◎――◎ 山崎渉
__∧＿∧_
|（ ＾＾ ）| ＜寝るぽ（＾＾）
|＼⌒⌒⌒＼
＼ |⌒⌒⌒~| 山崎渉
~￣￣￣￣
∧＿∧
（ ＾＾ ）＜ ぬるぽ（＾＾）
(⌒V⌒)
│ ＾ ＾ │＜これからも僕を応援して下さいね（＾＾）。
⊂| |つ
（＿）（＿） 山崎パン
∧＿∧ ∧＿∧
ﾋﾟｭ.ｰ （ ・３・） ( ＾＾ ） ＜これからも僕たちを応援して下さいね（＾＾）。
＝〔~∪￣￣￣∪￣￣〕
＝ ◎――――――◎ 山崎渉&ぼるじょあ
0205名無しさん＠３周年
05/03/09 14:05:25の時の、x, y, zの組み合わせを使いたいのですが、
どうすれば求められるのでしょうか。
0206名無しさん＠５周年
2006/01/28(土) 09:29:16(n=1,2,3,4,5)という問題です。
これはそれぞれフーリエ正弦級数展開、余弦級数展開せよということなのでしょうか。
前者はFourierTrigSeries[f(x),x,5]、
後者はFourierCosSeries[f(x),x,5]で良いのでしょうか。
正弦はTrigで出るらしいのですが余弦は適当なコマンドが見つからなかったので
Cosをつければ出るのかと思って自分で勝手に作ったのですが・・
ｎが５まで＝５次まで展開という意味なのでしょうか。
教えてくださいですいません。よろしくお願いします。
また展開した関数のグラフを元の関数と同じグラフ上にプロットする方法も
御教示いただけるとありがたいです。
別の板にも同趣旨の書き込みしましたがレスが全然来ないのでマルチ勘弁。
0207名無しさん＠５周年
2006/02/06(月) 15:02:04GUI環境は同じなのでしょうか？それとも行ベースのCUIになりますか？
他のマシン経由でTELNETあるいはSSH、RLOGINでログインして
XでのGUI環境が使えますか？
SMPのマシンの場合に、性能はシングルのマシンに比べて向上できる
のでしょうか？
0208名無しさん＠５周年
2006/02/06(月) 19:58:14このサイトはすごすぎます！！！！
http://www.1yen.tv/stall/intro.cgi?intro=kyonkyon
0209名無しさん＠５周年
2006/02/07(火) 13:13:350210名無しさん＠５周年
2006/02/08(水) 19:14:360211名無しさん＠５周年
2006/03/07(火) 21:50:58ただ、FORTRANほどパッケージはないがな。
0212名無しさん＠５周年
2006/05/12(金) 14:31:400213名無しさん＠５周年
2006/07/13(木) 05:33:06X client のフロントエンドがあるよ。
それより、Microsoft Windows のフロントエンドから
Linux 上の mathematica Kernel を使って計算する方法をどなたか教えて下さい。
rsh の設定はできていて、Linux 上でmathematica Kernel のプロセスは起動するんだけど、
計算結果がMicrosoft Windowsに返ってこないんです。
portブロックはLinux でも Windowsでもしてないし、Windowsファイアーウォールは無効にしたし、、、、
0215名無しさん＠５周年
2006/11/14(火) 23:27:430216名無しさん＠５周年
2007/05/10(木) 15:04:16matlabが6になったの何年前だと思ってんだよ、おめでてーな(ﾜﾗ
0217名無しさん＠５周年
2007/05/13(日) 23:58:34教えて頂けませんか？
0218名無しさん＠５周年
2007/06/16(土) 11:25:530219名無しさん＠５周年
2007/06/27(水) 09:09:33おしえて
0220名無しさん＠５周年
2007/11/27(火) 08:18:17として
f[0.5]
を計算しようとすると
out=f[0.5]
となってしまうのですが、
関数Ｊの形が複雑だとこうなっちゃうんですかね…
解決策ｷﾎﾞﾝﾇです。
教えてｴﾛｲ人
0221名無しさん＠５周年
2007/11/30(金) 01:42:26f[0.5]=D[J[0.5],0.5]
と解釈されて「0.5は変数じゃないから微分できないよ」
と怒られてるんじゃないのかな？
f[x]/.{x->0.5}
で置換しとけばいいと思うよ。
0222名無しさん＠５周年
2007/12/04(火) 17:04:50㌧ｸｽ!
ほんとに助かったありがとう
0223名無しさん＠５周年
2007/12/09(日) 21:13:27いったいどうやって自己解決したのだろうか．非常に気になる．
0224名無しさん＠５周年
2008/04/18(金) 00:29:580225名無しさん＠５周年
2008/10/05(日) 16:34:57光に関連した授業の問で全く分かりません。。。
どなたか力を貸していただけると嬉しいです。
矩形および円形開口のフラウンホーファー回折のパターンを
マスマティカなどで計算して、開口の大きさを変えた時に
パターンがどのように変化するのか調べよ。
まず、マスマティカに入力する式が分かりません。
よろしくお願いします。
0226名無しさん＠５周年
2008/11/08(土) 18:28:37当店の雰囲気に合った服装というのは勿論ですが、
何よりもその服装がお客様ご自身に合っているかどうか、
という点を最も重視してチェックさせて頂いております。
特別に決まりや基準があるわけではなくファッションの基本と同じですので、
難しくお考えにならずにご来店頂ければと思います。
0227名無しさん＠５周年
2008/11/13(木) 17:03:06初期条件、境界条件としてどんなのが使えますか？
0228名無しさん＠５周年
2009/02/14(土) 16:37:18http://www.wolfram.com/news/homeedition.html
295ドル！ 安い。貧乏人の君でも
アルバイトすれば買えるメドが立つだろう？
0229名無しさん＠５周年
2009/05/10(日) 09:11:50For[x2 = 1, x2 < nx2 + 1, x2 = x2 + 1,
im03ft[[y2, x2]] = If[im03f[[y2, x2]] > 30, 1, 0];
]
]
MathematicaでForループ内のIf文の30を変数にして、
Manipulateをつかいたいのですが、うまくできません。。。
どなたか知恵を御貸しください。よろしくお願いします。
0230名無しさん＠５周年
2009/07/29(水) 16:51:530231名無しさん＠５周年
2009/08/02(日) 11:46:47Mathematica はモデリングやシミュレーションには向かないよ。
0232名無しさん＠５周年
2009/08/31(月) 22:56:200233名無しさん＠５周年
2009/10/04(日) 03:22:36誰か出だしだけでも教えてもらえないでしょうか？
0234名無しさん＠５周年
2010/02/27(土) 22:54:550235名無しさん＠５周年
2010/02/28(日) 19:14:57おどろきだな
0236名無しさん＠５周年
2010/02/28(日) 19:48:440237名無しさん＠５周年
2010/04/19(月) 05:55:46たぶん、ネット全体にあらゆる情報が整備されたためだと思うけど。
0238名無しさん＠５周年
2010/04/19(月) 05:57:11Ｃ言語やフォートランで書いた方が、将来の役に立つだろうし。特にＣ言語つかうなら。
0239名無しさん＠５周年
2010/06/02(水) 14:30:57戦闘の数字のリストを作れ。またそのうち、１が何個あるか調べよ。
わからないのですみませんが誰か教えてくれませんか？お願いします。
0240名無しさん＠５周年
2010/06/09(水) 14:22:22ただし、東西南北をそれぞれ、｛１、０｝｛ー１、０｝｛０、ー１｝｛０、１｝で
表すことにする。
（ヒント：りすとdirectotions={{1,0},{-1,0},{0,-1},{0,1}}を作り、1、2、3、4の範囲の
乱数kを発生させ、directotionsのk番目の要素をとるとランダムな方角が得られる）
すいませんが練習だと思って解いてくれませんか？課題が出せないんでお願い
0241名無しさん＠５周年
2010/06/16(水) 13:53:56正整数nのケタ数を与える関数f2を書け。たとえば、f2[40172]＝５
0242名無しさん＠５周年
2010/06/16(水) 15:02:13問題間違えました。
０より大きく１より小さい既約分数で分母がn以下のものを
小さい方から並べてできるリストを出力する関数f5[n]を書け。たとえば、
f5[n]={1/4,1/3,1/2,1/3,1/4}
（このような数列をnに対するFarey数列といい、いろいろおもしろい性質が知られている。
たとえば、隣り合う項をa/b,c/dとすると、必ずab-cd=-1となっている。）
すいませんがお願いします。
■ このスレッドは過去ログ倉庫に格納されています