fdtd22.m

FDTD analysis in two dimension using MATLAB environment....

%1D simulation (FDTD(2,2))
clear;
NX=100;
m0=12.566E-7; %permability of free space
e0=8.854E-12; %permittivity of free space
er=1; %relative permittivity
e=er*e0;
m=m0;
c=1/sqrt(e*m); %velocity of light

Q=0.5; %Courant stability Criterion 
dz=0.05; %space step
dt=Q*dz/c; %time step

Nt=200; %Number of time steps

c1=dt/(e*dz);
c2=dt/(m*dz);

EX(1:NX)=0;
HY(1:NX-1)=0;
P(1:NX)=0;
% source constants
t0=8e-9;
T=1e-9;
zsource=1; 
% source
for n=1:Nt
f(n)=exp(-(n*dt-t0).^2./T^2);   
end

for n=1:Nt
   EX(zsource)=f(n); %source
   EX(2:NX-1)=EX(2:NX-1)-c1*(HY(2:NX-1)-HY(1:NX-2));
   %Observation point
   P(n)=EX(2);
      
   HY(1:NX-1)=HY(1:NX-1)-c2*(EX(2:NX)-EX(1:NX-1));
   
end
t=dt:dt:Nt*dt;
plot(t./1e-9,P,'r');