Plane Wave Addition

128 days ago by qwerty123xx

f(t,z, dw, w1, k1)=2*(cos((dw/(2*w1))*(k1)*z-(dw*t)/2))*cos(k1*z-w1*t) 
       
f1=f(dw=1, w1=50, k1=4000000*pi, t=1) 
       
f1 
       
2*cos(40000*pi*z - 1/2)*cos(4000000*pi*z - 50)
2*cos(40000*pi*z - 1/2)*cos(4000000*pi*z - 50)
plot(f1, 0, .00004, gridlines=True) 
       
plot(f1^2, 0, .00004, gridlines=True) 
       
plot(f(dw=1, w1=50, k1=4000000*pi, t=2), 0, .00004, gridlines=True) 
       
plot(f(dw=1, w1=50, k1=4000000*pi, t=2)^2, 0, .00004, gridlines=True) 
       
plot(f(dw=1, w1=50, k1=4000000*pi, t=3), 0, .00004, gridlines=True) 
       
plot(f(dw=1, w1=50, k1=4000000*pi, t=3)^2, 0, .00004, gridlines=True) 
       
a = animate([f(dw=1, w1=50, k1=4000000*pi, t=k) for k in srange(0,2*pi,0.3)], xmin=0, xmax=.00004, figsize=[10,10]) 
       
       
Animation with 21 frames
Animation with 21 frames
a.show() 
       
b = animate([f(dw=1, w1=50, k1=4000000*pi, t=k)^2 for k in srange(0,2*pi,0.3)], xmin=0, xmax=.00004, figsize=[10,10]) 
       
b.show()