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個文件信息