〓 Mathematica 2 〓
■ このスレッドは過去ログ倉庫に格納されています
0001１３２人目の素数さん
NGNG, - ' ,..､､.ヾ{{フ'⌒`ヽ､
／ ,:', -‐‐` ´ '´⌒ヽ ヾ:､
. ,' ,'´ ,ィ ,ｨ ,' , ｀ヽ', ',-<
,' .i /|. /.| { i, i, }. }_,,）)
! | ! .,'-.{ ! !|; |`､.}ﾞ!.! |. ! ヽ.
', ', |Vｧ=､ﾞ､ ｀ﾞ､!-_:ﾄ,ﾘ', l ! | ﾞ',
ヽ､', l:!Kノ}. f:_.)iﾞi: ﾘ ! l ル
| l!iヽﾞ- ' , .!__:ﾉ ﾞ ,ﾘ l ﾘ'´
. ',|!!､ r‐┐ ｀ ノ'. /,ｲ ／￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣￣
'i!ﾞ､ヽ､ ﾞｰ' _, ィ,:',:''´ ＜ Mathematica に関する話題はここに書いてね！
ﾞ:､ｨ､jヾｰ::: 'iﾍ .ノ',ﾘ. ＼＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿
,､- '´ ヽ､ﾞ､ { ｀＞"､
/＼＼ ', } //｀ヽ
★前スレ
〓Mathematica〓
http://science2.2ch.net/test/read.cgi/math/985023298/
0083KingMathematician ◆5lHaaEvFNc
NGNG出来ないはずはない。
Plot[x^5-2x^4+3x^3+x^2-x+1,{x,-3,3}]
Plot[Exp[x]-3Sin[x]-5,{x,-3,3}]
Plot[Exp[x]+1/(1-x^2),{x,-3,3}]
NSolve[x^5-2x^4+3x^3+x^2-x+1==0,x]
NSolve[Exp[x]-3Sin[x]-5==0,x]
NSolve[Exp[x]+1/(1-x^2)==0,x]
(現在手許にMathematicaが無いのでこれで大丈夫かどうか…。)
0084KingMathematician ◆5lHaaEvFNc
NGNG0085１３２人目の素数さん
NGNG転置のコマンドってTransposeのほかにあるんすか？
Transpose使ったら、列ベクトルとして出力されずに、行列になっちゃうんですよね。
行列としてではなく、列ベクトルで表示したいんですよ。
基底の組が二つなら、列ベクトルが二つ出力されるような。
0086KingMathematician ◆5lHaaEvFNc
NGNG行列でいいんじゃないの？
0087１３２人目の素数さん
NGNGやっぱりそうですか・・・俺もそう思うんですよね。
そんなコマンドどう考えてもないような気がして。
0088KingMathematician ◆5lHaaEvFNc
NGNGあるいは、行列を得たうえ、
MatrixForm関数で出力してみるとか…。
0089１３２人目の素数さん
NGNG0090１３２人目の素数さん
NGNG振り子が振れるごとに針が動くような。
お願いします(ToT)
0091１３２人目の素数さん
NGNGもっと問題点を絞ってくれよう。
振り子の運動方程式がわからないのか、方程式の解き方がわからないのか、
グラフィックの描き方があわからないのか、アニメーションにしたいのか・・・
0092１３２人目の素数さん
NGNGすいません。グラフィックが書きたいです。
できれば、振り子が一回振れるごとに秒針が一秒進んでいくような時計を作りたいと思っています。
どなたか宜しくお願いしますm(_ _)m
他人の力に頼りすぎというのは分かってますが、どこから始めていいのか分からないです。
0093１３２人目の素数さん
NGNG0094１３２人目の素数さん
NGNGNSolve[Exp[x]+1/(1-x^2)==0,x]
で、計算させると次のような注意をされてしまいます。
Solve::tdep:
本質的に非代数的な方法で解かれる変数が方程式に含まれているようです
解決策はありますか？
Plotの方は上手く行ったみたいなのですが・・・。
お願いします。
0095１３２人目の素数さん
NGNGFindRoot[Exp[x] - 3Sin[x] - 5 == 0, {x, 2}]
ヘルプを引け。See alsoぐらい読め。
0096１３２人目の素数さん
NGNG遂に、怒られてしまった。
でも、ありがとうございます。
0098１３２人目の素数さん
NGNGGraphicsも綴れない上にコピペすらできないとは
さすが教えて君ですね。
0099１３２人目の素数さん
NGNG0101１３２人目の素数さん
NGNG0102１３２人目の素数さん
NGNG010392
NGNG0104１３２人目の素数さん
NGNGそういう問題じゃないよ。
要求するグラフィックのレベルがわからんし、激しく面倒。
掲示板でちょこっと張れるレベルじゃない。
PolygonとかLineとかつかって地道に描け。
あるいは外部のドローソフトやモデラーでつくって取り込むか。
シンプルな2Dでよければ、Arrow, Circle, Rectangleあたりを組あわせて、
時計を描いて、振り子の運動方程式にしたがって振り子部分を動かす。
0105１３２人目の素数さん
NGNG分からないことがあったら、聞きますので、その時はよろしくお願いします。
0106１３２人目の素数さん
NGNGコピペで書けますけど・・・。
微分積分はライフワークですよね！
教えて君に救いの手を・・・。
0108１３２人目の素数さん
NGNGあれ？できた！！できたぞー！！
0109１３２人目の素数さん
NGNGバカですいません。
0110KingMathematician ◆5lHaaEvFNc
NGNG振り子は、運動がある球面に束縛されていると考えると良いだろう。
当然、滑らかな束縛である。
0111１３２人目の素数さん
NGNG式で具体的に書くとどうなりますでしょうか？
教えて君ですいません。
0112KingMathematician ◆5lHaaEvFNc
NGNGとりあえず、重力加速度は一定の値gをとることにしよう。
振り子の質量をmとして、rを振り子の長さとして、振り子の鉛直下向きからの角度をθとすると、
振り子に働く張力と重力の合成は、
mgsin(θ)という大きさで、
方向は振り子の「軸」に向かう接線であり、
向きは下に向かう方向である。(θ=0,πの場合は0だから方向も向きもない。)
とりあえず、振り子の中心を、三次元空間の原点にもってきて、
鉛直上向きを三次元方向の正の向きとすると、
(x,y,z)(ただし、x^2+y^2+z^2=r^2)における振り子の運動方程式(振り子の座標が(x,y,z)である。)は、
m(x'',y'',z'')=-mgsin(θ)(cos(φ)z/√(x^2+y^2+z^2),sin(φ)z/√(x^2+y^2+z^2),√(x^2+y^2)/√(x^2+y^2+z^2))
但し、φ=Arctan(y/x)(x>0),Arctan(y/x)+π(x<0),π/2(y>0,x=0),3π/2(y<0,x=0)であり、
θ=Arctan(z/√(x^2+y^2))(x^2+y^2>0),0(z=-r),π(z=r)である。
検算よろしく。
0113１３２人目の素数さん
NGNGありがとうございます。２Ｄで書こうと思ってるので、z座標は考えなくていいですよね。
本当にありがとうございました。
0114１３２人目の素数さん
NGNG0115１３２人目の素数さん
NGNG0116１３２人目の素数さん
NGNGttp://up.isp.2ch.net/up/7c3ec30d5063.nb
こう書いてるのに正しいファイルが出力されない…
どなたか原因分かる方いらっしゃいませんか？
Mathematica4の日本語化したものです
0117１３２人目の素数さん
NGNG> こう書いてるのに正しいファイルが出力されない…
プログラムは正常に動作しているようだが、
君の脳内「正しいファイル」というのを他人にもわかるように
説明してくれないか。
いったい何がやりたいわけ？
0118１３２人目の素数さん
NGNGが有って
その和を取りたいのですが
Σfi
とやるとf*1＋f*2…と認識されるのですが
どうすれば良いのでしょうか…
0119KingMathematician ◆5lHaaEvFNc
NGNGf1,f2,f3,などもまたシンボルとなる。
第i項をfiなどとは書けない。
これは、素直にf1+f2+f3+…
と書くのが早いだろう。
関数がとても多い場合は、
f[i,x]のように関数を定義し直した方がいいと思う。
他にも方法はあるかもしれない。
0120１３２人目の素数さん
NGNG0121１３２人目の素数さん
NGNGf[i,x]と、擬似的な２変数関数にするのは上手い方法ですね。
0122１３２人目の素数さん
NGNGf[i][x]と定義するのが良い。
0123１３２人目の素数さん
NGNG0124１３２人目の素数さん
NGNG計算の途中式を吐かせることはできるんでしょうか
0125１３２人目の素数さん
NGNG0126１３２人目の素数さん
NGNGsin(2x^2)というときに
(2x^2)'*cos(2x^2)
＝4xcos(2x^2)
などという風に途中の計算を入れるということはできるのでしょうか
0127１３２人目の素数さん
NGNGx = n;
d2 = d1
これじゃだめかな？
0128１３２人目の素数さん
NGNGNumberFormのExponentStepみたいなオプションを指定して、
指数表記させないことってできるのですか？
ぐぐっていろいろ調べたけど、どこにも書いてなくて....
ご存知の方、教えてください。
環境：Ver 5.0.1（OS X）
0129１３２人目の素数さん
NGNGここに色々あるから見てみると良いと思う
0130１３２人目の素数さん
NGNGうぅ〜ん、やはりでていないですね。
日本語サイトはだいたい見たとこばかりだし...
どうもありがとう。
0131１３２人目の素数さん
NGNGSolve関数などを用いて、導いた複数の解をそのまま次の演算に用いたい場合、
どうすればいいのでしょうか？
例えば
Solve[{2x+3y+5z==12,
5x-3y-4z==4,
3x+2y-z==-2},{x, y, z}]
などとして方程式を解かせ、その直後にこの解を下のような式に代入させ計算したいとします
x^2+y^3+z
どういう風にやればいいのでしょうか？初歩的ですいませんが、教えてください。お願いします。
0132１３２人目の素数さん
NGNGParametricPlot3Dで書くとどうもギザギザで…
0133KingOfKingMathematician ◆H06dC8bpwA
NGNGPlotPointsとかいうオプションなかったっけ？
0134１３２人目の素数さん
NGNG組込み関数の途中計算を知ることはできない。
Mathematicaの式として定義されているものならTraceで見れる。
MyD[x_^n_, x_] := n x^(n - 1)
MyD[f_ g_, x_] := MyD[f, x]g + f MyD[g, x]
MyD[f_, x_] /; FreeQ[f, x] := 0
MyD[2 x^2, x] Cos[2 x^2] // Trace
>>128
自分でFormatを定義せよ。
Format[MyForm[expr_]] := FromDigits[RealDigits[expr]]
Plot[10^x, {x,0,10}, FormatType -> MyForm]
>>131
x^2 + y^3 + z /. Solve[...]
0135１３２人目の素数さん
NGNG<<Graphics`
Sphere[1,30,30]//Graphics3D//Show
0137KingOfKingMathematician ◆H06dcQwwvA
NGNGH06dF0OtH. : #px:wSH,c
H06ddWXDNY : #um1m+｢MN
H06dcuQgY2 : #B45L0K\y
H06deD4siU : #6Oｦj*ｾ?ｻ
H06d2xYhm6 : #Fｻn4wJ+.
H06dQfJLXY : #3)H0x'U/
H06dgGdMtc : #W5IBF@Gz
H06dWwDj/. : #@=c\X8｣x
H06d4z7h8. : #8,:M0Fwf
H06d/CP0JY : #=!7｣*~mp
H06dpnDLQ6 : #xTfjQF3b
H06dhKAsgc : #,2C1n{5-
H06dEu.odw : #WYmH+?Sq
0138KingOfKingMathematician ◆H06dC8bpwA
NGNG0139１３２人目の素数さん
NGNG鬱だ
0140１３２人目の素数さん
NGNG0141１３２人目の素数さん
NGNG0142128
NGNGFormatで書式を指定する方法ではPlotLabel、FrameLabelを設定していると
エラーになってしまいます...
ということで、いろいろ悩んでFrameTicksのオプションを保存して書き換える
方法で解決しました。MathematicaのFAQに書いてありました。
fticks = FullOptions[gl1, FrameTicks];
newfticks = fticks /. {
x_, x_, {0.00625,
0.}, {GrayLevel[0.], AbsoluteThickness[0.25]}} :> {
x, FromDigits[RealDigits[x]], {
0.00625, 0.}, {GrayLevel[0.], AbsoluteThickness[0.25]}};
Show[gl1, FrameTicks -> newfticks]
もっと簡潔な方法ってありますか？
0143KingOfKingMathematician ◆H06dcQwwvA
NGNGH06dABjF.Q : #S\_I//x｢
H06dligu62 : #E?ｦｻy4E6
H06dUsR8lg : #｣%'ｼ|｢N2
H06dRCIVrg : #F?Z3vｦa7
H06dhLHyjw : #0\SE8.99
H06dgb.cI6 : #sｻ.Z[vb0
H06d/75ZQQ : #*Kl]U67h
H06dS8zXQQ : #S:g5xxr]
0144１３２人目の素数さん
NGNG型毎に違う出力をするFormatを書く。
Format[MyForm[expr_Real]] := FromDigits[RealDigits[expr]]
Format[MyForm[expr_]] := OutputForm[expr]
Plot[10^x, {x, 0, 10}, PlotLabel -> "10^x", Frame -> True,
FrameLabel -> {Bottom, Left, Top, Right}, FormatType -> MyForm]
0145128
NGNGありがとうございます。
なるほど、型ごとの指定というのが当然ながら可能なのですね。
大変勉強になります。
ヘルプ、オンライン等で確認しましたが、その他expr_*で指定できる型に
ついて参照することができませんでした。探し方が悪いのでしょうか。
ヘルプならどこを参照すればよろしいでしょうか。
0146１３２人目の素数さん
NGNGMathematicaの言語としての基礎だ。
Principles of Mathematicaを全部読めといいたいところだが、
とりあえず、必ず「Expressionsの章を読んでから」、Patternsの章、
特に、2.3.4を読め。
0147KingOfKingMathematician ◆H06dC8bpwA
NGNG0148１３２人目の素数さん
NGNG何をやっとるのか知らんがスレ違いでうざい。
0149KingOfKingMathematician ◆H06dyzvgzA
NGNGトリップ ザ ベスト
H06dWILLhA : #/{\@%YwX
H06djy9xBA : #SgHdO'H%
H06dYXOYLA : #*｢A@?NVF
H06dhKnt9A : #[Aｼsudｾl
H06dWifa1A : #{SfbN(6ｦ
H06dyzvgzA : #QAiEｼEp- ←使用中
0150KingOfKingMathematician ◆H06dC8bpwA
NGNG0151１３２人目の素数さん
NGNGp[x_] = a[x + i];
q[x_] = b[x + n - i];
(*Cauchyの積級数*)
Sum[Sum[p[x] q[x], {i, 0, n}], {n, 0, Infinity}]
0152１３２人目の素数さん
NGNGしようとしているのですが、どうもうまくいきません。
バージョンは4.1です。
<< Utilities`Notation`
Symbolize[expr_x]
Symbolize[expr_y]
Symbolize[expr_z]
p={px,py,pz}
q={qx,qy,qz}
p.q
（xyzは下付けの添字）
をやってみたのですが、どうもシンボル化ができていないみたいで、必ず強制終了
することになりうまくいきません。expr_を使わず変数を一つずつSymboliseに入れてみても
もだめでした。よくは分からないんですが<< Utilities`Notation`をNeeds["Utilities`Notation`"]
にしても変化なしです。サポートページにキーボードでなくパレットから、とあった
ので基本的にSymbolizeはパレットから持ってきています。ただ、一応キーボードで
やってみてもだめでした。
何をどうすればいいのでしょうか？
0153１３２人目の素数さん
NGNGパレットのSymbolizeをクリック。
ノートブックにフォーカスを写して
expr_<Ctrl+->x<Shift+Enter>
これで定義される。
ただし、<A+B>はキーAを押しながらキーBを押すことをあらわす。
a<Ctrl+->x // FullForm
でa_Subscript_xとなることを確認。
0154152
NGNGありがとうございます。できました！
ただ、152の計算で定義したあと一度目のShift+Enterでは
入力ミスではないですか。q_Subscript_ｙはp_Subscript_ｙに似ています
・
・
・
と似たようなエラーがいくつか出たんですが、消してもう一度Shift+Enter
を押したらちゃんとできました。
よく分かりませんが、使えるよう勉強させてもらいます。
0155１３２人目の素数さん
NGNG> 入力ミスではないですか。q_Subscript_ｙはp_Subscript_ｙに似ています
入力ミスじゃないんだから無視しろ。
0156１３２人目の素数さん
NGNG0157155
NGNG>>152
そのメッセージはwarningであってerrorじゃない。親切に入力ミスの可能性を
指摘してくれているだけで、変数の定義自体は正常に行われているので、
入力ミスでない今回の場合は無視すれば良い。
0158152
NGNGアドバイスありがとうございます。
そんな親切な機能があったのですね、分かりました。
0159１３２人目の素数さん
NGNGマネーサプライが増えた時と ぴグー効果が起こった時 の状態を描きたいんですが
まったくわかりません
とにかく
plot してみたんだけど グラフが ショボイっていうか軸により過ぎててわかんない
どうやって調整するんでしょうか
あと それぞれのグラフを色分けしたいんですが タグがまったくわからん
ひとつでも間違うとちんぷんかんぷん
ISLM分析 にPを付け加えた 3次元分析でも それぞれの定数が変化した時に起こる
需要曲線の傾き変動をきれいに江が北印ですが
だれか
ここにタグ貼ってください よろしくお願いします
0160１３２人目の素数さん
NGNGという基本タグ
をどのように拡張して行けば それぞれの線をきれいに色分けして きれいなグラフで
シフトがわかるようになるのかということです
0161１３２人目の素数さん
NGNGどのようにシフトするか分析するために
三次現プロットしなければ生けないんですがこのタグも出来たらわかるかたいませんかね
0162１３２人目の素数さん
NGNG例えば コブダグラス型の生産関数です
ほんとは k=1/3 l=2/3 なんですが
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]
0163１３２人目の素数さん
NGNGr = 4 - 0.5y
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]
0164１３２人目の素数さん
NGNGg = Plot3D[r, {p, 0, 10}, {y, 0, 10}, ViewPoint -> {1.5, 1, 0.7}]
ContourPlot[r, {p, 0, 10}, {y, 0, 10}, PlotPoints -> 30]
isもどき 1/pなんた゜か゜ なんだかエラーになってしまうので変えた
0165１３２人目の素数さん
NGNG0166１３２人目の素数さん
NGNG夜露死苦って素数なんだぜ！
分かりにくく答えが出しにくい青春時代を良く表してるよね！
0167１３２人目の素数さん
NGNGhttp://www.catkk.co.jp/mathcd/index.html
これ使うと勉強効率あがるかしら？
0168１３２人目の素数さん
NGNG0169１３２人目の素数さん
NGNGr = 2
g = Plot3D[r, {x, 0, 10}, {y, 0, 10}, ViewPoint -> {1.5, 0.5, 0.7}]
s = y
g = Plot3D[s, {p, 0, 10}, {y, 0, 10}, ViewPoint -> {1.5, 1, 0.7}]
gjy = Show[r, s]
Show::"gcomb": "\!\(Show[\(\(2, \
y\)\)]\)のグラフィックオブジェクトの合成処理においてエラーが発生しました．"
0170１３２人目の素数さん
NGNGミクロのコブダグラス関数しか表示してるサイト無い
ISLM分析から接線をADであらわしてるサイトなんて無い
だれか うまくシミュレートできるタグわからないかな
0171１３２人目の素数さん
NGNGもう使いたく無い
Flash MX よりかは扱いやすいが
ヘルプは醜くて拉致があかない
ミクロ分析の解説本はで照るみたいだけど
マクロも出るの解説本でてないかな
0172１３２人目の素数さん
NGNGr = 2, s = y なんだから Show[r, s] = Show[2, y]だね。
Show[2,y]でなにが表示されると？
ここまでの書き込みを読んでいても思ったのだが、君もうちょっと頭使わないと。
0173KingOfKingMathematician ◆H06dyzvgzA
NGNGKingOfKingMathematicianの後に付けるのがおしゃれ。
H06dWILLhA : #/{\@%YwX
H06djy9xBA : #SgHdO'H%
H06dYXOYLA : #*｢A@?NVF
H06dhKnt9A : #[Aｼsudｾl
H06dWifa1A : #{SfbN(6ｦ
H06dyzvgzA : #QAiEｼEp- ←使用中
0174１３２人目の素数さん
NGNGできませんよ
あなたこそちゃんとタグの使い方マスターしてくださいね
Show::"gcomb": "\!\(Show[\(\(2, \
y\)\)]\)のグラフィックオブジェクトの合成処理においてエラーが発生しました．"
0175１３２人目の素数さん
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.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 で 数理的に説明するよりも
自分で勝手に三次元グラフを作ってしまったほうが手っ取り早いんですね
0176１３２人目の素数さん
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 で 数理的に説明するよりも
自分で勝手に三次元グラフを作ってしまったほうが手っ取り早いんですね
0177１３２人目の素数さん
NGNGこの式から見て Mの増加や
C(Y-{T+kY}) = a+b(Y-{T+kY}) となる累進課税時に
曲線の傾きがどうなるかは簡単にわかりますが
実際にグラフで視覚的に体験して見たいのです
0178１３２人目の素数さん
NGNGそれが需要曲線の変動とどうつながるか
三次元グラフを使って立体的に示したいんです
0179１３２人目の素数さん
NGNG0180KingOfKingMathematician ◆H06dC8bpwA
NGNG0181KingOfKingMathematician ◆H06dC8bpwA
NGNG0182KingOfKingMathematician ◆H06dC8bpwA
NGNG■ このスレッドは過去ログ倉庫に格納されています