FDTD学习资料

关于FDTD算法的讲解和实例，可供大家进行学习使用。

obj=ones（515151）;                           %存放空间媒质编号，包括散射体，迭代次数
Ntime=100;
freq=1e9;                                        %频率上限
Nmedia=2;media=zeros（Nmedia4）;                  %可存Nmedia种介质的电磁参数表
media=[1100;113.72e70];                    %[epslon_rmiu_rsigmasigma_m]

Tin=1/freq;                                      %入射波时间常数
Vin=1;t0=1e-9;
Phi=0;                                           %平面波入射角
WL=40;                                           %数值波长，wavel有40个网格
Ninwave=100;
%---------不变参数-------------
miu0=4e-7*pi;epslon0=8.854187818e-12;            %真空电磁参数
vc=299792458;wavel=vc/freq;                      %光速和波长设置
Z0=sqrt（miu0/epslon0）;                           %波阻抗
%------------------------------
%------导入模型到odj离散空间----------
%readmodel;
%[xsyszs]=size（mmd）;
%obj（101-floor（xs/2）:101-floor（xs/2）+xs-1101-floor（ys/2）:101-floor（ys/2）+ys-1101-floor（zs/2）:101-floor（zs/2）+zs-1）=mmd+1;
obj（24:2824:2824:28）=obj（24:2824:2824:28）+1;
%----------导入模型结束---------------

Izmin=11;Izmax=41;%总场区边界
Ismin=6;Ismax=46;%散射场区边界
Fzmin=9;Fzmax=43;%远场外推边界
TEM_flag=1;        %TE（＝1） or TM（＝2）
Ez=zeros（515151）;Ex=Ez;Ey=Ez;Hx=Ez;Hy=Ez;Hz=Ez;%存放电磁场的空间
Ein=zeros（1Ninwave）;Hin=Ein;%存放1-D FDTD入射电磁场的空间
CA=zeros（NmediaNmediaNmediaNmedia）;CB=CA;%电场迭代系数
CP=zeros（NmediaNmedia）;CQ=CP;%磁场迭代系数
delta=wavel/WL;
dt=delta/（2*vc）;
%---申请吸收边界暂存空间（参考葛德彪P63）----
%6个面，每个面2个切向吸收E分量，每个分量需4层暂存层空间（1:Pn 2:Pn-1 3:Qn 4:Qn-1 ）
Ab1x=zeros（51451）;Ab1z=Ab1x;      %1面xz分量4层
Ab2x=zeros（51451）;Ab2z=Ab2x;      %2面xz分量4层
Ab3x=zeros（51514）;Ab3y=Ab3x;      %3面xy分量4层
Ab4x=zeros（51514）;Ab4y=Ab4x;      %4面xy分量4层
Ab5z=zeros（45151）;Ab5y=Ab5z;      %5面zy分量4层
Ab6z=zeros（45151）;Ab6y=Ab6z;      %6面zy分量4层
%---吸收边界暂存空间 end----

%计算迭代系数
for i=1:Nmedia
    for j=1:Nmedia
        for k=1:Nmedia
            for l=1:Nmedia
                epslon_eff=0.25*（media（i1）+media（j1）+media（k1）+media（l1））;
                sigma_eff=0.25*（media（i3）+media（j3）+media（k3）+media（l3））;
                CA（ijkl）=（1-sigma_eff*dt/（2*epslon0*epslon_eff））/（1+sigma_eff*dt/（2*epslon0*epslon_eff））;
                CB（ijkl）=（dt/（epslon0*delta））/（epslon_eff+sigma_eff*dt/（2*epslon0））;
            end
        end
        miu_eff=0.5*（media（i2）+media（j2））;
        sigmam_eff=0.5*（media（i4）+media（j4））;
        CP（ij）=（1-sigmam_eff*dt/（2*miu0*miu_eff））/（1+sigmam_eff*dt/（2*miu0*miu_eff））;
        CQ（ij）=（dt/miu0/delta）/（miu_eff+sigmam_eff*dt/（2*miu0））;
    end
