level 2
o可不可以
楼主
Subscript[\[CurlyPhi], 1] = 15*\[Pi]/180
Subscript[v, 0] = 20
Subscript[s, 1] = UnitStep[\[Epsilon]]
Subscript[s, 2] = UnitStep[\[Epsilon] - Subscript[v, 0] + v/2]
Subscript[s, 3] = UnitStep[\[Epsilon] + v]
Subscript[k, 1][v_] = \[Epsilon]*Cos[Subscript[\[CurlyPhi], 1]]
Subscript[k, y][v_] = \[Epsilon]*Sin[Subscript[\[CurlyPhi], 1]]
Subscript[k, 2][v_] =
Sqrt[(\[Epsilon] - Subscript[v, 0] + v/2)^2 - Subscript[k, y][v]^2]
Subscript[k, 3][v_] = Sqrt[(\[Epsilon] + v)^2 - Subscript[k, y][v]^2]
Subscript[\[CurlyPhi], 2][v_] =
ArcTan[Subscript[k, y][v]/Subscript[k, 2][v]]
Subscript[\[CurlyPhi], 3][v_] =
ArcTan[Subscript[k, y][v]/Subscript[k, 3][v]]
b = 1
Subscript[t, 1][
v_] = (Subscript[s, 1]*Subscript[s, 2]*2*
Cos[Subscript[\[CurlyPhi], 1]]*Cos[Subscript[\[CurlyPhi], 2][v]]*
Exp[-i*Subscript[k, 3][v]*b])/(Subscript[s, 1]*Subscript[s, 2]*
Exp[-i*Subscript[\[CurlyPhi], 1]]*
Cos[Subscript[k, 2][v]*b + Subscript[\[CurlyPhi], 2][v]] +
Subscript[s, 2]*Subscript[s, 3]*
Exp[i*Subscript[\[CurlyPhi], 3][v]]*
Cos[Subscript[k, 2][v]*b - Subscript[\[CurlyPhi], 2][v]] -
i*Sin[Subscript[k, 2][v]*b] -
Subscript[s, 1]*Subscript[s, 3]*
Exp[i*(Subscript[\[CurlyPhi], 3][v] - Subscript[\[CurlyPhi],
1])]*i*Sin[Subscript[k, 2][v]*b])
Subscript[t, 2][
v_] = (Subscript[s, 1]*Subscript[s, 2]*2*
Cos[Subscript[\[CurlyPhi], 1]]*Cos[Subscript[\[CurlyPhi], 2][v]]*
Exp[i*Subscript[k, 3][v]*b])/(Subscript[s, 1]*Subscript[s, 2]*
Exp[i*Subscript[\[CurlyPhi], 1]]*
Cos[Subscript[k, 2][v]*b + Subscript[\[CurlyPhi], 2][v]] +
Subscript[s, 2]*Subscript[s, 3]*
Exp[-i*Subscript[\[CurlyPhi], 3][v]]*
Cos[Subscript[k, 2][v]*b - Subscript[\[CurlyPhi], 2][v]] -
i*Sin[Subscript[k, 2][v]*b] -
Subscript[s, 1]*Subscript[s, 3]*
Exp[-i*(Subscript[\[CurlyPhi], 3][v] - Subscript[\[CurlyPhi],
1])]*i*Sin[Subscript[k, 2][v]*b])
T[v_] = Subscript[t, 1][v]*Subscript[t, 2][v]
p = Table[
Plot[T[v], {v, 0, 20}, PlotRange -> {0.0, 1.0},
AspectRatio -> 1/0.8], {\[Epsilon], 5, 30, 5}]
Show[p]
2018年03月02日 05点03分
1
Subscript[v, 0] = 20
Subscript[s, 1] = UnitStep[\[Epsilon]]
Subscript[s, 2] = UnitStep[\[Epsilon] - Subscript[v, 0] + v/2]
Subscript[s, 3] = UnitStep[\[Epsilon] + v]
Subscript[k, 1][v_] = \[Epsilon]*Cos[Subscript[\[CurlyPhi], 1]]
Subscript[k, y][v_] = \[Epsilon]*Sin[Subscript[\[CurlyPhi], 1]]
Subscript[k, 2][v_] =
Sqrt[(\[Epsilon] - Subscript[v, 0] + v/2)^2 - Subscript[k, y][v]^2]
Subscript[k, 3][v_] = Sqrt[(\[Epsilon] + v)^2 - Subscript[k, y][v]^2]
Subscript[\[CurlyPhi], 2][v_] =
ArcTan[Subscript[k, y][v]/Subscript[k, 2][v]]
Subscript[\[CurlyPhi], 3][v_] =
ArcTan[Subscript[k, y][v]/Subscript[k, 3][v]]
b = 1
Subscript[t, 1][
v_] = (Subscript[s, 1]*Subscript[s, 2]*2*
Cos[Subscript[\[CurlyPhi], 1]]*Cos[Subscript[\[CurlyPhi], 2][v]]*
Exp[-i*Subscript[k, 3][v]*b])/(Subscript[s, 1]*Subscript[s, 2]*
Exp[-i*Subscript[\[CurlyPhi], 1]]*
Cos[Subscript[k, 2][v]*b + Subscript[\[CurlyPhi], 2][v]] +
Subscript[s, 2]*Subscript[s, 3]*
Exp[i*Subscript[\[CurlyPhi], 3][v]]*
Cos[Subscript[k, 2][v]*b - Subscript[\[CurlyPhi], 2][v]] -
i*Sin[Subscript[k, 2][v]*b] -
Subscript[s, 1]*Subscript[s, 3]*
Exp[i*(Subscript[\[CurlyPhi], 3][v] - Subscript[\[CurlyPhi],
1])]*i*Sin[Subscript[k, 2][v]*b])
Subscript[t, 2][
v_] = (Subscript[s, 1]*Subscript[s, 2]*2*
Cos[Subscript[\[CurlyPhi], 1]]*Cos[Subscript[\[CurlyPhi], 2][v]]*
Exp[i*Subscript[k, 3][v]*b])/(Subscript[s, 1]*Subscript[s, 2]*
Exp[i*Subscript[\[CurlyPhi], 1]]*
Cos[Subscript[k, 2][v]*b + Subscript[\[CurlyPhi], 2][v]] +
Subscript[s, 2]*Subscript[s, 3]*
Exp[-i*Subscript[\[CurlyPhi], 3][v]]*
Cos[Subscript[k, 2][v]*b - Subscript[\[CurlyPhi], 2][v]] -
i*Sin[Subscript[k, 2][v]*b] -
Subscript[s, 1]*Subscript[s, 3]*
Exp[-i*(Subscript[\[CurlyPhi], 3][v] - Subscript[\[CurlyPhi],
1])]*i*Sin[Subscript[k, 2][v]*b])
T[v_] = Subscript[t, 1][v]*Subscript[t, 2][v]
p = Table[
Plot[T[v], {v, 0, 20}, PlotRange -> {0.0, 1.0},
AspectRatio -> 1/0.8], {\[Epsilon], 5, 30, 5}]
Show[p]