In[1]:= IdentityTr[x_] := InverseFourier[Fourier[x]]
In[2]:= a = {0, 0, 0.33333333333333333, 0}
b = {0, 0, 0.333333333333333333, 0}
Precision[a[[3]]]
Precision[b[[3]]]
Out[2]= {0, 0, 0.333333, 0}
Out[3]= {0, 0, 0.33333333333333333, 0}
Out[4]= MachinePrecision
Out[5]= 17.5229
In[8]:= NestList[IdentityTr, a, 20]
NestList[IdentityTr, b, 20]
Out[8]= {{0, 0, 0.333333, 0}, {0., 0., 0.333333, 0.}, {0., 0.,
0.333333, 0.}, {0., 0., 0.333333, 0.}, {0., 0., 0.333333, 0.}, {0.,
0., 0.333333, 0.}, {0., 0., 0.333333, 0.}, {0., 0., 0.333333,
0.}, {0., 0., 0.333333, 0.}, {0., 0., 0.333333, 0.}, {0., 0.,
0.333333, 0.}, {0., 0., 0.333333, 0.}, {0., 0., 0.333333, 0.}, {0.,
0., 0.333333, 0.}, {0., 0., 0.333333, 0.}, {0., 0., 0.333333,
0.}, {0., 0., 0.333333, 0.}, {0., 0., 0.333333, 0.}, {0., 0.,
0.333333, 0.}, {0., 0., 0.333333, 0.}, {0., 0., 0.333333, 0.}}
Out[9]= {{0, 0, 0.33333333333333333, 0}, {0.*10^-18,
0.*10^-18 + 0.*10^-18 I, 0.33333333333333333,
0.*10^-18 + 0.*10^-18 I}, {0.*10^-17, 0.*10^-18 + 0.*10^-18 I,
0.3333333333333333, 0.*10^-18 + 0.*10^-18 I}, {0.*10^-17,
0.*10^-17 + 0.*10^-17 I, 0.333333333333333,
0.*10^-17 + 0.*10^-17 I}, {0.*10^-16, 0.*10^-16 + 0.*10^-16 I,
0.333333333333333, 0.*10^-16 + 0.*10^-16 I}, {0.*10^-15,
0.*10^-15 + 0.*10^-15 I, 0.33333333333333,
0.*10^-15 + 0.*10^-15 I}, {0.*10^-14, 0.*10^-15 + 0.*10^-15 I,
0.3333333333333, 0.*10^-15 + 0.*10^-15 I}, {0.*10^-14,
0.*10^-14 + 0.*10^-14 I, 0.333333333333,
0.*10^-14 + 0.*10^-14 I}, {0.*10^-13, 0.*10^-13 + 0.*10^-13 I,
0.33333333333, 0.*10^-13 + 0.*10^-13 I}, {0.*10^-12,
0.*10^-12 + 0.*10^-12 I, 0.33333333333,
0.*10^-12 + 0.*10^-12 I}, {0.*10^-11, 0.*10^-11 + 0.*10^-11 I,
0.3333333333, 0.*10^-11 + 0.*10^-11 I}, {0.*10^-10,
0.*10^-11 + 0.*10^-11 I, 0.333333333,
0.*10^-11 + 0.*10^-11 I}, {0.*10^-10, 0.*10^-10 + 0.*10^-10 I,
0.33333333, 0.*10^-10 + 0.*10^-10 I}, {0.*10^-9,
0.*10^-9 + 0.*10^-9 I, 0.33333333,
0.*10^-9 + 0.*10^-9 I}, {0.*10^-8, 0.*10^-8 + 0.*10^-8 I, 0.3333333,
0.*10^-8 + 0.*10^-8 I}, {0.*10^-7, 0.*10^-8 + 0.*10^-8 I, 0.333333,
0.*10^-8 + 0.*10^-8 I}, {0.*10^-7, 0.*10^-7 + 0.*10^-7 I, 0.33333,
0.*10^-7 + 0.*10^-7 I}, {0.*10^-6, 0.*10^-6 + 0.*10^-6 I, 0.3333,
0.*10^-6 + 0.*10^-6 I}, {0.*10^-5, 0.*10^-5 + 0.*10^-5 I, 0.3333,
0.*10^-5 + 0.*10^-5 I}, {0.*10^-4, 0.*10^-4 + 0.*10^-4 I, 0.333,
0.*10^-4 + 0.*10^-4 I}, {0.*10^-3, 0.*10^-4 + 0.*10^-4 I, 0.33,
0.*10^-4 + 0.*10^-4 I}}
2013年11月16日土曜日
Mathematica Fourier変換について
バグを疑い始めたので,プログラムを純化させてみる.
登録:
コメントの投稿 (Atom)
0 件のコメント:
コメントを投稿