usrsetup.m

非连续波导时间域有限差分法仿真（TD-FDM）

% problem parameters
% constants calculation

mu = 4*pi*1e-7;
epsilon = 8.854*1e-12;

speed = 1/(sqrt(mu*epsilon)); % speed of travelling wave

% geometry dimension (unit: m)
L = 0.025;
h = 0.007;
l = 0.05;
H = 0.035;

% origin coordinates values
xo = 0;
yo = 0;

% spatial discretization (unit: cm)
hx = 0.0025;
hy = 0.0035;

% time discretization (unit: sec)
dt = .8e-12;

% Gaussian pulse definition
fc = 0.9e9;     % fc=900MHz

% parameters for update iterations
tol = 1e-10;
MAX_STEP = 1000;
scheme = 0;     % update scheme, '1' for explicit, '0' for implicit

% number of sampling points from TD to FD
N = realpow(2,16);
startFreq = 300e6;
endFreq = 900e6;