資源簡(jiǎn)介
2維FDTD全波數(shù)值算法+摩爾吸收邊界,基于matlab模擬二維平面波的空間傳遞,并且設(shè)置金屬塊模擬波的反射過程
代碼片段和文件信息
clear?all;
u=4.*pi.*10.^(-7);
e=8.854187817.*10.^(-12);
cig=0;
c=3*10^8;
Imax=200;
Jmax=200;
Kmax=200;
Nmax=800;
d=0.001;
dt=d/2/c;
s=(c*dt-d)/(c*dt+d);???????????????????????????????%邊界計(jì)算公式系數(shù)(相鄰)
s1=(c*dt-sqrt(2)*d)/(c*dt+sqrt(2)*d);??????????????%邊界計(jì)算公式系數(shù)(斜對(duì)角)
Ca=1;
Cb=dt/e/d;
Da=1;
Db=dt/u/d;
Ex=zeros(ImaxJmax);
Ey=zeros(ImaxJmax);
Ez=zeros(ImaxJmax);
Hx=zeros(ImaxJmax);
Hy=zeros(ImaxJmax);
Hz=zeros(ImaxJmax);
A=zeros(ImaxJmax);?????????????????????????????????%三組記錄上一時(shí)刻的數(shù)組
B=zeros(ImaxJmax);
C=zeros(ImaxJmax);
for?n=1:Nmax
????
????
?????A=Ex;??????????????????????????????????????????%記錄上一時(shí)刻Ex????
????for?x=1:60??????????????????????????????????????%計(jì)算Ex(1-60)?
????????for?y=2:Jmax
?????????Ex(xy)=Ca*Ex(xy)+Cb*(Hz(xy)-Hz(xy-1));
?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????文件???????3481??2017-04-26?14:57??FDTD2withM.m
-----------?---------??----------?-----??----
?????????????????3481????????????????????1
評(píng)論
共有 條評(píng)論