トップページ⇒math
1001コメント517KB

〓 Mathematica 伍 〓

レス数が1000を超えています。これ以上書き込みはできません。
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/
0951１３２人目の素数さん2014/11/27(木) 23:07:47.44
>>950
FullSimplify[(m^2)^(3/2)/m^3] はどうなるかをを考えてみ
0952１３２人目の素数さん2014/11/28(金) 07:17:52.06
＞＞951
有難う御座います。考えます。
0953１３２人目の素数さん2014/11/28(金) 21:41:47.33
お返事有難う御座います。
＞FullSimplify[(m^2)^(3/2)/m^3] はどうなるかをを考えてみ
ボキャブラリーというか、ExpnadAllで展開させることしか、思い浮かばないです。
まず、このFullSimplify[(m^2)^(3/2)/m^3]を、どうしたら、“１”になるのか？わからないです。
±があるので、正確には、これ以上どうしょうも、ないのでしょうか？
（ルートｍ＾２の解は、＋ｍと-ｍの２つが考えられる。）
このコンプトン散乱の場合、±の-は、一応、無視したいです。
その場合でも、どうしたら、良いのか？わからないです。
0954１３２人目の素数さん2014/11/28(金) 21:44:09.26
不細工ですが
√ｍ＾２→ｍ　に、力ずくで、置き換えれば良いかもしれませんね。
0955１３２人目の素数さん2014/11/28(金) 23:29:22.97
>>953
高校数学やり直し
0956１３２人目の素数さん2014/11/29(土) 10:50:53.66
ワロタ
0957１３２人目の素数さん2014/11/29(土) 10:51:45.15
FullSimplify[(m^2)^(3/2)/m^3, m > 0]
0958１３２人目の素数さん2014/11/29(土) 12:19:06.14
こんにちは、
ご回答有難う御座います。
今度は、図に出来ません。
原因は、虚数だからでしょうか？
（しかし、式全体に“Ｉ”を掛けても、図にできません。）
目盛りが、一致してないからでしょうか？

http://www.geocities.jp/dirac_equation/index1.html
0959１３２人目の素数さん2014/11/29(土) 13:22:18.28
>>958
この計算の答えは、存在しない（間違っている）
ということでしょうか
0960１３２人目の素数さん2014/11/30(日) 20:52:35.08
>>958
>>955
「複素数」って知ってる？
0961１３２人目の素数さん2014/11/30(日) 21:48:59.39
そうか。考え付かなかった。複素共役から、それなりのグラフ化が可能かも、しれませんね。
0962１３２人目の素数さん2014/11/30(日) 23:18:21.15
>>959
Mathematicaが出力したものが「この計算の答え」じゃないなら何なんだ
>>945によれば、そもそも入力が物理的にはデタラメだということは理解してるんだろ？
0963１３２人目の素数さん2014/12/01(月) 00:12:38.31
ノート見てないから適当言うけど、m って質量じゃないの？
0964１３２人目の素数さん2014/12/23(火) 15:31:03.39
こんにちは
乱数で表示された数のうち，1つだけ数えたいのですが・・・
1回1回実行するのが面倒なのでなんとかして，1度に複数回こなせるようなものを考えたいんですけどできません。
どうすればいいんでしょうか・・・？
0965１３２人目の素数さん2014/12/23(火) 19:33:03.97
>>964
>乱数で表示された数のうち，1つだけ数えたい
意味不明。それを1回実行するMathematiaコードを書いてみて。
0966１３２人目の素数さん2014/12/23(火) 22:21:46.85
このスレはコンプトンをやってる人に個人的に教えてあげるスレですか？
0967１３２人目の素数さん2014/12/24(水) 10:00:12.32
>>965
>意味不明
そりゃそうですよね・・・
すいませんでした。
一回実行するコードはこんな感じです。

dice010402 = Table[Random[Integer, {1, 4}],{22}]
dice010602 = Table[Random[Integer, {1, 6}], {3}]
count1[x_] := Which[x == 1, 1, x != 1, 0]
Sum[count1[dice010402[[k]]], {k, 1, 22}]
Sum[count1[dice010602[[k]]], {k, 1, 3}]

これを実行すると一回分の結果しか出ませんが，それを一度に複数回こなせるようなものに改良したいんです。
どうすればいいのでしょうか？
0968１３２人目の素数さん2014/12/24(水) 17:42:18.10
>>967
こういったかんじ？
Tally@RandomInteger[{1, 4}, 22]
Tally@RandomInteger[{1, 6}, 3]
0969１３２人目の素数さん2014/12/25(木) 13:58:44.35ID:U2/O3mDM
>>968
そんなかんじです
それを改良したいんですけど…
0970１３２人目の素数さん2014/12/25(木) 14:59:44.26ID:y7AcYpT7
>>969
改良とは、例題での1の度数だけでなく、総ての目の度数を
しらべることではないの？

