利用FDTD的方法计算场值

利用FDTD计算场值

%***********************************************************************
%     3-D FDTD code with PEC boundaries
%***********************************************************************
%
%     Program author: Susan C. Hagness
%                     Department of Electrical and Computer Engineering
%                     University of Wisconsin-Madison
%                     1415 Engineering Drive
%                     Madison WI 53706-1691
%                     608-265-5739
%                     hagness@engr.wisc.edu
%
%     Date of this version:  February 2000
%
%     This MATLAB M-file implements the finite-difference time-domain
%     solution of Maxwell‘s curl equations over a three-dimensional
%     Cartesian space lattice comprised of uniform cubic grid cells.
%     
%     To illustrate the algorithm an air-filled rectangular cavity 
%     resonator is modeled.  The length width and height of the 
%     cavity are 10.0 cm （x-direction） 4.8 cm （y-direction） and 
%     2.0 cm （z-direction） respectively.
%
%     The computational domain is truncated using PEC boundary 
%     conditions:
%          ex（ijk）=0 on the j=1 j=jb k=1 and k=kb planes
%          ey（ijk）=0 on the i=1 i=ib k=1 and k=kb planes
%          ez（ijk）=0 on the i=1 i=ib j=1 and j=jb planes
%     These PEC boundaries form the outer lossless walls of the cavity.
%
%     The cavity is excited by an additive current source oriented
%     along the z-direction.  The source waveform is a differentiated 
%     Gaussian pulse given by 
%          J（t）=-J0*（t-t0）*exp（-（t-t0）^2/tau^2） 
%     where tau=50 ps.  The FWHM spectral bandwidth of this zero-dc-
%     content pulse is approximately 7 GHz. The grid resolution 
%     （dx = 2 mm） was chosen to provide at least 10 samples per 
%     wavelength up through 15 GHz.
%
%     To execute this M-file type “fdtd3D“ at the MATLAB prompt.
%     This M-file displays the FDTD-computed Ez fields at every other
%     time step and records those frames in a movie matrix M which 
%     is played at the end of the simulation using the “movie“ command.
%
%***********************************************************************

clear

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

cc=2.99792458e8;            %speed of light in free space
muz=4.0*pi*1.0e-7;          %permeability of free space
epsz=1.0/（cc*cc*muz）;       %permittivity of free space

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

ie=50;       %number of grid cells in x-direction
je=24;       %number of grid cells in y-direction
ke=10;       %number of grid cells in z-direction

ib=ie+1;     
jb=je+1;   
kb=ke+1;   

is=26;       %location of z-directed 

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2013-06-17 11:00  电磁场时域计算方法概论\
     文件        5483  2011-03-23 20:23  电磁场时域计算方法概论\stright.f90
     文件       10280  2011-03-23 20:23  电磁场时域计算方法概论\wire circle.f90
     目录           0  2013-06-17 11:00  电磁场时域计算方法概论\课件\
     目录           0  2013-06-17 11:00  电磁场时域计算方法概论\课件\\FDTD\
     文件     2203211  2010-04-29 08:09  电磁场时域计算方法概论\课件\\FDTD\abc.pdf
     文件      296497  2010-04-25 21:44  电磁场时域计算方法概论\课件\\FDTD\ch2 FDTD吸收边界MUR（讲稿）.pdf
     文件      303032  2010-05-24 09:09  电磁场时域计算方法概论\课件\\FDTD\ch3 FDTD近场远场外推（讲稿1234）.pdf
     文件      105494  2010-05-31 09:07  电磁场时域计算方法概论\课件\\FDTD\ch4 亚网格和共形网格技术（讲稿2）.pdf
     文件      552346  2010-05-27 06:24  电磁场时域计算方法概论\课件\\FDTD\ch4 亚网格和共形网格技术（讲稿）.pdf
     文件      419280  2010-05-31 09:07  电磁场时域计算方法概论\课件\\FDTD\ch5 周期结构和色散介质（讲稿）.pdf
     文件      936093  2010-06-07 10:11  电磁场时域计算方法概论\课件\\TDFEM2010.pdf
     目录           0  2013-06-17 11:00  电磁场时域计算方法概论\课件\\TDIE\
     文件      339345  2010-03-11 18:36  电磁场时域计算方法概论\课件\\TDIE\TDIE ch1讲课用.pdf
     文件      330720  2010-03-11 18:34  电磁场时域计算方法概论\课件\\TDIE\TDIE ch2（a直导线）讲课用.pdf
     目录           0  2013-06-17 11:00  电磁场时域计算方法概论\作业\
     目录           0  2013-06-17 11:00  电磁场时域计算方法概论\作业\2\
     文件       14943  2011-06-20 17:55  电磁场时域计算方法概论\作业\2\\fdtd-2D.m
     文件      406016  2011-06-20 17:48  电磁场时域计算方法概论\作业\2\\李科（1005121894）.doc
     目录           0  2013-06-17 11:00  电磁场时域计算方法概论\作业\1\
     文件        7828  2011-05-18 09:05  电磁场时域计算方法概论\作业\1\\fdtd.m
     文件      888832  2011-06-20 20:02  电磁场时域计算方法概论\作业\1\\柴水荣1005110280.doc