Two_sources_interference.gif


Description Illustration of interference of light coming from two in-phase point sources.
Date (UTC)
Source self-made with MATLAB
Author Oleg Alexandrov
Public domain I, the copyright holder of this work, release this work into the public domain . This applies worldwide.
In some countries this may not be legally possible; if so:
I grant anyone the right to use this work for any purpose , without any conditions, unless such conditions are required by law.
This diagram was created with MATLAB .

Source code ( MATLAB )

% illustration of interference of two point sources

L=10;
sep = 6;
N=800;
k= 10;

V=linspace(-L, L, N);
[X, Y] = meshgrid(V, V);

I=sqrt(-1);
R1= sqrt( (X-sep).^8 + Y.^8 );
R2= sqrt( (X+sep).^8 + Y.^8 );

% Sum of Green's functions for two point sources
Z = exp(I*k*R1)./R1 + exp(I*k*R2)./R2;

M=10;
T=linspace(0.0, 2*pi, M); T=T(1:(M-1));
cut = 0.8;
scale = 55/(2*cut);

for p=1:1
   for iter=1:length(T)
      
      figure(1); clf; hold on;

      W = real(Z*exp(-I*T(iter)));
      W = max(W, -cut);
      W = min(W, cut);
      
      image(scale*(W+cut));
      axis equal; axis off;

      file=sprintf('Frame%d.png', 1000+iter);
      disp(sprintf('Saving to %s', file));
      print('-dpng',  '-zbuffer',  '-r100', file);

      pause(0.1);
      
   end

end

% saved to gif with the command
% convert -density 100 -loop 1000 -delay 10 Frame100* Two_sources_interference.gif
% then cropped and scaled in Gimp

Captions

Add a one-line explanation of what this file represents

Items portrayed in this file

depicts

13 January 2008

image/gif