5ちゃんねる ★スマホ版★ ■掲示板に戻る■ 全部 1- 最新50  

■ このスレッドは過去ログ倉庫に格納されています

mathematicaによるシミュレーション教えてください

1 :名無しさん@1周年 :2000/09/26(火) 22:20
mathematicaを使って時間を変数としてもつ関数を
求めたのですが,この関数の例えば0.1秒刻みの値を
出力させるにはどうしたらよいのでしょうか?
できれば,Exelに出力させたいのですが。
ご存知の方,教えていただけませんか。

156 :あぼーん:あぼーん
あぼーん

157 :名無しさん@1周年:02/12/19 21:59
須藤元気がもっと強くなりますように。

158 :名無しさん@1周年:02/12/21 10:01
mathematicaで多重積分(n重責分)を表現するにはどうすればよろしいですか?
これができないものかと悩んでいます。。。

159 :名無しさん@1周年:02/12/21 14:02
Integrate[f[x,y],{x,0,1},{y,0,1}]
じゃだめ?

160 :158:02/12/21 17:22
2重積分ならわかるんですけどね。
一般にn重となるとどうしていいかわからなくて。。。。
∫…∫dx1dx2・・・・dxn
っていうのをmathmaticaを使って計算させてみたいんです。

161 :名無しさん@1周年:02/12/21 18:44
>>160
具体的にどういう計算をしたいか良くわからんが、
Ver4.1以降ならBoole関数使うと上手く行くかも。

162 :名無しさん@1周年:02/12/21 20:02
>160
つまり、引数が可変ということ?
だとすると、こんなのはどうでしょう(頭で考えただけで未検証)
1)積分変数とその区間の組を List で作る
2)そのリストに被積分関数を Prepend
3)Apply で List を Integrate



163 :162:02/12/21 23:45
もしも Integrate[f[r_vector], Domain] なんてのがご希望でしたら
すいません。

大人しくリーマン和かモンテカルロを構成するのかなぁ。。。

164 :名無しさん@1周年:02/12/22 22:23
あのお、マスマティカで経路積分がしたいんですが、よr

165 :名無しさん@1周年:02/12/22 22:45
>164

時間分割タイプの式をちらっとしか見たことありませんが
ちょっとしたソフトウエアプロジェクト並になると思います。
ググったら、2〜3の米国大学でシラバスのトピックで
見かけた内容ですね。

領域が変なタイプの場合の対応として
積分領域を含む単純タイプの新領域 +
  領域内外区別用の Heviside x 被積分関数
にしてみたり

有限分割で式を組み立てて、 Limit もしくは
分割数大で数値的な収束を見る

C とかでライブラリがありそうな気がしますので、
そのソースを解読し、 Mathematica ならではの
リスト処理をするのが早いと思います。

166 :名無しさん@1周年:02/12/24 23:35
>>164
無限重積分なんてできっこないだろ!


167 :名無しさん@1周年:02/12/26 12:24
単純にxをn回積分するのも簡単にできないのかなぁ?

168 :名無しさん@1周年:03/01/07 23:51
n個のデータをn個のファイルに保存するにはどうすればよいですか?

169 :名無しさん@1周年:03/01/08 11:35
ひとつひとつ、Exportするのはだめ?
多分、任意のn個のdataを自動的にファイルの名前にaaa_###.txt
としてセーブしたいんですよね?
めちゃめんどうです。

170 :169:03/01/08 13:03
パターンの認識は他のプログラム言語(perl、awk etc)
に比べてかなり面倒です。
data[[n]]としてわりあてて
Table[Export[data[[n]].txt,data[[n]],"Table"],{n,1,100}]
として発生しておき、実行するしかないでしょう。
やったことないのでちゃんと動くか分かりませんが…

171 :名無しさん@3周年:03/02/14 17:40
誤差逆伝搬法のアルゴリズムを教えてください.
各層3層でのさくっと書いてくださる方,お願いします.


172 :BP:03/02/19 12:53
>>173
バックプロパゲーションならネットでも簡単にひっかかるだろ
検索してみな

173 :名無しさん@3周年:03/04/16 23:58
「ますます出来るMathematica」
http://www.catkk.co.jp/mathcd/index.html

これ使うと勉強効率あがるかしら?

174 :名無しさん@3周年:03/04/25 16:50
VBAからmathematica動かしたいんですが、
詳しく書かれてるHPとかありませんかね?
wolframのHPだけだと、よくわかりません、、、