968ではそうしてるけど･･･(出力フォーマットは別として)
0971１３２人目の素数さん2014/12/30(火) 15:20:41.68ID:yc9gAwbH
Jane Styleから書き込みえできません。
0972１３２人目の素数さん2014/12/30(火) 15:23:19.58ID:yc9gAwbH
>>971
書き込み出来ました。

こんにちは、
以下を教えて下さい。
１．ガンマ行列のトレース計算をしました。
http://www.geocities.jp/dirac_equation/index1.html

この計算で正しいのでしょうか？
ps={p0,-I*p1,-I*p2,-I*p3};
qs={q0,-I*q1,-I*q2,-I*q3};
ks={k0,-I*k1,-I*k2,-I*k3};
js={j0,-I*j1,-I*j2,-I*j3};
のように、空間成分に-Iが付くのでしょうか？

２．電子電子散乱の計算で、
http://www.geocities.jp/dirac_equation/landau.pdf
の式（82.6）の上の２つの式が、式（82.6）に一致することを計算しました。
http://www.geocities.jp/dirac_equation/index1.html

２つの式の内、
y01 = 2/t^2*((ps.ks)*(ps.ks) + (ps.js)*(ps.js) + 2*m^2*(m^2 - (ps.qs)))
y02 = 2/(t*u)*((ps.ks) - 2*m^2) *(ps.ks)
ｙ02は、式（82.6）に一致しますが、ｙ01は、一致しません。
どこが間違っているのでしょうか？
0973１３２人目の素数さん2014/12/30(火) 21:38:48.59ID:J3Yd0htO
特殊相対論の勉強を頑に拒むのはなぜ？
0974１３２人目の素数さん2015/02/11(水) 16:44:01.60ID:Yq7k8Ekq
てｓ
0975１３２人目の素数さん2015/02/11(水) 16:46:02.28ID:Yq7k8Ekq
てｓ
0976１３２人目の素数さん2015/02/11(水) 16:58:23.97ID:Yq7k8Ekq
てｓ
0977１３２人目の素数さん2015/02/17(火) 10:37:44.40ID:jGAH4sK+
1つのノートブック内に3つのセルを作成しました．

セルA：　計算パラメータの定義を集めたもの．
セルB：　計算パラメータの変更に依存しない部分の計算を集めたもの．
セルC：　計算パラメータに依存する計算を集めたもの．

この状況で，各セルを適宜再実行しながら計算するコードをセルDに書きたいと考えています．
そこで「あるセルから他のセルを呼び出して実行」ができるとよいのですが，可能でしょうか？
0978１３２人目の素数さん2015/02/17(火) 17:35:27.68ID:EdiRDBQD
>>977
Block とか Module とかじゃいかんの？
0979１３２人目の素数さん2015/02/18(水) 18:18:02.68ID:Af217Zem
>>978
うまく出来そうです．
ありがとうございます．
0980１３２人目の素数さん2015/02/23(月) 02:56:52.76ID:Jtdu54/Z
Needs["Benchmarking`"];
BenchmarkReport[]

1.33だった
0981１３２人目の素数さん2015/02/23(月) 18:50:29.06ID:x7brCvuk
>>980
System: Microsoft Windows (64-bit)
Date: February 23, 2015
Mathematica Version: 10.0.2
Benchmark Result: 1.95
0982１３２人目の素数さん2015/02/23(月) 19:05:54.69ID:Jtdu54/Z
どんだけいいマシン使ってんだよｗ
比較用に列挙されてるやつの一番いいやつでも1.89じゃん

マシンのスペック書いて！
最近は定格で4GHz超えてるマシンも出てきたみたいやね
0983１３２人目の素数さん2015/02/23(月) 22:26:55.60ID:x7brCvuk
>>982
CPU: Intel(R) Core(TM) i7-3970X CPU @ 3.50GHz
RAM: 64.0 GB
0984１３２人目の素数さん2015/02/23(月) 22:31:26.42ID:x7brCvuk
CPU: Intel(R) Xeon(R) CPU E5-2690 v2 @ 3.00GHz
RAM: 256 GB

System: Microsoft Windows (64-bit)
Date: February 23, 2015
Mathematica Version: 10.0.1
Benchmark Result: 1.41



