level 3
可魯贝洛斯
楼主
w0 = 1; s = 0.1 ;
EAx = (x - s + 1)/2 Exp (-((x^2 + y^2)/w0^2));
EAy = I (-x + s + 1)/2 Exp (-((x^2 + y^2)/w0^2));
EBx = (I y + 1)/2 Exp (-((x^2 + y^2)/w0^2));
EBy = ( y + I)/2 Exp (-((x^2 + y^2)/w0^2));
e1xA = Abs[EAx]; e2xA = e1xA^2;
e1yA = Abs[EAy]; e2yA = e1yA^2;
e1xB = Abs[EBx]; e2xB = e1xB^2;
e1yB = Abs[EBy]; e2yB = e1yB^2;
S11A = e2xA - e2yA;
S11B = e2xB - e2yB;
S1 = S11A + S11B;
R1A = Conjugate[EAx];
R1B = Conjugate[EBx];
S2A = 2 Re[R1A EAy];
S2B = 2 Re[R1B EBy];
S2 = S2A + S2B;
t2 = ContourPlot[N[Arg[S1 + I S2]], {x, -0.6`, 0.6}, {y, -0.6`, 0.6},
ContourStyle -> Automatic, ContourShading -> True, PlotRange -> All,
PlotPoints -> 50, ColorFunction -> "Rainbow",
Contours -> {-\[Pi], -(1/8) (6 \[Pi]), -(1/8) (4 \[Pi]), -(1/
8) (2 \[Pi]), 0, (2 \[Pi])/8, (4 \[Pi])/8, (6 \[Pi])/8, (
8 \[Pi])/8}]
求帮忙看一下代码的问题
2023年03月27日 03点03分
1
EAx = (x - s + 1)/2 Exp (-((x^2 + y^2)/w0^2));
EAy = I (-x + s + 1)/2 Exp (-((x^2 + y^2)/w0^2));
EBx = (I y + 1)/2 Exp (-((x^2 + y^2)/w0^2));
EBy = ( y + I)/2 Exp (-((x^2 + y^2)/w0^2));
e1xA = Abs[EAx]; e2xA = e1xA^2;
e1yA = Abs[EAy]; e2yA = e1yA^2;
e1xB = Abs[EBx]; e2xB = e1xB^2;
e1yB = Abs[EBy]; e2yB = e1yB^2;
S11A = e2xA - e2yA;
S11B = e2xB - e2yB;
S1 = S11A + S11B;
R1A = Conjugate[EAx];
R1B = Conjugate[EBx];
S2A = 2 Re[R1A EAy];
S2B = 2 Re[R1B EBy];
S2 = S2A + S2B;
t2 = ContourPlot[N[Arg[S1 + I S2]], {x, -0.6`, 0.6}, {y, -0.6`, 0.6},
ContourStyle -> Automatic, ContourShading -> True, PlotRange -> All,
PlotPoints -> 50, ColorFunction -> "Rainbow",
Contours -> {-\[Pi], -(1/8) (6 \[Pi]), -(1/8) (4 \[Pi]), -(1/
8) (2 \[Pi]), 0, (2 \[Pi])/8, (4 \[Pi])/8, (6 \[Pi])/8, (
8 \[Pi])/8}]
求帮忙看一下代码的问题
