二维FDTD圆柱散射问题RCS计算matlab代码包含场迭代动画

用matlab实现的二维FDTD代码：包含平面波引入，场值迭代，2阶MUR边界条件，UPML边界条件，近远场外推，后处理RCS计算代码以及圆柱散射的解析解代码。 总之二维散射问题能用到的这个包里都有啦。应该算不叫全了。

%***********************************************************************
%     1-D FDTD code with absorbing boundary conditions
%***********************************************************************
%                     
%                     
%     J.Y.Zhu2013/08/20
%
%     This MATLAB M-file implements the finite-difference time-domain
%     solution of Maxwell‘s curl equations over a one-dimensional
%     Cartesian space lattice comprised of uniform square grid cells.
%
%     To illustrate the algorithm a Guassian pulse propagating 
%     in free space is modeled. 
%
%     The grid resolution （dx = 3 mm） was chosen to provide 20 samples
%     per wavelength at the max frequency we interested in （which in turn
%     provides approximately 10 samples per wavelength at the high end
%     of the hxcitation spectrum around 10 Gex）.
%
%     To execute this M-file type “fdtd_1d“ at the MATLAB prompt.
%     This M-file displays the FDTD-computed Hy and ex fields at 
%     every 5th time step and records those frames in a movie matrix 
%     M which is played at the end of the simulation using the “movie“ 
%     command.
%
%***********************************************************************
close all;
clear all;

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

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

freq = 5e9;                %max frequency of plane wave hxcitation
lambda=cc/freq;             %minimum wavelength of plane wave hxcitation
omega=2.0*pi*freq;          

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

Nx=140;           %number of grid cells in x-direction of fdtd region
Nx_SF_min = 20;   %the boundary of the total-field and scatter field
Nx_SF_max = Nx - 20;
Nxp1=Nx+1; 

dx=3e-3;        %space increment of square lattice
dt=dx/（2.0*cc）;   %time step
x = linspace（0141142）;
nmax=Nx*3;         %total number of time steps

%***********************************************************************
%     Material parameters
%***********************************************************************

media=2;

eps=[1.0 1.0];
sig=[0.0 0.0];
mur=[1.0 1.0];
sim=[0.0 0.0];

%***********************************************************************
%     Wave excitation
%***********************************************************************

tau=2/freq;
delay=0.8*tau;

%***********************************************************************
%     Field arrays
%***********************************************************************
hy=ze

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        6991  2013-09-05 01:09  fdtd_1d.m
     文件       43770  2013-09-06 20:33  forallcode_TM_2D_Cylinder_UPML_Mur2.m
     文件        3830  2013-09-05 06:50  phisbistatic.m
     文件        9080  2013-09-05 15:51  version1_filed_computation.m
     文件       13456  2013-09-05 01:08  version2_planewave.m
     文件         702  2013-08-24 19:18  GetCylinderRCS.m