175 :名無しさん@3周年:03/05/07 22:15
age

176 :200:03/06/11 19:57
200get
ついでにage

177 :_:03/06/11 20:14
http://homepage.mac.com/hiroyuki44/hankaku03.html

178 :超初心者:04/02/05 02:45
図を描くのに必要なのですが、
ImplicitPlotはどうやったら使えるようになるのでしょうか?
ヘルプにある以下のプログラムを打ち込んでも
下の行の式をそのまま返してくるのですが・・・

<< Graphics`ImplicitPlot`
ImplicitPlot[x^2 + 2 y^2 == 3, {x, -2, 2}]

きちんと各行shift+enterで入力決定しております。
よろしくお願いします。

179 :名無しさん@3周年:04/02/05 06:34
>>178
Packageがうまく組み込めていないかもと思われ。
Input menu->Get File Pathで、<<のあとにパスを入れてミソ
きれいな楕円だね?


180 :超初心者:04/02/05 15:09
>>179
お返事ありがとうございます。
申し訳ありませんが、具体的なプログラムを書いていただけますでしょうか?
<<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}]
ではダメでした。

181 :名無しさん@3周年:04/02/05 20:40
<<"C:\\Program Files\\Wolfram
Research\\Mathematica\\3.0\\AddOns\\StandardPackages\\Graphics\\ImplicitPlot.m"
で、shift+enter
ImplicitPlot[x^2 + 2 y^2 == 3, {x, -2, 2}]
で、また、shift+enter


182 :超初心者:04/02/06 02:09
>>181
大変ありがとうございます。
実は、何をやったわけではないのですが、何度か新規作成を繰り返していたら
なぜか出来るようになってしまいました。不思議です。

183 :名無しさん@3周年:04/03/05 21:10
Solveで連立線形方程式を解こうとしてるのですが、
解きたい変数の数が何十個とあります。

仕方ないので

vars = {v1,v2,…全部の変数を入力してる…,v30}
Solve[ eqns, vars ]

とやってるんですが、これではあんまり大変なので
なにか上手い方法はないでしょうか。


184 :名無しさん@3周年:04/03/07 14:49
>>183
var=Table[v[[i]],{i,1,30,1}];
にした少しは楽になるのでは?

>>182
WindowsでMathematica使うとよくそんなことあるの?(パッケージがうまくよめないとか)
Mac使ってますが、一度もそんなことないなぁ

185 :名無しさん@3周年:04/03/08 20:50
>>183
Solveじゃなくて行列使って計算するのじゃだめなのかな?

186 :名無しさん@3周年:04/04/10 09:36
mathematicaの入門とか、初心者向けの本、教えてください。

187 :名無しさん@3周年:04/04/12 10:07
>182,184
横レスですが、Mathematicaは操作の途中で生成した、
オブジェクト変数名を起動時からずっと記憶しますよね。
途中で重複した変数名が残ってしまって悪さをしたのでは?
それら途中履歴をクリアする、Clear文っていうのもあります。

>186
Mahtematica関係の本は入門書からとんでもない高度な話題のまで
いろいろありますが、自分のオススメは、

牧野書店 Mathematica 数値数式プログラミング
ISBN4-7952-0113-7 C3055



188 :名無しさん@3周年:04/04/19 21:16
初心者です。shift+enterを押しても、計算結果が出力されません。1+1を入力してもエラーがでます。
何が原因なんでしょ。

189 :名無しさん@3周年:04/04/22 14:48
>>188
1+1を1+1っていれているからじゃないか?

190 :名無しさん@3周年:04/04/28 17:21
>>188
まずは使っているOS, Mathematicaのバージョンなどを書いてください

191 :名無しさん@3周年:04/06/25 22:01
マスマティカで きれいに需要曲線の変化を描きたい
たとえば  ケインジアン交差図と投資関数から導かれる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]

192 :名無しさん@3周年:04/06/25 22:04
これを応用して

r = 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 名前: 132人目の素数さん 投稿日: 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\)\)]\)のグラフィックオブジェクトの合成処理においてエラーが発生しました."


なってしまう

193 :名無しさん@3周年:04/06/27 21:10
解決しました
http://www.geocities.jp/toritsumacro2004/


194 :名無しさん@3周年:04/10/18 19:58:19
パターン・マッチを使って関数を定義する際に,
変数とパラメータの区別ってどうしたらできますか.

