%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%TM波的傳播
%今日完成?tm波的傳播

clear?all;?clc;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%?真空常數(shù)設(shè)定
epsz=1/（4*pi*9*10^9）;??????????%真空介電常數(shù)
mu=4*pi*10^（-7）;???????????????%真空磁導(dǎo)率
Z=sqrt（mu/epsz）;???????????????%真空阻抗
epsilon=1;?????????????????????%真空相對介電常數(shù)
sigma=0;???????????????????????%真空電導(dǎo)率

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%?設(shè)定截止頻率
v=3e8;?????????????????????????%波速
f=0.3e9;???????????????????????%截止頻率
lamda=v/f;?????????????????????%波長
k=2*pi/lamda;??????????????????%波數(shù)

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%?設(shè)定網(wǎng)格
IE=150;????????????????????????%x向網(wǎng)格數(shù)
JE=150;????????????????????????%y向網(wǎng)格數(shù)
ISteps=450;????????????????????%迭代次數(shù)
ddx=lamda/20;??????????????????%網(wǎng)格尺寸
dt=ddx/（2*v）;??????????????????%時間步大小

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%?設(shè)定總場
ia=30;
ib=IE-ia;??????????????????????%x向總場邊界位置
ja=30;?????????????????????????
jb=JE-ja;??????????????????????%y向總場邊界位置

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%?PML層數(shù)
npml=10;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%??設(shè)定脈沖源參數(shù)
spread=8;??????????????????????%脈沖寬度
t0=25;?????????????????????????%脈沖高度
%x0=IE/2;
x0=11;
%y0=JE/2;???????????????????????%脈沖位置
%x0=50;
y0=50;
pulse=0;???????????????????????%記錄脈沖大小


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%?計算中的變量內(nèi)存分配
u=zeros（IStepsIEJE）;?
u1=zeros（ISteps）;
u2=zeros（IEJE）;
dz=zeros（IEJE）;???????????????%z向電場密度，歸一化后的
ez=zeros（IEJE）;???????????????%z向電場強(qiáng)度

em=zeros（IEJE）;??%%%%%%%%%%%%%
IE1=IE/3;
IE2=IE1+1;
IE3=IE/2;
IE4=IE3+1;

iz=zeros（IEJE）;???????????????%迭代中間變量
hx=zeros（IEJE）;???????????????%x向磁場變量
hy=zeros（IEJE）;???????????????%y向磁場變量
ihx=zeros（IEJE）;??????????????%x向hx計算的中間變量
ihy=zeros（IEJE）;??????????????%y向hy計算的中間變量
?????for?i=1:IE1;
????????for?j=1:JE;
????????????em（ij）=1;
????????end
?????end
?????for?i=IE2:IE3;
????????for?j=1:JE;
????????????em（ij）=2;
????????end
?????end
?????for?i=IE4:IE;
????????for?j=1:JE;
????????????em（ij）=3;
????????end
?????end
%%?PML初始化設(shè)置
gi2=ones（1IE）;??????????????????%x向系數(shù)
gi3=ones（1IE）;
fi1=zeros（1IE）;
fi2=ones（1IE）;
fi3=ones（1IE）;

gj2=ones（1JE）;??????????????????%y向系數(shù)
gj3=ones（1JE）;
fj1=zeros（1JE）;
fj2=ones（1JE）;
fj3=ones（1JE）;

%%?PML層中阻抗?jié)u變設(shè)置
for?i=1:npml+1;
????xnum=npml-i+1;
????xxn=xnum/npml;
????xn=0.333*（xxn^3）;
????gi2（i）=1/（1+xn）;
????gi2（IE-i+1）=1/（1+xn）;
????gi3（i）=（1-xn）/（1+xn）;
????gi3（IE-i+1）=（1-xn）/（1+xn）;
????
????xxn=（xnum-0.5）/npml;
????xn=0.25*（xxn^3）;
????fi1（i）=xn;
????fi1（IE-i）=xn;
????fi2（i）=1/（1+xn）;
????fi2（IE-i）=1/（1+xn）;
????fi3（i）=（1-xn）/（1+xn）;
????fi3（IE-i）=（1-xn）/（1+xn）;
end;

for?j=1:npml+1;
????xnum=npml-j+1;???
????xxn=xnum/npml;??????
????xn=0.33*（xxn^3）;???
????gj2（j）=1/（1+xn）;???
????gj2（JE-j+1）=1/（1+xn）;???
????gj3（j）=（1-xn）/（1+xn）;???
????gj3（JE-j+1）=（1-xn）/（1+xn）;???
???????
????xxn=（xnum-0.5）/npml;???
????xn=0.25*（xxn^3）;???
????fj1（j）=xn;???
????fj1（JE-j）=xn;???
????fj2（j）=1/（1+xn）;???
????fj2（JE-j）=1/（1+xn）;???
?

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????5618??2019-04-06?22:05??matlabfdtd\todaywork.m

?????文件????8876895??2018-10-30?08:21??matlabfdtd\外國博士正演.pdf

?????目錄??????????0??2019-04-06?22:34??matlabfdtd

-----------?---------??----------?-----??----

??????????????8882513????????????????????3