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

Mathematica 11.0 code

\[Lambda] = 1;
k0 = (2 \[Pi])/\[Lambda];
G[r_] := I/4 HankelH1[0, k0 Norm[r] ];
s = Table[{h,(Abs[NIntegrate[G[{Cos[\[Theta]], Sin[\[Theta]]}] - G[{Cos[\[Theta]], Sin[\[Theta]]} - {2*h, 0}], {\[Theta], 0, 2 \[Pi]}] ]/s0)^2}, {h, 3, 0, -0.01}];
p1 = Table[
   GraphicsRow[{
DensityPlot[Abs[G[{x, y}] - G[{x + 2*h, y}] ]^2, {x, -3, 7}, {y, -5, 5}, PlotPoints -> 100, Frame -> False, ColorFunction -> "AvocadoColors", PlotLabel -> "|E\!\(\*SuperscriptBox[\(|\), \(2\)]\)", 
      PlotRange -> {0, 0.02}, LabelStyle -> {Black, Bold}, Epilog -> {Gray, Polygon[{{-5, -5}, {-h, -5}, {-h, 5}, {-5, 5}}]}]
     ,
     ListPlot[s, PlotRange -> {0, 2}, Joined -> True, Axes -> False, Frame -> True, LabelStyle -> {Black, Bold}, PlotStyle -> {Black, Thick}, FrameLabel -> {"Distance from the mirror / \[Lambda]", 
        "P/\!\(\*SubscriptBox[\(P\), \(\[Infinity]\)]\)"}, Epilog -> {PointSize[0.03], Red, Point[{h, Evaluate[(1/s0 Abs[NIntegrate[G[{Cos[\[Theta]], Sin[\[Theta]]}] - G[{Cos[\[Theta]], Sin[\[Theta]]} - {2*h, 0}], {\[Theta], 0, 2 \[Pi]}] ])^2]}]}]
     }, ImageSize -> Large]
   , {h, 3, 0, -0.05}];
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