明らかに得手不得手が出てるな
0985１３２人目の素数さん2015/02/23(月) 23:56:40.43ID:2NtpOD86
>>983
自分のは
Intel Core i7
Processor Speed: 2.7 GHz
RAM 16GB
だけどなんで64GBもメモリを積めるんですか？
0986１３２人目の素数さん2015/02/24(火) 00:22:09.26ID:1T2TwbB0
>>985
ﾏｻﾞﾎﾞが違うと思われ
俺が使ってるやつ↓
http://www.asus.com/jp/Motherboards/P9X79/
0987１３２人目の素数さん2015/02/24(火) 02:02:33.70ID:XH+kD/gH
メモリ256GBは初めて聞いた
CPUクロック数の上昇もそろそろ限界に達してるから
世の中のHPCの動向がそうであるように
MathematicaもそろそろGPUへの対応をした方がいいね
CUDAコード出力するだけじゃなくて組み込みレベルで。
0988１３２人目の素数さん2015/02/24(火) 20:16:28.04ID:9BvsOY5B
>>984
Mathematicaのデフォは4coreがmaxでなかったっけ？、
だとしたら10coreに上げてもclockを下げた分不利…

core extension導入でcore数を増やせるらしいけど、
ベンチマーク上は微増らしい、つか明示的に並列ｐ
書かんと速度上がらんぽい(コア管理がユーサー側

これじゃMathematicaを使う意味ないじゃんと想う
0989１３２人目の素数さん2015/02/25(水) 00:49:54.91ID:4QJ9DHPX
一部GPU対応しているはずだけどね
ちなみに近頃のHPC
ttp://www.hpc.co.jp/HPC7000-XI8192R8P_features.html
0990１３２人目の素数さん2015/02/25(水) 09:02:28.35ID:GacQMgDm
三年百六日。
0991１３２人目の素数さん2015/02/25(水) 20:51:50.18ID:GA/9W6kV
DateDifference[{2011, 11, 11}, DateValue[{"Year", "Month", "Day"}]]

->1202 days(*平均1ﾚｽ/day以下*)
0992１３２人目の素数さん2015/02/26(木) 00:18:27.64ID:UjfUV2xI
今はMathematica10だからスレッドもpart10になるぐらいのスピードにしないとな
1バージョン1スレッド
0993１３２人目の素数さん2015/02/26(木) 02:32:45.39ID:UjfUV2xI
このスレにはいないと思うが、これできる人いたらお願い。
http://wc2014.2ch.net/test/read.cgi/math/1424876583/
0994１３２人目の素数さん2015/02/26(木) 09:02:39.00ID:GgAkER3I
三年百七日。
0995１３２人目の素数さん2015/02/26(木) 09:54:15.69ID:xhFlfLHc
>>納n=1,m](n*C(m,n)*P(m,n)*P(l-m,m-n)/P(l,m))=m^2/l

微妙な処はあるが、Mathematica返しだと、
Simplify[Sum[
n Binomial[m, n] FactorialPower[m, n] FactorialPower[l - m,
m - n], {n, 1, m}]/FactorialPower[l, m] // FunctionExpand,
Assumptions -> Element[m, Integers]]

->m^2/l

だが？
0996１３２人目の素数さん2015/02/26(木) 17:44:44.33ID:GaHfmGDi
>>994
こういう階乗を含む級数は
Racahの公式 (http://journals.aps.org/pr/abstract/10.1103/PhysRev.62.438, (112)式) を使えば計算できる
0997１３２人目の素数さん2015/02/26(木) 23:09:26.80ID:UjfUV2xI
>>995
おお！サンクス！
まさかMathematicaにこの変形ができるとは思ってなかったわ！
どんなプロセスを経てるのかはブラックボックスだが。
4色定理の証明にコンピュータプログラムが認められたのは
ソースがオープンで内容について議論できたからだろうね。
クローズドソースなMathematicaについては証明とはちょっと呼べない。
もちろん、等式の正しさが限りなく100%に近づいたけどね。

>>996
リンク先みてみたんだけど、登録が必要で見れなかった。
ちょっと調べてみたけど、量子物理学系の話かな？
数学の公式のようなものは見当たらなかった。





本スレで「ヴァンデルモンドの恒等式」で証明できるというアドバイスを
もらったので、ちょっと考えてみるつもり。
0998１３２人目の素数さん2015/02/27(金) 00:18:28.42ID:EbhLAWFK
Racahの公式とヴァンデルモンドの恒等式は同じもの
これを少し変形すれば>>993の等式は証明できる
0999１３２人目の素数さん2015/02/27(金) 10:38:50.75ID:EbhLAWFK
梅
1000１３２人目の素数さん2015/02/27(金) 10:39:25.89ID:EbhLAWFK
〓 Mathematica 七 〓 [転載禁止](c)2ch.net
http://wc2014.2ch.net/test/read.cgi/math/1425001084/
1001１００１Over 1000Thread
このスレッドは１０００を超えました。
もう書けないので、新しいスレッドを立ててくださいです。。。
レス数が1000を超えています。これ以上書き込みはできません。