求助,公式作图一直不出图,是代码的问题吗?
mathematica吧
全部回复
仅看楼主
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
level 3
\[Lambda]A = 0.84/1000; \[Lambda]B = 1.06 /1000; kA = (
2 \[Pi])/\[Lambda]A; kB = (
2 \[Pi])/\[Lambda]B; w0 = 1; s = 0.1 ; z = 200; f = 200;
uA = 1/w0^2 + I kA (1/(2 z) - 1/(2 f));
v1A = (I kA x)/(2 z); v2A = (I kA y)/(2 z);
gA = (I \[Pi] )/(2 \[Lambda]A uA^2 z); MA =
Exp[v1A^2/uA + v2A^2/uA - (I kA)/(2 z) (x^2 + y^2)];
uB = 1/w0^2 + I kB (1/(2 z) - 1/(2 f));
v1B = (I kB x)/(2 z); v2B = (I kB y)/(2 z);
gB = (I \[Pi] )/(2 \[Lambda]B uB^2 z); MB =
Exp[v1B^2/uB + v2B^2/uB - (I kB )/(2 z) (x^2 + y^2)];
EAx = gA MA[v1A + uA (1 - s)];
EAy = I gA MA[(1 + s) uA - v1A];
EBx = gB MB (uB + I v2B);
EBy = gB MB (I uB + v2B);
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分 2
1楼和2楼全是括号的问题。把Mathematica里中括号和小括号的作用好好复习下,再一点点检查过去。
2023年04月01日 03点04分
@xzcyr 谢谢,目前已经改正了。
2023年04月10日 08点04分
level 6
运行成功,你查一下帮助看看Exp怎么用,应该不是(),是[ ]
2023年03月28日 11点03分 5
非常感谢![真棒]
2023年03月28日 11点03分
大佬求看看二楼的程序哪儿有问题
2023年03月28日 11点03分
1