clear;clc;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
IE=100;
JE=100;% 网格数量
L=2000; % 迭代次数
f=1.6e9; % 频率
dx=.01; % 网格尺寸
dt=dx/2/3e8; % 时间间隔
ic=3;%平面波源位置
e0=8.8e-12; % 真空介电常数
m0=pi*4e-7; % 真空磁导率
epsilon=1; % 相对介电常数
sigma=0; % 电导率
ia=7;
ib=IE-ia+1;
ja=7;
jb=JE-ja+1;
for i=1:IE
for j=1:JE
ga(i,j)=1/(epsilon+sigma*dt/e0); % 求和参量
gb(i,j)=sigma*dt/e0; % 求和参量
end
end
ez_inc=zeros(1,IE);
hx_inc=zeros(1,JE);
ez_low1=0;
本内容试读结束,登录后可阅读更多
下载后可阅读完整内容,剩余3页未读,立即下载