% FDTD仿真二维TE波目标双站RCS；
% 编写者：tingo;
% Email:**********************.
%%
clear;
clf;
clc;
close all;
mu0=4*pi*(1e-7); %磁导率；
epsilon0=(1/36/pi)*(1e-9); %介电常数；
Vvacuum=1/sqrt(mu0*epsilon0); % 真空中的波传播速度；
Zc=sqrt(mu0/epsilon0);
namda=1e-2; % 1GHz对应的波长；
dx=namda/40;
dy=dx;
nx=150;
ny=nx;
Lx=nx*dx;
三二三事变Ly=Lx;
a=2*namda;
dt=dx/Vvacuum/2; % 时间步长；
itmax=1000; % 仿真时间步；
%% PML层电磁参数设置；
Rmax=1e-3; % 电磁波垂直入射的反射系数；
iebc=4; % 左右两边PML层网格数目；
jebc=4; % 上下两边PML层网格数目；
orderbc=4; % pML层中电导率分布阶数；
iebcf=nx+2*iebc;
iebcb=nx+2*iebc;
%% 为储存电磁数据的矩阵变量分配内存；
%========FDTD 区域=========%
Ex=zeros(nx,ny+1);
Ey=zeros(nx+1,ny);
Hz=zeros(nx,ny);
caex=zeros(nx,ny+1);
cbex=zeros(nx,ny+1);
caey=zeros(nx+1,ny);
cbey=zeros(nx+1,ny);
cahz=zeros(nx,ny);
cbhz=zeros(nx,ny);
%=========PML层==========%
% front;
Exbcf=zeros(iebcf,jebc);
Eybcf=zeros(iebcf+1,jebc);
Hzxbcf=zeros(iebcf,jebc);
Hzybcf=zeros(iebcf,jebc);
caexbcf=zeros(iebcf,jebc); %Exbcf计算系数1,
cbexbcf=zeros(iebcf,jebc); %Exbcf计算系数2；
caeybcf=zeros(iebcf+1,jebc);
cbeybcf=zeros(iebcf+1,jebc);
cahzxbcf=zeros(iebcf,jebc);
cbhzybcf=zeros(iebcf,jebc);
% back;
Exbcb=zeros(iebcb,jebc); % jebc+1
Eybcb=zeros(iebcb+1,jebc);
Hzxbcb=zeros(iebcb,jebc);
Hzybcb=zeros(iebcb,jebc);
caexbcb=zeros(iebcb,jebc);
cbexbcb=zeros(iebcb,jebc);
caeybcb=zeros(iebcb+1,jebc);
cbeybcb=zeros(iebcb+1,jebc);
cahzxbcb=zeros(iebcb,jebc);
cbhzxbcb=zeros(iebcb,jebc);
cahzybcb=zeros(iebcb,jebc);
cbhzybcb=zeros(iebcb,jebc);
% left
Exbcl=zeros(iebc,ny+1);
Eybcl=zeros(iebc,ny);
Hzxbcl=zeros(iebc,ny);
Hzybcl=zeros(iebc,ny);
红领巾小五年规划
caexbcl=zeros(iebc,ny+1);
比例尺
cbexbcl=zeros(iebc,ny+1);
caeybcl=zeros(iebc,ny);
cbeybcl=zeros(iebc,ny);
cahzxbcl=zeros(iebc,ny);
cbhzxbcl=zeros(iebc,ny);
cahzybcl=zeros(iebc,ny);
cbhzybcl=zeros(iebc,ny);
% right;
东升镇高级中学
Exbcr=zeros(iebc,ny+1);
Eybcr=zeros(iebc,ny);
Hzxbcr=zeros(iebc,ny);
Hzybcr=zeros(iebc,ny);
caexbcr=zeros(iebc,ny+1);
cbexbcr=zeros(iebc,ny+1);
caeybcr=zeros(iebc,ny);
cbeybcr=zeros(iebc,ny);
cahzxbcr=zeros(iebc,ny);
cbhzxbcr=zeros(iebc,ny);
cahzybcr=zeros(iebc,ny);
cbhzybcr=zeros(iebc,ny);
%% 仿真区域媒质参数设定；
nmedium=2; % 媒质数目；
中国式养老困局epsilonr=[1,1]; % 相对介电常数；
mur=[1,1]; % 相对磁导率
sigma=[0,6e7]; % 电导率；
rou=[0,0]; % 磁阻率；
朱时华ca=zeros(1,nmedium);
cb=zeros(1,nmedium);
cp=zeros(1,nmedium);
cq=zeros(1,nmedium);
for i=1:nmedium
epsilon=epsilon0*epsilonr(i);
mu=mu0*mur(i);
eaf=sigma(i)*dt/2/epsilon;
haf=rou(i)*dt/2/mu;
ca(i)=(1-eaf)/(1+eaf);
cb(i)=dt/epsilon/(1+eaf);
cp(i)=(1-haf)/(1+haf);
cq(i)=dt/mu/(1+haf);
end
%% 自由空间FDTD迭代计算系数；
% Ex;
caex(:,:)=ca(1);
cbex(:,:)=cb(1);
豫公网安备41160202000603
站长QQ:729038198
关于我们
投诉建议