〓 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/
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ウリナラきげんをやられるとむかつくよ
0652１３２人目の素数さん
2014/02/11(火) 13:01:47.71下記HPには、３個のmathematicaプログラムが入っています。
Compton44は、４×４行列のγ行列を４個使って、コンプトン散乱を計算しています。
Compton164は、16×16行列のγ行列を４個使って、コンプトン散乱を計算しています。
この２つは、クライン仁科の式を導出しています。
さて問題は、Compton164です。これは、16×16行列のγ行列を８個使って、コンプトン散乱を計算しています。直接項と交換項は、導出できたのですが、干渉項の値が違っています。
16×16行列のγ行列の種類や係数等を、上手く使うと、クライン仁科の式を導出できると思うのですが、どのようなγ行列の種類や係数を使えば、良いでしょうか？
http://www.geocities.jp/dirac_equation/index8.htm
0653１３２人目の素数さん
2014/02/11(火) 13:04:48.57すいません。説明不足です。
ファイルは、
http://www.geocities.jp/dirac_equation/index8.htm
の中の
compton00
に入っています。
0654１３２人目の素数さん
2014/02/11(火) 15:27:26.10さて問題は、Compton168でした。
0655１３２人目の素数さん
2014/02/11(火) 17:31:29.38すいません。大きな間違いに、気づきました。
652〜654は、なかったことにしてください。
失礼しました。
0656１３２人目の素数さん
2014/02/11(火) 18:40:23.47すいません。やっぱり間違っていなかったです。
質問を書き直します。
下記HPには、３個のmathematicaプログラムが入っています。
Compton44は、４×４行列のγ行列を４個使って、コンプトン散乱を計算しています。
Compton164は、16×16行列のγ行列を４個使って、コンプトン散乱を計算しています。
この２つは、クライン仁科の式を導出しています。
問題は、Compton168です。これは、16×16行列のγ行列を８個使って、コンプトン散乱を計算しています。直接項と交換項は、導出できたのですが、干渉項の値が違っています。
16×16行列のγ行列の種類や係数等を、上手く使うと、クライン仁科の式を導出できると思うのですが、どのようなγ行列の種類や係数を使えば、良いでしょうか？
ファイルは、
http://www.geocities.jp/dirac_equation/index8.htm
の中の
compton00
に入っています。
0657１３２人目の素数さん
2014/02/11(火) 23:21:08.19そっちにちゃんと質問を書いて、答えも投稿できるようにして
お願いしますの一言とURLだけ貼れよ
0658１３２人目の素数さん
2014/02/12(水) 13:46:01.07質問への答は「違うべくして違うものを計算しているのだから一致しないのが正しい」
もはや板違いだが、まず
・何次元の時空での物理を問題にしているのか
・その時空では、電磁場の成分はいくつなのか、特に光子の偏極ベクトルについてはどうか
をはっきりさせて考察すべき。
0659１３２人目の素数さん
2014/02/12(水) 20:41:00.78＞質問への答は「違うべくして違うものを計算しているのだから一致しないのが正しい」
一般的には、そうかもしれません。しかし、実際に計算してみますと、１６×１６のγ行列８個で、直接項と交換項は、導出できました。多分、残りの干渉項も、導出できるような気がしています。
存在する、すべての１６×１６のγ行列８個で計算して、クライン仁科の式を導出できない場合は、納得します。
0660１３２人目の素数さん
2014/02/14(金) 01:00:31.52ver7をたま〜に使う標準的なユーザです。
ver9 にupdateしてまで使うメリットってありますか？
0661１３２人目の素数さん
2014/02/14(金) 03:57:00.53ver 5.2でも十分なんじゃねーの
0662１３２人目の素数さん
2014/02/14(金) 09:06:55.560663１３２人目の素数さん
2014/02/14(金) 13:55:12.330664１３２人目の素数さん
2014/02/14(金) 18:49:33.006以降を使うと戻れない
0665１３２人目の素数さん
2014/02/14(金) 22:21:26.894.0と6.0はどうちがうのですか？あなたが6.0を
持っているというだけですか？
自分のは4.5でした
■ このスレッドは過去ログ倉庫に格納されています