matlab非均匀时域有限差分,三维时域有限差分法MATLAB代码

%***********************************************************************

% 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.【矩形谐振腔 10 * 4.8 * 2.0 cm】

%

% The computational domain is truncated using PEC boundary

% conditions:

% ex(i,j,k)=0 on the j=1, j=YB, k=1, and k=ZB planes

% ey(i,j,k)=0 on the i=1, i=XB, k=1, and k=ZB planes

% ez(i,j,k)=0 on the i=1, i=XB, j=1, and j=YB 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.【为了一