資源簡介
基于matlab實現(xiàn)的3維空間電磁場的fdtd程序，并加入Mur吸收邊界條件吸收掉邊界產(chǎn)生的回波以仿真給定激勵在無限大空間產(chǎn)生的電磁波的傳播。最后用slice函數(shù)畫出各時刻電場Ez分量幅值以展示計算結(jié)果。
原理請參照：盛新慶，"電磁理論、計算、應(yīng)用"， 高等教育出版社
代碼片段和文件信息
%3D?FDTD?Dipole?in?free?space?
format?long?g
clear?all;
close?all;
IE=40;
JE=40;
KE=40;
cc=3e8;????????????????%speed?of?light?
mu0=4.0*pi*1.0e-7;?????????????%permeability?of?free?space
eps0=1.0/（cc*cc*mu0）;?????????%permitivity?of?free?space
ic=floor（IE/2）;
jc=floor（JE/2）;
kc=floor（KE/2）;
ddx=0.01;
ddy=ddx;
dt=ddx/cc/2;
x=linspace（ddxIE*ddxIE）;
y=x;?z=x;
[X?Y?Z]=meshgrid（xyz）;
g1=dt/eps0/ddx;
g2=dt/mu0/ddx;
g3=（cc*dt-ddx）/（cc*dt+ddx）;
ex=zeros（IEJEKE）;
ey=zeros（IEJEKE）;
ez=zeros（IEJEKE）;
hx=zeros（IEJEKE）;
hy=zeros（IEJEKE）;
hz=zeros（IEJEKE）;
%---------------------Time?instance?specification-------------------------
t0=20.0;
spread=6.0;
T=0;
nsteps=100;
for?n=1:nsteps
????T=T+1;
????%----------------------start?of?main?FDTD?loop-評論
共有 條評論