(Clear ky;
Clear t1;
Clear t2;
Clear \[CurlyPhi];
Clear \[CurlyPhi]1;
Clear \[CurlyPhi]2;
Clear k1;
Clear \[CurlyPhi]3;
Clear b;
Clear a;
Clear c;
Clear t;
f1 = "C:\\Users\\Administrator\\Desktop\\巨可爱\\\[Phi].xls";
OpenWrite[f1];
kongge = " ";
v1 := 3.0;
v2 := -3.0;
bb1 := 3.0;
bb2 := -3.0;
b1 := 0.50;
b2 := 0.50;
w1 := 0.50;
w2 := 0.50;
s0 := Sign[\[Epsilon]];
s1 := Sign[\[Epsilon] - v1];
s3 := Sign[\[Epsilon] - v2];
kx0 := \[Epsilon]*Cos[\[CurlyPhi]];
k0 := Sqrt[(\[Epsilon])^2 - ky^2];
k1 := Sqrt[((\[Epsilon] - v1))^2 - ky^2];
k2 := Sqrt[(\[Epsilon])^2 - (ky + bb1)^2];
k3 := Sqrt[((\[Epsilon] - v2))^2 - ky^2];
k4 := Sqrt[(\[Epsilon])^2 - (ky + bb2)^2];
\[CurlyPhi]1 := ArcSin[ky/Abs[\[Epsilon] - v1]];
\[CurlyPhi]2 := ArcSin[(ky + bb1)/\[Epsilon]];
\[CurlyPhi]3 := ArcSin[ky/Abs[\[Epsilon] - v2]];
\[CurlyPhi]4 := ArcSin[(ky + bb2)/\[Epsilon]];
\[Phi]1 :=
ArcTan[Tan[
k1*b1]*(s0*s1 -
Sin[\[CurlyPhi]]*Sin[\[CurlyPhi]1])/(Cos[\[CurlyPhi]]*
Cos[\[CurlyPhi]1])];
\[Phi]2 :=
ArcTan[Tan[
k2*w1]*(1 -
Sin[\[CurlyPhi]]*Sin[\[CurlyPhi]2])/(Cos[\[CurlyPhi]]*
Cos[\[CurlyPhi]2])];
\[Phi]3 :=
ArcTan[Tan[
k3*b2]*(s0*s3 -
Sin[\[CurlyPhi]]*Sin[\[CurlyPhi]3])/(Cos[\[CurlyPhi]]*
Cos[\[CurlyPhi]3])];
\[Phi]4 :=
ArcTan[Tan[
k4*w2]*(1 -
Sin[\[CurlyPhi]]*Sin[\[CurlyPhi]4])/(Cos[\[CurlyPhi]]*
Cos[\[CurlyPhi]4])];
mm = -1267;
de = 0.01;
\[Epsilon] := -1367.0*de;
\[CurlyPhi] := 10.0*3.1415926/180.0;
For[n = -1367,
n <= mm,
\[Epsilon] = n*de,
(Print[\[Epsilon]];
phi1 :=
Simplify[\[Phi]3 /. {ky -> \[Epsilon]*Sin[\[CurlyPhi]]}] // N;
phi2 :=
Simplify[\[Phi]4 /. {ky -> \[Epsilon]*Sin[\[CurlyPhi]]}] // N;
Print[phi1];
Print[phi2];
n = n + 1;
Write[f1, \[Epsilon], OutputForm[kongge], Re[phi1],
OutputForm[kongge], Re[phi2]];
);
];
Close[f1];
Print[OK];
)
出现的问题是：打印的数值有，但导入文件excle就会出现数据错乱，如图。