I wrote a few lines that numerically solve Maxwell's equations.
The result is a moving wave that looks like a single pulse.
This looks strange to me because I expect waves to move in oscillator fashion, perhaps like cos(x).
Why doesn't the wave appear to oscillate?
How does this kind of Gaussian wave work? Can it be broken into solutions?
Do similar solutions exist in 2D? Are they still Gaussian?
%%Super basic FDTD
clear all;
gsz=100;
ez=zeros(gsz);hy=zeros(gsz);
imp=377.0;
%
steps = 100;
for i=1:steps
for j=1:(gsz-1)
hy(j)=hy(j)+(ez(j+1)-ez(j))/imp;
end
for j=2:gsz
ez(j)=ez(j)+(hy(j)-hy(j-1))*imp;
end
ez(gsz)= ez(gsz)+exp(-((i-30)^2) /100);
subplot(2,1,1);
plot(ez);
subplot(2,1,2);
plot(hy);
M(i)=getframe;
end
