Краткое описание

Mathematica 11.0 code

\[Lambda] = 1;
\[Sigma] = 0.3 \[Lambda];
k0 = (2 \[Pi])/\[Lambda];
zR = \[Sigma]^2 k0;
w0 = Sqrt[(\[Lambda] zR)/\[Pi]];
w[z_] := w0 Sqrt[1 + z^2/zR^2];
R[z_] := (zR^2 + z^2)/z;
E2[x_, y_, z_] := w0/w[z] E^(I ArcTan[zR, z]) E^-((x^2 + y^2)/(w[z])^2) E^(I k0 (z + (x^2 + y^2)/(2 R[z])));
f = E2[x, 0, z];
psi = f;
p1 = Table[
   Show[
    RegionPlot[True, {z, -5, 5}, {x, 0, 2}, PlotPoints -> 200, ColorFunctionScaling -> False, ColorFunction -> Function[{z, x}, Hue[Arg[psi*E^(-I \[Phi])]/(2 Pi), 1, Abs[psi]^1]], Frame -> False]
    ,
    RegionPlot[True, {z, -5, 5}, {x, -0.2, 0}, PlotPoints -> 200, ColorFunctionScaling -> False, ColorFunction -> Function[{z, x}, Hue[Arg[E^(I (k0 z - \[Pi]/2))*E^(-I \[Phi])]/(2 Pi), 1, 1]], Frame -> False]
    , PlotRange -> {{-5, 5}, {-0.2, 2}}
    ]
   , {\[Phi], 0, 2 \[Pi], \[Pi]/10}];
ListAnimate[p1]

Лицензирование

Я, владелец авторских прав на это произведение, добровольно публикую его на условиях следующей лицензии:

Этот файл доступен на условиях Creative Commons CC0 1.0 Универсальной передачи в общественное достояние (Universal Public Domain Dedication).

Лица, связанные с работой над этим произведением, решили передать данное произведение в общественное достояние, отказавшись от всех прав на произведение по всему миру в рамках закона об авторских правах (а также связанных и смежных прав), в той степени, которую допускает закон. Вы можете копировать, изменять, распространять, исполнять данное произведение в любых целях, в том числе в коммерческих, без получения на это разрешения автора. http://creativecommons.org/publicdomain/zero/1.0/deed.enCC0Creative Commons Zero, Public Domain Dedicationfalsefalse