look.m

三维FDTD,matlal编程,mur边界条件,平面波光源

clear
clc
nm=1e-9;
um=1e-6;
opengl neverselect

load coupling.mat
% r=[linspace(0,  0,32),linspace(0,1,32),linspace(1,1,32),linspace(1,0.5,32)];
% g=[linspace(0,  0,32),linspace(0,1,32),linspace(1,0,32),linspace(0,0,32)];
% b=[linspace(0.5,1,32),linspace(1,1,32),linspace(1,0,32),linspace(0,0,32)];
% mapp=[r',g',b'];

% mapp2=textread('dkbluered','','headerlines',2,'delimiter',' ');
% mapp=mapp2(:,1:end-1);

Ex=double(Exr+Exi*j);
Hy=double(Hyr+Hyi*j);
Ey=double(Eyr+Eyi*j);
Hx=double(Hxr+Hxi*j);
[nx,ny,nz]=size(Ex);
figure

% label=0;
label=1;


if label==0
    Sz=.5*real(Ex.*conj(Hy)-Ey.*conj(Hx));
pcolor(squeeze(real(Hy(:,ny/2+1,1:end))));

% colormap(mapp);
colormap jet
caxis([-5,5])
% colormap hot
% caxis([0,.5])

set(gca,'PlotBoxAspectRatio',[nz*dz,nx*dx,1])
xlabel('z/\Deltax')
ylabel('x/\Deltay')
shading interp
title('Ey ')
colorbar
set(gca,'TickDir','out')
% figure
% 
% plot(10*log(squeeze((Sz(47,ny/2+1,11:end)))))
%==========================================================================
else if label==1
pcolor(squeeze(real(Hy(44,:,1:end))))
colormap jet
caxis([-1,1])
% colormap(mapp)
% caxis([0,.5])

set(gca,'PlotBoxAspectRatio',[nz*dz,ny*dy,1])
xlabel('z/\Deltax')
ylabel('x/\Deltay')
shading interp
title('Ey ')
colorbar
set(gca,'TickDir','out')
%==========================================================================
else
pcolor(squeeze(abs(Hy(:,:,42))).^2)
% caxis([0,0.0005])
colormap(hot);

set(gca,'PlotBoxAspectRatio',[ny*dy,nx*dx,1])
xlabel('x/\Deltax')
ylabel('y/\Deltay')
% rectangle('position',[6,6,71,10],'edgecolor','white','linestyle','--')
% rectangle('position',[41-6,6,13,10],'edgecolor','red','linestyle','--')
shading interp
title('|Ey|^2')
colorbar
set(gca,'TickDir','out')
% 
% figure
% plot(squeeze(abs(Hx(41,:,80))))
end
end