VERSCHRÄNKTE WELT - DER MIKROKOSMOS DER QUANTEN (3)

 

Die Vorlesungsfolien können Sie hier als pdf herunterladen...

Python-Skript Interferenz am Doppelspalt

import numpy as np

from matplotlib import pyplot as plt

# Eingabebereich
####################################################################
gbereich = 200
spaltabstand = 0.5
amplitude = 4
frequenz = 10
####################################################################

xmin, xmax = 0, 20
ymin, ymax = -10, 10

xpoints, ypoints = gbereich, gbereich

x = np.linspace(xmin, xmax, xpoints)
y = np.linspace(ymin, ymax, ypoints)

xx, yy = np.meshgrid(x, y, sparse = False)

points = np.concatenate([xx.reshape(-1, 1), yy.reshape(-1, 1)], axis=-1)
source1=np.array([0, spaltabstand])
source2=np.array([0, -spaltabstand])
points1 = points - source1
points2 = points - source2

A1= amplitude
A2= amplitude
k = frequenz

wave1 = A1*(np.sin( k * (points1[:, 0]**2 + points1[:, 1]**2)**0.5))
wave2 = A2*(np.sin( k * (points2[:, 0]**2 + points2[:, 1]**2)**0.5))

# Überlagerung der beiden Wellen
A = (wave1 + wave2)

intensity = A

# Plotten
plt.figure(figsize=(8, 8))
plt.xlim(xmin, xmax)
plt.ylim(ymin, ymax)
plt.scatter(points[:, 0], points[:, 1], c = intensity) 
plt.scatter(*source1, c='red')
plt.scatter(*source2, c='red')