end
FE=dt/epslon0/delta;FH=dt/miu0/delta; %入射波1-D FDTD迭代系数
EBin=zeros（14）;%1-D FDTD吸收边界变量
%====================================prepare work end=====================================================
%
%=========================================Main loop=======================================================
for Nt=1:Ntime
    Nt
    A（Nt）=Ez（26926）
    %---------creat add_in wave using 1-D FDTD----------------
    for i=1:Ninwave-1
        Hin（i）=Hin（i）-FH*（Ein（i+1）-Ein（i））;
    end
    for i=2:Ninwave
        Ein（i）=Ein（i）-FE*（Hi

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

     文件      23261  2007-11-22 21:01  fdtd资料集合2\40972112fdtd\fdtd\fdtd.asv

     文件      24171  2007-11-22 21:29  fdtd资料集合2\40972112fdtd\fdtd\fdtd.m

     文件     139709  2007-11-17 13:06  fdtd资料集合2\40972112fdtd\fdtd\feiji.txt

     文件        702  2007-11-21 18:39  fdtd资料集合2\40972112fdtd\fdtd\readmodel.m

     文件      24171  2007-11-22 21:29  fdtd资料集合2\fdtd\fdtd.m

     文件       5475  2010-10-25 10:34  fdtd资料集合2\fdtd\fdtd1D.m

     文件      20220  2010-10-25 10:34  fdtd资料集合2\fdtd\fdtd2D.m

     文件       7860  2010-11-03 11:31  fdtd资料集合2\fdtd\fdtd3D.m

     文件     139709  2007-11-17 13:06  fdtd资料集合2\fdtd\feiji.txt

     文件      15628  2010-11-23 19:16  fdtd资料集合2\fdtd\mile.m

     文件       1400  2010-11-18 23:42  fdtd资料集合2\fdtd\per.fig

     文件       2686  2010-11-18 23:15  fdtd资料集合2\fdtd\per.m

     文件       3401  2002-01-01 11:46  fdtd资料集合2\fdtd\zhende2wei.m

     文件       5082  2006-12-14 11:06  fdtd资料集合2\fdtd-matlab\2474_fdtd_1D_demo_v2[1].m

     文件       1589  2006-12-14 11:08  fdtd资料集合2\fdtd-matlab\demo_3dFDTD（1）.m

     文件        807  2006-12-14 11:07  fdtd资料集合2\fdtd-matlab\drv_fdtd[1].m

     文件       4014  2006-12-14 11:15  fdtd资料集合2\fdtd-matlab\FDTD1（1）.m

     文件      19909  2006-12-11 11:46  fdtd资料集合2\fdtd-matlab\fdtd2D.m

     文件       4094  2006-12-14 11:15  fdtd资料集合2\fdtd-matlab\fdtd2de（1）.m

     文件       8540  2006-12-14 11:07  fdtd资料集合2\fdtd-matlab\FDTD3D_Main[1].m

     文件       1437  2006-12-14 11:07  fdtd资料集合2\fdtd-matlab\fdtd_1d_3[1].m

     文件       9222  2006-12-14 11:06  fdtd资料集合2\fdtd-matlab\FDTD_1D_abc[1].m

     文件       5314  2006-12-14 11:15  fdtd资料集合2\fdtd-matlab\FDTD_1D_Erand（1）.m

     文件       4262  2006-12-14 10:33  fdtd资料集合2\fdtd-matlab\FDTD_1D_ill.m

     文件       4262  2006-12-14 11:06  fdtd资料集合2\fdtd-matlab\FDTD_1D_ill[1].m

     文件       2591  2006-12-14 11:06  fdtd资料集合2\fdtd-matlab\fdtd_1D_TEM[1].m

     文件       8964  2006-12-14 10:30  fdtd资料集合2\fdtd-matlab\fdtd_2d_demo_v1_1.m

     文件       2100  2006-12-14 11:08  fdtd资料集合2\fdtd-matlab\fdtd_2d_TM_3（1）.m

     文件       2100  2006-12-13 10:02  fdtd资料集合2\fdtd-matlab\fdtd_2d_TM_3.m

     文件       2286  2006-12-14 11:13  fdtd资料集合2\fdtd-matlab\Lossy_FDTD（1）.m

............此处省略15个文件信息