例えば,x_^y_ にマッチさせて,
x と y のどれが変数なのかによって,
返す結果を変えたいということなんですけど.

195 :名無しさん@3周年:04/10/18 20:04:25
パターン・マッチを使って関数を定義する際に,
変数とパラメータの区別ってどうしたらできますか.

例えば,x_^y_ にマッチさせて,
x と y のどれが変数なのかによって,
返す結果を変えたいということなんですけど.

196 :名無しさん@3周年:04/11/07 13:53:03
1シートあたりのデータ数の最高数(エクセルだと65000くらい)はいくつですか?

197 :名無しさん@3周年:04/11/14 03:13:01
最近来なくなっちゃったみたいだけど、mathematica愛好家って良い香具師だな

198 :名無しさん@3周年:04/12/08 18:41:09
楕円のシュミレーションなんですけど、楕円の回転のさせ方おしえてください。


199 :名無しさん@3周年:04/12/18 11:05:38
バージョン5.1って出てたのか・・・

なんか、面白そうな機能がついてるな。
http://www.wolfram.com/products/mathematica/newin51/equationtrekker.html

200 :名無しさん@3周年:05/01/13 08:50:31
     ∧_∧  ∧_∧
ピュ.ー (  ・3・) (  ^^ ) <これからも僕たちを応援して下さいね(^^)。
  =〔~∪ ̄ ̄ ̄∪ ̄ ̄〕
  = ◎――――――◎                      山崎渉&ぼるじょあ

201 :名無しさん@3周年:05/02/04 00:23:35
<< Graphics`ImplicitPlot`
ImplicitPlot3D[x^2 + y^2 + z^2, {x, -3, 3}, {y, -3, 3}, {z, -3, 3},
Contours -> {1.5}]

この様な式が本に載っていたのでそのまま書いているのですが、
Options[ImplicitPlot]のタグImplicitPlotはプロテクトされています。
と出てしまいグラフ化できません。
どう言ったことが原因なのでしょうか?

202 :ぼるじょあ ◆yBEncckFOU :05/02/20 09:25:04
                                         
     ∧_∧  ∧_∧                             
ピュ.ー (  ・3・) (  ^^ ) <これからも僕たちを応援して下さいね(^^)。
  =〔~∪ ̄ ̄ ̄∪ ̄ ̄〕                             
  = ◎――――――◎                      山崎渉&ぼるじょあ
                                          

203 :名無しさん@3周年:05/02/21 01:05:24
ううむ分からん.山崎渉は分かるが,ぼるじょあって何だ?
良く見かけるが・・・

204 :山.崎 渉:05/02/22 19:37:24
...これからも僕を応援して下さいね(^^)。   
  
━―━―━―━―━―━―━―━―━[JR山崎駅(^^)]━―━―━―━―━―━―━―━―━―
         
     ∧_∧
ピュ.ー (  ^^ ) <これからも僕を応援して下さいね(^^)。                         
  =〔~∪ ̄ ̄〕                                            
  = ◎――◎                      山崎渉                       
                                
 __∧_∧_                                                 
 |(  ^^ )| <寝るぽ(^^)      
 |\⌒⌒⌒\                                
 \ |⌒⌒⌒~|         山崎渉             
   ~ ̄ ̄ ̄ ̄                            
                            
   ∧_∧                                       
  (  ^^ )< ぬるぽ(^^)      
                                                       
    (⌒V⌒)                    
   │ ^ ^ │<これからも僕を応援して下さいね(^^)。   
  ⊂|    |つ                                
   (_)(_)                      山崎パン 
                                         
     ∧_∧  ∧_∧
ピュ.ー (  ・3・) (  ^^ ) <これからも僕たちを応援して下さいね(^^)。
  =〔~∪ ̄ ̄ ̄∪ ̄ ̄〕                          
  = ◎――――――◎                      山崎渉&ぼるじょあ

205 :名無しさん@3周年:05/03/09 14:05:25
x + y + z =1 ( x>0 , y>0, z>0 )

の時の、x, y, zの組み合わせを使いたいのですが、
どうすれば求められるのでしょうか。

53 KB
■ このスレッドは過去ログ倉庫に格納されています

★スマホ版★ 掲示板に戻る 全部 前100 次100 最新50

read.cgi ver 05.04.00 2017/10/04 Walang Kapalit ★
FOX ★ DSO(Dynamic Shared Object)