%***********************************************************************
%?????1-D?FDTD?code?with?absorbing?boundary?conditions
%***********************************************************************
%?????????????????????
%?????????????????????
%?????J.Y.Zhu2013/08/20
%
%?????This?MATLAB?M-file?implements?the?finite-difference?time-domain
%?????solution?of?Maxwell‘s?curl?equations?over?a?one-dimensional
%?????Cartesian?space?lattice?comprised?of?uniform?square?grid?cells.
%
%?????To?illustrate?the?algorithm?a?Guassian?pulse?propagating?
%?????in?free?space?is?modeled.?
%
%?????The?grid?resolution?（dx?=?3?mm）?was?chosen?to?provide?20?samples
%?????per?wavelength?at?the?max?frequency?we?interested?in?（which?in?turn
%?????provides?approximately?10?samples?per?wavelength?at?the?high?end
%?????of?the?hxcitation?spectrum?around?10?Gex）.
%
%?????To?execute?this?M-file?type?“fdtd_1d“?at?the?MATLAB?prompt.
%?????This?M-file?displays?the?FDTD-computed?Hy?and?ex?fields?at?
%?????every?5th?time?step?and?records?those?frames?in?a?movie?matrix?
%?????M?which?is?played?at?the?end?of?the?simulation?using?the?“movie“?
%?????command.
%
%***********************************************************************
close?all;
clear?all;

%***********************************************************************
%?????Fundamental?constants
%***********************************************************************

muz=4.0*pi*1.0e-7;???????????????????????%permeability?of?free?space
epsz=1.0/（35950207149.4727056*pi）;???????%permittivity?of?free?space
cc=1/sqrt（muz*epsz）;?????????????????????%speed?of?light?in?free?space
z0?=?sqrt（muz/epsz）;

freq?=?5e9;????????????????%max?frequency?of?plane?wave?hxcitation
lambda=cc/freq;?????????????%minimum?wavelength?of?plane?wave?hxcitation
omega=2.0*pi*freq;??????????

%***********************************************************************
%?????Grid?parameters
%***********************************************************************

Nx=140;???????????%number?of?grid?cells?in?x-direction?of?fdtd?region
Nx_SF_min?=?20;???%the?boundary?of?the?total-field?and?scatter?field
Nx_SF_max?=?Nx?-?20;
Nxp1=Nx+1;?

dx=3e-3;????????%space?increment?of?square?lattice
dt=dx/（2.0*cc）;???%time?step
x?=?linspace（0141142）;
nmax=Nx*3;?????????%total?number?of?time?steps

%***********************************************************************
%?????Material?parameters
%***********************************************************************

media=2;

eps=[1.0?1.0];
sig=[0.0?0.0];
mur=[1.0?1.0];
sim=[0.0?0.0];

%***********************************************************************
%?????Wave?excitation
%***********************************************************************

tau=2/freq;
delay=0.8*tau;

%***********************************************************************
%?????Field?arrays
%***********************************************************************
hy=ze

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????6991??2013-09-05?01:09??fdtd_1d.m
?????文件???????43770??2013-09-06?20:33??forallcode_TM_2D_Cylinder_UPML_Mur2.m
?????文件????????3830??2013-09-05?06:50??phisbistatic.m
?????文件????????9080??2013-09-05?15:51??version1_filed_computation.m
?????文件???????13456??2013-09-05?01:08??version2_planewave.m
?????文件?????????702??2013-08-24?19:18??GetCylinderRCS.m