Построить интерференционную картину

Question

Есть код - интерференция света. Здесь в качестве длины волны значения 5. Но хотелось бы взять видимый спектр значения. Как можно реализовать?

from math import sqrt,sin,pi
from numpy import empty
from pylab import imshow,gray,show

wavelength = 5
k = 2*pi/wavelength
xi0 = 1.0
separation = 20.0      # Separation of centers in cm
side = 100.0           # Side of the square in cm
points = 500           # Number of grid points along each side
spacing = side/points  # Spacing of points in cm


# Calculate the positions of the centers of the circles
x1 = side/2 + separation/2
y1 = side/2
x2 = side/2 - separation/2
y2 = side/2

# Make an array to store the heights
xi = empty([points,points],float)

# Calculate the values in the array
for i in range(points):
    y = spacing*i
    for j in range(points):
        x = spacing*j
        r1 = sqrt((x-x1)**2+(y-y1)**2)
        r2 = sqrt((x-x2)**2+(y-y2)**2)
        xi[i,j] = xi0*sin(k*r1) + xi0*sin(k*r2)

# Make the plot
imshow(xi,origin="lower",extent=[0,side,0,side])
gray()
show()

Answer 1

Единицы измерения в коде имеют условный характер (какими назовете единицами, такие единицы и будут - хоть сантиметрами назовите, хоть метрами, программе все равно, главное чтобы везде одна и та же единица измерения была). Можно выбрать, что на самом деле единица измерения - метр, поставить длину волны 550 нанометров (550 * 10**-9), все другие размеры заменить более-менее пропорционально (в том числе сторону рисунка), тогда будет нормальная интерференционная картина:

wavelength = 550E-9 # 550 нм (550 * 10**-9 м)
k = 2*pi/wavelength
xi0 = 1.0
separation = wavelength*4  # между источниками 4 длины волны (0.0022 мм)
side = 1E-5  # Сторона рисунка - 0.01 мм
points = 500
spacing = side/points

Нужно понимать, что это не отображение интерференционной картины на стену (при проекции на стену из-за удаления от источников картина становится "больше"), а как бы "сечение" пространства плоскостью, проходящей через два источника, поэтому все величины (в том числе и сторона рисунка) такие маленькие.

БЛОГ НА HUSL

Построить интерференционную картину

Ответы (1 шт):