三维电磁场FDTD程序(3D-FDTD matlab)PEC边界.docx

三维电磁场FDTD程序（3D-FDTD matlab）PEC边界% FDTD Main Function Jobs to Workers %%***********************************************************************%???? 3-D FDTD code with PEC boundaries%***********************************************************************%%???? This MATLAB M-file implements the finite-difference time-domain%???? solution of Maxwells 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 X cm (x-direction), Y cm (y-direction), and %???? Z cm (z-direction), respectively.%%???? The computational domain is truncated using PEC boundary %???? conditions:%????????? ex(i,j,k)=0 on the j=1, j=jb, k=1, and k=kb planes%????????? ey(i,j,k)=0 on the i=1, i=ib, k=1, and k=kb planes%????????? ez(i,j,k)=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.%%***********************************************************************function [Ex,Ey,Ez]=FDTD3D_Main(handles)globalSimRunStop% if ~isdir(C:\MATLAB7\work\cavity\figures)%???? mkdir C:\MATLAB7\work\cavity\figures% end%********************************