〓 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/
0554１３２人目の素数さん
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.82■ このスレッドは過去ログ倉庫に格納されています