电磁波的加上Mur吸收边界条件的FDTD算法

最新推荐文章于 2024-04-16 21:55:06 发布

weixin_38005545

最新推荐文章于 2024-04-16 21:55:06 发布

阅读量3.4k

点赞数 6

文章标签： Mur吸收边界，FDTD算法，matlab实现

本文链接：https://blog.csdn.net/weixin_38005545/article/details/102977253

版权

本文介绍了电磁波的时域有限差分FDTD算法在模拟开域问题中的应用，特别是如何结合Mur吸收边界条件来减少反射。通过matlab程序展示了TM波开域仿真，并附有仿真后的电磁波传播图形和三维空间显示效果。

摘要由CSDN通过智能技术生成

电磁波的时域有限差分FDTD算法

从物理观点来说：模拟开域问题，只有在微波实验室墙壁上傅一层吸波材料使其波在此界面无反射，形成微波暗室

从计算观点：模拟开域问题，在截断边界上FDTD的E和H分量的计算公式中，至少有一个环绕分量位于截断边界之外，需要特殊处理，既是需要引入吸收边界条件。

仿真结果实现

TM波开域加上Mur吸收边界条件的仿真

matlab程序如下

c=3*10^8;
e0=8.854*10^(-12);
u0=4*pi*10^(-7);
lamda=0.01;%波长
deta_x=0.001;
deta_y=0.001;%空间网格长度
x_max=400;
y_max=400;%网格数
deta_t=1/c/sqrt((1/deta_x)^2+(1/deta_x)^2)%时间步长
step=600;%迭代次数,也就是迭代时间数
source_x=200;
source_y=200;%电流源位置
field_x=180;
field_y=180;%观察点位之可变
frequency=15*10^9;%波长0.02
omiga=2*pi*frequency;
t0=0.0667*10^(-9);%高斯时间中心
tt=t0;%高斯源
result=zeros(step,1);
%常量置零
hx1=zeros(x_max+1,y_max);%存Hx(n-1/2)
hx2=zeros(x_max+1,y_max);;%存Hx(n+1/2)
hy1=zeros(x_max,y_max+1);%存Hy(n-1/2)