模拟光束在二维三维介质中传输，用FDTD法，Matlab语言编写

模拟光束在二维三维介质中传输，用FDTD法，Matlab语言编写

function FDTD_2D_Yee_2D_TM_final

%constants
c_0 = 3E8;                   % Speed of light in free space
mu_0 = 4.0*pi*1.0e-7;        % Permeability of free space
eps_0 = 1.0/（c_0*c_0*mu_0）;  % Permittivty of free space

%to be set
Nx = 100;                    % Number of cells in x-direction
Ny = 100;                    % Number of cells in y-direction
Nt = 500;                    % Number of time steps

c_ref = c_0;                 % Reference velocity
eps_ref = eps_0;
mu_ref = mu_0;

f_0 = 10e9;                  % Excitation frequency
f_ref = f_0;                 % Reference frequency 

omega_0 = 2.0*pi*f_0;        % Excitation circular frequency
omega_ref = omega_0;

lambda_ref = c_ref/f_ref;    % Reference wavelength

Dx_ref = lambda_ref/20;      % Reference cells width
Dy_ref = lambda_ref/20;

X = Nx*Dx_ref;
Y = Ny*Dy_ref;

r = 0.5;                     % Normalized factor
Dt_ref = r/c_ref*Dx_ref;     % Reference time step
Dt = Dt_ref;

% initiali

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       8353  2006-11-29 19:32  1121204768158150.m

----------- ---------  ---------- -----  ----

                 8353                    1