Subscript[B, L] = 3.0
Subscript[k, L] = Sqrt[\[Epsilon]^2 - Subscript[B, L]^2]
Subscript[b, L] = 0.5
Subscript[k, w] = \[Epsilon]
Subscript[b, R] = 0.5
Subscript[k, R][\[Epsilon]_] = Sqrt[\[Epsilon]^2 - Subscript[B, R]^2]
w = 0.5
Subscript[T,
L][\[Epsilon]_] = (1 + (Sin[
Subscript[k, L]*Subscript[b, L]])^2*(\[Epsilon]*
Subscript[B, L]/(Subscript[k, w]*Subscript[k, L]))^2)^-1
Subscript[T,
R][\[Epsilon]_] = (1 + (Sin[
Subscript[k, R][\[Epsilon]]*Subscript[b, R]])^2*(\[Epsilon]*
Subscript[B,
R]/(Subscript[k, w]*Subscript[k, R][\[Epsilon]]))^2)^-1
Subscript[T, M][\[Epsilon]_] =
Subscript[T, L][\[Epsilon]]*
Subscript[T,
R][\[Epsilon]]/(1 -
Sqrt[1 - Subscript[T, L][\[Epsilon]]]*
Sqrt[1 - Subscript[T, R][\[Epsilon]]])^2
Subscript[D,
S][\[Epsilon]_] = (2*(1 - Subscript[T, L][\[Epsilon]])^(1/4)*(1 -
Subscript[T, R][\[Epsilon]])^(1/4))/(1 -
Sqrt[1 - Subscript[T, L][\[Epsilon]]]*
Sqrt[1 - Subscript[T, R][\[Epsilon]]])
Subscript[\[Phi], L][\[Epsilon]_] =
ArcTan[(Subscript[k, w]^2 + Subscript[k, L]^2 + Subscript[B, L]^2)*
Tan[Subscript[k, L]*Subscript[b, L]]/(2*Subscript[k, L]*Subscript[
k, w])]
Subscript[\[Phi], R][\[Epsilon]_] =
ArcTan[(Subscript[k, w]^2 + Subscript[k, R][\[Epsilon]]^2 +
Subscript[B, R]^2)*
Tan[Subscript[k, R][\[Epsilon]]*Subscript[b, R]]/(2*
Subscript[k, R][\[Epsilon]]*Subscript[k, w])]
\[Phi][\[Epsilon]_] =
Subscript[\[Phi], L][\[Epsilon]] +
Subscript[\[Phi], R][\[Epsilon]] + 2*Subscript[k, w]*w
T[\[Epsilon]_] =
Subscript[T,
M][\[Epsilon]]/(1 + (Subscript[D, S][\[Epsilon]]*
Cos[\[Phi][\[Epsilon]]/2])^2)
p = Table[
Plot[T[\[Epsilon]], {\[Epsilon], 0.5, 3.0}, PlotRange -> {0.0, 1.0},
AspectRatio -> 2/1.25], {Subscript[B, R], 0.6, 3, 0.6}]
Show[p]
2018年03月02日 05点03分
5