〓 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/
0497１３２人目の素数さん
2013/07/15(月) NY:AN:NY.ANできるならそのとおりにやればいい
0498１３２人目の素数さん
2013/07/15(月) NY:AN:NY.ANID:cF5MhH6Y!0499１３２人目の素数さん
2013/07/15(月) NY:AN:NY.ANかな。
0500496
2013/07/15(月) NY:AN:NY.AN>>499
すみません、左辺も右辺も、どう書けばいいのかわかりません・・・。
フォルダ内にtime1.datとtime2.datとtime3.datがあり、それぞれをリストの形で読み込みたい、という質問です。
つまり
t1=Import["フォルダ名\\time1.dat,"Table""];
t2=Import["フォルダ名\\time2.dat,"Table""];
t3=Import["フォルダ名\\time3.dat,"Table""];
という作業を（実際はもっと多くの連番ファイルがあるので）Ｆｏｒ文で簡潔におこないたいのです。
教えてくださっているのは、右辺の部分でしょうか？
0501１３２人目の素数さん
2013/07/15(月) NY:AN:NY.ANt[i]=Import["フォルダ名\\time"<>i<>".dat","Table"]
]
ていう感じか？動くかどうかは確認していないが。
0502１３２人目の素数さん
2013/07/15(月) NY:AN:NY.AN○For[i=1,i<=3,i++
0503496
2013/07/15(月) NY:AN:NY.ANありがとうございます！
>>502を踏まえて、
For[i=1,i<=3,i++,
t[i]=Import["フォルダ名\\time"<>i<>".dat","Table"]
]
とやってみたり、i++のあとの,を;にしてみたりしたのですが、
StringJoin::string: String expected at position 2 in フォルダ名\time<>1<>.dat. >>
Import::chtype: First argument フォルダ名\time<>1<>.dat is not a valid file, directory, or URL specification. >>
というエラーが出てしまいました。
また、右辺に関してもt[1]やt[2]・・・というように配列に入れていくのではなく、t1やt2という名前の変数として定義したいのですがどうしたらよいのでしょうか。（t0095のように、tのあとに0を付けることもできるようにしたいためです）
0504１３２人目の素数さん
2013/07/16(火) NY:AN:NY.AN0505１３２人目の素数さん
2013/07/16(火) NY:AN:NY.AN0506１３２人目の素数さん
2013/07/16(火) NY:AN:NY.ANまず、499の書いている通り、
<>i<>
でなく
<>ToString[i]<>
でないとだめだと思う。
次に、たとえば、t00095とかにしたいなら
NumberForm[95, 5, NumberPadding -> {"0", ""}]が使えると思う。
よくわからないがTostringしないと"t"とつながらなかった。
最後に
ToExpression["a345=123"]とかつかえば、
a345に123が代入される。
0507１３２人目の素数さん
2013/07/16(火) NY:AN:NY.AN>Tostringしないと"t"とつながらなかった。
Headをみればわかりますよ。
NumberForm[95, 5, NumberPadding -> {"0", ""}]//Head
->NumberForm(*これでは<>:StringJoinは機能しない*)
ToString@NumberForm[95, 5, NumberPadding -> {"0", ""}]//Head
->String(*これなら<>:StringJoinは機能できる*)
0508１３２人目の素数さん
2013/07/16(火) NY:AN:NY.ANドキュメントセンターのグラフのこと？
事前にデフォルトオプションを変更して、
SetOptions[Plot, Exclusions -> None];
Plot[SawtoothWave[x], {x, -2, 2}]
とするか、個別にオプションを指定して
、
Plot[SawtoothWave[x], {x, -2, 2}, Exclusions -> None]
としないとつながらないよ。
0509505
2013/07/16(火) NY:AN:NY.AN0510１３２人目の素数さん
2013/07/17(水) NY:AN:NY.AN当時のやっつけの記述ですので、洗練はされてません。
(*filePath = "DataFolderにPassが通っていないときに記述"*);
fileName =
Table[
StringJoin["AgtABCD_image", StringTake["0" <> ToString[n], -2], ".gif"],
{n, 0, 63}(*データ数:64*)
];
imageList = Map[Import[#(*, Path -> filePath*)] &, fileName];
(*fileNameにリテラル値をいれる場合の例*)
literalName =
Table[
StringJoin["literal_image", StringTake["0" <> ToString[n], -2], ".gif"],
{n, 0, 63}
];
fileName = literalName;
※StringTake["0" <> ToString[n], -2]部は、今だとNumberFormの方がスマート？かも。
0511１３２人目の素数さん
2013/08/17(土) NY:AN:NY.AN0512１３２人目の素数さん
2013/08/25(日) NY:AN:NY.AN0513１３２人目の素数さん
2013/08/25(日) NY:AN:NY.ANWolframもアドビみたいに
ライセンスで金とるだけか
0514１３２人目の素数さん
2013/09/15(日) 19:19:21.980515計算
2013/09/25(水) 22:26:18.76約４時間程度と思います。
自分のPCの速度を見極めたいです。
http://www.geocities.jp/dirac_equation/index3.htm
0516１３２人目の素数さん
2013/09/26(木) 20:54:18.390517計算
2013/09/26(木) 21:03:18.141000秒程度のはずです。
http://www.geocities.jp/dirac_equation/index4.htm
0518計算
2013/09/26(木) 21:26:45.13下記計算を、昨日の２３時ころからやっておりますが、未だに終わりません。
PDF参照
本来なら、約６０秒×２５６／（６０秒×６０分）＝４時間 程度で終わるはずです。
なぜでしょうか？
http://www.geocities.jp/dirac_equation/index3.htm
の計算結果です。
http://www.geocities.jp/dirac_equation/time.pdf
0519１３２人目の素数さん
2013/09/26(木) 21:37:05.960520計算
2013/09/26(木) 22:02:30.42メモリー不足ですと、下記のように62.343750 秒と表示されても、実際はもっと計算時間がかかっているのでしょうか？
x= 0,y= 0
62.343750, Null <
0521１３２人目の素数さん
2013/09/27(金) 15:50:14.81ガンマ行列の性質を使って整理してから g[4〜15]=0 に注意すれば Mathematica に頼るまでもない計算のようだが。
0522計算
2013/09/27(金) 23:31:26.45そうでしょうか？このガンマ行列は、通常の４×４ではありません。256×256が１６個あります。またg[4〜15]=0は、ここで、特別に示しただけです。０でないケースを計算して、どうなるか？を検討するのが目的です。
ところで、
やっぱり、おかしいです。
下記計算を、一昨日の２３時ころからやっておりますが、未だに終わりません。
本来なら、約６０秒×２５６／（６０秒×６０分）＝４時間 程度で終わるはずです。
なぜでしょうか？
http://www.geocities.jp/dirac_equation/index3.htm
0523１３２人目の素数さん
2013/09/27(金) 23:59:17.47D次元（今はD=16）の場合
・γ^μγ^ν+γ^νγ^μ
・γ^μγ_μ
・tr(γ^μγ^ν)
と、それから
・γ^νγ^μγ_ν
・tr(γ^μγ^νγ^ργ^σ)
がそれぞれどうなるか御存知だろうか？それらを使えば、件のトレースはずっと簡単に表される。
>ガンマ行列の性質を使って整理
とはそういうこと。
0524計算
2013/09/28(土) 00:16:20.99手計算はメンドクサイし、間違いも発生します。
すでに、コードは完成しておりますので、 Mathematica に頼ります。
このコードを完成させるのに、むちゃ苦労しました。
0525１３２人目の素数さん
2013/09/28(土) 00:26:35.49「完成している」コードの動作が「おかしい」とはこれ如何に
0526計算
2013/09/28(土) 00:27:21.12256×256の１６個のガンマ行列を使って、下記を得るのが目標です。
256×256の４個のガンマ行列なら、簡単です。
http://www.geocities.jp/dirac_equation/fig.pdf
0527計算
2013/09/28(土) 00:31:26.98そこを質問しているのです。でも多分メモリー不足だと思います。
１つの結果を出すのに、１週間も時間がかかるのに、８GBのメモリーですから、、
0528１３２人目の素数さん
2013/09/28(土) 00:46:13.960529１３２人目の素数さん
2013/09/28(土) 02:05:15.77「コードは完成してる(はずだ)」けど、「正しく実行されてない」と疑ってるってことでしょ。
原因はマシン性能の可能性があるから、もっとパワーのあるマシンで実行してみて欲しいと。
0530137
2013/09/28(土) 09:17:35.90質問１
y = {2, 5, 1, 3, 4}
のような １〜５の整数の５個の乱数のリストを作るには、どのようなコマンドを使用すれば良いでしょうか？
質問２
z = x + x^2 + 4 r + t^2
の場合、式の項数は、４個ですが、どのようなコマンドを使用すれば４個と表示するでしょうか？
0531137
2013/09/28(土) 10:03:09.770532137
2013/09/28(土) 10:47:59.04繰り返し計算をｎ回させ、Xに１から１０の乱数を発生させて、その乱数が２以下の場合に、ｘ＝ｙとして、その値ｙを表示させて、繰り返し計算をそこで、ストップさせることは可能でしょうか？
下記は、失敗例です。どこを修正すれば良いでしょうか？
For[n = 1, n <= 10, n++,
v = RandomSample[Range[10]];
x = v[[n]];
Print["n=", n, ",x=", x];
If[x < 1, y = x & Print[y]];
];
0533137
2013/09/28(土) 13:24:27.210534計算
2013/09/28(土) 13:32:05.57もっと速く計算する方法はないでしょうか？
例えば、マッキントッシュ用のmathematicとパソコンを使用すれば、速く計算可能でしょうか？その場合、どのくらい速くなるのでしょうか？
その他の方法はないでしょうか？
（※；但し、メモリーは８GBで不足してます。購入中です。）
http://www.geocities.jp/dirac_equation/index6.htm
0535１３２人目の素数さん
2013/09/28(土) 20:23:26.69q0とかも、256×256の行列でしょうか？
やっていることは、256×256の行列をいくつか用意して、
それらの間で、何らかの単純な計算をしているということでいいのでしょうか？
逆行列を求めるとかでない限り、そんなに時間がかかるプログラムにはならないと思います。
Cとか、自分でメモリの管理ができる言語で作り直すほうが、手っ取り早いのでは？
0536計算
2013/09/28(土) 20:48:51.29例えば、C言語で、２回微分させるコードを書くのに、５０行必要としました。
Mathematicaなら、１行で済みます。
この計算は、一般相対論の計算より、かなり複雑なので、C言語で書き直すのは、難しすぎて、出来ないです。
メモリの管理なのですが、マックパソコンでは、Mathematicaのメモリー管理が出来るのではないでしょうか？また、windowsよりも、速いのではないでしょうか？
0537１３２人目の素数さん
2013/09/28(土) 21:25:51.91速さはWindowとかMacとかによるものではない
0538計算
2013/09/29(日) 09:30:05.17Windows 8 - 128GB
Windows 8 Professional - 512GB
Windows 8 Enterprise - 512GB
まで、増設できるので、windowsPCのメモリー増設で対応できそうな気がします。
http://bto-pc.jp/btopc-com/repair/ram-max-beginner-pc.html
0539１３２人目の素数さん
2013/09/29(日) 17:31:05.37まず>>523を考えた方がいい。
0540計算
2013/09/29(日) 18:27:47.85当然、523の計算はやっております。ただし、それはここでアップしたコードの計算後です。
このやり方は、教科書の記載通りです。この方法で、４×４ の４個のガンマ行列を使い、以下の正しい答えを導いてます。
http://www.geocities.jp/dirac_equation/fig.pdf
現在、１６GBに増設して、計算中です。
これで、メモリー不足になって計算できない場合は、以下のPCの購入も検討するつもりです。
http://www.dospara.co.jp/5shopping/detail_prime/2/56/3208/0/
0541１３２人目の素数さん
2013/09/29(日) 20:52:48.60>当然、523の計算はやっております。
それを知ってて
http://www.geocities.jp/dirac_equation/index6.htm
のxとyの和をいちいち成分で足すとかどういうマゾ趣味か、とは思うなぁ
0542１３２人目の素数さん
2013/09/29(日) 21:41:14.200543１３２人目の素数さん
2013/09/30(月) 09:10:34.45例えば Σ[k=1,n]k=n(n+1)/2 を知ってたら 1+2+3+4+…+99999+100000 を一つずつ足し算するようなことはせんだろ。計算氏がやってるのはそういうことだ。
0544計算
2013/10/01(火) 21:32:41.36そうかもしれませんが、この方法で、
http://www.geocities.jp/dirac_equation/fig.pdf
の正しい結果を導出しております。
但し、ご指摘通り、遠回りをしている可能性はあります。しかし、この方法は教科書通りです。
ところで、計算ができない理由が解りました。
２５６回の繰り返し計算で、項数が多い式を
y1 = y1 + s1;
で、作って、後で
Print[FullSimplify[ExpandAll[y1]]];
で、一気に整理しようとしたから、メモリーがパンクして計算できませんでした。
１回の繰り返し計算ごとに
Y1=FullSimplify[ExpandAll[y1+s1]];
で、整理するように、すれば、計算できるはずです。（只今、計算中）
0545１３２人目の素数さん
2013/10/01(火) 22:44:55.40もしかして、
>D次元（今はD=16）の場合
>・γ^μγ^ν+γ^νγ^μ
>・γ^μγ_μ
>・tr(γ^μγ^ν)
>と、それから
>・γ^νγ^μγ_ν
>・tr(γ^μγ^νγ^ργ^σ)
>がそれぞれどうなるか
の結果を知らないのか？知ってたら書いてみて。
0546計算
2013/10/01(火) 22:56:02.51y1 = 1024*(-64*k2^2*m^2 - 64*k3^2*m^2 + 64*m^4 - 48*m^2*p0^2 + 48*m^2*p1^2 - 16*k2*m^2*p2 + 48*m^2*p2^2 - 16*k3*m^2*p3 + 48*m^2*p3^2 + 49*k2^2*p0*q0 + 49*k3^2*p0*q0 - 63*m^2*p0*q0 + 49*p0^3*q0 - 49*p0*p1^2*q0 -
49*p0*p2^2*q0 - 49*p0*p3^2*q0 - 49*k2^2*p1*q1 - 49*k3^2*p1*q1 + 63*m^2*p1*q1 - 49*p0^2*p1*q1 + 49*p1^3*q1 + 49*p1*p2^2*q1 + 49*p1*p3^2*q1 + 112*k2*m^2*q2 - 98*k2*p0^2*q2 + 98*k2*p1^2*q2 + 49*k2^2*p2*q2 -
49*k3^2*p2*q2 + 63*m^2*p2*q2 - 49*p0^2*p2*q2 + 49*p1^2*p2*q2 + 98*k2*p2^2*q2 + 49*p2^3*q2 + 98*k2*k3*p3*q2 + 98*k2*p3^2*q2 + 49*p2*p3^2*q2 +
7*(7*k3^2*p3 + p3*(-7*k2^2 + 9*m^2 + 7*(-p0^2 + p1^2 + p2^2 + p3^2)) + 2*k3*(8*m^2 + 7*(-p0^2 + p1^2 + p2*(k2 + p2) + p3^2)))*q3 + k1^2*(-64*m^2 + 49*(p0*q0 + p1*q1 - p2*q2 - p3*q3)) +
k0^2*(64*m^2 + 49*(p0*q0 + p1*q1 + p2*q2 + p3*q3)) + 2*k1*(-8*m^2*(p1 - 7*q1) + 49*((-p0^2)*q1 + p1^2*q1 + (p2*(k2 + p2) + p3*(k3 + p3))*q1 + p1*(k2*q2 + k3*q3))) -
2*k0*(-8*m^2*(p0 - 7*q0) + 49*((-p0^2)*q0 + (p1*(k1 + p1) + p2*(k2 + p2) + p3*(k3 + p3))*q0 + p0*(k1*q1 + k2*q2 + k3*q3))));
纏めると
y1=-24064m^4+26112m^2s+25088m^2u-25088su;
になります。
0547１３２人目の素数さん
2013/10/10(木) 06:53:26.430548１３２人目の素数さん
2013/10/10(木) 20:08:33.09, '´ _. -‐'''"二ニﾆ=-`ヽ、
／ ／:::::; -‐''" `ｰノ
/ ／:::::／ ＼
/ /::::::/ | | | |
| |:::::/ / | | | | | |
| |::/ / / | | || | | ,ﾊ .| ,ﾊ|
| |/ / / /| ,ﾊﾉ| /|ﾉﾚ,ﾆ|ﾙ'
| | | / / ﾚ',二､ﾚ′ ,ｨｲ|ﾞ/
. | ＼ ∠イ ,ｲイ| ,`-' |
| l＾,人| ` `-' ゝ | このスレは馬と鹿と豚さんばかりね。
| ` -'＼ ｰ' 人
| /(l ＿＿／ ヽ、
| （:::::`‐-､__ |::::`､ ﾋﾆﾆヽ、
| ／ `‐-､::::::::::`‐-､::::＼ /,ﾆﾆ､＼
| |::::::::::::::::::|` -､:::::::,ﾍ￣|'､ ﾋﾆ二、 ＼
. | /::::::::::::::::::|::::::::＼/:::Ｏ`､::＼ | '､ ＼
| /:::::::::::::::::::/:::::::::::::::::::::::::::::'､::::＼ﾉ ヽ、 |
| |:::::／:::::::::/:::::::::::::::::::::::::::::::::::'､',::::'､ /:＼__/‐､
| |／:::::::::::/::::::::::::::::::::::::::::::::::Ｏ::| '､::| く:::::::::::::￣|
| /_..-'´￣`ｰ-､:::::::::::::::::::::::::::::::::::|/:/`‐'::＼;;;;;;;＿|
| |／::::::::::::::::::::::＼:::::::::::::::::::::::::::::|::/::::|::::/:::::::::::/
| /:::::::::::::::::::::::::::::::::|:::::::::::::::::::::Ｏ::|::|::::::|:::::::::::::::/
0549１３２人目の素数さん
2013/10/11(金) 10:22:32.09例えば A, B, C を互いに非可換としてその他の数や変数は可換扱いにしたいです．
このとき
Simplify[10*A*20*B-10*B*20*A]
のような計算があったとして，
結果が 200 (A**B-B**A) になるようにしたいのですが…
0550549
2013/10/11(金) 13:02:22.60http://math.ucsd.edu/~ncalg/
0551１３２人目の素数さん
2013/11/08(金) 09:12:06.47例えば｛a,b,c,d,e}（a〜eは数値）といったリストとA（数値）があって、
a,b,c,d,eによる線形結合で、ただし係数は整数、Aに一番近い値を求めたいです。
このときA＝１a+２b+３c+４d+５e
みたいにどういう形なのかも具体的に知りたいです
こういうこ機能のある関数はありますか？
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.82■ このスレッドは過去ログ倉庫に格納されています