波干涉相长时的演示程序

此程序用来演示干涉相长现象，两列波相遇时，中间绿色箭头表示亮条纹。

程序演示

# Written by Ruth Chabay
scene = canvas(width = 300,height = 300,forward = vector(0,-1,-2),background = color.white,range = 0.35)
lamb = 0.1     
c = 3e8
omega = 2*pi*c/lamb
d = 2*lamb ## slit spacing  

Evec1 = []
Evec2 = []

dist_to_screen = 4*lamb    
scene.center = vector(dist_to_screen*.65,-d/2,0)
scene.caption = "点击屏幕即可查看电磁波"

ds = lamb/20
dt = lamb/c/100
E0 = lamb/3

screen = curve(pos = [vector(dist_to_screen,0,0),vector(dist_to_screen,0,2*d)], color=vector(0.6,0.6,0.6), radius=0.001)
slit1 = vector(0, 0, -d/2) 
slit2 = vector(0, 0, d/2)  

max1 = vector(dist_to_screen,0,0)
max2 = vector(dist_to_screen,0,0.24)
min1 = vector(dist_to_screen,0,0.105) 

## vectors from slits to max2:
r1 = max2 - slit1
r2 = max2 - slit2

# vectors from slits to loc. on screen
dr1 = ds*norm(r1)
dr2 = ds*norm(r2)

rr1 = slit1 + vector(0,0,0) ## current loc along wave 1
rr2 = slit2 + vector(0,0,0) ## current loc along wave 2 

i1 = None
i2 = None

## create first wave
for i in range(120):
    ea = arrow(pos=rr1, axis=vector(0,E0*cos(2*pi*mag(rr1-slit1)/lamb),0), color=color.red, shaftwidth=lamb/40)
    ba = arrow(pos=rr1, axis=vector(0,0,0), color=vector(.4,.4,1), shaftwidth=lamb/40, visible=False)
    ea.B = ba
    if abs(ea.pos.x - dist_to_screen) < 0.002 and i1 == None:
        i1 = ea
        i1.visible = 0
    else:
        Evec1.append(ea)
    rr1 = rr1 + dr1
    
## create second wave
ct = 0
while ct < 100:
    ea = arrow(pos=rr2, axis=vector(0,E0*cos(2*pi*mag(rr2-slit2)/lamb),0), color=vector(1.,.6,0), shaftwidth=lamb/40)
    ba = arrow(pos=rr2, axis=vector(0,0,0), color=color.cyan, shaftwidth=lamb/40, visible=False)
    ea.B = ba
    if abs(ea.pos.x - dist_to_screen) < 0.002 and i2 == None:
        i2 = ea
    else:
        Evec2.append(ea)
    rr2 = rr2 + dr2
    ct = ct + 1

i2.color = color.green
i2.B.color = color.magenta

def clicked(evt):
    for a in Evec1:
        a.B.visible = not a.B.visible
    for a in Evec2:
        a.B.visible = not a.B.visible
    i2.B.visible = not i2.B.visible

scene.bind('click', clicked)

t=0
while True:
    rate(50)
    t = t+dt
    for ea in Evec1:
       ea.axis = vector(0,E0*cos(omega*t - 2*pi*mag(ea.pos-slit1)/lamb),0)
       ea.B.axis = -cross(ea.axis,norm(r1))*.7       
    for ea in Evec2:
       ea.axis = vector(0,E0*cos(omega*t - 2*pi*mag(ea.pos-slit2)/lamb),0)
       ea.B.axis = -cross(ea.axis,norm(r2))*.7
    sum = E0*cos(omega*t - 2*pi*mag(i1.pos-slit1)/lamb)+E0*cos(omega*t - 2*pi*mag(i2.pos-slit2)/lamb) # superposition
    i2.axis = vector(0,sum,0) 
    i2.B.axis = -cross(i2.axis,norm(r2))*.7