Question

用Mathematica做的动态图，怎么出不来结果呢

用Mathematica做的动态图，怎么出不来结果呢
t2 = ((E^(2 I k + I k1) (1 - 2 E^(I k) Cos[k]) -
E^(I k1) (E^(4 I k) - 2 E^(3 I k) Cos[k])) (E^(
2 I k1 (-1 + L) + I k1 L) -
E^(I k1 (-1 + L) + 2 I k1 L) \[Alpha] -
E^(2 I k1 (-1 + L) + I k1 L) \[Beta] -
2 E^(I k1 (-1 + L) + 2 I k1 L) Cos[k] +
E^(2 I k1 L) (-E^(I k1 L) + E^(I k1 (-1 + L)) \[Alpha] +
E^(I k1 L) \[Beta] + 2 E^(I k1 (-1 + L)) Cos[k])))/((-E^(
I k1 L + I k (2 + L)) +
2 E^(I k1 L + I k (1 + L)) Cos[k]) (-E^(
4 I k + I k1) - (E^(4 I k) - 2 E^(3 I k) Cos[k]) (-E^(I k) +
E^(I k + I k1) \[Alpha] + E^(I k) \[Beta] +
2 E^(I k + I k1) Cos[k])) (E^(2 I k1 (-1 + L) + I k1 L) -
E^(I k1 (-1 + L) + 2 I k1 L) \[Alpha] -
E^(2 I k1 (-1 + L) + I k1 L) \[Beta] -
2 E^(I k1 (-1 + L) + 2 I k1 L) Cos[k] +
E^(2 I k1 L) (-E^(I k1 L) + E^(I k1 (-1 + L)) \[Alpha] +
E^(I k1 L) \[Beta] + 2 E^(I k1 (-1 + L)) Cos[k])) + (E^(
I k1 (-1 + L) + I k1 L +
I k (1 + L)) - (-E^(I k1 L) + E^(I k1 (-1 + L)) \[Alpha] +
E^(I k1 L) \[Beta] +
2 E^(I k1 (-1 + L)) Cos[k]) (-E^(I k1 L + I k (2 + L)) +
2 E^(I k1 L + I k (1 + L)) Cos[k])) (E^(
4 I k + I k1) + (E^(4 I k) -
2 E^(3 I k) Cos[k]) (-E^(I k + 2 I k1) +
E^(I k + I k1) \[Alpha] + E^(I k + 2 I k1) \[Beta] +
2 E^(I k + I k1) Cos[k]) +
E^(2 I k1 L) (-E^(
4 I k +
I k1) - (E^(4 I k) - 2 E^(3 I k) Cos[k]) (-E^(I k) +
E^(I k + I k1) \[Alpha] + E^(I k) \[Beta] +
2 E^(I k + I k1) Cos[k]))));
FullSimplify[t2]
\[Beta] = 1 - Sqrt[2]; \[Alpha] = 1;
k1 = ArcCos[(\[Alpha] + 2*Cos[k])/(2*Sqrt[2])];
Animate[Plot[Abs[FullSimplify[t2]], {k, 0, \[Pi]}], {L, 1, 10, 1}]

Answer 1

代码修改一下，那个 t2 定义成函数最好，另外既然是画图，就不需要进行化简了，很费时间，去掉 FullSimplify。
t2[k_, L_] := ((E^(2 I k + I k1) (1 - 2 E^(I k) Cos[k]) -
E^(I k1) (E^(4 I k) -
2 E^(3 I k) Cos[k])) (E^(2 I k1 (-1 + L) + I k1 L) -
E^(I k1 (-1 + L) + 2 I k1 L) \[Alpha] -
E^(2 I k1 (-1 + L) + I k1 L) \[Beta] -
2 E^(I k1 (-1 + L) + 2 I k1 L) Cos[k] +
E^(2 I k1 L) (-E^(I k1 L) + E^(I k1 (-1 + L)) \[Alpha] +
E^(I k1 L) \[Beta] +
2 E^(I k1 (-1 + L)) Cos[k])))/((-E^(I k1 L + I k (2 + L)) +
2 E^(I k1 L + I k (1 + L)) Cos[
k]) (-E^(4 I k + I k1) - (E^(4 I k) -
2 E^(3 I k) Cos[k]) (-E^(I k) + E^(I k + I k1) \[Alpha] +
E^(I k) \[Beta] +
2 E^(I k + I k1) Cos[k])) (E^(2 I k1 (-1 + L) + I k1 L) -
E^(I k1 (-1 + L) + 2 I k1 L) \[Alpha] -
E^(2 I k1 (-1 + L) + I k1 L) \[Beta] -
2 E^(I k1 (-1 + L) + 2 I k1 L) Cos[k] +
E^(2 I k1 L) (-E^(I k1 L) + E^(I k1 (-1 + L)) \[Alpha] +
E^(I k1 L) \[Beta] +
2 E^(I k1 (-1 + L)) Cos[k])) + (E^(I k1 (-1 + L) + I k1 L +
I k (1 + L)) - (-E^(I k1 L) + E^(I k1 (-1 + L)) \[Alpha] +
E^(I k1 L) \[Beta] +
2 E^(I k1 (-1 + L)) Cos[k]) (-E^(I k1 L + I k (2 + L)) +
2 E^(I k1 L + I k (1 + L)) Cos[k])) (E^(4 I k +
I k1) + (E^(4 I k) -
2 E^(3 I k) Cos[k]) (-E^(I k + 2 I k1) +
E^(I k + I k1) \[Alpha] + E^(I k + 2 I k1) \[Beta] +
2 E^(I k + I k1) Cos[k]) +
E^(2 I k1 L) (-E^(4 I k + I k1) - (E^(4 I k) -
2 E^(3 I k) Cos[k]) (-E^(I k) +
E^(I k + I k1) \[Alpha] + E^(I k) \[Beta] +
2 E^(I k + I k1) Cos[k]))));
\[Beta] = 1 - Sqrt[2]; \[Alpha] = 1;
k1 = ArcCos[(\[Alpha] + 2*Cos[k])/(2*Sqrt[2])];
Animate[Plot[Abs[t2[k, L]], {k, 0, \[Pi]}], {L, 1, 10, 1}]