〓 Mathematica 伍 〓
■ このスレッドは過去ログ倉庫に格納されています
0001１３２人目の素数さん
2011/11/11(金) 09:02:28.82, - ' ,..､､.ヾ{{フ'⌒`ヽ､
／ ,:', -‐‐` ´ '´⌒ヽ ヾ:､
. ,' ,'´ ,ィ ,ｨ ,' , ｀ヽ', ',-<
,' .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/
〓 Mathematica 2 〓
http://science3.2ch.net/test/read.cgi/math/1078534285/
〓 Mathematica 3 〓
http://science6.2ch.net/test/read.cgi/math/1121413040/
〓 Mathematica 四 〓
http://kamome.2ch.net/test/read.cgi/math/1197270001/ (鯖飛びでログ消滅)
http://mimizun.com/log/2ch/math/1197270001/
〓 Mathematica 5 〓
http://kamome.2ch.net/test/read.cgi/math/1285859504/
0552１３２人目の素数さん
2013/11/09(土) 01:02:05.66http://reference.wolfram.com/mathematica/ref/LatticeReduce.html
の例題のアプリケーションを見ると良いと思う。
0553551
2013/11/09(土) 13:11:47.000554１３２人目の素数さん
2013/11/10(日) 01:48:52.81>>552のリンク先のLatticeReduceのドキュメントだけど、
アプリケーションの2番目の例題で
自明な初期関係：
a={
{1, 0, 0, 0, -a0},
{0, 1, 0, 0, -a1},
{0, 0, 1, 0, -a2},
{0, 0, 0, 1, -a3},
{0, 0, 0, 0, -a4}
};
ってなってるけど、間違ってるよね？
3つの例を統一的に表現できるコード書いてみた。
t = 3^(1/3);
v = {1, t, t^2, t^3, t^4};
n = 7;
c = Round[10^n v];
A = Join[IdentityMatrix[Length[v]], -Transpose[{c}], 2]
B = LatticeReduce[A]
A.Append[c, 1]
B.Append[c, 1]
B[[;; , ;; -2]] // MatrixForm
r = B[[;; , ;; -2]].v
N[r]
0555１３２人目の素数さん
2013/11/10(日) 01:59:57.67x*π+y*e=123.456
これを上記のコードにあてはめて
v = {Pi, E, 123.456};
n = 10;
c = Round[10^n v];
A = Join[IdentityMatrix[Length[v]], -Transpose[{c}], 2]
B = LatticeReduce[A]
A.Append[c, 1]
B.Append[c, 1]
B[[;; , ;; -2]] // MatrixForm
r = B[[;; , ;; -2]].v
N[r]
を実行してみると
-6747π+3256e+100*123.456≒0
が得られて、変形すると
(6747/100)π-(3256/100)e≒123.456
となるんだけど、これは求めたい解じゃないよね？
つまり、x,yは整数じゃないといけないんだよね？
となるとLatticeReduceで解けなくない？
0556１３２人目の素数さん
2013/11/10(日) 15:06:57.820557551
2013/11/11(月) 01:07:35.82ご指摘の通り係数は整数のものがいいので、LatticeReduceは使うのはやめました
この問題に関しては効率は悪いですがとりあえず別のプログラムを組んで解決しました
でもこういう処理するいい方法がないのか気になりますね
0558551
2013/11/11(月) 01:10:45.510559551
2013/11/11(月) 01:15:59.880560１３２人目の素数さん
2013/11/11(月) 19:19:11.00ttp://twicsy.com/i/AehpGb
ハート形の立体なんておもしろいなぁとか思いつつ、
記載されてる数式を眺めてたんだけど、
「320」ってなんだっ！？と思ったわけ。
明らかに必要ないからね。
んで、いろいろ調べてるうちに以下に辿り着いた。
ttp://demonstrations.wolfram.com/EquationsForValentines/
Taubinっていうのがオリジナルでそれを定数倍したのが
Trott（ウルフラム社員）になってることがわかる。
つまり、本来は同じ結果が出なきゃいけないのに定数倍することによって
じゃっかん滑らかになったからってそれを大々的にアピールするとかｗｗｗ
ウルフラムって馬鹿なの？自社製品を宣伝するんじゃなくて、
自社製品のバグを宣伝してどうすんだよｗｗｗ
0561556
2013/11/11(月) 22:34:17.50>あと整数の範囲は自然数です
>あとa,b,c,d,e,Aは全て正の数です
としても、最適解でなかったり、等価な解を拾えなかったり。
Minimize[{Abs[1 a + 2 b + 3 c + 4 d + 5 e - 100],
Element[a|b|c|d|e,Integers]
&& (a>0 && b>0 && c>0 && d>0 && e>0)},
{a,b,c,d,e}]
駄式ですが最初に見つかった解を返して終了してんだろな。
結局関数を探すよりコードを自分で書いた方が早いと思う。
0562１３２人目の素数さん
2013/11/12(火) 00:29:55.20プログラムの例等って有るのでしょうか。
0563１３２人目の素数さん
2013/11/12(火) 23:33:34.60>>555の問題解かしたらちゃんと答え出たぞ。
Minimize[{Abs[Pi x + E y - 123.456],
Element[x | y, Integers] && (x > 0 && y > 0)}, {x, y}]
{0.091152, {x -> 35, y -> 5}}
これが大域最適解じゃない可能性があると言っている？
0564１３２人目の素数さん
2013/11/13(水) 22:10:27.49問題条件は>>558，559で修正されました。
領域設定可能なMin関数は幾つかあります
が、大域解の保証もされてないってこと。
0565１３２人目の素数さん
2013/11/14(木) 20:25:39.28＞問題条件は>>558，559で修正されました。
だから何？
＞領域設定可能なMin関数は幾つかあります
Minimize以外に何がある？
＞が、大域解の保証もされてないってこと。
「f および cons が線形あるいは多項式の場合，Minimizeは常に大域的な最小値を求める．」
ってあるから>>563はだめだけど以下なら大域最適解なんじゃない？
Minimize[{(Pi x + E y - 123.456)^2,
Element[x | y, Integers] && (x > 0 && y > 0)}, {x, y}]
{0.00830869, {x -> 35, y -> 5}}
0566１３２人目の素数さん
2013/11/14(木) 22:31:17.76かみつかないでよ･･･(いちいち覚えてないです、これで終わり)。
1)だから何？
問題修正で符号がすべて正になりました(正負で相殺できない)。
ではPi x-E y=0の最適近似値はどこでしょう(答え：桁数次第
2)Minimize以外･･･
MinimizeでHelpすればMinValue,ArgMin,NMinimizeとかがでます、
つか、やたら関数が増えて、説明をみても違いが分からん。
3)＞が、大域解の保証もされてないってこと。
Minimize[{Abs[Pi x - E y - 123.456], (*Eの符号が反転ですが線形*)
Element[x | y, Integers] && (x > 0 && y > 0)}, {x, y}]
どうなりますか、{x -> 41, y -> 2}ならはずれ
0567１３２人目の素数さん
2013/11/15(金) 01:04:07.19(Pi x + E y - 123.456)^2は多項式だから、マニュアル通りに解釈すると、
結果は大域最適解のはずなんじゃないかということ。
「{x -> 41, y -> 2}ならはずれ」ということらしいが、
これよりも良い解の例を挙げてくれる？
Minimize[{(Pi x - E y - 123.456)^2,
Element[x | y, Integers] && (x > 0 && y > 0)}, {x, y}]
{0.00761516, {x -> 41, y -> 2}}
0568狸 ◆BvcplLXSGo
2013/11/15(金) 01:53:51.67○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●
●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○
○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●
●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○
○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●
●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○
○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●
●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○
○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●
●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○
○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●
●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○●○
0569１３２人目の素数さん
2013/11/16(土) 22:12:04.25-0.0872649
Pi x - E y - 123.456 /. {x -> 54, y -> 17}
-0.0207878
|
Pi x - E y - 123.456 /. {x -> 1412253, y -> 1632134}
-2.8871*10^-8
これでええかマニュアルくん
0570１３２人目の素数さん
2013/11/17(日) 10:11:49.10thx.
納得した。
またマニュアルのミスか。
ところで後半2つの解の算出方法は？
あと、解の精度を無限に高められることを直感的に理解する方法ってある？
なんかディオファントス近似とかディリクレの定理とかが関係してそうな気がするんだけど。。。
0571１３２人目の素数さん
2013/11/17(日) 22:56:09.39｢Minimizeが近似数を含む式に使われた場合は，自動的にNMinimizeが呼ばれる。｣
→NMinimizeに問題(近似数でなくすとMinimizeは解を放棄する)
>ところで後半2つの解の算出方法は？
Brute Force Attack...
>解の精度を無限に高められることを直感的に理解する方法
知らない。。。
この場合は超越数が係数でxは無限まで動けるため、無限
走るほどに最適解が生まれる。と思った。。。
0572１３２人目の素数さん
2013/11/22(金) 12:36:23.96http://ignt.hk/1aAiWv5
ついにタダ！
0573１３２人目の素数さん
2013/11/22(金) 12:49:58.630574１３２人目の素数さん
2013/11/23(土) 01:22:17.16iphone のsimフリー版の日本での販売開始
といい、最近，世の中急に慌ただしくなってきな
感じがするね。
0575１３２人目の素数さん
2013/11/23(土) 07:51:09.54「日本のユーザがたくさん支払ってくれるので
世界の子供たちにMathematicaを只で使ってもらえるのです」
0576１３２人目の素数さん
2013/11/23(土) 09:59:40.74今さらmathematicaでもないだろ
0577１３２人目の素数さん
2013/11/23(土) 23:14:57.27あったら、普通にPCの速度でMathematicaが動かせるじゃないか。
0578１３２人目の素数さん
2013/11/24(日) 07:57:21.61>Paspberry Pi のPCエミュレータ
Windows版があるようだ
0579１３２人目の素数さん
2013/11/24(日) 09:26:01.49OSX+QEMUでRaspberryPiをエミュレートする at Curb_Press
http://ignt.hk/IcackH
0580１３２人目の素数さん
2013/11/24(日) 09:27:31.39hirax.net::「Wolfram CDF PlayerをMathematicaとして使う方法」をRubyでもっと簡単にしてみた
http://ignt.hk/18ZVaoz
0581１３２人目の素数さん
2013/11/24(日) 12:11:24.130582１３２人目の素数さん
2013/11/25(月) 16:03:27.00mathematica動かしてみたが、遅くて話にならなかった、、、。
rasberrypiでより遅いかも、、、。
やり方まずいのかな、、、。qemuとか初めてだし、音もでないよ。
0583１３２人目の素数さん
2013/11/25(月) 16:15:32.80エミュは遅すぎ。
ラズパイ複数のクラスタって出来るのかなぁ？
0584１３２人目の素数さん
2013/11/25(月) 16:17:29.78どっちも古すぎだし、Macで動かないorz
0585１３２人目の素数さん
2013/11/25(月) 18:11:25.06大学生なら学校が提携してて無料でインストールできる場合があるので
学校に聞いてみるといい。
あとは大学の生協でアカデミックを買うかだな。
0586１３２人目の素数さん
2013/11/25(月) 18:39:02.98いつまでもボッタクリしてろや。そのうち完全に只になる。
0587１３２人目の素数さん
2013/11/25(月) 20:42:12.08買うのは無理
生活保護の貧乏人だもん
0588１３２人目の素数さん
2013/11/25(月) 20:46:12.220589１３２人目の素数さん
2013/11/25(月) 21:24:42.97日本価格で買うぐらいなら
Maxima,Sage,Octave,REDUCEで頑張る
0590１３２人目の素数さん
2013/11/25(月) 21:48:33.73でも起動するのが面倒でQEMUになる
0591１３２人目の素数さん
2013/11/29(金) 20:37:41.060592狸 ◆2VB8wsVUoo
2013/11/29(金) 22:55:42.32ケケケ狸
0593１３２人目の素数さん
2013/11/30(土) 16:57:31.72MacBook AirでVNC接続のRaspberry Pi
簡単なものならサクサク動くよ
0594１３２人目の素数さん
2013/11/30(土) 16:59:47.30iPhone 4Sからssh接続
グラフィックは無いけど、何処でもmathematica
0595１３２人目の素数さん
2013/11/30(土) 20:26:28.28593みたいにVNCでPCから接続するのが一番良いな。
RaspberryPiは常時起動していても、電気代はたいしたことない
だろうし。
0596１３２人目の素数さん
2013/12/01(日) 11:29:43.820597１３２人目の素数さん
2013/12/01(日) 16:38:15.340598狸 ◆2VB8wsVUoo
2013/12/01(日) 16:39:21.37■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□
□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■
■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□
□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■
■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□
□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■
■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□
□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■
■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□
□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■
■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□
□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■
0599１３２人目の素数さん
2013/12/01(日) 18:23:00.30TTEST(data1;data2;,2;3)（不正確な書き方だが意味わかるよね）をやると0.2846271788になるんだが、
対応する計算をmathematicav9でやると、
TTest[{{1, 2, 3}, {2, 3, 5}}, 0, AlternativeHypothesis -> "Unequal"]（だよね？）が0.274577になるんだがこれであってるの？
ちなみにOpenOfficeでTTEST(data1;data2;,2;2)は0.2745766291になる
0600狸 ◆2VB8wsVUoo
2013/12/01(日) 18:38:48.86■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□
□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■
■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□
□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■
■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□
□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■
■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□
□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■
■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□
□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■
■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□
□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■□■
0601１３２人目の素数さん
2013/12/03(火) 17:20:18.74俺の環境ではディスプレイドライバをオフにすると正常な起動速度になる．
刺さってるカードがGeforceGTX680 (2枚) なのだが，CUDA関係で起動が重くなってるのかもしれない．
0602１３２人目の素数さん
2013/12/16(月) 14:11:33.50Windows8/Mathematica8 です
0603１３２人目の素数さん
2013/12/16(月) 19:54:20.68新規ノートブックのことだったら、知らん、てかどうでもいいだろ。
0604１３２人目の素数さん
2013/12/16(月) 20:23:51.10毎回使いやすい位置に移動するのが面倒な人もいるんだよ
そんなことも考え付かない想像力のなさでは数学者になるのは難しい
0605１３２人目の素数さん
2013/12/17(火) 03:13:20.270606１３２人目の素数さん
2013/12/17(火) 16:36:48.670607１３２人目の素数さん
2013/12/17(火) 19:00:27.880608１３２人目の素数さん
2013/12/17(火) 21:45:31.78いまはMacBookAirを持ち歩いているんで
0609１３２人目の素数さん
2013/12/17(火) 23:30:34.12年間使用料もあるな
個人で使うには馬鹿らしいので
もう、Raspberry PiとCDFの汎用化でいいや
cloud.sagemathがiPadでも使えて便利すぎ
ネット無しならiCASのREDUCEがかなり使える。
0610１３２人目の素数さん
2013/12/18(水) 16:11:46.83>iPad版
日本版だけ高価になるかも
0611１３２人目の素数さん
2013/12/18(水) 16:54:33.95>日本版だけ高価になるかも
間違いない。きっとそうだ。
日本語化なんてしなくていいってのに、無能な日本人？スタッフがごり押しするから、
結局日本では高くて売れなくなってんだよ。
0612１３２人目の素数さん
2013/12/18(水) 17:23:00.33ぶっちゃけmathは高機能電卓でワープロ。
工学には至極便利だが理学では突き詰めるとソース非公開ブラックボックスが問題だ。
また、商用の壁が大きい。
Sageの様な多数の商用システムとの統合は出来ない。
mathは十年程度で消滅の予感
UNIXがlinuxに侵食されるごとくに。
0613１３２人目の素数さん
2013/12/18(水) 17:52:46.03MatlabはOctaveでもScilabでもRでも同じだがMathematicaをMaximaで置き換えることはできない
有料ソフトで今後１０年生き残れる会社はAdobe SystemsとWolfram researchしかないだろう
それほどMathematicaの完成度は高い
0614ビートたけし
2013/12/18(水) 22:19:02.98英語でつかうにかぎるよ
0615１３２人目の素数さん
2013/12/19(木) 03:49:22.62どこがだよｗ
バグだらけで困ってるよ！
そのくせ問い合わせても内部アルゴリズムは公開できないとか言うからタチが悪い
matlabやmapleに比べてその思想には共感できるが、
もう少し完成度を上げてもらわないと困る
バージョン10に期待
いつでるんだろうね？いつもこの時期なのに・・・
0616１３２人目の素数さん
2013/12/19(木) 17:08:48.9210の機能確認ならRaspberry Piで出来るよ
mathの問題はバグった時、そこで袋小路
オプソなら自分で必要に応じて修正も拡張も出来るね
0617１３２人目の素数さん
2013/12/20(金) 01:08:21.81これで十分。
ver 4でもいいくらい。
0618１３２人目の素数さん
2013/12/20(金) 13:35:31.74想像力のなさでは数学者になるのは難しい
0619１３２人目の素数さん
2013/12/20(金) 19:24:57.680620１３２人目の素数さん
2014/01/02(木) 08:27:00.980621１３２人目の素数さん
2014/01/12(日) 19:17:16.43下記HPの計算を行うと、mathematica4.1では計算して、答えを表示しますが、mathematica9.0では計算しません。なぜでしょうか？
http://www.geocities.jp/dirac_equation/index6.htm
0622１３２人目の素数さん
2014/01/12(日) 19:27:16.89http://www.geocities.jp/dirac_equation/question.pdf
0623１３２人目の素数さん
2014/01/12(日) 20:34:35.390624１３２人目の素数さん
2014/01/13(月) 08:06:32.67http://www.geocities.jp/dirac_equation/index7.htm
0625１３２人目の素数さん
2014/01/13(月) 20:53:14.66やりたいことは、ｙ３に含まれるたとえば、「p0^2-p1^2-p2^2-p3^2」を 「m^2」
に変換して、短くまとめたいだけです。（やりたいことは、単純です）
（１） p0^2-p1^2-p2^2-p3^2;→変換→ m^2
（２） q0^2-q1^2-q2^2-q3^2;;→変換→ m^2
（３） k0^2-k1^2-k2^2-k3^2;;→変換→ 0
（４） j0^2-j1^2-j2^2-j3^2;;→変換→ 0
（５） p0*q0-p1*q1-p2*q2-p3*q3;;→変換→ m^2-t/2
（６） j0*p0-j1*p1-j2*p2-j3*p3;;→変換→ (1/2)*(m^2-u)
（７） j0*q0-j1*q1-j2*q2-j3*q3;;→変換→ (1/2)*(s-m^2)
（８） k0*j0-k1*j1-k2*j2-k3*j3;;→変換→ -(t/2)
（９） k0*p0-k1*p1-k2*p2-k3*p3;;→変換→ (1/2)*(s-m^2)
（１０） k0*q0-k1*q1-k2*q2-k3*q3;;→変換→ (1/2)*(m^2-u)*)
最終的に、
ｙ３＝-512 (40 m＾４+42 s u-43 m＾２(s+u))
という形にまとめたいです。
Mathematicaで、どのようコードを作れば良いでしょうか？
下記コードは、不安定で、mathematicaのバージョンにより計算したり、しなかったりします。
（mathematica4.1では計算して、9.0ではしません。mathematicaのバグではないでしょうか？）
http://www.geocities.jp/dirac_equation/index7.htm
0626１３２人目の素数さん
2014/01/13(月) 21:50:33.58・エネルギー運動量保存則を使う
・都合のいい座標系（例えば始／終状態の重心系）を取る
・そもそも成分にばらさずに計算する方法を考える
で、>>545は如何？
0627１３２人目の素数さん
2014/01/14(火) 20:51:38.40＞・そもそも成分にばらさずに計算する方法を考える
＞で、>>545は如何？
私にとっては、その方法の方が難しいです。その方法は、別途考えたいと思います。
今は、ここまで、コードを作ったので、現在の方法（成分にバラす方法）で、計算します。しかし、手詰まりです。理由は、Mathematica4.1では計算が出来て、9.0では出来ないからです。コードは、これで合っていて、計算するはずなのですが、、、、
式を変換させるだけの単純なことなのですが、どうしたら良いか？本当に困っています。
v =｛7, 8, 4, 10, 2, 5, 6, 1, 3, 9｝;とかを使っている理由は、（１）から（１０）の式を、ランダムに、変換させていただけです。
この部分は、元々、ランダムに、１から１０のリストを作っています。
ちょっと、今から、散歩に行って、いいアイデアが出ないか？歩きながら考えます。
0628１３２人目の素数さん
2014/01/15(水) 00:43:10.19>成分にバラす方法
だったら
>・エネルギー運動量保存則を使う
>・都合のいい座標系（例えば始／終状態の重心系）を取る
に従って
Simplify[y3 //. {p1 -> 0, p2 -> 0,
k0 -> p3, k1 -> 0, k2 -> 0, k3 -> -p3,
q0 -> p0, q1 -> p3*Sqrt[1 - z^2], q2 -> 0, q3 -> p3*z,
j0 -> p3, j1 -> -p3*Sqrt[1 - z^2], j2 -> 0, j3 -> -p3*z,
p0 -> (s + m^2)/(2 Sqrt[s]), p3 -> (s - m^2)/(2 Sqrt[s]),
z -> 1 + t/(2 p3^2), t -> 2 m^2 - s - u}]
とでもやれば
-512 (40 m^4 + 42 s u - 43 m^2 (s + u))
は出てくるけど。
0629１３２人目の素数さん
2014/01/15(水) 23:05:42.50MATLABの代わりになるでしょうか
0630１３２人目の素数さん
2014/01/16(木) 01:23:18.660631１３２人目の素数さん
2014/01/16(木) 06:48:26.730632１３２人目の素数さん
2014/01/16(木) 18:10:22.710633１３２人目の素数さん
2014/01/16(木) 21:32:40.04代理店からメールがきた。
アイコンは8,9と赤できてるから変えてほしいね。
俺が一番好きなアイコンは4だったか4.2だったかの
カラフルなやつ、めっちゃかっこよかった。
0634１３２人目の素数さん
2014/01/18(土) 12:54:16.61お返事有難う御座います。
＞-512 (40 m^4 + 42 s u - 43 m^2 (s + u))
＞は出てくるけど。
じぇじぇじぇじぇじぇじぇじぇじぇじぇじぇじぇじぇです。
すごいです。やっぱり世の中、賢い方がおられるのですね。
この方法は、スマートで、一発で、回答が得られます。
有難う御座いました。
申し訳御座いませんが、コンプトン散乱の最後の以下の計算につきましても、ご教示頂きましたら幸いです。
y1=(m^2*(w - w0)^2 + 2*m*w*(w - w0)*w0 + w*w0*(w^2 + w0^2))/(2*w0^4)
という式の
(w0-w)/(w0*w)
の部分を
1/m*(1-Cos[theta])
に変換したいのですが、どのようにすれば良いでしょうか？
下記の方法では、置き換えられません。
y1=(m^2*(w - w0)^2 + 2*m*w*(w - w0)*w0 + w*w0*(w^2 + w0^2))/(2*w0^4)/.(w0-w)/(w0*w)->1/m*(1-Cos[theta])
y2=(w^2*(-1 + w/w0 + w0/w + Cos[theta]^2))/(2*w0^2)
0635１３２人目の素数さん
2014/01/18(土) 21:17:43.34《下記の方法》はパターンがあっていません。むしろこの変換は
補助式による記号ｍの消去ではないでしょうか？
本来は式の組からの変数消去はEliminate[]ですが、手っ取り早く
書くと、
y1 /. Solve[(w0 - w)/(w0*w)
== 1/m*(1 - Cos[theta]), m] // Simplify
ですかね･･･
ただ w/w0の比で纏めるには、細工がいりますが、
簡単な整理ですのでマニュアルの方が早いです。
0636１３２人目の素数さん
2014/01/18(土) 22:42:25.08＞ただ w/w0の比で纏めるには、細工がいりますが、
＞簡単な整理ですのでマニュアルの方が早いです。
マニュアルを読んで纏めるようにします。
話が遡りますが、質問が中途半端でしたので、ご参考に、コンプトン散乱を計算した資料とコードを添付します。
ｐ１〜４は、「相対論的量子力学」ランダウから抜粋し、
ｐ５〜１０は、「輻射の量子論」ハイトラーから抜粋しました。
P9の図を描くため、ｐ９の式（４０）を導出したいです。
compton
の“ｃ”をクリック願います。
http://www.geocities.jp/dirac_equation/index8.htm
0637１３２人目の素数さん
2014/01/18(土) 23:39:01.55y1 /. Solve[(w0 - w)/(w0*w)
== 1/m*(1 - Cos[theta]), m] // Simplify // InputForm
は、
{(w*(w^2 - w*w0 + w0^2 + w*w0*Cos[theta]^2))/(2*w0^3)}
だから･･･あとは手操作(マニュアル)でやてってことだ。。。
ランダウか、教程はリフシッツがフォローしてたけど、場の
古典論以降ずれてる、何故今更にハイトラーなのかも判らん。
0638１３２人目の素数さん
2014/01/19(日) 10:21:53.07＞だから･･･あとは手操作(マニュアル)でやてってことだ。。。
失礼しました。マニュアルを読んで、後は、オートで計算できるように
してください。と勘違いしました。文章を読み返すと、そのように
解釈できませんが、、、早とちりでした。
しかし、全自動計算をやりたいです。この部分だけ、手計算になるからです。
（５年以上前からの悲願の部分です。）
初期条件だけ入力すれば、アウトプットの図が出るようにしたいです。
要するに、8 (m^4− s u + m^2 (3 s + u)等の計算結果から、P9の第１０図の
コンプトン散乱の実験結果を、アウトプットとして出したいです。
添付の資料から、他の計算方法を使って、全自動で計算することは出来ない
でしょうか？
＞ランダウか、教程はリフシッツがフォローしてたけど、場の
＞古典論以降ずれてる、
ランダウの本に記載されていることは、古典論ではなく、立派な
「場の量子論」の計算だと思います。
最近の教科書に記載されている計算式と基本的には変わっていないと
理解しています。（ゲージ理論等は反映されていないですが）
コプトン散乱の計算について、「ずれている」部分をご教示願います。
＞何故今更にハイトラーなのかも判らん。
P9の第１０図のようなコンプトン散乱の実験結果が記載された資料が
ないからです。ご存じでしたら、ご教示頂きましたら幸いです。
http://oshiete.goo.ne.jp/qa/8289193.html
0639１３２人目の素数さん
2014/01/19(日) 11:50:21.94＞y1 /. Solve[(w0 - w)/(w0*w)
＞== 1/m*(1 - Cos[theta]), m] // Simplify
＞ですかね･･･
本日は、休みなので、落ち着いて、計算してみますと、ご教示頂きました方法で
マニュアル操作することもなく、全自動で計算できました。
助かりました。本当に有難う御座いました。
0640１３２人目の素数さん
2014/01/19(日) 18:16:01.35青→赤になるようにグラフを色付けしてるんですが、もっといい方法あります？
あとこの方法だとカラーバーの凡例がAutomaticで表示できないのですが、
こちらもいい方法ありますか。
0641１３２人目の素数さん
2014/01/21(火) 00:23:09.34昔使った記憶しか残っていない。
他のソフトでは計算できないようなベッセル関数を
含むめんどくさい関数の積分も簡単にできた。
ファイルやフォルダもプログラムコード使って操作が
出来てなかなかの優れものと感心した覚えがある。
でも、これだけは言える。Mathematicaは価格が高すぎる。
買う気にならない。
0642１３２人目の素数さん
2014/01/21(火) 16:49:26.68スチューデント版
買え
0643１３２人目の素数さん
2014/01/21(火) 17:51:09.390644１３２人目の素数さん
2014/01/26(日) 21:09:27.43図を、繰り返し表示するにはどうしたらよいでしょうか？
以下では、表示しませんが、以下のようにして、５回plotを表示したいです。plot以外のコマンドでも、結構です。
For[n=1,n5,n++,
Plot[Sin[x],{x,0,6Pi}]
];
0645１３２人目の素数さん
2014/01/26(日) 22:08:48.900646１３２人目の素数さん
2014/01/27(月) 05:30:31.69Print@Plot[Sin[x], {x, 0, 6 Pi}]
0647１３２人目の素数さん
2014/01/27(月) 19:33:59.51ご回答有難う御座います。
よく解りました。この方法を使って、表示するようにします。
0648１３２人目の素数さん
2014/02/07(金) 11:29:43.02Spikeyちゃんだよ
http://mathworld.wolfram.com/images/gifs/SpikeyVersions.gif
http://mathworld.wolfram.com/images/eps-gif/Spiky_700.gif
http://mathworld.wolfram.com/Spikey.html
0649１３２人目の素数さん
2014/02/10(月) 23:21:19.09いちいちプログラムをすべて打ち込むの面倒なんだが
Excelとかもっと簡易的な計算ツールがあるけど
mathematicaの他に比べても優れている点、魅力って何よ？
0650１３２人目の素数さん
2014/02/11(火) 10:54:20.92きほんてきにじょうほうはなげないよ
0651１３２人目の素数さん
2014/02/11(火) 10:58:57.71ウリナラきげんをやられるとむかつくよ
■ このスレッドは過去ログ倉庫に格納されています