仿真在介质层有损耗的传输的电磁场入射反射折射的3D图_matlab

clear cc = 3.0e8; %speed of light in the vacuum eps0 = 8.8546e-12; mu_0= 4.0*pi*1.0e-7; % Permeability of free space mu_0*mu_r f0=1e12;% 1e6, 1e12; sample_rate=30; dt=1/f0/sample_rate; nmax=50; sigma=100; % 10 epsr1=1-1i*sigma/(2*pi*f0*eps0); % epsr1=(1-4000i); epsl =eps0 *epsr1; mu_r=1;%2.25; mu=mu_0*mu_r; aimp0 = sqrt(mu_0/eps0); aimp0over1 = sqrt(epsr1/mu_r); R01=(1-aimp0over1)/(1+aimp0over1); T01=2/(1+aimp0over1); R12=(aimp0over1-1)/(1+aimp0over1); T12=2*aimp0over1/(1+aimp0over1); R10=R12; T10=T12; k0=2*pi*f0*sqrt(mu_0*eps0); k1=2*pi*f0*sqrt(mu*epsl); kx0=0.6*k0; kz0=0.8*k0; kx1=kx0; kz1=sqrt(k1*k1-kx1*kx1); size_x=50; size_z=150; boundary1=50; boundary2=100; dx=dt*cc; dz=dx; Ey_i=zeros(size_x,size_z); Ey_r=zeros(size_x,size_z); Ey_t=zeros(size_x,size_z); Ey=zeros(size_x,size_z); for n=1:nmax for x=1:size_x for z=1:boundary1 Ey_i(x,z)=exp(1i*(2*pi*f0*n*dt-kx0*dx*x-kz0*dx*(z-boundary1))); Ey_r(x,z)=R01*exp(1i*(2*pi*f0*n*dt-kx0*dx*x+kz0*dx*(z-boundary1))); Ey_t(x,z)=T01*R12*T10*exp(1i*(2*pi*f0*n*dt-kx0*dx*x+kz0*dx*(z-boundary1))); Ey(x,z)=Ey_i(x,z)+Ey_r(x,z)+Ey_t(x,z); end for z=boundary1+1:boundary2 Ey_i(x,z)=T01*exp(1i*(2*pi*f0*n*dt-kx0*dx*x-kz1*dz*(z-boundary1))); Ey_r(x,z)=T01*R12*exp(1i*(2*pi*f0*n*dt-kx0*dx*x-kz1*dz*(boundary2-boundary1)+kz1*dz*(z-boundary2))); Ey(x,z)=Ey_i(x,z)+Ey_r(x,z); end for z=boundary2+1:size_z Ey(x,z)=T01*T12*exp(1i*(2*pi*f0*n*dt-kx0*dx*x-kz1*dz*(boundary2-boundary1)-kz0*dz*(z-boundary2))); end end pcolor(real(Ey)) shading interp caxis([-0.8 0.8]) colorbar pause(0.1) end