【电磁波】一维FDTD方法电磁波传播及反射透射仿真【含Matlab源码 3871期】

⛄一、获取代码方式

获取代码方式1：
完整代码已上传我的资源：【电磁波】基于matlab一维FDTD方法电磁波传播及反射透射仿真【含Matlab源码 3871期】
点击上面蓝色字体，直接付费下载，即可。

获取代码方式2：
付费专栏Matlab物理应用（初级版）

备注：
点击上面蓝色字体付费专栏Matlab物理应用（初级版），扫描上面二维码，付费29.9元订阅海神之光博客付费专栏Matlab物理应用（初级版），凭支付凭证，私信博主，可免费获得1份本博客上传CSDN资源代码（有效期为订阅日起，三天内有效）；
点击CSDN资源下载链接：1份本博客上传CSDN资源代码

⛄二、部分源代码

clear;close all;

cc=2.99792458e8;
muz=4.0pi1.0e-7;
epsz=1.0/(ccccmuz);
freq=1e+9;
Tc=0.2e-9;
Tau=1/freq;
lambda=cc/freq;
omega=2.0pifreq;
display=2;
Source_type=5;
M_factor=1;
nmax=380M_factor;
nstep=50;
ie=100M_factor;
je=100M_factor;
ib=ie+1;
jb=je+1;
icenter=65M_factor; jcenter=55M_factor;
Dim_Tumor=10M_factor;
imp_x=50M_factor;
imp_y=50M_factor;
d_inter=20M_factor;
det_x=20M_factor;
det_y1=50M_factor;
det_y2=30M_factor;
det_y3=70M_factor;
for Tumor_case=1:2
switch Tumor_case
case 1
width=5;
case 2
width=-1;
end
dx=1e-3/M_factor;
dt=dx/(2.0cc);
iebc=8M_factor;
jebc=8M_factor;
rmax=0.000001;
orderbc=2;
ibbc=iebc+1;
jbbc=jebc+1;
iefbc=ie+2iebc;
jefbc=je+2jebc;
ibfbc=iefbc+1;
jbfbc=jefbc+1;

% Material parameters
media=2;
eps=[1.0 1.0];
sig=[0.0 3.73*1e7];
mur=[1.0 1.0];
sim=[0.0 0.0];

% Wave excitation
source=zeros(1,nmax);
switch Source_type
case 5
for n=1:nmax
t(n)=dtn;
source(n)=-sqrt(exp(1))(2pi/Tau)(t(n)-Tc)exp(-(1/2)(2pi(t(n)-Tc)/Tau)^2); % normalized
end
case 8
n1=int2str((Tc-Tau/2)/dt);
n2=int2str((Tc-Tau/6)/dt);
n3=int2str((Tc+Tau/6)/dt);
n4=int2str((Tc+Tau/2)/dt);
for n=60:100
t(n)=dtn;
source(n)=1;
end
for n=140:180
t(n)=dtn;
source(n)=-1;
end
for n=1:nmax
t(n)=dt*n;
end
end

⛄三、运行结果

⛄四、matlab版本及参考文献

1 matlab版本
2014a

2 参考文献
[1]张鹏，陈国平，王宇. 基于特征线法和特征参数法的管道瞬变流计算[J]. 河海大学学报（自然科学版）, 2014, 42(6): 602-607.

3 备注
简介此部分摘自互联网，仅供参考，若侵权，联系删除