二维FDTD TE波圆柱仿真(国外英语资料).doc

二维FDTD TE波圆柱仿真(国外英语资料)%二维FDTD TE波圆柱仿真清除所有；关闭所有；中图分类号；%定义常数% -------------------------C = 3.0e8；1.2566e-6μ=每股收益= 8.8542e-12；F = 1e9；%频率λ= C / F；%波w_max = 2×π×2×F；cylinder_circ =λ* 8；%圆柱的周长，改变圆柱的尺cylinder_rad = cylinder_circ／2／PI；%圆柱的半径定义FDTD网格%% ----------------------del_s =λ/ 20；%每最小波长20个采样del_t = 0.5 * del_s / C；%迭代时间步长dim_s = 5；s_range = dim_s×λ；s_range =细胞（s_range / del_s）* del_s；%计算区域的长度s_cells = s_range / del_s；细胞= s_cells；%划分的网格数节点=细胞+ 1；%采样点数%定义时间脉冲源% ----------------------------------------dev_larger = 2；DEV = 1 / w_max * dev_larger；死亡= 4；意味着死亡；T = linspace（0.9×DEV，1000）；术语=（t-mean）；脉冲=（1 / sqrt（2×π）/ dev ^ 3期）；脉冲=脉冲* *（（1 / 2 /开发单位2）*长期；pulsenorm = max（脉冲）；P = s_range / 2；dead_s = 1.5；dev_s = s_range / 4 / dead_s；空间= linspace（0，s_range、节点）；%在空间划分采样点数锥度= 1 / sqrt（2×π）/ dev_s×exp（- 1 *（空间- P）。^ 2 / 2 / dev_s ^ 2）；锥度=锥度/最大（锥度）；数字（1）；情节（空间、锥度、C）；标题（“脉冲源”）；xlabel（X轴”）；ylabel（pluse”）；TE波的分量初始化%% ----------------------------------------前=零（节点，节点）；0（节点，节点）；赫兹=零（节点，节点）；%加入金属圆柱% --------------------------------center_s =圆（节点/ 2）；cn = center_s * del_s；PEC =人（节点、节点）；% *节点1节点的全矩阵对于k = 1：节点，对于j = 1：节点，拉德= sqrt（（K×del_s -CN）^ 2 +（J×del_s -CN）^ 2）；%加入圆柱也可以改为方柱如果（RAD = cylinder_rad），PEC（K，j）= 0；结束；结束；结束；数字（2）；轮廓（PEC）；%等高绘图标题（“金属圆柱”）；%计算参数设置% --------------------------------完成= 1；n = 0；f = 0；c_mu = del_t /亩/ del_s；c_eps = del_t /每股收益/ del_s；c_mur =（C del_t - del_s）/（C del_t + del_s）；帧= 80；数字（3）；轴（[ 0节点0节点0 1 ]）；fdtd_m =电影（框架）；而（完成= 0），%初始化ex_o =前；ey_o = Ey；hz_o =赫兹；%电场值的中心差分hterm = c_eps *（Hz（2：nodes-1,1：节点）赫兹（1：nodes-2,1：节点））；前（2 nodes-1,1：节点）=前（2：nodes-1,1：节点）+ hterm；hterm = c_eps *（Hz（1节点，1：nodes-2）赫兹（1节点，2：nodes-1））；安永（1节点，2：nodes-1）= Ey（1节点，2：nodes-1）+ hterm；%加入脉冲源T = N * del_t；术语=（t-mean）；脉冲=（1 / sqrt（2×π）/ dev ^ 3）*项；脉冲，脉冲* exp（（1 / 2 / dev ^ 2）*项^ 2）/ pulsenorm；源=脉冲*的（1，节点）*锥度；%的DBD校正% %可选择高斯脉冲t = 80；%项=（永久性的）；%脉冲=范围（4 * * * *术语2 /（10公顷）；%源=脉冲*的（1，节点）；前（2,1：节点）=前（2,1：节点）+源；% 4条边界线的处理%左边前（1,1：节点）= ex_o（2,1：节点）+ c_mur *（前（2,1：节点）-前（1,1：节点））；%上边安永（1：节点、节点）= ey_o（1：节点，n