〓 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/
0682１３２人目の素数さん
2014/02/20(木) 00:59:02.600683１３２人目の素数さん
2014/02/20(木) 01:33:38.71Download版なので何もない
メールだけ
0684１３２人目の素数さん
2014/02/20(木) 14:14:51.34Fedexをはじめて受け取ったのがMathematicaのバージョンアップだった
0685１３２人目の素数さん
2014/02/20(木) 20:45:43.16伝統もオワコンになるかもシレン
0686１３２人目の素数さん
2014/02/24(月) 22:18:31.31ここに、８個の物があるとします。これを、A,B,C,Dの４種類に、対称性があるように、並べたいのですが、mathematicaのコマンドを使って、並べることは可能でしょうか？
例として
{A,A,A,A,A,A,A,A}
{A,B,A,B,A,B,A,B}
{A,A,A,A,B,B,B,B}
{A,A,B,B,C,C,D,D}
などです。
{B,D,B,C,A,C,A,A}
など、規則性がなく、でたらめに並べるのは、駄目です。
A,B,C,Dでなくても、１，２，３，４でも、青、赤、黄、緑の４種類でも結構です。
0687１３２人目の素数さん
2014/02/24(月) 22:41:36.52普通は全部列挙して「対称性」があるヤツだけ抜き出せばよろし
0688１３２人目の素数さん
2014/02/25(火) 20:30:12.68厳密には、そうかもしれません。
対称性があれば、なんでも、構いませんので、mathematicaで、並べるコマンドがあれば、教えてください。
＞普通は全部列挙して「対称性」があるヤツだけ抜き出せばよろし
それはそうですが、何卒、お願いします。
実は、８個ではなく、１６個で、やりたいのです。問題例を簡単にするため、８個としました。
0689１３２人目の素数さん
2014/02/25(火) 22:14:47.69Sort
0690１３２人目の素数さん
2014/02/25(火) 22:37:51.70>{A,B,A,B,A,B,A,B}
これを対称性があるというのか?
0691１３２人目の素数さん
2014/02/25(火) 23:30:03.92それを逆順にしてコピーしたのと２つ並べれば全部できるけどな．
0692１３２人目の素数さん
2014/02/26(水) 05:30:49.52691を応用して
左右や、上下や回転の関数を作って
1x16から4x4までの行列(ベクトル)に適用して
それをリストにして、重複するやつを排除したら？
0693１３２人目の素数さん
2014/02/26(水) 19:59:57.46イメージは、解りました。
ご教示頂いた方法で、作るかないですか？
たとえばSubsets[Range とかに、似たようなコマンドで、
作れないでしょうか？
0694１３２人目の素数さん
2014/02/26(水) 20:26:49.86まず「対称性がある」と「対称性がない」の識別条件を不足なく全部列挙してみ。
Mathematicaへの実装を考えるのはそれから。
0695１３２人目の素数さん
2014/02/26(水) 20:49:15.010696１３２人目の素数さん
2014/02/26(水) 21:32:26.770697１３２人目の素数さん
2014/02/26(水) 22:09:26.57しらんがな...
Tuples
0698１３２人目の素数さん
2014/02/27(木) 13:16:11.66「ゆとり」かね
0699１３２人目の素数さん
2014/02/27(木) 17:27:52.78>ここに、８個の物があるとします。これを、A,B,C,Dの４種類に、対称性があるように、並べたいのですが、mathematicaのコマンドを使って、並べることは可能でしょうか？
>例として
>{A,A,A,A,A,A,A,A}
>{A,B,A,B,A,B,A,B}
>{A,A,A,A,B,B,B,B}
>{A,A,B,B,C,C,D,D}
>などです。
まず
こいつのいう対称性というのがなあ
対称性という言葉を覚えたてなのか
0700１３２人目の素数さん
2014/02/27(木) 17:33:16.57if (俺の中で対称) {
選択;
} else {
捨てる;
}
}
このアルゴリズムを書くだけだから本人には簡単じゃん
他人には無理だけど
0701１３２人目の素数さん
2014/02/27(木) 18:14:10.71出題者の期待は完全に裏切ってるけどw
0702１３２人目の素数さん
2014/02/27(木) 20:11:12.24うそでしょう。
>このアルゴリズムを書くだけだから本人には簡単じゃん
たぶん、出来ません。
理由は、次のとおり、
＞ if (俺の中で対称) {
頭の中では、想像できても、実際にコードを考えると、定義が難しくて、出来ないでしょう。
＞for (全パターン) {
１６個となると、メモリーがパンクするでしょう。
0703１３２人目の素数さん
2014/02/28(金) 14:19:12.09メモリーが増え続けるはずが無いプログラムのはずなんだが、Mathematicaの
（ mathkernelじゃない）使用メモリーがどんどん増えて行き、長く表示させておくと、
必ずクラッシュする。メモリー増加のせいでクラッシュすると思うんだが、
どうにかならないかな？
0704１３２人目の素数さん
2014/02/28(金) 18:49:48.950705703
2014/02/28(金) 20:11:06.89MathKernelの使用メモリーが増えていないんだから、
プログラムではメモリー消費が増加しづけていない、と考えていいんだよね？
0706703
2014/02/28(金) 20:21:21.77でも Mathematicaの使用メモリーの増加が尋常じゃない、、、。
Magitianというユーティリティでメモリー開放すると、ぐっと減る、、、。
放っておいても mathematicaが不使用メモリーを開放してくれるのかな、、、。
でもクラッシュするからな、、、。
0707１３２人目の素数さん
2014/02/28(金) 21:02:30.58CDF Playerで動かすとどうかな?
0708１３２人目の素数さん
2014/02/28(金) 21:44:09.92結局、「出来ない。」の結論ですね。
>{A,B,A,B,A,B,A,B}
は、ローカル対称性です。
対称性は、{B,D,B,C,A,C,A,A}など、規則性がなく、でたらめに並べるのは、NGで、規則性があれば、「対称性あり」と定義します。
＞コマンド一発でなんとかなると思うのは、やっぱり＞「ゆとり」かね
そこまで、言うなら、創造性あるプログラムが書けますか？匿名で、能書きたれるだけって、空し過ぎませんか？→言うからには、責任を持って、プログラム書いてください！無責任な発言に、本気で怒っています！
0709１３２人目の素数さん
2014/02/28(金) 22:08:36.61> 規則性があれば、「対称性あり」と定義します。
{3,1,4,1,1,1,2,2}
{1,3,1,4,1,3,1,3}
{2,3,4,4,2,2,2,4}
{3,3,4,3,2,4,1,1}
…
は、ある（乱数を使わない）規則に従って生成されたものだが、「対称性あり」なのだね？
0710１３２人目の素数さん
2014/02/28(金) 22:13:06.00誤：{3,3,4,3,2,4,1,1}
正：{3,3,4,3,2,3,1,1}
だった。
0711１３２人目の素数さん
2014/02/28(金) 22:15:36.80他人にアルゴリズムを実装させようという発想について
0712１３２人目の素数さん
2014/02/28(金) 23:12:01.70>ローカル対称性
でたらめな名前を付けるな
自分勝手な対称性だろ
0713１３２人目の素数さん
2014/03/01(土) 06:58:15.38普通、プログラムで使用する乱数ですら規則に基づいて作られる。
君の言う規則が何をさすのか分からない以上、答えようがない。
0714１３２人目の素数さん
2014/03/01(土) 09:47:03.45なんかのアイデアが出るかもしれないのにな。
0715１３２人目の素数さん
2014/03/01(土) 10:25:55.67できない理由を、「対称性の定義が不明確であること。」にしていると思います。
明らかに、問題のすり替えです。
その証拠に、「対称性を、あなたが定義して、16個の物を、A,B,C,Dの４種類に、
mathematicaのコマンドを使って、並べてください。」という質問に変更しても、
あなたは、プログラムを作ることは、できないでしょう。
たぶん、出来ない理由を、再度、問題のすり替えによって、回答するだけでしょう。
0716１３２人目の素数さん
2014/03/01(土) 10:59:25.430717１３２人目の素数さん
2014/03/01(土) 11:40:48.68問題としては >>691,692 で十分でしょ
どこかよその掲示板の
もうちょっと実務的、実用的な言語で質問して
あとでMathematicaで書き直せばいいじゃん
ここで質問しても、もうまともには相手にされないよ
0718１３２人目の素数さん
2014/03/01(土) 13:54:14.27そもそも「対称性が自発的に破れる」という発想が必要なんじゃないかな。
0719１３２人目の素数さん
2014/03/02(日) 00:58:22.26>できない理由を、「対称性の定義が不明確であること。」にしていると思います。
「>>686が用いている『対称性』という用語の定義を>>686自身が明示していないこと」が原因。
「私が見て対称性があると思う」では数学的な定義になっていない。
>>708の
>規則性があれば、「対称性あり」と定義します。
も「私が見て規則性があると思う」でしかなく、定義になっていない。
とにかく>>694に答えてみてよ。
0720１３２人目の素数さん
2014/03/02(日) 01:24:22.83ところで、ホントにやりたいことは何？
0721１３２人目の素数さん
2014/03/02(日) 02:11:08.15対称性とやらが定義されていないのに，どうやって実装するのか逆に聞きたい．
君のこれまでの発言の中の「対称性」って言葉を，
全て「うんこ」という言葉に置き換えたら，
このままでは実装できないことに気づくんじゃない？
ちゃんと定義されていれば「うんこ」でも実装できるだろうよ．
0722１３２人目の素数さん
2014/03/02(日) 09:31:45.98サウンドの機能を使って、「ビンゴ」と、音を出すことは可能でしょうか？
0723１３２人目の素数さん
2014/03/02(日) 10:26:00.880724１３２人目の素数さん
2014/03/02(日) 10:38:37.47Speakは？
さもなきゃSystemDialogInput["RecordSound"]
で録音して使えば？
しかし俺のは日本語版なんだが、 Speak["ビンゴ"]としてもしゃべらない！
日本語化されてない！
0725１３２人目の素数さん
2014/03/02(日) 11:01:46.58Speak[bingo]で、なんとか、しゃべってくれました。
0726１３２人目の素数さん
2014/03/15(土) 10:43:35.80下記HPのPDFは、電子の対消滅の計算式です。
（素粒子物理 牧二郎先生、林浩一先生著P１３４，１３５抜粋）
この計算を、下記HPのmathematica eμ計算の通り、作成したのですが、エラーが発生して、答が出ません。（1/0で無限になります。）
どこを修正すれば、良いでしょうか？
http://www.geocities.jp/dirac_equation/index10.htm
0727１３２人目の素数さん
2014/03/15(土) 14:39:01.880728１３２人目の素数さん
2014/03/15(土) 18:15:57.380729１３２人目の素数さん
2014/03/16(日) 21:24:18.45Array[a,{4,4}]/IdentityMatrix[4] と
Array[a,{4,4}]/Array[b,{4,4}] を試せば
0730１３２人目の素数さん
2014/03/21(金) 19:53:42.72＞Array[a,{4,4}]/IdentityMatrix[4] と
＞Array[a,{4,4}]/Array[b,{4,4}] を試せば
試しました。前者は、1/0となり、エラーになります。
HPの計算も結局、IdentityMatrix[4]で、割っています。
IdentityMatrix[4] を、Array[b,{4,4}]で置き換えて、後で、
Array[b,{4,4}]を、IdentityMatrix[4] に、更に置き換えれば
計算できるのでしょうか？
0731１３２人目の素数さん
2014/03/22(土) 15:40:21.53>Array[a,{4,4}]/Array[b,{4,4}]
の結果を観察すれば、どうして
>Array[a,{4,4}]/IdentityMatrix[4]
でエラーが起きるのかわかるでしょ。それで、貴方がやりたいのはそういう計算なのか、ということが問題。
0732１３２人目の素数さん
2014/03/22(土) 17:12:45.41T1 = Simplify[
y1 //. {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}];
の修正が必要です。Mをどこに追加すれば、良いでしょか？
0733１３２人目の素数さん
2014/03/22(土) 22:20:54.40そういう質問が出るということは、貴方は自分が何を計算しているのか分かっていないということだ。
0734１３２人目の素数さん
2014/03/29(土) 10:38:45.82http://www.geocities.jp/dirac_equation/index10.htm
に再計算した結果をアップしておりますが、正しい答が得られません。
下記式の意味を理解してないからです。
この式の意味は、どのような本、資料を見れば、解るでしょうか？
T1 = Simplify[
y1 //. {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}];
0735１３２人目の素数さん
2014/03/29(土) 23:14:40.83持ってる本に書いてあるだろうが、まず
p0, p1, p2, p3, q0, q1, q2, q3, k0, k1, k2, k3, j0, j1, j2, j3
は各々どういう量なのか説明してみて
0736１３２人目の素数さん
2014/03/30(日) 22:03:52.32＞は各々どういう量なのか説明してみて
p0, p1, p2, p3,は、e(-)の4元運動量で、p0がエネルギー で、p1, p2, p3（空間成分）が運動量 です。
q0, q1, q2, q3は、e(+)の4元運動量、k0, k1, k2, k3は、μ(-)の4元運動量、j0, j1, j2, j3は、μ(+)の4元運動量です。
ｍは電子e、Mはμの質量です。
0737１３２人目の素数さん
2014/03/31(月) 00:43:20.45それがわかってるなら
>>626の
>・都合のいい座標系（例えば始／終状態の重心系）を取る
に従って適切な座標系を設定して運動量の成分の値を決めるだけ
0738１３２人目の素数さん
2014/03/31(月) 20:07:01.41＞に従って適切な座標系を設定して運動量の成分の値を決めるだけ
考えます。以下は、教えて頂いた式を並べたものです。
p0 -> (s + m^2)/(2 Sqrt[s]),
p1 -> 0,
p2 -> 0,
p3 -> (s - m^2)/(2 Sqrt[s]),
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,
z -> 1 + t/(2 p3^2), t -> 2 m^2 - s - u}];
（１） 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)
0739１３２人目の素数さん
2014/04/05(土) 15:12:47.46重心系、実験室系、マンデルシュタム変数に関して記述している文献が少なくて、難儀しております。
p0 -> (s + m^2)/(2 Sqrt[s]),
p1 -> 0,
p2 -> 0,
p3 -> (s - m^2)/(2 Sqrt[s]),
---------------------------------------
q0 -> p0,
q1 -> p3*Sqrt[1 - z^2],
q2 -> 0,
q3 -> p3*z,
---------------------------------------
k0 -> p3,
k1 -> 0,
k2 -> 0,
k3 -> -p3,
---------------------------------------
j0 -> p3,
j1 -> -p3*Sqrt[1 - z^2],
j2 -> 0,
j3 -> -p3*z,
---------------------------------------
z -> 1 + t/(2 p3^2),
t -> 2 m^2 - s - u}];
は、下記HPの式（6.26）から導出できるのでしょうか？
http://osksn2.hep.sci.osaka-u.ac.jp/~naga/kogi/konan-class04/ch6-qed.pdf#search='%E5%AE%9F%E9%A8%93%E5%AE%A4%E7%B3%BB+%E9%87%8D%E5%BF%83%E7%B3%BB+%E3%83%9E%E3%83%B3%E3%83%87%E3%83%AB%E3%82%B7%E3%83%A5%E3%82%BF%E3%83%A0'
0740１３２人目の素数さん
2014/04/07(月) 03:44:03.14>・都合のいい座標系（例えば始／終状態の重心系）を取る
について、試したことを書いてみて
0741１３２人目の素数さん
2014/04/07(月) 07:32:46.87Mathematicaブックはver4までだね。
4.2？の追加が出たのが、日本語での製本版の最後だったかと。
>>681
ペーパーレスもそうだけど、市販で1万前後をバンドルするのは、
それなりのコストと考えたんだろうな。
あと、Helpの形式が今みたいなドキュメントセンター？に変わったのが、
次のver5からだと記憶してるから、それも理由だろうな。
>>724
画像を読み込んでのOCR？のTextRecognizeなんて機能もあるけど、
日本語はサンプル以外、使えたもんじゃないよなｗ
ちょっと前からgeocitiesで示してる人いるけど、
もしかしてMathematica cafeで締め出された人じゃねーの？ｗ
0742１３２人目の素数さん
2014/04/07(月) 07:51:39.02俺が知ってる範囲で教えよう。
あくまで、推察だけどなｗ
代理店がボッてるというのは誤解。国内での販売価格はW社が決めてる。
納入価格は代理店によって違うかもしれないけど。
おそらく、W社としては日本じゃ大学orキャンパスむけみたいな
大規模なライセンスを維持、拡大できればいいと思ってる程度なんだろう。
比較的新しい、ファイナンスとモデラーが日本語化されない、
ネタにもならないというのが証左だろう。
個人でも買える単体ライセンス（HomeなりProなりStdなり？）が
日本で高いのは、日本語化のコストの問題。
それと、Mathematicaに限らないが、こういうソフトが主に使われる
アカデミックな場でのライセンス意識の低さを考慮してると思う。
0743１３２人目の素数さん
2014/04/07(月) 08:05:11.18Internationalでライセンスの区分があり、値段が違った（Domesticの方が安い）。
たぶん、EU諸国もUSより高かったはず。
アメリカで買った人が日本で新バージョンにアップグレードする時は、
アップグレード費用に加えて、USとInternationalの差額を追加する必要があった気がする。
ver3くらいから、日本語をふくめて各種言語のアドオン？も登場する。
その国では、現地言語をふくめたパッケージしか販売しないことで、
値段差への言い訳ができたんだろうな。
現在の価格差に関しては、現地語翻訳なりマーケティング？や、
先述のライセンス意識の低さによる設定だろうけど、
基本的にMathematicaは何処で買うか、なので日本でも英語版は買える。
多分、フランス語版さえも。
でも、「何処で買うのか」の違いなので、日本で英語版を注文しても、
Home Editionなら日本語版と同じ6万くらい？になる訳だ。
むしろ、インストーラに英語しか含まれてないステータスにされる。
日本語がサッパリな外人でもなければ、普通に日本語を買った方が手間がない。
オンラインアクティベーションは、俺の知る限り言語と住まいはチェックしてないはず。
0744１３２人目の素数さん
2014/04/08(火) 21:07:50.14ご回答有難う御座います。
＞>・都合のいい座標系（例えば始／終状態の重心系）を取る
＞について、試したことを書いてみて
何も試しておりません。重心系、実験室系、マンデルスタム係数に関する知識が足りなくて、考えられない状態です。
下記を読んでも理解できないです。また素粒子物理 牧二郎先生、林浩一先生著には、その記載がありません。
私には、まず重心系、実験室系、マンデルスタム係数に関する知識が必要なのですが、素粒子物理のどの本を見れば載っているでしょうか？
理解されている方から見れば、「何で解らんのか？」と思われるかもしれませんが、よろしくお願い致します。
http://osksn2.hep.sci.osaka-u.ac.jp/~naga/kogi/konan-class04/ch6-qed.pdf#search='%E5%AE%9F%E9%A8%93%E5%AE%A4%E7%B3%BB+%E9%87%8D%E5%BF%83%E7%B3%BB+%E3%83%9E%E3%83%B3%E3%83%87%E3%83%AB%E3%82%B7%E3%83%A5%E3%82%BF%E3%83%A0'
0745１３２人目の素数さん
2014/04/08(火) 21:54:27.02>重心系、実験室系、マンデルスタム係数に関する知識
重心系：扱っている粒子の（空間）運動量の合計がゼロであるような座標系
実験室系：↑であるように実験を設定すれば重心系と同じ
マンデルスタム変数：http://www.geocities.jp/dirac_equation/desnhishoumetu.pdf の(8.55)式
で終わり。
0746１３２人目の素数さん
2014/04/08(火) 22:21:29.12考えてみます。
ん、ん、ん、でも、この知識から、以下を導出するのは、難しくないですか？
p0 -> (s + m^2)/(2 Sqrt[s]),
p1 -> 0,
p2 -> 0,
p3 -> (s - m^2)/(2 Sqrt[s]),
、、、、、以下続く。
しかし、ゆっくり考えてみます。
0747１３２人目の素数さん
2014/04/09(水) 00:39:00.69変形べっする関数を含む積分を解いていたら、
解の中に BesselI^(1,0)[...] というのがでてきました。
この BesslI と引数の間の肩の部分は何を意味するのでしょうか？
0748１３２人目の素数さん
2014/04/09(水) 00:50:02.77偏導関数では
0749１３２人目の素数さん
2014/04/09(水) 01:03:27.53早速のお答えありがとうございます！
このシンタックスの意味をどう検索すればよいかさえわからず困り果てていたところでした。
本当に助かりました。
0750１３２人目の素数さん
2014/04/09(水) 02:08:20.260751１３２人目の素数さん
2014/04/10(木) 17:16:37.87また質問させてください。
Integrate で積分範囲を複素数全体に指定することはできないでしょうか。
実数であれば
Integrate[f[x],{x,-Infinity,Infinity}]
とするだけなわけですが、複素数も同じように手軽にできればと思っています。
現状では絶対値と偏角、あるいは実部と虚部に分けるなどの
変数変換を手動で書いて対応しているのですが、
これが省略できるとかなり作業が楽になるのです。
0752１３２人目の素数さん
2014/04/10(木) 18:09:51.84>Integrate[f[x],{x,-Infinity,Infinity}]
に比べて
Integrate[f[z]/.z->x+I*y,{x,-Infinity,Infinity},{y,-Infinity,Infinity}]
がとても面倒とも思わないけど、具体的にはどんな？
0753１３２人目の素数さん
2014/04/10(木) 19:26:39.11/. はありがたみが理解できないままスルーしてしまっていました。。。
今までやってきたことはというと、例えば
Integrate[Exp[-Abs[z1-z2]^2] /. z1->x1+I*y1, {x1,-Infinity,Infinity}, {y1,-Infinity,Infinity]
と書けば済む問題があれば、
x1=Re[z1], x2=Re[z2], y1=Im[z1], y2=Im[z2] ととらえて
Integrate[Exp[-(x1-x2)^2-(y1-y2)^2],{x1,-Infinity,Infinity},{y1,-Infinity,Infinity}, Assumptions->Element[x2,Reals]&&Element[y2,Reals]]
と書き直したり、極座標系に変数変換した形で
r1=Abs[x1], r2=Abs[x2], t=Arg[x1-x2] とみなして
Integrate[
r1 Exp[-(r1^2 + r2^2 - 2 r1 r2 Cos[t])^2], {t, -Pi, Pi}, {r2, 0,
Infinity}, Assumptions -> Element[r1, Reals]]
にしたりとか、やたらと手動で作業してました。
勉強不足でお恥ずかしい限りです。
0754１３２人目の素数さん
2014/04/26(土) 10:32:08.50(Showとdraw2d、FindRootとmnewton程度の)共通的な方法で
連想配列的にデータを渡すにはどうすればいいでしょうか？
いまは、とりあえず、データの固まりの定義(r1、r2)
データの指定(r)
受け渡し(/.、ev)
をこんな感じにしているのですが
r1 = r -> {a -> 1, b -> 2};
r2 = r -> {a -> 3, b -> 4};
r = r /. r1;
{a, b} /. r
r1:r=[a=1,b=2]$
r2:r=[a=3,b=4]$
r:ev(r,r1)$
ev([a,b],r);
0755１３２人目の素数さん
2014/04/27(日) 11:18:36.42PDF 2200ページにもなりますた！！！！
0756１３２人目の素数さん
2014/04/27(日) 23:27:47.419^9^9^9
やってみろ
0757１３２人目の素数さん
2014/04/28(月) 17:47:11.44, '´ _. -‐'''"二ニﾆ=-`ヽ、
／ ／:::::; -‐''" `ｰノ
/ ／:::::／ ＼
/ /::::::/ | | | |
| |:::::/ / | | | | | |
| |::/ / / | | || | | ,ﾊ .| ,ﾊ|
| |/ / / /| ,ﾊﾉ| /|ﾉﾚ,ﾆ|ﾙ'
| | | / / ﾚ',二､ﾚ′ ,ｨｲ|ﾞ/
. | ＼ ∠イ ,ｲイ| ,`-' |
| l＾,人| ` `-' ゝ |
| ` -'＼ ｰ' 人 私は死なないわよ。
| /(l ＿＿／ ヽ、 でも最近一寸太ったかしら。
| （:::::`‐-､__ |::::`､ ﾋﾆﾆヽ、 Windows ver.10 で
| ／ `‐-､::::::::::`‐-､::::＼ /,ﾆﾆ､＼ 元の痩せた姿にしてよね。
| |::::::::::::::::::|` -､:::::::,ﾍ￣|'､ ﾋﾆ二、 ＼
. | /::::::::::::::::::|::::::::＼/:::Ｏ`､::＼ | '､ ＼
| /:::::::::::::::::::/:::::::::::::::::::::::::::::'､::::＼ﾉ ヽ、 |
| |:::::／:::::::::/:::::::::::::::::::::::::::::::::::'､',::::'､ /:＼__/‐､
| |／:::::::::::/::::::::::::::::::::::::::::::::::Ｏ::| '､::| く:::::::::::::￣|
| /_..-'´￣`ｰ-､:::::::::::::::::::::::::::::::::::|/:/`‐'::＼;;;;;;;＿|
| |／::::::::::::::::::::::＼:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/
| /:::::::::::::::::::::::::::::::::|:::::::::::::::::::::Ｏ::|::|::::::|:::::::::::::::/
0758１３２人目の素数さん
2014/05/02(金) 18:41:01.62v='a
v::1
をMathematicaでできるでしょうか？
aが未定義のうちは
v=a、v=`a、v=Symbol["a"]
など(いずれか)として
Set[Evaluate[v],1]
でaに値が入りますが
一度aが定義されてしまうと
1=1
になってしまうためうまく行きません
0759１３２人目の素数さん
2014/05/02(金) 20:41:06.13maximaを知らないので何とも。どういう動作をさせたいのでしょうか？
0760１３２人目の素数さん
2014/05/02(金) 20:48:56.76/.を使っておこないます
例
x^2+2*x+1/.x->1
こたえ4
0761１３２人目の素数さん
2014/05/02(金) 21:52:00.16あまり意味のある例ではないと思いますが
Do[Set[Evaluate[v],{}],{v,{`a,`b,`c}}]
(記述の雰囲気だけです。動きません。)
と書いたら
{a,b,c}=Table[{},{3}]
となるような感じで
maximaの場合
for v in ['a,'b,'c] do (v::[]);
と書くと
[a,b,c]:makelist([],3);
になります
(758の2行目は v:'a でした)
0762１３２人目の素数さん
2014/05/02(金) 22:12:56.40>Do[Set[Evaluate[v],{}],{v,{`a,`b,`c}}]
をやったら
{a,b,c} は {{},{},{}} を返すんじゃないの？
0763１３２人目の素数さん
2014/05/02(金) 23:23:17.100764１３２人目の素数さん
2014/05/02(金) 23:39:02.74それが動くように見えるのは、
最初に書いた通りa,b,cが未定義の状態だけです
そのあとで
Do[Set[Evaluate[v],{1}],{v,{`a,`b,`c}}]
とやると失敗します
{{},{},{}}={{1},{1},{1}}
になるからです
0765１３２人目の素数さん
2014/05/04(日) 23:13:10.130766１３２人目の素数さん
2014/05/12(月) 23:33:12.320767１３２人目の素数さん
2014/05/13(火) 19:02:25.8510
が来るぞ
0768１３２人目の素数さん
2014/05/24(土) 21:18:49.41a,bは定数として、xだけが変数でグラフを書きたいんだけど
どうしたらいい？
上記のだとグラフが一切表示されないんだけど
0769１３２人目の素数さん
2014/05/24(土) 21:58:55.700770１３２人目の素数さん
2014/05/24(土) 22:02:28.27a, bを引数とする関数にすればいいんじゃないの？
それをg(a,b)として、a,bを外から与えれば
グラフがかけるという風にすればいいでしょ。
0771１３２人目の素数さん
2014/05/24(土) 22:30:07.77？
プログラミングにあまり詳しくないから
使う関数名と簡単な例を用いて教えてくれ
0772１３２人目の素数さん
2014/05/24(土) 22:44:19.730773１３２人目の素数さん
2014/05/25(日) 15:28:38.230774１３２人目の素数さん
2014/05/26(月) 10:16:46.34Manipulateが最適解と思うが…
768は釣りだから相手にするな。
0775１３２人目の素数さん
2014/05/26(月) 20:22:22.44釣りじゃねえよハゲ
0776１３２人目の素数さん
2014/05/26(月) 20:33:36.61report[board_] := Sow[board]
solutionQ[board_] := And[Not[MemberQ[board, 0, {2}]], report[board]; True]
candidates[board_, position_] := Complement[Range[9],
board[[position[[1]], ;;]],
board[[;; , position[[2]]]],
Flatten[Take[board, Sequence @@ (# - {2, 0} & /@ (3 Ceiling[position/3]))]]]
children[board_] := With[{position = First[Position[board, 0]]}, ReplacePart[board, position -> #] & /@ candidates[board, position]]
descendantSolvableQ[board_] := Or @@ (solvableQ /@ children[board])
solvableQ[board_] := Or[solutionQ[board], descendantSolvableQ[board]]
0777１３２人目の素数さん
2014/05/26(月) 20:37:28.97board = {
{1, 0, 0, 0, 0, 7, 0, 9, 0},
{0, 3, 0, 0, 2, 0, 0, 0, 8},
{0, 0, 9, 6, 0, 0, 5, 0, 0},
{0, 0, 5, 3, 0, 0, 9, 0, 0},
{0, 1, 0, 0, 8, 0, 0, 0, 2},
{6, 0, 0, 0, 0, 4, 0, 0, 0},
{3, 0, 0, 0, 0, 0, 0, 1, 0},
{0, 4, 0, 0, 0, 0, 0, 0, 7},
{0, 0, 7, 0, 0, 0, 3, 0, 0}};
AbsoluteTiming[output = Reap[solvableQ[board]]]
{1.070002, {True, {{{{1, 6, 2, 8, 5, 7, 4, 9, 3}, {5, 3, 4, 1, 2, 9,
6, 7, 8}, {7, 8, 9, 6, 4, 3, 5, 2, 1}, {4, 7, 5, 3, 1, 2, 9, 8,
6}, {9, 1, 3, 5, 8, 6, 7, 4, 2}, {6, 2, 8, 7, 9, 4, 1, 3,
5}, {3, 5, 6, 4, 7, 8, 2, 1, 9}, {2, 4, 1, 9, 3, 5, 8, 6,
7}, {8, 9, 7, 2, 6, 1, 3, 5, 4}}}}}}
TableForm[output[[2]]]
（ｒｙ
0778１３２人目の素数さん
2014/05/26(月) 20:40:19.990779１３２人目の素数さん
2014/05/28(水) 03:43:31.02まずPECE法が調べてもさっぱりなんだけども
予測子修正子法（PECE）でP(q,Dq)の値を与えた時
(Dq,DDq)を返すにはどうしたらいいかな？
0780１３２人目の素数さん
2014/05/28(水) 19:05:05.86ドラえもんではありません
0781狸 ◆2VB8wsVUoo
2014/05/29(木) 13:04:53.43>論理性が欠如していようがなんだろうが、芳雄は『後世に語り継がれる秀逸な結果』を残しただろうが
>なのに何が科学者の敵だ
>芳雄は過程はどうであれ結果を残した、お前は科学者をなめるんじゃない
>お前は芳雄を妨害して芳雄の研究成果に悪い影響を与えている、お前こそ研究者の敵だろうが
>今からでもいいから素直になって芳雄に謝ってこい、それぐらいはできるだろうが
>
■ このスレッドは過去ログ倉庫に格納されています